Using increment of diversity to predict mitochondrial proteins of malaria parasite: integrating pseudo-amino acid composition and structural alphabet

Due to the complexity of Plasmodium falciparum (PF) genome, predicting mitochondrial proteins of PF is more difficult than other species. In this study, using the n-peptide composition of reduced amino acid alphabet (RAAA) obtained from structural alphabet named Protein Blocks as feature parameter, the increment of diversity (ID) is firstly developed to predict mitochondrial proteins. By choosing the 1-peptide compositions on the N-terminal regions with 20 residues as the only input vector, the prediction performance achieves 86.86% accuracy with 0.69 Mathew’s correlation coefficient (MCC) by the jackknife test. Moreover, by combining with the hydropathy distribution along protein sequence and several reduced amino acid alphabets, we achieved maximum MCC 0.82 with accuracy 92% in the jackknife test by using the developed ID model. When evaluating on an independent dataset our method performs better than existing methods. The results indicate that the ID is a simple and efficient prediction method for mitochondrial proteins of malaria parasite.     

Function, structure and regulation of cyanobacterial and chloroplast ATP synthase

The proton-translocating ATP synthase from chloroplasts and cyanobacteria forms ATP upon photosynthetic electron transport by using the proton gradient across the thylakoid membrane. Both enzymes contain nine different subunits and from the similarity in gene organisation and the high degree of amino acid sequence homology of the subunits it appears that these ATP synthases might have a common ancestor. Both enzymes need to be activated by membrane energisation in order to perform catalytic activity but, in contrast to the chloroplast ATP synthase, that from the studied cyanobacteria (with the exception of Spirulina platensis ) shows no effect of the redox state on activation. Functionally, the cyanobacterial enzyme corresponds to the reduced form of the chloroplast ATP synthase. In the chloroplast enzyme a stretch of 9 amino acids, including two cysteines in the γ-subunit, is involved in this redox effect and this stretch is absent in cyanobacteria. With γ-mutants from the cyanobacterium Synechocystis 6803 the role of this stretch is studied. When active, both the cyanobacterial and the reduced chloroplast ATP synthase transport 4 protons per ATP synthesised and hydrolysed. This ratio may depend on the environment of the enzyme such as protein and lipid composition and pH.     

Structural classification and properties of ketoacyl synthases

Ketoacyl synthases (KSs) catalyze condensing reactions combining acyl-CoA or acyl-acyl carrier protein (acyl-ACP) with malonyl-CoA to form 3-ketoacyl-CoA or with malonyl-ACP to form 3-ketoacyl-ACP. In each case, the resulting acyl chain is two carbon atoms longer than before, and CO2 and either CoA or ACP are formed. KSs also join other activated molecules in the polyketide synthesis cycle. Our classification of KSs by their primary and tertiary structures instead of by their substrates and the reactions that they catalyze enhances insights into this enzyme group. KSs fall into five families separated by their characteristic primary structures, each having members with the same catalytic residues, mechanisms, and tertiary structures. KS1 members, overwhelmingly named 3-ketoacyl-ACP synthase III or its variants, are produced predominantly by bacteria. Members of KS2 are mainly produced by plants, and they are usually long-chain fatty acid elongases/condensing enzymes and 3-ketoacyl-CoA synthases. KS3, a very large family, is composed of bacterial and eukaryotic 3-ketoacyl-ACP synthases I and II, often found in multidomain fatty acid and polyketide synthases. Most of the chalcone synthases, stilbene synthases, and naringenin-chalcone synthases in KS4 are from eukaryota. KS5 members are all from eukaryota, most are produced by animals, and they are mainly fatty acid elongases. All families except KS3 are split into subfamilies whose members have statistically significant differences in their primary structures. KS1 through KS4 appear to be part of the same clan. KS sequences, tertiary structures, and family classifications are available on the continuously updated ThYme (Thioester-active enzYme) database.     

Critical aspartic acid residues in pseudouridine synthases.

