Disulfide bonds and the stability of globular proteins.

An understanding of the forces that contribute to stability is pivotal in solving the protein-folding problem. Classical theory suggests that disulfide bonds stabilize proteins by reducing the entropy of the denatured state. More recent theories have attempted to expand this idea, suggesting that in addition to configurational entropic effects, enthalpic and native-state effects occur and cannot be neglected. Experimental thermodynamic evidence is examined from two sources: (1) the disruption of naturally occurring disulfides, and (2) the insertion of novel disulfides. The data confirm that enthalpic and native-state effects are often significant. The experimental changes in free energy are compared to those predicted by different theories. The differences between theory and experiment are large near 300 K and do not lend support to any of the current theories regarding the stabilization of proteins by disulfide bonds. This observation is a result of not only deficiencies in the theoretical models but also from difficulties in determining the effects of disulfide bonds on protein stability against the backdrop of numerous subtle stabilizing factors (in both the native and denatured states), which they may also affect.     

Structural principles of alpha/beta barrel proteins: the packing of the interior of the sheet

Alpha/beta barrel structures very similar to that first observed in triose phosphate isomerase are now known to occur in 14 enzymes. To understand the origin of this fold, we analyzed in three of these proteins the geometry of the eight-stranded beta-sheets and the packing of the residues at the center of the barrel. The packing in this region is seen in its simplest form in glycolate oxidase. It consists of 12 residues arranged in three layers. Each layer contains four side chains. The packing of RubisCO and TIM can be understood in terms of distortions of this simple pattern, caused by residues with small side chains at some of the positions inside the barrel. Two classes of packing are found. In one class, to which RubisCO and TIM belong, the central layer is formed by a residue from the first, third, fifth, and seventh strands; the upper and lower layers are formed by residues from the second, fourth, sixth, and eighth strands. In the second class, to which GAO belongs, this is reversed: it is side chains from the even-numbered strands that form the central layer, and side chains from the odd-numbered strands that form the outer layers. Our results suggest that not all proteins with this fold are related by evolution, but that they represent a common favorable solution to the structural problems involved in the creation of a closed beta barrel.     

Dynamic polyhedral models of globular proteins

We have devised several mechanical models of globular proteins by approximating them to various polyhedra (dodecahedron, truncated octahedron, icosahedron, truncated icosahedron). The models comprise hollow blocks linked together in a flexible chain. Between blocks there is a set of several reversible, weak magnetic interactions such that when the chain is agitated, it will fold into a stable polyhedral structure about the size of a hand. Folding may be followed in real time with a video camera. Key to the success of the folding process is the lightness of the chain. Several side chains may also be added to the blocks such that they come together to create a polyhedral core when the chain folds. The models have a number of similarities to globular proteins: each chain folds into a unique, but dynamic, three-dimensional structure; the instructions that determine this structure are built into the configuration of blocks; and it is difficult to predict this structure given the unfolded block configuration. Furthermore, the chains fold quickly, generally in less than a minute, several pathways are involved, and these pathways progress through elements of "native" structure. In particular, the models emphasize the importance of restricted conformational mobility in assisting the chain to fold, and also in eliminating undesirable interactions. Because of these similarities to globular proteins, we believe that the polyhedral models will, with continued development, be helpful in understanding the protein folding process, while at the same time acting as valuable educational visual aids. They might also inspire the construction of new types of microscopic, self-assembling devices.     

Analysis of anisotropic side-chain packing in proteins and application to high-resolution structure prediction

pi-pi, Cation-pi, and hydrophobic packing interactions contribute specificity to protein folding and stability to the native state. As a step towards developing improved models of these interactions in proteins, we compare the side-chain packing arrangements in native proteins to those found in compact decoys produced by the Rosetta de novo structure prediction method. We find enrichments in the native distributions for T-shaped and parallel offset arrangements of aromatic residue pairs, in parallel stacked arrangements of cation-aromatic pairs, in parallel stacked pairs involving proline residues, and in parallel offset arrangements for aliphatic residue pairs. We then investigate the extent to which the distinctive features of native packing can be explained using Lennard-Jones and electrostatics models. Finally, we derive orientation-dependent pi-pi, cation-pi and hydrophobic interaction potentials based on the differences between the native and compact decoy distributions and investigate their efficacy for high-resolution protein structure prediction. Surprisingly, the orientation-dependent potential derived from the packing arrangements of aliphatic side-chain pairs distinguishes the native structure from compact decoys better than the orientation-dependent potentials describing pi-pi and cation-pi interactions.     

The energetics of specific binding of AT-hooks from HMGA1 to target DNA

