Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.     

Interior and surface of monomeric proteins

The solvent-accessible surface area (As) of 46 monomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. The As of these proteins can be determined to within 1 to 2% and that of their individual residues to within 10 to 20%. The As values of proteins are correlated with their molecular weight (Mr) in the range 4000 to 35,000: the power law As = 6.3 M0.73 predicts protein As values to within 4% on average. The average water-accessible surface is found to be 57% non-polar, 24% polar and 19% charged, with 5% root-mean-square variations. The molecular surface buried inside the protein is 58% non-polar, 39% polar and 4% charged. The buried surface contains more uncharged polar groups (mostly peptides) than the surface that remains accessible, but many fewer charged groups. On average, 15% of residues in small proteins and 32% in larger ones may be classed as "buried residues", having less than 5% of their surface accessible to the solvent. The accessibilities of most other residues are evenly distributed in the range 5 to 50%. Although the fraction of buried residues increases with molecular weight, the amino acid compositions of the protein interior and surface show no systematic variation with molecular weight, except for small proteins that are often very rich in buried cysteines. From amino acid compositions of protein surfaces and interiors we calculate an effective coefficient of partition for each type of residue, and derive an implied set of transfer free energy values. This is compared with other sets of partition coefficients derived directly from experimental data. The extent to which groups of residues (charged, polar and non-polar) are buried within proteins correlates well with their hydrophobicity derived from amino acid transfer experiments. Within these three groups, the correlation is low.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

Amino acid side-chain partition energies and distribution of residues in soluble proteins

Energies required to transfer amino acid side chains from water to less polar environments were calculated from results of several studies and compared with several statistical analyses of residue distributions in soluble proteins. An analysis that divides proteins into layers parallel with their surfaces is more informative than those that simply classify residues as exposed or buried. Most residues appear to be distributed as a function of the distance from the protein-water interface in a manner consistent with partition energies calculated from partitioning of amino acids between water and octanol phases and from solubilities of amino acids in water, ethanol, and methanol. Lys, Arg, Tyr, and Trp residues tend to concentrate near the water-protein interface where their apolar side-chain components are more buried than their polar side-chain components. Residue distributions calculated in this manner do not correlate well with side-chain solvation energies calculated from vapor pressures of side-chain analogs over a water phase. Results of statistical studies that classify residues as exposed to solvent or buried inside the protein interior appear to depend on the method used to classify residues. Data from some of these studies correlate better with solvation energies, but other data correlate better with partition energies. Most other statistical methods that have been used to evaluate effects of water on residue distributions yield results that correlate better with partition energies than with solvation energies.     

Prediction of catalytic residues in enzymes based on known tertiary structure,stability profile,and sequence conservation

The catalytic or functionally important residues of a protein are known to exist in evolutionarily constrained regions. However, the patterns of residue conservation alone are sometimes not very informative, depending on the homologous sequences available for a given query protein. Here, we present an integrated method to locate the catalytic residues in an enzyme from its sequence and structure. Mutations of functional residues usually decrease the activity, but concurrently often increase stability. Also, catalytic residues tend to occupy partially buried sites in holes or clefts on the molecular surface. After confirming these general tendencies by carrying out statistical analyses on 49 representative enzymes, these data together with amino acid conservation were evaluated. This novel method exhibited better sensitivity in the prediction accuracy than traditional methods that consider only the residue conservation. We applied it to some so-called "hypothetical" proteins, with known structures but undefined functions. The relationships among the catalytic, conserved, and destabilizing residues in enzymatic proteins are discussed.     

Thermodynamic stability of histone H3 is a necessary but not sufficient driving force for its evolutionary conservation

