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

Molecular simulation is an extremely useful, but computationally very expensive tool for studies of chemical and biomolecular systems Here, we present a new implementation of our molecular simulation toolkit GROMACS which now both achieves extremely high performance on single processors from algorithmic optimizations and hand-coded routines and simultaneously scales very well on parallel machines The code encompasses a minimal-communication domain decomposition algorithm, full dynamic load balancing, a state-of-the-art parallel constraint solver, and efficient virtual site algorithms that allow removal of hydrogen atom degrees of freedom to enable integration time steps up to 5 fs for atomistic simulations also in parallel To improve the scaling properties of the common particle mesh Ewald electrostatics algorithms, we have in addition used a Multiple-Program, Multiple-Data approach, with separate node domains responsible for direct and reciprocal space interactions Not only does this combination of a

/pdf/gromacs-4-algorithms-for-highly-efficient-load-balanced-and-16xnd83poh.pdf

GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

/pdf/gromacs-fast-flexible-and-free-15aazpel5u.pdf

GROMACS: Fast, flexible, and free

http://kth.diva-portal.org/smash/get/diva2:1303357/FULLTEXT01

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Copernicus, a hybrid dataflow and peer-to-peer scientific computing platform for efficient large-scale ensemble sampling

Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mechanical work. They are responsible for a spectrum of physiological processes in living organisms, including electrical signaling and cell-cycle progression. Although the mechanism of voltage-sensing has been well characterized for some membrane proteins, including voltage-gated ion channels, even the location of the voltage-sensing elements remains unknown for others. Moreover, the detection of these elements by using experimental techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estimation of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis experiments will significantly reduce the number of potential voltage-sensitive elements to be tested.

/pdf/determining-the-molecular-basis-of-voltage-sensitivity-in-2dwg43es54.pdf

Determining the molecular basis of voltage sensitivity in membrane proteins.

The effect of the inhomogeneous environment upon solvent molecules close to a macromolecular surface is evaluated from a molecular-dynamics simulation of a protein, myoglobin, in water solution. The simulation is analyzed in terms of a mean-field potential from the protein upon the water molecules and spatially varying translational diffusion coefficients for solvent molecules in directions parallel and perpendicular to the protein surface. The diffusion coefficients can be obtained from the slope of the average-square displacements vs time, as well as from the integral of the velocity autocorrelation functions. It is shown that the former procedure gives a lot of ambiguities due to the variation of the slope of the curve with time. The latter, however, after analytic correction for the contribution from algebraic long-time tails, furnish a much more reliable alternative.

Solvent diffusion outside macromolecular surfaces

Pentameric ligand-gated ion channels (pLGICs) are allosteric receptors that mediate rapid electrochemical signal transduction in the animal nervous system through the opening of an ion pore upon binding of neurotransmitters. Orthologs have been found and characterized in prokaryotes and they display highly similar structure–function relationships to eukaryotic pLGICs; however, they often encode greater architectural diversity involving additional amino-terminal domains (NTDs). Here we report structural, functional, and normal-mode analysis of two conformational states of a multidomain pLGIC, called DeCLIC, from a Desulfofustis deltaproteobacterium, including a periplasmic NTD fused to the conventional ligand-binding domain (LBD). X-ray structure determination revealed an NTD consisting of two jelly-roll domains interacting across each subunit interface. Binding of Ca2+ at the LBD subunit interface was associated with a closed transmembrane pore, with resolved monovalent cations intracellular to the hydrophobic gate. Accordingly, DeCLIC-injected oocytes conducted currents only upon depletion of extracellular Ca2+; these were insensitive to quaternary ammonium block. Furthermore, DeCLIC crystallized in the absence of Ca2+ with a wide-open pore and remodeled periplasmic domains, including increased contacts between the NTD and classic LBD agonist-binding sites. Functional, structural, and dynamical properties of DeCLIC paralleled those of sTeLIC, a pLGIC from another symbiotic prokaryote. Based on these DeCLIC structures, we would reclassify the previous structure of bacterial ELIC (the first high-resolution structure of a pLGIC) as a “locally closed” conformation. Taken together, structures of DeCLIC in multiple conformations illustrate dramatic conformational state transitions and diverse regulatory mechanisms available to ion channels in pLGICs, particularly involving Ca2+ modulation and periplasmic NTDs.

/pdf/structural-basis-for-allosteric-transitions-of-a-multidomain-2cbk1zfkb6.pdf

Erik Lindahl

Papers

Copernicus, a hybrid dataflow and peer-to-peer scientific computing platform for efficient large-scale ensemble sampling

Determining the molecular basis of voltage sensitivity in membrane proteins.

Solvent diffusion outside macromolecular surfaces

Structural basis for allosteric transitions of a multidomain pentameric ligand-gated ion channel

Author response: Accelerated cryo-EM structure determination with parallelisation using GPUs in RELION-2