On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations

Improving the replica-exchange molecular-dynamics method for efficient sampling in the temperature space

Replica-exchange molecular dynamics (REMD) is a popular sampling method in the molecular simulation. By frequently exchanging the replicas at different temperatures, the molecule can jump out of the minima and sample efficiently in the conformational space. Although REMD has been shown to be practical in a lot of applications, it does have a critical limitation. All the replicas at all the temperatures must be simulated for a period between the replica-exchange steps. This may be problematic for the reaction with high free energy barriers. In that case, too many replicas are required in the simulation. To reduce the calculation quantity and improve its performance, in this paper we propose a modified REMD method. During the simulation, each replica at each temperature can stay in either the active or inactive state and only switch between the states at the exchange step. In the active state, the replica moves freely in the canonical ensemble by the normal molecular dynamics, and in the inactive state, the replica is frozen temporarily until the next exchange step. The number of the replicas in the active states (active replicas) depends on the number of CPUs in the computer. Using the additional inactive replicas, one can perform an REMD simulation in a wider temperature space. The practical applications show that the modified REMD method is reliable. With the same number of active replicas, this REMD method can produce a more reasonable free energy surface around the free energy minima than the standard REMD method.

What causes tumbling of altro-α-CD derivatives? Insight from computer simulations

Spontaneous tumbling of altro-α-CD derivatives stems from the solvent and the side chain. Simulation results provide a theoretical basis for design of novel rotaxane-based molecular reels.

Trypsinogen activation as observed in accelerated molecular dynamics simulations

Serine proteases are involved in many fundamental physiological processes, and control of their activity mainly results from the fact that they are synthetized in an inactive form that becomes active upon cleavage. Three decades ago Martin Karplus's group performed the first molecular dynamics simulations of trypsin, the most studied member of the serine protease family, to address the transition from the zymogen to its active form. Based on the computational power available at the time, only high frequency fluctuations, but not the transition steps, could be observed. By performing accelerated molecular dynamics (aMD) simulations, an interesting approach that increases the configurational sampling of atomistic simulations, we were able to observe the N-terminal tail insertion, a crucial step of the transition mechanism. Our results also support the hypothesis that the hydrophobic effect is the main force guiding the insertion step, although substantial enthalpic contributions are important in the activation mechanism. As the N-terminal tail insertion is a conserved step in the activation of serine proteases, these results afford new perspective on the underlying thermodynamics of the transition from the zymogen to the active enzyme.

Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics

BACKGROUND

Accelerated molecular dynamics (aMD) has been proven to be a powerful biasing method for enhanced sampling of biomolecular conformations on general-purpose computational platforms. Biologically important long timescale events that are beyond the reach of standard molecular dynamics can be accessed without losing the detailed atomistic description of the system in aMD. Over other biasing methods, aMD offers the advantages of tuning the level of acceleration to access the desired timescale without any advance knowledge of the reaction coordinate.

SCOPE OF REVIEW

Recent advances in the implementation of aMD and its applications to small peptides and biological macromolecules are reviewed here along with a brief account of all the aMD variants introduced in the last decade.

MAJOR CONCLUSIONS

In comparison to the original implementation of aMD, the recent variant in which all the rotatable dihedral angles are accelerated (RaMD) exhibits faster convergence rates and significant improvement in statistical accuracy of retrieved thermodynamic properties. RaMD in conjunction with accelerating diffusive degrees of freedom, i.e. dual boosting, has been rigorously tested for the most difficult conformational sampling problem, protein folding. It has been shown that RaMD with dual boosting is capable of efficiently sampling multiple folding and unfolding events in small fast folding proteins.

GENERAL SIGNIFICANCE

RaMD with the dual boost approach opens exciting possibilities for sampling multiple timescales in biomolecules. While equilibrium properties can be recovered satisfactorily from aMD-based methods, directly obtaining dynamics and kinetic rates for larger systems presents a future challenge. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Adaptive accelerated ReaxFF reactive dynamics with validation from simulating hydrogen combustion

We develop here the methodology for dramatically accelerating the ReaxFF reactive force field based reactive molecular dynamics (RMD) simulations through use of the bond boost concept (BB), which we validate here for describing hydrogen combustion. The bond order, undercoordination, and overcoordination concepts of ReaxFF ensure that the BB correctly adapts to the instantaneous configurations in the reactive system to automatically identify the reactions appropriate to receive the bond boost. We refer to this as adaptive Accelerated ReaxFF Reactive Dynamics or aARRDyn. To validate the aARRDyn methodology, we determined the detailed sequence of reactions for hydrogen combustion with and without the BB. We validate that the kinetics and reaction mechanisms (that is the detailed sequences of reactive intermediates and their subsequent transformation to others) for H2 oxidation obtained from aARRDyn agrees well with the brute force reactive molecular dynamics (BF-RMD) at 2498 K. Using aARRDyn, we then extend our simulations to the whole range of combustion temperatures from ignition (798 K) to flame temperature (2998K), and demonstrate that, over this full temperature range, the reaction rates predicted by aARRDyn agree well with the BF-RMD values, extrapolated to lower temperatures. For the aARRDyn simulation at 798 K we find that the time period for half the H2 to form H2O product is ∼538 s, whereas the computational cost was just 1289 ps, a speed increase of ∼0.42 trillion (10(12)) over BF-RMD. In carrying out these RMD simulations we found that the ReaxFF-COH2008 version of the ReaxFF force field was not accurate for such intermediates as H3O. Consequently we reoptimized the fit to a quantum mechanics (QM) level, leading to the ReaxFF-OH2014 force field that was used in the simulations.

