Local and nonlocal interactions in globular proteins and mechanisms of alcohol denaturation.

How important are helical propensities in determining the conformations of globular proteins? Using the two-dimensional lattice model and two monomer types, H (hydrophobic) and P (polar), we explore both nonlocal interactions, through an HH contact energy, as developed in earlier work, and local interactions, through a helix energy, σ. By computer enumeration, the partition functions for short chains are obtained without approximation for the full range of both types of energy. When nonlocal interactions dominate, some sequences undergo coil-globule collapse to a unique native structure. When local interactions dominate, all sequences undergo helix–coil transitions. For two different conformational properties, the closest correspondence between the lattice model and proteins in the Protein Data Bank is obtained if the model local interactions are made small compared to the HH contact interaction, suggesting that helical propensities may be only weak determinants of globular protein structures in water. For some HP sequences, varying σ/ leads to additional sharp transitions (sometimes several) and to “conformational switching” between unique conformations. This behavior resembles the transitions of globular proteins in water to helical states in alcohols. In particular, comparison with experiments shows that whereas urea as a denaturant is best modeled as weakening both local and nonlocal interactions, trifluoroethanol is best modeled as mainly weakening HH interactions and slightly enhancing local helical interactions.     

Dynamic structures of globular proteins with respect to correlative movements of residues calculated in the normal mode analysis

Dynamic structures of globular proteins are studied on the basis of correlative movements of residues around their native conformations, which are computed by means of the normal mode analysis. To describe the dynamic structures of a protein, the core regions moving with strong positive or negative correlations to other regions of the polypeptide chain are detected from the correlation maps of the movements of residues. Such core regions are different, according to the definition, from the regions defined from a geometrical point of view, such as secondary structures, domains, modules, and so on. The core regions are actually detected for four proteins, myoglobin, Bence-Jones protein, flavodoxin, and hen egg-white lysozyme, with different folding types from each other. The results show that some of them coincide with the secondary structures, domains, or modules, but others do not. Then, the dynamic structure of each protein is discussed in terms of the dynamic cores detected, as compared with the secondary structures, domains, and modules.     

Global analysis of the thermal and chemical denaturation of the N-terminal domain of the ribosomal protein L9 in H2O and D2O. Determination of the thermodynamic parameters, deltaH(o), deltaS(o), and deltaC(o)p and evaluation of solvent isotope effects.

The stability of the N-terminal domain of the ribosomal protein L9, NTL9, from Bacillus stearothermophilus has been monitored by circular dichroism at various temperatures and chemical denaturant concentrations in H2O and D2O. The basic thermodynamic parameters for the unfolding reaction, deltaH(o), deltaS(o), and deltaC(o)p, were determined by global analysis of temperature and denaturant effects on stability. The data were well fit by a model that assumes stability varies linearly with denaturant concentration and that uses the Gibbs-Helmholtz equation to model changes in stability with temperature. The results obtained from the global analysis are consistent with information obtained from individual thermal and chemical denaturations. NTL9 has a maximum stability of 3.78 +/- 0.25 kcal mol(-1) at 14 degrees C. DeltaH(o)(25 degrees C) for protein unfolding equals 9.9 +/- 0.7 kcal mol(-1) and TdeltaS(o)++(25 degrees C) equals 6.2 +/- 0.6 kcal mol(-1). DeltaC(o)p equals 0.53 +/- 0.06 kcal mol(-1) deg(-1). There is a small increase in stability when D2O is substituted for H2O. Based on the results from global analysis, NTL9 is 1.06 +/- 0.60 kcal mol(-1) more stable in D2O at 25 degrees C and Tm is increased by 5.8 +/- 3.6 degrees C in D2O. Based on the results from individual denaturation experiments, NTL9 is 0.68 +/- 0.68 kcal mol(-1) more stable in D2O at 25 degrees C and Tm is increased by 3.5 +/- 2.1 degrees C in D2O. Within experimental error there are no changes in deltaH(o) (25 degrees C) when D2O is substituted for H2O.     

Estimation of the maximum change in stability of globular proteins upon mutation of a hydrophobic residue to another of smaller size.

