Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.

Methyl Jasmonate Induces Lauric Acid ω-Hydroxylase Activity and Accumulation of CYP94A1 Transcripts but Does Not Affect Epoxide Hydrolase Activities in Vicia sativa Seedlings1

Treatment of etiolated Vicia sativa seedlings by the plant hormone methyl jasmonate (MetJA) led to an increase of cytochrome P450 content. Seedlings that were treated for 48 h in a 1 mm solution of MetJA stimulated ω-hydroxylation of 12:0 (lauric acid) 14-fold compared with the control (153 versus 11 pmol min−1 mg−1 protein, respectively). Induction was dose dependent. The increase of activity (2.7-fold) was already detectable after 3 h of treatment. Activity increased as a function of time and reached a steady level after 24 h. Northern-blot analysis revealed that the transcripts coding for CYP94A1, a fatty acid ω-hydroxylase, had already accumulated after 1 h of exposure to MetJA and was maximal between 3 and 6 h. Under the same conditions, a study of the enzymatic hydrolysis of 9,10-epoxystearic acid showed that both microsomal and soluble epoxide hydrolase activities were not affected by MetJA treatment.

Designing safer chemicals: predicting the rates of metabolism of halogenated alkanes.

A computational model is presented that can be used as a tool in the design of safer chemicals. This model predicts the rate of hydrogen-atom abstraction by cytochrome P450 enzymes. Excellent correlations between biotransformation rates and the calculated activation energies (delta Hact) of the cytochrome P450-mediated hydrogen-atom abstractions were obtained for the in vitro biotransformation of six halogenated alkanes (1-fluoro-1,1,2,2-tetrachloroethane, 1,1-difluoro-1,2,2-trichloroethane, 1,1,1-trifluro-2,2-dichloroethane, 1,1,1,2-tetrafluoro-2-chloroethane, 1,1,1,2,2,-pentafluoroethane, and 2-bromo-2-chloro-1,1,1-trifluoroethane) with both rat and human enzyme preparations: In(rate, rat liver microsomes) = 44.99 - 1.79(delta Hact), r2 = 0.86; In(rate, human CYP2E1) = 46.99 - 1.77(delta Hact), r2 = 0.97 (rates are in nmol of product per min per nmol of cytochrome P450 and energies are in kcal/mol). Correlations were also obtained for five inhalation anesthetics (enflurane, sevoflurane, desflurane, methoxyflurane, and isoflurane) for both in vivo and in vitro metabolism by humans: In[F(-)]peak plasma = 42.87 - 1.57(delta Hact), r2 = 0.86. To our knowledge, these are the first in vivo human metabolic rates to be quantitatively predicted. Furthermore, this is one of the first examples where computational predictions and in vivo and in vitro data have been shown to agree in any species. The model presented herein provides an archetype for the methodology that may be used in the future design of safer chemicals, particularly hydrochlorofluorocarbons and inhalation anesthetics.

Benzoic acid 2-hydroxylase, a soluble oxygenase from tobacco, catalyzes salicylic acid biosynthesis.

Benzoic acid 2-hydroxylase (BA2H) catalyzes the biosynthesis of salicylic acid from benzoic acid. The enzyme has been partially purified and characterized as a soluble protein of 160 kDa. High-efficiency in vivo labeling of salicylic acid with 18O2 suggested that BA2H is an oxygenase that specifically hydroxylates the ortho position of benzoic acid. The enzyme was strongly induced by either tobacco mosaic virus inoculation or benzoic acid infiltration of tobacco leaves and it was inhibited by CO and other inhibitors of cytochrome P450 hydroxylases. The BA2H activity was immunodepleted by antibodies raised against SU2, a soluble cytochrome P450 from Streptomyces griseolus. The anti-SU2 antibodies immunoprecipitated a radiolabeled polypeptide of around 160 kDa from the soluble protein extracts of L-[35S]-methionine-fed tobacco leaves. Purified BA2H showed CO-difference spectra with a maximum at 457 nm. These data suggest that BA2H belongs to a novel class of soluble, high molecular weight cytochrome P450 enzymes.

4-Hydroxylation of estradiol by human uterine myometrium and myoma microsomes: implications for the mechanism of uterine tumorigenesis.

