Anisotropic Network Model: Difference between revisions

Revision as of 18:57, 23 November 2014

Latest revision as of 06:45, 12 May 2025

(25 intermediate revisions by 13 users not shown)

Line 1:

{{No footnotes|date=May 2025}}

The '''Anisotropic Network Model''' (ANM) is a simple yet powerful tool made for [[Normal Mode]] Analysis of proteins, which has been successfully applied for exploring the relation between function and dynamics for many proteins. It is essentially an Elastic Network Model for the Cα atoms with a step function for the dependence of the force constants on the inter-particle distance.

[[File:Elastic network model.png|thumb|250x250px|Anisotrpic Network Model use an elastic mass-and-spring network to represent biological macromolecule ([[Elastic Network Model]])]]

The '''Anisotropic Network Model''' (ANM)''' is a simple yet powerful tool made for [[Normalnormal Modemode]] Analysisanalysis of proteins, which has been successfully applied for exploring the relation between function and dynamics for many proteins[[protein]]s. It is essentially an Elastic Network Model for the [[Alpha and beta carbon|Cα atoms]] with a step function for the dependence of the force constants on the inter-particle distance.

== Theory ==

The Anisotropic Network Model was introduced in 2000 (Atilgan et al., 2001; Doruker et al., 2000), inspired by the pioneering work of Tirion (1996), succeeded by the development of the [[Gaussian network model]] (GNM) (Bahar et al., 1997; Haliloglu et al., 1997), and by the work of Hinsen (1998) who first demonstrated the validity of performing EN NMA at residue level. <br /> It represents the biological macromolecule as an elastic mass-and-spring network (Figure.1.), to explain the internal motions of a protein subject to a harmonic potential. In the network each node is the Cα atom of the residue and the springs represent the interactions between the nodes. The overall potential is the sum of harmonic potentials between interacting nodes. To describe the internal motions of the spring connecting the two atoms, there is only one [[degree of freedom]]. Qualitatively, this corresponds to the compression and expansion of the spring in a direction given by the locations of the two atoms. In other words, ANM is an extension of the Gaussian Network Model to three coordinates per atom, thus accounting for directionality.

The network includes all interactions within a cutoff distance, which is the only predetermined parameter in the model. Information about the orientation of each interaction with respect to the global coordinates system is considered within the Forceforce constant matrix (''H'') and allows prediction of anisotropic motions. Consider a sub-system consisting of nodes ''i'' and ''j'', let ''r''_''i'' = (''x''_''i'' ''y''_''i'' ''z''_''i'') and let ''r''_''j'' = (''x''_''j'' ''y''_''j'' ''z''_''j'') be the instantaneous positions of atoms ''i'' and ''j''. The equilibrium distance between the atoms is represented by ''s''_''ij''^O and the instantaneous distance is given by ''s''_''ij''. For the spring between ''i'' and ''j'', the harmonic potential in terms of the unknown spring constant ''γ'', is given by:<br />

<gallery>

File:elastic network model.png|Figure 1: Elastic Network Model

</gallery>

: <math>V_{ij} = {\gamma\over 2}{( s_{ij} - {s_{ij}}^o)}^2</math>

The network includes all interactions within a cutoff distance, which is the only predetermined parameter in the model. Information about the orientation of each interaction with respect to the global coordinates system is considered within the Force constant matrix (H) and allows prediction of anisotropic motions. Consider a sub-system consisting of nodes i and j, let r_i = (x_i y_i z_i) and let r_j = (x_j y_j z_j) be the instantaneous positions of atoms i and j. The equilibrium distance between the atoms is represented by s_ij^O and the instantaneous distance is given by s_ij. For the spring between i and j, the harmonic potential in terms of the unknown spring constant γ, is given by:<br />

The second derivatives of the potential, ''V''_''ij'' with respect to the components of ''r''_''i'' are evaluated at the equilibrium position, i.e. ''s''_''ij''^O = &nbsp;''s''_''ij'', are

<math>V_{ij} = {\gamma\over 2}{( s_{ij} - {s_{ij}}^o)}^2</math>

: <math>{\partial^2 V_{ij}\over\partial x_i^2} = {\partial^2 y_jV_{ij}\over\partial x_j^2} = {-\gamma\over {s_{ij}}^2} {(x_j - x_i)}{(y_j-y_i)}^2 </math> <br />

The second derivatives of the potential, V_ij with respect to the components of r_i are evaluated at the equilibrium position, i.e. s_ij^O = s_ij, are

: <math>{\partial^2 V_{ij}\over\partial{ x_i}^2 \, \partial y_j} = {-\gamma\over {s_{ij}}^2} {(x_j - x_i)}^2{(y_j-y_i)} </math>

The above is a direct outcome of one of the key underlying assumptions of ANM – that a given crystal structure is an energetic minimum and does not require energy minimization.

<math>{\partial^2 V_{ij}\over\partial x_i\partial y_j} = {-\gamma\over {s_{ij}}^2} {(x_j - x_i)}{(y_j-y_i)} </math> <br />

The force constant of the system can be described by the [[Hessian Matrix]] – (second partial derivative of potential V):

The force constant of the system can be described by the [[Hessian Matrixmatrix]] – (second partial derivative of potential ''V''):

<math>\Eta = \begin{bmatrix} {H_{ii}} & {H_{ij}} \\ {H_{ji}} & {H_{jj}} \end{bmatrix} </math> <br />

: <math>\Eta = \begin{bmatrix} {H_{ii}} & {H_{ij}} \\ {H_{ji}} & {H_{jj}} \end{bmatrix}. </math> <br />

Each element Hi,j is a 3×3 matrix which holds the anisotropic information regarding the orientation of nodes i,j. Each such sub matrix (or the "super element" of the Hessian) is defined as:

Each element Hi''H''_''i'',''j'' is a 3×33&nbsp;×&nbsp;3 matrix which holds the anisotropic information regarding the orientation of nodes ''i'',''j''. Each such sub matrix (or the "super element" of the Hessian) is defined as:

<math>H_{ij} = \begin{bmatrix} {\partial^2 V_{ij}\over\partial x_i\partial x_j} & {\partial^2 V_{ij}\over\partial x_i\partial y_j} & {\partial^2 V_{ij}\over\partial x_i\partial z_j} \\ {\partial^2 V_{ij}\over\partial y_i\partial x_j} & {\partial^2 V_{ij}\over\partial y_i\partial y_j} & {\partial^2 V_{ij}\over\partial y_i\partial z_j} \\ {\partial^2 V_{ij}\over\partial z_i\partial x_j} & {\partial^2 V_{ij}\over\partial z_i\partial y_j} & {\partial^2 V_{ij}\over\partial z_i\partial z_j}\end{bmatrix} </math> <br / >

: <math>H_{ij} = \begin{bmatrix} {\partial^2 V_{ij}\over\partial x_i \, \partial x_j} & {\partial^2 V_{ij}\over\partial x_i \, \partial y_j} & {\partial^2 V_{ij}\over\partial x_i \, \partial z_j} \\ {\partial^2 V_{ij}\over\partial y_i \, \partial x_j} & {\partial^2 V_{ij}\over\partial y_i \, \partial y_j} & {\partial^2 V_{ij}\over\partial y_i \, \partial z_j} \\ {\partial^2 V_{ij} \over\partial z_i \, \partial x_j} & {\partial^2 V_{ij}\over\partial z_i \, \partial y_j} & {\partial^2 V_{ij}\over\partial z_i \, \partial z_j}\end{bmatrix}. </math> <br / >

Using the definition of the potential, the Hessian can be expanded as,

Using the definition of the potential, the Hessian can be expanded as,

<math>H_{ij} = {-\gamma\over {s_{ij}}^2} \begin{bmatrix} {(x_j - x_i)(x_j - x_i)} & {(x_j - x_i)(y_j - y_i)} & {(x_j - x_i)(z_j - z_i)} \\ {(y_j - y_i)(x_j - x_i)} & {(y_j - y_i)(y_j - y_i)} & {(y_j - y_i)(z_j - z_i)} \\{(z_j - z_i)(x_j - x_i)} & {(z_j - z_i)(y_j - y_i)} & {(z_j - z_i)(z_j - z_i)} \end{bmatrix}</math>

: <math>H_{ij} = {-\gamma\over {s_{ij}}^2} \begin{bmatrix} {(x_j - x_i)(x_j - x_i)} & {(x_j - x_i)(y_j - y_i)} & {(x_j - x_i)(z_j - z_i)} \\ {(y_j - y_i)(x_j - x_i)} & {(y_j - y_i)(y_j - y_i)} & {(y_j - y_i)(z_j - z_i)} \\{(z_j - z_i)(x_j - x_i)} & {(z_j - z_i)(y_j - y_i)} & {(z_j - z_i)(z_j - z_i)} \end{bmatrix}</math>

which can then be written as,

which can then be written as,

<math>H_{ij} = {-\gamma\over {s_{ij}}^2} \begin{bmatrix} x_j - x_i\\y_j - y_i\\z_j-z_i \end{bmatrix} \begin{bmatrix} x_j - x_i & y_j - y_i & z_j-z_i \end{bmatrix}</math>

: <math>H_{ij} = {-\gamma\over {s_{ij}}^2} \begin{bmatrix} x_j - x_i\\y_j - y_i\\z_j-z_i \end{bmatrix} \begin{bmatrix} x_j - x_i & y_j - y_i & z_j-z_i \end{bmatrix}</math>

Here, the force constant matrix, or the hessian matrix H holds information about the orientation of the nodes, but not about the type of the interaction (such is whether the interaction is covalent or non-covalent, hydrophobic or non-hydrophobic, etc.). In addition, the distance between the interacting nodes is not considered directly. To account for the distance between the interactions we can weight each interaction between nodes i, j by the distance, sp. The new off-diagonal elements of the Hessian matrix take the below form, where p is an empirical parameter:<br />

<math>H_{ij} = {- 1 \over {s_{ij}}^{p+2}} \begin{bmatrix} {(X_j - X_i)(X_j - X_i)} & {(X_j - X_i)(Y_j - Y_i)} & {(X_j - X_i)(Z_j - Z_i)}\\{(Y_j - Y_i)(X_j - X_i)} & {(Y_j - Y_i)(Y_j - Y_i)} & {(Y_j - Y_i)(Z_j - Z_i)}\\{(Z_j - Z_i)(X_j - X_i)} & {(Z_j - Z_i)(Y_j - Y_i)} & {(Z_j - Z_i)(Z_j - Z_i)} \end{bmatrix}</math>

Here, the force constant matrix, or the hessian matrix ''H'' holds information about the orientation of the nodes, but not about the type of the interaction (such isas whether the interaction is covalent or non-covalent, hydrophobic or non-hydrophobic, etc.). In addition, the distance between the interacting nodes is not considered directly. To account for the distance between the interactions we can weight each interaction between nodes ''i'', ''j'' by the distance, &nbsp;''sp''. The new off-diagonal elements of the Hessian matrix take the below form, where ''p'' is an empirical parameter:<br />

The counterpart of the [[Kirchhoff matrix]] Γ of the GNM is simply (1/γ) Η in the ANM. Its decomposition yields 3N - 6 non-zero [[eigenvalues]], and 3N - 6 eigenvectors that reflect the respective frequencies and shapes of the individual modes. The inverse of Η, which holds the desired information about fluctuations is composed of N x N super-elements, each of which scales with the 3 x 3 matrix of correlations between the components of pairs of fluctuation vectors. The Hessian, however is not invertible, as its rank is 3N-6 (6 variables responsible to a rigid body motion). To obtain a pseudo inverse, a solution to the eigenvalue problem is obtained:

: <math>H_{ij} = {- 1 \over {s_{ij}}^{p+2}} \begin{bmatrix} {(X_j - X_i)(X_j - X_i)} & {(X_j - X_i)(Y_j - Y_i)} & {(X_j - X_i)(Z_j - Z_i)}\\{(Y_j - Y_i)(X_j - X_i)} & {(Y_j - Y_i)(Y_j - Y_i)} & {(Y_j - Y_i)(Z_j - Z_i)}\\{(Z_j - Z_i)(X_j - X_i)} & {(Z_j - Z_i)(Y_j - Y_i)} & {(Z_j - Z_i)(Z_j - Z_i)} \end{bmatrix}</math>

<math>H = U\Lambda{U^T}</math>

The counterpart of the [[Kirchhoff matrix]] Γ of the GNM is simply (1/''γ'') &nbsp;''Η'' in the ANM. Its decomposition yields 3N - 3''N''&nbsp;−&nbsp;6 non-zero [[eigenvalues]], and 3N - 3''N''&nbsp;−&nbsp;6 eigenvectors that reflect the respective frequencies and shapes of the individual modes. The inverse of ''Η'', which holds the desired information about fluctuations is composed of ''N x ''&nbsp;×&nbsp;''N'' super-elements, each of which scales with the 3 x &nbsp;×&nbsp;3 matrix of correlations between the components of pairs of fluctuation vectors. The Hessian, however is not invertible, as its rank is 3N-6 (6 variables responsible to a rigid body motion). In other words, the eigen values corresponding to the rigid motion are 0, resulting in the determinant being 0, making the matrix not invertible. To obtain a pseudo inverse, a solution to the eigenvalue problem is obtained:

The pseudo-inverse is composed of the 3N-6 eigenvectors and their respective non-zero eigen values. Where λi are the eigenvalues of H sorted by their size from small to large and Ui the corresponding eigenvectors. The eigenvectors (the columns of the matrix U) describe the vibrational direction and the relative amplitude in the different modes.

: <math>H = U\Lambda{ U^T} </math>

The pseudo-inverse is composed of the 3N-3''N''&nbsp;−&nbsp;6 eigenvectors and their respective non-zero eigen values., Wherewhere λi''λ''_''i'' are the eigenvalues of ''H'' sorted by their size from small to large and Ui''U''_''i'' the corresponding eigenvectors. The eigenvectors (the columns of the matrix ''U'') describe the vibrational direction and the relative amplitude in the different modes.

== Comparing ANM and GNM ==

Line 48 ⟶ 49:

ANM and GNM are both based on an elastic network model. The GNM has proven itself to accurately describe the vibrational dynamics of proteins and their complexes in numerous studies. Whereas the GNM is limited to the evaluation of the [[mean squared displacement]]s and cross-correlations between fluctuations, the motion being projected to a mode space of N dimensions, the ANM approach permits us to evaluate directional preferences and thus provides 3-D descriptions of the 3N - 6 internal modes.

It has been observed that GNM fluctuation predictions agree better with experiments than those computed with ANM. The higher performance of GNM can beenbe attributed to its underlying potential, which takes account of orientational deformations, in addition to distance changes.

== Evaluation of the Modelmodel ==

ANM has been evaluated on a large set of proteins to establish the optimal model parameters that achieve the highest correlation with experimental data and its limits of accuracy and applicability. The ANM is evaluated by comparing the fluctuations predicted from theory and those experimentally observed (B-factors deposited in the PDB). During evaluation, the following observations have been made about the models behavior.

- * ANM shows insensitivity to the choice of cutoff distance within a certain range, like GNM.<br />

- * Weighting the interactions by distance improves the correlation.<br />

- ANM shows insensitivity to the choice of cutoff distance within a certain range, like GNM.<br />

- * Residue fluctuations in globular proteins are shown to be more accurately predicted, than those in non-globular proteins.<br />

- Weighting the interactions by distance improves the correlation.<br />

- * Significant improvement in agreement with experiments is observed with increase in the resolution of the examined structure.<br />

- Residue fluctuations in globular proteins are shown to be more accurately predicted, than those in non-globular proteins.<br />

- * While understanding how the accuracy of the predicted fluctuations is related to solvent accessibilities, the predictions for buried residues are shown to be in significantly better agreement with the experimental data as compared to the solvent-exposed ones.<br />

- Significant improvement in agreement with experiments is observed with increase in the resolution of the examined structure.<br />

- * Polar/charged residues are more accurately predicted than hydrophobic ones, a possible consequence of the involvement of surface hydrophobic residues in crystal contacts.

- While understanding how the accuracy of the predicted fluctuations is related to solvent accessibilities, the predictions for buried residues are shown to be in significantly better agreement with the experimental data as compared to the solvent-exposed ones.<br />

- Polar/charged residues are more accurately predicted than hydrophobic ones, a possible consequence of the involvement of surface hydrophobic residues in crystal contacts.

== Applications of ANM ==

Line 73:

- [[Membrane proteins]] including potassium channels, by Shrivastava and Bahar, 2006.<br />

- [[Rhodopsin]], by Rader et al., 2004.<br />

- [[Nicotinic acetylcholine receptor]], by Hung et al., 2005; Taly et al., 2005.<br and a few more./>

- [[Auxiliary Activity family 9]] and [[Auxiliary Activity family 10]] family of lytic polysaccharide monooxygenases by Arora et al.,2019 [https://www.sciencedirect.com/science/article/pii/S1093326318306776?via%3Dihub] and a few more.

== ANM Webweb Serversservers ==

The ANM web server developed by Eyal E, Yang LW, Bahar I. in 2006, presents a web-based interface for performing ANM calculations, the main strengths of which are the rapid computing ability and the user-friendly graphical capabilities for analyzing and interpreting the outputs.<br />

-* Anisotropic Network Model web server. [http://ignmtest.ccbb.pitt.edu/cgi-bin/anm/anm1.cgi]<br />

-* ANM server. [https://web.archive.org/web/20070624170153/http://gor.bb.iastate.edu/anm/anm.htm]

== References ==

1. # "Anisotropy of fluctuation dynamics of proteins with an elastic network model", A.R. Atilgan et al., ''Biophys. J.'' 80, 505 (2001).<br />

2. # "Anisotropic network model: systematic evaluation and a new web interface", Eyal E, Yang LW, Bahar I. ''Bioinformatics''. 22, 2619–2627, (2006)<br />

1. Anisotropy of fluctuation dynamics of proteins with an elastic network model, A.R. Atilgan et al., Biophys. J. 80, 505 (2001).<br />

3. # "Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor", Doruker, P, Atilgan, AR & Bahar, I, ''Proteins'', 15, 512-524512–524, ( 2000).<br />

2. Anisotropic network model: systematic evaluation and a new web interface, Eyal E, Yang LW, Bahar I. Bioinformatics. 22, 2619–2627, (2006)<br />

4. # Hinsen, K. (1998) "Analysis of ___domain motions by approximate normal mode calculations", ''Proteins'', 33, 417-429417–429. {{PMID |11159421<br />}}

3. Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor, Doruker, P, Atilgan, AR & Bahar, I, Proteins, 15, 512-524, ( 2000).<br />

5. # Bahar,I. et al. (1997) "Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential". ''Fold Des'', 2, 173-181<br />173–181

4. Hinsen, K. (1998) Analysis of ___domain motions by approximate normal mode calculations, Proteins, 33, 417-429. PMID 11159421<br />

6. # Chennubhotla,C. et al. (2005) "Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies". ''Phys Biol'', 2, S173-S180pp.&nbsp;173–180.<br />

5. Bahar,I. et al. (1997) Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. Fold Des, 2, 173-181<br />

7. # Cui,Q. and Bahar,I. (2006) ''Normal Mode Analysis: Theory and Applications to Biological and Chemical Systems''. Chapman & Hall/CRC, Boca Raton, FL.

6. Chennubhotla,C. et al. (2005) Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies. Phys Biol, 2, S173-S180.<br />

# Arora et al. (2019) "Structural dynamics of lytic polysaccharide monoxygenases reveals a highly flexible substrate binding region". J ''Mol Graph Model'', 88, 1–10. [https://www.sciencedirect.com/science/article/pii/S1093326318306776?via%3Dihub]

7. Cui,Q. and Bahar,I. (2006) Normal Mode Analysis: Theory and Applications to Biological and Chemical Systems. Chapman & Hall/CRC, Boca Raton, FL.

== See also ==

- Gaussian Network Model* [[Gaussian network model]]

- Gaussian Network Model [[Gaussian network model]]

[[Category:Molecular modelling]]