Although the hydrophobic effect is generally considered to be one of the most important forces in stabilizing the folded structure of a globular protein molecule, there is a lack of consensus on the precise magnitude of this effect. The magnitude of the hydrophobic effect is most directly measured by observing the change in stability of a protein molecule when an internal hydrophobic residue is mutated to another of smaller size. Results of such measurements have, however, been confusing because they vary greatly and are generally considerably larger than expected from the transfer free energies of corresponding small molecules. In this article, a thermodynamic argument is presented to show (1) that the variation is mainly due to that in the flexibility of the protein molecule at the site of mutation, (2) that the maximum destabilization occurs when the protein at the site of mutation is rigid, in which case the value of the destabilization is approximately given by the work of cavity formation in water, and (3) that the transfer free energy approximately gives the minimum of the range of variations. The best numerical agreements between the small molecule and the protein systems are obtained when the data from the small molecule system are expressed as the molarity-based standard free energies without other corrections.     

Changes in mechanical characteristics of immobilized crystals and films of globular proteins during substitution of D2O for H2O

The influence of substitution of the isotopic composition of the medium on the mechanical properties of immobilized crystals and films from bovine pancreatic ribonuclease and hen egg white lysozyme was investigated. The order of magnitude of the observed effects indicates that the contribution of the electrostatic interaction to the observed isotopic effect may be considered inessential. The absence of aggregation in the H2O and D2O medium under experimental conditions is demonstrated by the method of the low angle dispersion of X-rays. The observed effects of D2O on the mechanical behavior of crystals and films of proteins may be accounted for by the strengthening of molecular interactions in the samples.     

Information and redundancy in the burial folding code of globular proteins within a wide range of shapes and sizes

Recent ab initio folding simulations for a limited number of small proteins have corroborated a previous suggestion that atomic burial information obtainable from sequence could be sufficient for tertiary structure determination when combined to sequence‐independent geometrical constraints. Here, we use simulations parameterized by native burials to investigate the required amount of information in a diverse set of globular proteins comprising different structural classes and a wide size range. Burial information is provided by a potential term pushing each atom towards one among a small number L of equiprobable concentric layers. An upper bound for the required information is provided by the minimal number of layers L^min still compatible with correct folding behavior. We obtain L^min between 3 and 5 for seven small to medium proteins with 50 ≤ N_r ≤ 110 residues while for a larger protein with N_r = 141 we find that L ≥ 6 is required to maintain native stability. We additionally estimate the usable redundancy for a given L ≥ L^min from the burial entropy associated to the largest folding‐compatible fraction of “superfluous” atoms, for which the burial term can be turned off or target layers can be chosen randomly. The estimated redundancy for small proteins with L = 4 is close to 0.8. Our results are consistent with the above‐average quality of burial predictions used in previous simulations and indicate that the fraction of approachable proteins could increase significantly with even a mild, plausible, improvement on sequence‐dependent burial prediction or on sequence‐independent constraints that augment the detectable redundancy during simulations. Proteins 2016; 84:515–531. © 2016 Wiley Periodicals, Inc.     

Physical-chemical determinants of turn conformations in globular proteins

Globular proteins are assemblies of alpha-helices and beta-strands, interconnected by reverse turns and longer loops. Most short turns can be classified readily into a limited repertoire of discrete backbone conformations, but the physical-chemical determinants of these distinct conformational basins remain an open question. We investigated this question by exhaustive analysis of all backbone conformations accessible to short chain segments bracketed by either an alpha-helix or a beta-strand (i.e., alpha-segment-alpha, beta-segment-beta, alpha-segment-beta, and beta-segment-alpha) in a nine-state model. We find that each of these four secondary structure environments imposes its own unique steric and hydrogen-bonding constraints on the intervening segment, resulting in a limited repertoire of conformations. In greater detail, an exhaustive set of conformations was generated for short backbone segments having reverse-turn chain topology and bracketed between elements of secondary structure. This set was filtered, and only clash-free, hydrogen-bond-satisfied conformers having reverse-turn topology were retained. The filtered set includes authentic turn conformations, observed in proteins of known structure, but little else. In particular, over 99% of the alternative conformations failed to satisfy at least one criterion and were excluded from the filtered set. Furthermore, almost all of the remaining alternative conformations have close tolerances that would be too tight to accommodate side chains longer than a single beta-carbon. These results provide a molecular explanation for the observation that reverse turns between elements of regular secondary can be classified into a small number of discrete conformations.     

Secondary H/D isotope effect on hydrogen-bonded hydroxyl groups as a tool for recognizing distance constraints in conformational analysis of oligosaccharides

