Quantifying the effect of burial of amino acid residues on protein stability

The average contribution of individual residue to folding stability and its dependence on buried accessible surface area (ASA) are obtained by two different approaches. One is based on experimental mutation data, and the other uses a new knowledge-based atom-atom potential of mean force. We show that the contribution of a residue has a significant correlation with buried ASA and the regression slopes of 20 amino acid residues (called the buriability) are all positive (pro-burial). The buriability parameter provides a quantitative measure of the driving force for the burial of a residue. The large buriability gap observed between hydrophobic and hydrophilic residues is responsible for the burial of hydrophobic residues in soluble proteins. Possible factors that contribute to the buriability gap are discussed.     

Prediction of protein accessible surface areas by support vector regression

A novel support vector regression (SVR) approach is proposed to predict protein accessible surface areas (ASAs) from their primary structures. In this work, we predict the real values of ASA in squared angstroms for residues instead of relative solvent accessibility. Based on protein residues, the mean and median absolute errors are 26.0 A(2) and 18.87 A(2), respectively. The correlation coefficient between the predicted and observed ASAs is 0.66. Cysteine is the best predicted amino acid (mean absolute error is 13.8 A(2) and median absolute error is 8.37 A(2)), while arginine is the least predicted amino acid (mean absolute error is 42.7 A(2) and median absolute error is 36.31 A(2)). Our work suggests that the SVR approach can be directly applied to the ASA prediction where data preclassification has been used.     

Conformational behavior of ionic self-complementary peptides

Several de novo designed ionic peptides that are able to undergo conformational change under the influence of temperature and pH were studied. These peptides have two distinct surfaces with regular repeats of alternating hydrophilic and hydrophobic side chains. This permits extensive ionic and hydrophobic interactions resulting in the formation of stable beta-sheet assemblies. The other defining characteristic of this type of peptide is a cluster of negatively charged aspartic or glutamic acid residues located toward the N-terminus and positively charged arginine or lysine residues located toward the C-terminus. This arrangement of charge balances the alpha-helical dipole moment (C --> N), resulting in a strong tendency to form stable alpha-helices as well. Therefore, these peptides can form both stable alpha-helices and beta-sheets. They are also able to undergo abrupt structural transformations between these structures induced by temperature and pH changes. The amino acid sequence of these peptides permits both stable beta-sheet and alpha-helix formation, resulting in a balance between these two forms as governed by the environment. Some segments in proteins may also undergo conformational changes in response to environmental changes. Analyzing the plasticity and dynamics of this type of peptide may provide insight into amyloid formation. Since these peptides have dynamic secondary structure, they will serve to refine our general understanding of protein structure.     

IR study of self-assembly of capsular exopolymers from Pseudomonas sp. NCIMB 2021 on hydrophilic and hydrophobic surfaces

Capsular exopolymers (EPS) of the bacterium Pseudomonas sp. NCIMB 2021 are allowed to self-assemble on hydrophilic and hydrophobic gold surfaces. Tapping mode atomic force microscopy confirms the differences in the surface topography between EPS adsorbed on both surfaces. Fourier-transform IR spectroscopy indicates that the EPS surface coverage is much greater on the hydrophobic surface. Furthermore, an increased contribution is observed from hydrophobic (i.e., methyl and tyrosyl residues) and electrostatic (i.e., carboxylate residues) groups at the hydrophobic surface, but there is relatively less neutral polymer compared to the hydrophilic surface. The behavior of this EPS is in agreement with the behavior of cells of Pseudomonas sp. NCIMB 2021 at hydrophilic and hydrophobic surfaces.     

Look-up tables for protein solvent accessibility prediction and nearest neighbor effect analysis

