Desmond (software)

Last updated
Desmond
Developer(s) D. E. Shaw Research
Operating system Linux
Platform x86, x86-64, computer clusters
Available inEnglish
Type Computational chemistry
License Proprietary freeware, commercial software
Website www.deshawresearch.com/resources_desmond.html , schrodinger.com/desmond

Desmond is a software package developed at D. E. Shaw Research to perform high-speed molecular dynamics simulations of biological systems on conventional computer clusters. [1] [2] [3] [4] The code uses novel parallel algorithms [5] and numerical methods [6] to achieve high performance on platforms containing multiple processors, [7] but may also be executed on a single computer.

Contents

The core and source code are available at no cost for non-commercial use by universities and other not-for-profit research institutions, and have been used in the Folding@home distributed computing project. Desmond is available as commercial software through Schrödinger, Inc.

Molecular dynamics program

Desmond supports algorithms typically used to perform fast and accurate molecular dynamics. Long-range electrostatic energy and forces can be calculated using particle mesh Ewald-based methods. [8] [9] Constraints can be enforced using the M-SHAKE algorithm. These methods can be used together with time-scale splitting (RESPA-based) integration schemes.

Desmond can compute energies and forces [10] for many standard fixed-charged force fields used in biomolecular simulations, and is also compatible with polarizable force fields based on the Drude formalism. A variety of integrators and support for various ensembles have been implemented in the code, including methods for temperature control (Andersen thermostat, Nosé-Hoover, and Langevin) and pressure control (Berendsen, Martyna-Tobias-Klein, and Langevin). The code also supports methods for restraining atomic positions and molecular configurations; allows simulations to be carried out using a variety of periodic cell configurations; and has facilities for accurate checkpointing and restart.

Desmond can also be used to perform absolute and relative free energy calculations (e.g., free energy perturbation). Other simulation methods (such as replica exchange) are supported through a plug-in-based infrastructure, which also allows users to develop their own simulation algorithms and models.

Desmond is also available in a graphics processing unit (GPU) accelerated version that is about 60-80 times faster than the central processing unit (CPU) version.

Along with the molecular dynamics program, the Desmond software also includes tools for minimizing and energy analysis, both of which can be run efficiently in a parallel environment.

Force fields parameters can be assigned using a template-based parameter assignment tool called Viparr [11] . It currently supports several versions of the CHARMM, Amber and OPLS force fields, and a range of different water models.

Desmond is integrated with a molecular modeling environment (Maestro, developed by Schrödinger, Inc.) for setting up simulations of biological and chemical systems, and is compatible with Visual Molecular Dynamics (VMD) for trajectory viewing and analysis.

See also

Related Research Articles

<span class="mw-page-title-main">Computational chemistry</span> Branch of chemistry

Computational chemistry is a branch of chemistry that uses computer simulations to assist in solving chemical problems. It uses methods of theoretical chemistry incorporated into computer programs to calculate the structures and properties of molecules, groups of molecules, and solids. The importance of this subject stems from the fact that, with the exception of some relatively recent findings related to the hydrogen molecular ion, achieving an accurate quantum mechanical depiction of chemical systems analytically, or in a closed form, is not feasible. The complexity inherent in the many-body problem exacerbates the challenge of providing detailed descriptions of quantum mechanical systems. While computational results normally complement information obtained by chemical experiments, it can occasionally predict unobserved chemical phenomena.

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

Ewald summation, named after Paul Peter Ewald, is a method for computing long-range interactions in periodic systems. It was first developed as the method for calculating the electrostatic energies of ionic crystals, and is now commonly used for calculating long-range interactions in computational chemistry. Ewald summation is a special case of the Poisson summation formula, replacing the summation of interaction energies in real space with an equivalent summation in Fourier space. In this method, the long-range interaction is divided into two parts: a short-range contribution, and a long-range contribution which does not have a singularity. The short-range contribution is calculated in real space, whereas the long-range contribution is calculated using a Fourier transform. The advantage of this method is the rapid convergence of the energy compared with that of a direct summation. This means that the method has high accuracy and reasonable speed when computing long-range interactions, and it is thus the de facto standard method for calculating long-range interactions in periodic systems. The method requires charge neutrality of the molecular system to accurately calculate the total Coulombic interaction. A study of the truncation errors introduced in the energy and force calculations of disordered point-charge systems is provided by Kolafa and Perram.

Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.