The interaction of the second and third AT-hooks of HMGA1 (formerly HMGI/Y), which bind selectively in the minor groove of an AT-rich DNA sequence, was studied at different temperatures and ionic strengths by spectropolarimetry, spectrofluorimetry, isothermal titration calorimetry and differential scanning calorimetry. The data show that binding of the ten amino acid core element of the two AT-hooks, which penetrates deep into the minor groove, is entropically driven: both the entropy and enthalpy of association of the peptides to the target DNA are positive up to 50 degrees C. The seven amino acid extension of the core in the second AT-hook, which extends out from the minor groove and loops over the phosphodiester backbone, adds a substantial negative enthalpic component into the binding of the 17 residue DBD2 peptide to DNA that corresponds in magnitude to the enthalpy of formation of two hydrogen bonds. The ionic strength dependence of the association constant allowed an estimation of the electrostatic component of binding and, by subtraction, the contribution of the non-electrostatic component, which results from dehydration of the contacting surfaces and makes up almost 70% of the total energy of complex formation. The exceptionally large positive entropy and enthalpy of association of the core AT-hook peptides with target DNA suggest that the water, which is removed from the minor groove of DNA upon binding, is in a highly ordered state. Acetylation of the lysine residue in the second AT-hook, which corresponds to Lys65 of HMGA1, has little effect on the DNA binding; so it appears that repression of the hIFNbeta gene, which follows this modification, is not a direct result of the abrogation of DNA binding.     

A continuum model to study the internal fluctuations in globular proteins

The salient features of the differential equation model to study protein dynamics are presented with results for 19 proteins.     

Analysis and prediction of inter-strand packing distances between beta-sheets of globular proteins

Any two beta-strands belonging to two different beta-sheets in a protein structure are considered to pack interactively if each beta-strand has at least one residue that undergoes a loss of one tenth or more of its solvent contact surface area upon packing. A data set of protein 3-D structures (determined at 2.5 A resolution or better), corresponding to 428 protein chains, contains 1986 non-identical pairs of beta-strands involved in interactive packing. The inter-axial distance between these is significantly correlated to the weighted sum of the volumes of the interacting residues at the packing interface. This correlation can be used to predict the changes in the inter-sheet distances in equivalent beta-sheets in homologous proteins and, therefore, is of value in comparative modelling of proteins.     

Structural mobility and transformations in globular proteins

Journal of Biosciences - Although globular proteins are endowed with well defined three-dimensional structures, they exhibit substantial mobility within the framework of the given three-dimensional...     

Hierarchic organization of globular proteins: Experimental studies on thermolysin

Previous studies from this laboratory have shown that the thermolysin fragment 121–316, comprising entirely the“all-α” COOH-terminal structural domain 158–316, as well as fragment 206–316 (fragment FII) are able to refold into a native-like, stable structure independently from the rest of the protein molecule. The present report describes conformational properties of fragments 228–316 and 255–316 obtained by chemical and enzymatic cleavage of fragment FII, respectively. These subfragments are able to acquire a stable conformation of native-like characteristics, as judged by quantitative analysis of secondary structure from far-ultra-violet circular dichroism spectra and immunochemical properties using rabbit anti-thermolysin antibodies. Melting curves of the secondary structure of the fragments show cooperativity with a temperature of half-denaturationT _mof 65–66°C. The results of this study provide evidence that it is possible to isolate stable supersecondary structures (folding units) of globular proteins and correlate well with predictions of subdomains of the COOH-terminal structural domain 158–316 of thermolysin.     

Families and the structural relatedness among globular proteins.

Protein structures come in families. Are families “closely knit” or “loosely knit” entities? We describe a measure of relatedness among polymer conformations. Based on weighted distance maps, this measure differs from existing measures mainly in two respects: (1) it is computationally fast, and (2) it can compare any two proteins, regardless of their relative chain lengths or degree of similarity. It does not require finding relative alignments. The measure is used here to determine the dissimilarities between all 12, 403 possible pairs of 158 diverse protein structures from the Brookhaven Protein Data Bank (PDB). Combined with minimal spanning trees and hierarchical clustering methods, this measure is used to define structural families. It is also useful for rapidly searching a dataset of protein structures for specific substructural motifs. By using an analogy to distributions of Euclidean distances, we find that protein families are not tightly knit entities.     

Structural and energetic consequences of mutations in a solvated hydrophobic cavity

