Recent advances in helix-coil theory

Effects of Charge Sequence Pattern and Lysine-to-Arginine Substitution on the Structural Stability of Bioinspired Polyampholytes

A comprehensive study focusing on the combined influence of the charge sequence pattern and the type of positively charged amino acids on the formation of secondary structures in sequence-specific polyampholytes is presented. The sequences of interest consisting exclusively of ionizable amino acids (lysine, K; arginine, R; and glutamic acid, E) are (EKEK)5, (EKKE)5, (ERER)5, (ERRE)5, and (EKER)5. The stability of the secondary structure was examined at three pH values in the presence of urea and NaCl. The results presented here underscore the combined prominent effects of the charge sequence pattern and the type of positively charged monomers on secondary structure formation. Additionally, (ERRE)5 readily aggregated across a wide range of pH. In contrast, sequences with the same charge pattern, (EKKE)5, as well as the sequences with the equivalent amino acid content, (ERER)5, exhibited no aggregate formation under equivalent pH and concentration conditions.

Chemical Protein Unfolding - A Simple Cooperative Model

Chemical unfolding with guanidineHCl or urea is a common method to study the conformational stability of proteins. The analysis of unfolding isotherms is usually performed with an empirical linear extrapolation method (LEM). A large positive free energy is assigned to the native protein, which is usually considered to be a minimum of the free energy. The method thus contradicts common expectations. Here, we present a multistate cooperative model that addresses specifically the binding of the denaturant to the protein and the cooperativity of the protein unfolding equilibrium. The model is based on a molecular statistical-mechanical partition function of the ensemble, but simple solutions for the calculation of the binding isotherm and the associated free energy are presented. The model is applied to 23 published unfolding isotherms of small and large proteins. For a given denaturant, the binding constant depends on temperature and pH but shows little protein specificity. Chemical unfolding is less cooperative than thermal unfolding. The cooperativity parameter σ is at least 2 orders of magnitude larger than that of thermal unfolding. The multistate cooperative model predicts zero free energy for the native protein, which becomes strongly negative beyond the midpoint concentration of unfolding. The free energy to unfold a cooperative unit corresponds exactly to the diffusive energy of the denaturant concentration gradient necessary for unfolding. The temperature dependence of unfolding isotherms yields the denaturant-induced unfolding entropy and, in turn, the unfolding enthalpy. The multistate cooperative model provides molecular insight and is as simple to apply as the LEM but avoids the conceptual difficulties of the latter.

Understanding the helical stability of charged peptides

Cationic helical peptides play a crucial role in applications such as anti-microbial and anticancer activity. The activity of these peptides directly correlates with their helicity. In this study, we have performed extensive all-atom molecular dynamics simulations of 25 Lysine-Leucine co-polypeptide sequences of varying charge density ( λ ) and patterns. Our findings showed that, an increase in the charge density on the peptide leads to a gradual decrease in the helicity up to a critical charge density λ c . Beyond λ c _, a complete helix to coil transition was observed. The decrease in the helicity is correlated with the increased number of water molecules in first solvation shell, solvent-exposed surface area, and a higher value of the radius of gyration of the peptide.

Molecular understanding of calorimetric protein unfolding experiments

Testing and predicting protein stability gained importance because proteins, including antibodies, became pharmacologically relevant in viral and cancer therapies. Isothermal scanning calorimetry is the principle method to study protein stability. Here, we use the excellent experimental heat capacity C_p(T) data from the literature for a critical inspection of protein unfolding as well as for the test of a new cooperative model. In the relevant literature, experimental temperature profiles of enthalpy, H_cal(T), entropy, S_cal(T), and free energy, G_cal(T) are missing. First, we therefore calculate the experimental H_cal(T), S_cal(T), and G_cal(T) from published C_p(T) thermograms. Considering only the unfolding transition proper, the heat capacity and all thermodynamic functions are zero in the region of the native protein. In particular, the free energy of the folded proteins is also zero and G_cal(T) displays a trapezoidal temperature profile when cold denaturation is included. Second, we simulate the DSC-measured thermodynamic properties with a new molecular model based on statistical-mechanical thermodynamics. The model quantifies the protein cooperativity and predicts the aggregate thermodynamic variables of the system with molecular parameters only. The new model provides a perfect simulation of all thermodynamic properties, including the observed trapezoidal G_cal(T) temperature profile. Importantly, the new cooperative model can be applied to a broad range of protein sizes, including antibodies. It predicts not only heat and cold denaturation but also provides estimates of the unfolding kinetics and allows a comparison with molecular dynamics calculations and quasielastic neutron scattering experiments.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Energetics and structure of alanine-rich α-helices via adaptive steered molecular dynamics

The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.

A Two-State Model Describes the Temperature-Dependent Conformational Equilibrium in the Alanine-Rich Domains in Elastin