<span class="mw-page-title-main">Water model</span> Aspect of computational chemistry

In computational chemistry, a water model is used to simulate and thermodynamically calculate water clusters, liquid water, and aqueous solutions with explicit solvent. The models are determined from quantum mechanics, molecular mechanics, experimental results, and these combinations. To imitate a specific nature of molecules, many types of models have been developed. In general, these can be classified by the following three points; (i) the number of interaction points called site, (ii) whether the model is rigid or flexible, (iii) whether the model includes polarization effects.

In computational chemistry, a constraint algorithm is a method for satisfying the Newtonian motion of a rigid body which consists of mass points. A restraint algorithm is used to ensure that the distance between mass points is maintained. The general steps involved are: (i) choose novel unconstrained coordinates, (ii) introduce explicit constraint forces, (iii) minimize constraint forces implicitly by the technique of Lagrange multipliers or projection methods.

In computer simulations of mechanical systems, energy drift is the gradual change in the total energy of a closed system over time. According to the laws of mechanics, the energy should be a constant of motion and should not change. However, in simulations the energy might fluctuate on a short time scale and increase or decrease on a very long time scale due to numerical integration artifacts that arise with the use of a finite time step Δt. This is somewhat similar to the flying ice cube problem, whereby numerical errors in handling equipartition of energy can change vibrational energy into translational energy.

<span class="mw-page-title-main">CP2K</span>

CP2K is a freely available (GPL) quantum chemistry and solid state physics program package, written in Fortran 2008, to perform atomistic simulations of solid state, liquid, molecular, periodic, material, crystal, and biological systems. It provides a general framework for different methods: density functional theory (DFT) using a mixed Gaussian and plane waves approach (GPW) via LDA, GGA, MP2, or RPA levels of theory, classical pair and many-body potentials, semi-empirical and tight-binding Hamiltonians, as well as Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid schemes relying on the Gaussian Expansion of the Electrostatic Potential (GEEP). The Gaussian and Augmented Plane Waves method (GAPW) as an extension of the GPW method allows for all-electron calculations. CP2K can do simulations of molecular dynamics, metadynamics, Monte Carlo, Ehrenfest dynamics, vibrational analysis, core level spectroscopy, energy minimization, and transition state optimization using NEB or dimer method.

D. E. Shaw Research (DESRES) is a privately held biochemistry research company based in New York City. Under the scientific direction of David E. Shaw, the group's chief scientist, D. E. Shaw Research develops technologies for molecular dynamics simulations and applies such simulations to basic scientific research in structural biology and biochemistry, and to the process of computer-aided drug design.

<span class="mw-page-title-main">Anton (computer)</span> Supercomputer designed and built by D. E. Shaw Research

Anton is a massively parallel supercomputer designed and built by D. E. Shaw Research in New York, first running in 2008. It is a special-purpose system for molecular dynamics (MD) simulations of proteins and other biological macromolecules. An Anton machine consists of a substantial number of application-specific integrated circuits (ASICs), interconnected by a specialized high-speed, three-dimensional torus network.

Particle–Particle–Particle–Mesh (P3M) is a Fourier-based Ewald summation method to calculate potentials in N-body simulations.

<span class="mw-page-title-main">Metadynamics</span> Scientific computer simulation method

Metadynamics is a computer simulation method in computational physics, chemistry and biology. It is used to estimate the free energy and other state functions of a system, where ergodicity is hindered by the form of the system's energy landscape. It was first suggested by Alessandro Laio and Michele Parrinello in 2002 and is usually applied within molecular dynamics simulations. MTD closely resembles a number of newer methods such as adaptively biased molecular dynamics, adaptive reaction coordinate forces and local elevation umbrella sampling. More recently, both the original and well-tempered metadynamics were derived in the context of importance sampling and shown to be a special case of the adaptive biasing potential setting. MTD is related to the Wang–Landau sampling.

