Alpha-helix stability in proteins. II. Factors that influence stability at an internal position

Amino acid propensities for secondary structures and its variation across protein structures using exhaustive PDB data

Amino acid propensities for protein secondary structures are vital for protein structure prediction, understanding folding, and design, and have been studied using various theoretical and experimental methods. Traditional assessments of average propensities using statistical methods have been done on relatively smaller dataset for only a few secondary structures. They also involve averaging out the environmental factors and lack insights into consistency of preferences across diverse protein structures. While a few studies have explored variations in propensities across protein structural classes and folds, exploration of such variations across protein structures remains to be carried out. In this work, we have revised the average propensities for all six different secondary structures, namely α-helix, β-strand, 310-helix, π-helix, turn and coil, analyzing the most exhaustive dataset available till date using two robust secondary structure assignment algorithms, DSSP and STRIDE. The propensities evaluated here can serve as a standard reference. Moreover, we present here, for the first time, the propensities within individual protein structures and investigated how the preferences of residues and more interestingly, of their groups formed based on their structural features, vary across different unique structures. We devised a novel approach- the minimal set analysis, based on the propensity distribution of residues, which along with the group propensities led us to the conclusion that a residue's preference for a specific secondary structure is primarily dictated by its side chain's structural features. The findings in this study provide a more insightful picture of residues propensities and can be useful in protein folding and design studies.

MEnTaT: A machine-learning approach for the identification of mutations to increase protein stability

Enhancing protein thermal stability is important for biomedical and industrial applications as well as in the research laboratory. Here, we describe a simple machine-learning method which identifies amino acid substitutions that contribute to thermal stability based on comparison of the amino acid sequences of homologous proteins derived from bacteria that grow at different temperatures. A key feature of the method is that it compares the sequences based not simply on the amino acid identity, but rather on the structural and physicochemical properties of the side chain. The method accurately identified stabilizing substitutions in three well-studied systems and was validated prospectively by experimentally testing predicted stabilizing substitutions in a polyamine oxidase. In each case, the method outperformed the widely used bioinformatic consensus approach. The method can also provide insight into fundamental aspects of protein structure, for example, by identifying how many sequence positions in a given protein are relevant to temperature adaptation.

Circular permutation at azurin's active site slows down its folding

Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

Propensities of Some Amino Acid Pairings in α-Helices Vary with Length

The results of secondary structure prediction methods are widely used in applications in biotechnology and bioinformatics. However, the accuracy limit of these methods could be improved up to 92%. One approach to achieve this goal is to harvest information from the primary structure of the peptide. This study aims to contribute to this goal by investigating the variations in propensity of amino acid pairings to α-helices in globular proteins depending on helix length. (n):(n + 4) residue pairings were determined using a comprehensive peptide data set according to backbone hydrogen bond criterion which states that backbone hydrogen bond is the dominant driving force of protein folding. Helix length is limited to 13 to 26 residues. Findings of this study show that propensities of ALA:GLY and GLY:GLU pairings to α-helix in globular protein increase with increasing helix length but of ALA:ALA and ALA:VAL decrease. While the frequencies of ILE:ALA, LEU:ALA, LEU:GLN, LEU:GLU, LEU:LEU, MET:ILE and VAL:LEU pairings remain roughly constant with length, the 25 residue pairings have varying propensities in narrow helix lengths. The remaining pairings have no prominent propensity to α-helices.

Advances in Cryochemistry: Mechanisms, Reactions and Applications

It is counterintuitive that chemical reactions can be accelerated by freezing, but this amazing phenomenon was discovered as early as the 1960s. In frozen systems, the increase in reaction rate is caused by various mechanisms and the freeze concentration effect is the main reason for the observed acceleration. Some accelerated reactions have great application value in the chemistry synthesis and environmental fields; at the same time, certain reactions accelerated at low temperature during the storage of food, medicine, and biological products should cause concern. The study of reactions accelerated by freezing will overturn common sense and provide a new strategy for researchers in the chemistry field. In this review, we mainly introduce various mechanisms for accelerating reactions induced by freezing and summarize a variety of accelerated cryochemical reactions and their applications.

Learning the changes of barnase mutants thermostability from structural fluctuations obtained using anisotropic network modeling

In biotechnology applications, rational design of new proteins with improved physico-chemical properties includes a number of important tasks. One of the greatest practical and fundamental challenges is the design of highly thermostable protein enzymes that maintain catalytic activity at high temperatures. This problem may be solved by introducing mutations into the wild-type enzyme protein. In this work, to predict the impact of such mutations in barnase protein we applied the anisotropic network modeling approach, revealing atomic fluctuations in structural regions that are changed in mutants compared to the wild-type protein. A regression model was constructed based on these structural features that can allow one to predict the thermal stability of new barnase mutants. Moreover, the analysis of regression model provides a mechanistic explanation of how the structural features can contribute to the thermal stability of barnase mutants.

Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields

Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.

