G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. Their structure and function are strongly influenced by their hydration shells. Although...G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. Their structure and function are strongly influenced by their hydration shells. Although extensively studied through various experimental and computational methods, hydration patterns near DNA remain under debate due to the chemically and topologically heterogeneous nature of the exposed surface. In this work, we employed all-atom molecular dynamics(MD) simulation to study the hydration patterns of GQ DNA. The Drude oscillator model was used in MD simulation as a computationally efficient method for modeling electronic polarization in DNA ion solutions. Hydration structure was analyzed in terms of radial distribution functions and high-density three-dimensional hydration sites. Analysis of hydration dynamics focused on self-diffusion rates and orientation time correlation at different structural regions of GQ DNA.The results show highly heterogeneous hydration patterns in both structure and dynamics;for example, there are several insular high-density sites in the inner channel, and ‘spine of water’ in the groove. For water inside the loop, anomalous diffusion is present over a long time scale, but for water around the phosphate group and groove, diffusion becomes normal after ~30 ps. These essentially correspond to deeply buried structural water and strong interaction with DNA, respectively.展开更多
G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. The topology of GQs is associated with the sequences and lengths of DNA, the types of li...G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. The topology of GQs is associated with the sequences and lengths of DNA, the types of linking loops, and the associated metal cations. However, our understanding on the basic physical properties of the formation process and the stability of GQs is rather limited. In this work, we employed ab initio, molecular dynamics(MD), and steered MD(SMD)simulations to study the interaction between loop bases and ions, and the effect on the stability of G-quadruplex DNA, the Drude oscillator model was used in MD and SMD simulations as a computationally efficient manner method for modeling electronic polarization in DNA ion solutions. We observed that the binding energy between DNA bases and ions(K^(+)/Na^(+))is about the base stacking free energies indicates that there will be a competition among the binding of M^(+)-base, H-bonds between bases, and the base-stacking while ions were bound in loop of GQs. Our SMD simulations indicated that the side loop inclined to form the base stacking while the loop sequence was Thy or Ade, and the cross-link loop upon the G-tetrads was not easy to form the base stacking. The base stacking side loop complex K+was found to have a good stabilization synergy. Although a stronger interaction was observed to exist between Cyt and K+, such an interaction was unable to promote the stability of the loop with the sequence Cyt.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 11705160 and 11647074)。
文摘G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. Their structure and function are strongly influenced by their hydration shells. Although extensively studied through various experimental and computational methods, hydration patterns near DNA remain under debate due to the chemically and topologically heterogeneous nature of the exposed surface. In this work, we employed all-atom molecular dynamics(MD) simulation to study the hydration patterns of GQ DNA. The Drude oscillator model was used in MD simulation as a computationally efficient method for modeling electronic polarization in DNA ion solutions. Hydration structure was analyzed in terms of radial distribution functions and high-density three-dimensional hydration sites. Analysis of hydration dynamics focused on self-diffusion rates and orientation time correlation at different structural regions of GQ DNA.The results show highly heterogeneous hydration patterns in both structure and dynamics;for example, there are several insular high-density sites in the inner channel, and ‘spine of water’ in the groove. For water inside the loop, anomalous diffusion is present over a long time scale, but for water around the phosphate group and groove, diffusion becomes normal after ~30 ps. These essentially correspond to deeply buried structural water and strong interaction with DNA, respectively.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11705160 and 11647074)。
文摘G-quadruplexes(GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. The topology of GQs is associated with the sequences and lengths of DNA, the types of linking loops, and the associated metal cations. However, our understanding on the basic physical properties of the formation process and the stability of GQs is rather limited. In this work, we employed ab initio, molecular dynamics(MD), and steered MD(SMD)simulations to study the interaction between loop bases and ions, and the effect on the stability of G-quadruplex DNA, the Drude oscillator model was used in MD and SMD simulations as a computationally efficient manner method for modeling electronic polarization in DNA ion solutions. We observed that the binding energy between DNA bases and ions(K^(+)/Na^(+))is about the base stacking free energies indicates that there will be a competition among the binding of M^(+)-base, H-bonds between bases, and the base-stacking while ions were bound in loop of GQs. Our SMD simulations indicated that the side loop inclined to form the base stacking while the loop sequence was Thy or Ade, and the cross-link loop upon the G-tetrads was not easy to form the base stacking. The base stacking side loop complex K+was found to have a good stabilization synergy. Although a stronger interaction was observed to exist between Cyt and K+, such an interaction was unable to promote the stability of the loop with the sequence Cyt.
