Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning

Overview

Journal J Chem Theory Comput

Publisher American Chemical Society

Specialties Biochemistry
Chemistry

Date 2024 May 14

PMID 38742922

Authors

Xudong Yang

Chengwen Liu

Pengyu Ren

Affiliations

Soon will be listed here.

Abstract

Conformational dynamics play a crucial role in determining the behavior of the biomolecules. Polarizable force fields, such as AMOEBA, can accurately capture electrostatic interactions underlying the conformational space. However, applying a polarizable force field in molecular dynamics (MD) simulations can be computationally expensive, especially in studying long-time-scale dynamics. To overcome this challenge, we incorporated the AMOEBA potential with Milestoning, an enhanced sampling method in this work. This integration allows us to efficiently sample the rare and important conformational states of a biomolecule by using many short and independent molecular dynamics trajectories with the AMOEBA force field. We applied this method to investigate the conformational dynamics of alanine dipeptide, DNA, and RNA A-B form conversion. Well-converged thermodynamic and kinetic properties were obtained, including the free energy difference, mean first passage time, and critical transitions between states. Our results demonstrate the power of integrating polarizable force fields with enhanced sampling methods in quantifying the thermodynamic and kinetic properties of biomolecules at the atomic level.

References

Jose D, Porschke D . The dynamics of the B-A transition of natural DNA double helices. J Am Chem Soc. 2005; 127(46):16120-8. DOI: 10.1021/ja053691d. View

Zuckerman D, Chong L . Weighted Ensemble Simulation: Review of Methodology, Applications, and Software. Annu Rev Biophys. 2017; 46:43-57. PMC: 5896317. DOI: 10.1146/annurev-biophys-070816-033834. View

Li S, Olson W, Lu X . Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures. Nucleic Acids Res. 2019; 47(W1):W26-W34. PMC: 6602438. DOI: 10.1093/nar/gkz394. View

Porschke D . Kinetics of the B-A transition of DNA: analysis of potential contributions to a reaction barrier. Eur Biophys J. 2018; 47(4):325-332. PMC: 5982448. DOI: 10.1007/s00249-018-1276-4. View

Romo T, Leioatts N, Grossfield A . Lightweight object oriented structure analysis: tools for building tools to analyze molecular dynamics simulations. J Comput Chem. 2014; 35(32):2305-18. PMC: 4227929. DOI: 10.1002/jcc.23753. View

Templeton C, Elber R . Why Does RNA Collapse? The Importance of Water in a Simulation Study of Helix-Junction-Helix Systems. J Am Chem Soc. 2018; 140(49):16948-16951. PMC: 6470002. DOI: 10.1021/jacs.8b11111. View

Knee K, Dixit S, Aitken C, Ponomarev S, Beveridge D, Mukerji I . Spectroscopic and molecular dynamics evidence for a sequential mechanism for the A-to-B transition in DNA. Biophys J. 2008; 95(1):257-72. PMC: 2426640. DOI: 10.1529/biophysj.107.117606. View

Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C . Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins. 2006; 65(3):712-25. PMC: 4805110. DOI: 10.1002/prot.21123. View

Manavalan B, Kuwajima K, Lee J . PFDB: A standardized protein folding database with temperature correction. Sci Rep. 2019; 9(1):1588. PMC: 6367381. DOI: 10.1038/s41598-018-36992-y. View

10.

DICKERSON R . Definitions and nomenclature of nucleic acid structure components. Nucleic Acids Res. 1989; 17(5):1797-803. PMC: 317523. DOI: 10.1093/nar/17.5.1797. View

11.

Ng H, Dickerson R . Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG. Nucleic Acids Res. 2002; 30(18):4061-7. PMC: 137100. DOI: 10.1093/nar/gkf515. View

12.

Kulkarni M, Mukherjee A . Sequence dependent free energy profiles of localized B- to A-form transition of DNA in water. J Chem Phys. 2013; 139(15):155102. DOI: 10.1063/1.4825175. View

13.

Invernizzi M, Parrinello M . Rethinking Metadynamics: From Bias Potentials to Probability Distributions. J Phys Chem Lett. 2020; 11(7):2731-2736. DOI: 10.1021/acs.jpclett.0c00497. View

14.

Noy A, Perez A, Laughton C, Orozco M . Theoretical study of large conformational transitions in DNA: the B<-->A conformational change in water and ethanol/water. Nucleic Acids Res. 2007; 35(10):3330-8. PMC: 1904281. DOI: 10.1093/nar/gkl1135. View

15.

Zgarbova M, Jurecka P, Sponer J, Otyepka M . A- to B-DNA Transition in AMBER Force Fields and Its Coupling to Sugar Pucker. J Chem Theory Comput. 2017; 14(1):319-328. DOI: 10.1021/acs.jctc.7b00926. View

16.

Rackers J, Wang Z, Lu C, Laury M, Lagardere L, Schnieders M . Tinker 8: Software Tools for Molecular Design. J Chem Theory Comput. 2018; 14(10):5273-5289. PMC: 6335969. DOI: 10.1021/acs.jctc.8b00529. View

17.

Kovaleva N, Strelnikov I, Zubova E . Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer. ACS Omega. 2021; 5(51):32995-33006. PMC: 7774075. DOI: 10.1021/acsomega.0c04247. View

18.

Ivanov P, Wadhams M, Roby R, Holland M, Weedn V, Parsons T . Mitochondrial DNA sequence heteroplasmy in the Grand Duke of Russia Georgij Romanov establishes the authenticity of the remains of Tsar Nicholas II. Nat Genet. 1996; 12(4):417-20. DOI: 10.1038/ng0496-417. View

19.

Shi Y, Xia Z, Zhang J, Best R, Wu C, Ponder J . The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins. J Chem Theory Comput. 2013; 9(9):4046-4063. PMC: 3806652. DOI: 10.1021/ct4003702. View

20.

Hollingsworth S, Dror R . Molecular Dynamics Simulation for All. Neuron. 2018; 99(6):1129-1143. PMC: 6209097. DOI: 10.1016/j.neuron.2018.08.011. View