Uncovering dehydration in cytochrome c refolding from urea- and guanidine hydrochloride-denatured unfolded state by high pressure spectroscopy

To investigate the dehydration associated with protein folding, the partial molar volume changes for protein unfolding (ΔV_u) in cytochrome c (Cyt c) were determined using high pressure absorption spectroscopy. ΔV_u values for the unfolding to urea- and guanidine hydrochloride (GdnHCl)-denatured Cyt c were estimated to be 56±5 and 29±1 mL mol⁻¹, respectively. Considering that the volume change for hydration of hydrophobic groups is positive and that Cyt c has a covalently bonded heme, a positive ΔV_u reflects the primary contribution of the hydration of heme. Because of the marked tendency of guanidium ions to interact with hydrophobic groups, a smaller number of water molecules were hydrated with hydrophobic groups in GdnHCl-denatured Cyt c than in urea-denatured Cyt c, resulting in the smaller positive ΔV_u. On the other hand, urea is a relatively weak denaturant and urea-denatured Cyt c is not completely hydrated, which retains the partially folded structures. To unfold such partial structures, we introduced a mutation near the heme binding site, His26, to Gln, resulting in a negatively shifted ΔV_u (4±2 mL mol⁻¹) in urea-denatured Cyt c. The formation of the more solvated and less structured state in the urea-denatured mutant enhanced hydration to the hydrophilic groups in the unfolding process. Therefore, we confirmed the hydration of amino acid residues in the protein unfolding of Cyt c by estimating ΔV_u, which allows us to discuss the hydrated structures in the denatured states of proteins.

Comparison of Peptide Ion Conformers Arising from Non-Helical and Helical Peptides Using Ion Mobility Spectrometry and Gas-Phase Hydrogen/Deuterium Exchange

The dominant gas-phase conformer of [M+3H]³⁺ ions of the model peptide acetyl-PSSSSKSSSSKSSSSKSSSSK has been examined with ion mobility spectrometry (IMS), gas-phase hydrogen deuterium exchange (HDX), and mass spectrometry (MS) techniques. The [M+3H]³⁺ peptide ions are observed predominantly as a relatively compact conformer type. Upon subjecting these ions to electron transfer dissociation (ETD), the level of protection for each amino acid residue in the peptide sequence is assessed. The overall per-residue deuterium uptake is observed to be relatively more efficient for the neutral residues than for the model peptide acetyl-PAAAAKAAAAKAAAAKAAAAK. In comparison, the N-terminal and C-terminal regions of the serine peptide show greater relative protection compared with interior residues. Molecular dynamics (MD) simulations have been used to generate candidate structures for collision cross section and HDX reactivity matching. Hydrogen accessibility scoring (HAS) for select structural candidates from MD simulations has been used to suggest conformer types that could contribute to the observed HDX patterns. The results are discussed with respect to recent studies employing extensive MD simulations of gas-phase structure establishment of a peptide system. Graphical Abstract ᅟ.

Mutations and functional analysis of 14-3-3 stress response protein from Triticum aestivum: An evolutionary analysis through in silico structural biochemistry approach

Wheat (Triticum aestivum), having high nutritional values is one of the staple food of most of the countries in the world. The productivity of the crop decreases drastically when it encounters various abiotic stresses, most common of which are heat, drought, flood and salinity. There is a crucial role of stress response proteins for the survival of the crops in stress conditions. So the study of wheat stress response proteins is of great importance to raise wheat production in different stress conditions. In this study, we analysed 14-3-3 protein, a stress response protein that is expressed in three major stresses, for example heat, drought and salinity and helps the plants to survive in those conditions. Effect of mutations in the 14-3-3 sequence was predicted using its domain, secondary structure and multiple sequence alignment of amino acid sequences from wheat and its related species. The functional diversity of the protein in different species was correlated with mutations, change in secondary structure and the evolutionary relatedness of the protein in different species. This is the first novel work for analysing the mutational effect on the structure and function of a stress response protein (14-3-3) from Triticum aestivum and its related species.

Structural Capacitance in Protein Evolution and Human Diseases