An isotopomer-selected NOE (ISNOE) method for the unequivocal identification of mutually hydrogen-bond-linked hydroxyl groups is described. It relies on the fact that the OH group's signal patterns obtained for a partially deuterated sample originate from both isotopomers of the partner hydroxyl, whereas a NOE for this group can originate from cross-relaxation with the protio isotopomer of this hydroxyl only. Hence, the isotopically shifted component of this group's signal does not appear in a ROE difference spectrum obtained with selective excitation of the partner hydroxyl. This method is also applicable in those cases when only one of two mutually hydrogen-bonded groups exhibits resolvable isotope shifts. Furthermore, it is shown that isotope shifts may occur even for pairs of OH groups that are not mutually hydrogen-bonded, if these participate in hydrogen bonds with other hydroxyls and thereby affect conformational equilibria. The ISNOE experiment enables one to distinguish between these two sources of isotope shifts. Since the OO distance for hydrogen-bonded hydroxyls in sugars is known to lie between 2.7 and 3.0 Å , the hydrogen bonds established by ISNOE can be used in conformational analysis as reliable, motionally non-averaged distance constraints for the conformations containing these bonds.     

'Molten-globule state': a compact form of globular proteins with mobile side-chains

In rat lacrimal glands, Forskolin induces a dose-dependent [³H]protein release. This effect can be potentiated by papaverine. As for the other inducers whose effects on protein secretion are assumed to be cAMP-mediated, Forskolin secretion time course shows a latency. Isoproterenol decreases the Forskolin EC_5- at least 60--times. On the other hand, Forskolin enhances the efficacy of isoproterenol without affecting its potency. As a whole, the data collected show that isoproterenol-induced [³H]protein secretion in rat lacrimal glands involved adenylate cyclase activation by coupling with β-adrenergic receptors.     

Prediction of probable pathways of folding in globular proteins

A method is described for the prediction of probable folding pathways of globular proteins, based on the analysis of distance maps. It is applicable to proteins of unknown spatial structure but known amino acid sequence as well as to proteins of known structure. It is based on an objective procedure for the determination of the boundary of compact regions that contain high densities of interresidue contacts on the distance map of a globular protein. The procedure can be used both with contact maps derived from a known three-dimensional protein structure and with predicted contact maps computed by means of a statistical procedure from the amino acid sequence alone. The computed contact map can also be used to predict the location of compact short-range structures, viz. -helices and -turns, thereby complementing other statistical predictive procedures. The method provides an objective basis for the derivation of a theoretically predicted pathway of protein folding, proposed by us earlier [Tanaka and Scheraga (1977) Macromolecules10, 291–304; Némethy and Scheraga (1979) Proc. Natl. Acad. Sci., U.S.A.76, 6050–6054].     

The effect of deuterium oxide on the stability of the collagen model peptides H‐(Pro‐Pro‐Gly)10‐OH,H‐(Gly‐Pro‐4(R)Hyp)9‐OH,and Type I collagen

The collagen triple helix has a larger accessible surface area per molecular mass than globular proteins, and therefore potentially more water interaction sites. The effect of deuterium oxide on the stability of collagen model peptides and Type I collagen molecules was analyzed by circular dichroism and differential scanning calorimetry. The transition temperatures (T_m) of the protonated peptide (Pro‐Pro‐Gly)₁₀ were 25.4 and 28.7°C in H₂O and D₂O, respectively. The increase of the T_m of (Pro‐Pro‐Gly)₁₀ measured calorimetrically at 1.0°C min^?1 in a low pH solution from the protonated to the deuterated solvent was 5.1°C. The increases of the T_m for (Gly‐Pro‐4(R)Hyp)₉ and pepsin‐extracted Type I collagen were measured as 4.2 and 2.2°C, respectively. These results indicated that the increase in the T_m in the presence of D₂O is comparable to that of globular proteins, and much less than reported previously for collagen model peptides [Gough and Bhatnagar, J Biomol Struct Dyn 1999, 17, 481–491]. These experimental results suggest that the interaction of water molecules with collagen is similar to the interaction of water with globular proteins, when the ratio of collagen to water is very small and collagen is monomerically dispersed in the solvent. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 93–101, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Simulating the minimum core for hydrophobic collapse in globular proteins.

