Correlation of co-ordinated amino acid changes at the two-domain interface of cysteine proteases with protein stability

Using experimental evolution to probe molecular mechanisms of protein function

Directed evolution is a powerful tool for engineering protein function. The process of directed evolution involves iterative rounds of sequence diversification followed by assaying activity of variants and selection. The range of sequence variants and linked activities generated in the course of an evolution are a rich information source for investigating relationships between sequence and function. Key residue positions determining protein function, combinatorial contributors to activity and even potential functional mechanisms have been revealed in directed evolutions. The recent application of high throughput sequencing substantially increases the information that can be retrieved from directed evolution experiments. Combined with computational analysis this additional sequence information has allowed high-resolution analysis of individual residue contributions to activity. These developments promise to significantly enhance the depth of insight that experimental evolution provides into mechanisms of protein function.

Prediction of contacts from correlated sequence substitutions

Recent work has led to a substantial improvement in the accuracy of predictions of contacts between amino acids using evolutionary information derived from multiple sequence alignments. Where large numbers of diverse sequence relatives are available and can be aligned to the sequence of a protein of unknown structure it is now possible to generate high-resolution models without recourse to the structure of a template. In this review we describe these exciting new techniques and critically assess the state-of-the-art in contact prediction in the light of these. While concentrating on methods, we also discuss applications to protein and RNA structure prediction as well as potential future developments.

Prediction of protein contacts from correlated sequence substitutions

Recent work has led to a substantial improvement in the accuracy of predictions of contacts between amino acids using evolutionary information derived from multiple sequence alignments. Where large numbers of diverse sequence relatives are available and can be aligned to the sequence of a protein of unknown structure, it is now possible to generate high-resolution models without recourse to the structure of a template. In this review, we describe these exciting new techniques and critically assess the state of the art in contact prediction in light of them. We discuss areas for immediate research and development as well as potential future developments.

Using evolutionary information to find specificity-determining and co-evolving residues

Intricate networks of protein interactions rely on the ability of a protein to recognize its targets: other proteins, ligands, and sites on DNA and RNA. To recognize other molecules, it was suggested that a protein uses a small set of specificity-determining residues (SDRs). How can one find these residues in proteins and distinguish them from other functionally important amino acids? A number of bioinformatics methods to predict SDRs have been developed in recent years. These methods use genomic information and multiple sequence alignments to identify positions exhibiting a specific pattern of conservation and variability. The challenge is to delineate the evolutionary pattern of SDRs from that of the active site residues and the residues responsible for formation of the protein's structure. The phylogenetic history of a protein family makes such analysis particularly hard. Here we present two methods for finding the SDRs and the co-evolving residues (CERs) in proteins. We use a Monte Carlo approach for statistical inference, allowing us to reveal specific evolutionary patterns of SDRs and CERs. We apply these methods to study specific recognition in the bacterial two-component system and in the class Ia aminoacyl-tRNA synthetases. Our results agree well with structural information and the experimental analyses of these systems. Our results point at the complex and distinct patterns characteristic of the evolution of specificity in these systems.

Structural insights into the substrate specificity and activity of ervatamins, the papain-like cysteine proteases from a tropical plant, Ervatamia coronaria

Multiple proteases of the same family are quite often present in the same species in biological systems. These multiple proteases, despite having high homology in their primary and tertiary structures, show deviations in properties such as stability, activity, and specificity. It is of interest, therefore, to compare the structures of these multiple proteases in a single species to identify the structural changes, if any, that may be responsible for such deviations. Ervatamin-A, ervatamin-B and ervatamin-C are three such papain-like cysteine proteases found in the latex of the tropical plant Ervatamia coronaria, and are known not only for their high stability over a wide range of temperature and pH, but also for variations in activity and specificity among themselves and among other members of the family. Here we report the crystal structures of ervatamin-A and ervatamin-C, complexed with an irreversible inhibitor 1-[l-N-(trans-epoxysuccinyl)leucyl]amino-4-guanidinobutane (E-64), together with enzyme kinetics and molecular dynamic simulation studies. A comparison of these results with the earlier structures helps in a correlation of the structural features with the corresponding functional properties. The specificity constants (k(cat)/K(m)) for the ervatamins indicate that all of these enzymes have specificity for a branched hydrophobic residue at the P2 position of the peptide substrates, with different degrees of efficiency. A single amino acid change, as compared to ervatamin-C, in the S2 pocket of ervatamin-A (Ala67-->Tyr) results in a 57-fold increase in its k(cat)/K(m) value for a substrate having a Val at the P2 position. Our studies indicate a higher enzymatic activity of ervatamin-A, which has been subsequently explained at the molecular level from the three-dimensional structure of the enzyme and in the context of its helix polarizibility and active site plasticity.