Canonical mechanisms of protein evolution include the duplication and diversification of pre-existing folds through genetic alterations that include point mutations, insertions, deletions, and copy number amplifications, as well as post-translational modifications that modify processes such as folding efficiency and cellular localization. Following a survey of the human mutation database, we have identified an additional mechanism that we term "structural capacitance," which results in the de novo generation of microstructure in previously disordered regions. We suggest that the potential for structural capacitance confers select proteins with the capacity to evolve over rapid timescales, facilitating saltatory evolution as opposed to gradualistic canonical Darwinian mechanisms. Our results implicate the elements of protein microstructure generated by this distinct mechanism in the pathogenesis of a wide variety of human diseases. The benefits of rapidly furnishing the potential for evolutionary change conferred by structural capacitance are consequently counterbalanced by this accompanying risk. The phenomenon of structural capacitance has implications ranging from the ancestral diversification of protein folds to the engineering of synthetic proteins with enhanced evolvability.

Thermodynamic Scale of β-Amino Acid Residue Propensities for an α-Helix-like Conformation

A thiol-thioester exchange system has been used to measure the propensities of diverse β-amino acid residues to participate in an α-helix-like conformation. These measurements depend on formation of a parallel coiled-coil tertiary structure when two peptide segments become linked by thioester formation. One peptide segment contains a "guest" site that accommodates diverse β residues and is distal to the coiled-coil interface. We find that helix propensity is influenced by side chain placement within the β residue [β3 (side chain adjacent to nitrogen) slightly favored relative to β2 (side chain adjacent to carbonyl)]. The previously recognized helix stabilization resulting from five-membered ring incorporation is quantified. These results are significant because so few quantitative thermodynamic measurements have been reported for α/β-peptide folding.

Helix Propensities of Amino Acid Residues via Thioester Exchange

We describe the use of thioester exchange equilibria to measure the propensities of amino acid residues to participate in helical secondary structure at room temperature in the absence of denaturants. Thermally or chemically induced unfolding has previously been employed to measure α-helix propensities among proteinogenic α-amino acid residues, and quantitative comparison with precedents indicates that the thioester exchange system is reliable for residues that lack side chain charge. This system allows the measurement of α-helix propensities for d-α-amino acid residues and propensities of residues with nonproteinogenic backbones, such as those derived from a β-amino acid, to participate in an α-helix-like secondary structure.

Discovering functional, non-proteinogenic amino acid containing, peptides using genetic code reprogramming

The protein synthesis machinery of the cell, the ribosome and associated factors, is able to accurately follow the canonical genetic code, that which maps RNA sequence to protein sequence, to assemble functional proteins from the twenty or so proteinogenic amino acids. A number of innovative methods have arisen to take advantage of this accurate, and efficient, machinery to direct the assembly of non-proteinogenic amino acids. We review and compare these routes to 'reprogram the genetic code' including in vitro translation, engineered aminoacyl tRNA synthetases, and RNA 'flexizymes'. These studies show that the ribosome is highly tolerant of unnatural amino acids, with hundreds of unusual substrates of varying structure and chemistries being incorporated into protein chains. We also discuss how these methods have been coupled to selection techniques, such as phage display and mRNA display, opening up an exciting new avenue for the production of proteins and peptides with properties and functions beyond that which is possible using proteins composed entirely of the proteinogenic amino acids.

Stabilization of G protein-coupled receptors by point mutations

G protein-coupled receptors (GPCRs) are flexible integral membrane proteins involved in transmembrane signaling. Their involvement in many physiological processes makes them interesting targets for drug development. Determination of the structure of these receptors will help to design more specific drugs, however, their structural characterization has so far been hampered by the low expression and their inherent instability in detergents which made protein engineering indispensable for structural and biophysical characterization. Several approaches to stabilize the receptors in a particular conformation have led to breakthroughs in GPCR structure determination. These include truncations of the flexible regions, stabilization by antibodies and nanobodies, fusion partners, high affinity and covalently bound ligands as well as conformational stabilization by mutagenesis. In this review we focus on stabilization of GPCRs by insertion of point mutations, which lead to increased conformational and thermal stability as well as improved expression levels. We summarize existing mutagenesis strategies with different coverage of GPCR sequence space and depth of information, design and transferability of mutations and the molecular basis for stabilization. We also discuss whether mutations alter the structure and pharmacological properties of GPCRs.

Pairwise amino acid secondary structural propensities

We investigate the propensities for amino acids to form a specific secondary structure when they are paired with other amino acids. Our investigations use molecular dynamics (MD) computer simulations, and we compare the results to those from the Protein Data Bank (PDB). Proper comparison requires weighting of the MD results in a manner consistent with the relative frequency of appearance in the PDB of each possible pair of amino acids. We find that the propensity for an amino acid to assume a secondary structure varies dramatically depending on the amino acid that is before or after it in the primary sequence. This cooperative effect means that when selecting amino acids to facilitate the formation of a secondary structure in peptide engineering experiments, the adjacent amino acids must be considered.

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

