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

An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms. Timings and accuracies are presented for three large crystalline ionic systems.

Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems

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

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

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

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

/pdf/scalable-molecular-dynamics-with-namd-51o1z7rwp8.pdf

Scalable molecular dynamics with NAMD

A new dynamical model for the interpretation of nuclear magnetic resonance relaxation data is presented. It is based on a normal mode description, treating the low frequencies associated with collective motions as adjustable parameters to optimize agreement between calculated and experimental relaxation order parameters. This model provides a compact representation of many aspects of internal dynamics and characterizes motions affecting different spin pairs in a correlated way. Furthermore, it links together dynamical characteristics of different types of NMR observables and allows one to assess vibrational thermodynamic properties. Sample applications are given for a 25-residue zinc-finger peptide.

Collective NMR relaxation model applied to protein dynamics.

The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.

/pdf/competitive-interaction-of-monovalent-cations-with-dna-from-h21ojwin9k.pdf

Competitive interaction of monovalent cations with DNA from 3D-RISM.

Molecular dynamics simulations of crystals can enlighten interpretation of experimental X-ray crystallography data and elucidate structural dynamics and heterogeneity in biomolecular crystals. Furthermore, because of the direct comparison against experimental data, they can inform assessment of molecular dynamics methods and force fields. We present microsecond scale results for triclinic hen egg-white lysozyme in a supercell consisting of 12 independent unit cells using four contemporary force fields (Amber ff99SB, ff14ipq, ff14SB, and CHARMM 36) in crystalline and solvated states (for ff14SB only). We find the crystal simulations consistent across multiple runs of the same force field and robust to various solvent equilibration schemes. However, convergence is slow compared with solvent simulations. All the tested force fields reproduce experimental structural and dynamic properties well, but Amber ff14SB maintains structure and reproduces fluctuations closest to the experimental model: its average backbone structure differs from the deposited structure by 0.37A; by contrast, the average backbone structure in solution differs from the deposited by 0.65A. All the simulations are affected by a small progressive deterioration of the crystal lattice, presumably due to imperfect modeling of hydrogen bonding and other crystal contact interactions; this artifact is smallest in ff14SB, with average lattice positions deviating by 0.20A from ideal. Side-chain disorder is surprisingly low with fewer than 30% of the nonglycine or alanine residues exhibiting significantly populated alternate rotamers. Our results provide helpful insight into the methodology of biomolecular crystal simulations and indicate directions for future work to obtain more accurate energy models for molecular dynamics.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/pro.2713

Molecular dynamics simulation of triclinic lysozyme in a crystal lattice.

Most theoretical models for NMR relaxation in liquids assume that overall rotational motion can be described as rotational diffusion with a single diffusion tensor. Such models cannot handle motions (such as between “closed” and “open” states of an enzyme, or between conformers of a partially disordered system) where the shape of the molecule (and hence its rotational diffusion behavior) fluctuates. We provide here a formalism for dealing with such problems. The model involves jumps between discrete conformers with different overall diffusion tensors, and a master (rate) equation to describe the transitions between these conformers. Numerical examples are given for a two-site jump model where global and local motions are concerted, showing how the rate of conformational transitions (relative to the rate of rotational diffusion) affects the observed relaxation parameters.

Influence of the coupling of interdomain and overall motions on NMR relaxation.

Spin‐restricted and spin‐unrestricted Xα multiple scattering calculations are reported for iron porphine with D4h symmetry in the absence of axial ligands. A 3A2g ground state is predicted when the radius of the porphyrin ‘‘hole’’ is less than 2.0 A, whereas a quintet ground state is expected for larger Fe–N distances. A large asymmetry is found in the occupations of the three Fe 4p orbitals. This has important implications for the interpretation of Mossbauer quadrupole coupling constants. Magnetic susceptibility and NMR results can be understood in terms of a model that involves spin‐orbit mixing of an excited 3Eg state with the ground 3A2g state; an essential element in explaining the NMR data is that the proton hyperfine constants are not the same in these two states. Tentative assignments are given for observed absorptions in the optical and near‐infrared regions. It is shown that spin‐unrestricted wave functions are particularly suited for the qualitative interpretation of electronic structure of metallophorphyrins.

David A. Case

Papers

Collective NMR relaxation model applied to protein dynamics.

Competitive interaction of monovalent cations with DNA from 3D-RISM.

Molecular dynamics simulation of triclinic lysozyme in a crystal lattice.

Influence of the coupling of interdomain and overall motions on NMR relaxation.

Xα multiple scattering calculations on iron(II) porphine