Computer-Aided Saturation Mutagenesis of Arabidopsis thaliana Ent-Copalyl Diphosphate Synthase.

Ent-copalyl diphosphate synthase controls the biosynthesis of gibberellin plant hormones, which in turn coordinate the expression of numerous enzymes. Some gibberellin-dependent genes encode enzymes coordinating the biosynthesis of tanshinones: diterpene derivatives with broad medical applications. New biotechnological approaches, such as metabolic engineering using naturally occurring or mutated enzymes, have been proposed to meet the growing demand for tanshinones which is currently met by the Chinese medicinal plant Salvia miltiorrhiza Bunge. These mutants may be prepared by directed evolution, saturation mutagenesis or rational enzyme design. In the presented paper, 15,257 non-synonymous variants of Arabidopsis thaliana ent-copalyl diphosphate synthase were obtained using the SNAP2 tool. The obtained forms were screened to isolate variants with potentially improved biological functions. A group of 455 mutants with potentially improved stability was isolated and subjected to further screening on the basis of ligand-substrate affinity, and both secondary structure and active site structure stability. Finally, a group of six single mutants was obtained, which were used to construct double mutants with potentially improved stability and ligand affinity. The potential influence of single mutations on protein stability and ligand affinity was evaluated by double mutant cycle analysis. Finally, the procedure was validated by in silico assessment of the experimentally verified enzyme mutants with reduced enzymatic activity.

Structure of a Kunitz-type potato cathepsin D inhibitor.

Potato cathepsin D inhibitor (PDI) is a glycoprotein of 188 amino acids which can inhibit both the aspartic protease cathepsin D and the serine protease trypsin. Here we report the first X-ray structure of PDI at a resolution of 2.1 Å showing that PDI adopts a ß-trefoil fold, which is typical of the Kunitz-family protease inhibitors, with the inhibitory loops protruding from the core. Possible reactive-site loops including one involving a unique disulphide and another involving a protruding 310 helix are identified and docking studies indicate the mode of action of this unusual bi-functional inhibitor.

Enzyme renaturation to higher activity driven by the sol-gel transition: Carbonic anhydrase.

We describe a so-far unknown route for renaturing denatured enzymes, namely subjecting the denatured enzyme to an oxide sol-gel transition. The phenomenon was revealed in a detailed study of denatured carbonic anhydrase which was subjected to an alumina sol-gel transition, up to the thermally stabilizing entrapment in the final xerogel. Remarkably, not only that the killed enzyme regained its activity during the sol-gel process, but its activity increased to 180% of the native enzyme. To the best of our knowledge, this is the first report of enhanced activity following by renaturing (a "Phoenix effect"). Kinetic study which revealed a five-orders of magnitude (!) increase in the Arrhenius prefactor upon entrapment compared to solution. Circular dichroism analysis, differential scanning calorimetry, zeta potential analyses as well as synchronous fluorescence measurements, all of which were used to characterize the phenomenon, are consistent with a proposed mechanism which is based on the specific orienting interactions of the active site of the enzyme with respect to the alumina interface and its pores network.

Computational redesign of metalloenzymes for catalyzing new reactions.

The ability to design novel activities in existing metalloenzyme active sites is a stringent test of our understanding of enzyme mechanisms, sheds light on enzyme evolution, and would have many practical applications. Here, we describe a computational method in the context of the macromolecular modeling suite Rosetta to repurpose active sites containing metal ions for reactions of choice. The required inputs for the method are a model of the transition state(s) for the reaction and a set of crystallographic structures of proteins containing metal ions. The coordination geometry associated with the metal ion (Zn(2+), for example) is automatically detected and the transition state model is aligned to the open metal coordination site(s) in the protein. Additional interactions to the transition state model are made using RosettaMatch and the surrounding amino acid side chain identities are optimized for transition state stabilization using RosettaDesign. Validation of the design is performed using docking and molecular dynamics simulations, and candidate designs are generated for experimental validation. Computational metalloenzyme repurposing is complementary to directed evolution approaches for enzyme engineering and allows large jumps in sequence space to make concerted sequence and structural changes for introducing novel enzymatic activities and specificities.

