Hydrophobicity of transmembrane proteins: spatially profiling the distribution

A hallmark of soluble globular protein tertiary structure is a hydrophobic core and a protein exterior populated predominantly by hydrophilic residues. Recent hydrophobic moment profiling of the spatial distribution of 30 globular proteins of diverse size and structure had revealed features of this distribution that were comparable. Analogous profiling of the hydrophobicity distribution of the alpha-helical buried bundles of several transmembrane proteins, as the lipid/protein interface is approached from within the bilayer, reveals spatial hydrophobicity profiles that contrast with those obtained for the soluble proteins. The calculations, which enable relative changes of hydrophobicity to be simply identified over the entire spatial extent of the multimer within the lipid bilayer, show the accumulated zero-order moments of the bundles to be mainly inverted with respect to that found for the soluble proteins. This indicates a statistical increase in the average residue hydrophobic content as the lipid bilayer is approached. This result differs from that of a relatively recent calculation and qualitatively agrees with earlier calculations involving lipid exposed and buried residues of the alpha-helices of transmembrane proteins. Spatial profiling, over the entire spatial extent of the multimer with scaled values of residue hydrophobicity, provides information that is not available from calculations using lipid exposure alone.     

Spatial assignment of amino acid residues in globular proteins: An approach from information theory

In the native folded state of globular proteins, amino acid residues place themselves at various positions from the centroid of the molecule. Applying information theory on 19 protein crystals the spatial preferences have been found out from the frequencies of occurrence of residues within various concentric ellipsoidal zones of proteins. The intrinsic spatial preferences of individual residues are related to their physical and chemical properties. The directing power of the individual residues on the chain path and the spatial information contained by doublets of residues have been found out. The derived information is used to predict the spatial/zonal preference of residues in carp myogen using the knowledge of amino acid sequence. The implication of packing densities in different spatial zones are discussed.     

A New Spatial Representation of Amino Acid Residues in Globular Proteins

Abstract

For the globular proteins with known three-dimensional structures, an ellipsoid model of each protein was constructed with least volume and its dimensions were derived. The spatial arrangements were made for the C_α and side chain atoms of that protein within that ellipsoid. This new spatial representation shows the residue position from the centroid, as well as the depth from the surface. The average spatial parameters were then calculated. The correlations between these new spatial parameters and the existing parameters of the amino acid residues were then derived.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

General architecture of the alpha-helical globule

A model is presented for the arrangement of alpha-helices in globular proteins. In the model, helices are placed on certain ribs of "quasi-spherical" polyhedra. The polyhedra are chosen so as to allow the close packing of helices around a hydrophobic core and to stress the collective interactions of the individual helices. The model predicts a small set of stable architectures for alpha-helices in globular proteins and describes the geometries of the helix packings. Some of the predicted helix arrangements have already been observed in known protein structures; others are new. An analysis of the three-dimensional structures of all proteins for which co-ordinates are available shows that the model closely approximates the arrangements and packing of helices actually observed. The average deviations of the real helix axes from those in the model polyhedra is +/- 20 degrees in orientation and +/- 2 A in position (1 A = 0.1 nm). We also show that for proteins that are not homologous, but whose helix arrangements are described by the same polyhedron, the root-mean-square difference in the position of the C alpha atoms in the helices is 1.6 to 3.0 A.     

Mammalian hemodynamics: a new similarity principle

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

The spatial distribution of physical, chemical, energetic and conformational properties of amino acid residues in globular proteins

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

Prediction of exposure degree diagram and sites of limited proteolysis in globular proteins as an approach to computer-aided design of protein bioregulators with prolonged action

In order to prolong the lifetime of protein bioregulators in blood it is possible to engineer analogs with protected sites of limited proteolysis. To determine the sites, primarily accessible to trypsin-like proteases, a computer procedure has been developed, including a prediction algorithm, to produce the residue diagram of a globular protein and a discriminant algorithm to determine the sites most liable to proteolysis. The accuracy of prediction of amino acid residue exposure is characterised by correlation coefficients between experimental and theoretical exposure values, the coefficients being about 0.7 as calculated for 10 globular proteins. The classification of Arg and Lys residues into two groups, susceptible or insusceptible to protease, has an error percentage of about 25.     

Amino acid composition and thermal stability of proteins

This study investigates the relationship between the thermal stability of a globular protein and its amino acid composition. The method deals with the relationship between the amino acid compositions and melting points in a set of proteins by computing single-residue and group correlations. Groups of residues are shown to stabilize or destabilize the molecule against temperature. The stabilizing group consists of polar-charged residues and nonpolar residues possessing high surrounding hydrophobicity. The polar-uncharged residues destabilize the molecule against temperature, serine being the most destabilizing residue. A very high cooperativity exists among the stabilizing nonpolar residues suggesting that their characteristic clustering inside the globule may enhance the thermostability of a protein. In small globular proteins which act as single cooperative units, the melting temperature remains mainly a function of amino acid composition, whereas in complex molecules it depends on other factors also.     

Prediction of natively unfolded regions in protein chain

