Speed limit of protein folding evidenced in secondary structure dynamics

Novel Covalent Modifier-Induced Local Conformational Changes within the Intrinsically Disordered Region of the Androgen Receptor

Intrinsically disordered regions (IDRs) of transcription factors play an important biological role in liquid condensate formation and gene regulation. It is thus desirable to investigate the druggability of IDRs and how small-molecule binders can alter their conformational stability. For the androgen receptor (AR), certain covalent ligands induce important changes, such as the neutralization of the condensate. To understand the specificity of ligand-IDR interaction and potential implications for the mechanism of neutralizing liquid-liquid phase separation (LLPS), we modeled and performed computer simulations of ligand-bound peptide segments obtained from the human AR. We analyzed how different covalent ligands affect local secondary structure, protein contact map, and protein-ligand contacts for these protein systems. We find that effective neutralizers make specific interactions (such as those between cyanopyrazole and tryptophan) that alter the helical propensity of the peptide segments. These findings on the mechanism of action can be useful for designing molecules that influence IDR structure and condensate of the AR in the future.

Transactivation domain of Adenovirus Early Region 1A (E1A): Investigating folding dynamics and aggregation

Transactivation domain of Adenovirus Early region 1A (E1A) oncoprotein is an intrinsically disordered molecular hub protein. It is involved in binding to different domains of human cell transcriptional co-activators such as retinoblastoma (pRb), CREB-binding protein (CBP), and its paralogue p300. The conserved region 1 (TAD) of E1A is known to undergo structural transitions and folds upon interaction with transcriptional adaptor zinc finger 2 (TAZ2). Previous reports on Taz2-E1A studies have suggested the formation of helical conformations of E1A-TAD. However, the folding behavior of the TAD region in isolation has not been studied in detail. Here, we have elucidated the folding behavior of E1A peptide at varied temperatures and solution conditions. Further, we have studied the effects of macromolecular crowding on E1A-TAD peptide. Additionally, we have also predicted the molecular recognition features of E1A using MoRF predictors. The predicted MoRFs are consistent with its structural transitions observed during TAZ2 interactions for transcriptional regulation in literature. Also, as a general rule of MoRFs, E1A undergoes helical transitions in alcohol and osmolyte solution. Finally, we studied the aggregation behavior of E1A, where we observed that the E1A could form amyloid-like aggregates that are cytotoxic to mammalian cells.

Dissecting Multiple Pathways in the Relaxation Dynamics of Helix <==> Coil Transitions with Optimum Dimensionality Reduction

Fast kinetic experiments with dramatically improved time resolution have contributed significantly to understanding the fundamental processes in protein folding pathways involving the formation of a-helices and b-hairpin, contact formation, and overall collapse of the peptide chain. Interpretation of experimental results through application of a simple statistical mechanical model was key to this understanding. Atomistic description of all events observed in the experimental findings was challenging. Recent advancements in theory, more sophisticated algorithms, and a true long-term trajectory made way for an atomically detailed description of kinetics, examining folding pathways, validating experimental results, and reporting new findings for a wide range of molecular processes in biophysical chemistry. This review describes how optimum dimensionality reduction theory can construct a simplified coarse-grained model with low dimensionality involving a kinetic matrix that captures novel insights into folding pathways. A set of metastable states derived from molecular dynamics analysis generate an optimally reduced dimensionality rate matrix following transition pathway analysis. Analysis of the actual long-term simulation trajectory extracts a relaxation time directly comparable to the experimental results and confirms the validity of the combined approach. The application of the theory is discussed and illustrated using several examples of helix <==> coil transition pathways. This paper focuses primarily on a combined approach of time-resolved experiments and long-term molecular dynamics simulation from our ongoing work.

Protons Are Fast and Smart; Proteins Are Slow and Dumb: On the Relationship of Electrospray Ionization Charge States and Conformations

