Investigating biological systems using first principles Car-Parrinello molecular dynamics simulations

Current Opinion in Structural Biology

Volumn 17, Issue 2, 2007, Pages 149-156

  Dal Peraro, Matteo ^a   Ruggerone, Paolo ^b   Raugei, Simone ^c   Gervasio, Francesco Luigi ^d   Carloni, Paolo ^c  

^a UNIVERSITY OF PENNSYLVANIA   (United States)

^b UNIVERSITY OF CAGLIARI   (Italy)

^c SISSA   (Italy)

^d ETH ZURICH   (Switzerland)

Author keywords

[No Author keywords available]

Indexed keywords

ANTINEOPLASTIC AGENT; CISPLATIN;

BIOMOLECULAR ELECTRONICS; DNA ADDUCT; DRUG BINDING; DRUG DNA INTERACTION; DRUG TARGETING; ELECTRON TRANSPORT; HUMAN; MOLECULAR BIOLOGY; MOLECULAR DYNAMICS; MOLECULAR INTERACTION; MOLECULAR MECHANICS; PRIORITY JOURNAL; REVIEW; SIMULATION; SYSTEM ANALYSIS; SYSTEMS THEORY;

ANTINEOPLASTIC AGENTS; ELECTROCHEMISTRY; ELECTRON TRANSPORT; ENZYMES; LIGANDS; METALLOPROTEINS; MODELS, BIOLOGICAL; MODELS, MOLECULAR; OXIDATION-REDUCTION; PLATINUM COMPOUNDS; SYSTEMS BIOLOGY; THERMODYNAMICS;

EID: 34147185256     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.sbi.2007.03.018     Document Type: Review

References (72)

