Tailoring Heme-Thiolate Proteins into Efficient Biocatalysts with High Specificity and Selectivity

Cytochrome P450 monooxygenases, one of the most important classes of heme-thiolate proteins, have attracted considerable interest in the biochemical community because of its catalytic versatility, substrate diversity and great number in the superfamily. Although P450s are capable of catalyzing numerous difficult oxidation reactions, the relatively low stability, low turnover rates and the need of electron-donating cofactors have limited their practical biotechnological and pharmaceutical applications as isolated enzymes. The goal of this study is to tailor such heme-thiolate proteins into efficient biocatalysts with high specificity and selectivity by protein engineering and to better understand the structure-function relationship in cytochromes P450. In the effort to engineer P450cam, the prototype member of the P450 superfamily, into an efficient peroxygenase that utilizes hydrogen peroxide via the “peroxide-shunt” pathway, site-directed mutagenesis has been used to elucidate the critical roles of hydrophobic residues in the active site. Various biophysical, biochemical and spectroscopic techniques have been utilized to investigate the wild-type and mutant proteins. Three important P450cam variants were obtained showing distinct structural and functional features. In P450camV247H mutant, which exhibited almost identical spectral properties with the wild-type, it is demonstrated that a single amino acid switch turned the monooxygenase into an efficient preoxidase by increasing the peroxidase activity nearly one thousand folds. In order to tune the distal pocket of P450cam with polar residues, Leu 246 was replaced with a basic residue, lysine, resulting in a mutant with spectral features identical to P420, the inactive species of P450. But this inactive-species-like mutant showed catalytic activities without the facilitation of any cofactors. By substituting Gly 248 with a histidine, a novel Cys-Fe-His ligation set was obtained in P450cam which represented the very rare case of His ligation in heme-thiolate proteins. In addition to serving as a convenient model for hemoprotein structural studies, the G248H mutant also provided evidence about the nature of the axial ligand in cytochrome P420 and other engineered hemoproteins with thiolate ligations. Furthermore, attempts have been made to replace the proximal ligand in sperm whale myoglobin to construct a heme-thiolate protein model by mimicking the protein environment of cytochrome P450cam and chloroperoxidase.

Structural Basis of RNA Recognition by Mycobacterium tuberculosis MazF Homologues

The MazEF toxin-antitoxin (TA) system consists of the antitoxin MazE and the toxin MazF. MazF is a sequence-specific endoribonuclease that upon activation causes cellular growth arrest and increass the level of persisters. Moreover, MazF-induced cells are in a quasi-dormant state that cells remain metabolically active while stop dividing. The quasi-dormancy is similar to the nonreplicating state of M. tuberculosis during latent tuberculosis, thus suggesting the role of mazEF in M. tuberculosis dormancy and persistence. M. tuberculosis has nine mazEF TA modules, each with different RNA cleavage specificities and implicated in selective gene expression during stress conditions. To date only the Bacillus subtilis MazF-RNA complex structure has been determined. As M. tuberculosis MazF homologues recognize distinct RNA sequences, their molecular mechanisms of substrate specificity remain unclear. By taking advantage of X-ray crystallography, we have determined structures of two M. tuberculosis MazF-RNA complexes, MazF-mt1 (Rv2801c) and MazF-mt3 (Rv1991c) in complex with an uncleavable RNA substrate. These structures have provided the molecular basis of sequence-specific RNA recognition and cleavage by MazF toxins.

Both MazF-mt1-RNA and MazF-mt3-RNA complexes showed similar structural organization with one molecule of RNA bound to a MazF-mt1 or MazF-mt3 dimer and occupying the same pocket within the MazF dimer interface. Similar to B. subtilis MazF-RNA complex, MazF-mt1 and MazF-mt3 displayed a conserved active site architecture, where two highly conserved residues, Arg and Thr, form hydrogen bonds with the scissile phosphate group in the cleavage site of the bound RNA. The MazF-mt1-RNA complex also showed specific interactions with its three-base RNA recognition element. Compared with the B. subtilis MazF-RNA complex, our structures showed that residues involved in sequence-specific recognition of target RNA vary between the MazF homologues, therefore explaining the molecular basis for their different RNA recognition sequences. In addition, local conformational changes of the loops in the RNA binding site of MazF-mt1 appear to play a role in MazF targeting different RNA lengths and sequences. In contrast, the MazF-mt3-RNA complex is in a non-optimal RNA binding state with a symmetry-related MazF-mt3 molecule found to make interactions with the bound RNA in the crystal. The crystal-packing interactions were further examined by isothermal titration calorimetry (ITC) studies on selected MazF-mt3 mutants. Our attempts to utilize a MazF-mt3 mutant bearing mutations involved in crystal contacts all crystallized with few nucleotides, which are still found to interact with a symmetry mate. However, these different crystal forms revealed the conformational flexibility of loops in the RNA binding interface of MazF-mt3, suggesting their role in RNA binding and recognition, which will require further studies on additional MazF-mt3-RNA complex interactions.

In conclusion, the structures of the MazF-mt1-RNA and MazF-mt3-RNA complexes provide the first structural information on any M. tuberculosis MazF homologues. Supplemented with structure-guided mutational studies on MazF toxicity in vivo, this study has addressed the structural basis of different RNA cleavage specificities among MazF homologues. Our work will guide future studies on the function of other M. tuberculosis MazF and MazE-MazF homologues, and will help delineate their physiological roles in M. tuberculosis stress responses and pathogenesis.

Development and Characterization of Monovalent and Bivalent RNA Aptamers Targeting the Common Pathway of Coagulation

Anticoagulant agents are commonly used drugs to reduce blood coagulation in acute and chronic clinical settings. Many of these drugs target the common pathway of coagulation because it is critical for thrombin generation and disruption of this portion of the pathway has profound effects on the hemostatic process. Currently available drugs for these indications struggle with balancing desired activity with immunogenicity and poor reversibility or irreversibility in the event of hemorrhage. While improvements are being made with the current drugs, new drugs with better therapeutic indices are needed for surgical intervention and chronic indications to prevent thrombosis from occurring.

A class of therapeutics known as aptamers may be able to meet the need for safer anticoagulant agents. Aptamer are short single-stranded RNA oligonucleotides that adopt specific secondary and tertiary structures based upon their sequence. They can be generated to both enzymes and cofactors because they derive their inhibitory activity by blocking protein-protein interactions, rather than active site inhibition. They inhibit their target proteins with a high level of specificity and bind with high affinity to their target. Additionally, they can be reversed using two different antidote approaches, specific oligonucleotide antidotes, or with cationic, “universal” antidotes. The reversal of their activity is both rapid and durable.

The ability of aptamers to be generated to cofactors has been conclusively proven by generating an aptamer targeting the common pathway coagulation cofactor, Factor V (FV). We developed two aptamers with anticoagulant ability that bind to both FV and FVa, the active cofactor. Both aptamers were truncated to smaller functional sizes and had specific point mutant aptamers developed for use as controls. The anticoagulant activity of both aptamer-mutant pairs was characterized using plasma-based clotting assays and whole blood assays. The mechanism of action resulting in anticoagulant activity was assessed for one aptamer. The aptamer was found to block FVa docking to membrane surfaces, a mechanism not previously observed in any of our other anticoagulant aptamers.

To explore development of aptamers as anticoagulant agents targeting the common pathway for surgical interventions, we fused two anticoagulant aptamers targeting Factor X and prothrombin into a single molecule. The bivalent aptamer was truncated to a minimal size while maintaining robust anticoagulant activity. Characterization of the bivalent aptamer in plasma-based clotting assays indicated we had generated a very robust anticoagulant therapeutic. Furthermore, we were able to simultaneously reverse the activity of both aptamers with a single oligonucleotide antidote. This rapid and complete reversal of anticoagulant activity is not available in the antithrombotic agents currently used in surgery.

Structural Aspects of the Distinct Biochemical Properties of Glutaredoxin 1 and Glutaredoxin 2 from Saccharomyces cerevisiae

Glutaredoxins (Grxs) are small (9-12 kDa) heat-stable proteins that are ubiquitously distributed. In Saccharomyces cerevisiae, seven Grx enzymes have been identified. Two of them (yGrx1 and yGrx2) are dithiolic, possessing a conserved Cys-Pro-Tyr-Cys motif. Here, we show that yGrx2 has a specific activity 15 times higher than that of yGrx1, although these two oxidoreductases share 64% identity and 85% similarity with respect to their amino acid sequences. Further characterization of the enzymatic activities through two-substrate kinetics analysis revealed that yGrx2 possesses a lower Km for glutathione and a higher turnover than yGrx1. To better comprehend these biochemical differences, the pK(a) of the N-terminal active-site cysteines (Cys27) of these two proteins and of the yGrx2-C30S mutant were determined. Since the pK(a) values of the yGrx1 and yGix2 Cys27 residues are very similar, these parameters cannot account for the difference observed between their specific activities. Therefore, crystal structures of yGrx2 in the oxidized form and with a glutathionyl mixed disulfide were determined at resolutions of 2.05 and 1.91 angstrom, respectively. Comparisons of yGrx2 structures with the recently determined structures of yGrx1 provided insights into their remarkable functional divergence. We hypothesize that the substitutions of Ser23 and Gln52 in yGrx1 by Ala23 and Glu52 in yGrx2 modify the capability of the active-site C-terminal cysteine to attack the mixed disulfide between the N-terminal active-site cysteine and the glutathione molecule. Mutagenesis studies supported this hypothesis. The observed structural and functional differences between yGrx1 and yGrx2 may reflect variations in substrate specificity. (C) 2008 Elsevier Ltd. All rights reserved.

S-nitrosylation homeostasis in plants: key enzymes and regulatory pathways

The nitrosylated form of glutathione (GSNO) has been acknowledged to be the most important nitrosylating agent of the plant cell, and the tuning of its intracellular concentration is of pivotal importance for photosynthetic life. During my time as a PhD student, I focused my attention on the enzymatic systems involved in the degradation of GSNO. Hence, we decided to study the structural and catalytic features of alcohol dehydrogenases (GSNOR and ADH1) from the model land plant Arabidopsis thaliana (At). These enzymes displayed a very similar 3D structure except for their active site which might explain the extreme catalytic specialization of the two enzymes. They share NAD(H) as a cofactor, but only AtGSNOR was able to catalyze the reduction of GSNO whilst being ineffective in oxidizing ethanol. Moreover, our study on the enzyme from the unicellular green alga Chlamydomonas reinhardtii (Cr) revealed how this S-nitrosoglutathione reductase (GSNOR) specifically use NADH to catalyze GSNO reduction and how its activity responds to thiol-based post-translational modifications. Contextually, the presence of NADPH-dependent GSNO-degrading systems in algal protein extract was highlighted and resulted to be relatively efficient in this model organism. This activity could be ascribed to several proteins whose contribution has not been defined yet. Intriguingly, protein extract from GSNOR null mutants of Arabidopsis displayed an increased NADPH-dependent ability to degrade GSNO and our quantitative proteome profiling on the gsnor mutant revealed the overexpression of two class 4 aldo-keto reductases (AKR), specifically AtAKR4C8 and AtAKR4C9. Later, all four class 4 AKRs showed to possess a NADPH-dependent GSNO-degrading activity. Finally, we initiated a preliminary analysis to determine the kinetic parameters of several plant proteins, including GSNOR, AKR4Cs, and thioredoxins. These data suggested GSNOR to be the most effective enzyme in catalyzing GSNO reduction because of its extremely high catalytic proficiency compared to NADPH-dependent systems.

Molecular And Biochemical Characterization Of Caffeine Synthase And Purine Alkaloid Concentration In Guarana Fruit.

Guarana seeds have the highest caffeine concentration among plants accumulating purine alkaloids, but in contrast with coffee and tea, practically nothing is known about caffeine metabolism in this Amazonian plant. In this study, the levels of purine alkaloids in tissues of five guarana cultivars were determined. Theobromine was the main alkaloid that accumulated in leaves, stems, inflorescences and pericarps of fruit, while caffeine accumulated in the seeds and reached levels from 3.3% to 5.8%. In all tissues analysed, the alkaloid concentration, whether theobromine or caffeine, was higher in young/immature tissues, then decreasing with plant development/maturation. Caffeine synthase activity was highest in seeds of immature fruit. A nucleotide sequence (PcCS) was assembled with sequences retrieved from the EST database REALGENE using sequences of caffeine synthase from coffee and tea, whose expression was also highest in seeds from immature fruit. The PcCS has 1083bp and the protein sequence has greater similarity and identity with the caffeine synthase from cocoa (BTS1) and tea (TCS1). A recombinant PcCS allowed functional characterization of the enzyme as a bifunctional CS, able to catalyse the methylation of 7-methylxanthine to theobromine (3,7-dimethylxanthine), and theobromine to caffeine (1,3,7-trimethylxanthine), respectively. Among several substrates tested, PcCS showed higher affinity for theobromine, differing from all other caffeine synthases described so far, which have higher affinity for paraxanthine. When compared to previous knowledge on the protein structure of coffee caffeine synthase, the unique substrate affinity of PcCS is probably explained by the amino acid residues found in the active site of the predicted protein.

Synthesis And Antiproliferative Activity Of 8-hydroxyquinoline Derivatives Containing A 1,2,3-triazole Moiety.

Twelve novel 8-hydroxyquinoline derivatives were synthesized with good yields by performing copper-catalyzed Huisgen 1,3-dipolar cycloaddition (click reaction) between an 8-O-alkylated-quinoline containing a terminal alkyne and various aromatic or protected sugar azides. These compounds were evaluated in vitro for their antiproliferative activity on various cancer cell types. Protected sugar derivative 16 was the most active compound in the series, exhibiting potent antiproliferative activity and high selectivity toward ovarian cancer cells (OVCAR-03, GI50 < 0.25 μg mL(-1)); this derivative was more active than the reference drug doxorubicin (OVCAR-03, GI50 = 0.43 μg mL(-1)). In structure-activity relationship (SAR) studies, the physico-chemical parameters of the compounds were evaluated and docking calculations were performed for the α-glucosidase active site to predict the possible mechanism of action of this series of compounds.

Digestive Peptidase Evolution In Holometabolous Insects Led To A Divergent Group Of Enzymes In Lepidoptera.

Trypsins and chymotrypsins are well-studied serine peptidases that cleave peptide bonds at the carboxyl side of basic and hydrophobic l-amino acids, respectively. These enzymes are largely responsible for the digestion of proteins. Three primary processes regulate the activity of these peptidases: secretion, precursor (zymogen) activation and substrate-binding site recognition. Here, we present a detailed phylogenetic analysis of trypsins and chymotrypsins in three orders of holometabolous insects and reveal divergent characteristics of Lepidoptera enzymes in comparison with those of Coleoptera and Diptera. In particular, trypsin subsite S1 was more hydrophilic in Lepidoptera than in Coleoptera and Diptera, whereas subsites S2-S4 were more hydrophobic, suggesting different substrate preferences. Furthermore, Lepidoptera displayed a lineage-specific trypsin group belonging only to the Noctuidae family. Evidence for facilitated trypsin auto-activation events were also observed in all the insect orders studied, with the characteristic zymogen activation motif complementary to the trypsin active site. In contrast, insect chymotrypsins did not seem to have a peculiar evolutionary history with respect to their mammal counterparts. Overall, our findings suggest that the need for fast digestion allowed holometabolous insects to evolve divergent groups of peptidases with high auto-activation rates, and highlight that the evolution of trypsins led to a most diverse group of enzymes in Lepidoptera.

Synthesis, screening for antiacetylcholinesterase activity and binding mode prediction of a new series of [3-(disubstituted-phosphate)-4,4,4-trifluoro-butyl]-carbamic acid ethyl esters

A series of nine new [3-(disubstituted-phosphate)-4,4,4-trifluoro-butyl]-carbamic acid ethyl esters (phosphate-carbamate compounds) was obtained through the reaction of (4,4,4-trifluoro-3-hydroxybut-1-yl)-carbamic acid ethyl esters with phosphorus oxychloride followed by the addition of alcohols. The products were characterized by ¹H, 13C, 31P, and 19F NMR spectroscopy, GC-MS, and elemental analysis. All the synthesized compounds were screened for acetylcholinesterase (AChE) inhibitory activity using the Ellman method. All compounds containing phosphate and carbamate pharmacophores in their structures showed enzyme inhibition, being the compound bearing the diethoxy phosphate group (2b) the most active compound. Molecular modeling studies were performed to investigate the detailed interactions between AChE active site and small-molecule inhibitor candidates, providing valuable structural insights into AChE inhibition.

Hypervalent organochalcogenanes as inhibitors of protein tyrosine phosphatases

A series of organochalcogenanes was synthesized and evaluated as protein tyrosine phosphatases (PTPs) inhibitors. The results indicate that organochalcogenanes inactivate the PTPs in a time- and concentration-dependent fashion, most likely through covalent modification of the active site sulfur-moiety by the chalcogen atom. Consequently, organochalcogenanes represent a new class of mechanism-based probes to modulate the PTP-mediated cellular processes.

Electrochemical and spectroscopic characterization of screen-printed gold-based electrodes modified with self-assembled monolayers and Tc85 protein

This work investigates the formation of self-assembled monolayers (SAMs) of cystamine and cystamine-glutaraldehyde on a screen-printed electrode, and the immobilization of the Tc85 protein (from Trypanosoma cruzi) on these monolayers. The methods used included infrared techniques, cyclic voltammetry, and electrochemical impedance spectroscopy. The electrochemical studies were performed at pH 6.9 in 0.1 mol L(-1) phosphate buffer solution containing Fe(CN)(6)(-3/-4) redox species. The surface coverage (0) of the electrode was 0.10 (cystamine), 0.35 (cystamine-glutaraldehyde) and 0.84 (Tc85). Interpretation of electrochemical impedance spectroscopy results was based on a charge-transfer reaction involving Fe(CN)(6)(-3/-4) species at high frequencies, followed by a diffusion through the monolayers at lower frequencies. Estimates of the electrode surface coverage, active site radius, and distance between two adjacent sites assumed that charge transfer occurred at the active sites, and that there was a planar diffusion of redox species to these sites. (C) 2009 Elsevier B.V. All rights reserved.

