Molecular dynamics simulation of the phosphorylation‐induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho‐protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau_225?250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β‐sheet pattern (²²⁶VAVVR²³⁰ and ²⁴⁴QTAPVP²⁴⁹), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α‐helix (²³⁹AKSRLQT²⁴⁵). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.     

Molecular dynamics simulations to investigate the stability and aggregation behaviour of the amyloid-forming peptide VQIVYK from tau protein

The formation of paired helical filaments arising from the short hexapeptide in the third repeat of tau protein, ³⁰⁶VQIVYK³¹¹, is critical for tau polymerisation. The atomic structure of the VQIVYK oligomer has revealed a dry, tightly self-complementing structure between the neighbouring β-sheet layers, termed as ‘steric zipper’. In this study, several molecular dynamics simulations with all-atom explicit water were conducted to investigate the structural stability and aggregation behaviour of the VQIVYK peptide with various sizes and its single alanine replacement mutations. Our results indicate that the van der Waals interaction between side chains of Q2, the π–π stacking interaction between aromatic rings of Y5, and the electrostatic interaction between K6 and the C-terminus play an important role in stabilising the VQIVYK oligomers within the same β-sheet layer, while hydrophobic steric zipper involving V1, I3 and Y5 is responsible for holding the neighbouring β-sheet layers together. The twisted angles of the VQIVYK oligomers were also analysed and shown to be size dependent. The present results not only provide atomic insights into amyloid formation, but are also helpful for designing new or modified capping peptides and inhibitors to prevent fibril formation of the VQIVYK peptide from tau protein.     

A molecular dynamics simulation study of segment B1 of protein G

The immunoglobulin binding protein, segment B1 of protein G, has been studied experimentally as a paradigm for protein folding. This protein consists of 56 residues, includes both β sheet and α helix and contains neither disulfide bonds nor proline residues. We report an all-atom molecular dynamics study of the native manifold of the protein in explicit solvent. A 2-ns simulation starting from the nuclear magnetic resonance (NMR) structure and a 1-ns control simulation starting from the x-ray structure were performed. The difference between average structures calculated over the equilibrium portion of trajectories is smaller than the difference between their starting conformations. These simulation averages are structurally similar to the x-ray structure and differ in systematic ways from the NMR-determined structure. Partitioning of the fluctuations into fast (<20 ps) and slow (<20 ps) components indicates that the β sheet displays greater long-time mobility than does the α helix. Clore and Gronenborn [J. Mol. Biol. 223:853–856, 1992] detected two long-residence water molecules by NMR in a solution structure of segment B1 of protein G. Both molecules were found in the fully exposed regions and were proposed to be stabilized by bifurcated hydrogen bonds to the protein backbone. One of these long-residence water molecules, found near an exposed loop region, is identified in both of our simulations, and is seen to be involved in the formation of a stable water-mediated hydrogen bond bridge. The second water molecule, located near the middle of the α helix, is not seen with an exceptional residence time in either as a result of the conformation being closer to the x-ray structure in this region of the protein. Proteins 29:193–202, 1997. © 1997 Wiley-Liss, Inc.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Folding propensity of anoplin: A molecular dynamics study of the native peptide and four mutated isoforms

Anoplin, a cationic decapeptide amide GLLKRIKTLL‐NH₂ derived from venom sac of the solitary wasp Anoplius samariensis has been investigated through Molecular Dynamics. The wild‐type (WT) and four isoforms were simulated both in water and in the membrane‐mimicking solvent trifluoroethanol (TFE). In water all the investigated species, found to be in rapid equilibrium between different conformational states, can be considered as unfolded. On the other hand, in TFE all the systems enhance their rigidity and, in general, show α‐helix as the main folded conformation. Interestingly, a semi‐quantitative thermodynamic analysis has suggested that the folding driving force is not always the same being in some cases (e.g., the WT Anoplin) of entropic nature and in other cases of energetic nature. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 692–701, 2015.     

