Modulating the proteinase inhibitory profile of a plant cystatin by single mutations at positively selected amino acid sites

Cysteine proteinase inhibitors of the cystatin superfamily have several important functions in plants, including the inhibition of exogenous cysteine proteinases during herbivory or infection. Here we used a maximum-likelihood approach to assess whether plant cystatins, like other proteins implicated in host-pest interactions, have been subject to positive selection during the course of their evolution. Several amino acid sites were identified as being positively selected in cystatins from either Poaceae (monocots) and Solanaceae (dicots). These hypervariable sites were located at strategic positions on the protein: on each side of the conserved glycine residues in the N-terminal trunk, within the first and second inhibitory loops entering the active site of target enzymes, and surrounding the larfav motif, a sequence of unknown function conserved among plant cystatins. Supporting the assumption that positively selected, hypervariable sites are indicative of amino acid sites implicated in functional diversity, mutants of the 8th cystatin unit of tomato multicystatin including alternative residues at positively selected sites in the N-terminal trunk exhibited highly variable affinities for the cysteine proteases papain, cathepsin B and cathepsin H. Overall, these observations support the hypothesis that plant cystatins have been under selective pressure to evolve in response to predatory challenges by herbivorous enemies. They also indicate the potential of site-directed mutagenesis at positively selected sites for the generation of cystatins with improved binding properties.     

Altering the specificity of subtilisin Bacillus lentus through the introduction of positive charge at single amino acid sites

The use of methanethiosulfonates as thiol-specific modifying reagents in the strategy of combined site-directed mutagenesis and chemical modification allows virtually unlimited opportunities for creating new protein surface environments. As a consequence of our interest in electrostatic manipulation as a means of tailoring enzyme activity and specificity, we have recently adopted this approach for the controlled incorporation of multiple negative charges at single sites in the representative serine protease, subtilisin Bacillus lentus (SBL). We now describe the use of this strategy to introduce multiple positive charges. A series of mono-, di- and triammonium methanethiosulfonates were synthesized and used to modify cysteine mutants of SBL at positions 62 in the S2 site, 156 and 166 in the S1 site and 217 in the S1' site. Kinetic parameters for these chemically modified mutants (CMM) enzymes were determined at pH 8.6. The presence of up to three positive charges in the S1, S1' and S2 subsites of SBL resulted in up to 77-fold lowered activity, possibly due to interference with the histidinium ion formed in the transition state of the hydrolytic reactions catalyzed.     

Inhibition of Plasmodium falciparum cysteine proteases by the sugarcane cystatin CaneCPI-4

Malaria is a disease caused by Plasmodium parasites that affects hundreds of millions of people. Plasmodium proteases are involved in invasion, erythrocyte egress and degradation of host proteins. Falcipains are well-studied cysteine peptidases located in P. falciparum food vacuoles that participate in hemoglobin degradation. Cystatins are natural cysteine protease inhibitors that are implicated in a wide range of regulatory processes. Here, we report that a cystatin from sugarcane, CaneCPI-4, is selectively internalized into P. falciparum infected erythrocytes and is not processed by the parasite proteolytic machinery. Furthermore, we demonstrated the inhibition of P. falciparum cysteine proteases by CaneCPI-4, suggesting that it can exert inhibitory functions inside the parasites. The inhibition of the proteolytic activity of parasite cells is specific to this cystatin, as the addition of an anti-CaneCPI-4 antibody completely abolished the inhibition. We extended the studies to recombinant falcipain-2 and falcipain-3 and demonstrated that CaneCPI-4 strongly inhibits these enzymes, with IC₅₀ values of 12 nM and 42 nM, respectively. We also demonstrated that CaneCPI-4 decreased the hemozoin formation in the parasites, affecting the parasitemia. Taken together, this study identified a natural molecule as a potential antimalarial that specifically targets falcipains and also contributes to a better understanding of macromolecule acquisition by Plasmodium falciparum infected RBCs.     

Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene.

Several codon-based models for the evolution of protein-coding DNA sequences are developed that account for varying selection intensity among amino acid sites. The "neutral model" assumes two categories of sites at which amino acid replacements are either neutral or deleterious. The "positive-selection model" assumes an additional category of positively selected sites at which nonsynonymous substitutions occur at a higher rate than synonymous ones. This model is also used to identify target sites for positive selection. The models are applied to a data set of the V3 region of the HIV-1 envelope gene, sequenced at different years after the infection of one patient. The results provide strong support for variable selection intensity among amino acid sites The neutral model is rejected in favor of the positive-selection model, indicating the operation of positive selection in the region. Positively selected sites are found in both the V3 region and the flanking regions.     

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