We developed dictionaries of two-, three-, and five-residue patterns in proteins and computed the average solvent accessibility of the central residues in their native proteins. These dictionaries serve as a look-up table for making subsequent predictions of solvent accessibility of amino acid residues. We find that predictions made in this way are very close to those made using more sophisticated methods of solvent accessibility prediction. We also analyzed the effect of immediate neighbors on the solvent accessibility of residues. This helps us in understanding how the same residue type may have different accessible surface areas in different proteins and in different positions of the same protein. We observe that certain residues have a tendency to increase or decrease the solvent accessibility of their neighboring residues in C- or N-terminal positions. Interestingly, the C-terminal and N-terminal neighbor residues are found to have asymmetric roles in modifying solvent accessibility of residues. As expected, similar neighbors enhance the hydrophobic or hydrophilic character of residues. Detailed look-up tables are provided on the web at www.netasa.org/look-up/.     

Surface, subunit interfaces and interior of oligomeric proteins

The solvent-accessible surface area (As) of 23 oligomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. As is correlated with the protein molecular weight, and a power law predicts its value to within 5% on average. The accessible surface of the average oligomer is similar to that of monomeric proteins in its hydropathy and amino acid composition. The distribution of the 20 amino acid types between the protein surface and its interior is also the same as in monomers. Interfaces, i.e. surfaces involved in subunit contacts, differ from the rest of the subunit surface. They are enriched in hydrophobic side-chains, yet they contain a number of charged groups, especially from Arg residues, which are the most abundant residues at interfaces except for Leu. Buried Arg residues are involved in H-bonds between subunits. We counted H-bonds at interfaces and found that several have none, others have one H-bond per 200 A2 of interface area on average (1 A = 0.1 nm). A majority of interface H-bonds involve charged donor or acceptor groups, which should make their contribution to the free energy of dissociation significant, even when they are few. The smaller interfaces cover about 700 A2 of the subunit surface. The larger ones cover 3000 to 10,000 A2, up to 40% of the subunit surface area in catalase. The lower value corresponds to an estimate of the accessible surface area loss required for stabilizing subunit association through the hydrophobic effect alone. Oligomers with small interfaces have globular subunits with accessible surface areas similar to those of monomeric proteins. We suggest that these oligomers assemble from preformed monomers with little change in conformation. In oligomers with large interfaces, isolated subunits should be unstable given their excessively large accessible surface, and assembly is expected to require major structural changes.     

Molecular dynamics simulations of a protein on hydrophobic and hydrophilic surfaces.

Molecular dynamics simulations have been used to investigate the behavior of the peripheral membrane protein, cytochrome c, covalently tethered to hydrophobic (methyl-terminated) and hydrophilic (thiol-terminated) self-assembled monolayers (SAMs). The simulations predict that the protein will undergo minor structural changes when it is tethered to either surface, and the structures differ qualitatively on the two surfaces: the protein is less spherical on the hydrophilic SAM where the polar surface residues reach out to interact with the SAM surface. The protein is completely excluded from the hydrophobic SAM but partially dissolves in the hydrophilic SAM. Consequently, the surface of the thiol-terminated SAM is considerably less ordered than that of the methyl-terminated SAM, although a comparable, high degree of order is maintained in the bulk of both SAMs: the chains exhibit collective tilts in the nearest-neighbor direction at angles of 20 degrees and 17 degrees with respect to the surface normal in the hydrophobic and the hydrophilic SAMs, respectively. On the hydrophobic SAM the protein is oriented so that the heme plane is more nearly parallel to the surface, whereas on the hydrophilic surface it is more nearly perpendicular. The secondary structure of the protein, dominated by alpha helices, is not significantly affected, but the structure of the loops as well as the helix packing is slightly modified by the surfaces.     

A novel method reveals that solvent water favors polyproline II over beta-strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA)

