Coevolution in defining the functional specificity

Covariation between sites can arise due to a common evolutionary history. At the same time, structure and function of proteins play significant role in evolvability of different sites that are not directly connected with the common ancestry. The nature of forces which cause residues to coevolve is still not thoroughly understood, it is especially not clear how coevolutionary processes are related to functional diversification within protein families. We analyzed both functional and structural factors that might cause covariation of specificity determinants and showed that they more often participate in coevolutionary relationships with each other and other sites compared with functional sites and those sites that are not under strong functional constraints. We also found that protein sites with higher number of coevolutionary connections with other sites have a tendency to evolve slower. Our results indicate that in some cases coevolutionary connections exist between specificity sites that are located far away in space but are under similar functional constraints. Such correlated changes and compensations can be realized through the stepwise coevolutionary processes which in turn can shed light on the mechanisms of functional diversification.     

The interplay of sequence and stereochemistry in defining conformation in proteins and polypeptides

Sequential specification of conformation in proteins and polypeptides is a triangular interplay involving the system of linked peptides, the sequences in side chains, and water as solvent. Stereochemistry in side chain linkages is obviously important in the interaction between all of the players, but no specification of its explicit role, if any, in linking sequence with conformation has been made. Flory and coworkers made a puzzling observation in 1967 that, when mutated from poly-L to alternating-L,D stereochemical structure, polypeptides will suffer a reduction in overall dimension or characteristic ratio by an astonishing factor of 10 and to a value even lower than that predicted for free rotation (Miller, W. G.; Brant, D. A.; Flory, P. J. J Mol Biol 1967, 23, 67-80). Enquiring into this longstanding puzzle, Durani and coworkers found that the stereochemical modification will also abolish conformational sensitivity in polypeptide structure to solvent, because electrostatic interactions in the system of linked peptides are transformed from a condition of mutual conflict to one of harmony (Ramakrishnan, V.; Ranbhor, R.; Kumar, A.; Durani, S. J Phys Chem B 2006, 110, 9314-9323). Thus, poly-L stereochemistry could be the fulcrum linking sequences with phi,psis in protein and polypeptide structures, via dielectric arbitrations in a conflicting type of interpeptide electrostatics, in agreement with the electrostatic screening model of Avbelj and Moult (Avbelj, F.; Moult, J. Biochemistry 1995, 34, 755-764).     

Thermal denaturation of oligomeric proteins: structural and functional modifications, the sequence of steps

On the basis of literature and own experimental data changes in the structural and functional properties of some oligomeric proteins of the blood system (hemoglobin, lactate dehydrogenase, superoxide dismutase, catalase) exposed to the influence of temperature in a broad range were analyzed. The many-phase character of the temperature modification of protein molecules with different values of functional and kinetic parameters for each of revealed stages was discovered. At a critical temperature and at higher values, the dissociation of oligomeric proteins into separate subunits was shown to occur along with the typical "loosening" of the protein globule. It was shown that low-molecular components (subunits) can subsequently associate with one another and with oligometic proteins, which leads to irreversible denaturation and to the unusual physicochemical behavior of protein molecules. Schemes of processes underlying the temperature modifications of the proteins studied were elaborated.     

Finding functional sites in structural genomics proteins

Assigning function to structures is an important aspect of structural genomics projects, since they frequently provide structures for uncharacterized proteins. Similarities uncovered by structure alignment can suggest a similar function, even in the absence of sequence similarity. For proteins adopting novel folds or those with many functions, this strategy can fail, but functional clues can still come from comparison of local functional sites involving a few key residues. Here we assess the general applicability of functional site comparison through the study of 157 proteins solved by structural genomics initiatives. For 17, the method bolsters confidence in predictions made based on overall fold similarity. For another 12 with new folds, it suggests functions, including a putative phosphotyrosine binding site in the Archaeal protein Mth1187 and an active site for a ribose isomerase. The approach is applied weekly to all new structures, providing a resource for those interested in using structure to infer function.     

Metalloproteomics: high-throughput structural and functional annotation of proteins in structural genomics