The engineering of a protein containing an alternative local residue packing for a set of side-chains has proven to be a major challenge because compositional, volumetric and steric constraints must be respected. Homologous proteins should provide examples of alternative groups of residues leading to a similar functional result. The functional significance of a pair of co-ordinated changes that are observed in the cysteine proteases family has been investigated by comparing the effect of individual or double changes on secretion, stability and activity of papain. The two changes are not independent. Detrimental effects of single mutations at one of the two positions can be partly suppressed by the co-ordinated mutation that reproduces naturally occurring contacts, indicating that these changes are concerted. Single mutations at the other position produce milder effects, suggesting a pathway for evolution.     

The specificity of TIMP-2 for matrix metalloproteinases can be modified by single amino acid mutations.

Residues 1-127 of human TIMP-2 (N-TIMP-2), comprising three of the disulfide-bonded loops of the TIMP-2 molecule, is a discrete protein domain that folds independently of the C-terminal domain. This domain has been shown to be necessary and sufficient for metalloproteinase inhibition and contains the major sites of interaction with the catalytic N-terminal domain of active matrix metalloproteinases (MMPs). Residues identified as being involved in the interaction with MMPs by NMR chemical shift perturbation studies and TIMP/MMP crystal structures have been altered by site-directed mutagenesis. We show, by measurement of association rates and apparent inhibition constants, that the specificity of these N-TIMP-2 mutants for a range of MMPs can be altered by single site mutations in either the TIMP "ridge" (Cys1-Cys3 and Ser68-Cys72) or the flexible AB loop (Ser31-Ile41). This work demonstrates that it is possible to engineer TIMPs with altered specificity and suggests that this form of protein engineering may be useful in the treatment of diseases such as arthritis and cancer where the selective inhibition of key MMPs is desirable.     

A single amino acid change in the plant alternative oxidase alters the specificity of organic acid activation.

The alternative oxidase is a quinol oxidase of the respiratory chain of plants and some fungi and protists. Its activity is regulated by redox-sensitive disulphide bond formation between neighbouring subunits and direct interaction with certain alpha-ketoacids. To investigate these regulatory mechanisms, we undertook site-directed mutagenesis of soybean and Arabidopsis alternative oxidase cDNAs, and expressed them in tobacco plants and Escherichia coli, respectively. The homologous C99 and C127 residues of GmAOX3 and AtAOX1a, respectively, were changed to serine. In the plant system, this substitution prevented oxidative inactivation of alternative oxidase and rendered the protein insensitive to pyruvate activation, in agreement with the recent results from other laboratories [Rhoads et al. (1998) J. Biol. Chem. 273, 30750-30756; Vanlerberghe et al. (1998) Plant Cell 10, 1551-1560]. However, the mutated protein is instead activated specifically by succinate. Measurements of AtAOX1a activity in bacterial membranes lacking succinate dehydrogenase confirmed that the stimulation of the mutant protein's activity by succinate did not involve its metabolism. Examples of alternative oxidase proteins with the C to S substitution occur in nature and these oxidases are expected to be activated under most conditions in vivo, with implications for the efficiency of respiration in the tissues which express them.     

Tailoring the peptide-binding specificity of an RNA by combinations of specificity-altering mutations

In this study, the ability to tailor the peptide-binding specificity of an RNA was investigated. First, variants of the Rev-response element (RRE) RNA with different specificities toward the natural binding partner, Rev, and two RRE-binding aptamers, the RSG-1.2 and the Kl peptides, were identified. Next, hybrid RRE mutants with combinations of two sets of specificity-altering substitutions were tested for peptide-binding specificity. It was shown that in most cases the results of the combination of individual mutations were of an additive nature, therefore providing a way to manipulate the peptide-binding specificity of an RNA in a predictable manner.     

Maximum likelihood analysis of mammalian p53 indicates the presence of positively selected sites and higher tumorigenic mutations in purifying sites

The tumor suppressor gene TP53 (p53) maintains genome stability. Mutation or loss of p53 is found in most cancers. Analysis of evolutionary constrains and p53 mutations reveal important sites for concomitant functional studies. In this study, phylogenetic analyses of the coding sequences of p53 from 26 mammals were carried out by applying a maximum likelihood method. The results display two branches under adaptive evolution in mammals. Moreover, each codon of p53 was analyzed by the PAML method for presence of positively selected sites. PAML identified several statistically significant amino acids that undergo positive selection. The data indicates that amino acids responsible for the core functions of p53 are highly conserved, while positively selected sites are predominantly located in the N- and C-terminus of p53. Further analysis of evolutionary pressure and mutations showed the occurrence of more frequent tumorigenic mutations in purifying sites of p53.     