Elastin is the insoluble elastomeric protein that provides extensibility and resilience to vertebrate tissues. Limited high-resolution structural data for elastin are notably complex. To access this information, this protein is considered in the simplified context of its two general domain types, that is, hydrophobic (HP) and crosslinking (CL). The question of elastin's structure-function has directed the focus of nearly all previous studies in the literature to the unique repeating sequences characteristic of this protein, found primarily in the HP domains. The CL domains were assumed to play a very limited role in biological elasticity due in part to the significant α-helical character that was (incorrectly) predicted for these regions. In this study, the conformational heterogeneity of alanines in native elastin's CL domains is examined in the context of helix-coil transition theory (HCTT) using solid-state nuclear magnetic resonance (SSNMR) spectroscopy in tandem with strategic isotopic labeling. Helix and coil populations are observed at all temperatures, but the former increases significantly at lower temperatures. Below the glass transition temperature (T_g), two major populations of alanines in the CL regions are resolved by two-dimensional SSNMR; one-dimensional methods are used for characterization in nativelike conditions. The spectra of ¹³CO-Ala in the CL regions are simulated using an HCTT-based statistical mechanical representation. Below T_g, longer segments with significant helical probabilities are consistent with the experimental data. At higher temperatures, the SSNMR lineshapes are best fit with a distribution of shorter (Ala)_n segments, most in random coil. These results are used to refine a structure-function model for elastin in the context of HCTT, redirecting attention to the CL domains and their role in elasticity.

Thermal and Chemical Unfolding of a Monoclonal IgG1 Antibody: Application of the Multistate Zimm-Bragg Theory

The thermal unfolding of a recombinant monoclonal antibody IgG1 (mAb) was measured with differential scanning calorimetry (DSC). The DSC thermograms reveal a pretransition at 72°C with an unfolding enthalpy of ΔH_cal ∼200-300 kcal/mol and a main transition at 85°C with an enthalpy of ∼900-1000 kcal/mol. In contrast to small single-domain proteins, mAb unfolding is a complex reaction that is analyzed with the multistate Zimm-Bragg theory. For the investigated mAb, unfolding is characterized by a cooperativity parameter σ ∼6 × 10^-5 and a Gibbs free energy of unfolding of g_nu ∼100 cal/mol per amino acid. The enthalpy of unfolding provides the number of amino acid residues ν participating in the unfolding reaction. On average, ν∼220 ± 50 amino acids are involved in the pretransition and ν∼850 ± 30 in the main transition, accounting for ∼90% of all amino acids. Thermal unfolding was further studied in the presence of guanidineHCl. The chemical denaturant reduces the unfolding enthalpy ΔH_cal and lowers the midpoint temperature T_m. Both parameters depend linearly on the concentration of denaturant. The guanidineHCl concentrations needed to unfold mAb at 25°C are predicted to be 2-3 M for the pretransition and 5-7 M for the main transition, varying with pH. GuanidineHCl binds to mAb with an exothermic binding enthalpy, which partially compensates the endothermic mAb unfolding enthalpy. The number of guanidineHCl molecules bound upon unfolding is deduced from the DSC thermograms. The bound guanidineHCl-to-unfolded amino acid ratio is 0.79 for the pretransition and 0.55 for the main transition. The pretransition binds more denaturant molecules and is more sensitive to unfolding than the main transition. The current study shows the strength of the Zimm-Bragg theory for the quantitative description of unfolding events of large, therapeutic proteins, such as a monoclonal antibody.

Structural features of split and unsplit βαβ-units

In this study, 1064 nonhomologous "unsplit", "one-strand split" and "two-strand split" right-handed βαβ-units having standard α-helices and loops up to seven residues in length have been analyzed. It was found that the α-helices in these kinds of βαβ-units have different distributions of the hydrophobic and hydrophilic amino acid residues along the chain. In the unsplit βαβ-units, most α-helices have hydrophobic residues in positions N4-N7-N8-N11 or N6-N7-N10, where N1 is the first N-terminal residue. In the one-strand split βαβ-units, most α-helices have hydrophobic residues in positions N4-N7-N8-N11 and those in two-strand split βαβ-units in positions N4-N5-N8-N12. On the other hand, in all kinds of βαβ-units, there are commonly occurring hydrophobic stripes of type C4-C7-C8 at the C-terminal parts of the α-helices. As a rule, the C- and N-terminal hydrophobic stripes overlap and the extent of their overlapping determine the length of α-helices.

Thermal and Chemical Unfolding of Lysozyme. Multistate Zimm-Bragg Theory Versus Two-State Model