Estradiol is converted to catechol estrogens via 2- and 4-hydroxylation by cytochrome P450 enzymes. 4-Hydroxyestradiol elicits biological activities distinct from estradiol, most notably an oxidant stress response induced by free radicals generated by metabolic redox cycling reactions. In this study, we have examined 2- and 4-hydroxylation of estradiol by microsomes of human uterine myometrium and of associated myomata. In all eight cases studied, estradiol 4-hydroxylation by myoma has been substantially elevated relative to surrounding myometrial tissue (minimum, 2-fold; mean, 5-fold). Estradiol 2-hydroxylation in myomata occurs at much lower rates than 4-hydroxylation (ratio of 4-hydroxyestradiol/2-hydroxyestradiol, 7.9 +/- 1.4) and does not significantly differ from rates in surrounding myometrial tissue. Rates of myometrial 2-hydroxylation of estradiol were also not significantly different from values in patients without myomata. We have used various inhibitors to establish that 4-hydroxylation is catalyzed by a completely different cytochrome P450 than 2-hydroxylation. In myoma, alpha-naphthoflavone and a set of ethynyl polycyclic hydrocarbon inhibitors (5 microM) each inhibited 4-hydroxylation more efficiently (up to 90%) than 2-hydroxylation (up to 40%), indicating > 10-fold differences in Ki (<0.5 microM vs. > 5 microM). These activities were clearly distinguished from the selective 2-hydroxylation of estradiol in placenta by aromatase reported previously (low Km, inhibition by Fadrozole hydrochloride or ICI D1033). 4-Hydroxylation was also selectively inhibited relative to 2-hydroxylation by antibodies raised against cytochrome P450 IB1 (rat) (53 vs. 17%). These data indicate that specific 4-hydroxylation of estradiol in human uterine tissues is catalyzed by a form(s) of cytochrome P450 related to P450 IB1, which contribute(s) little to 2-hydroxylation. This enzyme(s) is therefore a marker for uterine myomata and may play a role in the etiology of the tumor.

Molecular Characterization of an Arabidopsis Gene Encoding Hydroperoxide Lyase, a Cytochrome P-450 That Is Wound Inducible1

Hydroperoxide lyase (HPL) cleaves lipid hydroperoxides to produce volatile flavor molecules and also potential signal molecules. We have characterized a gene from Arabidopsis that is homologous to a recently cloned HPL from green pepper (Capsicum annuum). The deduced protein sequence indicates that this gene encodes a cytochrome P-450 with a structure similar to that of allene oxide synthase. The gene was cloned into an expression vector and expressed in Escherichia coli to demonstrate HPL activity. Significant HPL activity was evident when 13S-hydroperoxy-9(Z),11(E),15(Z)-octadecatrienoic acid was used as the substrate, whereas activity with 13S-hydroperoxy-9(Z),11(E)-octadecadienoic acid was approximately 10-fold lower. Analysis of headspace volatiles by gas chromatography-mass spectrometry, after addition of the substrate to E. coli extracts expressing the protein, confirmed enzyme-activity data, since cis-3-hexenal was produced by the enzymatic activity of the encoded protein, whereas hexanal production was limited. Molecular characterization of this gene indicates that it is expressed at high levels in floral tissue and is wound inducible but, unlike allene oxide synthase, it is not induced by treatment with methyl jasmonate.

Translation of cytochrome f is autoregulated through the 5′ untranslated region of petA mRNA in Chlamydomonas chloroplasts

A process that we refer to as control by epistasy of synthesis (CES process) occurs during chloroplast protein biogenesis in Chlamydomonas reinhardtii: the synthesis of some chloroplast-encoded subunits, the CES subunits, is strongly attenuated when some other subunits from the same complex, the dominant subunits, are missing. Herein we investigate the molecular basis of the CES process for the biogenesis of the cytochrome b6f complex and show that negative autoregulation of cytochrome f translation occurs in the absence of other complex subunits. This autoregulation is mediated by an interaction, either direct or indirect, between the 5′ untranslated region of petA mRNA, which encodes cytochrome f, and the C-terminal domain of the unassembled protein. This model for the regulation of cytochrome f translation explains both the decreased rate of cytochrome f synthesis in vivo in the absence of its assembly partners and its increase in synthesis when significant accumulation of the C-terminal domain of the protein is prevented. When expressed from a chimeric mRNA containing the atpA 5′ untranslated region, cytochrome f no longer showed an assembly-dependent regulation of translation. Conversely, the level of antibiotic resistance conferred by a chimeric petA-aadA-rbcL gene was shown to depend on the state of assembly of cytochrome b6f complexes and on the accumulation of the C-terminal domain of cytochrome f. We discuss the possible ubiquity of the CES process in organellar protein biogenesis.

A gene cluster for macrolide antibiotic biosynthesis in Streptomyces venezuelae: Architecture of metabolic diversity

