Molecular Dynamics Simulation Study of Parallel Telomeric DNA Quadruplexes at Different Ionic Strengths: Evaluation of Water and Ion Models

Warning

This publication doesn't include Faculty of Arts. It includes Central European Institute of Technology. Official publication website can be found on muni.cz.
Authors

REBIČ Matúš LAAKSONEN Aatto ŠPONER Jiří ULIČNÝ Jozef MOCCI Francesca

Year of publication 2016
Type Article in Periodical
Magazine / Source Journal of Physical Chemistry B
MU Faculty or unit

Central European Institute of Technology

Citation
web http://pubs.acs.org/doi/pdf/10.1021/acs.jpcb.6b06485
Doi http://dx.doi.org/10.1021/acs.jpcb.6b06485
Field Physical chemistry and theoretical chemistry
Keywords AMBER FORCE-FIELD; NUCLEIC-ACIDS; BIOMOLECULAR SIMULATIONS; LIQUID WATER; CONFORMATIONAL DYNAMICS; COMPUTER-SIMULATIONS; POTENTIAL FUNCTIONS; FREE-ENERGY; RNA; PARAMETERS
Description Most molecular dynamics (MD) simulations of DNA quadruplexes have been performed under minimal salt conditions using the Aqvist potential parameters for the cation with the TIP3P water model. Recently, this combination of parameters has been reported to be problematic for the stability of quadruplex DNA, especially caused by the ion interactions inside or near the quadruplex channel. Here, we verify how the choice of ion parameters and water model can affect the quadruplex structural stability and the interactions with the ions outside the channel. We have performed a series of MD simulations of the human full-parallel telomeric quadruplex by neutralizing its negative charge with K+ ions. Three combinations of different cation potential parameters and water models have been used: (a) Aqvist ion parameters, TIP3P water model; (b) Joung and Cheatham ion parameters, TIP3P water model; and (c) Joung and Cheatham ion parameters, TTP4P(ew) water model. For the combinations (b) and (c), the effect of the ionic strength has been evaluated by adding increasing amounts of KCl salt (50, 100, and 200 mM). Two independent simulations using the Aqvist parameters with the TIP3P model show that this combination is clearly less suited for the studied quadruplex with K+ counterions. In both simulations, one ion escapes from the channel, followed by significant deformation of the structure, leading to deviating conformation compared to that in the reference crystallographic data. For the other combinations of ion and water potentials, no tendency is observed for the channel ions to escape from the quadruplex channel. In addition, the internal mobility of the three loops, torsion angles, and counterion affinity have been investigated at varied salt concentrations. In summary, the selection of ion and water models is crucial as it can affect both the structure and dynamics as well as the interactions of the quadruplex with its counterions. The results obtained with the TIP4P(ew) model are found to be closest to the experimental data at all of the studied ion concentrations.
Related projects:

You are running an old browser version. We recommend updating your browser to its latest version.