Determining the forces that conserve amino acid positions in proteins across species is a fundamental pursuit of molecular evolution. Evolutionary conservation is driven by either a protein's function or its thermodynamic stability. Highly conserved histone proteins offer a platform to evaluate these driving forces. While the conservation of histone H3 and H4 "tail" domains and surface residues are driven by functional importance, the driving force behind the conservation of buried histone residues has not been examined. Using a computational approach, we determined the thermodynamically preferred amino acids at each buried position in H3 and H4. In agreement with what is normally observed in proteins, we find a significant correlation between thermodynamic stability and evolutionary conservation in the buried residues in H4. In striking contrast, we find that thermodynamic stability of buried H3 residues does not correlate with evolutionary conservation. Given that these H3 residues are not post-translationally modified and only regulate H3-H3 and H3-H4 stabilizing interactions, our data imply an unknown function responsible for driving conservation of these buried H3 residues.     

Amino acid architecture and the distribution of polar atoms on the surfaces of proteins

We propose that a necessary condition for a protein to be soluble is the absence of large hydrophobic patches on its solvent-accessible surface, which can cause aggregation to occur. We note that the polar nature of the backbone of all amino acids guarantees a minimum polar content and hence can interrupt such patches. As a result, a carefully conserved detailed atomic placement of residues on the protein surface is not necessary for solubility. In order to demonstrate this, we construct a measure based on the average hydrophobicity of a simply defined patch. We use this measurement to compare surfaces that exhibit a clear difference in their solubility properties, namely, a) the solvent accessible surfaces for a set of homo-dimers and the surfaces buried in their interfaces and b) for a set of monomers the surfaces of fragments of secondary structure which are solvent accessible/inaccessible. Having demonstrated a difference in the first set of distributions, we characterize the solvent accessible surfaces of monomeric proteins. To test if cooperative behavior occurs between the atoms for these surfaces, we construct a set of randomized surfaces, which obey a very simple stereochemical constraint. We find that the observed and randomized distributions are much more similar than the previous sets we examined. This implies that while surfaces of soluble proteins must have sufficient polar content, the relative placement of atoms of one amino acid with respect to the atoms of neighboring amino acid need not be finely tuned, which provides an innate robustness for protein design and folding.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Tolerance to the substitution of buried apolar residues by charged residues in the homologous protein structures

Occurrence and accommodation of charged amino acid residues in proteins that are structurally equivalent to buried non-polar residues in homologues have been investigated. Using a dataset of 1,852 homologous pairs of crystal structures of proteins available at 2A or better resolution, 14,024 examples of apolar residues in the structurally conserved regions replaced by charged residues in homologues have been identified. Out of 2,530 cases of buried apolar residues, 1,677 of the equivalent charged residues in homologues are exposed and the rest of the charged residues are buried. These drastic substitutions are most often observed in homologous protein pairs with low sequence identity (<30%) and in large protein domains (>300 residues). Such buried charged residues in the large proteins are often located in the interface of sub-domains or in the interface of structural repeats, Beyond 7A of residue depth of buried apolar residues, or less than 4% of solvent accessibility, almost all the substituting charged residues are buried. It is also observed that acidic sidechains have higher preference to get buried than the positively charged residues. There is a preference for buried charged residues to get accommodated in the interior by forming hydrogen bonds with another sidechain than the main chain. The sidechains interacting with a buried charged residue are most often located in the structurally conserved regions of the alignment. About 50% of the observations involving hydrogen bond between buried charged sidechain and another sidechain correspond to salt bridges. Among the buried charged residues interacting with the main chain, positively charged sidechains form hydrogen bonds commonly with main chain carbonyls while the negatively charged residues are accommodated by hydrogen bonding with the main chain amides. These carbonyls and amides are usually located in the loops that are structurally variable among homologous proteins.     

Chemical characterization of lumazine protein from Photobacterium leiognathi: comparison with lumazine protein from Photobacterium phosphoreum