We present simple considerations of how differences in time scales of motions of protons, the lightest and fastest chemical moiety, and the much longer time scales associated with the dynamics of proteins, among the heaviest and slowest analytes, may allow many protein conformations from solution to be kinetically trapped during the process of electrospraying protein solutions into the gas phase. In solution, the quantum nature of protons leads them to change locations by tunneling, an instantaneous process; moreover, the Grotthuss mechanism suggests that these small particles can respond nearly instantaneously to the dynamic motions of proteins that occur on much longer time scales. A conformational change is accompanied by favorable or unfavorable variations in the free energy of the system, providing the impetus for solvent ↔ protein proton exchange. Thus, as thermal distributions of protein conformations interconvert, protonation states rapidly respond, as specific acidic and basic sites are exposed or protected. In the vacuum of the mass spectrometer, protons become immobilized in locations that are specific to the protein conformations from which they were incorporated. In this way, conformational states from solution are preserved upon electrospraying them into the gas phase. These ideas are consistent with the exquisite sensitivity of electrospray mass spectra to small changes of the local environment that alter protein structure in solution. We might remember this approximation for the protonation of proteins in solution with the colloquial expression-protons are fast and smart; proteins are slow and dumb.

Length Dependent Folding Kinetics of Alanine-Based Helical Peptides from Optimal Dimensionality Reduction

We present a computer simulation study of helix folding in alanine homopeptides (ALA)n of length n = 5, 8, 15, and 21 residues. Based on multi-microsecond molecular dynamics simulations at room temperature, we found helix populations and relaxation times increasing from about 6% and ~2 ns for ALA5 to about 60% and ~500 ns for ALA21, and folding free energies decreasing linearly with the increasing number of residues. The helix folding was analyzed with the Optimal Dimensionality Reduction method, yielding coarse-grained kinetic models that provided a detailed representation of the folding process. The shorter peptides, ALA5 and ALA8, tended to convert directly from coil to helix, while ALA15 and ALA21 traveled through several intermediates. Coarse-grained aggregate states representing the helix, coil, and intermediates were heterogeneous, encompassing multiple peptide conformations. The folding involved multiple pathways and interesting intermediate states were present on the folding paths, with partially formed helices, turns, and compact coils. Statistically, helix initiation was favored at both termini, and the helix was most stable in the central region. Importantly, we found the presence of underlying universal local dynamics in helical peptides with correlated transitions for neighboring hydrogen bonds. Overall, the structural and dynamical parameters extracted from the trajectories are in good agreement with experimental observables, providing microscopic insights into the complex helix folding kinetics.

The Sensitivity of the Pair-Angle Distribution Function to Protein Structure

The continued development of X-ray free-electron lasers and serial crystallography techniques has opened up new experimental frontiers. Nanoscale dynamical processes such as crystal growth can now be probed at unprecedented time and spatial resolutions. Pair-angle distribution function (PADF) analysis is a correlation-based technique that has the potential to extend the limits of current serial crystallography experiments, by relaxing the requirements for crystal order, size and number density per exposure. However, unlike traditional crystallographic methods, the PADF technique does not recover the electron density directly. Instead it encodes substantial information about local three-dimensional structure in the form of three- and four-body correlations. It is not yet known how protein structure maps into the many-body PADF correlations. In this paper, we explore the relationship between the PADF and protein conformation. We calculate correlations in reciprocal and real space for model systems exhibiting increasing degrees of order and secondary structural complexity, from disordered polypeptides, single alpha helices, helix bundles and finally a folded 100 kilodalton protein. These models systems inform us about the distinctive angular correlations generated by bonding, polypeptide chains, secondary structure and tertiary structure. They further indicate the potential to use angular correlations as a sensitive measure of conformation change that is complementary to existing structural analysis techniques.

The Membrane-Integrated Steric Chaperone Lif Facilitates Active Site Opening of Pseudomonas aeruginosa Lipase A

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipase-specific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA. © 2019 Wiley Periodicals, Inc.

Kinetic pathway analysis of an α-helix in two protonation states: Direct observation and optimal dimensionality reduction