A high-throughput method for measuring transition metal content based on quantitation of X-ray fluorescence signals was used to analyze 654 proteins selected as targets by the New York Structural GenomiX Research Consortium. Over 10% showed the presence of transition metal atoms in stoichiometric amounts; these totals as well as the abundance distribution are similar to those of the Protein Data Bank. Bioinformatics analysis of the identified metalloproteins in most cases supported the metalloprotein annotation; identification of the conserved metal binding motif was also shown to be useful in verifying structural models of the proteins. Metalloproteomics provides a rapid structural and functional annotation for these sequences and is shown to be approximately 95% accurate in predicting the presence or absence of stoichiometric metal content. The project's goal is to assay at least 1 member from each Pfam family; approximately 500 Pfam families have been characterized with respect to transition metal content so far.     

Hidden Markov model-derived structural alphabet for proteins: the learning of protein local shapes captures sequence specificity

Understanding and predicting protein structures depend on the complexity and the accuracy of the models used to represent them. We have recently set up a Hidden Markov Model to optimally compress protein three-dimensional conformations into a one-dimensional series of letters of a structural alphabet. Such a model learns simultaneously the shape of representative structural letters describing the local conformation and the logic of their connections, i.e. the transition matrix between the letters. Here, we move one step further and report some evidence that such a model of protein local architecture also captures some accurate amino acid features. All the letters have specific and distinct amino acid distributions. Moreover, we show that words of amino acids can have significant propensities for some letters. Perspectives point towards the prediction of the series of letters describing the structure of a protein from its amino acid sequence.     

Plasmodium falciparum GPI mannosyltransferase-III has novel signature sequence and is functional

The glycosylphosphatidylinositol (GPI) anchors of Plasmodium falciparum are indispensable for parasite survival since merozoite surface proteins-1, -2, -4, -5, and -10, crucial for erythrocyte invasion, are GPI-anchored. Therefore, the GPI biosynthetic pathway can offer potential targets for novel anti-malarial drugs. Here, we characterized the putative P. falciparum PIG-B gene (PfPIGB) that encodes mannosyltransferase-III of GPI biosynthesis. PfPIGB mRNA is transcribed in a developmental stage specific manner. A protein corresponding to the expected size of PfPIG-B is expressed by the parasite and is localized in the endoplasmic reticulum. Treatment of parasites with PfPIG-B specific siRNA caused reduction in GPI synthesis, affecting the PIG-B specific GPI intermediate. These data demonstrate that PfPIG-B is functional and encodes mannosyltransferase-III of the parasite GPI biosynthesis. The parasite PfPIG-B is novel in that its signature sequence HKEHKI is unique and is only partially conserved as compared to HKEXRF signature motif of mammalian PIG-B enzymes.     

Role of the charge-charge interactions in defining stability and halophilicity of the CspB proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.     

Exposed loop domains of complexed 14-3-3 proteins contribute to structural diversity and functional specificity

The 14-3-3 family of proteins functions through protein:phosphoprotein interactions, the nature of which has been elucidated using x-ray crystallography. However, some key structural features in nonconserved regions have yet to be fully resolved, leaving open questions regarding the functional selectivity of 14-3-3 family members for diverse clients. In an effort to study surface accessible structural features in 14-3-3 containing macromolecular complexes and to illuminate important structure/function variations among the 14-3-3 isoforms, we determined the epitopes for three unique monoclonal antibodies (mAbs) developed against the Arabidopsis (Arabidopsis thaliana) G-box DNA:protein complex. The epitopes mapped to different loops in a phylogenetically important subset of the 13 14-3-3 family members. All three epitopes were on a common exposed face of complexed 14-3-3s. Two of the mAbs recognized linear sequences within loops 5 and 6, while the third mAb recognized 14-3-3 residues surrounding the pivotal medial Gly in the divalent cation-binding domain of loop 8, together with distal residue(s) in the putative dynamic 10th helix that has yet to be determined by crystallography. Gly at this loop 8 position is unique to nonepsilon 14-3-3 isoforms of the plant kingdom, suggesting that this region constitutes a plant-specific key functional 14-3-3 feature and highlighting that the loop 8 region is functionally significant. Mutagenesis of the medial amino acid in the loop 8 domain changed the flexibility of the C terminus and altered client peptide-binding selectivity, demonstrating the functional significance of the surface accessible, evolutionarily distinct loop 8 domain.     

