Generate the most probable trajectory between two states

Last modified on
Enter your e-mail address if you want to be notified as soon as the calculation completes.
A reminder of the URL for accessing the results will be part of the message.
Note that your results will be deleted after 15 days.
The title is for your own identification of the job.
Only alphanumerical characters, no space.
(Example: 1ANF_ca_profit.pdb)
The coordinate files should be in PDB format, with only a single structure (no multiple models).
The limit is 10000 residues
Only CA atoms for proteins
A selection of 3-4 atoms per residue for nucleic acids
(Example: 1OMP_ca.pdb)
Maximum is 10000 atoms, CA only for proteins.
The two PDB files (initial and final) should have EXACTLY the same number of atoms.
Note: They will be structurally aligned internally.
The length of the simulation T should be such that the total action Stot in the log file does not change very much if T is increased (i.e. the user should check that the asymptotic regime has been reached).
Typically 50-200 (in arbitrary units)
It should be of the order of 100.
See Figure below on the left for illustration
Elastic constants
A rule of thumb is that k × cutoff2 ≈ 10 kcal/Mol
See Figure below on the right for illustration

Works fine up to 1600 atoms

See below for more details on pseudo-φ.
Note that your results will be deleted after 15 days.

The program will output, in addition to the trajectory in PyMOL format, the coordinates of the transition state as well as the total action Stot and the energy of the transition state UL (measured with respect to the energy of the left -initial- state).

Q1 and Q2 will also be found in the logfile, where Q1 is the number of contacts as in the initial state, and Q2 the number of contacts as in the final state, for each point along the trajectory.
Here is an example (adenylate kinase), where the trajectory obtained with UMMS is also presented.

If needed, a Q1-specific - Q2-specific vs Q-common plot can be built upon request. This kind of plot reveals possible unfolding during the transition (see Okasaki et al., 2006).

The Elastic Energy during the transition could be calculated upon request. The attached plot demonstrates that MinActionPath generates a trajectory with less elastic energy than UMMS linearly interpolated one (see curve pink+green vs blue+red curve).
This again concernes adenylate kinase with elastic spring constants of 0.1 kcal/Mol/Å2 for both sides.

Finally, the question of the robustness of the trajectory with respect of the k1/k2 ratio, i.e. the elastic constants for both sides of the reaction, is of interest. Here is a Q1 vs Q2 plot for adenylate kinase and k2/k1 = 5, k1 = k2 and k1/k2 = 2.

An ASCII file ('something_overlap.out') will be generated. The first column corresponds to the frame number along the trajectory. The 106 following columns correspond, in ascending order, to the overlap between the set of difference vectors (with respect to the starting structure) and the N first normal modes (N=1,106, also calculated from the starting structure).
The 6 first columns must always have negligible values compared to the 7th.

A PDB file ('something_pseudophi.pdb') will be generated. It corresponds to the PDB file of the end structure where the B-factor column is replaced with the pseudo-'phi value'. This allows for direct visualisation using either PyMOL or VMD.
To make a long story short, here we monitor the elastic energy of each residue (through its Calpha) along the trajectory. The 'time' at which this local elastic energy is maximal corresponds to its own 'transition time'. This is what is stored in the B-column. Comparing this value to the global transition time will tell you whether this residue has a early or late transition. Other definitions of residue-specific transition time(s) are described in Zhu and Hummer (2009), Biophys. J. 97:2456-63.