We have shown that the ability of a protein to be in globular or in natively unfolded state (under native conditions) may be determined (besides low overall hydrophobicity and a large net charge) by such a property as the average environment density, the average number of residues enclosed at the given distance. A statistical scale of the average number of residues enclosed at the given distance for 20 types of amino acid residues in globular state has been created on the basis of 6626 protein structures. Using this scale for separation of 80 globular and 90 natively unfolded proteins we fail only in 11% of proteins (compared with 17% of errors which are observed if to use hydrophobicity scale). The present scale may be used both for prediction of form (folded or unfolded) of the native state of protein and for prediction of natively unfolded regions in protein chains. The results of comparison of our method of predicting natively unfolded regions with the other known methods show that our method has the highest fraction of correctly predicted natively unfolded regions (that is 87% and 77% if to make averaging over residues and over proteins correspondingly).     

Effect of spatial coupling of amino acid residues on the possibility of mutational substitution in globular proteins

Amino acid residues can be divided into similar groups by frequencies of interreplacements in the evolutionary pathway and by trends to spatial contacts at the tertiary structures of globular proteins. Each residue was compared to the cluster of spatial surrounding--the totality of residues spacially drawn together. 5210 clusters in 32 unhomologous proteins with established tertiary structure and 6447 clusters formed only by variables amino acid residues were analysed. Spatial contacts among residues were studied depending on the secondary structure and the amount of residues in a cluster. It was assumed that functionally admissible mutations may be defined, first of all, by the degree of neighboring of amino acid residues in the spatial surrounding.     

Establishment of Constraints on Amyloid Formation Imposed by Steric Exclusion of Globular Domains

In many disease-related and functional amyloids, the amyloid-forming regions of proteins are flanked by globular domains. When located in close vicinity of the amyloid regions along the chain, the globular domains can prevent the formation of amyloids because of the steric repulsion. Experimental tests of this effect are few in number and non-systematic, and their interpretation is hampered by polymorphism of amyloid structures. In this situation, modeling approaches that use such a clear-cut criterion as the steric tension can give us highly trustworthy results. In this work, we evaluated this steric effect by using molecular modeling and dynamics. As an example, we tested hybrid proteins containing an amyloid-forming fragment of Aβ peptide (17–42) linked to one or two globular domains of GFP. Searching for the shortest possible linker, we constructed models with pseudo-helical arrangements of the densely packed GFPs around the Aβ amyloid core. The molecular modeling showed that linkers of 7 and more residues allow fibrillogenesis of the Aβ-peptide flanked by GFP on one side and 18 and more residues when Aβ-peptide is flanked by GFPs on both sides. Furthermore, we were able to establish a more general relationship between the size of the globular domains and the length of the linkers by using analytical expressions and rigid body simulations. Our results will find use in planning and interpretation of experiments, improvement of the prediction of amyloidogenic regions in proteins, and design of new functional amyloids carrying globular domains.     

Variation of amino acid properties in all-beta globular and outer membrane protein structures

Discriminating outer membrane (OM) proteins from globular proteins is an important task. The structural analysis of β-strands dominating globular (all-β) proteins and OM proteins provides useful insight to distinguish between them. In this work, we analyze the characteristic features of the 20 amino acid residues in all-β and OM proteins. We set up numerical indices for several properties of amino acid residues, such as, conformational parameters, surrounding hydrophobicity, accessible surface area and reduction in accessibility, and inter-residue contacts. We found that all the aromatic residues prefer to be in β-strands of both globular and OM proteins. The surrounding hydrophobicity of aromatic and non-polar amino acid residues in globular proteins is significantly higher than that of OM proteins. The residues Trp, Arg, Phe and Gln show a remarkable difference of reduction in accessibility between all-β globular (βG) and OM proteins. The positively charged residues, Lys and Arg in the membrane part of OM proteins have more number of contacts than globular proteins. Further, the behavior of the 20 amino acid residues in β-strand segments of globular and OM proteins have been discussed. The parameters developed in this work can be used for identifying transmembrane β-strands in OM proteins and for discriminating βG proteins from OM proteins.     

Distribution of amino acid residues and residue-residue contacts in molecular chaperones.

The amino acid distribution and residue-residue contacts in molecular chaperones are different when compared to normal globular proteins. The study of molecular chaperones reveals a different surrounding environment to exist for the residues Cys, Trp, and His which may play an important role in determining the chaperone structures. Unlike globular proteins, it has been observed that a one-to-one correspondence between the amino acid distribution in a sequence and the structures of molecular chaperones. The preference of amino acid residues surrounding all 20 types of residues in secondary structures and their accessible surface areas have been analysed.     

Three-dimensional structure of porcine procarboxypeptidase B: a structural basis of its inactivity.