In a survey of microbial systems capable of generating unusual metabolite structural variability, Streptomyces venezuelae ATCC 15439 is notable in its ability to produce two distinct groups of macrolide antibiotics. Methymycin and neomethymycin are derived from the 12-membered ring macrolactone 10-deoxymethynolide, whereas narbomycin and pikromycin are derived from the 14-membered ring macrolactone, narbonolide. This report describes the cloning and characterization of the biosynthetic gene cluster for these antibiotics. Central to the cluster is a polyketide synthase locus (pikA) that encodes a six-module system comprised of four multifunctional proteins, in addition to a type II thioesterase (TEII). Immediately downstream is a set of genes for desosamine biosynthesis (des) and macrolide ring hydroxylation. The study suggests that Pik TEII plays a role in forming a metabolic branch through which polyketides of different chain length are generated, and the glycosyl transferase (encoded by desVII) has the ability to catalyze glycosylation of both the 12- and 14-membered ring macrolactones. Moreover, the pikC-encoded P450 hydroxylase provides yet another layer of structural variability by introducing regiochemical diversity into the macrolide ring systems. The data support the notion that the architecture of the pik gene cluster as well as the unusual substrate specificity of particular enzymes contributes to its ability to generate four macrolide antibiotics.

Replacement of the proximal heme thiolate ligand in chloroperoxidase with a histidine residue

Chloroperoxidase is a versatile heme enzyme which can cross over the catalytic boundaries of other oxidative hemoproteins and perform multiple functions. Chloroperoxidase, in addition to catalyzing classical peroxidative reactions, also acts as a P450 cytochrome and a potent catalase. The multiple functions of chloroperoxidase must be derived from its unique active site structure. Chloroperoxidase possesses a proximal cysteine thiolate heme iron ligand analogous to the P450 cytochromes; however, unlike the P450 enzymes, chloroperoxidase possesses a very polar environment distal to its heme prosthetic group and contains a glutamic acid residue in close proximity to the heme iron. The presence of a thiolate ligand in chloroperoxidase has long been thought to play an essential role in its chlorination and epoxidation activities; however, the research reported in this paper proves that hypothesis to be invalid. To explore the role of Cys-29, the amino acid residue supplying the thiolate ligand in chloroperoxidase, Cys-29 has been replaced with a histidine residue. Mutant clones of the chloroperoxidase genome have been expressed in a Caldariomyces fumago expression system by using gene replacement rather than gene insertion technology. C. fumago produces wild-type chloroperoxidase, thus requiring gene replacement of the wild type by the mutant gene. To the best of our knowledge, this is the first time that gene replacement has been reported for this type of fungus. The recombinant histidine mutants retain most of their chlorination, peroxidation, epoxidation, and catalase activities. These results downplay the importance of a thiolate ligand in chloroperoxidase and suggest that the distal environment of the heme active site plays the major role in maintaining the diverse activities of this enzyme.

Three-dimensional structure of NADPH–cytochrome P450 reductase: Prototype for FMN- and FAD-containing enzymes

Microsomal NADPH–cytochrome P450 reductase (CPR) is one of only two mammalian enzymes known to contain both FAD and FMN, the other being nitric-oxide synthase. CPR is a membrane-bound protein and catalyzes electron transfer from NADPH to all known microsomal cytochromes P450. The structure of rat liver CPR, expressed in Escherichia coli and solubilized by limited trypsinolysis, has been determined by x-ray crystallography at 2.6 Å resolution. The molecule is composed of four structural domains: (from the N- to C- termini) the FMN-binding domain, the connecting domain, and the FAD- and NADPH-binding domains. The FMN-binding domain is similar to the structure of flavodoxin, whereas the two C-terminal dinucleotide-binding domains are similar to those of ferredoxin–NADP+ reductase (FNR). The connecting domain, situated between the FMN-binding and FNR-like domains, is responsible for the relative orientation of the other domains, ensuring the proper alignment of the two flavins necessary for efficient electron transfer. The two flavin isoalloxazine rings are juxtaposed, with the closest distance between them being about 4 Å. The bowl-shaped surface near the FMN-binding site is likely the docking site of cytochrome c and the physiological redox partners, including cytochromes P450 and b5 and heme oxygenase.

Mitochondrial DNA rearrangements of Podospora anserina are under the control of the nuclear gene grisea

Podospora anserina is a filamentous fungus with a limited life span. Life span is controlled by nuclear and extranuclear genetic traits. Herein we report the nature of four alterations in the nuclear gene grisea that lead to an altered morphology, a defect in the formation of female gametangia, and an increased life span. Three sequence changes are located in the 5′ upstream region of the grisea ORF. One mutation is a G → A transition at the 5′ splice site of the single intron of the gene, leading to a RNA splicing defect. This loss-of-function affects the amplification of the first intron of the mitochondrial cytochrome c oxidase subunit I gene (COI) and the specific mitochondrial DNA rearrangements that occur during senescence of wild-type strains. Our results indicate that the nuclear gene grisea is part of a molecular machinery involved in the control of mitochondrial DNA reorganizations. These DNA instabilities accelerate but are not a prerequisite for the aging of P. anserina cultures.

Tracking molecular evolution of photosynthesis by characterization of a major photosynthesis gene cluster from Heliobacillus mobilis

