Integrity of H1 helix in prion protein revealed by molecular dynamic simulations to be especially vulnerable to changes in the relative orientation of H1 and its S1 flank

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

The cellular and pathologic prion protein

The cellular prion protein, PrP^C, is a small, cell surface glycoprotein with a function that is currently somewhat ill defined. It is also the key molecule involved in the family of neurodegenerative disorders called transmissible spongiform encephalopathies, which are also known as prion diseases. The misfolding of PrP^C to a conformationally altered isoform, designated PrP^TSE, is the main molecular process involved in pathogenesis and appears to precede many other pathologic and clinical manifestations of disease, including neuronal loss, astrogliosis, and cognitive loss. PrP^TSE is also believed to be the major component of the infectious "prion," the agent responsible for disease transmission, and preparations of this protein can cause prion disease when inoculated into a naïve host. Thus, understanding the biochemical and biophysical properties of both PrP^C and PrP^TSE, and ultimately the mechanisms of their interconversion, is critical if we are to understand prion disease biology. Although entire books could be devoted to research pertaining to the protein, herein we briefly review the state of knowledge of prion biochemistry, including consideration of prion protein structure, function, misfolding, and dysfunction.

Molecular dynamics studies on the buffalo prion protein

It was reported that buffalo is a low susceptibility species resisting to transmissible spongiform encephalopathies (TSEs) (same as rabbits, horses, and dogs). TSEs, also called prion diseases, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of species (except for rabbits, dogs, horses, and buffalo), manifesting as scrapie in sheep and goats; bovine spongiform encephalopathy (BSE or "mad-cow" disease) in cattle; chronic wasting disease in deer and elk; and Creutzfeldt-Jakob diseases, Gerstmann-Sträussler-Scheinker syndrome, fatal familial insomnia, and Kulu in humans etc. In molecular structures, these neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP(C)), predominantly with α-helices, into insoluble abnormally folded infectious prions (PrP(Sc)), rich in β-sheets. In this article, we studied the molecular structure and structural dynamics of buffalo PrP(C) (BufPrP(C)), in order to understand the reason why buffalo is resistant to prion diseases. We first did molecular modeling of a homology structure constructed by one mutation at residue 143 from the NMR structure of bovine and cattle PrP(124-227); immediately we found that for BufPrP(C)(124-227), there are five hydrogen bonds (HBs) at Asn143, but at this position, bovine/cattle do not have such HBs. Same as that of rabbits, dogs, or horses, our molecular dynamics studies also revealed there is a strong salt bridge (SB) ASP178-ARG164 (O-N) keeping the β2-α2 loop linked in buffalo. We also found there is a very strong HB SER170-TYR218 linking this loop with the C-terminal end of α-helix H3. Other information, such as (i) there is a very strong SB HIS187-ARG156 (N-O) linking α-helices H2 and H1 (if mutation H187R is made at position 187, then the hydrophobic core of PrP(C) will be exposed (L.H. Zhong (2010). Exposure of hydrophobic core in human prion protein pathogenic mutant H187R. Journal of Biomolecular Structure and Dynamics 28(3), 355-361)), (ii) at D178, there is a HB Y169-D178 and a polar contact R164-D178 for BufPrP(C) instead of a polar contact Q168-D178 for bovine PrP(C) (C.J. Cheng, & V. Daggett. (2014). Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH. Biomolecules 4(1), 181-201), (iii) BufPrP(C) owns three 310 helices at 125-127, 152-156, and in the β2-α2 loop, respectively, and (iv) in the β2-α2 loop, there is a strong π-π stacking and a strong π-cation F175-Y169-R164.(N)NH2, has been discovered.