Achieving stability and conformational specificity in designed proteins via binary patterning

We have developed a method to determine the optimal binary pattern (arrangement of hydrophobic and polar amino acids) of a target protein fold prior to amino acid sequence selection in protein design studies. A solvent accessible surface is generated for a target fold using its backbone coordinates and "generic" side-chains, which are constructs whose size and shape are similar to an average amino acid. Each position is classified as hydrophobic or polar according to the solvent exposure of its generic side-chain. The method was tested by analyzing a set of proteins in the Protein Data Bank and by experimentally constructing and analyzing a set of engrailed homeodomain variants whose binary patterns were systematically varied. Selection of the optimal binary pattern results in a designed protein that is monomeric, well-folded, and hyperthermophilic. Homeodomain variants with fewer hydrophobic residues are destabilized, while additional hydrophobic residues induce aggregation. Binary patterning, in conjunction with a force field that models folded state energies, appears sufficient to satisfy two basic goals of protein design: stability and conformational specificity.     

Understanding relationship between sequence and functional evolution in yeast proteins

The underlying relationship between functional variables and sequence evolutionary rates is often assessed by partial correlation analysis. However, this strategy is impeded by the difficulty of conducting meaningful statistical analysis using noisy biological data. A recent study suggested that the partial correlation analysis is misleading when data is noisy and that the principal component regression analysis is a better tool to analyze biological data. In this paper, we evaluate how these two statistical tools (partial correlation and principal component regression) perform when data are noisy. Contrary to the earlier conclusion, we found that these two tools perform comparably in most cases. Furthermore, when there is more than one ‘true’ independent variable, partial correlation analysis delivers a better representation of the data. Employing both tools may provide a more complete and complementary representation of the real data. In this light, and with new analyses, we suggest that protein length and gene dispensability play significant, independent roles in yeast protein evolution. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Ligand specificity and conformational stability of human fatty acid-binding proteins.

Fatty acid binding proteins (FABPs) are small cytosolic proteins with virtually identical backbone structures that facilitate the solubility and intracellular transport of fatty acids. At least eight different types of FABP occur, each with a specific tissue distribution and possibly with a distinct function. To define the functional characteristics of all eight human FABPs, viz. heart (H), brain (B), myelin (M), adipocyte (A), epidermal (E), intestinal (I), liver (L) and ileal lipid-binding protein (I-LBP), we studied their ligand specificity, their conformational stability and their immunological crossreactivity. Additionally, binding of bile acids to I-LBP was studied. The FABP types showed differences in fatty acid binding affinity. Generally, the affinity for palmitic acid was lower than for oleic and arachidonic acid. All FABP types, except E-FABP, I-FABP and I-LBP interacted with 1-anilinonaphtalene-8-sulphonic acid (ANS). Only L-FABP, I-FABP and M-FABP showed binding of 11-((5-dimethylaminonaphtalene-1-sulfonyl)amino)undecanoic acid (DAUDA). I-LBP showed increasing binding of bile acids in the order taurine-conjugated>glycine-conjugated>unconjugated bile acids. A hydroxylgroup of bile acids at position 7 decreased and at position 12 increased the binding affinity to I-LBP. The fatty acid-binding affinity and the conformation of FABP types were differentially affected in the presence of urea. Our results demonstrate significant differences in ligand binding, conformational stability and surface properties between different FABP types which may point to a specific function in certain cells and tissues. The preference of I-LBP (but not L-FABP) for conjugated bile acids is in accordance with a specific role in bile acid reabsorption in the ileum.     

Structure-based functional discovery of proteins: structural proteomics

The discovery of biochemical and cellular functions of unannotated gene products begins with a database search of proteins with structure/sequence homologues based on known genes. Very recently, a number of frontier groups in structural biology proposed a new paradigm to predict biological functions of an unknown protein on the basis of its three-dimensional structure on a genomic scale. Structural proteomics (genomics), a research area for structure-based functional discovery, aims to complete the protein-folding universe of all gene products in a cell. It would lead us to a complete understanding of a living organism from protein structure. Two major complementary experimental techniques, X-ray crystallography and NMR spectroscopy, combined with recently developed high throughput methods have played a central role in structural proteomics research; however, an integration of these methodologies together with comparative modeling and electron microscopy would speed up the goal for completing a full dictionary of protein folding space in the near future.     

