2010
DOI: 10.1063/1.3474951
|View full text |Cite
|
Sign up to set email alerts
|

Molecular dynamics simulation with weak coupling to heat and material baths

Abstract: A method for performing molecular dynamics simulation in the grand canonical ensemble is developed. The molecular dynamics, with coupling to an external bath, simulation method of [Berendsen et al., J. Chem. Phys. 81, 3684 (1984)] is extended for this purpose. Here the physical system of interest consists of real indistinguishable particles plus one fractional particle, whose potential energy of interaction with the rest of particles is scaled by a coupling parameter, ranging dynamically between zero and one. … Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
2
1

Citation Types

0
45
0

Year Published

2011
2011
2024
2024

Publication Types

Select...
10

Relationship

0
10

Authors

Journals

citations
Cited by 53 publications
(45 citation statements)
references
References 43 publications
0
45
0
Order By: Relevance
“…Finally, the generated output files from CHARMM-GUI were used for the next MD simulation step using GROMACS. Initially, the energy of the system was minimized using steepest descent method, followed by equilibration at room temperature of 300 K and 1 bar pressure, with the Berendsen weak coupling method to regulate the temperature and pressure of the system 77 . Finally, MD simulation production run was carried out for 50 ns to check the stability of the protein alone and in complex with the ligands.…”
Section: Methodsmentioning
confidence: 99%
“…Finally, the generated output files from CHARMM-GUI were used for the next MD simulation step using GROMACS. Initially, the energy of the system was minimized using steepest descent method, followed by equilibration at room temperature of 300 K and 1 bar pressure, with the Berendsen weak coupling method to regulate the temperature and pressure of the system 77 . Finally, MD simulation production run was carried out for 50 ns to check the stability of the protein alone and in complex with the ligands.…”
Section: Methodsmentioning
confidence: 99%
“…Finally 85 ns MD simulation for equilibration and production was carried out for the system with 3D periodic boundary in NPT ensemble imposed. Weak coupling methods [44] were applied with coupling constants of 0.1 ps and 0.5 ps respectively to maintain the systems at 298 K and 1 Bar. Bonds involving hydrogen were constained with LINCS algorithm for protein and SETTLE algorithm for water [45,46], thus permitted using a integrate step of 2 fs.…”
Section: Methodsmentioning
confidence: 99%
“…The ‘Standard Dynamics Cascade’ Protocol [21] of DS 2.5 was implemented with a time step of 1 fs. The initial minimization of the implicitly solvated protein–ligand complex was carried out in two steps on residues of a selected sphere of a 13 Å radius inscribing the ligand, binding site and its surrounding residues and loops, keeping about two thirds of the protein domain fixed throughout the whole dynamics simulation cascade.…”
Section: Materials and Methodologymentioning
confidence: 99%