RESUMO
The A1 family of eukaryotic aspartic proteinases (APs) forms one of the 16 AP families. Although one of the best characterized families, the recent increase in genome sequence data has revealed many fungal AP homologs with novel sequence characteristics. This study was performed to explore the fungal AP sequence space and to obtain an in-depth understanding of fungal AP evolution. Using a comprehensive phylogeny of approximately 700 AP sequences from the complete proteomes of 87 fungi and 20 nonfungal eukaryotes, 11 major clades of APs were defined of which clade I largely corresponds to the A1A subfamily of pepsin-archetype APs. Clade II largely corresponds to the A1B subfamily of nepenthesin-archetype APs. Remarkably, the nine other clades contain only fungal APs, thus indicating that fungal APs have undergone a large sequence diversification. The topology of the tree indicates that fungal APs have been subject to both "birth and death" evolution and "functional redundancy and diversification." This is substantiated by coclustering of certain functional sequence characteristics. A meta-analysis toward the identification of Cluster Determining Positions (CDPs) was performed in order to investigate the structural and biochemical basis for diversification. Seven CDPs contribute to the secondary structure of the enzyme. Three other CDPs are found in the vicinity of the substrate binding cleft. Tree topology, the large sequence variation among fungal APs, and the apparent functional diversification suggest that an amendment to update the current A1 AP classification based on a comprehensive phylogenetic clustering might contribute to refinement of the classification in the MEROPS peptidase database.
Assuntos
Ácido Aspártico Proteases/genética , Fungos/enzimologia , Fungos/genética , Filogenia , Sequência de Aminoácidos , Animais , Ácido Aspártico Proteases/química , Evolução Molecular , Fungos/química , Modelos Moleculares , Dados de Sequência Molecular , Proteoma/genética , Alinhamento de SequênciaRESUMO
The ascomycete plant pathogen Botrytis cinerea secretes aspartic proteinase (AP) activity. Functional analysis was carried out on five aspartic proteinase genes (Bcap1-5) reported previously. Single and double mutants lacking these five genes showed neither a reduced secreted proteolytic activity, nor a reduction in virulence and they showed no alteration in sensitivity to antifungal proteins purified from grape juice. Scrutiny of the B. cinerea genome revealed the presence of nine additional Bcap genes, denoted Bcap6-14. The product of the Bcap8 gene was found to constitute up to 23% of the total protein secreted by B. cinerea. Bcap8-deficient mutants secreted approximately 70% less AP activity but were just as virulent as the wild-type strain. Phylogenetic analysis showed that Bcap8 has orthologs in many basidiomycetes but only few ascomycetes including the biocontrol fungus Trichoderma harzanium. Potential functions of the 14 APs in B. cinerea are discussed based on their sequence characteristics, phylogeny and predicted localization.
Assuntos
Ácido Aspártico Proteases/metabolismo , Botrytis/enzimologia , Proteínas Fúngicas/metabolismo , Sequência de Aminoácidos , Antifúngicos/farmacologia , Ácido Aspártico Proteases/classificação , Ácido Aspártico Proteases/genética , Botrytis/efeitos dos fármacos , Botrytis/genética , Botrytis/patogenicidade , Clonagem Molecular , Citosol/enzimologia , Proteínas Fúngicas/classificação , Proteínas Fúngicas/genética , Deleção de Genes , Genes Fúngicos/genética , Dados de Sequência Molecular , Família Multigênica/genética , Filogenia , Doenças das Plantas/microbiologia , Proteínas de Plantas/farmacologia , Alinhamento de SequênciaRESUMO
The IA(3) polypeptide inhibitor from Saccharomyces cerevisiae interacts potently and selectively with its target, the S. cerevisiae vacuolar aspartic proteinase (ScPr). Upon encountering the enzyme, residues 2-32 of the intrinsically unstructured IA(3) polypeptide become ordered into an almost-perfect alpha-helix. In previous IA(3) mutagenesis studies, we identified important characteristics of the enzyme inhibitor interactions and generated a large dataset of variants with K(i) values determined experimentally at pH 3.1 and 4.7. Using this information, the three-dimensional structure of each variant was modelled in silico with the correct protonation for each experimental pH value. A set of descriptors of the inhibitor/ScPr interactions was then calculated and used to establish mathematical models relating the variant sequences to their inhibitory activities at each pH. Cross-validation, external-set validation and five separate selections of the training and test samples confirmed the robustness of the equations. A major contributor to the structure-activity relationship was the free energy of binding calculated by the FoldX program. The mathematical models were challenged further (i) by in silico alanine-scanning mutagenesis of residues 2-32 in IA(3) and relating binding energy to experimentally derived inhibition constants for selected representatives of these variants; and (ii) by predicting inhibitory-potencies for two novel IA(3)-variants. The predictions of the equations for these new IA(3)-variants with ScPr matched almost precisely the kinetic data determined experimentally. The models described represent valuable tools for the future design of novel inhibitor variants active against ScPr and other aspartic proteinases.