Skip to content


Open Access

Hybrid approaches for the detection of networks of critical residues involved in functional motions in protein families

BMC Bioinformatics201516(Suppl 3):A8

Published: 13 February 2015


Protein FamilyHybrid ApproachCombinatorial LibraryAnalysis SimulationProtein Dynamic


Currently there is great interest in identifying critical residues in proteins, to improve our understanding and allow for the engineering of protein families. Diverse approaches combine sequence information, structural data, dynamics analysis and functional description to determine the importance of amino acids with regards to protein function. In this work, we propose a hybrid approach for the identification of critical residues in proteins, combining the use of evolutionary information (co-evolution), cross-correlation of atomic fluctuations derived from Anisotropic Normal Mode Analysis simulations [1] (ANMA) and network analysis. Subsequently we have compared this method to existing approaches.


By combining the information of the covariance matrix derived from Statistical Coupling Analysis (SCA) [2] and the cross-correlation matrix of atomic fluctuations derived from ANMA, it was possible to identify a network of evolutionarily coupled residues involved in relevant motions in protein families. The outstanding sites revealed by our hybrid approach (ANMA.SCA) showed a high correspondence with experimental data, confirming the critical role of these sites in the functional mobility of proteins. In addition, our approach was found to be complementary to previous approaches. It maintained a good correspondence with approaches derived from extensive molecular dynamics, while being faster and less expensive in terms of computational resources [3].


The hybrid approach ANMA.SCA opens a wide range of possibilities in the study of functional motion within protein families. By means of detecting networks of critical sites and their topology it is able to reveal the hidden aspects of protein dynamics.



DA-M acknowledges the PhD fellowship (35083) from CONACYT and (IN-204714) DGAPA. EP-R was supported by a grant (IN-204714) from DGAPA and (155116) from CONACYT.

Authors’ Affiliations

Departamento de Ingeniería Celular y Biocatálisis, Instituto de Biotecnología, UNAM Av. Universidad 2001, Cuernavaca, Morelos, México
Unidad Multidisciplinaria de Docencia e Investigación, Sisal Facultad de Ciencias, UNAM, Sisal, Yucatán, México


  1. Eyal E, Yang LW, Bahar I: Anisotropic network model: systematic evaluation and a new web interface. Bioinformatics. 2006, 22 (21): 2619-2627. 10.1093/bioinformatics/btl448.View ArticlePubMedGoogle Scholar
  2. Socolich M, Lockless SW, Russ WP, Lee H, Gardner KH, Ranganathan R: Evolutionary information for specifying a protein fold. Nature. 2005, 437 (7058): 512-518. 10.1038/nature03991.View ArticlePubMedGoogle Scholar
  3. Armenta-Medina D, Perez-Rueda E, Segovia L: Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches. Proteins. 2011, 79 (5): 1662-1671. 10.1002/prot.22995.View ArticlePubMedGoogle Scholar


© Armenta-Medina and Perez-Rueda; licensee BioMed Central Ltd. 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.