Thermodynamically stable conformers of secondary structural elements make a stable tertiary/quaternary structure that performs its proper biological function efficiently. Formation mechanisms of secondary and tertiary/quaternary structural elements from the primary structure are driven by the kinetic properties of the respective systems. Here we have carried out thermodynamic and kinetic characterization of an alpha helical heteropeptide in two protonation states, created with the addition and removal of a proton involving a single histidine residue in the primary structure. Applying far-UV circular dichroism spectroscopy, the alpha helix is observed to be significantly more stable in the deprotonated state. Nanosecond laser temperature jump spectroscopy monitoring time-resolved tryptophan fluorescence on the protonated conformer is carried out to measure the kinetics of this system. The measured relaxation rates at a final temperature between 296K and 314 K generated a faster component of 20 ns-11 ns and a slower component of 314 ns-198 ns. Atomically detailed characterization of the helix-coil kinetic pathways is performed based on all-atom molecular dynamics trajectories of the two conformers. Application of clustering and kinetic coarse-graining with optimum dimensionality reduction produced description of the trajectories in terms of kinetic models with two to five states. These models include aggregate states corresponding to helix, coil, and intermediates. The "coil" state involves the largest number of conformations, consistent with the expected high entropy of this structural ensemble. The "helix" aggregate states are found to be mixed with the full helix and partially folded forms. The experimentally observed higher helix stability in the deprotonated form of the alpha helical heteropeptide is reflected in the nature of the "helix" aggregate state arising from the kinetic model. In the protonated form, the "coil" state exhibits the lowest free energy and longest lifetime, while in the deprotonated form, it is the "helix" that is found to be most stable. Overall, the coarse grained models suggest that the protonation of a single histidine residue in the primary structure induces significant changes in the free energy landscape and kinetic network of the studied helix-forming heteropeptide.

Helix-Coil Transition Courses Through Multiple Pathways and Intermediates: Fast Kinetic Measurements and Dimensionality Reduction

Nanosecond laser temperature jumps with tryptophan fluorescence detection and molecular dynamics simulation with kinetic dimensionality reduction were used to study the helix-coil transition in a 21-residue α-helical heteropeptide. Analysis of the temperature- dependent relaxation dynamics of this heteropeptide identified a distinct faster component of 20-35 ns, besides a slower component of 300-400 ns at temperatures between 296 and 280 K. To understand the mechanism of progression from a non-structured coil state to a structured helical state, we carried out a 12 μs molecular dynamics simulation of this peptide system. Clustering and optimal dimensionality reduction were applied to the molecular dynamics trajectory to generate low-dimensional coarse-grained models of the underlying kinetic network in terms of 2-5 metastable states. In accord with the generally accepted understanding of the multiple conformations and high entropy of the unfolded ensemble of states, the "coil" metastable set contains the largest number of structures. Interestingly, the helix metastable state was also found to be structurally heterogeneous, consisting of the completely helical form and several partly folded conformers that interconvert at a time scale faster that global folding. The intermediate states contain the fewest structures, have lowest populations, and have the shortest lifetimes. As the number of considered metastable states increases, more intermediates and more folding paths appear in the coarse-grained models. One of these intermediates corresponds to the transition state for folding, which involves an "off-center" helical region over residues 11-16. The kinetic network model is consistent with a statistical picture of folding following a simple reaction coordinate counting the helical population of individual residues. On the basis of simulations, we propose that the fast relaxation time should be assigned to cooperative folding/unfolding of segments of 1-4 neighboring residues.

Directly monitor protein rearrangement on a nanosecond-to-millisecond time-scale

In order to directly observe the refolding kinetics from a partially misfolded state to a native state in the bottom of the protein-folding funnel, we used a "caging" strategy to trap the β-sheet structure of ubiquitin in a misfolded conformation. We used molecular dynamics simulation to generate the cage-induced, misfolded structure and compared the structure of the misfolded ubiquitin with native ubiquitin. Using laser flash irradiation, the cage can be cleaved from the misfolded structure within one nanosecond, and we monitored the refolding kinetics of ubiquitin from this misfolded state to the native state by photoacoustic calorimetry and photothermal beam deflection techniques on nanosecond to millisecond timescales. Our results showed two refolding events in this refolding process. The fast event is shorter than 20 ns and corresponds to the instant collapse of ubiquitin upon cage release initiated by laser irradiation. The slow event is ~60 μs, derived from a structural rearrangement in β-sheet refolding. The event lasts 10 times longer than the timescale of β-hairpin formation for short peptides as monitored by temperature jump, suggesting that rearrangement of a β-sheet structure from a misfolded state to its native state requires more time than ab initio folding of a β-sheet.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.

Observation of α-Helical Hydrogen-Bond Cooperativity in an Intact Protein

The presence and extent of hydrogen-bonding (H-bonding) cooperativity in proteins remains a fundamental question, which in the past has been studied extensively, mostly by infrared and fluorescence measurements on model peptides. We demonstrate that such cooperativity can be studied in an intact protein by hydrogen/deuterium exchange NMR spectroscopy. The method is based on the fact that substitution of NH by ND in a backbone amide group slightly weakens the N-H···O═C hydrogen bond. Our results show that such substitution at position i in an α-helix impacts the (1)H and (15)N chemical shifts of the amide sites of residues i - 3 to i + 3. Quantum mechanical calculations indicate that the upfield shifts of (1)H and (15)N resonances at site i, observed upon H/D exchanges at sites i - 3, i + 1, i + 2, and i + 3, correspond to a decrease of the ith backbone amide electric dipole moment, which weakens its H-bonding and long-range electrostatic interactions with other backbone amides in the α-helix. These results provide new quantitative insights into the cooperativity of H-bonding in protein α-helices.