Thermal and chemical unfolding of lysozyme in the presence of the guanidine HCl denaturant is a model system to compare the conventional two-state model of protein unfolding with the multistate Zimm-Bragg theory. The two-state model is shown to be the noncooperative limit of the Zimm-Bragg theory. In particular, the Zimm-Bragg theory provides a molecular interpretation of the empirical linear extrapolation method (LEM) of the two-state model. Differential scanning calorimetry (DSC) experiments reported in the literature are analyzed with both methods. Lysozyme unfolding is associated with a large endothermic enthalpy that decreases significantly upon addition of guanidine HCl. In contrast, the Gibbs free energy of unfolding is small, negative, and independent of the guanidine HCl concentration, contradicting, in part, the conclusions of the LEM. The unfolding enthalpy is compensated by an even larger entropy term. The multistate Zimm-Bragg theory predicts a larger conformational enthalpy and a smaller Gibbs free energy than the two-state model. The Zimm-Bragg theory provides the protein cooperativity parameter, the average length of independently folding protein domains, and the Gibbs free energy of unfolding of individual amino acid residues. Guanidine HCl binding to lysozyme is exothermic and counteracts the endothermic unfolding enthalpy. The number of bound denaturant molecules is determined from the decrease in enthalpy and is extrapolated to the guanidine HCl-to-amino acid stoichiometry at complete lysozyme unfolding. Chemical unfolding isotherms measured with circular dichroism (CD) spectroscopy are analyzed with both models. The chemical Zimm-Bragg theory is a cooperative molecular model, yielding the guanidine HCl binding constant and the protein cooperativity parameter. It allows a quantitative comparison between thermal and chemical protein unfolding. The two reactions have almost identical changes in Gibbs free energy. However, thermal unfolding is significantly more cooperative than chemical unfolding. Finally, distinct differences are observed in thermal unfolding between DSC and CD spectroscopy.

Cooperative protein unfolding. A statistical-mechanical model for the action of denaturants

Knowledge of protein stability is of utmost importance in various fields of biotechnology. Protein stability can be assessed in solution by increasing the concentration of denaturant and recording the structural changes with spectroscopic or thermodynamic methods. The standard interpretation of the experimental data is to assume a 2-state equilibrium between completely folded and completely unfolded protein molecules. Here we propose a cooperative model based on the statistical-mechanical Zimm-Bragg theory. In this model protein unfolding is driven by the weak binding of a rather small number of denaturant molecules, inducing the cooperative unfolding with multiple dynamic intermediates. The modified Zimm-Bragg theory is applied to published thermodynamic and spectroscopic data leading to the following conclusions. (i) The binding constant K_D is correlated with the midpoint concentration, c₀, of the unfolding reaction according to c₀≅1/K_D. The better the binding of denaturant the lower is the concentration to achieve unfolding. (ii) The binding constant K_D agrees with direct thermodynamic measurements. A rather small number of bound denaturants suffices to induce the cooperative unfolding of the whole protein. (iii) Chemical unfolding occurs in the concentration range Δc_D=c_end-c_ini. The theory predicts the unfolding energy per amino acid residue as g_nu=RTK_D(c_end-c_ini). The Gibbs free energy of an osmotic gradient of the same size is ΔG_Diff=-RTln(c_end/c_ini). In all examples investigated ΔG_Diff exactly balances the unfolding energy g_nu. The total unfolding energy is thus close to zero. (iv) Protein cooperativity in chemical unfolding is rather low with cooperativity parameters σ≥3x10^-3. As a consequence, the theory predicts a dynamic mixture of conformations during the unfolding reaction. The probabilities of individual conformations are easily accessible via the partition function Z(c_D,σ).

Long-Lived Intermediates in a Cooperative Two-State Folding Transition

Biomolecular folding often occurs through a cooperative two-state reactant ↔ product transition; the term cooperative does not convey that intermediate structures are nonexistent but rather that these states are not observable by existing experimental techniques. Because of this, few intermediates have been studied and characterized. Recently, ion mobility spectrometry (IMS) measurements revealed that the oligomer polyproline-13 (Pro13, which in propanol (PrOH) favors the right-handed helical PPI structure having adjacent pyrrolidine rings in a cis configuration) folds through six sequential long-lived intermediates as it converts to the all-trans-configured PPII structure that is favored in aqueous solutions. Here, we examine the PPI_PrOH → PPII_aq folding transition for a HisPro13 sequence, i.e., Pro13 having a single histidine residue added to the N-terminus. Remarkably, the IMS measurements show that, upon addition of histidine, all of the IMS peaks associated with intermediate structures disappear. Instead, HisPro13 folds via a cooperative two-state transition, delayed by a significant induction period. The induction period is temperature dependent-shifting the transition to longer times at lower temperatures. Equilibrium studies show that the HisPro13 PPI_PrOH → PPII_aq transition is endothermic but favored entropically. From these clues, we propose a sequential folding mechanism and develop a model that suggests that ∼13-17 long-lived intermediates are likely responsible for the induction period. In this model, intermediates are separated by average individual activation barriers of ∼90 kJ·mol^-1, and are entropically favorable.

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp(3) atom, and (iii) reduced 1-4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric-isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14-C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14-C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric-isothermal MD simulations 2-4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490-1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Functionally conservative substitutions at cardiac troponin I S43/45