Interactions between heterologous helix-loop-helix proteins generate complexes that bind specifically to a common DNA sequence

A DNA binding and dimerization motif, with apparent amphipathic helices (the HLH motif), has recently been identified in various proteins, including two that bind to immunoglobulin enhancers (E12 and E47). We show here that various HLH proteins can bind as apparent heterodimers to a single DNA motif and also, albeit usually more weakly, as apparent homodimers. The HLH domain can mediate heterodimer formation between either daughterless, E12, or E47 (Class A) and achaete-scute T3 or MyoD (Class B) to form proteins with high affinity for the kappa E2 site in the immunoglobulin kappa chain enhancer. The achaete-scute T3 and MyoD proteins do not form kappa E2-binding heterodimers together, and no active complex with N-myc was evident. The formation of a heterodimer between the daughterless and achaete-scute T3 products may explain the similar phenotypes of mutants at these two loci and the genetic interactions between them. A role of E12 and E47 in mammalian development, analogous to that of daughterless in Drosophila, is likely.     

Correlation between the structural stability and aggregation propensity of proteins

Protein aggregation, being an outcome of improper protein folding, is largely dependent on the folding kinetics of a protein. Previous studies have reported a positive correlation between the stability of the secondary structural elements of a protein and their rate of folding/unfolding. In this in silico study, the secondary and tertiary structures of proteins a) that form inclusion bodies on overexpression in Escherichia coli, b) that form amyloid fibrils and c) that are soluble on overexpression in E. coli are analyzed for certain features that are known to be associated with structural stability. The study revealed that the soluble proteins seem to have a higher rate of folding (based on contact order) and a lower percentage of exposed hydrophobic residues as compared to the inclusion body forming or amyloidogenic proteins. The soluble proteins also seem to have a more favored helix and strand composition (based on the known secondary structural propensities of amino acids). The secondary structure analyses also reveal that the evolutionary pressure is directed against protein aggregation. This understanding of the positive correlation between structural stability and solubility, along with the other parameters known to influence aggregation, could be exploited in the design of mutations aimed at reducing the aggregation propensity of the proteins.     

Human hair keratin-associated proteins: sequence regularities and structural implications

In this paper, we undertake a sequence analysis of the human keratin-associated proteins (KAP). This analysis has revealed two fundamental pentapeptide quasi-repeats (A and B) of the form C-C-X-P-X and C-C-X-S/T-S/T, respectively. The A repeats are also commonly found in two subforms A1 and A2, -C-C-Q-P-X and C-C-R-P-X, respectively-similar to those found in sheep wool 30-40 years previously. Some high-sulphur and ultra-high sulphur proteins contain predominantly A repeats or B repeats but not regular combinations of them, whereas others are characterised by a contiguous pair of pentapeptide repeats that largely (though imperfectly) alternate to generate decapeptide motifs of the form AB, A1B or A2B. The A and B repeats sometimes occur in complex runs and can generate both 19- and 20-residue repeats of the form BABB' or BA1AA, respectively, where the prime indicates a motif truncated by one residue. Likewise, a 42-residue repeat with BA1BXAAAB (40 residues) separated by a di-serine (two residues) has been observed in an ultra-high sulphur protein from cuticle. To understand the possible conformations adopted by the A and B motifs, a search was initiated of the PDB structural database for a number of overlapping pentapeptide repeats. The total number of matches was 658 and these were found in 451 different proteins. From representative and unique structures the means and standard deviations were calculated for the Phi(i) and Psi(i) angles for the C-C-X-P-X and the C-C-X-S/T-S/T motifs. Molecular modelling has been employed to represent the "average" structure found from crystallographic and nmr data determined for each motif in other proteins. The conformation of consecutive A repeats with proline residues in the cis state is akin to a string of disulphide bond-stabilised pentapeptide knots between which there is relative freedom of rotation about the single bonds that link them. For B pentapeptides, however, the likelihood that a similar disulphide bond is formed appears much lower. This may give additional conformational flexibility to the chain and hence allow the A pentapeptides greater opportunity to interact appropriately with the IF via disulphide bonds, ionic interactions and/or hydrogen bonding.