The properties of lumazine proteins purified from the marine bioluminescent bacteria Photobacterium phosphoreum, a psychrophile, and Photobacterium leiognathi, a relatively thermophilic species, are compared. An accurate 1:1 stoichiometry of binding of the ligand 6,7-dimethyl-8-ribityllumazine to each lumazine protein is established by back-titration of the apoprotein with the authentic ligand, using both fluorescence and absorption measurements. Neither protein contains metal cofactors, organic phosphorus, or carbohydrate. Both proteins are anionic and hydrophilic. They each contain a single Trp residue and have blocked amino terminals but otherwise differ in amino acid composition and other properties (P. phosphoreum and P. leiognathi, respectively): Met (internal), 1, 2; Cys, 2, 1; Arg, 4, 7; pI, 4.78 and 4.83, 4.38 and 4.45; Mr, 19 750, 21 300. In the P. phosphoreum protein both Cys residues are accessible, but in the P. leiognathi protein the single Cys is "buried". Modification of this buried Cys and at least one Cys in the P. phosphoreum protein prevents binding of the ligand. The UV and visible absorption spectra of both lumazine proteins denatured in 6 M guanidine hydrochloride can be accurately modeled by using the number of equivalents of the lumazine derivative and blocked aromatic amino acid model compounds determined by chemical and spectrophotometric analyses for Trp, Tyr, and Phe.     

Analysis and functional prediction of reactive cysteine residues

Cys is much different from other common amino acids in proteins. Being one of the least abundant residues, Cys is often observed in functional sites in proteins. This residue is reactive, polarizable, and redox-active; has high affinity for metals; and is particularly responsive to the local environment. A better understanding of the basic properties of Cys is essential for interpretation of high-throughput data sets and for prediction and classification of functional Cys residues. We provide an overview of approaches used to study Cys residues, from methods for investigation of their basic properties, such as exposure and pK(a), to algorithms for functional prediction of different types of Cys in proteins.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

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.     

Conformational and functional analysis of the lipid binding protein Ag-NPA-1 from the parasitic nematode Ascaridia galli

Ag-NPA-1 (AgFABP), a 15 kDa lipid binding protein (LBP) from Ascaridia galli, is a member of the nematode polyprotein allergen/antigen (NPA) family. Spectroscopic analysis shows that Ag-NPA-1 is a highly ordered, alpha-helical protein and that ligand binding slightly increases the ordered secondary structure content. The conserved, single Trp residue (Trp17) and three Tyr residues determine the fluorescence properties of Ag-NPA-1. Analysis of the efficiency of the energy transfer between these chromophores shows a high degree of Tyr-Trp dipole-dipole coupling. Binding of fatty acids and retinol was accompanied by enhancement of the Trp emission, which allowed calculation of the affinity constants of the binary complexes. The distance between the single Trp of Ag-NPA-1 and the fluorescent fatty acid analogue 11-[(5-dimethylaminonaphthalene-1- sulfonyl)amino]undecanoic acid (DAUDA) from the protein binding site is 1.41 nm as estimated by fluorescence resonance energy transfer. A chemical modification of the Cys residues of Ag-NPA-1 (Cys66 and Cys122) with the thiol reactive probes 5-({[(2-iodoacetyl)amino]ethyl}amino) naphthalene-1-sulfonic acid (IAEDANS) and N,N'-dimethyl-N-(iodoacetyl)-N'-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)ethylenediamine (IANBD), followed by MALDI-TOF analysis showed that only Cys66 was labeled. The observed similar affinities for fatty acids of the modified and native Ag-NPA-1 suggest that Cys66 is not a part of the protein binding pocket but is located close to it. Ag-NPA-1 is one of the most abundant proteins in A. galli and it is distributed extracellularly mainly as shown by immunohistology and immunogold electron microscopy. This suggests that Ag-NPA-1 plays an important role in the transport of fatty acids and retinoids.     

Sequence and structure conservation in a protein core

In order to study structural aspects of sequence conservation in families of homologous proteins, we have analyzed structurally aligned sequences of 585 proteins grouped into 128 homologous families. The conservation of a residue in a family is defined as the average residue similarity in a given position of aligned sequences. The residue similarities were expressed in the form of log-odd substitution tables that take into account the environments of amino acids in three-dimensional structures. The protein core is defined as those residues that have less then 7% solvent accessibility. The density of a protein core is described in terms of atom packing, which is investigated as a criterion for residue substitution and conservation. Although there is no significant correlation between sequence conservation and average atom packing around nonpolar residues such as leucine, valine and isoleucine, a significant correlation is observed for polar residues in the protein core. This may be explained by the hydrogen bonds in which polar residues are involved; the better their protection from water access the more stable should be the structure in that position. Proteins 33:358–366, 1998. © 1998 Wiley-Liss, Inc.     

