Introduction to Molecular Dynamics and Enhanced Sampling Algorithms

2014 ◽  
pp. 24-53
2017 ◽  
Vol 38 (25) ◽  
pp. 2193-2206 ◽  
Author(s):  
Chigusa Kobayashi ◽  
Jaewoon Jung ◽  
Yasuhiro Matsunaga ◽  
Takaharu Mori ◽  
Tadashi Ando ◽  
...  

2015 ◽  
Vol 119 (46) ◽  
pp. 14594-14603 ◽  
Author(s):  
Ole Juul Andersen ◽  
Julie Grouleff ◽  
Perri Needham ◽  
Ross C. Walker ◽  
Frank Jensen

Author(s):  
Yu‐Peng Huang ◽  
Yijie Xia ◽  
Lijiang Yang ◽  
Jiachen Wei ◽  
Yi Isaac Yang ◽  
...  

2021 ◽  
Author(s):  
Fréderic Célerse ◽  
Theo Jaffrelot-Inizan ◽  
Louis Lagardère ◽  
Olivier Adjoua ◽  
Pierre Monmarché ◽  
...  

We detail a novel multi-level enhanced sampling strategy grounded on Gaussian accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs-accelerated implementation within the Tinker-HP molecular dynamics package. We then introduce the new "dual-water" mode and its use with the flexible AMOEBA polarizable force field. By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups thanks to the use of fast multiple--timestep integrators. To further reduce time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD—US/dual-water approach is tested on the 1D Potential of Mean Force (PMF) of the CD2-CD58 system (168000 atoms) allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added enabling AS-GaMD capabilities but also the introduction of the new Adaptive Sampling--US--GaMD (ASUS--GaMD) scheme. The highly parallel ASUS--GaMD setup decreases time to convergence by respectively 10 and 20 compared to GaMD--US and US.


2017 ◽  
Author(s):  
Irfan Alibay ◽  
Kepa K. Burusco ◽  
Neil J. Bruce ◽  
Richard A. Bryce

<p>Determining the conformations accessible to carbohydrate ligands in aqueous solution is important for understanding their biological action. In this work, we evaluate the conformational free energy surfaces of Lewis oligosaccharides in explicit aqueous solvent using a multidimensional variant of the swarm-enhanced sampling molecular dynamics (msesMD) method; we compare with multi-microsecond unbiased MD simulations, umbrella sampling and accelerated MD approaches. For the sialyl Lewis A tetrasaccharide, msesMD simulations in aqueous solution predict conformer landscapes in general agreement with the other biased methods and with triplicate unbiased 10 ms trajectories; these simulations find a predominance of closed conformer and a range of low occupancy open forms. The msesMD simulations also suggest closed-to-open transitions in the tetrasaccharide are facilitated by changes in ring puckering of its GlcNAc residue away from the <sup>4</sup>C<sub>1</sub> form, in line with previous work. For sialyl Lewis X tetrasaccharide, msesMD simulations predict a minor population of an open form in solution, corresponding to a rare lectin-bound pose observed crystallographically. Overall, from comparison with biased MD calculations, we find that triplicate 10 ms unbiased MD simulations may not be enough to fully sample glycan conformations in aqueous solution. However, the computational efficiency and intuitive approach of the msesMD method suggest potential for its application in glycomics as a tool for analysis of oligosaccharide conformation.</p>


2015 ◽  
Vol 11 (7) ◽  
pp. 3446-3454 ◽  
Author(s):  
Isidro Cortes-Ciriano ◽  
Guillaume Bouvier ◽  
Michael Nilges ◽  
Luca Maragliano ◽  
Thérèse E. Malliavin

Sign in / Sign up

Export Citation Format

Share Document