Coupled folding and binding of the disordered protein PUMA does not require particular residual structure

Many cellular proteins are 'disordered' in isolation. A subset of these intrinsically disordered proteins (IDPs) can, upon binding another molecule, fold to a well-defined three-dimensional structure. In the structurally heterogeneous, unbound ensemble of these IDPs, conformations are likely to exist that, in part, resemble the final bound form. It has been suggested that these conformations, displaying 'residual structure', could be important for the mechanism of such coupled folding and binding reactions. PUMA, of the BCL-2 family, is an IDP in isolation but will form a single, contiguous α-helix upon binding the folded protein MCL-1. Using the helix-breaking residue proline, we systematically target each potential turn of helix of unbound PUMA and assess the binding to MCL-1 using time-resolved stopped-flow techniques. All proline-containing mutants bound, and although binding was weaker than the wild-type protein, association rate constants were largely unaffected. We conclude that population of particular residual structure, containing a specific helical turn, is neither required for the binding nor for fast association of PUMA and MCL-1.

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.

Unraveling the role of the C-terminal helix turn helix of the coat-binding domain of bacteriophage P22 scaffolding protein

Many viruses encode scaffolding and coat proteins that co-assemble to form procapsids, which are transient precursor structures leading to progeny virions. In bacteriophage P22, the association of scaffolding and coat proteins is mediated mainly by ionic interactions. The coat protein-binding domain of scaffolding protein is a helix turn helix structure near the C terminus with a high number of charged surface residues. Residues Arg-293 and Lys-296 are particularly important for coat protein binding. The two helices contact each other through hydrophobic side chains. In this study, substitution of the residues of the interface between the helices, and the residues in the β-turn, by aspartic acid was used examine the importance of the conformation of the domain in coat binding. These replacements strongly affected the ability of the scaffolding protein to interact with coat protein. The severity of the defect in the association of scaffolding protein to coat protein was dependent on location, with substitutions at residues in the turn and helix 2 causing the most significant effects. Substituting aspartic acid for hydrophobic interface residues dramatically perturbs the stability of the structure, but similar substitutions in the turn had much less effect on the integrity of this domain, as determined by circular dichroism. We propose that the binding of scaffolding protein to coat protein is dependent on angle of the β-turn and the orientation of the charged surface on helix 2. Surprisingly, formation of the highly complex procapsid structure depends on a relatively simple interaction.

Dependence of α-helical and β-sheet amino acid propensities on the overall protein fold type

Background

A large number of studies have been carried out to obtain amino acid propensities for α-helices and β-sheets. The obtained propensities for α-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the β-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects β-sheet propensity.

Results

We calculated amino acid propensities for α-helices and β-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that α-helix propensities do not differ significantly by fold, but β-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for α-helix, but such is not the case for the β-sheet propensities. We also found some fold dependence on amino acid frequency in β-strands. Folds with a high Ser, Thr and Asn content at exposed sites in β-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = −0.90) and to have flat β-sheets. At buried sites in β-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = −0.93). "All-β" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "α/β" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The α-helix propensities are similar for all folds and for exposed and buried residues. However, β-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in β-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a β-strand or association between two β sheets.

Nucleation effects in peptide foldamers

Oligomers composed of β(3)-amino acid residues and a mixture of α- and β(3)-residues have emerged as proteolytically stable structural mimics of α-helices. An attractive feature of these oligomers is that they adopt defined conformations in short sequences. In this manuscript, we evaluate the impact of β(3)-residues as compared to their α-amino acid analogs in prenucleated helices. Our hydrogen-deuterium exchange results suggest that heterogeneous sequences composed of "αααβ" repeats are conformationally more rigid than the corresponding homogeneous α-peptide helices, with the macrocycle templating the helical conformation having a significant influence.

New Soft-Core Potential Function for Molecular Dynamics Based Alchemical Free Energy Calculations

The fields of rational drug design and protein engineering benefit from accurate free energy calculations based on molecular dynamics simulations. A thermodynamic integration scheme is often used to calculate changes in the free energy of a system by integrating the change of the system's Hamiltonian with respect to a coupling parameter. These methods exploit nonphysical pathways over thermodynamic cycles involving particle introduction and annihilation. Such alchemical transitions require the modification of the classical nonbonded potential energy terms by applying soft-core potential functions to avoid singularity points. In this work, we propose a novel formulation for a soft-core potential to be applied in nonequilibrium free energy calculations that alleviates singularities, numerical instabilities, and additional minima in the potential energy for all combinations of nonbonded interactions at all intermediate alchemical states. The method was validated by application to (a) the free energy calculations of a closed thermodynamic cycle, (b) the mutation influence on protein thermostability, (c) calculations of small ligand solvation free energies, and (d) the estimation of binding free energies of trypsin inhibitors. The results show that the novel soft-core function provides a robust and accurate general purpose solution to alchemical free energy calculations.