Characterization of cell-surface prion protein relative to its recombinant analogue: insights from molecular dynamics simulations of diglycosylated, membrane-bound human prion protein

The prion protein (PrP) is responsible for several fatal neurodegenerative diseases via conversion from its normal to disease-related isoform. The recombinant form of the protein is typically studied to investigate the conversion process. This constructs lacks the co- and post-translational modifications present in vivo , there the protein has two N-linked glycans and is bound to the outer leaflet of the plasma membrane via a glycosylphosphatidylinositol (GPI) anchor. The inherent flexibility and heterogeneity of the glycans, the plasticity of the GPI anchor, and the localization of the protein in a membrane make experimental structural characterization of biological constructs of cellular prion protein (PrP^C) challenging. Yet this characterization is central in determining not only the suitability of recombinant (rec)-PrP^C as a model for biological forms of the protein but also the potential role of co- and post-translational modifications on the disease process. Here, we present molecular dynamics simulations of three human prion protein constructs: (i) a protein-only construct modeling the recombinant form, (ii) a diglycosylated and soluble construct, and (iii) a diglycosylated and GPI-anchored construct bound to a lipid bilayer. We found that glycosylation and membrane anchoring do not significantly alter the structure or dynamics of PrP^C, but they do appreciably modify the accessibility of the polypeptide surface PrP^C. In addition, the simulations of membrane-bound PrP^C revealed likely recognition domains for the disease-initiating PrP^C:PrP^Sc (infectious and/or misfolded form of the prion protein) binding event and a potential mechanism for the observed inefficiency of conversion associated with differentially glycosylated PrP species.     

Tau-tubulin kinase 1 (TTBK1), a neuron-specific tau kinase candidate, is involved in tau phosphorylation and aggregation

Neurofibrillary tangles, which are major pathological hallmarks of Alzheimer's disease (AD), are composed of paired helical filaments (PHFs) containing hyperphosphorylated tau. Specific kinases regulate tau phosphorylation and are closely linked to the pathogenesis of AD. We have characterized a human tau-tubulin kinase 1 (TTBK1) gene located on chromosome 6p21.1. TTBK1 is a serine/threonine/tyrosine kinase that is conserved among species and belongs to the casein kinase 1 superfamily. It is specifically expressed in the brain, especially in the cytoplasm of cortical and hippocampal neurons. TTBK1 phosphorylates tau proteins in both a Mg2+- and a Mn2+-dependent manner. Phosphopeptide mapping and immunoblotting analysis confirmed a direct tau phosphorylation by TTBK1 at Ser198, Ser199, Ser202 and Ser422, which are also phosphorylated in PHFs. TTBK1 also induces tau aggregation in human neuronal cells in a dose-dependent manner. We conclude that TTBK1 is a neuron-specific dual kinase involved in tau phosphorylation at AD-related sites and is also associated with tau aggregation.     

Computational diagnosis of protein conformational diseases: short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia

The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins.     

Solvent sensitivity of protein aggregation in Cu,Zn superoxide dismutase: a molecular dynamics simulation study

Misfolding and aggregation of Cu, Zn Superoxide Dismutase (SOD1) is often found in amyotrophic lateral sclerosis (ALS) patients. The central apo SOD1 barrel was involved in protein maturation and pathological aggregation in ALS. In this work, we employed atomistic molecular dynamics (MD) simulations to study the conformational dynamics of SOD1^barrel monomer in different concentrations of trifluoroethanol (TFE). We find concentration dependence unusual structural and dynamical features, characterized by the local unfolding of SOD1^barrel. This partially unfolded structure is characterized by the exposure of hydrophobic core, is highly dynamic in nature, and is the precursor of aggregation seen in SOD1^barrel. Our computational studies supports the hypothesis of the formation of aggregation ‘building blocks’ by means of local unfolding of apo monomer as the mechanism of SOD1 fibrillar aggregation. The non-monotonic TFE concentration dependence of protein conformational changes was explored through simulation studies. Our results suggest that altered protein conformation and dynamics within its structure may underlie the aggregation of SOD1 in ALS.     

Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Molecular dynamics simulations reveal structural differences among wild-type NPC1 protein and its mutant forms

Abstract

NPC1 is a 25-exon gene located on the long arm of chromosome 18q11.2 and encodes NPC1, a transmembrane protein comprising 1278 amino acid residues. Mutations in the NPC1 gene can cause Niemann-Pick disease type C (NP-C), a rare autosomal-recessive neurovisceral disease. We assessed mutant protein folding using computer-based molecular dynamics (MD) simulations and molecular docking of the three most common NPC1 mutations, all of which result in changes in a cysteine-rich luminal loop region of the protein: a) I1061T is the most commonly detected variant in patients with NP-C worldwide; b) P1007A is the second most common variant, frequently detected in Portuguese, British and German patients; c) G992W occurs most often in patients of Acadian descent. Analyses of molecular structural information and related cellular physiological processes revealed that mutant NPC1 proteins exhibited altered function despite being far from the N-terminal domain cholesterol binding. MD simulations revealed that mutant I1061T protein shows remarkable instability in comparison the WT and also de other mutants, and interestingly this mutant has been identified as the most common variant. In the case of the mutant P1007A, it is presumed that this substitution promotes larger structural changes than proline due to their greater hydrophobic properties.

Structural changes related to the G992W mutation may affect the physicochemical space of G992W variant protein because tryptophan induces hydrophobic interactions. Cholesterol docking studies focused on binding recognition showed differences in the binding positions of variants versus the wild-type protein that go some way to explaining the molecular pathogenesis.

Communicated by Ramaswamy H. Sarma     

Kinetic studies of folding of the B-domain of staphylococcal protein A with molecular dynamics and a united-residue (UNRES) model of polypeptide chains

Langevin dynamics is used with our physics-based united-residue (UNRES) force field to study the folding pathways of the B-domain of staphylococcal protein A (1BDD (alpha; 46 residues)). With 400 trajectories of protein A started from the extended state (to gather meaningful statistics), and simulated for more than 35 ns each, 380 of them folded to the native structure. The simulations were carried out at the optimal folding temperature of protein A with this force field. To the best of our knowledge, this is the first simulation study of protein-folding kinetics with a physics-based force field in which reliable statistics can be gathered. In all the simulations, the C-terminal alpha-helix forms first. The ensemble of the native basin has an average RMSD value of 4 A from the native structure. There is a stable intermediate along the folding pathway, in which the N-terminal alpha-helix is unfolded; this intermediate appears on the way to the native structure in less than one-fourth of the folding pathways, while the remaining ones proceed directly to the native state. Non-native structures persist until the end of the simulations, but the native-like structures dominate. To express the kinetics of protein A folding quantitatively, two observables were used: (i) the average alpha-helix content (averaged over all trajectories within a given time window); and (ii) the fraction of conformations (averaged over all trajectories within a given time window) with Calpha RMSD values from the native structure less than 5 A (fraction of completely folded structures). The alpha-helix content grows quickly with time, and its variation fits well to a single-exponential term, suggesting fast two-state kinetics. On the other hand, the fraction of folded structures changes more slowly with time and fits to a sum of two exponentials, in agreement with the appearance of the intermediate, found when analyzing the folding pathways. This observation demonstrates that different qualitative and quantitative conclusions about folding kinetics can be drawn depending on which observable is monitored.     

Phosphorylation of tau at Ser214 mediates its interaction with 14-3-3 protein: implications for the mechanism of tau aggregation