Detection of pairwise residue proximity by covariation analysis for 3D-structure prediction of G-protein-coupled receptors

G-protein-coupled receptor (GPCR) is one of the most important targets for medicines. Homology modeling based on the crystal structure of bovine rhodopsin is currently the most frequently used method for GPCR targeted drug design. Information about residue-residue contacts and the structural specificity in the subfamily is essential for constructing more precise 3D structures, to distinguish the structural differences between the template and targets. In this study, we adopted the covariation analysis to extract information about residue-residue interactions from the amino acid sequence. In the opsin family, a large number of adjacent covarying residue pairs were detected. The detected residue pairs have a strong tendency to gather in some regions important for the structure and function. These results suggest that the covariation analysis is practically utilized to detect adjacent residue pairs and also to apply for predicting functional sites. Analyses of other GPCR subfamilies, olfactory receptor and chemokine receptor families, demonstrated that some adjacent covarying residue pairs were common. Thus, the covariation analysis has possibilities in the substantial improvement of the 3D-structure modeling of GPCRs and in the detection of functional sites such as the ligand-binding sites.

CRASP: a program for analysis of coordinated substitutions in multiple alignments of protein sequences

Recent results suggest that during evolution certain substitutions at protein sites may occur in a coordinated manner due to interactions between amino acid residues. Information on these coordinated substitutions may be useful for analysis of protein structure and function. CRASP is an Internet-available software tool for the detection and analysis of coordinated substitutions in multiple alignments of protein sequences. The approach is based on estimation of the correlation coefficient between the values of a physicochemical parameter at a pair of positions of sequence alignment. The program enables the user to detect and analyze pairwise relationships between amino acid substitutions at protein sequence positions, estimate the contribution of the coordinated substitutions to the evolutionary invariance or variability in integral protein physicochemical characteristics such as the net charge of protein residues and hydrophobic core volume. The CRASP program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/crasp/.

Cysteine protease isoforms from Trypanosoma cruzi, cruzipain 2 and cruzain, present different substrate preference and susceptibility to inhibitors

Cysteine-proteinases from parasitic protozoa have been recently characterized as factors of virulence and pathogenicity in several human and veterinary diseases. In Chagas' disease, the chronic infection caused by Trypanosoma cruzi, structure-functional studies on cysteine proteases were thus far limited to the parasite's major isoform, a cathepsin L-like lysosomal protease designated as cruzipain, cruzain or GP57/51. Encoded by a large gene family, cruzipain is efficiently targeted by synthetic inhibitors, which prevent parasite intracellular growth and differentiation. We have previously demonstrated that the multicopy cruzipain gene family includes polymorphic sequences, which could encode functionally different isoforms. We report here a comparative kinetic study between cruzain, the archetype of the cruzipain family, and an isoform, termed cruzipain 2, which is expressed preferentially by the mammalian stages of T. cruzi. Heterologous expression of the catalytic domain of cruzipain 2 in Saccharomyces cerevisae yielded an enzyme that differs markedly from cruzain with respect to pH stability, substrate specificity and sensitivity to inhibition by natural and synthetic inhibitors of cysteine proteases. We suggest that the structural-functional diversification imparted by genetic polymorphism of cruzipain genes may have contributed to T. cruzi adaptation to vertebrate hosts.

Analysis of covariation in an SH3 domain sequence alignment: applications in tertiary contact prediction and the design of compensating hydrophobic core substitutions

