Contribution of the 30/36 hydrophobic contact at the C-terminus of the alpha-helix to the stability of the ubiquitin molecule

Almost 50 Years of Monomeric Extracellular Ubiquitin (eUb)

Monomeric ubiquitin (Ub) is a 76-amino-acid highly conserved protein found in eukaryotes. The biological activity of Ub first described in the 1970s was extracellular, but it quickly gained relevance due to its intracellular role, i.e., post-translational modification of intracellular proteins (ubiquitination) that regulate numerous eukaryotic cellular processes. In the following years, the extracellular role of Ub was relegated to the background, until a correlation between higher survival rate and increased serum Ub concentrations in patients with sepsis and burns was observed. Although the mechanism of action (MoA) of extracellular ubiquitin (eUb) is not yet well understood, further studies have shown that it may ameliorate the inflammatory response in tissue injury and multiple sclerosis diseases. These observations, compounded with the high stability and low immunogenicity of eUb due to its high conservation in eukaryotes, have made this small protein a relevant candidate for biotherapeutic development. Here, we review the in vitro and in vivo effects of eUb on immunologic, cardiovascular, and nervous systems, and discuss the potential MoAs of eUb as an anti-inflammatory, antimicrobial, and cardio- and brain-protective agent.

Mapping side chain interactions at protein helix termini

Background

Interactions that involve one or more amino acid side chains near the ends of protein helices stabilize helix termini and shape the geometry of the adjacent loops, making a substantial contribution to overall protein structure. Previous work has identified key helix-terminal motifs, such as Asx/ST N-caps, the capping box, and hydrophobic and electrostatic interactions, but important questions remain, including: 1) What loop backbone geometries are favoured by each motif? 2) To what extent are multi-amino acid motifs likely to represent genuine cooperative interactions? 3) Can new motifs be identified in a large, recent dataset using the latest bioinformatics tools?

Results

Three analytical tools are applied here to answer these questions. First, helix-terminal structures are partitioned by loop backbone geometry using a new 3D clustering algorithm. Next, Cascade Detection, a motif detection algorithm recently published by the author, is applied to each cluster to determine which sequence motifs are overrepresented in each geometry. Finally, the results for each motif are presented in a CapMap, a 3D conformational heatmap that displays the distribution of the motif’s overrepresentation across loop geometries, enabling the rapid isolation and characterization of the associated side chain interaction. This work identifies a library of geometry-specific side chain interactions that provides a new, detailed picture of loop structure near the helix terminus. Highlights include determinations of the favoured loop geometries for the Asx/ST N-cap motifs, capping boxes, “big” boxes, and other hydrophobic, electrostatic, H-bond, and pi stacking interactions, many of which have not been described before.

Conclusions

This work demonstrates that the combination of structural clustering and motif detection in the sequence space can efficiently identify side chain motifs and map them to the loop geometries which they support. Protein designers should find this study useful, because it identifies side chain interactions which are good candidates for inclusion in synthetic helix-terminal loops with specific desired geometries, since they are used in nature to support these geometries. The techniques described here can also be applied to map side chain interactions associated with other structural components of proteins such as beta and gamma turns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0671-4) contains supplementary material, which is available to authorized users.

Thermal coefficients of the methyl groups within ubiquitin

Physiological processes such as protein folding and molecular recognition are intricately linked to their dynamic signature, which is reflected in their thermal coefficient. In addition, the local conformational entropy is directly related to the degrees of freedom, which each residue possesses within its conformational space. Therefore, the temperature dependence of the local conformational entropy may provide insight into understanding how local dynamics may affect the stability of proteins. Here, we analyze the temperature dependence of internal methyl group dynamics derived from the cross-correlated relaxation between dipolar couplings of two CH bonds within ubiquitin. Spanning a temperature range from 275 to 308 K, internal methyl group dynamics tend to increase with increasing temperature, which translates to a general increase in local conformational entropy. With this data measured over multiple temperatures, the thermal coefficient of the methyl group order parameter, the characteristic thermal coefficient, and the local heat capacity were obtained. By analyzing the distribution of methyl group thermal coefficients within ubiquitin, we found that the N-terminal region has relatively high thermostability. These results indicate that methyl groups contribute quite appreciably to the total heat capacity of ubiquitin through the regulation of local conformational entropy.