3D-Pharmacophore mapping of thymidine-based inhibitors of TMPK as potential antituberculosis agents

Tuberculosis (TB) is the primary cause of mortality among infectious diseases. Mycobacterium tuberculosis monophosphate kinase (TMPKmt) is essential to DNA replication. Thus, this enzyme represents a promising target for developing new drugs against TB. In the present study, the receptor-independent, RI, 4D-QSAR method has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 81 thymidine analogues, and two corresponding subsets, reported as inhibitors of TMPKmt. The resulting optimized models are not only statistically significant with r (2) ranging from 0.83 to 0.92 and q (2) from 0.78 to 0.88, but also are robustly predictive based on test set predictions. The most and the least potent inhibitors in their respective postulated active conformations, derived from each of the models, were docked in the active site of the TMPKmt crystal structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. Moreover, the QSAR models provide insights regarding a probable mechanism of action of the analogues.

Rational Design and 3D-Pharmacophore Mapping of 5 `-Thiourea-Substituted alpha-Thymidine Analogues as Mycobacterial TMPK Inhibitors

Thymidine monophosphate kinase (TMPK) has emerged as an attractive target for developing inhibitors of Mycobacterium tuberculosis growth. In this study the receptor-independent (RI) 4D-QSAR formalism has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 5`-thiourea-substituted alpha-thymidine inhibitors. Models were developed for the entire training set and for a subset of the training set consisting of the most potent inhibitors. The optimized (RI) 4D-QSAR models are statistically significant (r(2) = 0.90, q(2) = 0.83 entire set, r(2) = 0.86, q(2) = 0.80 high potency subset) and also possess good predictivity based on test set predictions. The most and least potent inhibitors, in their respective postulated active conformations derived from the models, were docked in the active site of the TMPK crystallographic structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. This model identifies new regions of the inhibitors that contain pharmacophore sites, such as the sugar-pyrimidine ring structure and the region of the 5`-arylthiourea moiety. These new regions of the ligands can be further explored and possibly exploited to identify new, novel, and, perhaps, better antituberculosis inhibitors of TMPKmt. Furthermore, the 3D-pharmacophores defined by these models can be used as a starting point for future receptor-dependent antituberculosis drug design as well as to elucidate candidate sites for substituent addition to optimize ADMET properties of analog inhibitors.

Structural complexes of human adenine phosphoribosyltransferase reveal novel features of the APRT catalytic mechanism

Adenine phosphoribosyltransferase (APRT) is an important enzyme component of the purine recycling pathway. Parasitic protozoa of the order Kinetoplastida are unable to synthesize purines de novo and use the salvage pathway for the synthesis of purine bases rendering this biosynthetic pathway an attractive target for antiparasitic drug design. The recombinant human adenine phosphoribosyltransferase (hAPRT) structure was resolved in the presence of AMP in the active site to 1.76 angstrom resolution and with the substrates PRPP and adenine simultaneously bound to the catalytic site to 1.83 angstrom resolution. An additional structure was solved containing one subunit of the dimer in the apo-form to 2.10 angstrom resolution. Comparisons of these three hAPRT structures with other `type I` PRTases revealed several important features of this class of enzymes. Our data indicate that the flexible loop structure adopts an open conformation before and after binding of both substrates adenine and PRPR Comparative analyses presented here provide structural evidence to propose the role of Glu 104 as the residue that abstracts the proton of adenine N9 atom before its nucleophilic attack on the PRPP anomeric carbon. This work leads to new insights to the understanding of the APRT catalytic mechanism.

Structural and functional characterisation of human sulfotransferases

The human aryl sulfotransferases HAST4 and HAST4v vary by only two amino acids but exhibit markedly different affinity towards the sulfonate acceptor p-nitrophenol and the sulfonate donor 3'-phosphoadenosine-5'-phosphosulfate (PAPS). To determine the importance of each of these amino acid differences, chimeric constructs were made of HAST4 and HAST4v. By attaching the last 120 amino acids of HAST-4v to HAST4 (changing Thr235 to Asn235) we have been able to produce a protein that has a K-m for PAPS similar to HAST4v. The reverse construct, HAST4v/4 produces a protein with a K-m for PAPS similar to HAST4. These data suggests that the COOH-terminal of sulfotransferases is involved in co-factor binding. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.