A phospho-null Ala substitution at protein kinase C (PKC)-targeted cardiac troponin I (cTnI) S43/45 reduces myocyte and cardiac contractile function. The goal of the current study was to test whether cTnIS43/45N is an alternative, functionally conservative substitution in cardiac myocytes. Partial and more extensive endogenous cTnI replacement was similar at 2 and 4 days after gene transfer, respectively, for epitope-tagged cTnI and cTnIS43/45N. This replacement did not significantly change thin filament stoichiometry. In functional studies, there were no significant changes in the amplitude and/or rates of contractile shortening and re-lengthening after this partial (2 days) and extensive (4 days) replacement with cTnIS43/45N. The cTnIS43/45N substitution also was not associated with adaptive changes in the myocyte Ca(2+) transient or in phosphorylation of the protein kinase A and C-targeted cTnIS23/24 site. These results provide evidence that cTnIS43/45N is a functionally conservative substitution, and may be appropriate for use as a phospho-null in rodent models designed for studies on PKC modulation of cardiac performance.

Thermal protein unfolding by differential scanning calorimetry and circular dichroism spectroscopy Two-state model versus sequential unfolding

Thermally-induced protein unfolding is commonly described with the two-state model. This model assumes only two types of protein molecules in solution, the native (N) and the denatured, unfolded (U) protein. In reality, protein unfolding is a multistep process, even if intermediate states are only sparsely populated. As an alternative approach we explore the Zimm-Bragg theory, originally developed for the α-helix-to-random coil transition of synthetic polypeptides. The theory includes intermediate structures with concentrations determined by the cooperativity of the unfolding reaction. We illustrate the differences between the two-state model and the Zimm-Bragg theory with measurements of apolipoprotein A-1 and lysozyme by differential scanning calorimetry (DSC) and CD spectroscopy. Nine further protein examples are taken from the literature. The Zimm-Bragg theory provides a perfect fit of the calorimetric unfolding transitions for all proteins investigated. In contrast, the transition curves and enthalpies predicted by the two-state model differ considerably from the experimental results. Apolipoprotein A-1 is ~50% α-helical at ambient temperature and its unfolding follows the classical α-helix-to-random coil equilibrium. The unfolding of proteins with little α-helix content, such as lysozyme, can also be analyzed with the Zimm-Bragg theory by introducing the concept of 'folded' and 'unfolded' peptide units assuming an average unfolding enthalpy per peptide unit. DSC is the method of choice to measure the unfolding enthalpy, , but CD spectroscopy in combination with the two-state model is often used to deduce the unfolding enthalpy. This can lead to erroneous result. Not only are different enthalpies required to describe the CD and DSC transition curves but these values deviate distinctly from the experimental result. In contrast, the Zimm-Bragg theory predicts the DSC and CD unfolding transitions with the same set of parameters.

Hydrocarbon stapled peptides as modulators of biological function

Peptide-based drug discovery has experienced a significant upturn within the past decade since the introduction of chemical modifications and unnatural amino acids has allowed for overcoming some of the drawbacks associated with peptide therapeutics. Strengthened by such features, modified peptides become capable of occupying a niche that emerges between the two major classes of today's therapeutics-small molecules (<500 Da) and biologics (>5000 Da). Stabilized α-helices have proven particularly successful at impairing disease-relevant PPIs previously considered "undruggable." Among those, hydrocarbon stapled α-helical peptides have emerged as a novel class of potential peptide therapeutics. This review provides a comprehensive overview of the development and applications of hydrocarbon stapled peptides discussing the benefits and limitations of this technique.

Ensemble Calculation for Intrinsically Disordered Proteins Using NMR Parameters

Intrinsically disordered proteins (IDPs) perform their function despite their lack of well-defined tertiary structure. Residual structure has been observed in IDPs, commonly described as transient/dynamic or expressed in terms of fractional populations. In order to understand how the protein primary sequence dictates the dynamic and structural properties of IDPs and in general to understand how IDPs function, atomic-level descriptions are needed. Nuclear magnetic resonance spectroscopy provides information about local and long-range structure in IDPs at amino acid specific resolution and can be used in combination with ensemble descriptions to represent the dynamic nature of IDPs. In this chapter we describe sample-and-select approaches for ensemble modelling of local structural propensities in IDPs with specific emphasis on validation of these ensembles.

Dipole-mediated rectification of intramolecular photoinduced charge separation and charge recombination

Controlling charge transfer at a molecular scale is critical for efficient light harvesting, energy conversion, and nanoelectronics. Dipole-polarization electrets, the electrostatic analogue of magnets, provide a means for "steering" electron transduction via the local electric fields generated by their permanent electric dipoles. Here, we describe the first demonstration of the utility of anthranilamides, moieties with ordered dipoles, for controlling intramolecular charge transfer. Donor-acceptor dyads, each containing a single anthranilamide moiety, distinctly rectify both the forward photoinduced electron transfer and the subsequent charge recombination. Changes in the observed charge-transfer kinetics as a function of media polarity were consistent with the anticipated effects of the anthranilamide molecular dipoles on the rectification. The regioselectivity of electron transfer and the molecular dynamics of the dyads further modulated the observed kinetics, particularly for charge recombination. These findings reveal the underlying complexity of dipole-induced effects on electron transfer and demonstrate unexplored paradigms for molecular rectifiers.