To investigate the nature of hydrophobic collapse considered to be the driving force in protein folding, we have simulated aqueous solutions of two model hydrophobic solutes, methane and isobutylene. Using a novel methodology for determining contacts, we can precisely follow hydrophobic aggregation as it proceeds through three stages: dispersed, transition, and collapsed. Theoretical modeling of the cluster formation observed by simulation indicates that this aggregation is cooperative and that the simulations favor the formation of a single cluster midway through the transition stage. This defines a minimum solute hydrophobic core volume. We compare this with protein hydrophobic core volumes determined from solved crystal structures. Our analysis shows that the solute core volume roughly estimates the minimum core size required for independent hydrophobic stabilization of a protein and defines a limiting concentration of nonpolar residues that can cause hydrophobic collapse. These results suggest that the physical forces driving aggregation of hydrophobic molecules in water is indeed responsible for protein folding.     

Distance-constraint approach to higher-order structures of globular proteins with empirically determined distances between amino acid residues

An analysis of higher-order structures of globular proteins by means of a distance-constraint approach is presented. Conformations are generated for each of 21 test proteins of small and medium sizes by optimizing an objective functionf=w _ij(d _ij–d _ij)², whered _ij is a distance between residuesi andj in a calculated conformation, d _ij is an assigned distance to the (ij) pair of residues which is determined based on the statistics of known three-dimensional structures of 14 proteins in the earlier study, andw _ij is a weighting factor. d _ij involves information about hydrophobicity and hydrophilicity of each amino acid residue and about connectivity of a polypeptide chain. In these calculations, only the amino acid sequence is used as input data specific to a calculated protein. With respect to higher-order structures regenerated in the optimized conformations, the following properties are analyzed: (a) N14 of a residue, defined as the number of residues surrounding the residue located within a sphere of radius of 14 Å; (b) root-mean-square differences of the global and local conformations from the corresponding X-ray conformations; (c) distance profiles in the short and medium ranges; and (d) distance maps. The effects of supplementary information about locations of secondary structures and disulfide bonds are also examined to discuss the potential ability of this methodology to predict the three-dimensional structures of globular proteins.     

Solvent denaturation of proteins as observed by resolution-enhanced Fourier transform infrared spectroscopy

Fourier self-deconvolution of Fourier transform infrared (FTIR) spectra and second derivative FTIR spectroscopy were applied to study solvent-induced conformational changes in globular proteins. For beta-lactoglobulin a total of three different denatured forms were identified in alkaline solution and in aqueous methanol-d1 and isopropanol-d1. In isopropanol-d1 solution a new conformation was identified which appears to resemble, but is not identical with, the beta-structure of native proteins. This conformation is characterized by absorption bands around 1615-1618 and 1684-1688 cm-1, and is also observed for concanavalin A and chymotrypsinogen A in aqueous isopropanol-d1 solution.     

Kinetics and thermodynamics of thermal denaturation in acyl carrier protein.

The denaturation of Escherichia coli acyl carrier protein (ACP) in buffers containing both monovalent and divalent cations was followed by variable-temperature NMR and differential scanning calorimetry. Both high concentrations of monovalent salts (Na+) and moderate concentrations of divalent salts (Ca2+) raise the denaturation temperature, but calorimetry indicates that a significant increase in the enthalpy of denaturation is obtained only with the addition of a divalent salt. NMR experiments in both low ionic strength monovalent buffers and low ionic strength monovalent buffers containing calcium ions show exchange between native and denatured forms to be slow on the NMR time scale. However, in high ionic strength monovalent buffers, where the temperature of denaturation is elevated as it is in the presence of Ca2+, the transition is fast on the NMR time scale. These results suggest that monovalent and divalent cations may act to stabilize ACP in different ways. Monovalent ions may nonspecifically balance the intrinsic negative charge of this protein in a way that is similar for native, denatured, and intermediate forms. Divalent cations provide stability by binding to specific sites present only in the native state.     

Stability and stabilization of globular proteins in solution