Local elevation is a technique used in computational chemistry or physics, mainly in the field of molecular simulation. It was developed in 1994 by Huber, Torda and van Gunsteren to enhance the searching of conformational space in molecular dynamics simulations and is available in the GROMOS software for molecular dynamics simulation. The method was, together with the conformational flooding method, the first to introduce memory dependence into molecular simulations. Many recent methods build on the principles of the local elevation technique, including the Engkvist-Karlström, adaptive biasing force, Wang–Landau, metadynamics, adaptively biased molecular dynamics, adaptive reaction coordinate forces, and local elevation umbrella sampling methods. The basic principle of the method is to add a memory-dependent potential energy term in the simulation so as to prevent the simulation to revisit already sampled configurations, which leads to the increased probability of discovering new configurations. The method can be seen as a continuous variant of the Tabu search method.

Path integral molecular dynamics (PIMD) is a method of incorporating quantum mechanics into molecular dynamics simulations using Feynman path integrals. In PIMD, one uses the Born–Oppenheimer approximation to separate the wavefunction into a nuclear part and an electronic part. The nuclei are treated quantum mechanically by mapping each quantum nucleus onto a classical system of several fictitious particles connected by springs governed by an effective Hamiltonian, which is derived from Feynman's path integral. The resulting classical system, although complex, can be solved relatively quickly. There are now a number of commonly used condensed matter computer simulation techniques that make use of the path integral formulation including Centroid Molecular Dynamics (CMD), Ring Polymer Molecular Dynamics (RPMD), and the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method. The same techniques are also used in path integral Monte Carlo (PIMC).

Adaptive sampling is a technique used in computational molecular biology to efficiently simulate protein folding when coupled with molecular dynamics simulations.

Surface hopping is a mixed quantum-classical technique that incorporates quantum mechanical effects into molecular dynamics simulations. Traditional molecular dynamics assume the Born-Oppenheimer approximation, where the lighter electrons adjust instantaneously to the motion of the nuclei. Though the Born-Oppenheimer approximation is applicable to a wide range of problems, there are several applications, such as photoexcited dynamics, electron transfer, and surface chemistry where this approximation falls apart. Surface hopping partially incorporates the non-adiabatic effects by including excited adiabatic surfaces in the calculations, and allowing for 'hops' between these surfaces, subject to certain criteria.

<span class="mw-page-title-main">Ali Alavi</span> British theoretical chemist

Ali Alavi FRS is a professor of theoretical chemistry in the Department of Chemistry at the University of Cambridge and a Director of the Max Planck Institute for Solid State Research in Stuttgart.

<span class="mw-page-title-main">NhaA family</span> Family of transport proteins

