![]() |
|||||||||
![]() |
|
||||||||
![]() 45485-AC5
|
![]() |
(1) |
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
The performance of the method is demonstrated on a long polymer chain and is shown to outperform replica-exchange Monte Carlo with only one trajectory. To implement REPSWA to polymers and biomolecules, the transformations are applied to the backbone dihedral angles, which are useful for determining the conformational space of the molecule. The comparison between REPSWA, hybrid Monte Carlo and parallel tempering for an all-atom 50-mer alkane chain is shown in Fig. 1 below.
Figure 1: Comparison between HMC, PT with 10 replicas, and REPSWA in their sampling efficiency for a 50-mer alkane chain. Top: N(i), number of barrier crossing events for each dihedral i from a 50-million MC step trajectory. Middle: Ramachandran plot of the central dihedral angles (24 and 25). Bottom: Value of the central dihedral (25) as a function of MC steps.
The figure shows, with one trajectory, REPSWA is able to visit more minima than either hybrid Monte Carlo or parallel tempering.
Earlier, we reported
a set of force field parameters for the
attachment of imidazole to a polyethyleneoxide spacer, as discussed
in the original proposal.
We have now tested this force field on the experimental crystal
structure of the material in order to ensure that the force
field parameters are able to maintain the crystal structure.
Fig. 2 shows a snapshot of the experimental
crystal structure, which was the starting structure for a
molecular dynamics calculation of a system containing
32 molecules (24
4 unit cells) at a temperature of 120 K.
Figure 2: (Left) Snapshot of the initial crystal structure (Right) Snapshot from a simulation at 120 K after 13 ps.
The figure also shows the final structure from the end of a 13 ps run in the NVT ensemble. We see that the crystal structure is well maintained. Fig. 3 shows the NH radial distribution function computed from this molecular dynamics trajectory together with that calculated from the original crystal structure. It can be seen that the peaks line up reasonably well, indicating that our force field is a sufficiently accurate one. The first peak in the NH distribution of the crystal structure is due to partially transferred protons not treated by the force field.
Figure 3: Carbon-carbon (top) and nitrogen-hydrogen (bottom) radial distribution functions from the simulation (solid black) and crystal structure (dashed blue).
Recently, we introduced the
Adiabatic Free Energy Dynamics (AFED) method
for computing free
energy profiles quickly and accurately using a dynamical adiabatic
separation between a set of collective variables or reaction
coordinates and the remaining degrees of freedom of a system.
This approach leads
to a significant gain in efficiency versus traditional
methods and is able to generate
multidimensional free-energy surfaces efficiently.
However, the need for coordinate transformations is a
significant drawback.
We have now extended this approach to circumvent
these via the introduction of a set of extended phase-space
variables, to which the adiabatic coupling and high-temperature
are applied. The figure below shows how a full
free energy surface in the radius of gyration, ,
and number of hydrogen bonds,
evolves (see Fig. 4 for illustration).
The surface is converged in 5 ns.
Figure 4:
Time evolution
of the free-energy
(
,
) for
N-acetyl-tryptophan-methylamide (NATMA) in the gas phase with
the CHARMM22 force field.