Anti-HLA-E mAb 3D12 mimics MEM-E/02 in binding to HLA-B and HLA-C alleles: Web-tools validate the immunogenic epitopes of HLA-E recognized by the antibodies

HLA-E shares several peptide sequences with HLA-class Ia molecules. Therefore, anti-HLA-E antibodies that recognize the shared sequences may bind to HLA-class Ia alleles. This hypothesis was validated with a murine anti-HLA-E monoclonal antibody (mAb) MEM-E/02, which reacted with microbeads coated with several HLA-B and HLA-C antigens. In this report, the hypothesis was reexamined with another mAb 3D12, considered to be specific for HLA-E. The antibody binding is evaluated by measuring mean fluorescence index [MFI] with Luminex Multiplex Flow-Cytometric technology. The peptide-inhibition experiments are carried out with synthetic shared peptides, most prevalent to HLA-E and HLA-Ia alleles. The results showed that mAb 3D12 simulated MEM-E/02 in recognizing several HLA-B and HLA-C antigens. Both 3D12 and MEM-E/02 did not bind to HLA-A, HLA-F and HLA-G molecules. As observed with MEM-E/02, binding of 3D12 to HLA-E is inhibited by the peptides sequences (115)QFAYDGKDY(123) and (137)DTAAQI(142). Decrease in binding of mAb 3D12 to HLA class Ia, after heat treatment of antigen coated microbeads, supports the contention that the epitope may be located at the outside of the "thermodynamically stable" α-helix conformations of HLA-E. Several sequence and structure-based web-tools were employed to validate the discontinuous epitopes recognized by the mAbs. The scores obtained by these web-tools distinguished the shared peptide sequences that inhibited the mAb binding to HLA-E. Furthermore, ElliPro web tool points out that both mAbs recognize the conformational discontinuous epitopes (the shared inhibitory peptide sequences) in the secondary structure of the HLA-E molecule. The study favors the contention that the domain of the shared inhibitory peptide sequences may be the most immunogenic site of HLA-E molecule. It also postulates and clarifies that amino acid substitution on or near the binding domains may account for the lack of cross reactivity of 3D12 and MEM-E/02 with HLA-A, HLA-F and HLA-G molecules.

Protein thermostability calculations using alchemical free energy simulations

Thermal stability of proteins is crucial for both biotechnological and therapeutic applications. Rational protein engineering therefore frequently aims at increasing thermal stability by introducing stabilizing mutations. The accurate prediction of the thermodynamic consequences caused by mutations, however, is highly challenging as thermal stability changes are caused by alterations in the free energy of folding. Growing computational power, however, increasingly allows us to use alchemical free energy simulations, such as free energy perturbation or thermodynamic integration, to calculate free energy differences with relatively high accuracy. In this article, we present an automated protocol for setting up alchemical free energy calculations for mutations of naturally occurring amino acids (except for proline) that allows an unprecedented, automated screening of large mutant libraries. To validate the developed protocol, we calculated thermodynamic stability differences for 109 mutations in the microbial Ribonuclease Barnase. The obtained quantitative agreement with experimental data illustrates the potential of the approach in protein engineering and design.

Sequence periodicity and secondary structure propensity in model proteins

We explore the question of whether local effects (originating from the amino acids intrinsic secondary structure propensities) or nonlocal effects (reflecting the sequence of amino acids as a whole) play a larger role in determining the fold of globular proteins. Earlier circular dichroism studies have shown that the pattern of polar, non polar amino acids (nonlocal effect) dominates over the amino acid intrinsic propensity (local effect) in determining the secondary structure of oligomeric peptides. In this article, we present a coarse grained computational model that allows us to quantitatively estimate the role of local and nonlocal factors in determining both the secondary and tertiary structure of small, globular proteins. The amino acid intrinsic secondary structure propensity is modeled by a dihedral potential term. This dihedral potential is parametrized to match with experimental measurements of secondary structure propensity. Similarly, the magnitude of the attraction between hydrophobic residues is parametrized to match the experimental transfer free energies of hydrophobic amino acids. Under these parametrization conditions, we systematically explore the degree of frustration a given polar, non polar pattern can tolerate when the secondary structure intrinsic propensities are in opposition to it. When the parameters are in the biophysically relevant range, we observe that the fold of small, globular proteins is determined by the pattern of polar, non polar amino acids regardless of their instrinsic secondary structure propensities. Our simulations shed new light on previous observations that tertiary interactions are more influential in determining protein structure than secondary structure propensity. The fact that this can be inferred using a simple polymer model that lacks most of the biochemical details points to the fundamental importance of binary patterning in governing folding.