We have analyzed sequence covariation in an alignment of 266 non-redundant SH3 domain sequences using chi-squared statistical methods. Artifactual covariations arising from close evolutionary relationships among certain sequence subgroups were eliminated using empirically derived sequence diversity thresholds. This covariation detection method was able to predict residue-residue contacts (side-chain centres of mass within 8 A) in the structure of the SH3 domain with an accuracy of 85 %, which is greater than that achieved in many previous covariation studies. In examining the positions involved most frequently in covariations, we discovered a dramatic over-representation of a subset of five hydrophobic core positions. This covariation information was used to design second and third site substitutions that could compensate for highly destabilizing hydrophobic core substitutions in the Fyn SH3 domain, thus providing experimental data to validate the covariation analysis. The testing of our covariation detection method on 15 other alignments showed that the accuracy of contact prediction is highly variable depending on which sequence alignment is used, and useful levels of prediction accuracy were obtained with only approximately one-third of alignments. The results presented here provide insight into the difficulties inherent in covariation analysis, and suggest that it may have limited usefulness in tertiary structure prediction. On the other hand, our ability to use covariation analysis to design stabilizing combinations of hydrophobic core substitutions attests to its potential utility for gaining deeper insight into the stability determinants and functional mechanisms of proteins with known three-dimensional structures.

Correlated mutations contain information about protein-protein interaction

Many proteins have evolved to form specific molecular complexes and the specificity of this interaction is essential for their function. The network of the necessary inter-residue contacts must consequently constrain the protein sequences to some extent. In other words, the sequence of an interacting protein must reflect the consequence of this process of adaptation. It is reasonable to assume that the sequence changes accumulated during the evolution of one of the interacting proteins must be compensated by changes in the other. Here we apply a method for detecting correlated changes in multiple sequence alignments to a set of interacting protein domains and show that positions where changes occur in a correlated fashion in the two interacting molecules tend to be close to the protein-protein interfaces. This leads to the possibility of developing a method for predicting contacting pairs of residues from the sequence alone. Such a method would not need the knowledge of the structure of the interacting proteins, and hence would be both radically different and more widely applicable than traditional docking methods. We indeed demonstrate here that the information about correlated sequence changes is sufficient to single out the right inter-domain docking solution amongst many wrong alternatives of two-domain proteins. The same approach is also used here in one case (haemoglobin) where we attempt to predict the interface of two different proteins rather than two protein domains. Finally, we report here a prediction about the inter-domain contact regions of the heat- shock protein Hsc70 based only on sequence information.

Role of the occluding loop in cathepsin B activity

Within the lysosomal cysteine protease family, cathepsin B is unique due to its ability to act both as an endopeptidase and a peptidyldipeptidase. This latter capacity to remove C-terminal dipeptides has been attributed to the presence of a 20-residue insertion, termed the occluding loop, that blocks the primed terminus of the active site cleft. Variants of human procathepsin B, where all or part of this element was deleted, were expressed in the yeast Pichia pastoris. A mutant, where the 12 central residues of the occluding loop were deleted, autoprocessed, albeit more slowly than the wild type proenzyme, to yield a mature form of the enzyme with endopeptidase activity comparable with the wild-type cathepsin B, but totally lacking exopeptidase activity. This deletion mutant showed a 40-fold higher affinity for the inhibitor cystatin C, suggesting that the occluding loop normally restricts access of this inhibitor to the active site. In addition, the binding affinity of the cathepsin B propeptide, which is a potent inhibitor of this enzyme, was 50-fold increased, consistent with the finding that the loop reorients on activation of the proenzyme. These results suggest that the endopeptidase activity of cathepsin B is an evolutionary remnant since, as a consequence of its membership in the papain family, the propeptide must be able to bind unobstructed through the full length of the active site cleft.

Contribution to activity of histidine-aromatic, amide-aromatic, and aromatic-aromatic interactions in the extended catalytic site of cysteine proteinases

Within the papain family of cysteine proteinases few other residues in addition to the catalytic triad, Cys25-His159-Asn175 (papain numbering) are completely conserved [Berti & Storer (1995) J. Mol. Biol. 246, 273-283]. One such residue is tryptophan 177 which participates in a Trp-His-type interaction with the catalytic His159. In all enzymes of this class for which a three-dimensional structure has been reported, an additional highly conserved tryptophan, Trp181, also interacts with Trp177 via an aromatic-aromatic interaction in which the planes of the indole rings are essentially perpendicular. Also, both indole rings participate as pseudo-hydrogen bond acceptors in interactions with the two side chain amide protons of Asn175. Clearly, the proximity of Trp177 and Trp181 to the catalytic triad residues His159 and Asn175 and their network of interactions points to potential contributions of these aromatic residues to catalysis. In this paper, using cathepsin S, a naturally occurring variant that has a phenylalanine residue at position 181, we report the kinetic characterization of mutants of residues 175, 177, and 181. The results are interpreted in terms of the side chain contributions to catalytic activity and thiolate-imidazolium ion-pair stability. For example, the side chain of Asn175 has a major influence on the ion-pair stability presumably through its hydrogen bond to His159. The magnitude of this effect is modulated by Trp177, which shields the His159-Asn175 hydrogen bond from solvent. The His159-Trp177 interaction also contributes significantly to ion-pair stability; however, Trp181 and its interactions with Asn175 and Trp177 do not influence ion-pair stability to a significant degree. The observation that certain mutations at positions 177 and 181 result in a reduction of kcat/Km but do not appear to influence ion-pair stability probably reflects the contributions of these residues to substrate binding.

