David van der Spoel

GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation

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 algorithms enable extremely long simulations of large systems but also it provides that simulation performance on quite modest numbers of standard cluster nodes.

GROMACS: Fast, flexible, and free

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.

GROMACS 4.5

MOTIVATION Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources. RESULTS Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations. AVAILABILITY GROMACS is an open source and free software available from http://www.gromacs.org. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Potential for biomolecular imaging with femtosecond X-ray pulses

Sample damage by X-rays and other radiation limits the resolution of structural studies on non-repetitive and non-reproducible structures such as individual biomolecules or cells. Cooling can slow sample deterioration, but cannot eliminate damage-induced sample movement during the time needed for conventional measurements. Analyses of the dynamics of damage formation suggest that the conventional damage barrier (about 200 X-ray photons per Å2 with X-rays of 12 keV energy or 1 Å wavelength) may be extended at very high dose rates and very short exposure times. Here we have used computer simulations to investigate the structural information that can be recovered from the scattering of intense femtosecond X-ray pulses by single protein molecules and small assemblies. Estimations of radiation damage as a function of photon energy, pulse length, integrated pulse intensity and sample size show that experiments using very high X-ray dose rates and ultrashort exposures may provide useful structural information before radiation damage destroys the sample. We predict that such ultrashort, high-intensity X-ray pulses from free-electron lasers that are currently under development, in combination with container-free sample handling methods based on spraying techniques, will provide a new approach to structural determinations with X-rays.

Femtosecond diffractive imaging with a soft-X-ray free-electron laser

Theory predicts1,2,3,4 that, with an ultrashort and extremely bright coherent X-ray pulse, a single diffraction pattern may be recorded from a large macromolecule, a virus or a cell before the sample explodes and turns into a plasma. Here we report the first experimental demonstration of this principle using the FLASH soft-X-ray free-electron laser. An intense 25 fs, 4×1013 W cm−2 pulse, containing 1012 photons at 32 nm wavelength, produced a coherent diffraction pattern from a nanostructured non-periodic object, before destroying it at 60,000 K. A novel X-ray camera assured single-photon detection sensitivity by filtering out parasitic scattering and plasma radiation. The reconstructed image, obtained directly from the coherent pattern by phase retrieval through oversampling5,6,7,8,9, shows no measurable damage, and is reconstructed at the diffraction-limited resolution. A three-dimensional data set may be assembled from such images when copies of a reproducible sample are exposed to the beam one by one10.

A systematic study of water models for molecular simulation: Derivation of water models optimized for use with a reaction field.

We have performed long molecular dynamics simulations of water using four popular water models, namely simple point charge (SPC), extended simple point charge (SPC/E), and the three point (TIP3P) and four point (TIP4P) transferable intermolecular potentials. System sizes of 216 and 820 molecules were used to study the dependence of properties on the system size. All systems were simulated at 300 K with and without reaction fields and with two different cutoff radii, in order to study the impact of the cutoff treatment on density, energy, dynamic, and dielectric properties. Furthermore we generated two special-purpose water models based on the SPC and TIP4P models, for use with a reaction field. The atomic charges and the Lennard-Jones C12 parameter were optimized to reproduce the correct energy and pressure using the weak coupling algorithm for parameters. Indeed, in simulations without parameter coupling of both new models the density and potential energy were found to be close to the experimental values...

Efficient docking of peptides to proteins without prior knowledge of the binding site

Reliability in docking of ligand molecules to proteins or other targets is an important challenge for molecular modeling. Applications of the docking technique include not only prediction of the binding mode of novel drugs, but also other problems like the study of protein‐protein interactions. Here we present a study on the reliability of the results obtained with the popular AutoDock program. We have performed systematical studies to test the ability of AutoDock to reproduce eight different protein/ligand complexes for which the structure was known, without prior knowledge of the binding site. More specifically, we look at factors influencing the accuracy of the final structure, such as the number of torsional degrees of freedom in the ligand. We conclude that the Autodock program package is able to select the correct complexes based on the energy without prior knowledge of the binding site. We named this application blind docking, as the docking algorithm is not able to “see” the binding site but can still find it. The success of blind docking represents an important finding in the era of structural genomics.

A temperature predictor for parallel tempering simulations

An algorithm is proposed that generates a set of temperatures for use in parallel tempering simulations (also known as temperature-replica exchange molecular dynamics simulations) of proteins to obtain a desired exchange probability Pdes. The input consists of the number of protein atoms and water molecules in the system, information about the use of constraints and virtual sites and the lower temperature limits. The temperatures generated yield probabilities which are very close to Pdes (correlation 97%), independent of force field and over a wide temperature range. To facilitate its use, the algorithm has been implemented as a web server at .

Dynamic properties of water/alcohol mixtures studied by computer simulation

We have studied mixtures of alcohol and water in an extensive series of 465 molecular-dynamics simulations with an aggregate length of 713 ns, in order to study excess properties of mixing, in particular the relation between mobility and viscosity. Methanol/water, ethanol/water, and 1-propanol/water mixtures were simulated using an alcohol content of 0–100 mass % in steps of 10%, using the OPLS (optimized potential for liquid simulations) force field for the alcohol molecules and the TIP4P (transferable intermolecular potential with four particles) water model. Computed densities and energies show very good agreement with experimental data for bulk simulations and the mixtures are satisfactory as well. The shear viscosity was computed using nonequilibrium molecular-dynamics simulations. Other properties studied include diffusion constants and rotational correlation times. We find the mobility to correlate well with the viscosity data, i.e., at intermediate alcohol concentrations the viscosity is maximal a...

Blind docking of drug-sized compounds to proteins with up to a thousand residues

Blind docking was introduced for the detection of possible binding sites and modes of peptide ligands by scanning the entire surface of protein targets. In the present study, the method is tested on a group of drug‐sized compounds and proteins with up to a thousand amino acid residues. Both proteins from complex structures and ligand‐free proteins were used as targets. Robustness, limitations and future perspectives of the method are discussed. It is concluded that blind docking can be used for unbiased mapping of the binding patterns of drug candidates.