This module defines a class and a function for anisotropic network model (ANM) calculations.
Class for Anisotropic Network Model (ANM) analysis of proteins ([PD00], [ARA01]).
See a usage example in Anisotropic Network Model (ANM).
[PD00] | Doruker P, Atilgan AR, Bahar I. Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: Application to a-amylase inhibitor. Proteins 2000 40:512-524. |
[ARA01] | Atilgan AR, Durrell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I. Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys. J. 2001 80:505-515. |
Add eigen vector and eigen value pair(s) to the instance. If eigen value is omitted, it will be set to 1. Inverse eigenvalues are set as variances.
Build Hessian matrix for given coordinate set.
Parameters: |
|
---|
Instances of Gamma classes and custom functions are accepted as gamma argument.
When Scipy is available, user can select to use sparse matrices for efficient usage of memory at the cost of computation speed.
Calculate normal modes. This method uses scipy.linalg.eigh() function to diagonalize the Hessian matrix. When Scipy is not found, numpy.linalg.eigh() is used.
Parameters: |
|
---|
Return a copy of eigenvectors array.
Return covariance matrix. If covariance matrix is not set or yet calculated, it will be calculated using available modes.
Return cutoff distance.
Return eigenvalues. For PCA and EDA models built using coordinate data in Å, unit of eigenvalues is Å2. For ANM, GNM, and RTB, on the other hand, eigenvalues are in arbitrary or relative units but they correlate with stiffness of the motion along associated eigenvector.
Return a copy of eigenvectors array.
Return spring constant (or the gamma function or Gamma instance).
Return a copy of the Hessian matrix.
Return a copy of the Kirchhoff matrix.
Return self.
Return title of the model.
Return variances. For PCA and EDA models built using coordinate data in Å, unit of variance is Å2. For ANM, GNM, and RTB, on the other hand, variance is the inverse of the eigenvalue, so it has arbitrary or relative units.
Return True if model is 3-dimensional.
Return number of atoms.
Return number of degrees of freedom.
Return number of modes in the instance (not necessarily maximum number of possible modes).
Set eigen vectors and eigen values. If eigen values are omitted, they will be set to 1. Inverse eigenvalues are set as variances.
Set Hessian matrix. A symmetric matrix is expected, i.e. not a lower- or upper-triangular matrix.
Set title of the model.