Requirement of Ala residues at g position in heptad sequence of alpha-helix-forming peptide for formation of fibrous structure

One feature of the alpha3-peptide, which has the amino acid sequence of (Leu-Glu-Thr-Leu-Ala-Lys-Ala)(3), that distinguishes it from many other alpha-helix-forming peptides is its ability to form fibrous assemblies that can be observed by transmission electron microscopy. In this study, the effects of Ala-->Gln substitution at the e (5th) or g (7th) position in the above heptad sequence of the alpha3-peptide on the formation of alpha-helix and fibrous assemblies were investigated by circular dichroism spectral measurement and atomic force microscopy. The 5Qalpha3-peptide obtained by Ala-->Gln substitution at the e position of the alpha3-peptide was found to form very short fibrils with long-elliptical shape, whereas the 7Qalpha3-peptide with Gln residues at the g position lost its ability to form such assemblies, in spite of alpha-helix formation in both peptides; the stabilities of both peptides decreased. These results indicate that Ala residues at the g position in the heptad sequence of the alpha3-peptide are key residues for the formation of fibrous assemblies, which may be due to hydrophobic interactions between alpha-helical bundle surfaces.

Transformation of an alpha-helix peptide into a beta-hairpin induced by addition of a fragment results in creation of a coexisting state

Intrinsic rules of determining the tertiary structure of a protein have been unknown partly because physicochemical factors that contribute to stabilization of a protein structure cannot be represented as a linear combination of local interactions. To clarify the rules on the nonlinear term caused by nonlocal interaction in a protein, we tried to transform a peptide that has a fully helical structure (Target Peptide or TP) into a peptide that has a beta-hairpin structure (Designed Peptide or DP) by adding seven residues to the C terminus of TP. According to analyses of nuclear magnetic resonance measurements, while the beta-hairpin structure is stabilized in some DPs, it is evident that the helical structure observed in TP is also persistent and even extended throughout the length of the molecule. As a result, we have produced a peptide molecule that contains both the alpha-helix and beta-hairpin conformation at an almost equally populated level. The helical structures contained in these DPs were more stable than the helix in TP, suggesting that stabilizing one conformation does not result in destabilizing the other conformation. These DPs can thus be regarded as an isolated peptide version of the chameleon sequence, which has the capability of changing the secondary structure depending on the context of the surrounding environment in a protein structure. The fact that the transformation of one secondary structure caused stabilization of both the original and the induced structure would shed light on the mechanism of protein folding.

Conformational entropy of alanine versus glycine in protein denatured states

The presence of a solvent-exposed alanine residue stabilizes a helix by 0.4-2 kcal.mol(-1) relative to glycine. Various factors have been suggested to account for the differences in helical propensity, from the higher conformational freedom of glycine sequences in the unfolded state to hydrophobic and van der Waals' stabilization of the alanine side chain in the helical state. We have performed all-atom molecular dynamics simulations with explicit solvent and exhaustive sampling of model peptides to address the backbone conformational entropy difference between Ala and Gly in the denatured state. The mutation of Ala to Gly leads to an increase in conformational entropy equivalent to approximately 0.4 kcal.mol(-1) in a fully flexible denatured, that is, unfolded, state. But, this energy is closely counterbalanced by the (measured) difference in free energy of transfer of the glycine and alanine side chains from the vapor phase to water so that the unfolded alanine- and glycine-containing peptides are approximately isoenergetic. The helix-stabilizing propensity of Ala relative to Gly thus mainly results from more favorable interactions of Ala in the folded helical structure. The small difference in energetics in the denatured states means that the Phi-values derived from Ala --> Gly scanning of helices are a very good measure of the extent of formation of structure in proteins with little residual structure in the denatured state.

Alpha-helix stabilization by alanine relative to glycine: roles of polar and apolar solvent exposures and of backbone entropy

