The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Silicatein α: Cathepsin L-like protein in sponge biosilica

Earth’s biota produces vast quantities of polymerized silica at ambient temperatures and pressures by mechanisms that are not understood. Silica spicules constitute 75% of the dry weight of the sponge Tethya aurantia, making this organism uniquely tractable for analyses of the proteins intimately associated with the biosilica. Each spicule contains a central protein filament, shown by x-ray diffraction to exhibit a highly regular, repeating structure. The protein filaments can be dissociated to yield three similar subunits, named silicatein α, β, and γ. The molecular weights and amino acid compositions of the three silicateins are similar, suggesting that they are members of a single protein family. The cDNA sequence of silicatein α, the most abundant of these subunits, reveals that this protein is highly similar to members of the cathepsin L and papain family of proteases. The cysteine at the active site in the proteases is replaced by serine in silicatein α, although the six cysteines that form disulfide bridges in the proteases are conserved. Silicatein α also contains unique tandem arrays of multiple hydroxyls. These structural features may help explain the mechanism of biosilicification and the recently discovered activity of the silicateins in promoting the condensation of silica and organically modified siloxane polymers (silicones) from the corresponding silicon alkoxides. They suggest the possibility of a dynamic role of the silicateins in silicification of the sponge spicule and offer the prospect of a new synthetic route to silica and siloxane polymers at low temperature and pressure and neutral pH.

The domain structure of protoporphyrinogen oxidase, the molecular target of diphenyl ether-type herbicides

Protoporphyrinogen oxidase (EC 1–3-3–4), the 60-kDa membrane-bound flavoenzyme that catalyzes the final reaction of the common branch of the heme and chlorophyll biosynthesis pathways in plants, is the molecular target of diphenyl ether-type herbicides. It is highly resistant to proteases (trypsin, endoproteinase Glu-C, or carboxypeptidases A, B, and Y), because the protein is folded into an extremely compact form. Trypsin maps of the native purified and membrane-bound yeast protoporphyrinogen oxidase show that this basic enzyme (pI > 8.5) was cleaved at a single site under nondenaturing conditions, generating two peptides with relative molecular masses of 30,000 and 35,000. The endoproteinase Glu-C also cleaved the protein into two peptides with similar masses, and there was no additional cleavage site under mild denaturing conditions. N-terminal peptide sequence analysis of the proteolytic (trypsin and endoproteinase Glu-C) peptides showed that both cleavage sites were located in putative connecting loop between the N-terminal domain (25 kDa) with the βαβ ADP-binding fold and the C-terminal domain (35 kDa), which possibly is involved in the binding of the isoalloxazine moiety of the FAD cofactor. The peptides remained strongly associated and fully active with the Km for protoporphyrinogen and the Ki for various inhibitors, diphenyl-ethers, or diphenyleneiodonium derivatives, identical to those measured for the native enzyme. However, the enzyme activity of the peptides was much more susceptible to thermal denaturation than that of the native protein. Only the C-terminal domain of protoporphyrinogen oxidase was labeled specifically in active site-directed photoaffinity-labeling experiments. Trypsin may have caused intramolecular transfer of the labeled group to reactive components of the N-terminal domain, resulting in nonspecific labeling. We suggest that the active site of protoporphyrinogen oxidase is in the C-terminal domain of the protein, at the interface between the C- and N-terminal domains.

Molecular recognition of angiogenesis inhibitors fumagillin and ovalicin by methionine aminopeptidase 2