In aqueous solution, the ensemble of conformations sampled by peptides and unfolded proteins is largely determined by their interaction with water. It has been a long-standing goal to capture these solute-water energetics accurately and efficiently in calculations. Historically, accessible surface area (ASA) has been used to estimate these energies, but this method breaks down when applied to amphipathic peptides and proteins. Here we introduce a novel method in which hydrophobic ASA is determined after first positioning water oxygens in hydrogen-bonded orientations proximate to all accessible peptide/protein backbone N and O atoms. This conditional hydrophobic accessible surface area is termed CHASA. The CHASA method was validated by predicting the polyproline-II (P(II)) and beta-strand conformational preferences of non-proline residues in the coil library (i.e., non-alpha-helix, non-beta-strand, non-beta-turn library derived from X-ray elucidated structures). Further, the method successfully rationalizes the previously unexplained solvation energies in polyalanyl peptides and compares favorably with published experimentally determined P(II) residue propensities. We dedicate this paper to Frederic M. Richards.     

An examination of the conservation of surface patch polarity for proteins

MOTIVATION: The solubility of a protein is crucial for its function and is therefore an evolutionary constraint. As the solubility of a protein is related to the distribution of polar and hydrophobic residues on its solvent accessible surface, such a constraint should provide a valuable insight into the evolution of protein surfaces. We examine how the surfaces of proteins have evolved by considering how the average hydrophobicities of patches of surface residues vary across homologous proteins. We derive distributions for the average hydrophobicity/philicity of surface patches at a residue-based level-which we refer to as the residue hydrophobic density. This is computed for a set of 28 monomeric proteins and their homologues. The resulting distributions are compared with a set of randomized sequences, with the same residue content. RESULTS: We find that the patches, involving typically more than 10 residues, maintain a more hydrophilic surface than one would expect from a random substitution model, indicating a cooperative behaviour for these surfaces residues in terms of this single variable. SUPPLEMENTARY INFORMATION: Additional plots for all of the proteins examined in this paper can be found at: http://www.ebi.ac.uk/~shanahan/PCon/index.html     

Solvent accessibility studies on polysaccharides

The solvent accessible surface area (ASA) of the polysaccharides, namely (i) carrageenan (1CAR); (ii) agarose (1AGA); (iii) guaran (GUR); (iv) capsular polysaccharide (1CAP); and (v) hyaluronan (1HUA), have been computed using the solvent accessibility technique of Lee and Richards. The results show that the average variation of ASA for the various atoms in the molecules lie in the range 1-30 A(2). Irrespective of position of sulfation, either at two or four in the sugar residues in 1CAR, the charged groups interact almost equally with the solvent. The ASA values for the chains A and B in 1AGA and 1CAR indicate that there are not much interchain interactions and the chains in both the molecules interact equally with the solvent. Residue-wise analysis indicates that the ASAs of residues vary alternately, high-low-high value pattern that is similar to that of the hydrophobic behaviour of beta-strands in proteins. The results also suggest that in these polysaccharides D-configuration residues have higher ASA than L-configuration residues.     

Adsorption of zein on surfaces with controlled wettability and thermal stability of adsorbed zein films

Adsorption characteristics of zein protein on hydrophobic and hydrophilic surfaces have been investigated to understand the orientation changes associated with the protein structure on a surface. The protein is adsorbed by a self-assembly procedure on a monolayer-modified gold surface. It is observed that zein shows higher affinity toward hydrophilic than hydrophobic surfaces on the basis of the initial adsorption rate followed by quartz crystal microbalance studies. Reflection absorption infrared (RAIR) spectroscopic studies reveal the orientation changes associated with the adsorbed zein films. Upon adsorption, the protein is found to be denatured and the transformation of alpha-helix to beta-sheet form is inferred. This transformation is pronounced when the protein is adsorbed on hydrophobic surfaces as compared to hydrophilic surfaces. Electrochemical techniques (cyclic voltammetry and impedance techniques) are very useful in assessing the permeability of zein film. It is observed that the zein moieties adsorbed on hydrophilic surfaces are highly impermeable in nature and act as a barrier for small molecules. The topographical features of the deposits before and after adsorption are analyzed by atomic force microscopy. The protein adsorbed on hydrophilic surface shows rod- and disclike features that are likely to be the base units for the growth of cylindrical structures of zein. The thermal stability of the adsorbed zein film has been followed by variable-temperature RAIR measurements.     

