Comparative structural studies of psychrophilic and mesophilic protein homologues by molecular dynamics simulation.

Comparative molecular dynamics simulations of psychrophilic type III antifreeze protein from the North-Atlantic ocean-pout Macrozoarces americanus and its corresponding mesophilic counterpart, the antifreeze-like domain of human sialic acid synthase, have been performed for 10 ns each at five different temperatures. Analyses of trajectories in terms of secondary structure content, solvent accessibility, intramolecular hydrogen bonds and protein-solvent interactions indicate distinct differences in these two proteins. The two proteins also follow dissimilar unfolding pathways. The overall flexibility calculated by the trace of the diagonalized covariance matrix displays similar flexibility of both the proteins near their growth temperatures. However at higher temperatures psychrophilic protein shows increased overall flexibility than its mesophilic counterpart. Principal component analysis also indicates that the essential subspaces explored by the simulations of two proteins at different temperatures are non-overlapping and they show significantly different directions of motion. However, there are significant overlaps within the trajectories and similar directions of motion of each protein especially at 298 K, 310 K and 373 K. Overall, the psychrophilic protein leads to increased conformational sampling of the phase space than its mesophilic counterpart. Our study may help in elucidating the molecular basis of thermostability of homologous proteins from two organisms living at different temperature conditions. Such an understanding is required for designing efficient proteins with characteristics for a particular application at desired working temperatures.

[1]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[2]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[3]  Elena Papaleo,et al.  Flexibility and enzymatic cold-adaptation: a comparative molecular dynamics investigation of the elastase family. , 2006, Biochimica et biophysica acta.

[4]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[5]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[6]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[7]  R A Sayle,et al.  RASMOL: biomolecular graphics for all. , 1995, Trends in biochemical sciences.

[8]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[9]  B. Brandsdal,et al.  Comparative unfolding studies of psychrophilic and mesophilic uracil DNA glycosylase: MD simulations show reduced thermal stability of the cold-adapted enzyme. , 2007, Journal of molecular graphics & modelling.

[10]  B. Brandsdal,et al.  Comparative molecular dynamics of mesophilic and psychrophilic protein homologues studied by 1.2 ns simulations. , 1999, Journal of biomolecular structure & dynamics.

[11]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[12]  Z. Jia,et al.  Quantitative and Qualitative Analysis of Type III Antifreeze Protein Structure and Function* , 1999, The Journal of Biological Chemistry.

[13]  M. G. Mateu,et al.  Thermodynamic stability of a cold‐adapted protein, type III antifreeze protein, and energetic contribution of salt bridges , 2006, Protein science : a publication of the Protein Society.

[14]  Z. Jia,et al.  Ice-binding surface of fish type III antifreeze. , 1999, Biophysical journal.

[15]  H. Berendsen,et al.  A comparison of techniques for calculating protein essential dynamics , 1997 .

[16]  B. Hess,et al.  Similarities between principal components of protein dynamics and random diffusion , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[18]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[19]  A. Amadei,et al.  On the convergence of the conformational coordinates basis set obtained by the essential dynamics analysis of proteins' molecular dynamics simulations , 1999, Proteins.

[20]  B. Hess Convergence of sampling in protein simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  H. Berendsen,et al.  The consistency of large concerted motions in proteins in molecular dynamics simulations. , 1996, Biophysical journal.

[22]  B. Brandsdal,et al.  Increased Flexibility as a Strategy for Cold Adaptation , 2005, Journal of Biological Chemistry.

[23]  E. Papaleo,et al.  Optimization of electrostatics as a strategy for cold-adaptation: a case study of cold- and warm-active elastases. , 2007, Journal of molecular graphics & modelling.

[24]  Elena Papaleo,et al.  Protein flexibility in psychrophilic and mesophilic trypsins. Evidence of evolutionary conservation of protein dynamics in trypsin‐like serine‐proteases , 2008, FEBS letters.

[25]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[26]  Temperature-induced unfolding pathway of a type III antifreeze protein: insight from molecular dynamics simulation. , 2008, Journal of molecular graphics & modelling.

[27]  Yoko Ito,et al.  Solution structure of the antifreeze‐like domain of human sialic acid synthase , 2006, Protein science : a publication of the Protein Society.

[28]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .