A novel Bayesian approach to drug design with simple and complex enzymes: Use with lens aldehyde dehydrogenase and the glyoxalase pathway

Purpose: Acquiring details of kinetic parameters of enzymes is crucial to biochemical understanding, drug development, and clinical diagnosis in ocular diseases. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. Methods: We have developed Bayesian utility functions to minimise kinetic parameter variance involving differentiation of model expressions and matrix inversion. These have been applied to the simple kinetics of the enzymes in the glyoxalase pathway (of importance in posttranslational modification of proteins in cataract), and the complex kinetics of lens aldehyde dehydrogenase (also of relevance to cataract). Results: Our successful application of Bayesian statistics has allowed us to identify a set of rules for designing optimum kinetic experiments iteratively. Most importantly, the distribution of points in the range is critical; it is not simply a matter of even or multiple increases. At least 60 % must be below the KM (or plural if more than one dissociation constant) and 40% above. This choice halves the variance found using a simple even spread across the range.With both the glyoxalase system and lens aldehyde dehydrogenase we have significantly improved the variance of kinetic parameter estimation while reducing the number and costs of experiments. Conclusions: We have developed an optimal and iterative method for selecting features of design such as substrate range, number of measurements and choice of intermediate points. Our novel approach minimises parameter error and costs, and maximises experimental efficiency. It is applicable to many areas of ocular drug design, including receptor-ligand binding and immunoglobulin binding, and should be an important tool in ocular drug discovery.

Molecular cloning and comparative analysis of four beta-galactosidase genes from Bifidobacterium bifidum NCIMB41171

Bifidobacterium bifidum NCIMB41171 carries four genes encoding different beta-galactosidases. One of them, named bbgIII, consisted of an open reading frame of 1,935 amino acid (a.a.) residues encoding a protein with a multidomain structure, commonly identified on cell wall bound enzymes, having a signal peptide, a membrane anchor, FIVAR domains, immunoglobulin Ig-like and discoidin-like domains. The other three genes, termed bbgI, bbgII and bbgIV, encoded proteins of 1,291, 689 and 1,052 a.a. residues, respectively, which were most probably intracellularly located. Two cases of protein evolution between strains of the same species were identified when the a.a. sequences of the BbgI and BbgIII were compared with homologous proteins from B. bifidum DSM20215. The homologous proteins were found to be differentiated at the C-terminal a.a. part either due to a single nucleotide insertion or to a whole DNA sequence insertion, respectively. The bbgIV gene was located in a gene organisation surrounded by divergently transcribed genes putatively for sugar transport (galactoside-symporter) and gene regulation (LacI-transcriptional regulator), a structure that was found to be highly conserved in B. longum, B. adolescentis and B. infantis, suggesting optimal organisation shared amongst those species.

Twists and turns: a tale of two shikimate-pathway enzymes

We are studying two enzymes from the shikimate pathway, dehydroquinate synthase (DHQS) and 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). Both enzymes have been the subject of numerous studies to elucidate their reaction mechanisms. Crystal structures of DHQS and EPSPS in the presence and absence of substrates, cofactors and/or inhibitors are now available. These structures reveal movements of domains, rearrangements of loops and changes in side-chain positions necessary for the formation of a catalytically competent active site. The potential for using complementary small-angle X-ray scattering (SAXS) studies to confirm the presence of these structural differences in solution has also been explored. Comparative analysis of crystal structures, in the presence and absence of ligands, has revealed structural features critical for substrate-binding and catalysis. We have also analysed these structures by generating GRID energy maps to detect favourable binding sites. The combination of X-ray crystallography, SAXS and computational techniques provides an enhanced analysis of structural features important for the function of these complex enzymes.

Platelets release novel thiol isomerase enzymes which are recruited to the cell surface following activation

The thiol isomerase enzymes protein disulphide isomerase (PDI) and endoplasmic reticulum protein 5 (ERp5) are released by resting and activated platelets. These re-associate with the cell surface where they modulate a range of platelet responses including adhesion, secretion and aggregation. Recent studies suggest the existence of yet uncharacterised platelet thiol isomerase proteins. This study aimed to identify which other thiol isomerase enzymes are present in human platelets. Through the use of immunoblotting, flow cytometry, cell-surface biotinylation and gene array analysis, we report the presence of five additional thiol isomerases in human and mouse platelets and megakaryocytes, namely; ERp57, ERp72, ERp44, ERp29 and TMX3. ERp72, ERp57, ERp44 and ERp29 are released by platelets and relocate to the cell surface following platelet activation. The transmembrane thiol isomerase TMX3 was also detected on the platelet surface but does not increase following activation. Extracellular PDI is also implicated in the regulation of coagulation by the modulation of tissue factor activity. ERp57 was identified within platelet-derived microparticle fractions, suggesting that ERp57 may also be involved in the regulation of coagulation as well as platelet function. These data collectively implicate the expanding family of platelet-surface thiol isomerases in the regulation of haemostasis.

Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

Background: Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites. Methodology and Principal Findings: Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed. Conclusions and Significance: Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.

Sequence and phylogenetic analysis of viper venom serine proteases

Snakebites are a major neglected tropical disease responsible for as many as 95000 deaths every year worldwide. Viper venom serine proteases disrupt haemostasis of prey and victims by affecting various stages of the blood coagulation system. A better understanding of their sequence, structure, function and phylogenetic relationships will improve the knowledge on the pathological conditions and aid in the development of novel therapeutics for treating snakebites. A large dataset for all available viper venom serine proteases was developed and analysed to study various features of these enzymes. Despite the large number of venom serine protease sequences available, only a small proportion of these have been functionally characterised. Although, they share some of the common features such as a C-terminal extension, GWG motif and disulphide linkages, they vary widely between each other in features such as isoelectric points, potential N-glycosylation sites and functional characteristics. Some of the serine proteases contain substitutions for one or more of the critical residues in catalytic triad or primary specificity pockets. Phylogenetic analysis clustered all the sequences in three major groups. The sequences with substitutions in catalytic triad or specificity pocket clustered together in separate groups. Our study provides the most complete information on viper venom serine proteases to date and improves the current knowledge on the sequence, structure, function and phylogenetic relationships of these enzymes. This collective analysis of venom serine proteases will help in understanding the complexity of envenomation and potential therapeutic avenues.

Metaproteomics analysis reveals the adaptation process for the chicken gut microbiota

The animal gastrointestinal tract houses a large microbial community, the gut microbiota, that confers many benefits to its host, such as protection from pathogens and provision of essential metabolites. Metagenomic approaches have defined the chicken fecal microbiota in other studies, but here, we wished to assess the correlation between the metagenome and the bacterial proteome in order to better understand the healthy chicken gut microbiota. Here, we performed high-throughput sequencing of 16S rRNA gene amplicons and metaproteomics analysis of fecal samples to determine microbial gut composition and protein expression. 16 rRNA gene sequencing analysis identified Clostridiales, Bacteroidaceae, and Lactobacillaceae species as the most abundant species in the gut. For metaproteomics analysis, peptides were generated by using the Fasp method and subsequently fractionated by strong anion exchanges. Metaproteomics analysis identified 3,673 proteins. Among the most frequently identified proteins, 380 proteins belonged to Lactobacillus spp., 155 belonged to Clostridium spp., and 66 belonged to Streptococcus spp. The most frequently identified proteins were heat shock chaperones, including 349 GroEL proteins, from many bacterial species, whereas the most abundant enzymes were pyruvate kinases, as judged by the number of peptides identified per protein (spectral counting). Gene ontology and KEGG pathway analyses revealed the functions and locations of the identified proteins. The findings of both metaproteomics and 16S rRNA sequencing analyses are discussed.

Synthesis and biological analysis of novel glycoside derivatives of L-AEP, as targeted antibacterial agents

To develop targeted methods for treating bacterial infections, the feasibility of using glycoside derivatives of the antibacterial compound L-R-aminoethylphosphonic acid (L-AEP) has been investigated. These derivatives are hypothesized to be taken up by bacterial cells via carbohydrate uptake mechanisms, and then hydrolysed in situ by bacterial borne glycosidase enzymes, to selectively afford L-AEP. Therefore the synthesis and analysis of ten glycoside derivatives of L-AEP, for selective targeting of specific bacteria, is reported. The ability of these derivatives to inhibit the growth of a panel of Gram-negative bacteria in two different media is discussed. β-Glycosides (12a) and (12b) that contained L-AEP linked to glucose or galactose via a carbamate linkage inhibited growth of a range of organisms with the best MICs being <0.75 mg/ml; for most species the inhibition was closely related to the hydrolysis of the equivalent chromogenic glycosides. This suggests that for (12a) and (12b), release of L-AEP was indeed dependent upon the presence of the respective glycosidase enzyme.

Differential proteome analysis of mature and germinated embryos of Araucaria angustifolia

Araucaria angustifolia is an endangered Brazilian native conifer tree. The aim of the present work was to identify differentially expressed proteins between mature and germinated embryos of A. angustifolia, using one and two dimensional gel electrophoresis approaches followed by protein identification by tandem mass spectrometry. The identities of 32 differentially expressed protein spots from two dimensional gel maps were successfully determined, including proteins and enzymes involved in storage mobilization such as the vicilin-like storage protein and proteases. A label free approach, based on spectral counts, resulted in detection of 10 and 14 mature and germinated enriched proteins, respectively. Identified proteins were mainly related to energetic metabolism pathways, translational processes. oxidative stress regulation and cellular signaling. The integrated use of both strategies permitted a comprehensive protein expression overview of changes in germinated embryos in relation to matures, providing insights into the this process in a recalcitrant seed species. Applications of the data generated on the monitoring and control of in vitro somatic embryos were discussed. Published by Elsevier Ltd.

Bothrops insularis venomics: A proteomic analysis supported by transcriptomic-generated sequence data

A joint transcriptomic and proteomic approach employing two-dimensional electrophoresis, liquid chromatography and mass spectrometry was carried out to identify peptides and proteins expressed by the venom gland of the snake Bothrops insularis, an endemic species of Queimada Grande Island, Brazil. Four protein families were mainly represented in processed spots, namely metalloproteinase, serine proteinase, phospholipase A(2) and lectin. Other represented families were growth factors, the developmental protein G10, a disintegrin and putative novel bradykinin-potentiating peptides. The enzymes were present in several isoforms. Most of the experimental data agreed with predicted values for isoelectric point and M(r) of proteins found in the transcriptome of the venom gland. The results also support the existence of posttranslational modifications and of proteolytic processing of precursor molecules which could lead to diverse multifunctional proteins. This study provides a preliminary reference map for proteins and peptides present in Bothrops insularis whole venom establishing the basis for comparative studies of other venom proteomes which could help the search for new drugs and the improvement of venom therapeutics. Altogether, our data point to the influence of transcriptional and post-translational events on the final venom composition and stress the need for a multivariate approach to snake venomics studies. (c) 2009 Elsevier B.V. All rights reserved.

Cysteine proteases of Trypanosoma (Megatrypanum) theileri: Cathepsin L-like gene sequences as targets for phylogenetic analysis, genotyping diagnosis

Although Trypanosoma theileri and allied trypanosomes are the most widespread trypanosomes in bovids little is known about proteolytic enzymes in these species. We have characterized genes encoding for cathepsin L-like (CATL) cysteine proteases from isolates of cattle, water buffalo and deer that largely diverged from homologues of other trypanosome species. Analysis of 78 CATL catalytic domain sequences from 22 T. theileri trypanosomes disclosed 6 genotypes tightly clustered together into the T. theileri clade. The CATL genes in these trypanosomes are organized in tandem arrays of similar to 1.7 kb located in 2 chromosomal bands of 600-720 kb. A diagnostic PCR assay targeting CATL sequences detected T. theileri of all genotypes from cattle, buffaloes and cervids and also from tabanid vectors. Expression of T. theileri cysteine proteases was demonstrated by proteolytic activity in gelatin gels and hydrolysis of Z-Phe-Arg-AMC substrate. Results from this work agree with previous data using ribosomal and spliced leader genes demonstrating that CATL gene sequences are useful for diagnosis, population genotyping and evolutionary studies of T. theileri trypanosomes. (c) 2010 Elsevier Ireland Ltd. All rights reserved.

Systematic structural studies of iron superoxide dismutases from human parasites and a statistical coupling analysis of metal binding specificity