Na+/H+ antiporter A (NhaA) family (TC# 2.A.33) contains a number of bacterial sodium-proton antiporter (SPAP) proteins. These are integral membrane proteins that catalyse the exchange of H+ for Na+ in a manner that is highly pH dependent. Homologues have been sequenced from a number of bacteria and archaea. Prokaryotes possess multiple paralogues. A representative list of the proteins that belong to the NhaA family can be found in the Transporter Classification Database.

<span class="mw-page-title-main">Mixed quantum-classical dynamics</span> Computational chemistry methods to simulate non-adiabatic processes

Mixed quantum-classical (MQC) dynamics is a class of computational theoretical chemistry methods tailored to simulate non-adiabatic (NA) processes in molecular and supramolecular chemistry. Such methods are characterized by:

  1. Propagation of nuclear dynamics through classical trajectories;
  2. Propagation of the electrons through quantum methods;
  3. A feedback algorithm between the electronic and nuclear subsystems to recover nonadiabatic information.

References

  1. Bowers, Kevin J.; Chow, Edmond; Xu, Huafeng; Dror, Ron O.; Eastwood, Michael P.; Gregersen, Brent A.; Klepeis, John L.; Kolossvary, Istvan; Moraes, Mark A.; Sacerdoti, Federico D.; Salmon, John K.; Shan, Yibing; Shaw, David E. (2006). "Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters" (PDF). ACM/IEEE SC 2006 Conference (SC'06). p. 43. doi:10.1109/SC.2006.54. ISBN   978-0-7695-2700-0. Archived from the original (PDF) on 2008-08-28. Retrieved 2009-01-16.
  2. Jensen, M. O.; Borhani, D. W.; Lindorff-Larsen, K.; Maragakis, P.; Jogini, V.; Eastwood, M. P.; Dror, R. O.; Shaw, D. E. (2010). "Principles of conduction and hydrophobic gating in K+ channels". Proceedings of the National Academy of Sciences. 107 (13): 5833–5838. Bibcode:2010PNAS..107.5833J. doi: 10.1073/pnas.0911691107 . PMC   2851896 . PMID   20231479.
  3. Dror, R. O.; Arlow, D. H.; Borhani, D. W.; Jensen, M. O.; Piana, S.; Shaw, D. E. (2009). "Identification of two distinct inactive conformations of the 2-adrenergic receptor reconciles structural and biochemical observations". Proceedings of the National Academy of Sciences. 106 (12): 4689–4694. Bibcode:2009PNAS..106.4689D. doi: 10.1073/pnas.0811065106 . PMC   2650503 . PMID   19258456.
  4. Shan, Y.; Seeliger, M. A.; Eastwood, M. P.; Frank, F.; Xu, H.; Jensen, M. O.; Dror, R. O.; Kuriyan, J.; Shaw, D. E. (2009). "A conserved protonation-dependent switch controls drug binding in the Abl kinase". Proceedings of the National Academy of Sciences. 106 (1): 139–144. Bibcode:2009PNAS..106..139S. doi: 10.1073/pnas.0811223106 . PMC   2610013 . PMID   19109437.
  5. Bowers, Kevin J.; Dror, Ron O.; Shaw, David E. (2006). "The midpoint method for parallelization of particle simulations". The Journal of Chemical Physics. 124 (18): 184109. Bibcode:2006JChPh.124r4109B. doi: 10.1063/1.2191489 . PMID   16709099.
  6. Lippert, Ross A.; Bowers, Kevin J.; Dror, Ron O.; Eastwood, Michael P.; Gregersen, Brent A.; Klepeis, John L.; Kolossvary, Istvan; Shaw, David E. (2007). "A common, avoidable source of error in molecular dynamics integrators". The Journal of Chemical Physics. 126 (4): 046101. Bibcode:2007JChPh.126d6101L. doi: 10.1063/1.2431176 . PMID   17286520. S2CID   38661350.
  7. Edmond Chow; Charles A. Rendleman; Kevin J. Bowers; Ron O. Dror; Douglas H. Hughes; Justin Gullingsrud; Federico D. Sacerdoti; David E. Shaw (July 2008). "Desmond Performance on a Cluster of Multicore Processors". D. E. Shaw Research Technical Report DESRES/TR--2008-01.{{cite journal}}: Cite journal requires |journal= (help)
  8. Bowers, K.J.; Lippert, R.A.; Dror, R.O.; Shaw, D.E. (2010). "Improved Twiddle Access for Fast Fourier Transforms". IEEE Transactions on Signal Processing. 58 (3): 1122–1130. Bibcode:2010ITSP...58.1122B. doi:10.1109/TSP.2009.2035984. S2CID   17240443.
  9. Shan, Yibing; Klepeis, John L.; Eastwood, Michael P.; Dror, Ron O.; Shaw, David E. (2005). "Gaussian split Ewald: A fast Ewald mesh method for molecular simulation". The Journal of Chemical Physics. 122 (5): 054101. Bibcode:2005JChPh.122e4101S. doi:10.1063/1.1839571. PMID   15740304. S2CID   35865319.
  10. Lindorff-Larsen, Kresten; Piana, Stefano; Palmo, Kim; Maragakis, Paul; Klepeis, John L.; Dror, Ron O.; Shaw, David E. (2010). "Improved side-chain torsion potentials for the Amber ff99SB protein force field". Proteins: Structure, Function, and Bioinformatics. 78 (8): 1950–8. doi:10.1002/prot.22711. PMC   2970904 . PMID   20408171.
  11. https://github.com/DEShawResearch/viparr