CAPS-DB: a structural classification of helix-capping motifs

The regions of the polypeptide chain immediately preceding or following an α-helix are known as Nt- and Ct cappings, respectively. Cappings play a central role stabilizing α-helices due to lack of intrahelical hydrogen bonds in the first and last turn. Sequence patterns of amino acid type preferences have been derived for cappings but the structural motifs associated to them are still unclassified. CAPS-DB is a database of clusters of structural patterns of different capping types. The clustering algorithm is based in the geometry and the (ϕ–ψ)-space conformation of these regions. CAPS-DB is a relational database that allows the user to search, browse, inspect and retrieve structural data associated to cappings. The contents of CAPS-DB might be of interest to a wide range of scientist covering different areas such as protein design and engineering, structural biology and bioinformatics. The database is accessible at: http://www.bioinsilico.org/CAPSDB.

Synthesis of a trigalacturonic acid analogue mimicking the expected transition state in the glycosidases

A trigalacturonic acid analogue carrying a cyclohexene framework in place of the central pyranose ring was synthesized as a molecular probe for the mechanistic investigation of endo-polygalacturonase 1 (endo-PG 1). Preliminary enzymatic studies revealed that this analogue inhibited endo-PG 1 activity by about 30% at 0.3mM concentration.

A computational approach for the rational design of stable proteins and enzymes: optimization of surface charge-charge interactions

The design of stable proteins and enzymes is not only of particular biotechnological importance, but also addresses some important fundamental questions. While there are a number of different options available for designing or engineering stable proteins, the field of computational design provides fast and universal methods for stabilizing proteins of interest. One of the successful computational design strategies focuses on stabilizing proteins through the optimization of charge-charge interactions on the protein surface. By optimizing surface interactions, it is possible to alleviate some of the challenges that accompany efforts to redesign the protein core. The rational design of surface charge-charge interactions also allows one to optimize only the interactions that are distant from binding sites or active sites, making it possible to increase stability without adversely affecting activity. The optimization of surface charge-charge interactions is discussed in detail along with the experimental evidence to demonstrate that this is a robust and universal approach to designing proteins with enhanced stability.

Molecular dynamics simulations of the native and partially folded states of ubiquitin: influence of methanol cosolvent, pH, and temperature on the protein structure and dynamics

A series of explicit-solvent molecular dynamics simulations of the protein ubiquitin are reported, which investigate the effect of environmental factors (presence of methanol cosolvent in the aqueous solution, neutral or low pH value, room or elevated temperature) on the structure, stability, and dynamics of the protein. The simulations are initiated either from the native structure of the protein or from a model of a partially folded state (A-state) that is known to exist at low pH in methanol-water mixtures. The main results of the simulations are: (1) The ubiquitin native structure is remarkably stable at neutral pH in water; (2) the addition of the methanol cosolvent enhances the stability of the secondary structure but weakens tertiary interactions within the protein; (3) this influence of methanol on the protein structure is enhanced at low pH, while the effect of lowering the pH in pure water is limited; and (4) the A-state of ubiquitin can be described as a set of relatively rigid secondary structure elements (a native-like beta-sheet and native-like alpha-helix plus two nonnative alpha-helices) connected by flexible linkers.

Computational spectroscopy of ubiquitin: Comparison between theory and experiments

Using the constrained molecular dynamics simulation method in combination with quantum chemistry calculation, Hessian matrix reconstruction, and fragmentation approximation methods, the authors have established computational schemes for numerical simulations of amide I IR absorption, vibrational circular dichroism (VCD), and two-dimensional (2D) IR photon echo spectra of the protein ubiquitin in water. Vibrational characteristic features of these spectra in the amide I vibration region are discussed. From the semiempirical quantum chemistry calculation results on an isolated ubiquitin, amide I local mode frequencies and vibrational coupling constants were fully determined. It turns out that the amide I local mode frequencies of ubiquitin in both gas phase and aqueous solution are highly heterogeneous and site dependent. To directly test the quantitative validity of thus obtained spectroscopic properties, they compared the experimentally measured amide I IR, 2D IR, and electronic circular dichroism spectra with experiments, and found good agreements between theory and experiments. However, the simulated VCD spectrum is just qualitatively similar to the experimentally measured one. This indicates that, due to delicate cancellations between the positive and negative VCD contributions, the prediction of protein VCD spectrum is critically relied on quantitative accuracy of the theoretical model for predicting amide I local mode frequencies. On the basis of the present comparative investigations, they found that the site dependency of amide I local mode frequency, i.e., diagonal heterogeneity of the vibrational Hamiltonian matrix in the amide I local mode basis, is important. It is believed that the present computational methods for simulating various vibrational and electronic spectra of proteins will be of use in further refining classical force fields and in addressing the structure-spectra relationships of proteins in solution.