Surface beta-strands in proteins: identification using an hydropathy technique

From a representative set of monomeric globular proteins with known three-dimensional structures, beta-strands with lengths > or = 5 amino acids have been identified and catalogued. By ascertaining the accessible surface areas of the constituent residues in these strands, and by checking whether the exposed/buried pattern is 80% or more similar to that in an idealized surface strand, a subset of structures can be delineated in which the beta-strands are all sited on the surface of the protein. The corresponding sequence data show that about 50% of the residues are apolar (Val, Ile, Leu, Phe, Tyr, Ala) and that the common occurrence of valine (14.3%), isoleucine (9.6%), and threonine (8.1%) is a characteristic feature. The frequencies of occurrence of those amino acids in the strands that face the aqueous environment and the interior have also been determined separately and show that most surface strands have a substructure of the form (apolar-X)(n), where X is approximately equally divided between apolar, charged, and hydrophilic residues. Using the frequency data thus obtained, allied with an algorithm to delineate potential surface beta-strands from characteristic hydropathy profiles, it is now possible to search through the sequences of proteins with unknown tertiary structures and make realistic predictions of the presence of this element of structure on the protein surface. In addition, new data are presented on the distribution of the various types of residues on the surface of proteins and in their interior. Significant differences were observed, not all of which have been identified previously. Furthermore, the distribution of the types of residue in a surface beta-strand was compared to that corresponding to the surfaces of all of the proteins in our database. Again, very characteristic differences were observed. These are helpful in recognizing the presence of surface beta-strands.     

Surface-accessible residues in the monomeric and assembled forms of a bacterial surface layer protein

The S-layer protein SbsB of the thermophilic, Gram-positive organism Bacillus stearothermophilus PV72/p2 forms a crystalline, porous array constituting the outermost component of the cell envelope. SbsB has a molecular mass of 98 kDa, and the corresponding S-layer exhibits an oblique lattice symmetry. To investigate the molecular structure and assembly of SbsB, we replaced 75 residues (mainly serine, threonine, and alanine), located throughout the primary sequence, with cysteine, which is not found in the wild-type protein. As determined by electron microscopy, 72 out of 75 mutants formed regularly-structured self-assembly products identical to wild-type, thereby proving that the replacement of most of the selected amino acids by cysteine does not dramatically alter the structure of the protein. The three defective mutants, which showed a greatly reduced ability to self-assemble, were, however, successfully incorporated into S-layers of wild-type protein. Monomeric SbsB mutants and SbsB mutants assembled into S-layers were subjected to a surface accessibility screen by targeted chemical modification with a 5-kDa hydrophilic cysteine-reactive polyethylene glycol conjugate. In the monomeric form of SbsB, 34 of the examined residues were not surface accessible, while 23 were classified as very accessible, and 18 were of intermediate surface accessibility. By contrast, in the assembled S-layers, 57 of the mutated residues were not accessible, six were very accessible, and 12 of intermediate accessibility. Together with other structural information, the results suggest a model for SbsB in which functional domains are segregated along the length of the polypeptide chain.     

Conformational transitions of human serum albumin depending on pH. Study using tritium markers]

Accessible surfaces of the HSA molecule in N-, F- and B-forms were studied in the present work by tritium labelling method which allowed to obtain detailed information on N-F- and N-B-transitions. In was shown that the F-form in comparison top the N-form is characterized by more high accessibility of Ser, Ala, Ile, Tyr, Phe, His, Arg, Pro, Val and Phe residues and in the B-form Tyr, Ser, Arg, Gly, Ile, Phe and Pro residues turn to be highly accessible. Full accessible surfaces of protein molecule at N-F- and N-B-transitions increase respectively from 39,000 to 70,400 A2 and from 39,000 to 47,000 A2. Basing on the prevailing increase of hydrophobic residues accessibility it is supposed that the molecule expansion testifies the separation of the subunits forming the molecule.     