Local motifs involved in the canonical structure of the ligand-binding domain in the nuclear receptor superfamily

Structural and sequence alignment analyses have revealed the existence of class-dependent and -independent local motifs involved in the overall fold of the ligand-binding domain (LBD) in the nuclear receptor (NR) superfamily. Of these local motifs, three local motifs, i.e., AF-2 fixed motifs, were involved in the agonist conformation of the activation function-2 (AF-2) region of the LBD. Receptor-agonist interactions increased the stability of these AF-2 fixed motifs in the agonist conformation. In contrast, perturbation of the AF-2 fixed motifs by a ligand or another protein molecule led the AF-2 architecture to adopt an antagonist conformation. Knowledge of this process should provide us with novel insights into the 'agonism' and 'antagonism' of NRs.

Stabilization and anomalous hydration of collagen fibril under heating

BACKGROUND

Type I collagen is the most common protein among higher vertebrates. It forms the basis of fibrous connective tissues (tendon, chord, skin, bones) and ensures mechanical stability and strength of these tissues. It is known, however, that separate triple-helical collagen macromolecules are unstable at physiological temperatures. We want to understand the mechanism of collagen stability at the intermolecular level. To this end, we study the collagen fibril, an intermediate level in the collagen hierarchy between triple-helical macromolecule and tendon.

METHODOLOGY/PRINCIPAL FINDING

When heating a native fibril sample, its Young's modulus decreases in temperature range 20-58°C due to partial denaturation of triple-helices, but it is approximately constant at 58-75°C, because of stabilization by inter-molecular interactions. The stabilization temperature range 58-75°C has two further important features: here the fibril absorbs water under heating and the internal friction displays a peak. We relate these experimental findings to restructuring of collagen triple-helices in fibril. A theoretical description of the experimental results is provided via a generalization of the standard Zimm-Bragg model for the helix-coil transition. It takes into account intermolecular interactions of collagen triple-helices in fibril and describes water adsorption via the Langmuir mechanism.

CONCLUSION/SIGNIFICANCE

We uncovered an inter-molecular mechanism that stabilizes the fibril made of unstable collagen macromolecules. This mechanism can be relevant for explaining stability of collagen.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Structural propensities and entropy effects in peptide helix-coil transitions

The helix-coil transition in peptides is a critical structural transition leading to functioning proteins. Peptide chains have a large number of possible configurations that must be accounted for in statistical mechanical investigations. Using hydrogen bond and local helix propensity interaction terms, we develop a method for obtaining and incorporating the degeneracy factor that allows the exact calculation of the partition function for a peptide as a function of chain length. The partition function is used in calculations for engineered peptide chains of various lengths that allow comparison with a variety of different types of experimentally measured quantities, such as fraction of helicity as a function of both temperature and chain length, heat capacity, and denaturation studies. When experimental sensitivity in helicity measurements is properly accounted for in the calculations, the calculated curves fit well with the experimental curves. We determine values of interaction energies for comparison with known biochemical interactions, as well as quantify the difference in the number of configurations available to an amino acid in a random coil configuration compared to a helical configuration.

Inhibition of amyloid aggregation by formation of helical assemblies

The formation of amyloid aggregates is responsible for a wide range of diseases, including Alzheimer's and Parkinson's disease. Although the amyloid-forming proteins have different structures and sequences, all undergo a conformational change to form amyloid aggregates that have a characteristic cross-β-structure. The mechanistic details of this process are poorly understood, but different strategies for the development of inhibitors of amyloid formation have been proposed. In most cases, chemically diverse compounds bind to an elongated form of the protein in a β-strand conformation and thereby exert their therapeutic effect. However, this approach could favor the formation of prefibrillar oligomeric species, which are thought to be toxic. Herein, we report an alternative approach in which a helical coiled-coil-based inhibitor peptide has been designed to engage a coiled-coil-based amyloid-forming model peptide in a stable coiled-coil arrangement, thereby preventing rearrangement into a β-sheet conformation and the subsequent formation of amyloid-like fibrils. Moreover, we show that the helix-forming peptide is able to disassemble mature amyloid-like fibrils.

Permanent electric dipole moments of carboxyamides in condensed media: what are the limitations of theory and experiment?

Electrostatic properties of proteins are crucial for their functionality. Carboxyamides are small polar groups that, as peptide bonds, are principal structural components of proteins that govern their electrostatic properties. We investigated the medium dependence of the molar polarization and of the permanent dipole moments of amides with different state of alkylation. The experimentally measured and theoretically calculated dipole moments manifested a solvent dependence that increased with the increase in the media polarity. We ascribed the observed enhancement of the amide polarization to the reaction fields in the solvated cavities. Chloroform, for example, caused about a 25% increase in the amide dipole moments determined for vacuum, as the experimental and theoretical results demonstrated. Another chlorinated solvent, 1,1,2,2-tetrachloroethane, however, caused an "abnormal" increase in the experimentally measured amide dipoles, which the theoretical approaches we used could not readily quantify. We showed and discussed alternatives for addressing such discrepancies between theory and experiment.

