MDVRY: a polarizable classical molecular dynamics package for biomolecules

Abstract

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: MDVRY Catalogue identifier: AEBY_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEBY_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 39 156 No. of bytes in distributed program, including test data, etc.: 277 197 Distribution format: tar.bz2 Programming language: C Computer: Linux machines with FFTW Fourier Transform package installed Operating system: Linux machines, SUSE & RedHat distributions Classification: 3, 16.13, 23 External 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.