Conservation of cis prolyl bonds in proteins during evolution

In proteins and peptides, the vast majority of peptide bonds occurs in trans conformation, but a considerable fraction (about 5%) of X-Pro bonds adopts the cis conformation. Here we study the conservation of cis prolyl residues in evolutionary related proteins. We find that overall, in contrast to local, protein sequence similarity is a clear indicator for the conformation of prolyl residues. We observe that cis prolyl residues are more often conserved than trans prolyl residues, and both are more conserved than the surrounding amino acids, which show the same extent of conservation as the whole protein. The pattern of amino acid exchanges differs between cis and trans prolyl residues. Also, the cis prolyl bond is maintained in proteins with sequence identity as low as 20%. This finding emphasizes the importance of cis peptide bonds in protein structure and function.     

The structure of hepadnaviral core antigens. Identification of free thiols and determination of the disulfide bonding pattern.

A set of wild-type and mutant human, woodchuck, and duck hepatitis viral core proteins have been prepared and used to study the free thiol groups and the disulfide bonding pattern present within the core particle. Human (HBcAg) and woodchuck (WHcAg) core proteins contain 4 cysteine residues, whereas duck (DHcAg) core protein contains a single cysteine residue. Each of the cysteines of HBcAg has been eliminated, either singly or in combinations, by a two-step mutagenesis procedure. All of the proteins were shown to have very similar physical and immunochemical properties. All assemble into essentially identical core particle structures. Therefore disulfide bonds are not essential for core particle formation. No intra-chain disulfide bonds occur. Cys107 is a free thiol buried within the particle structure, whereas Cys48 is present partly as a free sulfhydryl which is exposed at the surface of the particle. Cys61 is always and Cys48 is partly involved in interchain disulfide bonds with the identical residues of another monomer, whereas Cys183 is always involved in a disulfide bond with the Cys183 of a different monomer. WHcAg has the same pattern of bonding, whereas DHcAg lacks any disulfide bonds, and the single free sulfhydryl, Cys153 which is equivalent to Cys107 of HBcAg, is buried.     

Functional analyses of the Dof domain, a zinc finger DNA-binding domain, in a pumpkin DNA-binding protein AOBP

AOBP, a DNA-binding protein in pumpkin, contains a Dof domain that is composed of 52 amino acid residues and is highly conserved in several DNA-binding proteins of higher plants. The Dof domain has a significant resemblance to Cys2/Cys2 zinc finger DNA-binding domains of steroid hormone receptors and GATA1, but has a longer putative loop where an extra Cys residue is conserved. We show that the Dof domain in AOBP functions as a zinc finger DNA-binding domain and suggest that the Cys residue uniquely conserved in the putative loop might negatively regulate the binding to DNA.     

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Summary. Thermophilic proteins show substantially higher intrinsic thermal stability than their mesophilic counterparts. Amino acid composition is believed to alter the intrinsic stability of proteins. Several investigations and mutagenesis experiment have been carried out to understand the amino acid composition for the thermostability of proteins. This review presents some generalized features of amino acid composition found in thermophilic proteins, including an increase in residue hydrophobicity, a decrease in uncharged polar residues, an increase in charged residues, an increase in aromatic residues, certain amino acid coupling patterns and amino acid preferences for thermophilic proteins. The differences of amino acids composition between thermophilic and mesophilic proteins are related to some properties of amino acids. These features provide guidelines for engineering mesophilic protein to thermophilic protein. Authors’ addresses: Yuan-Jiang Pan, Institute of Chemical Biology and Pharmaceutical Chemistry, Zhejiang University, Zhejiang University Road 38, Hangzhou 310027, China; Wei-Fen Li, Microbiology Division, College of Animal Science, Zhejiang University, Hangzhou 310029, China     

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.