The contribution of entropy, enthalpy, and hydrophobic desolvation to cooperativity in repeat-protein folding

Cooperativity is a defining feature of protein folding, but its thermodynamic and structural origins are not completely understood. By constructing consensus ankyrin repeat protein arrays that have nearly identical sequences, we quantify cooperativity by resolving stability into intrinsic and interfacial components. Heteronuclear NMR and CD spectroscopy show that these constructs adopt ankyrin repeat structures. Applying a one-dimensional Ising model to a series of constructs chosen to maximize information content in unfolding transitions, we quantify stabilities of the terminal capping repeats, and resolve the effects of denaturant into intrinsic and interfacial components. Reversible thermal denaturation resolves interfacial and intrinsic free energies into enthalpic, entropic, and heat capacity terms. Intrinsic folding is entropically disfavored, whereas interfacial interaction is entropically favored and attends a decrease in heat capacity. These results suggest that helix formation and backbone ordering occurs upon intrinsic folding, whereas hydrophobic desolvation occurs upon interfacial interaction, contributing to cooperativity.

Direct assessment of the α-helix nucleation time

The nucleation event in α-helix formation is a fundamental process in protein folding. However, determining how quickly it takes place based on measurements of the relaxation dynamics of helical peptides is difficult because such relaxations invariably contain contributions from various structural transitions such as from helical to nonhelical states and helical to partial-helical conformations. Herein, we measure the temperature-jump (T-jump) relaxation kinetics of three model peptides that fold into a single-turn α-helix, using time-resolved infrared spectroscopy, aiming to provide a direct assessment of the helix nucleation rate. The α-helical structure of these peptides is stabilized by a covalent cross-linker formed between the side chains of two residues at the i and i + 4 positions. If we assume that this cross-linker mimics the structural constraint arising from a strong side chain-side chain interaction (e.g., a salt bridge) in proteins, these peptides would represent good models for studying the nucleation process of an α-helix in a protein environment. Indeed, we find that the T-jump induced relaxation rate of these peptides is approximately (0.6 μs)(-1) at room temperature, which is slower than that of commonly studied alanine-based helical peptides but faster than that of a naturally occurring α-helix whose folded state is stabilized by a series of side chain-side chain interactions. Taken together, our results put an upper limit of about 1 μs for the helix nucleation time at 20 °C and suggest that the subsequent propagation steps occur with a time constant of about 240 ns.

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.

Solid-phase synthesis of tris-benzamides as α-helix mimetics

Small molecules mimicking α-helices are of great interest since numerous protein-protein interactions use helical structures at the interface. With a goal of generating libraries of α-helix mimetics, an efficient solid-phase synthetic method was developed to produce tris-benzamides. The tris-benzamide scaffold was designed to place three side-chain functional groups found at the i, i + 4, and i + 7 positions of an α-helix, emulating one helical face. The synthetic strategy involves simple and iterative reactions of removal of an allyl ester, formation of an amide bond via an O → N acyl migration, and an O-alkylation. A small library of twenty tris-benzamides containing a variety of functional groups was prepared in high purity (83-99%) to demonstrate the versatility of the synthetic approach. This methodology allowed the facile and rapid construction of α-helix mimetics that would facilitate the identification of small molecules for target proteins.

Multivalent protein polymers with controlled chemical and physical properties

In this review, we describe our work on the design, characterization, and modification of a series of alanine-rich helical polypeptides with novel functions. Glycosylation of the polypeptides has permitted investigation of polymer architecture effects on multivalent interactions. One of the members of this polypeptide family exhibits polymorphological behavior that is easily manipulated via simple changes in solution pH and temperature. Polypeptide-based fibrils formed at acidic pH and high temperature were shown to direct the one-dimensional organization of gold nanoparticles via electrostatic interactions. As a precursor to fibrils, aggregates likely comprising alanine-rich cores form at low temperatures and acidic pH and reversibly dissociate into monomers upon deprotonation. PEGylation of these polypeptides does not alter the self-association or conformational behavior of the polypeptide, suggesting potential applications in the development of assembled delivery vehicles, as modification of the polypeptides should be a useful strategy for controlling assembly.

Two intracellular helices of G-protein coupling receptors could generally support oligomerization and coupling with transducers