Alteration of a single amino acid changes the substrate specificity of dihydroflavonol 4-reductase

Many plant species exhibit a reduced range of flower colors due to the lack of an essential gene or to the substrate specificity of a biosynthetic enzyme. Petunia does not produce orange flowers because dihydroflavonol 4-reductase (DFR) from this species, an enzyme involved in anthocyanin biosynthesis, inefficiently reduces dihydrokaempferol, the precursor to orange pelargonidin-type anthocyanins. The substrate specificity of DFR, however, has not been investigated at the molecular level. By analyzing chimeric DFRs of Petunia and Gerbera, we identified a region that determines the substrate specificity of DFR. Furthermore, by changing a single amino acid in this presumed substrate-binding region, we developed a DFR enzyme that preferentially reduces dihydrokaempferol. Our results imply that the substrate specificity of DFR can be altered by minor changes in DFR.     

Subsite specificity of trypanosomal cathepsin L-like cysteine proteases. Probing the S2 pocket with phenylalanine-derived amino acids.

The S2 subsite of mammalian cysteine proteinases of the papain family is essential for specificity. Among natural amino acids, all these enzymes prefer bulky hydrophobic residues such as phenylalanine at P2. This holds true for their trypanosomal counterparts: cruzain from Trypanosoma cruzi and congopain from T. congolense. A detailed analysis of the S2 specificity of parasitic proteases was performed to gain information that might be of interest for the design of more selective pseudopeptidyl inhibitors. Nonproteogenic phenylalanyl analogs (Xaa) have been introduced into position P2 of fluorogenic substrates dansyl-Xaa-Arg-Ala-Pro-Trp, and their kinetic constants (Km, kcat/Km) have been determined with congopain and cruzain, and related host cathepsins B and L. Trypanosomal cysteine proteases are poorly stereoselective towards D/L-Phe, the inversion of chirality modifying the efficiency of the reaction but not the Km. Congopain binds cyclohexylalanine better than aromatic Phe derivatives. Another characteristic feature of congopain compared to cruzain and cathepsins B and L was that it could accomodate a phenylglycyl residue (kcat/Km = 1300 mM-1.s-1), while lengthening of the side chain by a methylene group only slightly impaired the specificity constant towards trypanosomal cysteine proteases. Mono- and di-halogenation or nitration of Phe did not affect Km for cathepsin L-like enzymes, but the presence of constrained Phe derivatives prevented a correct fitting into the S2 subsite. A model of congopain has been built to study the fit of Phe analogs within the S2 pocket. Phe analogs adopted a positioning within the S2 pocket similar to that of the Tyr of the cruzain/Z-Tyr-Ala-fluoromethylketone complex. However, cyclohexylalanine has an energetically favorable chair-like conformation and can penetrate deeper into the subsite. Fitting of modeled Phe analogs were in good agreement with kinetic parameters. Furthermore, a linear relationship could be established with logP, supporting the suggestion that fitting into the S2 pocket of trypanosomal cysteine proteases depends on the hydrophobicity of Phe analogs.     

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.     

Enzymatic synthesis of the neuroexcitatory amino acid quisqualic acid by cysteine synthase

Purification of cysteine synthase from the leaves of Quisqualis indica var. villosa reveals the presence of two forms of this enzyme, separated by chromatography on DEAE-Sephadex A-50. Isoenzyme A was purified 10 000-fold and had a specific activity of 10.8 U/mg protein. Isoenzyme B was purified 460-fold with a specific activity of 0.49 U/mg protein. Both isoenzymes have the same M,s (58 000) and dissociate into identical subunits (M_r 29 000). The K_m value of isoenzyme A is 1.9 mM for O-acetyl-L-serine and 59 μM for sulphide, while that of isoenzyme B is 7.1 mM for O-acetyl-L-serine and 4.0 mM for 3,5-dioxo-1,2,4-oxadiazolidine. Both isoenzymes catalyse the formation of cysteine from O-acetyl-L-serine and hydrogen sulphide, but only isoenzyme B catalyses the formation of L-quisqualic acid. Other significant differences occur in the substrate specificity of the two isoenzymes. Some properties of the purified cysteine synthase isoenzymes are also described.     

Impact of plant nutrients on the relationship between a herbivorous insect and its symbiotic bacteria