Proteins are multifunctional: their amino acid sequences simultaneously determine folding, function and turnover. Correspondingly, evolution selected for compromises between rigidity (stability) and flexibility (folding/function/degradation), to the result that generally the free energy of stabilization of globular proteins in solution is the equivalent to only a few weak intermolecular interactions. Additional increments may come from extrinsic factors such as ligands or specific compatible solutes. Apart from the enthalpic effects, entropy may play a role by reducing the flexibility (cystine bridges, increased proline content), or by water release from residues buried upon folding and association. Additional quaternary interactions and closer packing are typical characteristics of proteins from thermophiles. In halophiles, protein stability and function are maintained by increased ion binding and glutamic acid content, both allowing the protein inventory to compete for water at high salt. Acidophiles and alkalophiles show neutral intracellular pH; proteins facing the outside extremes of pH possess anomalously high contents in ionizable amino acids. Global comparisons of the amino acid compositions and sequences of proteins from mesophiles and extremophiles did not result in general rules of protein stabilization, even after including complete genome sequences into the search. Obviously, proteins are individuals that optimize internal packing and external solvent interactions by very different mechanisms, each protein in its own way. Strategies deduced from specific ultrastable proteins allow stabilizing point mutations to be predicted.     

On the effect of hydrostatic pressure on the conformational stability of globular proteins

The model developed for cold denaturation (Graziano, PCCP 2010, 12, 14245‐14252) is extended to rationalize the dependence of protein conformational stability upon hydrostatic pressure, at room temperature. A pressure− volume work is associated with the process of cavity creation for the need to enlarge the liquid volume against hydrostatic pressure. This contribution destabilizes the native state that has a molecular volume slightly larger than the denatured state due to voids existing in the protein core. Therefore, there is a hydrostatic pressure value at which the pressure−volume contribution plus the conformational entropy loss of the polypeptide chain are able to overwhelm the stabilizing gain in translational entropy of water molecules, due to the decrease in water accessible surface area upon folding, causing denaturation. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 711–718, 2015.     

Native atomic burials, supplemented by physically motivated hydrogen bond constraints, contain sufficient information to determine the tertiary structure of small globular proteins

We investigate the possibility that atomic burials, as measured by their distances from the structural geometrical center, contain sufficient information to determine the tertiary structure of globular proteins. We report Monte Carlo simulated annealing results of all-atom hard-sphere models in continuous space for four small proteins: the all-beta WW-domain 1E0L, the alpha/beta protein-G 1IGD, the all-alpha engrailed homeo-domain 1ENH, and the alpha + beta engineered monomeric form of the Cro protein 1ORC. We used as energy function the sum over all atoms, labeled by i, of |R(i) - R(i) (*)|, where R(i) is the atomic distance from the center of coordinates, or central distance, and R(i) (*) is the "ideal" central distance obtained from the native structure. Hydrogen bonds were taken into consideration by the assignment of two ideal distances for backbone atoms forming hydrogen bonds in the native structure depending on the formation of a geometrically defined bond, independently of bond partner. Lowest energy final conformations turned out to be very similar to the native structure for the four proteins under investigation and a strong correlation was observed between energy and distance root mean square deviation (DRMS) from the native in the case of all-beta 1E0L and alpha/beta 1IGD. For all alpha 1ENH and alpha + beta 1ORC the overall correlation between energy and DRMS among final conformations was not as high because some trajectories resulted in high DRMS but low energy final conformations in which alpha-helices adopted a non-native mutual orientation. Comparison between central distances and actual accessible surface areas corroborated the implicit assumption of correlation between these two quantities. The Z-score obtained with this native-centric potential in the discrimination of native 1ORC from a set of random compact structures confirmed that it contains a much smaller amount of native information when compared to a traditional contact Go potential but indicated that simple sequence-dependent burial potentials still need some improvement in order to attain a similar discriminability. Taken together, our results suggest that central distances, in conjunction to physically motivated hydrogen bond constraints, contain sufficient information to determine the native conformation of these small proteins and that a solution to the folding problem for globular proteins could arise from sufficiently accurate burial predictions from sequence followed by minimization of a burial-dependent energy function.     

Hydrophobicity of amino acid subgroups in proteins

Protein folding studies often utilize areas and volumes to assess the hydrophobic contribution to conformational free energy (Richards, F.M. Annu. Rev. Biophys. Bioeng. 6:151-176, 1977). We have calculated the mean area buried upon folding for every chemical group in each residue within a set of X-ray elucidated proteins. These measurements, together with a standard state cavity size for each group, are documented in a table. It is observed that, on average, each type of group buries a constant fraction of its standard state area. The mean area buried by most, though not all, groups can be closely approximated by summing contributions from three characteristic parameters corresponding to three atom types: (1) carbon or sulfur, which turn out to be 86% buried, on average; (2) neutral oxygen or nitrogen, which are 40% buried, on average; and (3) charged oxygen or nitrogen, which are 32% buried, on average.