Accelerated adaptive integration method

Conformational changes that occur upon ligand binding may be too slow to observe on the time scales routinely accessible using molecular dynamics simulations. The adaptive integration method (AIM) leverages the notion that when a ligand is either fully coupled or decoupled, according to λ, barrier heights may change, making some conformational transitions more accessible at certain λ values. AIM adaptively changes the value of λ in a single simulation so that conformations sampled at one value of λ seed the conformational space sampled at another λ value. Adapting the value of λ throughout a simulation, however, does not resolve issues in sampling when barriers remain high regardless of the λ value. In this work, we introduce a new method, called Accelerated AIM (AcclAIM), in which the potential energy function is flattened at intermediate values of λ, promoting the exploration of conformational space as the ligand is decoupled from its receptor. We show, with both a simple model system (Bromocyclohexane) and the more complex biomolecule Thrombin, that AcclAIM is a promising approach to overcome high barriers in the calculation of free energies, without the need for any statistical reweighting or additional processors.

MOLECULAR SIMULATION OF α-TOCOPHEROL PASSING ACROSS DPPC LIPID USING POTENTIAL OF MEAN FORCE AND ACCELERATED MOLECULAR DYNAMICS METHOD

In this paper the process of α-tocopherol (TCP) passing across DPPC membrane was simulated using both the potential of mean force (PMF) and the accelerated molecular dynamics (aMD) methods, respectively. Energy properties, hydrogen bonds and orientation have been compared between these two methods and several conclusions have been obtained. The results indicate that TCP tends to stay at z = 1.2 nm of lipid bilayer. The binding free energy profiles of these two methods are alike. All these show that aMD could obtain comparable results as PMF method, while needs less computation time and resources. Therefore, aMD method could be used as an alternative method for prediction of transport properties of drug-lipid system.

Accelerated Molecular Dynamics Simulations with the AMOEBA Polarizable Force Field on Graphics Processing Units

The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling.

Substrate-dependent dynamics of UDP-galactopyranose mutase: Implications for drug design

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected tropical disease that represents one of the major health challenges of the Latin American countries. Successful efforts were made during the last few decades to control the transmission of this disease, but there is still no treatment for the 10 million adults in the chronic phase of the disease. In T. cruzi, as well as in other pathogens, the flavoenzyme UDP-galactopyranose mutase (UGM) catalyzes the conversion of UDP-galactopyranose to UDP-galactofuranose, a precursor of the cell surface β-galactofuranose that is involved in the virulence of the pathogen. The fact that UGM is not present in humans makes inhibition of this enzyme a good approach in the design of new Chagas therapeutics. By performing a series of computer simulations of T. cruzi UGM in the presence or absence of an active site ligand, we address the molecular details of the mechanism that controls the uptake and retention of the substrate. The simulations suggest a modular mechanism in which each moiety of the substrate controls the flexibility of a different protein loop. Furthermore, the calculations indicate that interactions with the substrate diphosphate moiety are especially important for stabilizing the closed active site. This hypothesis is supported with kinetics measurements of site-directed mutants of T. cruzi UGM. Our results extend our knowledge of UGM dynamics and offer new alternatives for the prospective design of drugs.

w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations

Free energy governs the equilibrium extent of many biological processes. High barriers separating free energy minima often limit the sampling in molecular dynamics (MD) simulations, leading to inaccurate free energies. Here, we demonstrate enhanced sampling and improved free energy calculations, relative to conventional MD, using windowed accelerated MD within a Hamiltonian replica exchange framework (w-REXAMD). We show that for a case in which multiple conformations are separated by large free energy barriers, w-REXAMD is a useful enhanced sampling technique, without any necessary reweighting.

Exploring the Photophysical Properties of Molecular Systems Using Excited State Accelerated ab Initio Molecular Dynamics

In the present work, we employ excited state accelerated ab initio molecular dynamics (A-AIMD) to efficiently study the excited state energy landscape and photophysical topology of a variety of molecular systems. In particular, we focus on two important challenges for the modeling of excited electronic states: (i) the identification and characterization of conical intersections and crossing seams, in order to predict different and often competing radiationless decay mechanisms, and (ii) the description of the solvent effect on the absorption and emission spectra of chemical species in solution. In particular, using as examples the Schiff bases formaldimine and salicylidenaniline, we show that A-AIMD can be readily employed to explore the conformational space around crossing seams in molecular systems with very different photochemistry. Using acetone in water as an example, we demonstrate that the enhanced configurational space sampling may be used to accurately and efficiently describe both the prominent features and line-shapes of absorption and emission spectra.