The structural and energetic consequences of modifications to the hydrophobic cavity of interleukin 1-beta (IL-1beta) are described. Previous reports demonstrated that the entirely hydrophobic cavity of IL-1beta contains positionally disordered water. To gain a better understanding of the nature of this cavity and the water therein, a number of mutant proteins were constructed by site-directed mutagenesis, designed to result in altered hydrophobicity of the cavity. These mutations involve the replacement of specific phenylalanine residues, which circumscribe the cavity, with tyrosine, tryptophan, leucine and isoleucine. Using differential scanning calorimetry to determine the relative stabilities of the wild-type and mutant proteins, we found all of the mutants to be destabilizing. X-ray crystallography was used to identify the structural consequences of the mutations. No clear correlation between the hydrophobicities of the specific side-chains introduced and the resulting stabilities was found.     

Influence of proline residues in transmembrane helix packing

Integral membrane proteins often contain proline residues in their alpha-helical transmembrane (TM) fragments, which may strongly influence their folding and association. Pro-scanning mutagenesis of the helical domain of glycophorin A (GpA) showed that replacement of the residues located at the center abrogates helix packing while substitution of the residues forming the ending helical turns allows dimer formation. Synthetic TM peptides revealed that a point mutation of one of the residues of the dimerization motif (L75P) located at the N-terminal helical turn of the GpA TM fragment, adopts a secondary structure and oligomeric state similar to the wild-type sequence in detergents. In addition, both glycosylation mapping in biological membranes and molecular dynamics showed that the presence of a proline residue at the lipid/water interface has as an effect the extension of the helical end. Thus, helix packing can be an important factor that determines appearance of proline in TM helices. Membrane proteins might accumulate proline residues at the two ends of their TM segments in order to modulate the exposition of key amino acid residues at the interface for molecular recognition events while allowing stable association and native folding.     

Core side-chain packing and backbone conformation in Lpp-56 coiled-coil mutants

Native proteins exhibit precise geometric packing of atoms in their hydrophobic interiors. Nonetheless, controversy remains about the role of core side-chain packing in specifying and stabilizing the folded structures of proteins. Here we investigate the role of core packing in determining the conformation and stability of the Lpp-56 trimerization domain. The X-ray crystal structures of Lpp-56 mutants with alanine substitutions at two and four interior core positions reveal trimeric coiled coils in which the twist of individual helices and the helix-helix spacing vary significantly to achieve the most favored superhelical packing arrangement. Introduction of each alanine "layer" into the hydrophobic core destabilizes the superhelix by 1.4 kcal mol(-1). Although the methyl groups of the alanine residues pack at their optimum van der Waals contacts in the coiled-coil trimer, they provide a smaller component of hydrophobic interactions than bulky hydrophobic side-chains to the thermodynamic stability. Thus, specific side-chain packing in the hydrophobic core of coiled coils are important determinants of protein main-chain conformation and stability.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Identification of compounds with binding affinity to proteins via magnetization transfer from bulk water*

A powerful screening by NMR methodology (WaterLOGSY), based on transfer of magnetization from bulk water, for the identification of compounds that interact with target biomolecules (proteins, RNA and DNA fragments) is described. The method exploits efficiently the large reservoir of H₂O magnetization. The high sensitivity of the technique reduces the amount of biomolecule and ligands needed for the screening, which constitutes an important requirement for high throughput screening by NMR of large libraries of compounds. Application of the method to a compound mixture against the cyclin-dependent kinase 2 (cdk2) protein is presented.     

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.     

Spatial optimization of electrostatic interactions between the ionized groups in globular proteins

A model approach is suggested to estimate the degree of spatial optimization of the electrostatic interactions in protein molecules. The method is tested on a set of 44 globular proteins, representative of the available crystallographic data. The theoretical model is based on macroscopic computation of the contribution of charge–charge interactions to the electrostatic term of the free energy for the native proteins and for a big number of virtual structures with randomly distributed on protein surface charge consetellations (generated by a Monte-Carlo technique). The statistical probability of occurrence of random structures with electrostatic energies lower than the energy of the native protein is suggested as a criterion for spatial optimization of the electrostatic interactions. The results support the hypothesis that the folding process optimizes the stabilizing effect of electrostatic interactions, but to very different degree for different proteins. A parallel analysis of ion pairs shows that the optimization of the electrostatic term in globular proteins has increasingly gone in the direction of rejecting the repulsive short contacts between charges of equal sign than of creating of more salt bridges (in comparison with the statistically expected number of shortrange ion pairs in the simulated random structures). It is observed that the decrease in the spatial optimization of the electrostatic interactions is usually compensated for by an appearance of disulfide bridges in the covalent structure of the examined proteins. © 1994 Wiley-Liss, Inc.