The pseudouridine synthases catalyze the isomerization of uridine to pseudouridine at particular positions in certain RNA molecules. Genomic data base searches and sequence alignments using the first four identified pseudouridine synthases led Koonin (Koonin, E. V. (1996) Nucleic Acids Res. 24, 2411-2415) and, independently, Santi and co-workers (Gustafsson, C., Reid, R., Greene, P. J., and Santi, D. V. (1996) Nucleic Acids Res. 24, 3756-3762) to group this class of enzyme into four families, which display no statistically significant global sequence similarity to each other. Upon further scrutiny (Huang, H. L., Pookanjanatavip, M., Gu, X. G., and Santi, D. V. (1998) Biochemistry 37, 344-351), the Santi group discovered that a single aspartic acid residue is the only amino acid present in all of the aligned sequences; they then demonstrated that this aspartic acid residue is catalytically essential in one pseudouridine synthase. To test the functional significance of the sequence alignments in light of the global dissimilarity between the pseudouridine synthase families, we changed the aspartic acid residue in representatives of two additional families to both alanine and cysteine: the mutant enzymes are catalytically inactive but retain the ability to bind tRNA substrate. We have also verified that the mutant enzymes do not release uracil from the substrate at a rate significant relative to turnover by the wild-type pseudouridine synthases. Our results clearly show that the aligned aspartic acid residue is critical for the catalytic activity of pseudouridine synthases from two additional families of these enzymes, supporting the predictive power of the sequence alignments and suggesting that the sequence motif containing the aligned aspartic acid residue might be a prerequisite for pseudouridine synthase function.     

Using K-minimum increment of diversity to predict secretory proteins of malaria parasite based on groupings of amino acids

Due to the complexity of Plasmodium falciparumis genome, predicting secretory proteins of P. falciparum is more difficult than other species. In this study, based on the measure of diversity definition, a new K-nearest neighbor method, K-minimum increment of diversity (K-MID), is introduced to predict secretory proteins. The prediction performance of the K-MID by using amino acids composition as the only input vector achieves 88.89% accuracy with 0.78 Mathew’s correlation coefficient (MCC). Further, the several reduced amino acids alphabets are applied to predict secretory proteins and the results show that the prediction results are improved to 90.67% accuracy with 0.83 MCC by using the 169 dipeptide compositions of the reduced amino acids alphabets obtained from Protein Blocks method.     

Insect chitin synthases: a review

Chitin is the most widespread amino polysaccharide in nature. The annual global amount of chitin is believed to be only one order of magnitude less than that of cellulose. It is a linear polymer composed of N-acetylglucosamines that are joined in a reaction catalyzed by the membrane-integral enzyme chitin synthase, a member of the family 2 of glycosyltransferases. The polymerization requires UDP–N-acetylglucosamines as a substrate and divalent cations as co-factors. Chitin formation can be divided into three distinct steps. In the first step, the enzymes‘ catalytic domain facing the cytoplasmic site forms the polymer. The second step involves the translocation of the nascent polymer across the membrane and its release into the extracellular space. The third step completes the process as single polymers spontaneously assemble to form crystalline microfibrils. In subsequent reactions the microfibrils combine with other sugars, proteins, glycoproteins and proteoglycans to form fungal septa and cell walls as well as arthropod cuticles and peritrophic matrices, notably in crustaceans and insects. In spite of the good effort by a hardy few, our present knowledge of the structure, topology and catalytic mechanism of chitin synthases is rather limited. Gaps remain in understanding chitin synthase biosynthesis, enzyme trafficking, regulation of enzyme activity, translocation of chitin chains across cell membranes, fibrillogenesis and the interaction of microfibrils with other components of the extracellular matrix. However, cumulating genomic data on chitin synthase genes and new experimental approaches allow increasingly clearer views of chitin synthase function and its regulation, and consequently chitin biosynthesis. In the present review, I will summarize recent advances in elucidating the structure, regulation and function of insect chitin synthases as they relate to what is known about fungal chitin synthases and other glycosyltransferases.     

Putrescine <Emphasis Type="Italic">N</Emphasis>-methyltransferases—a structure–function analysis

Putrescine N-methyltransferase (PMT) is a key enzyme of plant secondary metabolism at the start of the specific biosynthesis of nicotine, of tropane alkaloids, and of calystegines that are glycosidase inhibitors with nortropane structure. PMT is assumed to have developed from spermidine synthases (SPDS) participating in ubiquitous polyamine metabolism. In this study decisive differences between both enzyme families are elucidated. PMT sequences were known from four Solanaceae genera only, therefore additional eight PMT cDNA sequences were cloned from five Solanaceae and a Convolvulaceae. The encoded polypeptides displayed between 76% and 97% identity and typical amino acids different from plant spermidine synthase protein sequences. Heterologous expression of all enzymes proved catalytic activity exclusively as PMT and K _cat values between 0.16 s⁻¹ and 0.39 s⁻¹. The active site of PMT was initially inferred from a protein structure of spermidine synthase obtained by protein crystallisation. Those amino acids of the active site that were continuously different between PMTs and SPDS were mutated in one of the PMT sequences with the idea of changing PMT activity into spermidine synthase. Mutagenesis of active site residues unexpectedly resulted in a complete loss of catalytic activity. A protein model of PMT was based on the crystal structure of SPDS and suggests that overall protein folds are comparable. The respective cosubstrates S-adenosylmethionine and decarboxylated S-adenosylmethionine, however, appear to bind differentially to the active sites of both enzymes, and the substrate putrescine adopts a different position.     