Cloning, characterization and evaluation of potent inhibitors of Shigella sonnei acetohydroxyacid synthase catalytic subunit.

Acetohydroxyacid synthase (AHAS) is a thiamin diphosphate (ThDP)- and flavin adenine dinucleotide (FAD)-dependent plant and microbial enzyme that catalyzes the first common step in the biosynthesis of essential amino acids such as leucine, isoleucine and valine. To identify strong potent inhibitors against Shigella sonnei (S. sonnei) AHAS, we cloned and characterized the catalytic subunit of S. sonnei AHAS and found two potent chemicals (KHG20612, KHG25240) that inhibit 87-93% S. sonnei AHAS activity at an inhibitor concentration of 100uM. The purified S. sonnei AHAS had a size of 65kDa on SDS-PAGE. The enzyme kinetics revealed that the enzyme has a K(m) of 8.01mM and a specific activity of 0.117U/mg. The cofactor activation constant (K(s)) for ThDP and (K(c)) for Mg(++) were 0.01mM and 0.18mM, respectively. The dissociation constant (K(d)) for ThDP was found to be 0.14mM by tryptophan fluorescence quenching. The inhibition kinetics of inhibitor KHG20612 revealed an un-competitive inhibition mode with a K(ii) of 2.65mM and an IC(50) of 9.3µM, whereas KHG25240 was a non-competitive inhibitor with a K(ii of) 5.2mM, K(is) of 1.62mM and an IC(50) of 12.1µM. Based on the S. sonnei AHAS homology model structure, the docking of inhibitor KHG20612 is predicted to occur through hydrogen bonding with Met 257 at a 1.7Å distance with a low negative binding energy of -9.8kcal/mol. This current study provides an impetus for the development of a novel strong antibacterial agent targeting AHAS based on these potent inhibitor scaffolds.

Catalytic subunits atpα and atpß from the Pacific white shrimp Litopenaeus vannamei F(O)F (1) ATP-synthase complex: cDNA sequences, phylogenies, and mRNA quantification during hypoxia.

In the mitochondrial F(O)F(1) ATP-synthase/ATPase complex, subunits α and ß are part of the extrinsic portion that catalyses ATP synthesis. Since there are no reports about genes and proteins from these subunits in crustaceans, we analyzed the cDNA sequences of both subunits in the whiteleg shrimp Litopenaeus vannamei and their phylogenetic relationships. We also investigated the effect of hypoxia on shrimp by measuring changes in the mRNA amounts of atpα and atpß. Our results confirmed highly conserved regions for both subunits and underlined unique features among others. The ATPß deduced protein of shrimp was less conserved in size and sequence than ATPα. The relative mRNA amounts of atpα and atpß changed in shrimp pleopods; hypoxia at 1.5 mg/L caused an increase in atpß transcripts and a subsequent decrease when shrimp were re-oxygenated. Results confirm that changes in the mRNAs of the ATP-synthase subunits are part of the mechanisms allowing shrimp to deal with the metabolic adjustment displayed to tolerate hypoxia.

Prediction of the interaction of HIV-1 integrase and its dicaffeoylquinic acid inhibitor through molecular modeling approach.

An essential step in the life cycle of human immunodeficiency virus type 1 (HIV-1) is integration of the double-stranded retroviral DNA into the genome of the host cell. HIV-1 integrase, the enzyme that inserts the vital DNA into the host chromosome, is an attractive and rational target for anti-AIDS drug design because it is essential for HIV replication and there are no known counterparts in the host cell. Inhibitors of this enzyme have a great potential to complement the therapeutic use of HIV protease and reverse transcriptase inhibitors. Natural products have provided a source of new drug candidates for anti-AIDS therapy. Dicaffeoylquinic acids, isolated from traditional medicinal plants, are a novel class of integrase inhibitors. These compounds are potent inhibitors of HIV-1 replication in cultured cell lines and catalytic activities of integrase in vitro. They are therefore promising compounds for developing new anti-AIDS drugs. To understand how the inhibitors work and therefore design more potent and specific inhibitors, we have used molecular modeling techniques to investigate the binding modes of 3,4-dicaffeoylquinic acid. Our computational modeling study demonstrated that the inhibitor of this compound on HIV integrase is likely to proceed by two different but equivalent mechanisms with one bound to the active site region of the enzyme and another docked into the binding pocket located on the other side of the catalytic site. Our study will be of help to design new pharmaceuticals for the treatment of AIDS.

Conformational plasticity of the catalytic subunit of protein kinase CK2 and its consequences for regulation and drug design.

At the first glance CK2alpha, the catalytic subunit of protein kinase CK2, is a rigid molecule: in contrast to many eukaryotic protein kinases in CK2alpha the canonical regulatory key elements like the activation segment occur exclusively in their typical active conformations. This observation fits well to the constitutive activity of the enzyme, meaning, its independence from phosphorylation or other characteristic control factors. Most CK2alpha structures are based on the enzyme from Zea mays, supplemented by an increasing number of human CK2alpha structures. In the latter a surprising plasticity of important ATP-binding elements - the interdomain hinge region and the glycine-rich loop - was discovered. In fully active CK2alpha the hinge region is open and does not anchor the ATP ribose, but alternatively it can adopt a closed conformation, form hydrogen bonds to the ribose moiety and thus retract the gamma-phospho group from its functional position. In addition to this partially inactive state human CK2alpha was recently found in a fully inactive conformation. It is incompatible with ATP-binding due to a combination of a closed hinge and a collapse of the glycine-rich loop into the ATP cavity. These conformational transitions are apparently correlated with the occupation state of a remote docking site located at the interface to the non-catalytic subunit CK2beta: if CK2beta blocks this site, the fully active conformation of CK2alpha is stabilized, while the binding of certain small molecule seems to favour the partially and fully inactive states. This observation may be exploited to design effective and selective CK2 inhibitors.

[The regulation of peroxisomal matrix enzymes (alcohol oxidase and catalase) formation by the product of the gene Mth1 in methylotrophic yeast Pichia methanolica].

Two independent mutant strains of methylotrophic yeast Pichia methanolica (mth1 arg1 and mth2 arg4) from the initial line 616 (ade1 ade5) were investigated. The mutant strains possessed defects in genes MTH1 and MTH2 which resulted in the inability to assimilate methanol as a sole carbon source and the increased activity of alcohol oxidase (AO). The function of the AUG2 gene encoding one of the subunits of AO and CTA1, a probable homolog of peroxisomal catalase of Saccharomyces cereviseae, was investigated by analyses of the molecular forms of isoenzymes. It was shown that optimal conditions for the expression of the AUG2 gene on a medium supplemented with 3% of methanol leads to an increasing synthesis of peroxisomal catalase. The mutant mth1 possessed a dominant formation of AO isoform with electrophoretic mobility which is typical for isogenic form 9, the product of the AUG2 gene, and a decreased level of peroxisomal catalase. The restoration of growth of four spontaneous revertants of the mutant mth1 (Rmth1) on the methanol containing medium was accompanied by an increase in activity of AO isogenic form 9 and peroxisomal catalase. The obtained results confirmed the functional continuity of the structural gene AUG2 in mutant mth1. The correlation of activity of peroxisomal catalase and AO isogenic form 1 in different conditions evidenced the existence of common regulatory elements for genes AUG2 and CTA1 in methilotrophic yeast Pichia methanolica.

Identification of in vitro phosphorylation sites in the Arabidopsis thaliana somatic embryogenesis receptor-like kinases.