A DNA sequence has been obtained for a 35.6-kb genomic segment from Heliobacillus mobilis that contains a major cluster of photosynthesis genes. A total of 30 ORFs were identified, 20 of which encode enzymes for bacteriochlorophyll and carotenoid biosynthesis, reaction-center (RC) apoprotein, and cytochromes for cyclic electron transport. Donor side electron-transfer components to the RC include a putative RC-associated cytochrome c553 and a unique four-large-subunit cytochrome bc complex consisting of Rieske Fe-S protein (encoded by petC), cytochrome b6 (petB), subunit IV (petD), and a diheme cytochrome c (petX). Phylogenetic analysis of various photosynthesis gene products indicates a consistent grouping of oxygenic lineages that are distinct and descendent from anoxygenic lineages. In addition, H. mobilis was placed as the closest relative to cyanobacteria, which form a monophyletic origin to chloroplast-based photosynthetic lineages. The consensus of the photosynthesis gene trees also indicates that purple bacteria are the earliest emerging photosynthetic lineage. Our analysis also indicates that an ancient gene-duplication event giving rise to the paralogous bchI and bchD genes predates the divergence of all photosynthetic groups. In addition, our analysis of gene duplication of the photosystem I and photosystem II core polypeptides supports a “heterologous fusion model” for the origin and evolution of oxygenic photosynthesis.

Structure of a cytochrome P450–redox partner electron-transfer complex

The crystal structure of the complex between the heme- and FMN-binding domains of bacterial cytochrome P450BM-3, a prototype for the complex between eukaryotic microsomal P450s and P450 reductase, has been determined at 2.03 Å resolution. The flavodoxin-like flavin domain is positioned at the proximal face of the heme domain with the FMN 4.0 and 18.4 Å from the peptide that precedes the heme-binding loop and the heme iron, respectively. The heme-binding peptide represents the most efficient and coupled through-bond electron pathway to the heme iron. Substantial differences between the FMN-binding domains of P450BM-3 and microsomal P450 reductase, observed around the flavin-binding sites, are responsible for different redox properties of the FMN, which, in turn, control electron flow to the P450.

Heme transfer to the heme chaperone CcmE during cytochrome c maturation requires the CcmC protein, which may function independently of the ABC-transporter CcmAB

Cytochrome c maturation in Escherichia coli requires the ccm operon, which encodes eight membrane proteins (CcmABCDEFGH). CcmE is a periplasmic heme chaperone that binds heme covalently and transfers it onto apocytochrome c in the presence of CcmF, CcmG, and CcmH. In this work we addressed the functions of the ccmABCD gene products with respect to holo-CcmE formation and the subsequent ligation of heme to apocytochrome c. In the absence of the ccmABCD genes, heme is not bound to CcmE. We report that CcmC is functionally uncoupled from the ABC transporter subunits CcmA and CcmB, because it is the only Ccm protein that is strictly required for heme transfer and attachment to CcmE. Site-directed mutagenesis of conserved histidines inactivates the CcmC protein, which is in agreement with the hypothesis that this protein interacts directly with heme. We also present evidence that questions the role of CcmAB as a heme exporter; yet, the transported substrate remains unknown. CcmD was found to be involved in stabilizing the heme chaperone CcmE in the membrane. We propose a heme-trafficking pathway as part of a substantially revised model for cytochrome c maturation in E. coli.

Farnesyltransferase inhibitors induce cytochrome c release and caspase 3 activation preferentially in transformed cells

Farnesyltransferase inhibitors (FTIs) represent a new class of anticancer drugs that show promise in blocking the growth of tumors. Here, we report that FTIs are capable of inducing apoptosis of transformed but not untransformed cells. Treatment of v-K-ras-transformed normal rat kidney (KNRK) cells with FTIs leads to the induction of apoptotic cell morphology, chromatin condensation and DNA fragmentation. In addition, fluorescence-activated cell sorter analysis of FTI-treated KNRK cells shows a sub-G1 apoptotic peak (chromosome content of <2 N). This FTI-induced apoptosis is evident only when the cells are grown in low serum conditions (0.1% fetal calf serum) and is observed selectively with transformed KNRK cells and not with untransformed NRK cells. Further analysis of the mechanism underlying this apoptosis has shown that FTI treatment of KNRK cells results in the activation of caspase 3 but not caspase 1. Moreover, the addition of Z-DEVD-fmk, an agent that interferes with caspase 3 activity, can inhibit FTI-induced apoptosis in a dose-dependent manner. Introduction of the CASP-3 gene into MCF7 cells, which lack caspase 3 activity, results in a significant increase of FTI-induced apoptosis. Furthermore, FTI induces the release of cytochrome c into the cytosol. This release is an important feature of caspase 3-mediated apoptosis. These results suggest that FTIs induce apoptosis through the release of cytochrome c from the mitochondria resulting in caspase 3 activation.