1
|
Christiansen H, Errica F, Alesiani F. Self-tuning Hamiltonian Monte Carlo for accelerated sampling. J Chem Phys 2023; 159:234109. [PMID: 38108481 DOI: 10.1063/5.0177738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2023] [Accepted: 11/27/2023] [Indexed: 12/19/2023] Open
Abstract
The performance of Hamiltonian Monte Carlo simulations crucially depends on both the integration timestep and the number of integration steps. We present an adaptive general-purpose framework to automatically tune such parameters based on a local loss function that promotes the fast exploration of phase space. We show that a good correspondence between loss and autocorrelation time can be established, allowing for gradient-based optimization using a fully differentiable set-up. The loss is constructed in such a way that it also allows for gradient-driven learning of a distribution over the number of integration steps. Our approach is demonstrated for the one-dimensional harmonic oscillator and alanine dipeptide, a small protein commonly used as a test case for simulation methods. Through the application to the harmonic oscillator, we highlight the importance of not using a fixed timestep to avoid a rugged loss surface with many local minima, otherwise trapping the optimization. In the case of alanine dipeptide, by tuning the only free parameter of our loss definition, we find a good correspondence between it and the autocorrelation times, resulting in a >100 fold speedup in the optimization of simulation parameters compared to a grid search. For this system, we also extend the integrator to allow for atom-dependent timesteps, providing a further reduction of 25% in autocorrelation times.
Collapse
Affiliation(s)
- Henrik Christiansen
- NEC Laboratories Europe GmbH, Kurfürsten-Anlage 36, 69115 Heidelberg, Germany
| | - Federico Errica
- NEC Laboratories Europe GmbH, Kurfürsten-Anlage 36, 69115 Heidelberg, Germany
| | - Francesco Alesiani
- NEC Laboratories Europe GmbH, Kurfürsten-Anlage 36, 69115 Heidelberg, Germany
| |
Collapse
|
2
|
Barhaghi MS, Crawford B, Schwing G, Hardy DJ, Stone JE, Schwiebert L, Potoff J, Tajkhorshid E. py-MCMD: Python Software for Performing Hybrid Monte Carlo/Molecular Dynamics Simulations with GOMC and NAMD. J Chem Theory Comput 2022; 18:4983-4994. [PMID: 35621307 PMCID: PMC9760104 DOI: 10.1021/acs.jctc.1c00911] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
py-MCMD, an open-source Python software, provides a robust workflow layer that manages communication of relevant system information between the simulation engines NAMD and GOMC and generates coherent thermodynamic properties and trajectories for analysis. To validate the workflow and highlight its capabilities, hybrid Monte Carlo/molecular dynamics (MC/MD) simulations are performed for SPC/E water in the isobaric-isothermal (NPT) and grand canonical (GC) ensembles as well as with Gibbs ensemble Monte Carlo (GEMC). The hybrid MC/MD approach shows close agreement with reference MC simulations and has a computational efficiency that is 2 to 136 times greater than traditional Monte Carlo simulations. MC/MD simulations performed for water in a graphene slit pore illustrate significant gains in sampling efficiency when the coupled-decoupled configurational-bias MC (CD-CBMC) algorithm is used compared with simulations using a single unbiased random trial position. Simulations using CD-CBMC reach equilibrium with 25 times fewer cycles than simulations using a single unbiased random trial position, with a small increase in computational cost. In a more challenging application, hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to hydrate a buried binding pocket in bovine pancreatic trypsin inhibitor. Water occupancies produced by GCMC/MD simulations are in close agreement with crystallographically identified positions, and GCMC/MD simulations have a computational efficiency that is 5 times better than MD simulations. py-MCMD is available on GitHub at https://github.com/GOMC-WSU/py-MCMD.
Collapse
Affiliation(s)
- Mohammad Soroush Barhaghi
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Brad Crawford
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, United States
| | - Gregory Schwing
- Department of Computer Science, Wayne State University, Detroit, Michigan 48202, United States
| | - David J Hardy
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - John E Stone
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Loren Schwiebert
- Department of Computer Science, Wayne State University, Detroit, Michigan 48202, United States
| | - Jeffrey Potoff
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, United States
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
3
|
Schoeman D, Cloete R, Fielding BC. The Flexible, Extended Coil of the PDZ-Binding Motif of the Three Deadly Human Coronavirus E Proteins Plays a Role in Pathogenicity. Viruses 2022; 14:v14081707. [PMID: 36016329 PMCID: PMC9416557 DOI: 10.3390/v14081707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Revised: 07/22/2022] [Accepted: 07/29/2022] [Indexed: 02/04/2023] Open
Abstract
The less virulent human (h) coronaviruses (CoVs) 229E, NL63, OC43, and HKU1 cause mild, self-limiting respiratory tract infections, while the more virulent SARS-CoV-1, MERS-CoV, and SARS-CoV-2 have caused severe outbreaks. The CoV envelope (E) protein, an important contributor to the pathogenesis of severe hCoV infections, may provide insight into this disparate severity of the disease. We, therefore, generated full-length E protein models for SARS-CoV-1 and -2, MERS-CoV, HCoV-229E, and HCoV-NL63 and docked C-terminal peptides of each model to the PDZ domain of the human PALS1 protein. The PDZ-binding motif (PBM) of the SARS-CoV-1 and -2 and MERS-CoV models adopted a more flexible, extended coil, while the HCoV-229E and HCoV-NL63 models adopted a less flexible alpha helix. All the E peptides docked to PALS1 occupied the same binding site and the more virulent hCoV E peptides generally interacted more stably with PALS1 than the less virulent ones. We hypothesize that the increased flexibility of the PBM in the more virulent hCoVs facilitates more stable binding to various host proteins, thereby contributing to more severe disease. This is the first paper to model full-length 3D structures for both the more virulent and less virulent hCoV E proteins, providing novel insights for possible drug and/or vaccine development.
Collapse
Affiliation(s)
- Dewald Schoeman
- Molecular Biology and Virology Research Laboratory, Department of Medical Biosciences, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
| | - Ruben Cloete
- South African Medical Research Council Bioinformatics Unit, South African National Bioinformatics Institute, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
| | - Burtram C. Fielding
- Molecular Biology and Virology Research Laboratory, Department of Medical Biosciences, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
- Correspondence:
| |
Collapse
|
4
|
Basciu A, Callea L, Motta S, Bonvin AM, Bonati L, Vargiu AV. No dance, no partner! A tale of receptor flexibility in docking and virtual screening. VIRTUAL SCREENING AND DRUG DOCKING 2022. [DOI: 10.1016/bs.armc.2022.08.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
Peter EK, Manstein DJ, Shea JE, Schug A. CORE-MD II: A fast, adaptive, and accurate enhanced sampling method. J Chem Phys 2021; 155:104114. [PMID: 34525829 DOI: 10.1063/5.0063664] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we present a fast and adaptive correlation guided enhanced sampling method (CORE-MD II). The CORE-MD II technique relies, in part, on partitioning of the entire pathway into short trajectories that we refer to as instances. The sampling within each instance is accelerated by adaptive path-dependent metadynamics simulations. The second part of this approach involves kinetic Monte Carlo (kMC) sampling between the different states that have been accessed during each instance. Through the combination of the partition of the total simulation into short non-equilibrium simulations and the kMC sampling, the CORE-MD II method is capable of sampling protein folding without any a priori definitions of reaction pathways and additional parameters. In the validation simulations, we applied the CORE-MD II on the dialanine peptide and the folding of two peptides: TrpCage and TrpZip2. In a comparison with long time equilibrium Molecular Dynamics (MD), 1 µs replica exchange MD (REMD), and CORE-MD I simulations, we find that the level of convergence of the CORE-MD II method is improved by a factor of 8.8, while the CORE-MD II method reaches acceleration factors of ∼120. In the CORE-MD II simulation of TrpZip2, we observe the formation of the native state in contrast to the REMD and the CORE-MD I simulations. The method is broadly applicable for MD simulations and is not restricted to simulations of protein folding or even biomolecules but also applicable to simulations of protein aggregation, protein signaling, or even materials science simulations.
Collapse
Affiliation(s)
- Emanuel K Peter
- Institute for Biophysical Chemistry, Fritz-Hartmann-Centre for Medical Research, Hannover Medical School, Carl-Neuberg-Str. 1, Hannover 30625, Germany
| | - Dietmar J Manstein
- Institute for Biophysical Chemistry, Fritz-Hartmann-Centre for Medical Research, Hannover Medical School, Carl-Neuberg-Str. 1, Hannover 30625, Germany
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, Department of Physics, University of California, Santa Barbara, California 93106, USA
| | - Alexander Schug
- John von Neumann Institute for Computing and Jülich Supercomputing Centre, Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
6
|
Goel H, Yu W, Ustach VD, Aytenfisu AH, Sun D, MacKerell AD. Impact of electronic polarizability on protein-functional group interactions. Phys Chem Chem Phys 2020; 22:6848-6860. [PMID: 32195493 PMCID: PMC7194236 DOI: 10.1039/d0cp00088d] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Interactions of proteins with functional groups are key to their biological functions, making it essential that they be accurately modeled. To investigate the impact of the inclusion of explicit treatment of electronic polarizability in force fields on protein-functional group interactions, the additive CHARMM and Drude polarizable force field are compared in the context of the Site-Identification by Ligand Competitive Saturation (SILCS) simulation methodology from which functional group interaction patterns with five proteins for which experimental binding affinities of multiple ligands are available, were obtained. The explicit treatment of polarizability produces significant differences in the functional group interactions in the ligand binding sites including overall enhanced binding of functional groups to the proteins. This is associated with variations of the dipole moments of solutes representative of functional groups in the binding sites relative to aqueous solution with higher dipole moments systematically occurring in the latter, though exceptions occur with positively charged methylammonium. Such variation indicates the complex, heterogeneous nature of the electronic environments of ligand binding sites and emphasizes the inherent limitation of fixed charged, additive force fields for modeling ligand-protein interactions. These effects yield more defined orientation of the functional groups in the binding pockets and a small, but systematic improvement in the ability of the SILCS method to predict the binding orientation and relative affinities of ligands to their target proteins. Overall, these results indicate that the physical model associated with the explicit treatment of polarizability along with the presence of lone pairs in a force field leads to changes in the nature of the interactions of functional groups with proteins versus that occurring with additive force fields, suggesting the utility of polarizable force fields in obtaining a more realistic understanding of protein-ligand interactions.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| | - Vincent D Ustach
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| | - Asaminew H Aytenfisu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| | - Delin Sun
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| | - Alexander D MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St., Baltimore, Maryland 21201, USA.
| |
Collapse
|
7
|
Okamoto Y. Protein structure predictions by enhanced conformational sampling methods. Biophys Physicobiol 2019; 16:344-366. [PMID: 31984190 PMCID: PMC6976031 DOI: 10.2142/biophysico.16.0_344] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Accepted: 08/07/2019] [Indexed: 12/01/2022] Open
Abstract
In this Special Festschrift Issue for the celebration of Professor Nobuhiro Gō's 80th birthday, we review enhanced conformational sampling methods for protein structure predictions. We present several generalized-ensemble algorithms such as multicanonical algorithm, replica-exchange method, etc. and parallel Monte Carlo or molecular dynamics method with genetic crossover. Examples of the results of these methods applied to the predictions of protein tertiary structures are also presented.
Collapse
Affiliation(s)
- Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
- Structural Biology Research Center, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Aichi 464-8603, Japan
- Information Technology Center, Nagoya University, Nagoya, Aichi 464-8601, Japan
- JST-CREST, Nagoya, Aichi 464-8602, Japan
| |
Collapse
|
8
|
Das S, Shimshi M, Raz K, Nitoker Eliaz N, Mhashal AR, Ansbacher T, Major DT. EnzyDock: Protein–Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes. J Chem Theory Comput 2019; 15:5116-5134. [DOI: 10.1021/acs.jctc.9b00366] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Susanta Das
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Mor Shimshi
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Keren Raz
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | | | - Anil Ranu Mhashal
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Tamar Ansbacher
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
- Hadassah Academic College, 7 Hanevi’im Street, Jerusalem 9101001, Israel
| | - Dan T. Major
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
9
|
Yang KC, Cui FC, Shi C, Chen WD, Li YQ. Evolution of Conformation and Dynamics of Solvents in Hydration Shell along the Urea-induced Unfolding of Ubiquitin. CHINESE JOURNAL OF POLYMER SCIENCE 2019. [DOI: 10.1007/s10118-019-2238-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
|
10
|
Huang WT, Chen PS, Liu JJ, Chen YR, Chen YH. Dynamic configuration scheduling problem for stochastic medical resources. J Biomed Inform 2018; 80:96-105. [DOI: 10.1016/j.jbi.2018.03.005] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2017] [Revised: 12/27/2017] [Accepted: 03/12/2018] [Indexed: 11/28/2022]
|
11
|
In silico methods for design of biological therapeutics. Methods 2017; 131:33-65. [PMID: 28958951 DOI: 10.1016/j.ymeth.2017.09.008] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Revised: 09/21/2017] [Accepted: 09/23/2017] [Indexed: 12/18/2022] Open
Abstract
It has been twenty years since the first rationally designed small molecule drug was introduced into the market. Since then, we have progressed from designing small molecules to designing biotherapeutics. This class of therapeutics includes designed proteins, peptides and nucleic acids that could more effectively combat drug resistance and even act in cases where the disease is caused because of a molecular deficiency. Computational methods are crucial in this design exercise and this review discusses the various elements of designing biotherapeutic proteins and peptides. Many of the techniques discussed here, such as the deterministic and stochastic design methods, are generally used in protein design. We have devoted special attention to the design of antibodies and vaccines. In addition to the methods for designing these molecules, we have included a comprehensive list of all biotherapeutics approved for clinical use. Also included is an overview of methods that predict the binding affinity, cell penetration ability, half-life, solubility, immunogenicity and toxicity of the designed therapeutics. Biotherapeutics are only going to grow in clinical importance and are set to herald a new generation of disease management and cure.
Collapse
|
12
|
Urano R, Okamoto Y. Observation of helix associations for insertion of a retinal molecule and distortions of helix structures in bacteriorhodopsin. J Chem Phys 2016; 143:235101. [PMID: 26696075 DOI: 10.1063/1.4935964] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We applied a newly proposed prediction method for membrane protein structures to bacteriorhodopsin that has distorted transmembrane helices in the native structure. This method uses an implicit membrane model, which restricts sampling space during folding in a membrane region, and includes helix bending. Replica-exchange simulations were performed with seven transmembrane helices only without a retinal molecule. Obtained structures were classified into clusters of similar structures, which correspond to local-minimum free energy states. The two lowest free energy states corresponded to a native-like structure with the correct empty space for retinal and a structure with this empty space filled with a helix. Previous experiments of bacteriorhodopsin suggested that association of transmembrane helices enables them to make a room for insertion of a retinal. Our results are consistent with these results. Moreover, distortions of helices in the native-like structures were successfully reproduced. In the distortions, whereas the locations of kinks for all helices were similar to those of Protein Data Bank's data, the amount of bends was more similar for helices away from the retinal than for those close to the retinal in the native structure. This suggests a hypothesis that the amino-acid sequence specifies the location of kinks in transmembrane helices and that the amount of distortions depends on the interactions with the surrounding molecules such as neighboring helices, lipids, and retinal.
Collapse
Affiliation(s)
- Ryo Urano
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
| | - Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
| |
Collapse
|
13
|
Yahya R, Karjiban RA, Basri M, Rahman MBA, Girardi M. Monte Carlo simulation of mixed nonionic Brij surfactants in water. J Mol Model 2014; 20:2512. [PMID: 25381172 DOI: 10.1007/s00894-014-2512-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2014] [Accepted: 10/21/2014] [Indexed: 11/24/2022]
Abstract
Nonionic surfactants such as the Brij® series are important in the preparation of transdermal drug nanodelivery products using nanoemulsions because of their low toxicity and low irritancy. Here, Monte Carlo (MC) simulation was used to examine the physical behavior of the model deterministic system by using sampling procedures. Metropolis MC simulations were run on three mixtures of two different nonionic surfactants, Brij92 and Brij96, with different compositions in aqueous solution. The system was simulated in the canonical ensemble with constant temperature, volume and number of molecules. Hence, the acceptance ratio for single atom moves of the mixed surfactants increased as the concentration of surfactants increased from 0.494 to 0.591. The lowest total energy for the mixed surfactant systems was -99,039 kcal mol(-1) due to the interaction between all molecules in the system simulated. The physicochemical properties of models such as the radius of gyration and radial distribution function, were also determined. These observations indicate that the behavior and physicochemical of mixed surfactant and PKOEs nanoemulsion systems were described adequately during the simulation.
Collapse
Affiliation(s)
- Ruzanna Yahya
- Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, Selangor, Malaysia
| | | | | | | | | |
Collapse
|
14
|
Hédin F, Plattner N, Doll JD, Meuwly M. Spatial Averaging: Sampling Enhancement for Exploring Configurational Space of Atomic Clusters and Biomolecules. J Chem Theory Comput 2014; 10:4284-96. [DOI: 10.1021/ct500529w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Florent Hédin
- Department
of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Nuria Plattner
- Department
of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
| | - J. D. Doll
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Markus Meuwly
- Department
of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
15
|
Fang Y, Sanz-Serna JM, Skeel RD. Compressible generalized hybrid Monte Carlo. J Chem Phys 2014; 140:174108. [DOI: 10.1063/1.4874000] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Wüstner D, Sklenar H. Atomistic Monte Carlo simulation of lipid membranes. Int J Mol Sci 2014; 15:1767-803. [PMID: 24469314 PMCID: PMC3958820 DOI: 10.3390/ijms15021767] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2013] [Revised: 12/06/2013] [Accepted: 01/09/2014] [Indexed: 02/07/2023] Open
Abstract
Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.
Collapse
Affiliation(s)
- Daniel Wüstner
- Department of Biochemistry and Molecular Biology, University of Southern Denmark, Odense M DK-5230, Denmark.
| | - Heinz Sklenar
- Theoretical Biophysics Group, Max Delbrück Center for Molecular Medicine, Robert-Rössle-Str. 10, Berlin D-13125, Germany.
| |
Collapse
|
17
|
Small MC, Lopes P, Andrade RB, MacKerell AD. Impact of ribosomal modification on the binding of the antibiotic telithromycin using a combined grand canonical monte carlo/molecular dynamics simulation approach. PLoS Comput Biol 2013; 9:e1003113. [PMID: 23785274 PMCID: PMC3681621 DOI: 10.1371/journal.pcbi.1003113] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2013] [Accepted: 05/07/2013] [Indexed: 02/05/2023] Open
Abstract
Resistance to macrolide antibiotics is conferred by mutation of A2058 to G or methylation by Erm methyltransferases of the exocyclic N6 of A2058 (E. coli numbering) that forms the macrolide binding site in the 50S subunit of the ribosome. Ketolides such as telithromycin mitigate A2058G resistance yet remain susceptible to Erm-based resistance. Molecular details associated with macrolide resistance due to the A2058G mutation and methylation at N6 of A2058 by Erm methyltransferases were investigated using empirical force field-based simulations. To address the buried nature of the macrolide binding site, the number of waters within the pocket was allowed to fluctuate via the use of a Grand Canonical Monte Carlo (GCMC) methodology. The GCMC water insertion/deletion steps were alternated with Molecular Dynamics (MD) simulations to allow for relaxation of the entire system. From this GCMC/MD approach information on the interactions between telithromycin and the 50S ribosome was obtained. In the wild-type (WT) ribosome, the 2'-OH to A2058 N1 hydrogen bond samples short distances with a higher probability, while the effectiveness of telithromycin against the A2058G mutation is explained by a rearrangement of the hydrogen bonding pattern of the 2'-OH to 2058 that maintains the overall antibiotic-ribosome interactions. In both the WT and A2058G mutation there is significant flexibility in telithromycin's imidazole-pyridine side chain (ARM), indicating that entropic effects contribute to the binding affinity. Methylated ribosomes show lower sampling of short 2'-OH to 2058 distances and also demonstrate enhanced G2057-A2058 stacking leading to disrupted A752-U2609 Watson-Crick (WC) interactions as well as hydrogen bonding between telithromycin's ARM and U2609. This information will be of utility in the rational design of novel macrolide analogs with improved activity against methylated A2058 ribosomes.
Collapse
Affiliation(s)
- Meagan C. Small
- Department of Pharmaceutical Sciences, University of Maryland, Baltimore, Maryland, United States of America
| | - Pedro Lopes
- Department of Pharmaceutical Sciences, University of Maryland, Baltimore, Maryland, United States of America
| | - Rodrigo B. Andrade
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania, United States of America
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, University of Maryland, Baltimore, Maryland, United States of America
| |
Collapse
|
18
|
Boomsma W, Frellsen J, Harder T, Bottaro S, Johansson KE, Tian P, Stovgaard K, Andreetta C, Olsson S, Valentin JB, Antonov LD, Christensen AS, Borg M, Jensen JH, Lindorff-Larsen K, Ferkinghoff-Borg J, Hamelryck T. PHAISTOS: a framework for Markov chain Monte Carlo simulation and inference of protein structure. J Comput Chem 2013; 34:1697-705. [PMID: 23619610 DOI: 10.1002/jcc.23292] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2012] [Revised: 03/14/2013] [Accepted: 03/20/2013] [Indexed: 11/10/2022]
Abstract
We present a new software framework for Markov chain Monte Carlo sampling for simulation, prediction, and inference of protein structure. The software package contains implementations of recent advances in Monte Carlo methodology, such as efficient local updates and sampling from probabilistic models of local protein structure. These models form a probabilistic alternative to the widely used fragment and rotamer libraries. Combined with an easily extendible software architecture, this makes PHAISTOS well suited for Bayesian inference of protein structure from sequence and/or experimental data. Currently, two force-fields are available within the framework: PROFASI and OPLS-AA/L, the latter including the generalized Born surface area solvent model. A flexible command-line and configuration-file interface allows users quickly to set up simulations with the desired configuration. PHAISTOS is released under the GNU General Public License v3.0. Source code and documentation are freely available from http://phaistos.sourceforge.net. The software is implemented in C++ and has been tested on Linux and OSX platforms.
Collapse
Affiliation(s)
- Wouter Boomsma
- Department of Biology, University of Copenhagen, Copenhagen, 2200, Denmark
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
19
|
|
20
|
Ribeiro AAST, de Alencastro RB. Mixed monte carlo/molecular dynamics simulations in explicit solvent. J Comput Chem 2012; 33:901-5. [DOI: 10.1002/jcc.22925] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2011] [Revised: 11/04/2011] [Accepted: 12/09/2011] [Indexed: 11/06/2022]
|
21
|
Ramya L, Gautham N. Conformational space exploration of met- and Leu-enkephalin using the mols method, molecular dynamics, and Monte Carlo simulation-a comparative study. Biopolymers 2011; 97:165-76. [DOI: 10.1002/bip.21721] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2011] [Revised: 09/16/2011] [Accepted: 09/16/2011] [Indexed: 11/09/2022]
|
22
|
Abstract
Equilibrium sampling of biomolecules remains an unmet challenge after more than 30 years of atomistic simulation. Efforts to enhance sampling capability, which are reviewed here, range from the development of new algorithms to parallelization to novel uses of hardware. Special focus is placed on classifying algorithms--most of which are underpinned by a few key ideas--in order to understand their fundamental strengths and limitations. Although algorithms have proliferated, progress resulting from novel hardware use appears to be more clear-cut than from algorithms alone, due partly to the lack of widely used sampling measures.
Collapse
Affiliation(s)
- Daniel M Zuckerman
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15213, USA.
| |
Collapse
|
23
|
|
24
|
Acevedo O, Jorgensen WL. Advances in quantum and molecular mechanical (QM/MM) simulations for organic and enzymatic reactions. Acc Chem Res 2010; 43:142-51. [PMID: 19728702 DOI: 10.1021/ar900171c] [Citation(s) in RCA: 176] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Application of combined quantum and molecular mechanical (QM/MM) methods focuses on predicting activation barriers and the structures of stationary points for organic and enzymatic reactions. Characterization of the factors that stabilize transition structures in solution and in enzyme active sites provides a basis for design and optimization of catalysts. Continued technological advances allowed for expansion from prototypical cases to mechanistic studies featuring detailed enzyme and condensed-phase environments with full integration of the QM calculations and configurational sampling. This required improved algorithms featuring fast QM methods, advances in computing changes in free energies including free-energy perturbation (FEP) calculations, and enhanced configurational sampling. In particular, the present Account highlights development of the PDDG/PM3 semi-empirical QM method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo (MC) simulations, and a polynomial quadrature method for efficient modeling of proton-transfer reactions. The utility of this QM/MM/MC/FEP methodology is illustrated for a variety of organic reactions including substitution, decarboxylation, elimination, and pericyclic reactions. A comparison to experimental kinetic results on medium effects has verified the accuracy of the QM/MM approach in the full range of solvents from hydrocarbons to water to ionic liquids. Corresponding results from ab initio and density functional theory (DFT) methods with continuum-based treatments of solvation reveal deficiencies, particularly for protic solvents. Also summarized in this Account are three specific QM/MM applications to biomolecular systems: (1) a recent study that clarified the mechanism for the reaction of 2-pyrone derivatives catalyzed by macrophomate synthase as a tandem Michael-aldol sequence rather than a Diels-Alder reaction, (2) elucidation of the mechanism of action of fatty acid amide hydrolase (FAAH), an unusual Ser-Ser-Lys proteolytic enzyme, and (3) the construction of enzymes for Kemp elimination of 5-nitrobenzisoxazole that highlights the utility of QM/MM in the design of artificial enzymes.
Collapse
Affiliation(s)
- Orlando Acevedo
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849
| | - William L. Jorgensen
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107
| |
Collapse
|
25
|
|
26
|
Abstract
Diffusion on a low-dimensional free-energy surface is a remarkably successful model for the folding dynamics of small single-domain proteins. Complicating the interpretation of both simulations and experiments is the expectation that the effective diffusion coefficient D will in general depend on the position along the folding coordinate, and this dependence may vary for different coordinates. Here we explore the position dependence of D, its connection to protein internal friction, and the consequences for the interpretation of single-molecule experiments. We find a large decrease in D from unfolded to folded, for reaction coordinates that directly measure fluctuations in Cartesian configuration space, including those probed in single-molecule experiments. In contrast, D is almost independent of Q, the fraction of native amino acid contacts: Near the folded state, small fluctuations in position cause large fluctuations in Q, and vice versa for the unfolded state. In general, position-dependent free energies and diffusion coefficients for any two good reaction coordinates that separate reactant, product, and transition states, are related by a simple transformation, as we demonstrate. With this transformation, we obtain reaction coordinates with position-invariant D. The corresponding free-energy surfaces allow us to justify the assumptions used in estimating the speed limit for protein folding from experimental measurements of the reconfiguration time in the unfolded state, and also reveal intermediates hidden in the original free-energy projection. Lastly, we comment on the design of future single-molecule experiments that probe the position dependence of D directly.
Collapse
|
27
|
Brooks B, Brooks C, MacKerell A, Nilsson L, Petrella R, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner A, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor R, Post C, Pu J, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York D, Karplus M. CHARMM: the biomolecular simulation program. J Comput Chem 2009; 30:1545-614. [PMID: 19444816 PMCID: PMC2810661 DOI: 10.1002/jcc.21287] [Citation(s) in RCA: 6011] [Impact Index Per Article: 400.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.
Collapse
Affiliation(s)
- B.R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and
Blood Institute, National Institutes of Health, Bethesda, MD 20892
| | - C.L. Brooks
- Departments of Chemistry & Biophysics, University of
Michigan, Ann Arbor, MI 48109
| | - A.D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy,
University of Maryland, Baltimore, MD, 21201
| | - L. Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition,
SE-141 57, Huddinge, Sweden
| | - R.J. Petrella
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Department of Medicine, Harvard Medical School, Boston, MA
02115
| | - B. Roux
- Department of Biochemistry and Molecular Biology, University of
Chicago, Gordon Center for Integrative Science, Chicago, IL 60637
| | - Y. Won
- Department of Chemistry, Hanyang University, Seoul
133–792 Korea
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - M. Karplus
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Laboratoire de Chimie Biophysique, ISIS, Université de
Strasbourg, 67000 Strasbourg France
| |
Collapse
|
28
|
Schlick T. Monte Carlo, harmonic approximation, and coarse-graining approaches for enhanced sampling of biomolecular structure. F1000 BIOLOGY REPORTS 2009; 1:48. [PMID: 20948637 PMCID: PMC2924683 DOI: 10.3410/b1-48] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The rugged energy landscape of biomolecules and associated large-scale conformational changes have triggered the development of many innovative enhanced sampling methods, either based or not based on molecular dynamics (MD) simulations. Surveyed here are methods in the latter class - including Monte Carlo methods, harmonic approximations, and coarse graining - many of which yield valuable conformational insights into biomolecular structure and flexibility, despite altered kinetics. MD-based methods are surveyed in an upcoming issue of F1000 Biology Reports.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University 251 Mercer Street, New York, NY 10012 USA.
| |
Collapse
|
29
|
Velez-Vega C, Borrero EE, Escobedo FA. Kinetics and reaction coordinate for the isomerization of alanine dipeptide by a forward flux sampling protocol. J Chem Phys 2009; 130:225101. [DOI: 10.1063/1.3147465] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
30
|
Automatic Determination of Stepsize Parameters in Monte Carlo Simulation Tested on a Bromodomain-Binding Octapeptide. ALGORITHMS 2009. [DOI: 10.3390/a2010215] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
31
|
Vitalis A, Pappu RV. Methods for Monte Carlo simulations of biomacromolecules. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2009; 5:49-76. [PMID: 20428473 PMCID: PMC2860296 DOI: 10.1016/s1574-1400(09)00503-9] [Citation(s) in RCA: 94] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The state-of-the-art for Monte Carlo (MC) simulations of biomacromolecules is reviewed. Available methodologies for sampling conformational equilibria and associations of biomacromolecules in the canonical ensemble, given a continuum description of the solvent environment, are reviewed. Detailed sections are provided dealing with the choice of degrees of freedom, the efficiencies of MC algorithms and algorithmic peculiarities, as well as the optimization of simple movesets. The issue of introducing correlations into elementary MC moves, and the applicability of such methods to simulations of biomacromolecules is discussed. A brief discussion of multicanonical methods and an overview of recent simulation work highlighting the potential of MC methods are also provided. It is argued that MC simulations, while underutilized biomacromolecular simulation community, hold promise for simulations of complex systems and phenomena that span multiple length scales, especially when used in conjunction with implicit solvation models or other coarse graining strategies.
Collapse
Affiliation(s)
- Andreas Vitalis
- Department of Biomedical Engineering, Molecular Biophysics Program, Center for Computational Biology, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, MO 63130-4899, USA
| | | |
Collapse
|
32
|
Casey FP, Waterfall JJ, Gutenkunst RN, Myers CR, Sethna JP. Variational method for estimating the rate of convergence of Markov-chain Monte Carlo algorithms. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:046704. [PMID: 18999558 DOI: 10.1103/physreve.78.046704] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2008] [Indexed: 05/27/2023]
Abstract
We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces. We apply the method to one dimensional examples with Gaussian and quartic target densities, and we contrast the performance of the random walk Metropolis-Hastings algorithm with a "smart" variant that incorporates gradient information into the trial moves, a generalization of the Metropolis adjusted Langevin algorithm. We find that the variational method agrees quite closely with numerical simulations. We also see that the smart MCMC algorithm often fails to converge geometrically in the tails of the target density except in the simplest case we examine, and even then care must be taken to choose the appropriate scaling of the deterministic and random parts of the proposed moves. Again, this calls into question the utility of smart MCMC in more complex problems. Finally, we apply the same method to approximate the rate of convergence in multidimensional Gaussian problems with and without importance sampling. There we demonstrate the necessity of importance sampling for target densities which depend on variables with a wide range of scales.
Collapse
Affiliation(s)
- Fergal P Casey
- Complex and Adaptive Systems Laboratory, University College Dublin, Dublin 4, Ireland.
| | | | | | | | | |
Collapse
|
33
|
Predicting the folding pathway of engrailed homeodomain with a probabilistic roadmap enhanced reaction-path algorithm. Biophys J 2008; 94:1622-9. [PMID: 18024496 DOI: 10.1529/biophysj.107.119214] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
To predict a protein-folding pathway, we present an alternative to the time-consuming dynamic simulation of atomistic models. We replace the actual dynamic simulation with variational optimization of a reaction path connecting known initial and final protein conformations in such a way as to maximize an estimate of the reactive flux or minimize the mean first passage time at a given temperature, referred to as MaxFlux. We solve the MaxFlux global optimization problem with an efficient graph-theoretic approach, the probabilistic roadmap method (PRM). We employed CHARMM19 and the EEF1 implicit solvation model to describe the protein solution. The effectiveness of our MaxFlux-PRM is demonstrated in our promising simulation results on the folding pathway of the engrailed homeodomain. Our MaxFlux-PRM approach provides the direct evidence to support that the previously reported intermediate state is a genuine on-pathway intermediate, and the demand of CPU power is moderate.
Collapse
|
34
|
Lund M, Trulsson M, Persson B. Faunus: An object oriented framework for molecular simulation. SOURCE CODE FOR BIOLOGY AND MEDICINE 2008; 3:1. [PMID: 18241331 PMCID: PMC2266748 DOI: 10.1186/1751-0473-3-1] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/18/2007] [Accepted: 02/01/2008] [Indexed: 11/29/2022]
Abstract
Background We present a C++ class library for Monte Carlo simulation of molecular systems, including proteins in solution. The design is generic and highly modular, enabling multiple developers to easily implement additional features. The statistical mechanical methods are documented by extensive use of code comments that – subsequently – are collected to automatically build a web-based manual. Results We show how an object oriented design can be used to create an intuitively appealing coding framework for molecular simulation. This is exemplified in a minimalistic C++ program that can calculate protein protonation states. We further discuss performance issues related to high level coding abstraction. Conclusion C++ and the Standard Template Library (STL) provide a high-performance platform for generic molecular modeling. Automatic generation of code documentation from inline comments has proven particularly useful in that no separate manual needs to be maintained.
Collapse
Affiliation(s)
- Mikael Lund
- Institute of Organic Chemistry and Biochemistry, The Academy of Sciences of the Czech Republic, Flemingovo nam,2, CZ-16610 Prague 6, Czech Republic.
| | | | | |
Collapse
|
35
|
Abstract
A description of Monte Carlo methods for simulation of proteins is given. Advantages and disadvantages of the Monte Carlo approach are presented. The theoretical basis for calculating equilibrium properties of biological molecules by the Monte Carlo method is presented. Some of the standard and some of the more recent ways of performing Monte Carlo on proteins are presented. A discussion of the estimation of errors in properties calculated by Monte Carlo is given.
Collapse
Affiliation(s)
- David J Earl
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA, USA
| | | |
Collapse
|
36
|
Burghardt TP, Hu JY, Ajtai K. Myosin dynamics on the millisecond time scale. Biophys Chem 2007; 131:15-28. [PMID: 17913331 PMCID: PMC2505346 DOI: 10.1016/j.bpc.2007.08.008] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2007] [Revised: 08/27/2007] [Accepted: 08/27/2007] [Indexed: 11/20/2022]
Abstract
Myosin is a motor protein associating with actin and ATP. It translates along actin filaments against a force by transduction of free energy liberated with ATP hydrolysis. Various myosin crystal structures define time points during ATPase showing the protein undergoes large conformation change during transduction over a cycle with approximately 10 ms periodicity. The protein conformation trajectory between two intermediates in the cycle is surmised by non-equilibrium Monte Carlo simulation utilizing free-energy minimization. The trajectory shows myosin transduction of free energy to mechanical work giving evidence for: (i) a causal relationship between product release and work production in the native isoform that is correctly disrupted in a chemically modified protein, (ii) the molecular basis of ATP-sensitive tryptophan fluorescence enhancement and acrylamide quenching, (iii) an actin-binding site peptide containing the free-energy barrier to ATPase product release defining the rate limiting step and, (iv) a scenario for actin-activation of myosin ATPase.
Collapse
Affiliation(s)
- Thomas P Burghardt
- Department of Biochemistry and Molecular Biology, Mayo Clinic College of Medicine, Rochester, MN 55905, USA.
| | | | | |
Collapse
|
37
|
Hu J, Ma A, Dinner AR. Bias annealing: a method for obtaining transition paths de novo. J Chem Phys 2007; 125:114101. [PMID: 16999460 DOI: 10.1063/1.2335640] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Computational studies of dynamics in complex systems require means for generating reactive trajectories with minimum knowledge about the processes of interest. Here, we introduce a method for generating transition paths when an existing one is not already available. Starting from biased paths obtained from steered molecular dynamics, we use a Monte Carlo procedure in the space of whole trajectories to shift gradually to sampling an ensemble of unbiased paths. Application to basin-to-basin hopping in a two-dimensional model system and nucleotide-flipping by a DNA repair protein demonstrates that the method can efficiently yield unbiased reactive trajectories even when the initial steered dynamics differ significantly. The relation of the method to others and the physical basis for its success are discussed.
Collapse
Affiliation(s)
- Jie Hu
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | | | | |
Collapse
|
38
|
Ma A, Nag A, Dinner AR. Dynamic coupling between coordinates in a model for biomolecular isomerization. J Chem Phys 2006; 124:144911. [PMID: 16626249 DOI: 10.1063/1.2183768] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
To understand a complex reaction, it is necessary to project the dynamics of the system onto a low-dimensional subspace of physically meaningful coordinates. We recently introduced an automatic method for identifying coordinates that relate closely to stable-state commitment probabilities and successfully applied it to a model for biomolecular isomerization, the C(7eq)-->alpha(R) transition of the alanine dipeptide [A. Ma and A. R. Dinner, J. Phys. Chem. B 109, 6769 (2005)]. Here, we explore approximate means for estimating diffusion tensors for systems subject to restraints in one and two dimensions and then use the results together with an extension of Kramers theory for unimolecular reaction rates [A. Berezhkovskii and A. Szabo, J. Chem. Phys. 122, 014503 (2005)] to show explicitly that both the potential of mean force and the diffusion tensor are essential for describing the dynamics of the alanine dipeptide quantitatively. In particular, the signficance of off-diagonal elements of the diffusion tensor suggests that the coordinates of interest are coupled by the hydrodynamic-like response of the bath of remaining degrees of freedom.
Collapse
Affiliation(s)
- Ao Ma
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5640 S. Ellis Avenue, Chicago, Illinois 60637, USA
| | | | | |
Collapse
|
39
|
Cheluvaraja S, Meirovitch H. Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method. J Chem Phys 2006; 125:24905. [PMID: 16848609 DOI: 10.1063/1.2208608] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Hypothetical scanning (HS) is a method for calculating the absolute entropy S and free energy F from a sample generated by any simulation technique. With this approach each sample configuration is reconstructed with the help of transition probabilities (TPs) and their product leads to the configuration's probability, hence to the entropy. Recently a new way for calculating the TPs by Monte Carlo (MC) simulations has been suggested, where all system interactions are taken into account. Therefore, this method--called HSMC--is in principle exact where the only approximation is due to insufficient sampling. HSMC has been applied very successfully to liquid argon, TIP3P water, self-avoiding walks on a lattice, and peptides. Because molecular dynamics (MD) is considered to be significantly more efficient than MC for a compact polymer chain, in this paper HSMC is extended to MD simulations as applied to peptides. Like before, we study decaglycine in vacuum but for the first time also a peptide with side chains, (Val)(2)(Gly)(6)(Val)(2). The transition from MC to MD requires implementing essential changes in the reconstruction process of HSMD. Results are calculated for three microstates, helix, extended, and hairpin. HSMD leads to very stable differences in entropy TDeltaS between these microstates with small errors of 0.1-0.2 kcal/mol (T=100 K) for a wide range of calculation parameters with extremely high efficiency. Various aspects of HSMD and plans for future work are discussed.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260, USA
| | | |
Collapse
|
40
|
Bitetti-Putzer R, Dinner AR, Yang W, Karplus M. Conformational sampling via a self-regulating effective energy surface. J Chem Phys 2006; 124:174901. [PMID: 16689598 DOI: 10.1063/1.2171194] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The difficulty of efficiently sampling the phase space of complex systems with rough energy surfaces is well known. Typical solutions to the problem involve accelerating the crossing of barriers, but such methods often have the secondary problem that the low-energy states of interest are inadequately sampled, unless the parameters of the search algorithm are modified as the system evolves. A method is presented to improve the sampling with particular emphasis on the low-energy conformations, which make the most important contributions to the thermodynamics of the system. The algorithm proposed here samples the details of the minima, while easily surmounting barriers. This is achieved by introducing a self-regulating sampling variable which depends on the current state of the system. Two replicas of the system are introduced and the sampling variable is treated as a particle coupled to the physical system. The method is illustrated with a simple model system and is applied to the realistic example of barrier crossing in a protein-ligand complex.
Collapse
Affiliation(s)
- Ryan Bitetti-Putzer
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | | | | | | |
Collapse
|