Conformational Changes of Trialanine in Water Induced by Vibrational Relaxation of the Amide I Mode

Most of the protein-based diseases are caused by anomalies in the functionality and stability of these molecules. Experimental and theoretical studies of the conformational dynamics of proteins are becoming in this respect essential to understand the origin of these anomalies. However, a description of the conformational dynamics of proteins based on mechano-energetic principles still remains elusive because of the intrinsic high flexibility of the peptide chains, the participation of weak noncovalent interactions, and the role of the ubiquitous water solvent. In this work, the conformational dynamics of trialanine dissolved in water (D2O) is investigated through Molecular Dynamics (MD) simulations combined with instantaneous normal modes (INMs) analysis both at equilibrium and after the vibrational excitation of the C-terminal amide I mode. The conformational equilibrium between α and pPII conformers is found to be altered by the intramolecular relaxation of the amide I mode as a consequence of the different relaxation pathways of each conformer which modify the amount of vibrational energy stored in the torsional motions of the tripeptide, so the α → pPII and pPII → α conversion rates are increased differently. The selectivity of the process comes from the shifts of the vibrational frequencies with the conformational changes that modify the resonance conditions driving the intramolecular energy flows.

Role of side-chain interactions on the formation of α-helices in model peptides

The role played by side-chain interactions on the formation of α-helices is studied using extensive all-atom molecular dynamics simulations of polyalanine-like peptides in explicit TIP4P water. The peptide is described by the OPLS-AA force field except for the Lennard-Jones interaction between Cβ-Cβ atoms, which is modified systematically. We identify values of the Lennard-Jones parameter that promote α-helix formation. To rationalize these results, potentials of mean force (PMF) between methane-like molecules that mimic side chains in our polyalanine-like peptides are computed. These PMF exhibit a complex distance dependence where global and local minima are separated by an energy barrier. We show that α-helix propensity correlates with values of these PMF at distances corresponding to Cβ-Cβ of i-i+3 and other nearest neighbors in the α-helix. In particular, the set of Lennard-Jones parameters that promote α-helices is characterized by PMF that exhibit a global minimum at distances corresponding to i-i+3 neighbors in α-helices. Implications of these results are discussed.

Molecular modeling of peptides

This article presents a review of the field of molecular modeling of peptides. The main focus is on atomistic modeling with molecular mechanics potentials. The description of peptide conformations and solvation through potentials is discussed. Several important computer simulation methods are briefly introduced, including molecular dynamics, accelerated sampling approaches such as replica-exchange and metadynamics, free energy simulations and kinetic network models like Milestoning. Examples of recent applications for predictions of structure, kinetics, and interactions of peptides with complex environments are described. The reliability of current simulation methods is analyzed by comparison of computational predictions obtained using different models with each other and with experimental data. A brief discussion of coarse-grained modeling and future directions is also presented.

Dominance of misfolded intermediates in the dynamics of α-helix folding

Helices are the "hydrogen atoms" of biomolecular complexity; the DNA/RNA double hairpin and protein α-helix ubiquitously form the building blocks of life's constituents at the nanometer scale. Nevertheless, the formation processes of these structures, especially the dynamical pathways and rates, remain challenging to predict and control. Here, we present a general analytical method for constructing dynamical free-energy landscapes of helices. Such landscapes contain information about the thermodynamic stabilities of the possible macromolecular conformations, as well as about the dynamic connectivity, thus enabling the visualization and computation of folding pathways and timescales. We elucidate the methodology using the folding of polyalanine, and demonstrate that its α-helix folding kinetics is dominated by misfolded intermediates. At the physiological temperature of T = 298 K and midfolding time t = 250 ns, the fraction of structures in the native-state (α-helical) basin equals 22%, which is in good agreement with time-resolved experiments and massively distributed, ensemble-convergent molecular-dynamics simulations. We discuss the prominent role of β-strand-like intermediates in flight toward the native fold, and in relation to the primary conformational change precipitating aggregation in some neurodegenerative diseases.