Angiogenesis inhibitors are a novel class of promising therapeutic agents for treating cancer and other human diseases. Fumagillin and ovalicin compose a class of structurally related natural products that potently inhibit angiogenesis by blocking endothelial cell proliferation. A synthetic analog of fumagillin, TNP-470, is currently undergoing clinical trials for treatment of a variety of cancers. A common target for fumagillin and ovalicin recently was identified as the type 2 methionine aminopeptidase (MetAP2). These natural products bind MetAP2 covalently, inhibiting its enzymatic activity. The specificity of this binding is underscored by the lack of inhibition of the closely related type 1 enzyme, MetAP1. The molecular basis of the high affinity and specificity of these inhibitors for MetAP2 has remained undiscovered. To determine the structural elements of these inhibitors and MetAP2 that are involved in this interaction, we synthesized fumagillin analogs in which each of the potentially reactive epoxide groups was removed either individually or in combination. We found that the ring epoxide in fumagillin is involved in the covalent modification of MetAP2, whereas the side chain epoxide group is dispensable. By using a fumagillin analog tagged with fluorescein, His-231 in MetAP2 was identified as the residue that is covalently modified by fumagillin. Site-directed mutagenesis of His-231 demonstrated its importance for the catalytic activity of MetAP2 and confirmed that the same residue is covalently modified by fumagillin. These results, in agreement with a recent structural study, suggest that fumagillin and ovalicin inhibit MetAP2 by irreversible blockage of the active site.

Protein component of the ribozyme ribonuclease P alters substrate recognition by directly contacting precursor tRNA

The protein component of ribonuclease P (RNase P) binds to the RNA subunit, forming a functional ribonucleoprotein complex in vivo and enhancing the affinity of the precursor tRNA (pre-tRNA) substrate. Photocrosslinking experiments with pre-tRNA bound to RNase P reconstituted with the protein component of Bacillus subtilis ribonuclease P (P protein) site specifically modified with a crosslinking reagent indicate that: (i) the central cleft of P protein directly interacts with the single-stranded 5′ leader sequence of pre-tRNA, and (ii) the orientation and register of the pre-tRNA leader sequence in the central cleft places the protein component in close proximity to the active site. This unique mode of interaction suggests that the catalytic active site in RNase P occurs near the interface of RNA and protein. In contrast to other ribonucleoprotein complexes where the protein mainly stabilizes the active tertiary fold of the RNA, a critical function of the protein component of RNase P is to alter substrate specificity and enhance catalytic efficiency.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.

Studies on the enzymatic and transcriptional activity of the dimerization cofactor for hepatocyte nuclear factor 1

The relationship between the enzymatic and the transcriptional activity of the bifunctional protein pterin-4a-carbinolamine dehydratase/dimerization cofactor for hepatocyte nuclear factor 1 (DCoH) has been elucidated by site-directed mutagenesis. DCoH dimers harbor a binding site for hepatocyte nuclear factor 1 (HNF1), two active centers that bind pterins, and a saddle-shaped surface that resembles nucleic acid binding domains. Two domains of the protein have been selectively targeted to determine if a change in one activity affects the other. No strong correlation has been found, supporting the idea that carbinolamine dehydratase activity is not required for HNF1 binding in vitro or transcriptional coactivation in vivo. Double mutations in the active center, however, influence the in vivo transcriptional activity but not HNF1 binding. This finding suggests that some active center residues also are used during transcription, possibly for binding of another (macro)molecule. Several mutations in the saddle led to a surprising increase in transcription, therefore linking this domain to transcriptional regulation as well. The transcriptional function of DCoH therefore is composed of two parts, HNF1 binding and another contributing effect that involves the active site and, indirectly, the saddle.

B12-dependent ribonucleotide reductases from deeply rooted eubacteria are structurally related to the aerobic enzyme from Escherichia coli

The ribonucleotide reductases from three ancient eubacteria, the hyperthermophilic Thermotoga maritima (TM), the radioresistant Deinococcus radiodurans (DR), and the thermophilic photosynthetic Chloroflexus aurantiacus, were found to be coenzyme-B12 (class II) enzymes, similar to the earlier described reductases from the archaebacteria Thermoplasma acidophila and Pyrococcus furiosus. Reduction of CDP by the purified TM and DR enzymes requires adenosylcobalamin and DTT. dATP is a positive allosteric effector, but stimulation of the TM enzyme only occurs close to the temperature optimum of 80–90°C. The TM and DR genes were cloned by PCR from peptide sequence information. The TM gene was sequenced completely and expressed in Escherichia coli. The deduced amino acid sequences of the two eubacterial enzymes are homologous to those of the archaebacteria. They can also be aligned to the sequence of the large protein of the aerobic E. coli ribonucleotide reductase that belongs to a different class (class I), which is not dependent on B12. Structure determinations of the E. coli reductase complexed with substrate and allosteric effectors earlier demonstrated a 10-stranded β/α-barrel in the active site. From the conservation of substrate- and effector-binding residues we propose that the B12-dependent class II enzymes contain a similar barrel.

