New Insights into Protein (Un)Folding Dynamics

Reversing the typical pH stability profile of the Trp-cage

The Trp-cage, an 18-20 residue miniprotein, has emerged as a primary test system for evaluating computational fold prediction and folding rate determination efforts. As it turns out, a number of stabilizing interactions in the Trp-cage folded state have a strong pH dependence; all prior Trp-cage mutants have been destabilized under carboxylate-protonating conditions. Notable among the pH dependent stabilizing interactions within the Trp-cage are: (1) an Asp as the helix N-cap, (2) an H-bonded Asp9/Arg16 salt bridge, (3) an interaction between the chain termini which are in close spatial proximity, and (4) additional side chain interactions with Asp9. In the present study, we have prepared Trp-cage species that are significantly more stable at pH 2.5 (rather than 7) and quantitated the contribution of each interaction listed above. The Trp-cage structure remains constant with the pH change. The study has also provided measures of the stabilizing contribution of indole ring shielding from surface exposure and the destabilizing effects of an ionized Asp at the C-terminus of an α-helix.

From a Highly Disordered to a Metastable State: Uncovering Insights of α-Synuclein

α-Synuclein (αS) is a major constituent of Lewy bodies, the insoluble aggregates that are the hallmark of one of the most prevalent neurodegenerative disorders, Parkinson's disease (PD). The vast majority of experiments in vitro and in vivo provide extensive evidence that a disordered monomeric form is the predominant state of αS in water solution, and it undergoes a large-scale disorder-to-helix transition upon binding to vesicles of different types. Recently, another form, tetrameric, of αS with a stable helical structure was identified experimentally. It has been shown that a dynamic intracellular population of metastable αS tetramers and monomers coexists normally; and the tetramer plays an essential role in maintaining αS homeostasis. Therefore, it is of interest to know whether the tetramer can serve as a means of preventing or delaying the start of PD. Before answering this very important question, it is, first, necessary to find out, on an atomistic level, a correlation between tetramers and monomers; what mediates tetramer formation and what makes a tetramer stable. We address these questions here by investigating both monomeric and tetrameric forms of αS. In particular, by examining correlations between the motions of the side chains and the main chain, steric parameters along the amino-acid sequence, and one- and two-dimensional free-energy landscapes along the coarse-grained dihedral angles γ and δ and principal components, respectively, in monomeric and tetrameric αS, we were able to shed light on a fundamental relationship between monomers and tetramers, and the key residues involved in mediating formation of a tetramer. Also, the reasons for the stability of tetrameric αS and inability of monomeric αS to fold are elucidated here.

N-terminal diproline and charge group effects on the stabilization of helical conformation in alanine-based short peptides: CD studies with water and methanol as solvent

