Abstract
Finite element (FE) modeling approach has emerged as an efficient way of calculating the dynamic properties of supramolecular protein structures and their complexes. Its efficiency mainly stems from the fact that the complexity of three-dimensional shape of a molecular surface dominates the computational cost rather than the molecular size or the number of atoms. However, no critical evaluation of the method has been made yet particularly for its sensitivity to the parameters used in model construction. Here, we make a close investigation on the effect of FE model parameters by analyzing 135 representative protein structures whose normal modes calculated using all-atom normal mode analysis are publicly accessible online. Results demonstrate that it is more beneficial to use a contour surface of electron densities as the molecular surface, in general, rather than to employ a solvent excluded surface, and that the solution accuracy is almost insensitive to the model parameters unless we avoid extreme values leading to an inaccurate depiction of the characteristic shapes.













Similar content being viewed by others
References
McCammon JA, Gelin BR, Karplus M (1977) Dynamics of folded proteins. Nature 267:585–590
García AE, Harman JG (1996) Simulations of CRP:(cAMP)2 in noncrystalline environments show a subunit transition from the open to the closed conformation. Protein Sci 5:62–71
de Groot BL, Hayward S, van Aalten DM, Amadei A, Berendsen, HJC (1998) Domain motions in bacteriophage T4 lysozyme: a comparison between molecular dynamics and crystallographic data. Proteins 31:116–127
Karplus M, McCammon JA (2002) Molecular dynamics simulations of biomolecules. Nat Struct Biol 9:646–652
Cui Q, Bahar I (2005) Normal Mode Analysis: Theory and Applications to Biological and Chemical Systems. Chapman and Hall/CRC, Boca Raton, FL
Brooks B, Karplus M (1983) Harmonic dynamics of proteins: normal modes and fluctuations in bovine pancreatic trypsin inhibitor. Proc Natl Acad Sci USA 80:6571–6575
Levitt M, Sander C, Stern PS (1983) The normal modes of a protein: Native bovine pancreatic trypsin inhibitor. Int J Quant Chem Quant Biol Symp 10:181–199
Kitao A, Go N (1999) Investigating protein dynamics in collective coordinate space. Curr Opin Struct Biol 9:64–169
Tozzini V (2005) Coarse-grained models for proteins. Curr Opin Struct Biol 15:144–150
Tirion MM (1996) Large amplitude elastic motions in proteins from a single-parameter, atomic analysis. Phys Rev Lett 77:1905–1908
Ueda Y, Taketomi H, Gō N (1978) Studies on protein folding, unfolding, and fluctuations by computer simulation. II. A three-dimensional lattice model of lysozyme. Biopolymers 17:1531–1548
Bahar I, Atilgan AR, Erman B (1997) Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. Fold. Des. 2:173–181
Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I (2001) Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J 80:505–515
Kim MH, Kim MK (2014) Review: elastic network model for protein structural dynamics. JSM Enzymol. Protein Sci 1:1001
Na H, Song G (2014) Bridging between normal mode analysis and elastic network models. Proteins 82:1–12
Kim MK, Jernigan RL, Chirikjian GS (2002) Efficient generation of feasible pathways for protein cconformational transitions. Biophys J 83:1620–1630
Tama F, Miyashita O, Brooks CL III (2004) Normal mode based flexible fitting of high-resolution structure into low-resolution experimentlal data from cryo-EM. J Struct Biol 147:315–326
Bathe M (2008) A finite element framework for computation of protein normal modes and mechanical response. Proteins 70:1595–1609
Kim DN, Altschuler J, Strong C, McGill G, Bathe M (2011) Conformational dynamics data bank: a database for conformational dynamics of proteins and supramolecular protein assemblies. Nucleic Acids Res 39:451–455
Kim DN, Nguyen CT, Bathe M (2011) Conformational dynamics of supramolecular protein assemblies. J Struct Biol 173:261–270
Kim J, Kim JG, Yun G, Lee PS, Kim DN (2015) Towards modular analysis of supramolecular protein assemblies. J Chem Theory Comput 11:4260–4272
Gibbons MM, Klug WS (2007) Nonlinear finite-element analysis of nanoindentation of viral capsids. Phys. Rev. E 75:1–11
Sakalli I, Schöberl J, Knapp EW (2014) mFES: a robust molecular finite element solver for electrostatic energy compuations. J Chem Theory Comput 10:5095–5112
Na H, Song G (2015) Conventional NMA as a better standard for evaluating elastic network models. Proteins 83:259–267
Na H, Song G (2015) The performance of fine-grained and coarse-grained elastic network models and its dependence on various factors. Proteins 83:1273–1283
Connolly ML (1983) Analytical molecular surface calculation. J Appl Cryst 16:548–558
Sanner MF, Olson AJ, Spehner JC (1996) Reduced surface: an efficient way to compute molecular surfaces. Biopolymers 38:305–320
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE (2004) UCSF Chimera–A visualization system for exploratory research and analysis. J Comput Chem 25:1605–1612
Cignoni P, Callieri M, Corsini M, Dellepiane M, Ganovelli F, Ranzuglia G (2008) Meshlab: an open-source mesh processing tool. In Eurographics Italian Chapter Conference. Dblp, Salerno
Löhner R, Parikh P (1988) Generation of three-dimensional unstructured grids by the advancing-front method. Int J Numer Methods Fluids 8:1135–1149
Peraire J, Peiro J, Formaggia L, Morgan K, Zienkiewicz OC (1988) Finite element Euler computations in three dimensions. Int J Numer Methods Eng 26:2135–2159
Watson DF (1981) Computing the n-dimensional Delaunay tessellation with application to Voronoi polytopes. Comput J 24:167–172
Shenton D, Cendes Z (1985) Three-dimensional finite element mesh generation using Delaunay tesselation. IEEE Trans Magn 21:2535–2538
Fischer H, Polikarpov I, Craievich AF (2004) Average protein density is a molecular-weight-dependent function. Protein Sci 13:2825–2828
Sedeh RS, Bathe M, Bathe KJ (2010) The subspace iteration method in protein normal mode analysis. J Comput Chem 31:66–74
Garland M, Heckbert PS (1997). Surface simplification using quadric error metrics. In Proceedings of the 24th annual conference on Computer graphics and interactive techniques. ACM Press, New York.
Lorensen WE, Cline HE (1987) Marching cubes: a high resolution 3D surface construction algorithm. ACM Siggraph Comp Gr 21:163–169
Kondrashov DA, Van Wynsberghe AW, Bannen RM, Cui Q, Phillips GN (2007) Protein structural variation in computational models and crystallographic data. Structure 15:169–177
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Funding
This research was supported by the EDucation-research Integration through Simulation On the Net (EDISON) Program (Grant No. 2014M3C1A6038842), the Global Frontier R&D Program on Center for Wave Energy Control based on Metamaterials (Grant No. 2014M3A6B3063711), and the Basic Science Research Program (Grant No. 2016R1C1B2011098) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning.
Author contributions
The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Yun, G., Kim, J. & Kim, DN. A critical assessment of finite element modeling approach for protein dynamics. J Comput Aided Mol Des 31, 609–624 (2017). https://doi.org/10.1007/s10822-017-0027-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10822-017-0027-4