Detailed microscopic unfolding pathways of an α-helix and a β-hairpin: direct observation and molecular dynamics

We present a combined experimental and computational study of unfolding pathways of a model 21-residue α-helical heteropeptide (W1H5-21) and a 16-residue β-hairpin (GB41-56). Experimentally, we measured fluorescence energy transfer efficiency as a function of temperature, employing natural tryptophans as donors and dansylated lysines as acceptors. Secondary structural analysis was performed with circular dichroism and Fourier transform infrared spectroscopy. Our studies present markedly different unfolding pathways of the two elementary secondary structural elements. During thermal denaturation, the helical peptide exhibits an initial decrease in length, followed by an increase, while the hairpin undergoes a systematic increase in length. In the complementary computational part of the project, we performed microsecond length replica-exchange molecular dynamics simulations of the peptides in explicit solvent, yielding a detailed microscopic picture of the unfolding processes. For the α-helical peptide, we found a large heterogeneous population of intermediates that are primarily frayed single helices or helix-turn-helix motifs. Unfolding starts at the termini and proceeds through a stable helical region in the interior of the peptide but shifted off-center toward the C-terminus. The simulations explain the experimentally observed non-monotonic variation of helix length with temperature as due primarily to the presence of frayed-end single-helix intermediate structures. For the β-hairpin peptide, our simulations indicate that folding is initiated at the turn, followed by formation of the hairpin in zipper-like fashion, with Cα···Cα contacts propagating from the turn to termini and hairpin hydrogen bonds forming in parallel with these contacts. In the early stages of hairpin formation, the hydrophobic side-chain contacts are only partly populated. Intermediate structures with low numbers of β-hairpin hydrogen bonds have very low populations. This is in accord with the "broken zipper" model of Scheraga. The monotonic increase in length with temperature may be explained by the zipper-like breaking of the hairpin hydrogen bonds and backbone contacts.

Effects of side chains in helix nucleation differ from helix propagation

Helix-coil transition theory connects observable properties of the α-helix to an ensemble of microstates and provides a foundation for analyzing secondary structure formation in proteins. Classical models account for cooperative helix formation in terms of an energetically demanding nucleation event (described by the σ constant) followed by a more facile propagation reaction, with corresponding s constants that are sequence dependent. Extensive studies of folding and unfolding in model peptides have led to the determination of the propagation constants for amino acids. However, the role of individual side chains in helix nucleation has not been separately accessible, so the σ constant is treated as independent of sequence. We describe here a synthetic model that allows the assessment of the role of individual amino acids in helix nucleation. Studies with this model lead to the surprising conclusion that widely accepted scales of helical propensity are not predictive of helix nucleation. Residues known to be helix stabilizers or breakers in propagation have only a tenuous relationship to residues that favor or disfavor helix nucleation.

Side-chain-side-chain interactions and stability of the helical state

Understanding the driving forces that lead to the stability of the secondary motifs found in proteins, namely α-helix and β-sheet, is a major goal in structural biology. The thermodynamic stability of these repetitive units is a result of a delicate balance between many factors, which in addition to the peptide chain involves also the solvent. Despite the fact that the backbones of all amino acids are the same (except of that of proline), there are large differences in the propensity of the different amino acids to promote the helical structure. In this paper, we investigate by explicit-solvent molecular dynamics simulations the role of the side chains (modeled as coarse-grained single sites) in stabilizing α helices in an aqueous solution. Our model systems include four (six-mer-nine-mer) peptide lengths in which the magnitude of the effective attraction between the side chains is systematically increased. We find that these interactions between the side chains can induce (for the nine-mer almost completely) a transition from a coil to a helical state. This transition is found to be characterized by three states in which the intermediate state is a partially folded α-helical conformation. In the absence of any interactions between the side chains the free energy change for helix formation has a small positive value indicating that favorable contributions from the side chains are necessary to stabilize the helical conformation. Thus, the helix-coil transition is controlled by the effective potentials between the side-chain residues and the magnitude of the required attraction per residue, which is on the order of the thermal energy, reduces with the length of the peptide. Surprisingly, the plots of the population of the helical state (or the change in the free energy for helix formation) as a function of the total effective interactions between the side chains in the helical state for all peptide lengths fall on the same curve.

Non-exponential kinetics and a complete folding pathway of an α-helical heteropeptide: direct observation and comprehensive molecular dynamics