Amino acid substitutions affecting protein dynamics in eglin C do not affect heat capacity change upon unfolding

The heat capacity change upon unfolding (deltaC(p)) is a thermodynamic parameter that defines the temperature dependence of the thermodynamic stability of proteins; however, physical basis of the heat capacity change is not completely understood. Although empirical surface area-based calculations can predict heat capacity changes reasonably well, accumulating evidence suggests that changes in hydration of those surfaces is not the only parameter contributing to the observed heat capacity changes upon unfolding. Because packing density in the protein interior is similar to that observed in organic crystals, we hypothesized that changes in protein dynamics resulting in increased rigidity of the protein structure might contribute to the observed heat capacity change upon unfolding. Using differential scanning calorimetry we characterized the thermodynamic behavior of a serine protease inhibitor eglin C and two eglin C variants with altered native state dynamics, as determined by NMR. We found no evidence of changes in deltaC(p) in either of the variants, suggesting that changes in rigidity do not contribute to the heat capacity change upon unfolding in this model system.

Ubiquitin: a small protein folding paradigm

For the past twenty years, the small, 76-residue protein ubiquitin has been used as a model system to study protein structure, stability, folding and dynamics. In this time, ubiquitin has become a paradigm for both the experimental and computational folding communities. The folding energy landscape is now uniquely characterised with a plethora of information available on not only the native and denatured states, but partially structured states, alternatively folded states and locally unfolded states, in addition to the transition state ensemble. This Perspective focuses on the experimental characterisation of ubiquitin using a comprehensive range of biophysical techniques.

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.

Both helical propensity and side-chain hydrophobicity at a partially exposed site in α-helix contribute to the thermodynamic stability of ubiquitin

Improving helical propensity of residues was proposed as one of the approaches to increase protein stability. Here the contribution of the helix propensity and hydrophobicity of residues at partially buried positions of alpha-helix to the stability of a model protein-ubiquitin- is explored. Thermodynamic stabilities of 13 ubiquitin variants with substitutions at a partially buried helical residue were measured by differential scanning calorimetry. It was found that the dynamic range of stabilities for different amino acid residues at this partially buried position is 3 times larger than that expected based on helical propensity alone. Correlation analysis shows that both helical propensity and hydrophobicity are important in defining the relative stabilities of the studied ubiquitin variants. These results provide experimental evidence that partially buried positions are potentially useful sites for engineering proteins with enhanced thermostability.

Noncharged amino acid residues at the solvent-exposed positions in the middle and at the C terminus of the alpha-helix have the same helical propensity

It was established previously that helical propensities of different amino acid residues in the middle of alpha-helix in peptides and in proteins are very similar. The statistical analysis of the protein helices from the known three-dimensional structures shows no difference in the frequency of noncharged residues in the middle and at the C terminus. Yet, experimental studies show distinctive differences for the helical propensities of noncharged residues in the middle and in the C terminus in model peptides. Is this a general effect, and is it applicable to protein helices or is it specific to the model alanine-based peptides? To answer this question, the effects of substitutions at positions 28 (middle residue) and 32 (C2 position at the C terminus) of the alpha-helix of ubiquitin on the stability of this protein are measured by using differential scanning calorimetry. The two data sets produce similar values for intrinsic helix propensity, leading to a conclusion that noncharged amino acid residues at the solvent-exposed positions in the middle and at the C terminus of the alpha-helix have the same helical propensity. This conclusion is further supported with an excellent correlation between the helix propensity scale obtained for the two positions in ubiquitin with the experimental helix propensity scale established previously and with the statistical distribution of the residues in protein helices.