Isolation and characterization of a reducing polyketide synthase gene from the lichen-forming fungus Usnea longissima

The reducing polyketide synthases found in filamentous fungi are involved in the biosynthesis of many drugs and toxins. Lichens produce bioactive polyketides, but the roles of reducing polyketide synthases in lichens remain to be clearly elucidated. In this study, a reducing polyketide synthase gene (U1PKS3) was isolated and characterized from a cultured mycobiont of Usnea longissima. Complete sequence information regarding U1PKS3 (6,519 bp) was obtained by screening a fosmid genomic library. A U1PKS3 sequence analysis suggested that it contains features of a reducing fungal type I polyketide synthase with β-ketoacyl synthase (KS), acyltransferase (AT), dehydratase (DH), enoyl reductase (ER), ketoacyl reducatse (KR), and acyl carrier protein (ACP) domains. This domain structure was similar to the structure of ccRadsl, which is known to be involved in resorcylic acid lactone biosynthesis in Chaetomium chiversii. The results of phylogenetic analysis located U1PKS3 in the clade of reducing polyketide synthases. RT-PCR analysis results demonstrated that UIPKS3 had six intervening introns and that UIPKS3 expression was upregulated by glucose, sorbitol, inositol, and mannitol.     

Accuracy of sequence alignment and fold assessment using reduced amino acid alphabets

Reduced or simplified amino acid alphabets group the 20 naturally occurring amino acids into a smaller number of representative protein residues. To date, several reduced amino acid alphabets have been proposed, which have been derived and optimized by a variety of methods. The resulting reduced amino acid alphabets have been applied to pattern recognition, generation of consensus sequences from multiple alignments, protein folding, and protein structure prediction. In this work, amino acid substitution matrices and statistical potentials were derived based on several reduced amino acid alphabets and their performance assessed in a large benchmark for the tasks of sequence alignment and fold assessment of protein structure models, using as a reference frame the standard alphabet of 20 amino acids. The results showed that a large reduction in the total number of residue types does not necessarily translate into a significant loss of discriminative power for sequence alignment and fold assessment. Therefore, some definitions of a few residue types are able to encode most of the relevant sequence/structure information that is present in the 20 standard amino acids. Based on these results, we suggest that the use of reduced amino acid alphabets may allow to increasing the accuracy of current substitution matrices and statistical potentials for the prediction of protein structure of remote homologs.     

Structure–function studies on the complex iron–sulfur flavoprotein glutamate synthase: the key enzyme of ammonia assimilation

Glutamate synthases are complex iron–sulfur flavoproteins that participate in the essential ammonia assimilation pathway in microorganisms and plants. The recent determination of the 3-dimensional structures of the α subunit of the NADPH-dependent glutamate synthase form and of the ferredoxin-dependent enzyme of Synechocystis sp. PCC 6803 provides a framework for the interpretation of the functional properties of these enzymes, and highlights protein segments most likely involved in control and coordination of the partial catalytic activities of glutamate synthases, which take place at sites distant from each other in space. In this review, we focus on the current knowledge on structure–function relationships in glutamate synthases, and we discuss open questions on the mechanisms of control of the enzyme reaction and of electron transfer among the enzyme flavin cofactors and iron–sulfur clusters.     

Clustering of protein families into functional subtypes using Relative Complexity Measure with reduced amino acid alphabets

Background  

Phylogenetic analysis can be used to divide a protein family into subfamilies in the absence of experimental information. Most phylogenetic analysis methods utilize multiple alignment of sequences and are based on an evolutionary model. However, multiple alignment is not an automated procedure and requires human intervention to maintain alignment integrity and to produce phylogenies consistent with the functional splits in underlying sequences. To address this problem, we propose to use the alignment-free Relative Complexity Measure (RCM) combined with reduced amino acid alphabets to cluster protein families into functional subtypes purely on sequence criteria. Comparison with an alignment-based approach was also carried out to test the quality of the clustering.     