We have performed a combined experimental and computational study of the folding of a 21-residue α-helix-forming heteropeptide (WH21). Temperature jump kinetics with improved dynamic range at several temperatures revealed non-exponential relaxation that could be well described with two time constants of 20 and 300 ns at 298 K. In the computational part, we performed multi-microsecond molecular dynamics simulations of WH21 in explicit water, using the AMBER03 and OPLS/AA potentials. The simulations were in good agreement with available experimental data on helix content and relaxation times. On the basis of 70 individual transitions, we identified the main pathways of helix unfolding. Three paths were found in both force fields, with unfolding progressing through (1) N-terminus, C-terminus, and center; (2) C-terminus, N-terminus, and center; and (3) C-terminus, center, and N-terminus. An additional fourth path starting in the central region and expanding to the termini was detected only with AMBER03. Intermediate structures sampled along the pathway included a central helix with frayed termini, an off-center helix, and a helical hairpin. The simulations suggest that the short relaxation should be assigned to partly cooperative fluctuations of several neighboring hydrogen bonds. Overall, by a combination of ultrafast kinetic measurements and detailed microscopic description through comprehensive molecular dynamics, we have obtained important new insights into the helix folding process.

A membrane-proximal, C-terminal α-helix is required for plasma membrane localization and function of the G Protein-coupled receptor (GPCR) TGR5

The C terminus of G protein-coupled receptors (GPCRs) is important for G protein-coupling and activation; in addition, sorting motifs have been identified in the C termini of several GPCRs that facilitate correct trafficking from the endoplasmic reticulum to the plasma membrane. The C terminus of the GPCR TGR5 lacks any known sorting motif such that other factors must determine its trafficking. Here, we investigate deletion and substitution variants of the membrane-proximal C terminus of TGR5 with respect to plasma membrane localization and function using immunofluorescence staining, flow cytometry, and luciferase assays. Peptides of the membrane-proximal C-terminal variants are subjected to molecular dynamics simulations and analyzed with respect to their secondary structure. Our results reveal that TGR5 plasma membrane localization and responsiveness to extracellular ligands is fostered by a long (≥ 9 residues) α-helical stretch at the C terminus, whereas the presence of β-strands or only a short α-helical stretch leads to retention in the endoplasmic reticulum and a loss of function. As a proof-of-principle, chimeras of TGR5 containing the membrane-proximal amino acids of the β2 adrenergic receptor (β2AR), the sphingosine 1-phosphate receptor-1 (S1P1), or the κ-type opioid receptor (κOR) were generated. These TGR5β2AR, TGR5S1P1, or TGR5κOR chimeras were correctly sorted to the plasma membrane. As the exchanged amino acids of the β2AR, the S1P1, or the κOR form α-helices in crystal structures but lack significant sequence identity to the respective TGR5 sequence, we conclude that the secondary structure of the TGR5 membrane-proximal C terminus is the determining factor for plasma membrane localization and responsiveness towards extracellular ligands.

Predicting the Thermodynamics and Kinetics of Helix Formation in a Cyclic Peptide Model

The peptide Ac-(cyclo-2,6)-R[KAAAD]-NH2 (cyc-RKAAAD) is a short cyclic peptide known to adopt a remarkably stable single turn α-helix in water. Due to its simplicity and the availability of thermodynamic and kinetic experimental data, cyc-RKAAAD poses as an ideal model for evaluating the aptness of current molecular dynamics (MD) simulation methodologies to accurately sample conformations that reproduce experimentally observed properties. In this work, we extensively sample the conformational space of cyc-RKAAAD using microsecond-timescale MD simulations. We characterize the peptide conformational preferences in terms of secondary structure propensities and, using Cartesian-coordinate principal component analysis (cPCA), construct its free energy landscape, thus obtaining a detailed weighted discrimination between the helical and nonhelical subensembles. The cPCA state discrimination, together with a Markov model built from it, allowed us to estimate the free energy of unfolding (-0.57 kJ/mol) and the relaxation time (∼0.435 μs) at 298.15 K, which are in excellent agreement with the experimentally reported values (-0.22 kJ/mol and 0.42 μs, Serrano, A. L.; Tucker, M. J.; Gai, F. J. Phys. Chem. B, 2011, 115, 7472-7478.). Additionally, we present simulations conducted using two enhanced sampling methods: replica-exchange molecular dynamics (REMD) and bias-exchange metadynamics (BE-MetaD). We compare the free energy landscape obtained by these two methods with the results from MD simulations and discuss the sampling and computational gains achieved. Overall, the results obtained attest to the suitability of modern simulation methods to explore the conformational behavior of peptide systems with a high level of realism.

