4.5.2. Updated nucleic acid analysis — MDAnalysis.analysis.nucleicacids

Author

Alia Lescoulie

Year

2022

copyright

GNU Public Licence v3

The module provides classes for analyzing nucleic acids structures. This is an updated, higher performance version of previous nucleic acid tools. For applications see [Denning2011] and [Denning2012].

4.5.2.1. Distances

class MDAnalysis.analysis.nucleicacids.NucPairDist(selection1: List[AtomGroup], selection2: List[AtomGroup], **kwargs)[source]

Atom Pair distance calculation base class.

Takes two lists of AtomGroup and computes the distances between them over a trajectory. Used as a superclass for the other nucleic acid distances classes. The distance will be measured between atoms sharing an index in the two lists of AtomGroup.

Parameters

  • selection1 (List[AtomGroup]) – list of AtomGroup containing an atom of each nucleic acid being analyzed.

  • selection1 – list of AtomGroup containing an atom of each nucleic acid being analyzed.

  • kwargs (dict) – arguments for AnalysisBase

results

first index is selection second index is time

Type

numpy.ndarray

results.times

times used in analysis

Type

numpy.ndarray

Raises

ValueError – if the selections given are not the same length

run(start=None, stop=None, step=None, frames=None, verbose=None)

Perform the calculation

Parameters

  • start (int, optional) – start frame of analysis

  • stop (int, optional) – stop frame of analysis

  • step (int, optional) – number of frames to skip between each analysed frame

  • frames (array_like, optional) – array of integers or booleans to slice trajectory

  • verbose (bool, optional) – Turn on verbosity

Changed in version 2.1.0: Added ability to iterate through trajectory by passing a list of frame indices

class MDAnalysis.analysis.nucleicacids.WatsonCrickDist(strand1: List[Residue], strand2: List[Residue], n1_name: str = 'N1', n3_name: str = 'N3', g_name: str = 'G', a_name: str = 'A', u_name: str = 'U', t_name: str = 'T', c_name: str = 'C', **kwargs)[source]

Watson-Crick basepair distance for selected residues over a trajectory.

Takes two lists of Residue objects and calculates the Watson-Crick distance between them over the trajectory. Bases are matched by their index in the lists given as arguments.

Parameters

  • strand1 (List[Residue]) – First list of bases

  • strand2 (List[Residue]) – Second list of bases

  • n1_name (str (optional)) – Name of Nitrogen 1 of nucleic acids by default assigned to N1

  • n3_name (str (optional)) – Name of Nitrogen 3 of nucleic acids by default assigned to N3

  • g_name (str (optional)) – Name of Guanine in topology by default assigned to G

  • a_name (str (optional)) – Name of Adenine in topology by default assigned to G

  • u_name (str (optional)) – Name of Uracil in topology by default assigned to U

  • t_name (str (optional)) – Name of Thymine in topology by default assigned to T

  • c_name (str (optional)) – Name of Cytosine in topology by default assigned to C

  • **kwargs (dict) – arguments for AnalysisBase

results

first index is selection second index is time

Type

numpy.ndarray

results.times

times used in analysis

Type

numpy.ndarray

Raises

  • ValueError – if the residues given are not amino acids

  • ValueError – if the selections given are not the same length

run(start=None, stop=None, step=None, frames=None, verbose=None)

Perform the calculation

Parameters

  • start (int, optional) – start frame of analysis

  • stop (int, optional) – stop frame of analysis

  • step (int, optional) – number of frames to skip between each analysed frame

  • frames (array_like, optional) – array of integers or booleans to slice trajectory

  • verbose (bool, optional) – Turn on verbosity

Changed in version 2.1.0: Added ability to iterate through trajectory by passing a list of frame indices