The Arabidopsis thaliana somatic embryogenesis receptor-like kinase (SERK) family consists of five leucine-rich repeat receptor-like kinases (LRR-RLKs) with diverse functions such as brassinosteroid insensitive 1 (BRI1)-mediated brassinosteroid perception, development and innate immunity. The autophosphorylation activity of the kinase domains of the five SERK proteins was compared and the phosphorylated residues were identified by LC-MS/MS. Differences in autophosphorylation that ranged from high activity of SERK1, intermediate activities for SERK2 and SERK3 to low activity for SERK5 were noted. In the SERK1 kinase the C-terminally located residue Ser-562 controls full autophosphorylation activity. Activation loop phosphorylation, including that of residue Thr-462 previously shown to be required for SERK1 kinase activity, was not affected. In vivo SERK1 phosphorylation was induced by brassinosteroids. Immunoprecipitation of CFP-tagged SERK1 from plant extracts followed by MS/MS identified Ser-303, Thr-337, Thr-459, Thr-462, Thr-463, Thr-468, and Ser-612 or Thr-613 or Tyr-614 as in vivo phosphorylation sites of SERK1. Transphosphorylation of SERK1 by the kinase domain of the main brassinosteroid receptor BRI1 occurred only on Ser-299 and Thr-462. This suggests both intra- and intermolecular control of SERK1 kinase activity. Conversely, BRI1 was transphosphorylated by the kinase domain of SERK1 on Ser-887. BRI1 kinase activity was not required for interaction with the SERK1 receptor in a pull down assay.

Unconventional serine proteases: variations on the catalytic Ser/His/Asp triad configuration.

Serine proteases comprise nearly one-third of all known proteases identified to date and play crucial roles in a wide variety of cellular as well as extracellular functions, including the process of blood clotting, protein digestion, cell signaling, inflammation, and protein processing. Their hallmark is that they contain the so-called "classical" catalytic Ser/His/Asp triad. Although the classical serine proteases are the most widespread in nature, there exist a variety of "nonclassical" serine proteases where variations to the catalytic triad are observed. Such variations include the triads Ser/His/Glu, Ser/His/His, and Ser/Glu/Asp, and include the dyads Ser/Lys and Ser/His. Other variations are seen with certain serine and threonine peptidases of the Ntn hydrolase superfamily that carry out catalysis with a single active site residue. This work discusses the structure and function of these novel serine proteases and threonine proteases and how their catalytic machinery differs from the prototypic serine protease class.

Expression of glucose-6-phosphatase system genes in murine cortex and hypothalamus.

The glucose-6-phosphatase (G6Pase) system participates in the regulation of glucose homeostasis by converting glucose-6-phosphate (G6P) into glucose and inorganic phosphates. We have used an RT-PCR-based cloning and sequencing approach to study the expression of components of the G6Pase system in the hypothalamus and cortex tissues of the ob/ob mouse. We observed the expression of hepatic G6Pase catalytic subunit, G6PC, in both tissues, although increased template inputs were required for its detection. Conversely, expression of both the mouse homologue of the previously-described brain-specific G6P translocase T1 (G6PT1) variant and of the hepatic G6PT1 isoform was easily detectable in hypothalamus and cortex tissues. Of the proposed G6Pase catalytic subunit homologues, the expression of murine ubiquitous G6Pase catalytic subunit-related protein (UGRP, G6PC3) was also easily detectable in both tissues. However, islet-specific G6Pase catalytic subunit-related protein (IGRP, G6PC2) was expressed in a tissue-specific manner, and was detectable only in hypothalamus tissue at increased template inputs. We conclude that cells within ob/ob mouse hypothalamus and cortex tissues express genes with either established or proposed roles in G6P hydrolysis.

Functional domain organization of the potato alpha-glucan, water dikinase (GWD): evidence for separate site catalysis as revealed by limited proteolysis and deletion mutants.

The potato tuber (Solanum tuberosum) GWD (alpha-glucan, water dikinase) catalyses the phosphorylation of starch by a dikinase-type reaction mechanism in which the beta-phosphate of ATP is transferred to the glucosyl residue of amylopectin. GWD shows sequence similarity to bacterial pyruvate, water dikinase and PPDK (pyruvate, phosphate dikinase). In the present study, we examine the structure-function relationship of GWD. Analysis of proteolytic fragments of GWD, in conjunction with peptide microsequencing and the generation of deletion mutants, indicates that GWD is comprised of five discrete domains of 37, 24, 21, 36 and 38 kDa. The catalytic histidine, which mediates the phosphoryl group transfer from ATP to starch, is located on the 36 kDa fragment, whereas the 38 kDa C-terminal fragment contains the ATP-binding site. Binding of the glucan molecule appears to be confined to regions containing the three N-terminal domains. Deletion mutants were generated to investigate the functional interdependency of the putative ATP- and glucan-binding domains. A truncated form of GWD expressing the 36 and 38 kDa C-terminal domains was found to catalyse the E+ATP-->E-P+AMP+P(i) (where P(i) stands for orthophosphate) partial reaction, but not the E-P+glucan-->E+glucan-P partial reaction. CD experiments provided evidence for large structural changes on autophosphorylation of GWD, indicating that GWD employs a swivelling-domain mechanism for enzymic phosphotransfer similar to that seen for PPDK.