Amyloid beta(1-42) in aqueous environments: effects of ionic strength and E22Q (Dutch) mutation

Development of extracellular plaques characteristic of Alzheimer's disease is related to aggregation of amyloid peptides. The Aβ-42 peptide is the most aggregation prone species, and some missense mutant forms increase this aggregation ability. Due to its poor solubility as monomer in aqueous solutions, Aβ-42 conformational transitions in water have been largely investigated by molecular dynamics. Here we report an all-atom molecular dynamics analysis of the Aβ-42 peptide in aqueous environment using as starting conformation a structure obtained in an isotropic, low-polarity medium, representing a plausible model for the membrane-bound species. While previous studies commonly show that Aβ-42 is largely unstructured in aqueous solution, here we report that this peptide can adopt partially folded structures. Importance of ionic strength has been also investigated, showing that at physiological ionic strength condition a loop stabilizing electrostatic interaction involving Lys28 builds up. In addition, besides stable α-helix structures, we observe the appearance of 310 helix, similar to what was reported experimentally for the Aβ-40 species. The effect of E22Q (Dutch) mutation in high ionic strength condition has been explored. We show that this mutation has a dramatic impact on the Aβ-42 structure. Instead of a partially folded, but extended, conformation obtained with the wild type, the E22Q assumes a two-helix collapsed one due to the clustering of hydrophobic residues.

Testing the diffusing boundary model for the helix-coil transition in peptides

The dynamics of peptide α-helices have been studied extensively for many years, and the kinetic mechanism of the helix-coil dynamics has been discussed controversially. Recent experimental results have suggested that equilibrium helix-coil dynamics are governed by movement of the helix/coil boundary along the peptide chain, which leads to slower unfolding kinetics in the helix center compared with the helix ends and position-independent helix formation kinetics. We tested this diffusion of boundary model in helical peptides of different lengths by triplet-triplet energy transfer measurements and compared the data with simulations based on a kinetic linear Ising model. The results show that boundary diffusion in helical peptides can be described by a classical, Einstein-type, 1D diffusion process with a diffusion coefficient of 2.7⋅10(7) (amino acids)(2)/s or 6.1⋅10(-9) cm(2)/s. In helices with a length longer than about 40 aa, helix unfolding by coil nucleation in a helical region occurs frequently in addition to boundary diffusion. Boundary diffusion is slowed down by helix-stabilizing capping motifs at the helix ends in agreement with predictions from the kinetic linear Ising model. We further tested local and nonlocal effects of amino acid replacements on helix-coil dynamics. Single amino acid replacements locally affect folding and unfolding dynamics with a ϕf-value of 0.35, which shows that interactions leading to different helix propensities for different amino acids are already partially present in the transition state for helix formation. Nonlocal effects of amino acid replacements only influence helix unfolding (ϕf = 0) in agreement with a diffusing boundary mechanism.

Automated Optimization of Potential Parameters

An algorithm and software to refine parameters of empirical energy functions according to condensed phase experimental measurements are discussed. The algorithm is based on sensitivity analysis and local minimization of the differences between experiment and simulation as a function of potential parameters. It is illustrated for a toy problem of alanine dipeptide and is applied to folding of the peptide WAAAH. The helix fraction is highly sensitive to the potential parameters while the slope of the melting curve is not. The sensitivity variations make it difficult to satisfy both observations simultaneously. We conjecture that there is no set of parameters that reproduces experimental melting curves of short peptides that are modeled with the usual functional form of a force field.

On the foldability of tryptophan-containing tetra- and pentapeptides: an exhaustive molecular dynamics study

