Molecular dynamics time complexity

Figure 6 HIV-1 protease Ab initio molecular dynamics of complexes C(I), C(II) and C. (a) Location of the proton (b)-(d) (Left) 051 H distances plotted 2is a function of time and (right) final (thick line) and starting (thin line) structures of complexes C(I) (b), C(II) (c) and C (d). In (d)(left) only the last 0.9 ps are shown for the S2dce of clarity. (Reproduced with permission from ref. [100], Copyright 2000, Wiley.)...

West, A.M.A., R. Elber, and D. ShaUoway, Extending molecular dynamics time scales with milestoning Example of complex kinetics in a solvated peptide. Journal of Chemical Physics, 2007, 126(14) 1451104. 

R. Kimmich, NMR Diffusometiy Molecular Dynamics in Complex Systems Probed over Many Decades of Time, EPJ Web Conf., 2012, 30, 05003. 

Fig. 7.15 The variation in torsion angles can be effectively represented as a series of dials, where the time corresponds to the distance from the centre of the dial. Data from a molecular dynamics simulation of an intermolecular complex between the enzyme dihydrofolate reductase and a triazine inhibitor [Leach and Klein 1995].

Once the model of a ligand-receptor complex is built, its stability should be evaluated. Simple molecular mechanics optimization of the putative ligand-receptor complex leads only to the identification of the closest local minimum. However, molecular mechanics optimization of molecules lacks two crucial properties of real molecular systems temperature and, consequently, motion. Molecular dynamics studies the time-dependent evolution of coordinates of complex multimolecular systems as a function of inter- and intramolecular interactions (see Chapter 3). Because simulations are usually performed at nonnal temperature (—300 K), relatively low energy barriers, on the order of kT (0.6 kcal), can... 

Molecular dynamics simulations are capable of addressing the self-assembly process at a rudimentary, but often impressive, level. These calculations can be used to study the secondary structure (and some tertiary structure) of large complex molecules. Present computers and codes can handle massive calculations but cannot eliminate concerns that boundary conditions may affect the result. Eventually, continued improvements in computer hardware will provide this added capacity in serial computers development of parallel computer codes is likely to accomplish the goal more quickly. In addition, the development of realistic, time-efficient potentials will accelerate the useful application of dynamic simulation to the self-assembly process. In addition, principles are needed to guide the selec-... 

We have performed also a reaction field DFT/Molecular Dynamics simulation of this system. We found that after an initial time, when the complex oscillates within the cage at R(N-H) 2.0 a.u. and R(N-C1) 6.0 a.u., a small temperature variation is enough for allowing the complex to overcome the small energetic barrier and, with time, the distance between Cl" and the NH4 fragments starts to increase. Extrapolating to a real solution environment, the two fragments will be completely surrounded by water molecules, i.e. in a solution at infinite dilution the two ions are fully solvated. 

Our calculations show that the isomerization of the silyl-alkyl complex to form a V-allyl complex affords a significant stabilization as summarized in Figure 11. TheiV toil3 isomerization of 9a to the anti silyl-allyi complex, 10a-anti, results in a 9.5 kcal/mol stabilization and isomerization to the syn isomer, lOa-syn, results in a 7.2 kcal/mol stabilization. Isomerization of 9b to the anti silyl-allyi complex, 10b-anti, results in a 6.1 kcal/mol stabilization and isomerization to the syn isomer, lOb-.vyn results in a 5.4 kcal/mol stabilization. High temperature (500 °C) molecular dynamics simulations initiated at the V complex, 9a, reveal that the rj1 to T 3 isomerization has a minimal barrier and occurs in the sub-pico time frame. The inter-conversion between the syn and anti isomers has not been examined since both isomers are stereochemically equivalent, however, we expect the barrier to be small. 

---