Mutational analysis of a monoterpene synthase reaction: altered catalysis through directed mutagenesis of (-)-pinene synthase from Abies grandis

Two monoterpene synthases, (-)-pinene synthase and (-)-camphene synthase, from grand fir (Abies grandis) produce different product mixtures despite having highly homologous amino acid sequences and, presumably, very similar three-dimensional structures. The major product of (-)-camphene synthase, (-)-camphene, and the major products of (-)-pinene synthase, (-)-alpha-pinene, and (-)-beta-pinene, arise through distinct mechanistic variations of the electrophilic reaction cascade that is common to terpenoid synthases. Structural modeling followed by directed mutagenesis in (-)-pinene synthase was used to replace selected amino acid residues with the corresponding residues from (-)-camphene synthase in an effort to identify the amino acids responsible for the catalytic differences. This approach produced an enzyme in which more than half of the product was channeled through an alternative pathway. It was also shown that several (-)-pinene synthase to (-)-camphene synthase amino acid substitutions were necessary before catalysis was significantly altered. The data support a model in which the collective action of many key amino acids, located both in and distant from the active site pocket, regulate the course of the electrophilic reaction cascade.     

Role of cysteine residues in pseudouridine synthases of different families.

The pseudouridine synthases catalyze the isomerization of uridine to pseudouridine in RNA molecules. An attractive mechanism was proposed based on that of thymidylate synthase, in which the thiol(ate) group of a cysteine side chain serves as the nucleophile in a Michael addition to C6 of the isomerized uridine. Such a role for cysteine in the pseudouridine synthase TruA (also named Psi synthase I) has been discredited by site-directed mutagenesis, but sequence alignments have led to the conclusion that there are four distinct "families" of pseudouridine synthases that share no statistically significant global sequence similarity. It was, therefore, necessary to probe the role of cysteine residues in pseudouridine synthases of the families that do not include TruA. We examined the enzymes RluA and TruB, which are members of different families than TruA and each other. Substitution of cysteine for amino acids with nonnucleophilic side chains did not significantly alter the catalytic activity of either pseudouridine synthase. We conclude, therefore, that neither TruB nor RluA require thiol(ate) groups to effect catalysis, excluding their participation in a Michael addition to C6 of uridine, although not eliminating that mechanism (with an alternate nucleophile) from future consideration.     

Lipases provide a new mechanistic model for polyhydroxybutyrate (PHB) synthases: characterization of the functional residues in Chromatium vinosum PHB synthase

Polyhydroxybutyrate (PHB) synthases catalyze the conversion of beta-hydroxybutyryl coenzyme A (HBCoA) to PHB. These enzymes require an active site cysteine nucleophile for covalent catalysis. A protein BLASTp search using the Class III Chromatium vinosum synthase sequence reveals high homology to prokaryotic lipases whose crystal structures are known. The homology is very convincing in the alpha-beta-elbow (with the active site nucleophile)-alpha-beta structure, residues 131-175 of the synthase. A conserved histidine of the Class III PHB synthases aligns with the active site histidine of the lipases using the ClustalW algorithm. This is intriguing as this histidine is approximately 200 amino acids removed in sequence space from the catalytic nucleophile. Different threading algorithms suggest that the Class III synthases belong to the alpha/beta hydrolase superfamily which includes prokaryotic lipases. Mutagenesis studies were carried out on C. vinosum synthase C149, H331, H303, D302, and C130 residues. These studies reveal that H331 is the general base catalyst that activates the nucleophile, C149, for covalent catalysis. The model indicates that C130 is not involved in catalysis as previously proposed [Müh, U., Sinskey, A. J., Kirby, D. P., Lane, W. S., and Stubbe, J. (1999) Biochemistry 38, 826-837]. Studies with D302 mutants suggest D302 functions as a general base catalyst in activation of the 3-hydroxyl of HBCoA (or a hydroxybutyrate acyl enzyme) for nucleophilic attack on the covalently linked thiol ester intermediate. The relationship of the lipase model to previous models based on fatty acid synthases is discussed.     

Monoterpene and sesquiterpene synthases and the origin of terpene skeletal diversity in plants

The multitude of terpene carbon skeletons in plants is formed by enzymes known as terpene synthases. This review covers the monoterpene and sesquiterpene synthases presenting an up-to-date list of enzymes reported and evidence for their ability to form multiple products. The reaction mechanisms of these enzyme classes are described, and information on how terpene synthase proteins mediate catalysis is summarized. Correlations between specific amino acid motifs and terpene synthase function are described, including an analysis of the relationships between active site sequence and cyclization type and a discussion of whether specific protein features might facilitate multiple product formation.     

