Integrating protein structural dynamics and evolutionary analysis with Bio3D

BMC Bioinformatics

Lars Skjærven¹  Xin-Qiu Yao²  Guido Scarabelli²  Barry J Grant² 

[1] Structural and Computational Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany

[2] Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA

关键词: Evolution;    Principal component analysis;    Molecular dynamics;    Normal mode analysis;    Allostery;    Protein dynamics;    Protein structure;   

Others  :  1084471
DOI  :  10.1186/s12859-014-0399-6

 received in 2014-10-07, accepted in 2014-11-26,  发布年份 2014

【 摘 要 】

Background

Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution.

Results

Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case.

Conclusions

The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/ webcite.

【 授权许可】

   
2014 Skjærven et al.; licensee BioMed Central.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

【 参考文献 】

• [1]Teilum K, Olsen JG, Kragelund BB: Functional aspects of protein flexibility. Cell Mol Life Sci 2009, 66:2231-2247.
• [2]Henzler-Wildman K, Kern D: Dynamic personalities of proteins. Nature 2007, 450:964-972.
• [3]Henzler-Wildman KA, Thai V, Lei M, Ott M, Wolf-Watz M, Fenn T, Pozharski E, Wilson MA, Petsko GA, Karplus M, Hübner CG, Kern D: Intrinsic motions along an enzymatic reaction trajectory. Nature 2007, 450:838-844.
• [4]Vale RD, Milligan RA: The way things move: looking under the hood of molecular motor proteins. Science 2000, 288:88-95.
• [5]Yébenes H, Mesa P, Muñoz IG, Montoya G, Valpuesta JM: Chaperonins: two rings for folding. Trends Biochem Sci 2011, 36:424-432.
• [6]Smock RG, Gierasch LM: Sending signals dynamically. Science 2009, 324:198-203.
• [7]Karplus M, McCammon JA: Molecular dynamics simulations of biomolecules. Nat Struct Biol 2002, 9:646-652.
• [8]Lee EH, Hsin J, Sotomayor M, Comellas G, Schulten K: Discovery through the computational microscope. Structure 2009, 17:1295-1306.
• [9]Skjaerven L, Hollup SM, Reuter N: Normal mode analysis for proteins. J Mol Struct (THEOCHEM) 2009, 898:42-48.
• [10]Suhre K, Sanejouand Y-H: EN: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement. Nucleic Acids Res 2004, 32:W610-W614.
• [11]Krüger DM, Ahmed A, Gohlke H: NMSim web server: integrated approach for normal mode-based geometric simulations of biologically relevant conformational transitions in proteins. Nucleic Acids Res 2012, 40:W310-W316.
• [12]Eyal E, Yang L-W, Bahar I: Anisotropic network model: systematic evaluation and a new web interface. Bioinformatics 2006, 22:2619-2627.
• [13]Hinsen K: The molecular modeling toolkit: a new approach to molecular simulations. J Comput Chem 2000, 21:79-85.
• [14]Bakan A, Meireles LM, Bahar I: ProDy: protein dynamics inferred from theory and experiments. Bioinformatics 2011, 27:1575-1577.
• [15]Zimmermann MT, Kloczkowski A, Jernigan RL: MAVENs: motion analysis and visualization of elastic networks and structural ensembles. BMC Bioinformatics 2011, 12:264. BioMed Central Full Text
• [16]Grant B, Rodrigues A, ElSawy KM, McCammon JA, Caves LSD: Bio3d: an R package for the comparative analysis of protein structures. Bioinformatics 2006, 22:2695-2696.
• [17]Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215:403-410.
• [18]Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res 2011, 39:W29-W37.
• [19]Atilgan AR, Durell 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.
• [20]Yang L, Song G, Jernigan RL: Protein elastic network models and the ranges of cooperativity. Proc Natl Acad Sci U S A 2009, 106:12347-12352.
• [21]Hinsen K, Petrescu AJ, Dellerue S, Bellissent-Funel MC, Kneller GR: Harmonicity in slow protein dynamics. Chem Phys 2000, 261:25-37.
• [22]Moritsugu K, Smith JC: Coarse-grained biomolecular simulation with REACH: realistic extension algorithm via covariance Hessian. Biophys J 2007, 93:3460-3469.
• [23]Dehouck Y, Mikhailov AS: Effective harmonic potentials: insights into the internal cooperativity and sequence-specificity of protein dynamics. PLoS Comput Biol 2013, 9:e1003209.
• [24]Tirion M: Large amplitude elastic motions in proteins from a single-parameter, atomic analysis. Phys Rev Lett 1996, 77:1905-1908.
• [25]Hinsen K: Analysis of domain motions by approximate normal mode calculations. Proteins 1998, 33:417-429.
• [26]Fuglebakk E, Echave J, Reuter N: Measuring and comparing structural fluctuation patterns in large protein datasets. Bioinformatics 2012, 28:2431-2440.
• [27]Caves LSD, Evanseck JD, Karplus M: Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin. Protein Sci 1998, 7:649-666.
• [28]Gorfe AA, Grant BJ, McCammon JA: Mapping the nucleotide and isoform-dependent structural and dynamical features of ras proteins. Structure 2008, 16:885-896.
• [29]Grant BJ, McCammon JA, Caves LSD, Cross RA: Multivariate analysis of conserved sequence-structure relationships in kinesins: coupling of the active site and a tubulin-binding sub-domain. J Mol Biol 2007, 368:1231-1248.
• [30]Van Aalten DMF, de Groot BL, Findlay JBC, Berendsen HJC, Amadei A, VanAalten DMF, DeGroot BL: A comparison of techniques for calculating protein essential dynamics. J Comput Chem 1997, 18:169-181.
• [31]Amadei A, Ceruso MA, Di Nola A: On the convergence of the conformational coordinates basis set obtained by the essential dynamics analysis of proteins’ molecular dynamics simulations. Proteins 1999, 36:419-424.
• [32]Romo TD, Grossfield A: Validating and improving elastic network models with molecular dynamics simulations. Proteins 2011, 79:23-34.
• [33]Fuglebakk E, Reuter N, Hinsen K: Evaluation of protein elastic network models based on an analysis of collective motions. J Chem Theory Comput 2013, 9:5618-5628.
• [34]Kundu S, Melton JS, Sorensen DC, Phillips GN: Dynamics of proteins in crystals: comparison of experiment with simple models. Biophys J 2002, 83:723-732.
• [35]Romanowska J, Nowinski KS, Trylska J: Determining geometrically stable domains in molecular conformation sets. J Chem Theory Comput 2012, 8:2588-2599.
• [36]Ichiye T, Karplus M: Collective motions in proteins: a covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations. Proteins 1991, 11:205-217.
• [37]Lange OF, Grubmüller H: Generalized correlation for biomolecular dynamics. Proteins 2006, 62:1053-1061.
• [38]Sethi A, Eargle J, Black AA, Luthey-Schulten Z: Dynamical networks in tRNA:protein complexes. Proc Natl Acad Sci U S A 2009, 106:6620-6625.
• [39]Girvan M, Newman MEJ: Community structure in social and biological networks. Proc Natl Acad Sci U S A 2002, 99:7821-7826.
• [40]Schnell JR, Dyson HJ, Wright PE: Structure, dynamics, and catalytic function of dihydrofolate reductase. Annu Rev Biophys Biomol Struct 2004, 33:119-140.
• [41]Bhabha G, Ekiert DC, Jennewein M, Zmasek CM, Tuttle LM, Kroon G, Dyson HJ, Godzik A, Wilson I a, Wright PE: Divergent evolution of protein conformational dynamics in dihydrofolate reductase. Nat Struct Mol Biol 2013, 20:1243-1249.
• [42]Yao X-Q, Grant BJ: Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins. Biophys J 2013, 105:L08-L10.
• [43]Rasmussen SGF, DeVree BT, Zou Y, Kruse AC, Chung KY, Kobilka TS, Thian FS, Chae PS, Pardon E, Calinski D, Mathiesen JM, Shah STA, Lyons J a, Caffrey M, Gellman SH, Steyaert J, Skiniotis G, Weis WI, Sunahara RK, Kobilka BK: Crystal structure of the β2 adrenergic receptor-Gs protein complex. Nature 2011, 477:549-555.
• [44]Westfield GH, Rasmussen SGF, Su M, Dutta S, DeVree BT, Chung KY, Calinski D, Velez-Ruiz G, Oleskie AN, Pardon E, Chae PS, Liu T, Li S, Woods VL, Steyaert J, Kobilka BK, Sunahara RK, Skiniotis G: Structural flexibility of the G alpha s alpha-helical domain in the beta2-adrenoceptor Gs complex. Proc Natl Acad Sci U S A 2011, 108:16086-16091.
• [45]Hollup SM, Salensminde G, Reuter N: WEBnm@: a web application for normal mode analyses of proteins. BMC Bioinformatics 2005, 6:52. BioMed Central Full Text
• [46]Roe DR, Cheatham TE: PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J Chem Theory Comput 2013, 9:3084-3095.
• [47]Salomon-Ferrer R, Case D a, Walker RC: An overview of the Amber biomolecular simulation package. WIREs Comput Mol Sci 2013, 3:198-210.
• [48]Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5:R80. BioMed Central Full Text