Protein folding problem remains a formidable challenge as main chain, side chain and solvent interactions remain entangled and have been difficult to resolve. Alanine-based short peptides are promising models to dissect protein folding initiation and propagation structurally as well as energetically. The effect of N-terminal diproline and charged side chains is assessed on the stabilization of helical conformation in alanine-based short peptides using circular dichroism (CD) with water and methanol as solvent. A1 (Ac-Pro-Pro-Ala-Lys-Ala-Lys-Ala-Lys-Ala-NH₂ ) is designed to assess the effect of N-terminal homochiral diproline and lysine side chains to induce helical conformation. A2 (Ac-Pro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) and A3 (Ac-dPro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) with N-terminal homochiral and heterochiral diproline, respectively, are designed to assess the effect of Glu...Lys (i, i + 4) salt bridge interactions on the stabilization of helical conformation. The CD spectra of A1, A2 and A3 in water manifest different amplitudes of the observed polyproline II (PPII) signals, which indicate different conformational distributions of the polypeptide structure. The strong effect of solvent substitution from water to methanol is observed for the peptides, and CD spectra in methanol evidence A2 and A3 as helical folds. Temperature-dependent CD spectra of A1 and A2 in water depict an isodichroic point reflecting coexistence of two conformations, PPII and β-strand conformation, which is consistent with the previous studies. The results illuminate the effect of N-terminal diproline and charged side chains in dictating the preferences for extended-β, semi-extended PPII and helical conformation in alanine-based short peptides. The results of the present study will enhance our understanding on stabilization of helical conformation in short peptides and hence aid in the design of novel peptides with helical structures. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13-26) with hydrogen bonds and π-π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16-20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM-PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (-43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.

Rationally Designed Peptides and Peptidomimetics as Inhibitors of Amyloid-β (Aβ) Aggregation: Potential Therapeutics of Alzheimer's Disease

Alzheimer's disease (AD) is a progressive neurodegenerative disease with no clinically accepted treatment to cure or halt its progression. The worldwide effort to develop peptide-based inhibitors of amyloid-β (Aβ) aggregation can be considered an unplanned combinatorial experiment. An understanding of what has been done and achieved may advance our understanding of AD pathology and the discovery of effective therapeutic agents. We review here the history of such peptide-based inhibitors, including those based on the Aβ sequence and those not derived from that sequence, containing both natural and unnatural amino acid building blocks. Peptide-based aggregation inhibitors hold significant promise for future AD therapy owing to their high selectivity, effectiveness, low toxicity, good tolerance, low accumulation in tissues, high chemical and biological diversity, possibility of rational design, and highly developed methods for analyzing their mode of action, proteolytic stability (modified peptides), and blood-brain barrier (BBB) permeability.

Sequence-, structure-, and dynamics-based comparisons of structurally homologous CheY-like proteins

We recently introduced a physically based approach to sequence comparison, the property factor method (PFM). In the present work, we apply the PFM approach to the study of a challenging set of sequences-the bacterial chemotaxis protein CheY, the N-terminal receiver domain of the nitrogen regulation protein NT-NtrC, and the sporulation response regulator Spo0F. These are all response regulators involved in signal transduction. Despite functional similarity and structural homology, they exhibit low sequence identity. PFM sequence comparison demonstrates a statistically significant qualitative difference between the sequence of CheY and those of the other two proteins that is not found using conventional alignment methods. This difference is shown to be consonant with structural characteristics, using distance matrix comparisons. We also demonstrate that residues participating strongly in native contacts during unfolding are distributed differently in CheY than in the other two proteins. The PFM result is also in accord with dynamic simulation results of several types. Molecular dynamics simulations of all three proteins were carried out at several temperatures, and it is shown that the dynamics of CheY are predicted to differ from those of NT-NtrC and Spo0F. The predicted dynamic properties of the three proteins are in good agreement with experimentally determined B factors and with fluctuations predicted by the Gaussian network model. We pinpoint the differences between the PFM and traditional sequence comparisons and discuss the informatic basis for the ability of the PFM approach to detect physical differences between these sequences that are not apparent from traditional alignment-based comparison.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.

Protein States as Symmetry Transitions in the Correlation Matrices

Over the last few years, there has been significant progress in the knowledge on protein folding. However, some aspects of protein folding still need further attention. One of these is the exact relationship between the folded and unfolded states and the differences between them. Whereas the folded state is well known, at least from a structural point of view (just think of the thousands of structures in online databases), the unfolded state is more elusive. Also, these are dynamic states of matter, and this aspect cannot be overlooked. Molecular dynamics-derived correlation matrices are an invaluable source of information on the protein dynamics. Here, bulk eigenvalue spectra of the correlation matrices obtained from the Trp-cage dynamics in the folded and unfolded states have been analyzed. The associated modes represent localized vibrations and are significantly affected by the fine details of the structure and interactions. Therefore, these bulk modes can be used as probes of the protein local dynamics in different states. The results of these analyses show that the correlation matrices describing the folded and unfolded dynamics belong to different symmetry classes. This finding provides new support to the phase-transition models of protein folding.

Computational scrutiny of the effect of N-terminal proline and residue stereochemistry in the nucleation of α-helix fold

N-Terminal l- to d-residue mutation nucleate helical fold in Ac–^DAla–^LAla₃–NHMe (Ib, m2), Ac–^DPro–^LAla₃–NHMe (IIb, m1), and Ac–^DPro–^LPro–^LAla₂–NHMe (IIIb, m2) peptides.

Probing the role of electrostatics of polypeptide main-chain in protein folding by perturbing N-terminal residue stereochemistry: DFT study with oligoalanine models

Energetics of folding (ΔH_E→F, in kcal mol⁻¹) from the extended (E) structure to the folded (F) structure for Ia and Ib critically depend on the geometrical relationship between the backbone peptide units of the polypeptide structure.

Correlation Analysis of Trp-Cage Dynamics in Folded and Unfolded States

A fundamental and still debated problem is how folded structures of proteins are related to their unfolded state. Besides the classical view, in which a large number of conformations characterize the unfolded state while the folded one is dominated by a single structure, recently a reassessment of the denatured state has been suggested. A growing amount of evidence indicates that not only the folded but also the unfolded state is at least partially organized. Here, we try to answer the question of how different protein dynamics is in folded and unfolded states by performing all-atom molecular dynamics simulations on the model protein Trp-cage. Random matrix theory inspired analysis of the correlation matrices has been carried out. The spectra of these correlation matrices show that the low rank modes of Trp-cage dynamics are outside of the limit expected for a random system both in folded and in unfolded conditions. These findings shed light on the nature of the unfolded state of the proteins, suggesting that it is much less random than previously thought.

Intrinsic Localized Modes in Proteins

Protein dynamics is essential for proteins to function. Here we predicted the existence of rare, large nonlinear excitations, termed intrinsic localized modes (ILMs), of the main chain of proteins based on all-atom molecular dynamics simulations of two fast-folder proteins and of a rigid α/β protein at 300 K and at 380 K in solution. These nonlinear excitations arise from the anharmonicity of the protein dynamics. The ILMs were detected by computing the Shannon entropy of the protein main-chain fluctuations. In the non-native state (significantly explored at 380 K), the probability of their excitation was increased by a factor between 9 and 28 for the fast-folder proteins and by a factor 2 for the rigid protein. This enhancement in the non-native state was due to glycine, as demonstrated by simulations in which glycine was mutated to alanine. These ILMs might play a functional role in the flexible regions of proteins and in proteins in a non-native state (i.e. misfolded or unfolded states).

Equilibrium transitions between side-chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side-chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side-chain dihedral angles, we performed all-atom hard-sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo-chemical constraints and repulsive-only steric interactions between non-bonded atoms. We determine the relative populations of the different χ(1) and χ(2) dihedral angle combinations as a function of the backbone dihedral angles ϕ and ψ. We also propose, and test, a mechanism for inter-conversion between the different side-chain conformations. Specifically, we discover that some of the transitions between side-chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ(2), the Leu side-chain bond angles κ(1) and κ(2) must increase, whereas to transition in χ(1), the Ile bond angles λ(1) and λ(2) must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side-chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design.