BioSimSpace.Metadynamics.CollectiveVariable.RMSD#

class BioSimSpace.Metadynamics.CollectiveVariable.RMSD(system, reference, rmsd_indices, reference_index=None, hill_width=0.1000 nm, lower_bound=None, upper_bound=None, grid=None, alignment_type='optimal', pbc=True, property_map={})[source]#

A class for a root-mean-square deviation (RMSD) collective variable.

__init__(system, reference, rmsd_indices, reference_index=None, hill_width=0.1000 nm, lower_bound=None, upper_bound=None, grid=None, alignment_type='optimal', pbc=True, property_map={})[source]#

Constructor.

Parameters:
  • system (System) – The molecular system of interest.

  • reference (Molecule) – The reference molecule, against which the RMSD will be measured. This molecule should match with a single molecule from the system, i.e. contain the same residues as the matching molecule in the same order.

  • rmsd_indices ([int]) – The indices of the atoms within the reference that should be used when calculating the RMSD. Atoms from the reference that are not in rmsd_indices will instead be used for alignment.

  • reference_index (int) – The index of the molecule in the system that matches the reference. If none, then we will assume that the molecule with the closest number of residues is the match. If not all atoms from the reference are matched in the system, then an exception will be thrown.

  • hill_width (Length) – The width of the Gaussian hill used to sample this variable.

  • lower_bound (Bound) – A lower bound on the value of the collective variable.

  • upper_bound (Bound) – An upper bound on the value of the collective variable.

  • grid (Grid) – The grid on which the collective variable will be sampled. This can help speed up long metadynamics simulations where the number of Gaussian kernels can become prohibitive.

  • alignment_type (str) – The mannier in which RMSD alignment is performed. Options are “optimal” or “simple”.

  • pbc (bool) – Whether to use periodic boundary conditions when computing the collective variable.

  • property_map (dict) – A dictionary that maps system “properties” to their user defined values. This allows the user to refer to properties with their own naming scheme, e.g. { “charge” : “my-charge” }

Methods

__init__(system, reference, rmsd_indices[, ...])

Constructor.

getAlignmentType()

Return the RMSD alignment type.

getGrid()

Get the grid on which the collective variable is sampled.

getHillWidth()

Return the width of the Gaussian hill used to bias this collective variable.

getInitialValue()

Return the initial value of the collective variable.

getLowerBound()

Get the lower bound on the collective variable.

getPeriodicBoundaries()

Return whether to take account of periodic boundary conditions when computing the collective variable.

getReferenceIndex()

Return the index of the molecule in the system that corresponds to the reference.

getReferencePDB()

Return the reference PDB file as a list of strings.

getUpperBound()

Get the upper bound on the collective variable.

nComponents()

Return the number of components for the collective variable.

setAlignmentType(alignment_type)

Set the RMSD alignment type.

setGrid([grid])

Set a grid on which the collective variable will be sampled.

setHillWidth(hill_width)

Set the width of the Gaussian hills used to bias this collective variable.

setLowerBound([lower_bound])

Set a lower bound on the value of the collective variable.

setPeriodicBoundaries(pbc)

Set whether to use periodic_boundaries when calculating the collective variable.

setUpperBound([upper_bound])

Set an upper bound on the value of the collective variable.

getAlignmentType()[source]#

Return the RMSD alignment type.

Returns:

alignment_type – The RMSD alignment type.

Return type:

str

getGrid()#

Get the grid on which the collective variable is sampled.

Returns:

grid – The grid on which the collective variable is sampled.

Return type:

Grid

getHillWidth()[source]#

Return the width of the Gaussian hill used to bias this collective variable.

Returns:

hill_width – The width of the Gaussian hill.

Return type:

Length

getInitialValue()[source]#

Return the initial value of the collective variable.

Returns:

rmsd – The initial value of the collective variable.

Return type:

Length

getLowerBound()#

Get the lower bound on the collective variable.

Returns:

lower_bound – The lower bound on the value of the collective variable.

Return type:

Bound

getPeriodicBoundaries()[source]#

Return whether to take account of periodic boundary conditions when computing the collective variable.

Returns:

pbc – Whether to use periodic boundaries conditions.

Return type:

bool

getReferenceIndex()[source]#

Return the index of the molecule in the system that corresponds to the reference.

Returns:

reference_index – The index of the molecule in the system that corresponds to the reference.

Return type:

int

getReferencePDB()[source]#

Return the reference PDB file as a list of strings.

Returns:

pdb – The reference PDB file as list of strings.

Return type:

[str]

getUpperBound()#

Get the upper bound on the collective variable.

Returns:

upper_bound – The upper bound on the value of the collective variable.

Return type:

Bound

nComponents()#

Return the number of components for the collective variable.

setAlignmentType(alignment_type)[source]#

Set the RMSD alignment type. Options are “optimal” or “simple”.

Parameters:

alignment_type (str) – The RMSD alignment type.

setGrid(grid=None)#

Set a grid on which the collective variable will be sampled. Call with no arguments to clear the grid.

Parameters:

grid (Grid) – A grid for the collective variable.

setHillWidth(hill_width)[source]#

Set the width of the Gaussian hills used to bias this collective variable.

hill_widthLength

The width of the Gaussian hill.

setLowerBound(lower_bound=None)#

Set a lower bound on the value of the collective variable.

Parameters:

lower_bound (Bound) – A lower bound on the value of the collective variable.

setPeriodicBoundaries(pbc)[source]#

Set whether to use periodic_boundaries when calculating the collective variable.

Parameters:

pbc (bool) – Whether to use periodic boundaries conditions.

setUpperBound(upper_bound=None)#

Set an upper bound on the value of the collective variable.

Parameters:

upper_bound (Bound) – An upper bound on the value of the collective variable.