Structural background of cyclodextrin-protein interactions

Cyclodextrins are cyclic oligosaccharides with the shape of a hollow truncated cone. Their exterior is hydrophilic and their cavity is hydrophobic, which gives cyclodextrins the ability to accommodate hydrophobic molecules/moieties in the cavity. This special molecular arrangement accounts for the variety of beneficial effects cyclodextrins have on proteins, which is widely used in pharmacological applications. We have studied the interaction between beta-cyclodextrin and four non-carbohydrate-binding model proteins: ubiquitin, chymotrypsin inhibitor 2 (CI2), S6 and insulin SerB9Asp by NMR spectroscopy at varying structural detail. We demonstrate that the interaction of beta-cyclodextrin and our model proteins takes place at specific sites on the protein surface, and that solvent accessibility of those sites is a necessary but not compelling condition for the occurrence of an interaction. If this behaviour can be generalized, it might explain the wide range of different effects of cyclodextrins on different proteins: aggregation suppression (if residues responsible for aggregation are highly solvent accessible), protection against degradation (if point of attack of a protease is sterically 'masked' by cyclodextrin), alteration of function (if residues involved in function are 'masked' by cyclodextrin). The exact effect of cyclodextrins on a given protein will always be related to the particular structure of this protein.     

Solvent accessibility in native and isolated domain environments: general features and implications to interface predictability

A non-redundant database of 4536 structural domains, comprising more than 790,000 residues, has been used for the calculation of their solvent accessibility in the native protein environment and then in the isolated domain environment. Nearly 140,000 (18%) residues showed a change in accessible surface area in the above two conditions. General features of this change under these two circumstances have been pointed out. Propensities of these interfacing amino acid residues have been calculated and their variation for different secondary structure types has been analyzed. Actual amount of surface area lost by different secondary structures is higher in the case of helix and strands compared to coil and other conformations. Overall change in surface area in hydrophobic and uncharged residues is higher than that in charged residues. An attempt has been made to know the predictability of interface residues from sequence environments. This analysis and prediction results have significant implications towards determining interacting residues in proteins and for the prediction of protein-protein, protein-ligand, protein-DNA and similar interactions.     

Crystal structure analysis of amicyanin and apoamicyanin from Paracoccus denitrificans at 2.0 A and 1.8 A resolution.

The crystal structure of amicyanin, a cupredoxin isolated from Paracoccus denitrificans, has been determined by molecular replacement. The structure has been refined at 2.0 A resolution using energy-restrained least-squares procedures to a crystallographic residual of 15.7%. The copper-free protein, apoamicyanin, has also been refined to 1.8 A resolution with residual 15.5%. The protein is found to have a beta-sandwich topology with nine beta-strands forming two mixed beta-sheets. The secondary structure is very similar to that observed in the other classes of cupredoxins, such as plastocyanin and azurin. Amicyanin has approximately 20 residues at the N-terminus that have no equivalents in the other proteins; a portion of these residues forms the first beta-strand of the structure. The copper atom is located in a pocket between the beta-sheets and is found to have four coordinating ligands: two histidine nitrogens, one cysteine sulfur, and, at a longer distance, one methionine sulfur. The geometry of the copper coordination is very similar to that in the plant plastocyanins. Three of the four copper ligands are located in the loop between beta-strands eight and nine. This loop is shorter than that in the other cupredoxins, having only two residues each between the cysteine and histidine and the histidine and methionine ligands. The amicyanin and apoamicyanin structures are very similar; in particular, there is little difference in the positions of the coordinating ligands with or without copper. One of the copper ligands, a histidine, lies close to the protein surface and is surrounded on that surface by seven hydrophobic residues. This hydrophobic patch is thought to be important as an electron transfer site.