The energetics of alpha-helix formation are fairly well understood and the helix content of a given amino acid sequence can be calculated with reasonable accuracy from helix-coil transition theories that assign to the different residues specific effects on helix stability. In internal helical positions, alanine is regarded as the most stabilizing residue, whereas glycine, after proline, is the more destabilizing. The difference in stabilization afforded by alanine and glycine has been explained by invoking various physical reasons, including the hydrophobic effect and the entropy of folding. Herein, the contribution of these two effects and that of hydrophilic area burial is evaluated by analyzing Ala and Gly mutants implemented in three helices of apoflavodoxin. These data, combined with available data for similar mutations in other proteins (22 Ala/Gly mutations in alpha-helices have been considered), allow estimation of the difference in backbone entropy between alanine and glycine and evaluation of its contribution and that of apolar and polar area burial to the helical stabilization typically associated to Gly-->Ala substitutions. Alanine consistently stabilizes the helical conformation relative to glycine because it buries more apolar area upon folding and because its backbone entropy is lower. However, the relative contribution of polar area burial (which is shown to be destabilizing) and of backbone entropy critically depends on the approximation used to model the structure of the denatured state. In this respect, the excised-peptide model of the unfolded state, proposed by Creamer and coworkers (1995), predicts a major contribution of polar area burial, which is in good agreement with recent quantitations of the relative enthalpic contribution of Ala and Gly residues to alpha-helix formation.

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.

A model for the coupling of alpha-helix and tertiary contact formation

Peptides corresponding to excised alpha-helical segments of natural proteins can spontaneously form helices in solution. However, peptide helices are usually substantially less stable in solution than in the structural context of a folded protein, because of the additional interactions possible between helices in a protein. Such interactions can be thought of as coupling helix formation and tertiary contact formation. The relative energetic contributions of the two processes to the total energy of the folded state of a protein is a matter of current debate. To investigate this balance, an extended helix-coil model (XHC) that incorporates both effects has been constructed. The model treats helix formation with the Lifson-Roig formalism, which describes helix initiation and propagation through cooperative local interactions. The model postulates an additional parameter representing participation of a site in a tertiary contact. In the model, greater helix stability can be achieved through combinations of these short-range and long-range interactions. For instance, stronger tertiary contacts can compensate for helices with little intrinsic stability. By varying the strength of the nonlocal interactions, the model can exhibit behavior consistent with a variety of qualitative models describing the relative importance of secondary and tertiary structure. Moreover, the model is explicit in that it can be used to fit experimental data to individual peptide sequences, providing a means to quantify the two contributions on a common energetic basis.

Structure-based Prediction of bZIP Partnering Specificity

Predicting protein interaction specificity from sequence is an important goal in computational biology. We present a model for predicting the interaction preferences of coiled-coil peptides derived from bZIP transcription factors that performs very well when tested against experimental protein microarray data. We used only sequence information to build atomic-resolution structures for 1711 dimeric complexes, and evaluated these with a variety of functions based on physics, learned empirical weights or experimental coupling energies. A purely physical model, similar to those used for protein design studies, gave reasonable performance. The results were improved significantly when helix propensities were used in place of a structurally explicit model to represent the unfolded reference state. Further improvement resulted upon accounting for residue-residue interactions in competing states in a generic way. Purely physical structure-based methods had difficulty capturing core interactions accurately, especially those involving polar residues such as asparagine. When these terms were replaced with weights from a machine-learning approach, the resulting model was able to correctly order the stabilities of over 6000 pairs of complexes with greater than 90% accuracy. The final model is physically interpretable, and suggests specific pairs of residues that are important for bZIP interaction specificity. Our results illustrate the power and potential of structural modeling as a method for predicting protein interactions and highlight obstacles that must be overcome to reach quantitative accuracy using a de novo approach. Our method shows unprecedented performance in predicting protein-protein interaction specificity accurately using structural modeling and suggests that predicting coiled-coil interactions generally may be within reach.

Stabilization of the fibrous structure of an α-helix-forming peptide by sequence reversal

The alpha3-peptide, which comprises three repeats of the sequence Leu-Glu-Thr-Leu-Ala-Lys-Ala and forms an amphipathic alpha-helix, is unique among various alpha-helix-forming peptides in that it assembles into fibrous structures that can be observed by transmission electron microscopy. As part of our investigation of the structure-stability relationships of the alpha3-peptide, we synthesized the r3-peptide, whose amino acid sequence is the reverse of that of the alpha3-peptide, and we investigated the effects of sequence reversal on alpha-helix stability and the formation of fibrous structures. Unexpectedly, the r3-peptide formed a more-stable alpha-helix and longer fibers than did the alpha3-peptide. The stability of the r3-peptide helix decreased when the ionic strength of the buffer was increased and when the pH of the buffer was adjusted to 2 or 12. These results suggest that the r3-peptide underwent a "magnet-like" oligomerization and that an increase in the charge-distribution inequality may be the driving force for the formation of fibrous structures.

Thermodynamic and Structural Equivalence of Two HLA-B27 Subtypes Complexed with a Self-peptide