Procarboxypeptidase B is converted to enzymatically active carboxypeptidase B by limited proteolysis catalysed by trypsin, removing the long N-terminal activation segment of 95 amino acids. The three-dimensional crystal structure of procarboxypeptidase B from porcine pancreas has been determined at 2.3 A resolution and refined to a crystallographic R-factor of 0.169. The functional determinants of its enzymatic inactivity and of its activation by limited proteolysis have thus been unveiled. The activation segment folds in a globular region with an open sandwich antiparallel-alpha antiparallel-beta topology and in a C terminal alpha-helix which connects it to the enzyme moiety. The globular region (A7-A82) shields the preformed active site, and establishes specific interactions with residues important for substrate recognition. AspA41 forms a salt bridge with Arg145, which in active carboxypeptidase binds the C-terminal carboxyl group of substrate molecules. The connecting region occupies the putative extended substrate binding site. The scissile peptide bond cleaved by trypsin during activation is very exposed. Its cleavage leads to the release of the activation segment and to exposure of the substrate binding site. An open-sandwich folding has been observed in a number of other proteins and protein domains. One of them is the C-terminal fragment of L7/L12, a ribosomal protein from Escherichia coli that displays a topology similar to the activation domain of procarboxypeptidase.     

Solid model compounds and the thermodynamics of protein unfolding.

Analysis of thermodynamic data on the dissolution of solid cyclic dipeptides into water in terms of group additivity provides a rationale for the enthalpy and entropy convergence temperatures observed for small globular protein denaturation and the dissolution of model compounds into water. Convergence temperatures are temperatures at which the extrapolated enthalpy or entropy changes for a series of related compounds take on a common value. At these temperatures (TH* and TS*) the apolar contributions to the corresponding thermodynamic values (delta H degrees and delta S degrees) are shown to be zero. Other contributions such as hydrogen bonding and configurational effects can then be evaluated and their quantitative effects on the stability of globular proteins assessed. It is shown that the denaturational heat capacity is composed of a large positive contribution from the exposure of apolar groups and a significant negative contribution from the exposure of polar groups in agreement with previous results. The large apolar contribution suggests that a liquid hydrocarbon model of the hydrophobic effect does not accurately represent the apolar contribution to delta H degrees of denaturation. Rather, significant enthalpic stabilizing contributions are found to arise from peptide groups (hydrogen bonding). Combining the average structural features of globular proteins (i.e. number of residues, fraction of buried apolar groups and fraction of hydrogen bonds) with their specific group contributions permits a first-order prediction of the thermodynamic properties of proteins. The predicted values compare well with literature values for cytochrome c, myoglobin, ribonuclease A and lysozyme. The major thermodynamic features are described by the number of peptide and apolar groups in a given protein.     

Native protein sequences are designed to destabilize folding intermediates

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.     

Approaching a complete classification of protein secondary structure

A complete classification of types of the protein secondary structure is developed on the basis of computer analysis of the crystallographic structural data deposited in the protein Data Bank. The majority of amino acid residues fall into five conformation types. A conclusion is drawn that the number of sequence variants of torsion angles phi, psi in globular proteins is limited and is essentially less than the number of possible amino acid sequences for this chain length. Along with alpha-helix and beta-structure, the distribution analysis assigning every maximum of distribution of amino acid conformations on Ramachandran map to a certain type of the secondary structure exposed a third type of the secondary structure that was previously neglected. This type of the structure is extended left-handed helical conformation, designated as mobile (M-) conformation. A full set of M-conformation fragments that seems to play a major role in protein globule dynamics has been obtained, a small radius of correlation for the polypeptide chain in M-conformation is demonstrated. It explains a prevalence of short segments of mobile conformation revealed in globular proteins. For secondary structure types, the frequency of occurrence of amino acid residues has been computed.     

Protein sequences yield a proteomic code

Analysis of crystallized protein structures suggests that globular proteins are organized as consecutively connected units of 25-35 residues. These units are closed loops, that is returns of the polypeptide chain trajectory to a close contact with itself. This universal feature of apparently polymer-statistical nature is a basis for a principally novel view on the globular proteins as loop fold structures. The same unit size has been detected in protein sequences translated from complete prokaryotic genomes by positional autocorrelation analysis, which strongly indicates the evolutionary connection of the units. The units are further characterized by prototype sequences matching to their numerous derivatives in the translated genomes. The matches to five strongest prokaryotic prototypes and three prototypes of C. elegans are identified in the sequences of crystallized proteins, and their structures analyzed. Corresponding segments of the polypeptide chains in majority of cases form closed loops, though evolutionary fate of every prototype element is shown to be rather diverse. Then loop ends can be separated by a sequence-wise distant segments and stabilized by the spatial interactions in the context of the overall globular structure. The units belong to a presumably limited spectrum of the sequence prototypes, full repertoire of which would constitute a proteomic code.     

An analysis of side-chain conformation in proteins.

The crystal structures of a number of globular proteins are currently available. An analysis of the distribution of side-chains among different allowed conformations in these proteins has been carried out. The observed conformations of individual residues are discussed on the basis of well-known stereochemical criteria. The population distribution of side-chains in different allowed regions in conformational space can be explained largely on the basis of simple steric considerations. In addition to examining the conformational behaviour of individual residues, some population distributions of conformational angles of general interest involving groups of residues have also been analyzed.