For many G-protein coupling receptors (GPCRs), the upkeep of receptor dimers could depend on association with functional Gi α subunits. This is known for Y1, Y2 and Y4 neuropeptide Y receptors [presented in the companion paper (Estes et al., Amino Acids, doi: 10.1007/s00726-010-0642-z , 2010)]. Interactions with transducers use mainly intracellular domains of the receptors. Intracellular loops 1 and 2 in GPCRs are short and lack extensive helicity that could support transducer anchoring. Interaction with G-proteins is known to use the juxtamembrane Helix 8 in the fourth intracellular domain, for which we document a helix-stabilizing n/(n + 4) pattern of large hydrophobic sidechains. Another intracellular helix located in the C-terminal portion of the third intracellular loop does not display a strong stabilizing pattern, and is found in many studies to serve dynamically in association and activation of transducers and effectors. We show that these tracts share features across metazoan phyla not only in opsins and opsin-like receptors (including the Y receptors), but also in Taste-2 and Frizzled receptors. Similarities of these helices across GPCR groups could have both phylogenetic and functional roots.

Exactly solvable model for helix-coil-sheet transitions in protein systems

In view of the important role helix-sheet transitions play in protein aggregation, we introduce a simple model to study secondary structural transitions of helix-coil-sheet systems using a Potts model starting with an effective Hamiltonian. This energy function depends on four parameters that approximately describe entropic and enthalpic contributions to the stability of a polypeptide in helical and sheet conformations. The sheet structures involve long-range interactions between residues which are far in sequence, but are in contact in real space. Such contacts are included in the Hamiltonian. Using standard statistical mechanical techniques, the partition function is solved exactly using transfer matrices. Based on this model, we study thermodynamic properties of polypeptides, including phase transitions between helix, sheet, and coil structures.

Helix-coil energetics for helix formers and breakers reflect context and temperature: mutants of helically robust, guest-sensitive homopeptide hosts

The natural amino acids are primarily helix breakers at the low assignment temperatures characteristic of many studies, but recent genomic analyses of thermophilic proteins suggest that at high temperatures, some breakers may become strong helix formers. Moreover, the breaker/former inventory has not been previously characterized at the physiologically relevant temperature of 37 degrees C. The versatility of 13C==O NMR chemical shifts as helicity reporters allows construction of two mutant peptide series, tailored to expand the range of temperature assignments for helical propensities and derived from the core hosts tL-Ala9XxxAla9-tL and tL-AlaNva4XxxNva4Ala9-tL, Nva=norvaline. For three limiting guests Xxx, the helix former Nva and the breakers Gly and Pro, we report wXxx[T] assignments at seven temperatures from 2 to 80 degrees C, validating our reasoning and paving the way for assignment of a definitive wXxx[T] data-base.

Circular dichroism and UV resonance raman study of the impact of alcohols on the Gibbs free energy landscape of an alpha-helical peptide

We used CD and UV resonance Raman spectroscopy to study the impact of alcohols on the conformational equilibria and relative Gibbs free energy landscapes along the Ramachandran Psi-coordinate of a mainly poly-Ala peptide, AP with an AAAAA(AAARA)(3)A sequence. 2,2,2-Trifluoroethanol (TFE) most stabilizes the alpha-helix-like conformations, followed by ethanol, methanol, and pure water. The pi-bulge conformation is stabilized more than the alpha-helix, while the 3(10)-helix is destabilized due to the alcohol-increased hydrophobicity. Turns are also stabilized by alcohols. We also found that while TFE induces more alpha-helices, it favors multiple, shorter helix segments.

Residue patterning in helix interiorsThis paper is one of a selection of papers published in this special issue entitled “Canadian Society of Biochemistry, Molecular & Cellular Biology 52nd Annual Meeting — Protein Folding: Principles and Diseases” and has undergone the Journal's usual peer review process

The α-helix remains a focus of research because of its importance to protein folding and structure. Nevertheless, despite numerous empirical, computational, and theoretical studies, the fundamental structural properties governing their formation and stability are still unclear. We have examined the statistical occurrence of polar and apolar residue patterning in helical interiors in a large, non-redundant dataset, and compared these patterns with those found in other structural environments. While the familiar amphipathic distributions for both polar and apolar residues are evident, our analysis also finds significant differences between these distributions. Non-amphipathic signals can also be discerned within both distributions. Most interestingly, among various positional patterning, an analysis of immediate (i,i + 1) helical neighbours found: (i) clear neighbouring preferences, with high (low) occurrences of hydrophobics (hydrophilics) next to Gly, Pro, and short polar residues; (ii) high negative (positive) correlation between residue helical propensities and the degree of neighbouring hydrophobicity (hydrophilicity); and (iii) a preferred ordering among neighbours, implying an inherent helix directionality. Because (i,i + 1) helical pairs have limited side chain – side chain interactions, thermodynamic considerations cannot readily explain these observations, nor can evolutionary pressures that enhance tertiary interactions via amphipathicity, as this particular spacing does not segregate residues onto either the same or opposing helical faces. We suggest that the mechanism of helix formation may be partly responsible for these observations. In particular, the high negative correlation between residue helical propensities and neighbouring hydrophobicity suggests that hydrophobicity may play a more important role in helix formation than currently recognized.

Context-independent, temperature-dependent helical propensities for amino acid residues