Growth requirements of endothelial cells in culture: variations in serum and amino acid concentrations

Endothelial cell growth in vitro is limited to the availability of nutrients from commercially available media and added serum. Nutrients, such as amino acids, are chiefly derived from the cell culture medium, rather than from added serum, and optimal endothelial cell growth may be dependent on amino acid levels in the culture media. To test this hypothesis, porcine pulmonary artery-derived endothelial cells were exposed to culture medium 199 (M199), amino acid-deficient M199 (dM199), as well as dM199 supplemented with amino acids. Cell protein was similar in cells cultured for 3 d in M199 supplemented with 1, 3, 5 or 10% bovine calf serum, respectively. Addition of amino acid solutions (L-amino acids [Laa], DL-amino acids [DLaa], 2Laa, or Laa+glutamine) to dM199 demonstrated a cell dependence for optimal growth on the type of amino acids as well as on the total available nitrogen in the media. Compared with M199, dM199 supplemented with Laa only partially supported long-term growth of endothelial cells in culture. On the other hand, dM199 supplemented with either 2Laa, DLaa, or Laa+ glutamine was superior over M199 with regard to endothelial cell growth. The addition of Laa+glutamine to dM199 was most growth-supporting, with an increase of over 2.6-fold in total cell protein compared with cells cultured with M199. These results suggest that, in addition to the presence of essential amino acids, total available nitrogen in culture media may be a critical factor for optimal endothelial cell growth.

Structural and functional roles of asparagine 175 in the cysteine protease papain

The role of the asparagine residue in the Cys-His-Asn "catalytic triad" of cysteine proteases has been investigated by replacing Asn175 in papain by alanine and glutamine using site-directed mutagenesis. The mutants were expressed in yeast and kinetic parameters determined against the substrate carbobenzoxy-L-phenylalanyl-(7-amino-4-methylcoumarinyl)- L-arginine. At the optimal pH of 6.5, the specificity constant (k(cat)/KM)obs was reduced by factors of 3.4 and 150 for the Asn175-->Gln and Asn175-->Ala mutants, respectively. Most of this effect was the result of a decrease in k(cat), as neither mutation significantly affected KM. Substrate hydrolysis by these mutants is still much faster than the non-catalytic rate, and therefore Asn175 cannot be considered as an essential catalytic residue in the cysteine protease papain. Detailed analyses of the pH activity profiles for both mutants allow the evaluation of the role of the Asn175 side chain on the stability of the active site ion pair and on the intrinsic activity of the enzyme. Alteration of the side chain at position 175 was also found to increase aggregation and proteolytic susceptibility of the proenzyme and to affect the thermal stability of the mature enzyme, reflecting a contribution of the asparagine residue to the structural integrity of papain. The strict conservation of Asn175 in cysteine proteases might therefore result from a combination of functional and structural constraints.

Sequence conservation and correlation measures in protein structure prediction

The rapid elucidation of protein sequences has allowed multiple sequence alignments to be calculated for a wide variety of proteins. Such alignments reveal positions that exhibit amino acid conservation--either of specific chemical groups in active and binding sites or of the more chemically inert hydrophobic residues that contribute to the protein core. The latter can provide constraints on the position of the protein chain and any local periodicity can suggest the type of secondary structure. Conservation measures, however, cannot provide specific pairwise packing information (each conserved hydrophobic position might pack against any other). However, if correlated changes between positions were observed then specific pairs of residue could be identified as interacting and therefore probably spatially adjacent. Most 'observations' of correlated changes have been anecdotal and of the few systematic studies that have been made, most have mistakenly incorporated a strong bias towards selecting conserved positions. When the conservation effect is separated (as best as possible) then little correlation signal remains to help identify adjacent positions.