Using substitution matrices to estimate probability distributions for biological sequences.

Accurately estimating probabilities from observations is important for probabilistic-based approaches to problems in computational biology. In this paper we present a biologically-motivated method for estimating probability distributions over discrete alphabets from observations using a mixture model of common ancestors. The method is an extension of substitution matrix-based probability estimation methods. In contrast to previous such methods, our method has a simple Bayesian interpretation and has the advantage over Dirichlet mixtures that it is both effective and simple to compute for large alphabets. The method is applied to estimate amino acid probabilities based on observed counts in an alignment and is shown to perform comparably to previous methods. The method is also applied to estimate probability distributions over protein families and improves protein classification accuracy.     

Engineering of chimeric class II polyhydroxyalkanoate synthases

PHA synthase is a key enzyme involved in the biosynthesis of polyhydroxyalkanoates (PHAs). Using a combinatorial genetic strategy to create unique chimeric class II PHA synthases, we have obtained a number of novel chimeras which display improved catalytic properties. To engineer the chimeric PHA synthases, we constructed a synthetic phaC gene from Pseudomonas oleovorans (phaC1Po) that was devoid of an internal 540-bp fragment. Randomly amplified PCR products (created with primers based on conserved phaC sequences flanking the deleted internal fragment) were generated using genomic DNA isolated from soil and were substituted for the 540-bp internal region. The chimeric genes were expressed in a PHA-negative strain of Ralstonia eutropha, PHB(-)4 (DSM 541). Out of 1,478 recombinant clones screened for PHA production, we obtained five different chimeric phaC1Po genes that produced more PHA than the native phaC1Po. Chimeras S1-71, S4-8, S5-58, S3-69, and S3-44 exhibited 1.3-, 1.4-, 2.0-, 2.1-, and 3.0-fold-increased levels of in vivo activity, respectively. All of the mutants mediated the synthesis of PHAs with a slightly increased molar fraction of 3-hydroxyoctanoate; however, the weight-average molecular weights (Mw) of the PHAs in all cases remained almost the same. Based upon DNA sequence analyses, the various phaC fragments appear to have originated from Pseudomonas fluorescens and Pseudomonas aureofaciens. The amino acid sequence analyses showed that the chimeric proteins had 17 to 20 amino acid differences from the wild-type phaC1Po, and these differences were clustered in the same positions in the five chimeric clones. A threading model of PhaC1Po, developed based on homology of the enzyme to the Burkholderia glumae lipase, suggested that the amino acid substitutions found in the active chimeras were located mostly on the protein model surface. Thus, our combinatorial genetic engineering strategy proved to be broadly useful for improving the catalytic activities of PHA synthase enzymes.     

Cloning of an Arabidopsis thaliana cDNA coding for farnesyl diphosphate synthase by functional complementation in yeast

A cDNA encoding farnesyl diphosphate synthase, an enzyme that synthesizes C15 isoprenoid diphosphate from isopentenyl diphosphate and dimethylallyl diphosphate, was cloned from an Arabidopsis thaliana cDNA library by complementation of a mutant of Saccharomyces cerevisiae deficient in this enzyme. The A. thaliana cDNA was also able to complement the lethal phenotype of the erg20 deletion yeast mutant. As deduced from the full-length 1.22 kb cDNA nucleotide sequence, the polypeptide contains 343 amino acids and has a relative molecular mass of 39689. The predicted amino acid sequence presents about 50% identity with the yeast, rat and human FPP synthases. Southern blot analyses indicate that A. thaliana probably contains a single gene for farnesyl diphosphate synthase.     

On reduced amino acid alphabets for phylogenetic inference

We investigate the use of Markov models of evolution for reduced amino acid alphabets or bins of amino acids. The use of reduced amino acid alphabets can ameliorate effects of model misspecification and saturation. We present algorithms for 2 different ways of automating the construction of bins: minimizing criteria based on properties of rate matrices and minimizing criteria based on properties of alignments. By simulation, we show that in the absence of model misspecification, the loss of information due to binning is found to be insubstantial, and the use of Markov models at the binned level is found to be almost as effective as the more appropriate missing data approach. By applying these approaches to real data sets where compositional heterogeneity and/or saturation appear to be causing biased tree estimation, we find that binning can improve topological estimation in practice.