Na,K-ATPase alpha3 subunit in the goldfish retina during optic nerve regeneration.

The goldfish optic nerve can regenerate after injury. To understand the molecular mechanism of optic nerve regrowth, we identified genes whose expression is specifically up-regulated during the early stage of optic nerve regeneration. A cDNA library constructed from goldfish retina 5 days after transection was screened by differential hybridization with cDNA probes derived from axotomized or normal retina. Of six cDNA clones isolated, one clone was identified as the Na,K-ATPase catalytic subunit alpha3 isoform by high- sequence homology. In northern hybridization, the expression level of the mRNA was significantly increased at 2 days and peaked at 5-10 days, and then gradually decreased and returned to control level by 45 days after optic nerve transection. Both in situ hybridization and immunohistochemical staining have revealed the location of this transient retinal change after optic nerve transection. The increased expression was observed only in the ganglion cell layer and optic nerve fiber layer at 5-20 days after optic nerve transection. In an explant culture system, neurite outgrowth from the retina 7 days after optic nerve transection was spontaneously promoted. A low concentration of ouabain (50-100 nm ) completely blocked the spontaneous neurite outgrowth from the lesioned retina. Together, these data indicate that up-regulation of the Na,K-ATPase alpha3 subunit is involved in the regrowth of ganglion cell axons after axotomy.

Neutral inhibitors of the serine protease factor Xa.

A neutral inhibitor of the serine protease factor Xa was identified via a high-throughput screen of a commercial library. The initial lead 1 demonstrated reversible and competitive inhibition kinetics for factor Xa and possessed a high degree of selectivity versus other related serine proteases. Initial modeling efforts and the generation of a series of analogues of 1 are described.

Secreted phospholipase A(2) induces vascular endothelial cell migration.

Secreted phospholipase A(2) (sPLA(2)) regulates a variety of cellular functions. The present investigation was undertaken to elucidate the potential role of sPLA(2) in endothelial cell (EC) migration. Bovine aortic endothelial cells (BAECs) exposed to sPLA(2) placed in the lower compartment of a modified Boyden chamber displayed increased migration compared to cells exposed to vehicle. The effect of sPLA(2) on EC migration was time and dose dependent. Migration of BAECs was observed at 30 minutes, increased over 1 to 2 hours, and declined thereafter. At 2 hours of stimulation, sPLA(2) (0.01-2 micromol/L) induced 1.2- to 3-fold increased cell migration compared with media alone. Among the different sPLA(2)s tested, bee venom, Naja naja, and porcine and human pancreatic PLA(2)s all evoked a migratory response in ECs. Moreover, human synovial fluid, obtained from patients with arthritis and containing sPLA(2) activity, induced EC migration. Migration of ECs was significantly reduced after exposure to a catalytic site mutant of pancreatic sPLA(2) with decreased lipolytic activity as compared to wild-type sPLA(2). Similarly, pretreatment of human synovial fluid with p-bromophenacyl bromide, an irreversible inhibitor of sPLA(2), markedly decreased the ability of human synovial fluid to stimulate EC migration. Moreover, migration of ECs was stimulated on exposure to hydrolytic products of sPLA(2) activity including arachidonic acid, lysophosphatidic acid, and lysophosphatidylcholine. These findings suggest that sPLA(2) plays a physiologic role in induction of EC migration. Moreover, the effects of sPLA(2) on EC migration are mediated, at least in part, by its catalytic activity. (Blood. 2000;96:3809-3815)