Interaction of the herbicide glyphosate with its target enzyme 5-enolpyruvylshikimate 3-phosphate synthase in atomic detail

Biosynthesis of aromatic amino acids in plants, many bacteria, and microbes relies on the enzyme 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, a prime target for drugs and herbicides. We have identified the interaction of EPSP synthase with one of its two substrates (shikimate 3-phosphate) and with the widely used herbicide glyphosate by x-ray crystallography. The two-domain enzyme closes on ligand binding, thereby forming the active site in the interdomain cleft. Glyphosate appears to occupy the binding site of the second substrate of EPSP synthase (phosphoenol pyruvate), mimicking an intermediate state of the ternary enzyme⋅substrates complex. The elucidation of the active site of EPSP synthase and especially of the binding pattern of glyphosate provides a valuable roadmap for engineering new herbicides and herbicide-resistant crops, as well as new antibiotic and antiparasitic drugs.

A 1.2-Å snapshot of the final step of bacterial cell wall biosynthesis

The cell wall imparts structural strength and shape to bacteria. It is made up of polymeric glycan chains with peptide branches that are cross-linked to form the cell wall. The cross-linking reaction, catalyzed by transpeptidases, is the last step in cell wall biosynthesis. These enzymes are members of the family of penicillin-binding proteins, the targets of β-lactam antibiotics. We report herein the structure of a penicillin-binding protein complexed with a cephalosporin designed to probe the mechanism of the cross-linking reaction catalyzed by transpeptidases. The 1.2-Å resolution x-ray structure of this cephalosporin bound to the active site of the bifunctional serine type d-alanyl-d-alanine carboxypeptidase/transpeptidase (EC 3.4.16.4) from Streptomyces sp. strain R61 reveals how the two peptide strands from the polymeric substrates are sequestered in the active site of a transpeptidase. The structure of this complex provides a snapshot of the enzyme and the bound cell wall components poised for the final and critical cross-linking step of cell wall biosynthesis.

Tight Asp-85–Thr-89 association during the pump switch of bacteriorhodopsin

Unidirectional proton transport in bacteriorhodopsin is enforced by the switching machinery of the active site. Threonine 89 is located in this region, with its O—H group forming a hydrogen bond with Asp-85, the acceptor for proton transfer from the Schiff base of the retinal chromophore. Previous IR spectroscopy of [3-18O]threonine-labeled bacteriorhodopsin showed that the hydrogen bond of the O—D group of Thr-89 in D2O is strengthened in the K photocycle intermediate. Here, we show that the strength and orientation of this hydrogen bond remains unchanged in the L intermediate and through the M intermediate. Furthermore, a strong interaction between Asp-85 and the O—H (O—D) group of Thr-89 in M is indicated by a shift in the C⩵O stretching vibration of the former because of 18O substitution in the latter. Thus, the strong hydrogen bond between Asp-85 and Thr-89 in K persists through M, contrary to structural models based on x-ray crystallography of the photocycle intermediates. We propose that, upon photoisomerization of the chromophore, Thr-89 forms a tight, persistent complex with one of the side-chain oxygens of Asp-85 and is thereby precluded from participating in the switching process. On the other hand, the loss of hydrogen bonding at the other oxygen of Asp-85 in M may be related to the switching event.

Phylogenetic relationships among group II intron ORFs