Assigned from data sets measured in water at 2, 25, and 60 degrees C containing (13)C=O NMR chemical shifts and [theta](222) ellipticities, helical propensities are reported for the 20 genetically coded amino acids, as well as for norvaline and norleucine. These have been introduced by chemical synthesis at central sites within length-optimized, spaced, solubilized Ala(19) hosts. The resulting polyalanine-derived, quantitative propensity sets express for each residue its temperature-dependent but context-independent tendency to forego a coil state and join a preexisting helical conformation. At 2 degrees C their rank ordering is: P << G < H < C, T, N < S < Y, F, W < V, D < K < Q < I < R, M < L < E < A; at 60 degrees C the rank becomes: H, P < G < C < R, K < T, Y, F < N, V < S < Q < W, D < I, M < E < A < L. The DeltaDeltaG values, kcal/mol, relative to alanine, for the cluster T, N, S, Y, F, W, V, D, Q, imply that at 2 degrees C all are strong breakers: DeltaDeltaG(mean) = +0.63 +/- 0.11, but at 60 degrees C their breaking tendencies are dramatically attenuated and converge toward the mean: DeltaDeltaG(mean) = +0.25 +/- 0.07. Accurate modeling of helix-rich proteins found in thermophiles, mesophiles, and organisms that flourish near 0 degrees C thus requires appropriately matched propensity sets. Comparisons are offered between the temperature-dependent propensity assignments of this study and those previously assigned by the Scheraga group; the special problems that attend propensity assignments for charged residues are illustrated by lysine guest data; and comparisons of errors in helicity assignments from shifts and ellipticity data show that the former provide superior precision and accuracy.

Structure and stability of helices in square-well homopolymers

Recently, it has been demonstrated [Magee, Phys. Rev. Lett. 96, 207802 (2006)] that isolated square-well homopolymers can spontaneously break chiral symmetry and "freeze" into helical structures at sufficiently low temperatures. This behavior is interesting because the square-well homopolymer is itself achiral. In this work, we use event-driven molecular dynamics combined with an optimized parallel tempering scheme to study this polymer model over a wide range of parameters. We examine the conditions where the helix structure is stable and determine how the interaction parameters of the polymer govern the details of the helix structure. The width of the square well (proportional to lambda) is found to control the radius of the helix, which decreases with increasing well width until the polymer forms a coiled sphere for sufficiently large wells. The helices are found to be stable for only a "window" of molecular weights. If the polymer is too short, the helix will not form. If the polymer is too long, the helix is no longer the minimum energy structure, and other folded structures will form. The size of this window is governed by the chain stiffness, which in this model is a function of the ratio of the monomer size to the bond length. Outside this window, the polymer still freezes into a locked structure at low temperature; however, unless the chain is sufficiently stiff, this structure will not be unique and is similar to a glassy state.

Folding by numbers: primary sequence statistics and their use in studying protein folding

The exponential growth over the past several decades in the quantity of both primary sequence data available and the number of protein structures determined has provided a wealth of information describing the relationship between protein primary sequence and tertiary structure. This growing repository of data has served as a prime source for statistical analysis, where underlying relationships between patterns of amino acids and protein structure can be uncovered. Here, we survey the main statistical approaches that have been used for identifying patterns within protein sequences, and discuss sequence pattern research as it relates to both secondary and tertiary protein structure. Limitations to statistical analyses are discussed, and a context for their role within the field of protein folding is given. We conclude by describing a novel statistical study of residue patterning in β-strands, which finds that hydrophobic (i,i+2) pairing in β-strands occurs more often than expected at locations near strand termini. Interpretations involving β-sheet nucleation and growth are discussed.

Statistical thermodynamics of the collagen triple-helix/coil transition. Free energies for amino acid substitutions within the triple-helix

Collagen sequences frequently deviate from the most thermally stable (Gly-Pro-Hyp)(n) pattern, with many mutations causing osteogenesis imperfecta (or "brittle bone disease"). The effects of collagen mutations have been studied in short peptides. The analysis of this work is problematic, however, as triple-helices fray from their ends, making the coil/triple-helix equilibrium non-two-state. Here, I develop a statistical thermodynamic model to handle this equilibrium that is applicable to peptides that follow the (G-X-Y)(n) pattern, where Gly is present at every third position and where all three chains are identical. Parameters for substitutions at each position are included, as well as a penalty for initiating triple-helix formation. The model is applied to equilibrium experimental data at 37 degrees C to show that the extension of a triple-helix by a three residue unit stabilizes the triple-helix by 0.76 kcal/mol for Gly-Pro-Hyp and 0.33 kcal/mol for Gly-Pro-Pro. The replacement of Hyp by Arg, Asp, or Trp destabilizes the triple-helix by 1.5, 2.4, and 2.9 kcal/mol, respectively, where the substitution is present once in each chain. The model can thus be used to quantitatively interpret data on collagen peptides, giving free energies that can help rationalize mutations that affect collagen stability, and to design new collagen sequences.