Coordinate I/O tends to be the most time-consuming part of trajectory analysis. In addition, many types of analyses (for example two-dimensional RMSD and cluster analysis) require using coordinate frames multiple times. To simplify this, trajectory coordinates may be saved as a separate data set via the loadcrd command or createcrd action. Any action can then be performed on the COORDS data set with the crdaction command. The crdout command can be used to write coordinates to an output trajectory (similar to trajout).
Although COORDS data sets store everything internally with single-precision, they can still use a large amount of memory. Because of this there is a specialized type of COORDS data set called a TRAJ data set (trajectory), which functions exactly like a COORDS data set except all data is stored on disk. TRAJ data sets can be created with the loadtraj command. TRAJ data sets cannot be modified.
There are several analyses that can be performed using COORDS data sets, either as part of the normal analysis list or via the runanalysis command. Note that while these analyses can be run on specified COORDS data sets, if one is not specified a default COORDS data set will be created, made up of frames from trajin commands.
As an example of where this might be useful is in the calculation of atomic positional fluctuations. Previously this required two steps: one to generate an average structure, then a second to rms-fit to that average structure prior to calculating the fluctuations. This can now be done in one pass with the following input:
parm topology.parm7 loadcrd mdcrd.nc # Generate average structure PDB, @CA only crdaction mdcrd.nc average avg.pdb @CA # Load average structure PDB as reference parm avg.pdb reference avg.pdb parm avg.pdb # RMS-fit to average structure PDB crdaction mdcrd.nc rms reference @CA # Calculate atomic fluctuations for @CA only crdaction mdcrd.nc atomicfluct out fluct.dat bfactor @CA
The following COORDS data set commands are available:
|catcrd||Concatenate two or more COORDS sets.|
|combinecrd||Combine two or more COORDS sets.|
|crdaction||Run a single Action on a COORDS set.|
|crdout||Write a COORDS set to a file.|
|createcrd||(Action) Create a COORDS set during a run.|
|emin||Run simple energy minimization on a frame of a COORDS set.|
|graft||Graft part of one COORDS set onto another COORDS set.|
|loadcrd||Create or append to a COORDS set from a file.|
|loadtraj||Create special COORDS set where frames remain on disk.|
|permutedihedrals||Rotate specified dihedral(s) in given COORDS set by specific interval or to random values.|
|prepareforleap||Prepare a structure (usually loaded from a PDB) for processing with LEaP from Amber.|
|reference||Load a single trajectory frame as a reference.|
|rotatedihedral||Rotate specified dihedral to specified value or by given increment.|
|splitcoords||Split molecules in a COORDS set into a trajectory.|