MDVRY: a polarizable classical molecular dynamics package for biomolecules

The MDVRY classical molecular dynamics package is presented for the study of biomolecules in the gas and liquid phase. Electrostatic polarization has been implemented in the formalism of point induced dipoles following the model of Thole. Two schemes have been implemented for the calculation of induced dipoles, i.e. resolution of the self-consistent equations and a ‘Car-Parrinello’ dynamical approach. In this latter, the induced dipoles are calculated at each time step of the dynamics through the dynamics of additional degrees of freedom associated with the dipoles. This method saves computer time and allows to study polarized solvated proteins at a very low CPU cost. The program is written in C-language and runs on LINUX machines. A detailed manual of the code is given. The main features of the package are illustrated taking on examples of proteins in the gas phase or immersed in liquid water.

Program summary

Program title: MDVRYCatalogue identifier: AEBY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBY_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 39 156No. of bytes in distributed program, including test data, etc.: 277 197Distribution format: tar.bz2Programming language: CComputer: Linux machines with FFTW Fourier Transform package installedOperating system: Linux machines, SUSE & RedHat distributionsClassification: 3, 16.13, 23External routines: FFTW (http://www.fftw.org/)Nature of problem: Molecular Dynamics Software package.Solution method: Velocity Verlet algorithm. The implemented force field is composed of intra-molecular interactions and inter-molecular interactions (electrostatics, polarization, van der Waals). Polarization is accounted through induced point dipoles at each atomic site. Supplementary degrees of freedom are associated to the induced dipoles so that a modified Hamiltonian of the dynamics is written. This allows to calculate the induced dipoles with a very fast ‘Car-Parrinello’ type of dynamics.Running time: The test run provided takes approximately 6 minutes to run.