Superoxide dismutases (SODs) are a crucial class of enzymes in the combat against intracellular free radical damage. They eliminate superoxide radicals by converting them into hydrogen peroxide and oxygen. In spite of their very different life cycles and infection strategies, the human parasites Plasmodium falciparum, Trypanosoma cruzi and Trypanosoma brucei are known to be sensitive to oxidative stress. Thus the parasite Fe-SODs have become attractive targets for novel drug development. Here we report the crystal structures of FeSODs from the trypanosomes T. brucei at 2.0 angstrom and T. cruzi at 1.9 angstrom resolution, and that from P. falciparum at a higher resolution (2.0 angstrom) to that previously reported. The homodimeric enzymes are compared to the related human MnSOD with particular attention to structural aspects which are relevant for drug design. Although the structures possess a very similar overall fold, differences between the enzymes at the entrance to the channel which leads to the active site could be identified. These lead to a slightly broader and more positively charged cavity in the parasite enzymes. Furthermore, a statistical coupling analysis (SCA) for the whole Fe/MnSOD family reveals different patterns of residue coupling for Mn and Fe SODs, as well as for the dimeric and tetrameric states. In both cases, the statistically coupled residues lie adjacent to the conserved core surrounding the metal center and may be expected to be responsible for its fine tuning, leading to metal ion specificity.

The Crystal Complex of Phosphofructokinase-2 of Escherichia coli with Fructose-6-phosphate KINETIC AND STRUCTURAL ANALYSIS OF THE ALLOSTERIC ATP INHIBITION

Substrate inhibition by ATP is a regulatory feature of the phosphofructokinases isoenzymes from Escherichia coli (Pfk-1 and Pfk-2). Under gluconeogenic conditions, the loss of this regulation in Pfk-2 causes substrate cycling of fructose-6-phosphate (fructose-6-P) and futile consumption of ATP delaying growth. In the present work, we have broached the mechanism of ATP-induced inhibition of Pfk-2 from both structural and kinetic perspectives. The crystal structure of Pfk-2 in complex with fructose-6-P is reported to a resolution of 2 angstrom. The comparison of this structure with the previously reported inhibited form of the enzyme suggests a negative interplay between fructose-6-P binding and allosteric binding of MgATP. Initial velocity experiments show a linear increase of the apparent K(0.5) for fructose-6-P and a decrease in the apparent k(cat) as a function of MgATP concentration. These effects occur simultaneously with the induction of a sigmoidal kinetic behavior (n(H) of approximately 2). Differences and resemblances in the patterns of fructose-6-P binding and the mechanism of inhibition are discussed for Pfk-1 and Pfk-2, as an example of evolutionary convergence, because these enzymes do not share a common ancestor.

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R(free)=21.2%) at 1.70 angstrom resolution and 15.81% (R(free) = 19.2%) at 1.85 angstrom resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly beta-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity. (C) 2008 Elsevier Ltd. All rights reserved.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Tuberculosis (TB) is one of the most common infectious diseases known to man and responsible for millions of human deaths in the world. The increasing incidence of TB in developing countries, the proliferation of multidrug resistant strains, and the absence of resources for treatment have highlighted the need of developing new drugs against TB. The shikimate pathway leads to the biosynthesis of chorismate, a precursor of aromatic amino acids. This pathway is absent from mammals and shown to be essential for the survival of Mycobacterium tuberculosis, the causative agent of TB. Accordingly, enzymes of aromatic amino acid biosynthesis pathway represent promising targets for structure-based drug design. The first reaction in phenylalanine biosynthesis involves the conversion of chorismate to prephenate, catalyzed by chorismate mutase. The second reaction is catalyzed by prephenate dehydratase (PDT) and involves decarboxylation and dehydratation of prephenate to form phenylpyruvate, the precursor of phenylalanine. Here, we describe utilization of different techniques to infer the structure of M. tuberculosis PDT (MtbPDT) in solution. Small angle X-ray scattering and ultracentrifugation analysis showed that the protein oligomeric state is a tetramer and MtbPDT is a flat disk protein. Bioinformatics tools were used to infer the structure of MtbPDT A molecular model for MtbPDT is presented and molecular dynamics simulations indicate that MtbPDT i.s stable. Experimental and molecular modeling results were in agreement and provide evidence for a tetrameric state of MtbPDT in solution.