The microtubule associated protein tau is a major component of neurofibrillary tangles in Alzheimer disease brain, however the neuropathological processes behind the formation of neurofibrillary tangles are still unclear. Previously, 14-3-3 proteins were reported to bind with tau. 14-3-3 Proteins usually bind their targets through specific serine/threonine –phosphorylated motifs. Therefore, the interaction of tau with 14-3-3 mediated by phosphorylation was investigated. In this study, we show that the phosphorylation of tau by either protein kinase A (PKA) or protein kinase B (PKB) enhances the binding of tau with 14-3-3 in vitro . The affinity between tau and 14-3-3 is increased 12- to 14-fold by phosphorylation as determined by real time surface plasmon resonance studies. Mutational analyses revealed that Ser214 is critical for the phosphorylation-mediated interaction of tau with 14-3-3. Finally, in vitro aggregation assays demonstrated that phosphorylation by PKA/PKB inhibits the formation of aggregates/filaments of tau induced by 14-3-3. As the phosphorylation at Ser214 is up-regulated in fetal brain, tau's interaction with 14-3-3 may have a significant role in the organization of the microtubule cytoskeleton in development. Also as the phosphorylation at Ser214 is up-regulated in Alzheimer's disease brain, tau's interaction with 14-3-3 might be involved in the pathology of this disease.     

Partial unfolding and refolding for structure refinement: A unified approach of geometric simulations and molecular dynamics

The most successful protein structure prediction methods to date have been template‐based modeling (TBM) or homology modeling, which predicts protein structure based on experimental structures. These high accuracy predictions sometimes retain structural errors due to incorrect templates or a lack of accurate templates in the case of low sequence similarity, making these structures inadequate in drug‐design studies or molecular dynamics simulations. We have developed a new physics based approach to the protein refinement problem by mimicking the mechanism of chaperons that rehabilitate misfolded proteins. The template structure is unfolded by selectively (targeted) pulling on different portions of the protein using the geometric based technique FRODA, and then refolded using hierarchically restrained replica exchange molecular dynamics simulations (hr‐REMD). FRODA unfolding is used to create a diverse set of topologies for surveying near native‐like structures from a template and to provide a set of persistent contacts to be employed during re‐folding. We have tested our approach on 13 previous CASP targets and observed that this method of folding an ensemble of partially unfolded structures, through the hierarchical addition of contact restraints (that is, first local and then nonlocal interactions), leads to a refolding of the structure along with refinement in most cases (12/13). Although this approach yields refined models through advancement in sampling, the task of blind selection of the best refined models still needs to be solved. Overall, the method can be useful for improved sampling for low resolution models where certain of the portions of the structure are incorrectly modeled. Proteins 2015; 83:2279–2292. © 2015 Wiley Periodicals, Inc.     

Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation study

Substrate phosphorylation by cAMP-dependent-protein kinase A (protein kinase A, PKA) has been studied extensively. Phosphoryl transfer was found to be fast, whereas ADP release was found to be the slow, rate-limiting step. There is also evidence that ADP release may be preceded by a partially rate-limiting conformational change. However, the atomic details of the conformational change and the mode of ADP release are difficult to obtain experimentally. In this work, we studied ADP release from PKA by carrying out molecular dynamics simulations with different pulling forces applied to the ligand. The detailed ADP release pathway and the associated conformational changes were analyzed. The ADP release process was found to involve a swinging motion with the phosphate of ADP anchored to the Gly-rich loop, so that the more buried adenine base and ribose ring came out before the phosphate. In contrast to the common belief that a hinge-bending motion was responsible for the opening of the ligand-binding cleft, our simulations showed that the small lobe exhibited a large amplitude "rocking" motion when the ligand came out. The largest conformational change of the protein was observed at about the first quarter time point along the release pathway. Two prominent intermediate states were observed in the release process.     

Structural stability of disulfide mutants of basic pancreatic trypsin inhibitor: A molecular dynamics study

The structure and folding of basic pancreatic trypsin inhibitor (BPTI) has been studied extensively by experimental means. We report a computer simulation study of the structural stability of various disulfide mutants of BPTI, involving eight 250-psec molecular dynamics simulations of the proteins in water, with and without a phosphate counterion. The presence of the latter alters the relative stability of the single disulfide species [5–55] and [30–51]. This conclusion can explain results of mutational studies and the conservation of residues in homologues of BPTI, and suggests a possible role of ions in stabilizing one intermediate over another in unfolding or folding processes. © 1996 Wiley-Liss, Inc.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.