Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

The previously developed particle mesh Ewald method is reformulated in terms of efficient B‐spline interpolation of the structure factors This reformulation allows a natural extension of the method to potentials of the form 1/rp with p≥1 Furthermore, efficient calculation of the virial tensor follows Use of B‐splines in place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy We demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N) For biomolecular systems with many thousands of atoms this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 A or less

A smooth particle mesh Ewald method

CHARMM (Chemistry at HARvard Macromolecular Mechanics) is a highly flexible computer program which uses empirical energy functions to model macromolecular systems. The program can read or model build structures, energy minimize them by first- or second-derivative techniques, perform a normal mode or molecular dynamics simulation, and analyze the structural, equilibrium, and dynamic properties determined in these calculations. The operations that CHARMM can perform are described, and some implementation details are given. A set of parameters for the empirical energy function and a sample run are included.

CHARMM: A program for macromolecular energy, minimization, and dynamics calculations

6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

A simple atomic-level hydrophobicity scale reveals protein interfacial structure.

Globular proteins undergo structural transitions to denatured states when sufficient thermodynamic state or chemical perturbations are introduced to their native environment. Cold denaturation is a somewhat counterintuitive phenomenon whereby proteins lose their compact folded structure as a result of a temperature drop. The currently accepted explanation for cold denaturation is based on an associated favorable change in the contact free energy between water and nonpolar groups at colder temperatures which would weaken the hydrophobic interaction and is thought to eventually allow polymer entropy to disrupt protein tertiary structure. In this paper we explore how this environmental perturbation leads to changes in the protein hydration and local motions in apomyoglobin. We do this by analyzing changes in protein hydration and protein motion from molecular dynamics simulation trajectories initially at 310 K, followed by a temperature drop to 278 K. We observe an increase in the number of solvent contacts around the protein and, in particular, distinctly around nonpolar atoms. Further analysis shows that the fluctuations of some protein atoms increase with decreasing temperature. This is accompanied by an observed increase in the isothermal compressibility of the protein, indicating an increase in the protein interior interstitial space. Closer inspection reveals that atoms with increased compressibility and larger-than-expected fluctuations are localized within the protein core regions. These results provide insight into a description of the mechanism of cold denaturation. That is, the lower temperature leads to solvent-induced packing defects at the protein surface, and this more favorable water-protein interaction in turn destabilizes the overall protein structure.

Mechanistic elements of protein cold denaturation.

The extended RISM integral equation theory is employed to study molecular conformational equilibria in liquids. Solvent induced torsional free energysurfaces are computed for the solutes n‐butane and 1, 2‐dichloroethane at infinite dilution in apolar and polar model liquids, CCl4 and water. In addition, the intramolecular distribution of conformers for the neat molecular liquidsn‐butane and 1, 2‐dichloroethane are determined. The qualitative features of these free energies agree well with available simulation results. Quantitatively, the absolute solvent effects are consistently overestimated; however, the differences between polar and apolar systems are accurately predicted. The thermodynamic components of the solvent induced free energies for the dilute solutions are analyzed. For the systems examined here, the internal energy dominates the torsional free energies. In particular, for n‐butane in water, the increase in g a u c h e conformers is in accord with the notion of a hydrophobic interaction. However, the thermodynamic behavior found here opposes the usual belief that such association is an entropy driven process.

Molecular conformational equilibria in liquids

We present a dynamical analysis of the solvent in an aqueous solution of two nonpolar atomic solutes which are constrained to an interatomic distance corresponding to a solvent‐separated free energy minimum. The results are obtained from a molecular dynamics simulation using ST2 model water. Molecular mobility for solvent near the solutes is seen to be retarded, as evidenced in translational diffusion and rotational reorientation. These slower net motions are analogous to pure solvent dynamics at a temperature reduced by 10–15 °C. An analysis of intermolecular hydrogen bonding reveals that solvation shell molecules have correspondingly longer bond half‐lives compared to bulk molecules, by a factor of 1.5–2.0. The spectral densities for intermolecular vibrations are computed from translational and rotational velocity autocorrelation functions for shell and bulk motions. These densities are seen to correlate well with the local binding energy distributions.

Solvent molecular dynamics in regions of hydrophobic hydration

The dynamics of water next to hydrophobic groups is critical for several fundamental biochemical processes such as protein folding and amyloid fiber aggregation. Some biomolecular systems, like melittin or other membrane-associated proteins, exhibit extended hydrophobic surfaces. Due to the strain these surfaces impose on the hydrogen (H)-bond network, the water molecules shift from the clathrate-like arrangement observed around small solutes to an anticlathrate-like geometry with some dangling OH bonds pointing toward the surface. Here we examine the water reorientation dynamics next to a model hydrophobic surface through molecular dynamics simulations and analytic modeling. We show that the water OH bonds lying next to the hydrophobic surface fall into two subensembles with distinct dynamical reorientation properties. The first is the OH bonds tangent to the surface; these exhibit a behavior similar to the water OHs around small hydrophobic solutes, i.e. with a moderate reorientational slowdown explained by an excluded volume effect due to the surface. The second is the dangling OHs pointing toward the surface: these are not engaged in any H-bond, reorient much faster than in the bulk, and exhibit an unusual anisotropy decay which becomes negative for delays of a few picoseconds. The H-bond dynamics, i.e. the exchanges between the different configurations, and the resulting anisotropy decays are analyzed within the analytic extended jump model. We also show that a recent spectroscopy technique, two-dimensional time resolved vibrational spectroscopy (2D-IR), can be used to selectively follow the dynamics of dangling OHs, since these are spectrally distinct from H-bonded ones. By computing the first 2D-IR spectra of water next to a hydrophobic surface, we establish a connection between the spectral dynamics and the dynamical properties that we obtain directly from the simulations.

Peter J. Rossky

Papers

A simple atomic-level hydrophobicity scale reveals protein interfacial structure.

Mechanistic elements of protein cold denaturation.

Molecular conformational equilibria in liquids

Solvent molecular dynamics in regions of hydrophobic hydration

Water reorientation, hydrogen-bond dynamics and 2D-IR spectroscopy next to an extended hydrophobic surface