1. Car R., and Parrinello M. Unified approach for molecular-dynamics and density-functional theory. Phys Rev Lett 55 (1985) 2471-2474
2. VandeVondele J., Krack M., Mohamed F., Parrinello M., Chassaing T., and Hutter J. QUICKSTEP: fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput Phys Commun 167 (2005) 103-128
3. Warshel A., and Levitt M. Theoretical studies of enzymic reactions - dielectric, electrostatic and steric stabilization of carbonium-ion in reaction of lysozyme. J Mol Biol 103 (1976) 227-249
4. Laio A., VandeVondele J., and Rothlisberger U. A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations. J Chem Phys 116 (2002) 6941-6947
5. Eichinger M., Tavan P., Hutter J., and Parrinello M. A hybrid method for solutes in complex solvents: density functional theory combined with empirical force fields. J Chem Phys 110 (1999) 10452-10467
6. Mordasini T., Curioni A., and Andreoni W. Why do divalent metal ions either promote or inhibit enzymatic reactions? - The case of BamHI restriction endonuclease from combined quantum-classical simulations. J Biol Chem 278 (2003) 4381-4384
7. van Gunsteren WF, Billeter SR, Eising AA, Hünenberger PH,Krüger PKHC, Mark AE, Scott WRP, Tironi IG: Biomolecular simulation: The GROMOS96 manual and user guide. Zurich: vdf Hochschulverlag AG, 1996.
8. Dellago C., Bolhuis P.G., and Chandler D. On the calculation of reaction rate constants in the transition path ensemble. J Chem Phys 110 (1999) 6617-6625
9. Sprik M., and Ciccotti G. Free energy from constrained molecular dynamics. J Chem Phys 109 (1998) 7737-7744
10. Laio A., and Parrinello M. Escaping free-energy minima. Proc Natl Acad Sci USA 99 (2002) 12562-12566
11. Iannuzzi M., Laio A., and Parrinello M. Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. Phys Rev Lett 90 (2003) 238302
12. Raugei S., Cascella M., and Carloni P. A proficient enzyme: insights on the mechanism of orotidine monophosphate decarboxylase from computer simulations. J Am Chem Soc 126 (2004) 15730-15737
13. Cascella M., Micheletti C., Rothlisberger U., and Carloni P. Evolutionarily conserved functional mechanics across pepsin-like and retroviral aspartic proteases. J Am Chem Soc 127 (2005) 3734-3742. The combination of multinanosecond MD and MM/CPMD simulations indicates that mechanical large-scale fluctuations exhibited by proteases are crucial for the modulation of the catalytic efficiency.
14. Sulpizi M., Laio A., VandeVondele J., Cattaneo A., Rothlisberger U., and Carloni P. Reaction mechanism of caspases: insights from QM/MM Car-Parrinello simulations. Proteins 52 (2003) 212-224
15. Biarnes X., Nieto J., Planas A., and Rovira C. Substrate distortion in the Michaelis complex of Bacillus 1,3-1,4-beta-glucanase - insight from first principles molecular dynamics simulations. J Biol Chem 281 (2006) 1432-1441
16. Vidossich P., Piana S., Miani A., and Carloni P. Deuterium isotope effects in A:T and A:U base pairs: a computational NMR study. J Am Chem Soc 128 (2006) 7215-7221
17. + selective filter suggest that the backbone carbonyl groups are largely polarized and that there is sizeable charge transfer from the backbone to the cations. These effects might be relevant for an accurate description of the electrostatic repulsion among the carbonyl groups and the cations, which has been proposed to play a central role in the selectivity mechanism of the channel [67].
18. Dal Peraro M., Llarrull L.I., Rothlisberger U., Vila A.J., and Carloni P. Water-assisted reaction mechanism of monozinc β-lactamases. J Am Chem Soc 126 (2004) 12661-12668
19. 2+ ion in MβLs is discussed based on a direct and consistent comparison of monozinc and dizinc species.
20. De Vivo M., Ensing B., and Klein M.L. Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle. J Am Chem Soc 127 (2005) 11226-11227
21. De Vivo M., Ensing B., Dal Peraro M., Gomez G.A., Christianson D.W., and Klein M.L. Proton shuttles and phosphatase activity in soluble epoxide hydrolase. J Am Chem Soc 129 (2007) 387-394. The crucial role of an ordered hydrogen-bond network of water molecules in phosphoryl transfer reactions emerges from this set of calculations. This suggests the general importance of explicitly including solvation effects at the catalytic site for the correct description of a reaction mechanism.
22. + cations at the active site cooperate to activate this water via exchange between their solvation shells, sensibly improving the catalytic efficiency of formation of the inorganic phosphate.
23. Rovira C., Alfonso-Prieto M., Biarnes X., Carpena X., Fita I., and Loewen P.C. A first principles study of the binding of formic acid in catalase complementing high resolution X-ray structures. Chem Phys 323 (2006) 129-137
24. Alfonso-Prieto M., Borovik A., Carpena X., Mursnudov G., Melik-Adamyan W., Fita I., Rovira C., and Loewen P.C. The structures and electronic configuration of compound I intermediates of Helicobacter pylori and Penicillium vitale catalases determined by X-ray crystallography and QM/MM DFT calculations. J Am Chem Soc (2007) 10.1021/ja063660y
25. Raugei S., and Carloni P. Structure and function of vanadium haloperoxidases. J Phys Chem B Condens Matter Mater Surf Interfaces Biophys 110 (2006) 3747-3758. The rest state of the vanadium-dependent chloroperoxidases and the early stages of the catalytic cycle have been investigated by MM/CPMD metadynamics simulations. This is one of the first applications of metadynamics to a biological system within the MM/CPMD approach.
26. Boero M., Tateno M., Terakura K., and Oshiyama A. Double-metal-ion/single-metal-ion mechanisms of the cleavage reaction of ribozymes: first-principles molecular dynamics simulations of a fully hydrated model system. J Chem Theory Comput 1 (2005) 925-934
27. Crowder M.W., Spencer J., and Vila A.J. Metallo-β-lactamases: novel weaponry for antibiotic resistance in bacteria. Acc Chem Res 39 (2006) 721-728
28. Jamieson E.R., and Lippard S.J. Structure, recognition, and processing of cisplatin-DNA adducts. Chem Rev 99 (1999) 2467-2498
29. Reedijk J. New clues for platinum antitumor chemistry: kinetically controlled metal binding to DNA. Proc Natl Acad Sci USA 100 (2003) 3611-3616
30. Spiegel K., Rothlisberger U., and Carloni P. Cisplatin binding to DNA oligomers from hybrid Car-Parrinello/molecular dynamics simulations. J Phys Chem B 108 (2004) 2699-2707. Extensive MM/CPMD simulations of cisplatin-DNA adducts succeed in shedding light on questions related to structural features extracted from X-ray and NMR data. In particular, the platinated DNA has been simulated in complex with the high mobility group protein domain A, pinpointing electronic and structural perturbations induced by complexation.
31. Gelasco A., and Lippard S.J. NMR solution structure of a DNA dodecamer duplex containing a cis-diammineplatinum(II) d(GpG) intrastrand cross-link, the major adduct of the anticancer drug cisplatin. Biochemistry 37 (1998) 9230-9239
32. Magistrato A., Ruggerone P., Spiegel K., Carloni P., and Reedijk J. Binding of novel azole-bridged dinuclear platinum(II) anticancer drugs to DNA: insights from hybrid QM/MM molecular dynamics simulations. J Phys Chem B Condens Matter Mater Surf Interfaces Biophys 110 (2006) 3604-3613. MM/CPMD calculations have identified the electronic origin of the isomerization process involving a new platinum-derived drug, confirming the suggestion made on the basis of experimental data. The process enlarges the range of possible drug-DNA linkages and could cooperate in lowering drug resistance.
33. Allardyce C.S., Dyson P.J., Ellis D.J., and Heath S.L. [Ru(eta(6)-p-cymene)Cl-2(pta)] (pta=1,3,5-triaza-7-phosphatricyclo[3. 3. 1. 1]decane): a water soluble compound that exhibits pH dependent DNA binding providing selectivity for diseased cells. Chem Commun (2001) 1396-1397
34. Gossens C., Tavemelli I., and Rothlisberger U. Rational design of organo-ruthenium anticancer compounds. Chimia (Aarau) 59 (2005) 81-84
35. Novakova O., Chen H.M., Vrana O., Rodger A., Sadler P.J., and Brabec V. DNA interactions of monofunctional organometallic ruthenium(II) antitumor complexes in cell-free media. Biochemistry 42 (2003) 11544-11554
36. Vidossich P., and Carloni P. Binding of phosphinate and phosphonate inhibitors to aspartic proteases: a first-principles study. J Phys Chem B Condens Matter Mater Surf Interfaces Biophys 110 (2006) 1437-1442
37. 2 discrimination in myoglobin. J Phys Chem B Condens Matter Mater Surf Interfaces Biophys. 109 (2005) 3065-3070
38. Rovira C. Study of ligand-protein interactions by means of density functional theory and first-principles molecular dynamics. Methods Mol Biol 305 (2005) 517-554
39. 2 vs CO-heme: a time-dependent density functional study. J Am Chem Soc 125 (2003) 15710-15711
40. Kulik H.J., Cococcioni M., Scherlis D.A., and Marzari N. Density functional theory in transition-metal chemistry: a self-consistent Hubbard U approach. Phys Rev Lett 97 (2006) 103001. Transition metal complexes, usually problematic within DFT, benefit from a generalized-gradient approximation augmented by a Hubbard term (U), which is determined by a new self-consistent procedure based on linear response theory. This computationally efficient scheme provides remarkable improvements for spin energetics, state symmetries and reaction barriers, and represents a promising step towards an accurate theoretical description of transition-metal-based biochemistry.
41. Degtyarenko I., Nieminen R.M., and Rovira C. Structure and dynamics of dioxygen bound to cobalt and iron heme. Biophys J 91 (2006) 2024-2034
42. Schuster G.B., and Landman U. The mechanism of long-distance radical cation transport in duplex DNA: ion-gated hopping of polaron-like distortions. Top Curr Chem 236 (2004) 139-161
43. Gervasio F.L., Carloni P., and Parrinello M. Electronic structure of wet DNA. Phys Rev Lett 89 (2002) 108102
44. Hubsch A., Endres R.G., Cox D.L., and Singh R.R.P. Optical conductivity of wet DNA. Phys Rev Lett 94 (2005) 178102
45. Di Felice R., Calzolari A., Garbesi A., Alexandre S.S., and Soler J.M. Strain-dependence of the electronic properties in periodic quadruple helical G4-wires. J Phys Chem B Condens Matter Mater Surf Interfaces Biophys 109 (2005) 22301-22307
46. Gervasio F.L., Laio A., Parrinello M., and Boero M. Charge localization in DNA fibers. Phys Rev Lett 94 (2005) 158103
47. Gervasio F.L., Laio A., Iannuzzi M., and Parrinello M. Influence of DNA structure on the reactivity of the guanine radical cation. Chemistry 10 (2004) 4846-4852
48. Mundy C.J., Colvin M.E., and Quong A.A. Irradiated guanine: a Car-Parrinello molecular dynamics study of dehydrogenation in the presence of an OH radical. J Phys Chem A 106 (2002) 10063-10071
49. Marcus R.A. Theory of oxidation-reduction reactions involving electron transfer.1. J Chem Phys 24 (1956) 966-978
50. Marcus R.A. On theory of electron-transfer reactions. 6. Unified treatment for homogeneous and electrode reactions. J Chem Phys 43 (1965) 679
51. Warshel A. Dynamics of reactions in polar-solvents - semi-classical trajectory studies of electron-transfer and proton-transfer reactions. J Phys Chem 86 (1982) 2218-2224
52. Tavernelli I., Vuilleumier R., and Sprik M. Ab initio molecular dynamics for molecules with variable numbers of electrons. Phys Rev Lett 88 (2002) 213002
53. 2+ redox couple: a biased Born-Oppenheimer molecular dynamics investigation. J Chem Phys 124 (2006) 64507
54. Tateyama Y., Blumberger J., Sprik M., and Tavernelli I. Density-functional molecular-dynamics study of the redox reactions of two anionic, aqueous transition-metal complexes. J Chem Phys 122 (2005) 234505
55. Blumberger J., Bernasconi L., Tavernelli I., Vuilleumier R., and Sprik M. Electronic structure and solvation of copper and silver ions: a theoretical picture of a model aqueous redox reaction. J Am Chem Soc 126 (2004) 3928-3938
56. Cascella M., Magistrato A., Tavernelli I., Carloni P., and Rothlisberger U. Role of protein frame and solvent for the redox properties of azurin from Pseudomonas aeruginosa. Proc Natl Acad Sci USA 103 (2006) 19641-19646. MM/CPMD simulations coupled to a grand-canonical scheme have been performed to study the Cu(II)/Cu(I) half-reaction in azurin from P. aeruginosa. An accurate description at atomistic level of the environment surrounding the metal-binding site and thermal fluctuations of the protein structure are both essential for a correct quantitative description of the redox properties of the system. The redox potential is fine-tuned by long-range electrostatic interactions, whereas the reorganization energy is mainly due to solvent rearrangement at the protein surface.
57. Blumberger J., and Klein M.L. Reorganization free energies for long-range electron transfer in a porphyrin-binding four-helix bundle protein. J Am Chem Soc 128 (2006) 13854-13867. The authors have computed the reorganization free energy for both oxidation of a single Ru-porphyrin cofactor and electron self-exchange between two Ru-porphyrin cofactors binding to a solvated four-helix-bundle protein. Through accurate sampling of the classical electrostatic energy gap, they find that the fluctuations are well described by Gaussian statistics. Decomposing the outer-sphere reorganization free energy, they found that the solvent (aqueous ionic solution) is the primary outer-sphere medium for oxidation, whereas the protein contributes only for a smaller fraction.
58. Sulpizi M, Raugei S, VandeVondele J, Sprik M, Carloni P: Calculation of redox properties: understanding short and long range effects in rubredoxin. J Phys Chem B 2007, in press. The experimental redox potential difference between mesophilic and thermophilic rubredoxin is explained in terms of short-range contributions from a few residues close to the metal center. It is also proposed that the protein controls the reorganization free energy in a balance between shielding of the high dielectric solvent medium and preserving some water accessibility to the metal center for fast electron transfer to the surface. These results underline the importance of a molecular description of the solvent and of correct inclusion of the polarization effects.
59. 2: a self-interaction-corrected approach. Phys Rev B 71 (2005) 205210
60. Baletto F., Cavazzoni C., and Scandolo S. Surface trapped excess electrons on ice. Phys Rev Lett 95 (2005) 176801
61. VandeVondele J., and Sprik M. A molecular dynamics study of the hydroxyl radical in solution applying self-interaction-corrected density functional methods. Phys Chem Chem Phys 7 (2005) 1363-1367
62. Maurer P., Laio A., Hugossonz H., Colombo M.C., and Rothlisberger U. Automated parametrization of biomolecular force fields from QM/MM simulations through force-matching. J Chem Theory Comput 3 (2007) 628-639. A recipe to parameterize biomolecular systems 'on the fly' via first principles MM/CPMD simulations is presented in this paper. The generation of DFT-based parameters for simulating metal-containing proteins for which standard force fields poorly perform can largely improve the accuracy of current MD simulations.
63. Spiegel K, Magistrato A, Maurer P, Ruggerone P, Rothlisberger U, Carloni P, Reedijk J, Klein ML: Parameterization of azole-bridged dinuclear platinum anticancer drugs via a QM/MM force matching procedure. J Comput Chem 2007, in press.
64. Izvekov S., Parrinello M., Burnham C.J., and Voth G.A. Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching. J Chem Phys 120 (2004) 10896-10913
65. von Lilienfeld O.A., Tavernelli I., Rothlisberger U., and Sebastiani D. Optimization of effective atom centered potentials for London dispersion forces in density functional theory. Phys Rev Lett 93 (2004) 153004. Using an efficient algorithm, dispersion forces are included in DFT-based calculations at essentially no additional computational cost. The van der Waals forces are modeled by an atom-electron interaction, instead of approximating them via an atom-atom term. Thus, a remarkable improvement in the description of molecular properties within the DFT framework is achieved.
66. Laino T., Mohamed F., Laio A., and Parrinello M. An efficient real space multigrid QM/MM electrostatic coupling. J Chem Theory Comput 1 (2005) 1176-1184
67. Noskov S.Y., Berneche S., and Roux B. Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands. Nature 431 (2004) 830-833
68. Komeda S., Lutz M., Spek A.L., Chikuma M., and Reedijk J. New antitumor-active azole-bridged dinuclear platinum(II) complexes: synthesis, characterization, crystal structures, and cytotoxic studies. Inorg Chem 39 (2000) 4230-4236
69. Barnett R.N., Cleveland C.L., Joy A., Landman U., and Schuster G.B. Charge migration in DNA: ion-gated transport. Science 294 (2001) 567-571
70. Gervasio F.L., Boero M., and Parrinello M. Double proton coupled charge transfer in DNA. Angew Chem Int Ed Engl 45 (2006) 5606-5609. Based on CPMD and MM/CPMD calculations, the authors point to a key role for rare DNA tautomers in charge transfer processes in DNA.
71. Giese B., Amaudrut J., Kohler A.K., Spormann M., and Wessely S. Direct observation of hole transfer through DNA by hopping between adenine bases and by tunnelling. Nature 412 (2001) 318-320
72. Kühne T.D., Krack M., Mohamed F.R., and Parrinello M. Efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. Phys Rev Lett 98 (2007) 066401. The authors introduce an improved dynamics in which the CP equation of motion is replaced by an integration scheme (via an accurate predictor-corrector algorithm) to maintain the electrons as close as possible to the instantaneous ground state. The resulting dynamics is very accurate in terms of sampling distribution and deviation from the Born-Oppenheimer (BO) surface, and accelerated by at least one order of magnitude with respect to current CPMD and BOMD.