The interactions between herbivorous insects and their symbiotic micro-organisms can be influenced by the plant species on which the insects are reared, but the underlying mechanisms are not understood. Here, we identify plant nutrients, specifically amino acids, as a candidate factor affecting the impact of symbiotic bacteria on the performance of the phloem-feeding aphid Aphis fabae. Aphis fabae grew more slowly on the labiate plant Lamium purpureum than on an alternative host plant Vicia faba, and the negative effect of L. purpureum on aphid growth was consistently exacerbated by the bacterial secondary symbionts Regiella insecticola and Hamiltonella defensa, which attained high densities in L. purpureum-reared aphids. The amino acid content of the phloem sap of L. purpureum was very low; and A. fabae on chemically defined diets of low amino acid content also grew slowly and had elevated secondary symbiont densities. It is suggested that the phloem nutrient profile of L. purpureum promotes deleterious traits in the secondary symbionts and disturbs insect controls over bacterial abundance.     

Mutagenesis by single site-specific arylamine-DNA adducts. Induction of mutations at multiple sites

Two related carcinogen adducts, N-(deoxyguanosin-8-yl)-2-aminofluorene (AF) or N-(deoxyguanosin-8-yl)-N-acetyl-2-aminofluorene (AAF), were introduced into the lacZ' gene at base position 6253 of the minus strand of M13mp9 viral DNA. The construction of this site-specifically modified DNA was accomplished by first preparing a gapped heteroduplex missing 7 nucleotides at position 6251-6257 followed by ligation with an unmodified heptamer or with a heptamer containing either an AF or AAF adduct. These site-specifically modified templates were transfected into competent wild-type Escherichia coli cells (JM103) and a uvrA strain (SMH12). The mutation spectrum was determined by phenotypic selection of colorless plaques indicating a defective beta-galactosidase marker enzyme and by an in situ hybridization procedure to detect single base pair mismatches in the adduct region. DNA sequencing was used to characterize 179 of the mutants obtained. We found that both adducts were capable of inducing base substitution mutations at the adduct site and in the local region of the adduct. A specific frameshift (+1G) was also observed at a displaced site. All of the frameshift mutations occurred at the ligation site of the modified oligonucleotide. Control experiments with an unmodified oligonucleotide did not show an enhancement of mutations at this site, indicating that the adducts may have been responsible for these frameshifts. The mutations spectra induced by these adducts suggest that mutagenesis depends not only on adduct structure but also the sequence in which the adduct is located and the host cell type used for mutation expression.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Control of the false discovery rate applied to the detection of positively selected amino acid sites

In this article, we consider the probabilistic identification of amino acid positions that evolve under positive selection as a multiple hypothesis testing problem. The null hypothesis "H0,s: site s evolves under a negative selection or under a neutral process of evolution" is tested at each codon site of the alignment of homologous coding sequences. Standard hypothesis testing is based on the control of the expected proportion of falsely rejected null hypotheses or type-I error rate. As the number of tests increases, however, the power of an individual test may become unacceptably low. Recent advances in statistics have shown that the false discovery rate--in this case, the expected proportion of sites that do not evolve under positive selection among those that are estimated to evolve under this selection regime--is a quantity that can be controlled. Keeping the proportion of false positives low among the significant results generally leads to an increase in power. In this article, we show that controlling the false detection rate is relevant when searching for positively selected sites. We also compare this new approach to traditional methods using extensive simulations.     

Internalization of exogenous cystatin F supresses cysteine proteases and induces the accumulation of single-chain cathepsin L by multiple mechanisms

Cystatin F is an unusual member of the cystatin family of protease inhibitors, which is made as an inactive dimer and becomes activated by proteolysis in the endo/lysosome pathway of the immune cells that produce it. However a proportion is secreted and can be taken up and activated by other cells. We show here that cystatin F acquired in this way induces a dramatic accumulation of the single-chain form of cathepsin L (CatL). Cystatin F was observed in the same cellular compartments as CatL and was tightly complexed with CatL as determined by co-precipitation studies. The observed accumulation of single-chain CatL was partly due to cystatin F-mediated inhibition of the putative single-chain to two-chain CatL convertase AEP/legumain and partly to general suppression of cathepsin activity. Thus, cystatin F stabilizes CatL leading to the dramatic accumulation of an inactive complex composed either of the single-chain or two-chain form depending on the capacity of cystatin F to inhibit AEP. Cross-transfer of cystatin F from one cell to another may therefore attenuate potentially harmful effects of excessive CatL activity while paradoxically, inducing accumulation of CatL protein. Finally, we confirmed earlier data (Beers, C., Honey, K., Fink, S., Forbush, K., and Rudensky, A. (2003) J. Exp. Med. 197, 169-179) showing a loss of CatL activity, but not of CatL protein, in macrophages activated with IFNγ. However, we found equivalent loss of CatL activity in wild type and cystatin F-null macrophages suggesting that an inhibitory activity other than cystatin F quenches CatL activity in activated macrophages.     

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.