The F pocket of major histocompatibility complex (in humans HLA) class I molecules accommodates the C terminus of the bound peptide. Residues forming this pocket exhibit considerable polymorphism, and a single difference (Asp116 in HLA-B*2705 and His116 in HLA-B*2709 heavy chains) confers differential association of these two HLA-B27 subtypes to the autoimmune disease ankylosing spondylitis. As peptide presentation by HLA molecules is of central importance for immune responses, we performed thermodynamic (circular dichroism, differential scanning calorimetry, fluorescence polarization) and X-ray crystallographic analyses of both HLA-B27 subtypes complexed with the epidermal growth factor response factor 1-derived self-peptide TIS (RRLPIFSRL) to understand the impact of the Asp116His exchange on peptide display. This peptide is known to be presented in vivo by both subtypes, and as expected for a self-peptide, TIS-reactive cytotoxic T lymphocytes are absent in the respective individuals. The thermodynamic analyses reveal that both HLA-B27:TIS complexes exhibit comparable, relatively high thermostability (Tm approximately 60 degrees C) and undergo multi-step unfolding reactions, with dissociation of the peptide in the first step. As shown by X-ray crystallography, only subtle structural differences between the subtypes were observed regarding the architecture of their F pockets, including the presence of distinct networks of water molecules. However, no consistent structural differences were found between the peptide presentation modes. In contrast to other peptides displayed by the two HLA-subtypes which show either structural or dynamical differences in their peptide presentation modes, the TIS-complexed HLA-B*2705 and HLA-B*2709 subtypes are an example for thermodynamic and structural equivalence, in agreement with functional data.

Enthalpy of helix-coil transition: missing link in rationalizing the thermodynamics of helix-forming propensities of the amino acid residues

It is known that different amino acid residues have effects on the thermodynamic stability of an alpha-helix. The underlying mechanism for the thermodynamic helical propensity is not well understood. The major accepted hypothesis is the difference in the side-chain configurational entropy loss upon helix formation. However, the changes in the side-chain configurational entropy explain only part of the thermodynamic helical propensity, thus implying that there must be a difference in the enthalpy of helix-coil transition for different residues. This work provides an experimental test to this hypothesis. Direct calorimetric measurements of folding of a model host peptide in which the helix formation is induced by metal binding is applied to a wide range of residue types, both naturally occurring and nonnatural, at the guest site. Based on the calorimetric results for 12 peptides, it was found that indeed there is a difference in the enthalpy of helix-coil transition for different amino acid residues, and simple empirical rules that define these differences are presented. The obtained difference in the enthalpies of helix-coil transition complement the differences in configurational entropies and provide the complete thermodynamic characterization of the helix-forming tendencies.

Thermodynamics Of α‐Helix Formation

The alpha-helix was the first proposed and experimentally confirmed secondary structure. The elegant simplicity of the alpha-helical structure, stabilized by hydrogen bonding between the backbone carbonyl oxygen and the peptide amide four residues away, has captivated the scientific community. In proteins, alpha-helices are also stabilized by the so-called capping interactions that occur at both the C- and the N-termini of the helix. This chapter provides a brief historical overview of the thermodynamic studies of the energetics of helix formation, and reviews recent progress in our understanding of the thermodynamics of helix formation.

Structural and functional characteristics of the Val44Met insulin-like growth factor I missense mutation: correlation with effects on growth and development

We have previously described the phenotype resulting from a missense mutation in the IGF-I gene, which leads to expression of IGF-I with a methionine instead of a valine at position 44 (Val44Met IGF-I). This mutation caused severe growth and mental retardation as well as deafness evident at birth and growth retardation in childhood, but is relatively well tolerated in adulthood. We have conducted a biochemical and structural analysis of Val44Met IGF-I to provide a molecular basis for the phenotype observed. Val44Met IGF-I exhibits a 90-fold decrease in type 1 IGF receptor (IGF-1R) binding compared with wild-type human IGF-I and only poorly stimulates autophosphorylation of the IGF-1R. The ability of Val44Met IGF-I to signal via the extracellular signal-regulated kinase 1/2 and Akt/protein kinase B pathways and to stimulate DNA synthesis is correspondingly poorer. Binding or activation of both insulin receptor isoforms is not detectable even at micromolar concentrations. However, Val44Met IGF-I binds IGF-binding protein-2 (IGFBP-2), IGFBP-3, and IGFBP-6 with equal affinity to IGF-I, suggesting the maintenance of overall structure, particularly in the IGFBP binding domain. Structural analysis by nuclear magnetic resonance confirms retention of near-native structure with only local side-chain disruptions despite the significant loss of function. To our knowledge, our results provide the first structural study of a naturally occurring mutant human IGF-I associated with growth and developmental abnormalities and identifies Val44 as an essential residue involved in the IGF-IGF-1R interaction.