Subsections


structure.rmsd

Image structure

Synopsis

Determine the RMSD between structures.

Defaults

structure.rmsd(pipes=None, models=None, molecules=None, atom_id=None, atomic=False)

Keyword arguments

pipes: The data pipes to determine the RMSD for.

models: The list of models for each data pipe to determine the RMSD for. The number of elements must match the pipes argument. If no models are given, then all will be used.

molecules: The list of molecules for each data pipe to determine the RMSD for. The RMSD will only be calculated for atoms with identical residue name and number and atom name. The number of elements must match the pipes argument. If no molecules are given, then all will be used.

atom_id: The atom identification string of the coordinates of interest.

atomic: A flag which if True will allow for per-atom RMSDs to be additionally calculated.

Description

This allows the root mean squared deviation (RMSD) between all structures to be calculated. The RMSDs for individual structures to the mean structure will be calculated and reported. These values averaged to produce a global RMSD stored in the structural object of the current data pipe. If the `atomic' argument is set, per-atom RMSDs will additionally be calculated stored in spin containers.

Support for multiple structures is provided by the data pipes, model numbers and molecule names arguments. Each data pipe, model and molecule combination will be treated as a separate structure. As only atomic coordinates with the same residue name and number and atom name will be assembled, structures with slightly different atomic structures can be compared. If the list of models is not supplied, then all models of all data pipes will be used. If the optional molecules list is supplied, each molecule in the list will be considered as a separate structure for comparison between each other.

The atom ID string, which uses the same notation as the spin ID, can be used to restrict the coordinates compared to a subset of molecules, residues, or atoms. For example to only use backbone heavy atoms in a protein, set the atom ID to `@N,C,CA,O', assuming those are the names of the atoms in the 3D structural file.

Prompt examples

To determine the RMSD of all models in the current data pipe, simply type:

[numbers=none]
relax> structure.rmsd()

For the backbone heavy atom RMSD of all models in the current data pipe, simply type:

[numbers=none]
relax> structure.rmsd(atom_id='@N,C,CA,O')

To calculate the C-alpha backbone RMSDs of all models in the current data pipe, type:

[numbers=none]
relax> structure.rmsd(atom_id='CA', atomic=True)


The relax user manual (PDF), created 2020-08-26.