Hydrophobic interactions at the Ccap position of the C-capping motif of alpha-helices

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.

Removal of surface charge-charge interactions from ubiquitin leaves the protein folded and very stable

The contribution of solvent-exposed charged residues to protein stability was evaluated using ubiquitin as a model protein. We combined site-directed mutagenesis and specific chemical modifications to first replace all Arg residues with Lys, followed by carbomylation of Lys-amino groups. Under the conditions in which all carboxylic groups are protonated (at pH 2), the chemically modified protein is folded and very stable (DeltaG = 18 kJ/mol). These results indicate that surface charge-charge interactions are not an essential fundamental force for protein folding and stability.

Hydration of the peptide backbone largely defines the thermodynamic propensity scale of residues at the C' position of the C-capping box of alpha-helices

The C' position of the C-capping box is the second residue outside of the helix. Statistical analysis of residue distribution at the C' position in the alpha-helices' C-capping box showed that different amino acid residues occur with different probabilities, with the strongest preference being for glycine. To understand the physico-chemical basis for this preference, we studied the effects that 17 amino acid substitutions at the C' position in an alpha-helix of ubiquitin have on the stability of this protein. We determined the following rank order of amino acid residues at the C' position with respect to their effect on the stability: Gly>His>Asn>Arg>Lys>Gln>Ala>Phe>Met>Ser>Asp>Glu>Trp>Thr>Pro>Ile>Val. The effect of the amino acid substitutions on the structure also was evaluated by comparing the (1)H-(15)N heteronuclear sequential quantum correlation spectra and showed no significant changes in the structures of the most stable (Gly) and the least stable (Val) variants. The obtained changes in stability highly correlate (r = 0.85) with the statistical distribution of the residues at the C' position indicating that the measured thermodynamic propensities are unbiased by secondary interactions. We also found that the measured thermodynamic propensities correlate well with the amide hydrogen exchange data on short model peptides (r = 0.85) and the calculated hydration of the peptide backbone (r = 0.88). These results combined with the changes in enthalpy and entropy of unfolding of ubiquitin variants suggest that dehydration of the peptide backbone plays a significant role in defining the thermodynamic propensity scale at the C' position of the C-capping box in alpha-helices. This propensity scale is useful for protein secondary structure predictions and protein design.

Heat capacity changes upon burial of polar and nonpolar groups in proteins

In this paper we address the question of whether the burial of polar and nonpolar groups in the protein locale is indeed accompanied by the heat capacity changes, DeltaC(p), that have an opposite sign, negative for nonpolar groups and positive for polar groups. To accomplish this, we introduced amino acid substitutions at four fully buried positions of the ubiquitin molecule (Val5, Val17, Leu67, and Gln41). We substituted Val at positions 5 and 17 and Leu at position 67 with a polar residue, Asn. As a control, Ala was introduced at the same three positions. We also replaced the buried polar Gln41 with Val and Leu, nonpolar residues that have similar size and shape as Gln. As a control, Asn was introduced at Gln41 as well. The effects of these amino acid substitutions on the stability, and in particular, on the heat capacity change upon unfolding were measured using differential scanning calorimetry. The effect of the amino acid substitutions on the structure was also evaluated by comparing the (1)H-(15)N HSQC spectra of the ubiquitin variants. It was found that the Ala substitutions did not have a considerable effect on the heat capacity change upon unfolding. However, the substitutions of aliphatic side chains (Val or Leu) with a polar residue (Asn) lead to a significant (> 30%) decrease in the heat capacity change upon unfolding. The decrease in heat capacity changes does not appear to be the result of significant structural perturbations as seen from the HSQC spectra of the variants. The substitution of a buried polar residue (Gln41) to a nonpolar residue (Leu or Val) leads to a significant (> 25%) increase in heat capacity change upon unfolding. These results indicate that indeed the heat capacity change of burial of polar and nonpolar groups has an opposite sign. However, the observed changes in DeltaC(p) are several times larger than those predicted, based on the changes in water accessible surface area upon substitution.