Short peptides serve as minimal model systems to decipher the determinants of foldability due to their simplicity arising from their smaller size, their ability to echo protein-like structural characteristics, and their direct implication in force field validation. Here, we describe an effort to identify small peptides that can still form stable structures in aqueous solutions. We followed the in silico folding of a selected set of 8640 tryptophan-containing tetra- and pentapeptides through 15 210 molecular dynamics simulations amounting to a total of 272.46 μs using explicit representation of the solute and full treatment of the electrostatics. The evaluation and sorting of peptides is achieved through scoring functions, which include terms based on interatomic vector distances, atomic fluctuations, and rmsd matrices between successive frames of a trajectory. Highly scored peptides are studied further via successive simulation rounds of increasing simulation length and using different empirical force fields. Our method suggested only a handful of peptides with strong foldability prognosis. The discrepancies between the predictions of the various force fields for such short sequences are also extensively discussed. We conclude that the vast majority of such short peptides do not adopt significantly stable structures in water solutions, at least based on our computational predictions. The present work can be utilized in the rational design and engineering of bioactive peptides with desired molecular properties.

Consistent View of Protein Fluctuations from All-Atom Molecular Dynamics and Coarse-Grained Dynamics with Knowledge-Based Force-Field

It is widely recognized that atomistic Molecular Dynamics (MD), a classical simulation method, captures the essential physics of protein dynamics. That idea is supported by a theoretical study showing that various MD force-fields provide a consensus picture of protein fluctuations in aqueous solution [Rueda, M. et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 796-801]. However, atomistic MD cannot be applied to most biologically relevant processes due to its limitation to relatively short time scales. Much longer time scales can be accessed by properly designed coarse-grained models. We demonstrate that the aforementioned consensus view of protein dynamics from short (nanosecond) time scale MD simulations is fairly consistent with the dynamics of the coarse-grained protein model - the CABS model. The CABS model employs stochastic dynamics (a Monte Carlo method) and a knowledge-based force-field, which is not biased toward the native structure of a simulated protein. Since CABS-based dynamics allows for the simulation of entire folding (or multiple folding events) in a single run, integration of the CABS approach with all-atom MD promises a convenient (and computationally feasible) means for the long-time multiscale molecular modeling of protein systems with atomistic resolution.

Experiments and comprehensive simulations of the formation of a helical turn

We investigate the kinetics and thermodynamics of a helical turn formation in the peptide Ac-WAAAH-NH(2). NMR measurements indicate that this peptide has significant tendency to form a structure of a helical turn, while temperature dependent CD establishes the helix fraction at different temperatures. Molecular dynamics and milestoning simulations agree with experimental observables and suggest an atomically detailed picture for the turn formation. Using a network representation, two alternative mechanisms of folding are identified: (i) a direct co-operative mechanism from the unfolded to the folded state without intermediate formation of hydrogen bonds and (ii) an indirect mechanism with structural intermediates with two residues in a helical conformation. This picture is consistent with kinetic measurements that reveal two experimental time scales of sub-nanosecond and several nanoseconds.

Towards the physical basis of how intrinsic disorder mediates protein function

Intrinsically disordered proteins (IDPs) are an important class of functional proteins that is highly prevalent in biology and has broad association with human diseases. In contrast to structured proteins, free IDPs exist as heterogeneous and dynamical conformational ensembles under physiological conditions. Many concepts have been discussed on how such intrinsic disorder may provide crucial functional advantages, particularly in cellular signaling and regulation. Establishing the physical basis of these proposed phenomena requires not only detailed characterization of the disordered conformational ensembles, but also mechanistic understanding of the roles of various ensemble properties in IDP interaction and regulation. Here, we review the experimental and computational approaches that may be integrated to address many important challenges of establishing a "structural" basis of IDP function, and discuss some of the key emerging ideas on how the conformational ensembles of IDPs may mediate function, especially in coupled binding and folding interactions.

Residue-specific α-helix propensities from molecular simulation

Formation of α-helices is a fundamental process in protein folding and assembly. By studying helix formation in molecular simulations of a series of alanine-based peptides, we obtain the temperature-dependent α-helix propensities of all 20 naturally occurring residues with two recent additive force fields, Amber ff03w and Amber ff99SB(∗). Encouragingly, we find that the overall helix propensity of many residues is captured well by both energy functions, with Amber ff99SB(∗) being more accurate. Nonetheless, there are some residues that deviate considerably from experiment, which can be attributed to two aspects of the energy function: i), variations of the charge model used to determine the atomic partial charges, with residues whose backbone charges differ most from alanine tending to have the largest error; ii), side-chain torsion potentials, as illustrated by the effect of modifications to the torsion angles of I, L, D, N. We find that constrained refitting of residue charges for charged residues in Amber ff99SB(∗) significantly improves their helix propensity. The resulting parameters should more faithfully reproduce helix propensities in simulations of protein folding and disordered proteins.