Adaptive time stepping in biomolecular dynamics

Low-mass molecular dynamics simulation for configurational sampling enhancement: More evidence and theoretical explanation

It has been reported recently that classical, isothermal-isobaric molecular dynamics (NTP MD) simulations at a time step of 1.00 fs of the standard-mass time (Δt=1.00 fssmt) and a temperature of ≤340 K using uniformly reduced atomic masses by tenfold offers better configurational sampling than standard-mass NTP MD simulations at the same time step. However, it has long been reported that atomic masses can also be increased to improve configurational sampling because higher atomic masses permit the use of a longer time step. It is worth investigating whether standard-mass NTP MD simulations at Δt=2.00 or 3.16 fssmt can offer better or comparable configurational sampling than low-mass NTP MD simulations at Δt=1.00 fssmt. This article reports folding simulations of two β-hairpins showing that the configurational sampling efficiency of NTP MD simulations using atomic masses uniformly reduced by tenfold at Δt=1.00 fssmt is statistically equivalent to and better than those using standard masses at Δt=3.16 and 2.00 fssmt, respectively. The results confirm that, relative to those using standard masses at routine Δt=2.00 fssmt, the low-mass NTP MD simulations at Δt=1.00 fssmt are a simple and generic technique to enhance configurational sampling at temperatures of ≤340 K.

Discrete molecular dynamics: an efficient and versatile simulation method for fine protein characterization

Until now it has been impractical to observe protein folding in silico for proteins larger than 50 residues. Limitations of both force field accuracy and computational efficiency make the folding problem very challenging. Here we employ discrete molecular dynamics (DMD) simulations with an all-atom force field to fold fast-folding proteins. We extend the DMD force field by introducing long-range electrostatic interactions to model salt-bridges and a sequence-dependent semiempirical potential accounting for natural tendencies of certain amino acid sequences to form specific secondary structures. We enhance the computational performance by parallelizing the DMD algorithm. Using a small number of commodity computers, we achieve sampling quality and folding accuracy comparable to the explicit-solvent simulations performed on high-end hardware. We demonstrate that DMD can be used to observe equilibrium folding of villin headpiece and WW domain, study two-state folding kinetics, and sample near-native states in ab initio folding of proteins of ∼100 residues.

Evaluation of the regioselective delactonization of tri-sialic acid lactone by in-solution molecular dynamics simulation

An approximate model for the delactonization of tri-sialic acid lactone is presented with two water-layers that led to neutral hydrolysis of δ-lactone. The hydrolytic reactivity was studied with a 10-ns in-solution molecular dynamics simulation. The initial step of this hydrolysis involves a reactant water nucleophile complex via a proton transfer with another water molecule. Therefore, the probability of water molecules localized at the hydrolytic center correlates to the hydrolysis of δ-lactone. The stepwise delactonization of α2,8-(NeuAc)(3) lactone results/resulted from water concentration discrepancy near the carbonyl carbon of lactones in two water oxygen···carbonyl carbon shells, and the distances of OC···O(water) layers were 2.8 Å and 5.1 Å. Based on in-solution molecular dynamics study, the motion of water molecules over the re-face of the carbonyl groups was used for the quantitative description of the residence probability, p, whose value is 0.11 for lactone I and 0.33 for lactone II. The geometric criteria used to determine the residence statistics are (1) the distance of water-oxygen···carbonyl carbon in less than 5.1 Å and (2) the cone angle, θ, of carbonyl OC···O(water) in the range of 85-115°. As expected, a higher residence probability at lactone II led to its faster hydrolysis. Both the radial g(r) and angular p(θ) pair distribution functions of water oxygen and carbonyl groups of lactones ensure a better surrounding hydration encounter for lactone II. In contrast, water molecules around lactone I are deduced due to a steric hindrance by the turn structure of α2,8-(NeuAc)(3) lactone.

Methods for Monte Carlo simulations of biomacromolecules

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.