Structural stabilization of a rigid beta-sheet cluster of fucosylated proteinase inhibitor PMPC (Pars intercerebralis major peptide C) against thermal denaturation: An unfolding molecular dynamics simulation study.

Unfolding behavior of glycosylated- and unglycosylated proteinase inhibitor Pars intercerebralis major peptide C (PMPC) at 350 K were traced with molecular dynamics simulations using the CHARMM program. The fucosylated PMPC (FPMPC) possesses a nearly identical protein structure with PMPC, differing only by the presence of a single fucose residue linked to Thr9 in the PMPC. Attachment of a monomeric fucose residue to the Thr9 in PMPC resulted in a change of the denaturing process of FPMPC. Simulations showed that the unfolding of PMPC involved significant weakening of non-local interactions whereas fucosylation led FPMPC to preserve the non-local interactions, even in its denatured form. Even in simulations over 16 ns at 350 K, FPMPC remained relatively stable in a less denatured conformation. However, the conformation of PMPC transformed to a fully unfolded state within 5 ns in the simulation at 350 K. This difference was due to the formation of fucose-mediated hydrogen bonds and non-local contacts by the attached fucose residue of FPMPC. In the case of FPMPC, fucosyl residue was involved in maintaining a rigid beta-sheet cluster through interaction with the hydrogen bond network. These high-temperature unfolding MD simulations provide a theoretical basis for a previous experimental work in which FPMPC showed stable unfolding thermodynamics compared to unfucosylated PMPC, suggesting that single fucosylation induces conformational stabilization of PMPC by tertiary contacts.

[1]  S. N. Timasheff,et al.  Mechanism of the stabilization of ribonuclease a by sorbitol: Preferential hydration is greater for the denatured than for the native protein , 1997, Protein science : a publication of the Protein Society.

[2]  Leslie A Kuhn,et al.  Protein unfolding: Rigidity lost , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[3]  J. Lefèvre,et al.  Stabilization of proteins by glycosylation examined by NMR analysis of a fucosylated proteinase inhibitor , 1996, Nature Structural Biology.

[4]  C. Turano,et al.  Influence of the carbohydrate moiety on the stability of glycoproteins. , 1996, Biochemistry.

[5]  V. Daggett,et al.  Increasing temperature accelerates protein unfolding without changing the pathway of unfolding. , 2002, Journal of molecular biology.

[6]  Seunho Jung,et al.  Retardation of the unfolding process by single N‐glycosylation of ribonuclease A based on molecular dynamics simulations , 2008, Biopolymers.

[7]  J. Lefèvre,et al.  Solution structure of PMP-C: a new fold in the group of small serine proteinase inhibitors. , 1996, Journal of molecular biology.

[8]  V. Eijsink,et al.  Rational engineering of enzyme stability. , 2004, Journal of biotechnology.

[9]  H. Gohlke,et al.  Exploiting the Link between Protein Rigidity and Thermostability for Data‐Driven Protein Engineering , 2008 .

[10]  S. Elliott,et al.  Glycoengineering: the effect of glycosylation on the properties of therapeutic proteins. , 2005, Journal of pharmaceutical sciences.

[11]  S. Jackson,et al.  The folding pathway of ubiquitin from all-atom molecular dynamics simulations. , 2004, Biophysical chemistry.

[12]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[13]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[14]  S. Flitsch,et al.  Increased Thermal Stability of Site‐Selectively Glycosylated Dihydrofolate Reductase , 2005, Chembiochem : a European journal of chemical biology.

[15]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[16]  M. Kuttel,et al.  Water structure about the dimer and hexamer repeat units of amylose from molecular dynamics computer simulations , 2001, J. Comput. Chem..

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

[18]  P. Hünenberger,et al.  Trehalose–protein interaction in aqueous solution , 2004, Proteins.

[19]  Wei Gu,et al.  Molecular dynamics simulation of the unfolding of the human prion protein domain under low pH and high temperature conditions. , 2003, Biophysical chemistry.

[20]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .