Role of non-glycine residues in left-handed helical conformation for the conformational stability of human lysozyme.

To understand the role of non-Gly residues in the left-handed helical conformation for the conformational stability of a protein, the non-Gly to Gly and Ala mutations at six left-handed residues (R21, Y38, R50, Q58, H78, and N118) of the human lysozyme were examined. The thermodynamic parameters for denaturation were determined using a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. If a left-handed non-Gly had an unfavorable steric interaction between the side-chain Cbeta and backbone, the Gly mutation would be expected to stabilize more than the Ala mutation at the same position. For the mutant human lysozymes, however, there were few differences in the denaturation Gibbs energy (DeltaG) between the Gly and Ala mutants, except for the substitution at position 58. Analysis of the changes in stability (DeltaDeltaG) based on the structures of the wild-type and mutant proteins showed that the experimental DeltaDeltaG value of Q58G was approximately 7 kJ/mol higher than the estimated value without consideration of any local steric interaction. These results indicate that only Q58G increased the stability by elimination of local constraints. The residue 58 is located at the most rigid position in the left-handed non-Gly residues and is involved in its enzymatic function. It can be concluded that the left-handed non-Gly residues do not always have unfavorable strain energies as compared with Gly at the same position.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

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.     

Cooperative propagation of local stability changes from low-stability and high-stability regions in a SH3 domain

Site-directed mutagenesis has been used to produce local stability changes at two regions of the binding site surface of the alpha-spectrin SH3 domain (Spc-SH3) differing in their intrinsic stability. Mutations were made at residue 56, located at the solvent-exposed side of the short 3(10) helix, and at residue 21 in the tip of the flexible RT-loop. NMR chemical-shift analysis and X-ray crystallography indicated negligible changes produced by the mutations in the native structure limited to subtle rearrangements near the mutated residue and at flexible loops. Additionally, mutations do not alter importantly the SH3 binding site structure, although produce significant changes in its affinity for a proline-rich decapeptide. The changes in global stability measured by differential scanning calorimetry are consistent the local energy changes predicted by theoretical models, with the most significant effects observed for the Ala-Gly mutations. Propagation of the local stability changes throughout the domain structure has been studied at a per-residue level of resolution by NMR-detected amide hydrogen-deuterium exchange (HX). Stability propagation is remarkably efficient in this small domain, apparently due to its intrinsically low stability. Nevertheless, the HX-core of the domain is not fully cooperative, indicating the existence of co-operative subunits within the core, which is markedly polarized. An equilibrium phi-analysis of the changes in the apparent Gibbs energies of HX per residue produced by the mutations has allowed us to characterize structurally the conformational states leading to HX. Some of these states resemble notably the folding transition state of the Spc-SH3 domain, suggesting a great potential of this approach to explore the folding energy landscape of proteins. An energy perturbation propagates more effectively from a flexible region to the core than in the opposite direction, because the former affects a broader region of the energy landscape than the latter. This might be of importance in understanding the special thermodynamic signature of the SH3-peptide interaction and the relevance of the dual character of SH3 binding sites.     

Increasing and decreasing protein stability: effects of revertant substitutions on the thermal denaturation of phage lambda repressor

The thermal denaturations of five revertant lambda repressors containing single amino acid substitutions in their N-terminal domains have been studied by differential scanning calorimetry. Two substitutions slightly decrease stability, and the remaining three render the protein more stable than wild type. The Gly48----Asn and Gly48----Ser proteins are 4 degrees C more stable than wild type. These two substitutions replace an alpha helical residue, and in each case a poor helix forming residue, glycine, is replaced by a residue with a higher helical propensity. We also present data showing that one revertant, Tyr22----Phe, has reduced operator DNA binding affinity despite its enhanced stability.     

Crystal structure of a consensus-designed ankyrin repeat protein: implications for stability

Consensus-designed ankyrin repeat (AR) proteins are thermodynamically very stable. The structural analysis of the designed AR protein E3_5 revealed that this stability is due to a regular fold with highly conserved structural motifs and H-bonding networks. However, the designed AR protein E3_19 exhibits a significantly lower stability than E3_5 (9.6 vs. 14.8 kcal/mol), despite 88% sequence identity. To investigate the structural correlations of this stability difference between E3_5 and E3_19, we determined the crystal structure of E3_19 at 1.9 A resolution. E3_19 as well has a regular AR domain fold with the characteristic H-bonding patterns. All structural features of the E3_5 and E3_19 molecules appear to be virtually identical (RMSD(Calpha) approximately 0.7 A). However, clear differences are observed in the surface charge distribution of the two AR proteins. E3_19 features clusters of charged residues and more exposed hydrophobic residues than E3_5. The atomic coordinates of E3_19 have been deposited in the Protein Data Bank. PDB ID: 2BKG.     

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 α‐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 α‐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 α‐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.     

Conformations of helical Aib peptides containing a pair of l‐amino acid and d‐amino acid

A pair of l ‐leucine (l ‐Leu) and d ‐leucine (d ‐Leu) was incorporated into α‐aminoisobutyric acid (Aib) peptide segments. The dominant conformations of four hexapeptides, Boc‐l ‐Leu‐Aib‐Aib‐Aib‐Aib‐l ‐Leu‐OMe (1a), Boc‐d ‐Leu‐Aib‐Aib‐Aib‐Aib‐l ‐Leu‐OMe (1b), Boc‐Aib‐Aib‐l ‐Leu‐l ‐Leu‐Aib‐Aib‐OMe (2a), and Boc‐Aib‐Aib‐d ‐Leu‐l ‐Leu‐Aib‐Aib‐OMe (2b), were investigated by IR, ¹H NMR, CD spectra, and X‐ray crystallographic analysis. All peptides 1a,b and 2a,b formed 3₁₀‐helical structures in solution. X‐ray crystallographic analysis revealed that right‐handed (P) 3₁₀‐helices were present in 1a and 1b and a mixture of right‐handed (P) and left‐handed (M) 3₁₀‐helices was present in 2b in their crystalline states. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

Role of solvation barriers in protein kinetic stability

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).     

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (～ms) for substrate binding and conformational changes to the fast time (～ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X‐ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild‐type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.     

Natural selection for kinetic stability is a likely origin of correlations between mutational effects on protein energetics and frequencies of amino acid occurrences in sequence alignments

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.     

Structure, conformational stability, and enzymatic properties of acylphosphatase from the hyperthermophile Sulfolobus solfataricus

The structure of AcP from the hyperthermophilic archaeon Sulfolobus solfataricus has been determined by (1)H-NMR spectroscopy and X-ray crystallography. Solution and crystal structures (1.27 A resolution, R-factor 13.7%) were obtained on the full-length protein and on an N-truncated form lacking the first 12 residues, respectively. The overall Sso AcP fold, starting at residue 13, displays the same betaalphabetabetaalphabeta topology previously described for other members of the AcP family from mesophilic sources. The unstructured N-terminal tail may be crucial for the unusual aggregation mechanism of Sso AcP previously reported. Sso AcP catalytic activity is reduced at room temperature but rises at its working temperature to values comparable to those displayed by its mesophilic counterparts at 25-37 degrees C. Such a reduced activity can result from protein rigidity and from the active site stiffening due the presence of a salt bridge between the C-terminal carboxylate and the active site arginine. Sso AcP is characterized by a melting temperature, Tm, of 100.8 degrees C and an unfolding free energy, DeltaG(U-F)H2O, at 28 degrees C and 81 degrees C of 48.7 and 20.6 kJ mol(-1), respectively. The kinetic and structural data indicate that mesophilic and hyperthermophilic AcP's display similar enzymatic activities and conformational stabilities at their working conditions. Structural analysis of the factor responsible for Sso AcP thermostability with respect to mesophilic AcP's revealed the importance of a ion pair network stabilizing particularly the beta-sheet and the loop connecting the fourth and fifth strands, together with increased density packing, loop shortening and a higher alpha-helical propensity.     

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

Pfam family DUF1023 consists entirely of uncharacterized proteins generated by sequencing the genomes of Actinobacteria (Bateman A., et al., Nucleic Acids Res. 2004;32 Database issue:D138-141.) Utilizing sequence similarity detection methods, we infer homology between DUF1023 and alpha/beta hydrolases. DUF1023 proteins conserve the core secondary structures in alpha/beta hydrolase fold, and share similar catalytic machinery as that of alpha/beta hydrolases. We predict DUF1023 spatial structure and deduce that they function as hydrolases utilizing catalytic Ser-His-Asp triad with the serine as a nucleophile.     

Functional clues from the crystal structure of an orphan periplasmic ligand‐binding protein from Treponema pallidum

The spirochete Treponema pallidum is the causative agent of syphilis, a sexually transmitted infection of major global importance. Other closely related subspecies of Treponema also are the etiological agents of the endemic treponematoses, such as yaws, pinta, and bejel. The inability of T. pallidum and its close relatives to be cultured in vitro has prompted efforts to characterize T. pallidum's proteins structurally and biophysically, particularly those potentially relevant to treponemal membrane biology, with the goal of possibly revealing the functions of those proteins. This report describes the structure of the treponemal protein Tp0737; this polypeptide has a fold characteristic of a class of periplasmic ligand‐binding proteins associated with ABC‐type transporters. Although no ligand for the protein was observed in electron‐density maps, and thus the nature of the native ligand remains obscure, the structural data described herein provide a foundation for further efforts to elucidate the ligand and thus the function of this protein in T. pallidum.     

Crystal structure of a hypothetical protein,TM841 of Thermotoga maritima,reveals its function as a fatty acid-binding protein

We determined the three-dimensional (3D) crystal structure of protein TM841, a protein product from a hypothetical open-reading frame in the genome of the hyperthermophile bacterium Thermotoga maritima, to 2.0 A resolution. The protein belongs to a large protein family, DegV or COG1307 of unknown function. The 35 kDa protein consists of two separate domains, with low-level structural resemblance to domains from other proteins with known 3D structures. These structural homologies, however, provided no clues for the function of TM841. But the electron density maps revealed clear density for a bound fatty-acid molecule in a pocket between the two protein domains. The structure indicates that TM841 has the molecular function of fatty-acid binding and may play a role in the cellular functions of fatty acid transport or metabolism.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Comparative structural stability of subunits of the potato virus X coat protein in solution and in virus particles

Several optical methods and differential scanning calorimetry were used to study the structure and stability of free coat protein (CP) molecules and CP molecules in the virion of the potato virus X (PVX), a filamentous plant virus. All criteria suggest that PVX CP (hereinafter, CP) subunits in solution at room temperature display a certain preserved tertiary structure; however, this structure is very unstable and already denatures at 35°C. Very low concentrations of sodium dodecylsulfate or cetyltrimethylammonium bromide also disrupt the CP tertiary structure, three-five molecules of these detergents per one protein molecule being sufficient. However, the secondary structure of CP molecules does not change under the same conditions. Once included into the virion, CP subunits become considerably more stable towards increased temperature and detergents. This combination of a highly labile tertiary structure and a fairly stable secondary structure of free CP can be a structural basis for the recently discovered ability of PVX CP to assume two distinct functional states within the virion.     

The effect of specific proline residues on the kinetic stability of the triosephosphate isomerases of two trypanosomes

The effect of specific residues on the kinetic stability of two closely related triosephosphate isomerases (from Trypanosoma cruzi, TcTIM and Trypanosoma brucei, TbTIM) has been studied. Based on a comparison of their β‐turn occurrence, we engineered two chimerical enzymes where their super secondary β‐loop‐α motifs 2 ((βα)₂) were swapped. Differential scanning calorimetry (DSC) experiments showed that the (βα)₂ motif of TcTIM inserted into TbTIM (2Tc) increases the kinetic stability. On the other hand, the presence of the (βα)₂ motif of TbTIM inserted into TcTIM (2Tb) gave a chimerical protein difficult to purify in soluble form and with a significantly reduced kinetic stability. The comparison of the contact maps of the (βα)₂ of TbTIM and TcTIM showed differences in the contact pattern of residues 43 and 49. In TcTIM these residues are prolines, located at the N‐terminal of loop‐2 and the C‐terminal of α‐helix‐2. Twelve mutants were engineered involving residues 43 and 49 to study the effect over the unfolding activation energy barrier (E_A). A systematic analysis of DSC data showed a large decrease on the E_A of TcTIM (ΔE_A ranging from 468 to 678 kJ/mol) when the single and double proline mutations are present. The relevance of Pro43 to the kinetic stability is also revealed by mutation S43P, which increased the free energy of the transition state of TbTIM by 17.7 kJ/mol. Overall, the results indicate that protein kinetic stability can be severely affected by punctual mutations, disturbing the complex network of interactions that, in concerted action, determine protein stability. Proteins 2017; 85:571–579. © 2016 Wiley Periodicals, Inc.