Group II introns are widely believed to have been ancestors of spliceosomal introns, yet little is known about their own evolutionary history. In order to address the evolution of mobile group II introns, we have compiled 71 open reading frames (ORFs) related to group II intron reverse transcriptases and subjected their derived amino acid sequences to phylogenetic analysis. The phylogenetic tree was rooted with reverse transcriptases (RTs) of non-long terminal repeat retroelements, and the inferred phylogeny reveals two major clusters which we term the mitochondrial and chloroplast-like lineages. Bacterial ORFs are mainly positioned at the bases of the two lineages but with weak bootstrap support. The data give an overview of an apparently high degree of horizontal transfer of group II intron ORFs, mostly among related organisms but also between organelles and bacteria. The Zn domain (nuclease) and YADD motif (RT active site) were lost multiple times during evolution. Differences in domain structures suggest that the oldest ORFs were concise, while the ORF in the mitochondrial lineage subsequently expanded in three locations. The data are consistent with a bacterial origin for mobile group II introns.

Global analysis of proteasomal substrate specificity using positional-scanning libraries of covalent inhibitors

The proteasome is a large protease complex consisting of multiple catalytic subunits that function simultaneously to digest protein substrates. This complexity has made deciphering the role each subunit plays in the generation of specific protein fragments difficult. Positional scanning libraries of peptide vinyl sulfones were generated in which the amino acid located directly at the site of hydrolysis (P1 residue) was held constant and sequences distal to that residue (P2, P3, and P4 positions) were varied across all natural amino acids (except cysteine and methionine). Binding information for each of the individual catalytic subunits was obtained for each library under a variety of different conditions. The resulting specificity profiles indicated that substrate positions distal to P1 are critical for directing substrates to active subunits in the complex. Furthermore, specificity profiles of IFN-γ-regulated subunits closely matched those of their noninducible counterparts, suggesting that subunit swapping may modulate substrate processing by a mechanism that does require a change in the primary sequence specificity of individual catalytic subunits in the complex. Finally, specificity profiles were used to design specific inhibitors of a single active site in the complex. These reagents can be used to further establish the role of each subunit in substrate processing by the proteasome.

Crystal structure of cytochrome P450 14α-sterol demethylase (CYP51) from Mycobacterium tuberculosis in complex with azole inhibitors

Cytochrome P450 14α-sterol demethylases (CYP51) are essential enzymes in sterol biosynthesis in eukaryotes. CYP51 removes the 14α-methyl group from sterol precursors such as lanosterol, obtusifoliol, dihydrolanosterol, and 24(28)-methylene-24,25-dihydrolanosterol. Inhibitors of CYP51 include triazole antifungal agents fluconazole and itraconazole, drugs used in treatment of topical and systemic mycoses. The 2.1- and 2.2-Å crystal structures reported here for 4-phenylimidazole- and fluconazole-bound CYP51 from Mycobacterium tuberculosis (MTCYP51) are the first structures of an authentic P450 drug target. MTCYP51 exhibits the P450 fold with the exception of two striking differences—a bent I helix and an open conformation of BC loop—that define an active site-access channel running along the heme plane perpendicular to the direction observed for the substrate entry in P450BM3. Although a channel analogous to that in P450BM3 is evident also in MTCYP51, it is not open at the surface. The presence of two different channels, with one being open to the surface, suggests the possibility of conformationally regulated substrate-in/product-out openings in CYP51. Mapping mutations identified in Candida albicans azole-resistant isolates indicates that azole resistance in fungi develops in protein regions involved in orchestrating passage of CYP51 through different conformational stages along the catalytic cycle rather than in residues directly contacting fluconazole. These new structures provide a basis for rational design of new, more efficacious antifungal agents as well as insight into the molecular mechanism of P450 catalysis.

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 1010-fold: but how? Our structures of the yeast PPase product complex at 1.15 Å and fluoride-inhibited complex at 1.9 Å visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (Onu) directly attacks PPi with a pKa vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and Onu, in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.