Toward a generalized computational workflow for exploiting transient pockets as new targets for small molecule stabilizers: Application to the homogentisate 1,2-dioxygenase mutants at the base of rare disease Alkaptonuria

Alkaptonuria (AKU) is an inborn error of metabolism where mutation of homogentisate 1,2-dioxygenase (HGD) gene leads to a deleterious or misfolded product with subsequent loss of enzymatic degradation of homogentisic acid (HGA) whose accumulation in tissues causes ochronosis and degeneration. There is no licensed therapy for AKU. Many missense mutations have been individuated as responsible for quaternary structure disruption of the native hexameric HGD. A new approach to the treatment of AKU is here proposed aiming to totally or partially rescue enzyme activity by targeting of HGD with pharmacological chaperones, i.e. small molecules helping structural stability. Co-factor pockets from oligomeric proteins have already been successfully exploited as targets for such a strategy, but no similar sites are present at HGD surface; hence, transient pockets are here proposed as a target for pharmacological chaperones. Transient pockets are detected along the molecular dynamics trajectory of the protein and filtered down to a set of suitable sites for structural stabilization by mean of biochemical and pharmacological criteria. The result is a computational workflow relevant to other inborn errors of metabolism requiring rescue of oligomeric, misfolded enzymes.     

Structural basis of regiospecificity of a mononuclear iron enzyme in antibiotic fosfomycin biosynthesis

Hydroxypropylphosphonic acid epoxidase (HppE) is an unusual mononuclear iron enzyme that uses dioxygen to catalyze the oxidative epoxidation of (S)-2-hydroxypropylphosphonic acid (S-HPP) in the biosynthesis of the antibiotic fosfomycin. Additionally, the enzyme converts the R-enantiomer of the substrate (R-HPP) to 2-oxo-propylphosphonic acid. To probe the mechanism of HppE regiospecificity, we determined three X-ray structures: R-HPP with inert cobalt-containing enzyme (Co(II)-HppE) at 2.1 ? resolution; R-HPP with active iron-containing enzyme (Fe(II)-HppE) at 3.0 ? resolution; and S-HPP-Fe(II)-HppE in complex with dioxygen mimic NO at 2.9 ? resolution. These structures, along with previously determined structures of S-HPP-HppE, identify the dioxygen binding site on iron and elegantly illustrate how HppE is able to recognize both substrate enantiomers to catalyze two completely distinct reactions.     

Product bound structures of the soluble methane monooxygenase hydroxylase from Methylococcus capsulatus (Bath): protein motion in the alpha-subunit

The soluble methane monooxygenase hydroxylase (MMOH) alpha-subunit contains a series of cavities that delineate the route of substrate entrance to and product egress from the buried carboxylate-bridged diiron center. The presence of discrete cavities is a major structural difference between MMOH, which can hydroxylate methane, and toluene/o-xylene monooxygenase hydroxylase (ToMOH), which cannot. To understand better the functions of the cavities and to investigate how an enzyme designed for methane hydroxylation can also accommodate larger substrates such as octane, methylcubane, and trans-1-methyl-2-phenylcyclopropane, MMOH crystals were soaked with an assortment of different alcohols and their X-ray structures were solved to 1.8-2.4 A resolution. The product analogues localize to cavities 1-3 and delineate a path of product exit and/or substrate entrance from the active site to the surface of the protein. The binding of the alcohols to a position bridging the two iron atoms in cavity 1 extends and validates previous crystallographic, spectroscopic, and computational work indicating this site to be where substrates are hydroxylated and products form. The presence of these alcohols induces perturbations in the amino acid side-chain gates linking pairs of cavities, allowing for the formation of a channel similar to one observed in ToMOH. Upon binding of 6-bromohexan-1-ol, the pi helix formed by residues 202-211 in helix E of the alpha-subunit is extended through residue 216, changing the orientations of several amino acid residues in the active site cavity. This remarkable secondary structure rearrangement in the four-helix bundle has several mechanistic implications for substrate accommodation and the function of the effector protein, MMOB.     

Synthesis and characterisation of a ligand that forms a stable tetrahedral intermediate in the active site of the Aureobacterium species (-) gamma-lactamase

The crystal structure of a (-) gamma-lactamase from an Aureobacterium species showed a molecule bound covalently to the active site serine residue. This enzyme complex represented the first structure of a stably bound tetrahedral intermediate for an alpha/beta hydrolase fold enzyme. The structural elucidation of tetrahedral intermediates is important for the understanding of enzymatic mechanism, substrate recognition and enzyme inhibition. In this paper, we report the synthesis and subsequent characterisation of (3aR,7aS)-3a,4,7,7a-tetrahydrobenzo-[1,3]-dioxol-2-one (BD1), the molecule modelled into the Aureobacterium (-) gamma-lactamase active site. This molecule has been confirmed to be an inhibitor and to be displaced from the enzyme by the racemic gamma-lactam substrate.     

Hydroxylase activity of Met471Cys tyramine beta-monooxygenase

A series of mutations was targeted at the methionine residue, Met471, coordinating the Cu(M) site of tyramine beta-monooxygenase (TbetaM). The methionine ligand at Cu(M) is believed to be key to dioxygen activation and the hydroxylation chemistry of the copper monooxygenases. The reactivity and copper binding properties of three TbetaM mutants, Met471Asp, Met471Cys, and Met471His, were examined. All three mutants show similar metal binding affinities to wild type TbetaM in the oxidized enzyme forms. EPR spectroscopy suggests that the Cu(II) coordination geometry is identical to that of the WT enzyme. However, substrate hydroxylation was observed for the reaction of tyramine solely with Met471Cys TbetaM. Met471Cys TbetaM provides the first example of an active mutant directed at the Cu(M) site of this class of hydroxylases. The reactivity and altered kinetics of the Met471Cys mutant further highlight the central role of the methionine residue in the enzyme mechanism. The sole ability of the cysteine residue to support activity among the series of alternate amino acids investigated is relevant to theoretical and biomimetic investigations of dioxygen activation at mononuclear copper centers.     

Hydrophobic core flexibility modulates enzyme activity in HIV-1 protease

Human immunodeficiency virus Type-1 (HIV-1) protease is crucial for viral maturation and infectivity. Studies of protease dynamics suggest that the rearrangement of the hydrophobic core is essential for enzyme activity. Many mutations in the hydrophobic core are also associated with drug resistance and may modulate the core flexibility. To test the role of flexibility in protease activity, pairs of cysteines were introduced at the interfaces of flexible regions remote from the active site. Disulfide bond formation was confirmed by crystal structures and by alkylation of free cysteines and mass spectrometry. Oxidized and reduced crystal structures of these variants show the overall structure of the protease is retained. However, cross-linking the cysteines led to drastic loss in enzyme activity, which was regained upon reducing the disulfide cross-links. Molecular dynamics simulations showed that altered dynamics propagated throughout the enzyme from the engineered disulfide. Thus, altered flexibility within the hydrophobic core can modulate HIV-1 protease activity, supporting the hypothesis that drug resistant mutations distal from the active site can alter the balance between substrate turnover and inhibitor binding by modulating enzyme activity.     

Immobilization of the enzyme beta-lactamase on biotin-derivatized poly(L-lysine)-g-poly(ethylene glycol)-coated sensor chips: a study on oriented attachment and surface activity by enzyme kinetics and in situ optical sensing

Understanding the conformation, orientation, and specific activity of proteins bound to surfaces is crucial for the development and optimization of highly specific and sensitive biosensors. In this study, the very efficient enzyme beta-lactamase is used as a model protein. The wild-type form was genetically engineered by site-directed mutagenesis to introduce single cysteine residues on the surface of the enzyme. The cysteine thiol group is subsequently biotinylated with a dithiothreitol (DTT)-cleavable biotinylation reagent. beta-Lactamase is then immobilized site-specifically via the biotin group on neutral avidin-covered surfaces with the aim to control the orientation of the enzyme molecule at the surface and study its effect on enzymatic activity using Nitrocefin as the substrate. The DTT-cleavable spacer allows the release of the specifically bound enzyme from the surface. Immobilization of the enzyme is performed on a monolayer of the polycationic, biotinylated polymer PLL-g-PEG/PEG-biotin assembled on niobium oxide (Nb2O5) surfaces via neutral avidin as the docking site. Two different assembly protocols, the sequential adsorption of avidin and biotinylated beta-lactamase and the immobilization of preformed complexes of beta-lactamase and avidin, are compared in terms of immobilization efficiency. In situ optical waveguide lightmode spectroscopy and colorimetric analysis of enzymatic activity were used to distinguish between specific and unspecific enzyme adsorption, to sense quantitatively the amount of immobilized enzyme, and to determine Michaelis-Menten kinetics. All tested enzyme variants turned out to be active upon immobilization at the polymeric surface. However, the efficiency of immobilized enzymes relative to the soluble enzymes was reduced about sevenfold, mainly because of impaired substrate (Nitrocefin) diffusion or restricted accessibility of the active site. No significant effect of different enzyme orientations could be detected, probably because the enzymes were attached to the surface through long, flexible PEG chain linkers.     

Aggregation behavior of giant amphiphiles prepared by cofactor reconstitution

We report on biohybrid surfactants, termed "giant amphiphiles", in which a protein or an enzyme acts as the polar head group and a synthetic polymer as the apolar tail. It is demonstrated that the modification of horseradish peroxidase (HRP) and myoglobin (Mb) with an apolar polymer chain through the cofactor reconstitution method yields giant amphiphiles that form spherical aggregates (vesicles) in aqueous solution. Both HRP and Mb retain their original functionality when modified with a single polystyrene chain, but reconstitution has an effect on their activities. In the case of HRP the enzymatic activity decreases and for Mb the stability of the dioxygen myoglobin (oxy-Mb) complex is reduced, which is probably the result of a disturbed binding of the heme in the apo-protein or a reduced access of the substrate to the active site of the enzyme or protein.     

Structure and Mechanism Leading to Formation of the Cysteine Sulfinate Product Complex of a Biomimetic Cysteine Dioxygenase Model

Cysteine dioxygenase is a unique nonheme iron enzyme that is involved in the metabolism of cysteine in the body. It contains an iron active site with an unusual 3‐His ligation to the protein, which contrasts with the structural features of common nonheme iron dioxygenases. Recently, some of us reported a truly biomimetic model for this enzyme, namely a trispyrazolylborato iron(II) cysteinato complex, which not only has a structure very similar to the enzyme–substrate complex but also represents a functional model: Treatment of the model with dioxygen leads to cysteine dioxygenation, as shown by isolating the cysteine part of the product in the course of the work‐up. However, little is known on the conversion mechanism and, so far, not even the structure of the actual product complex had been characterised, which is also unknown in case of the enzyme. In a multidisciplinary approach including density functional theory calculations and X‐ray absorption spectroscopy, we have now determined the structure of the actual sulfinato complex for the first time. The Cys‐SO₂^? functional group was found to be bound in an η²‐O,O‐coordination mode, which, based on the excellent resemblance between model and enzyme, also provides the first support for a corresponding binding mode within the enzymatic product complex. Indeed, this is again confirmed by theory, which had predicted a η²‐O,O‐binding mode for synthetic as well as the natural enzyme.     

The role of copper in particulate methane monooxygenase from Methylosinus trichosporium OB3b

The effect of metal ions on particulate methane monooxygenase (pMMO) was studied. The pMMO activity in the membranes was partially inhibited by ethylenediaminetetraacetic acid (EDTA), but remained more than 70% of the as-isolated membranes. The activity of the EDTA-treated membranes was strongly influenced by the addition of metal ions. Among the metal ions, copper ion stimulated the activity, indicating that copper was needed for the activity. When duroquinol and dioxygen were introduced to the EDTA-treated membranes, the electron spin resonance signal of copper did not change, suggesting that the copper cluster did not play as an electron transport and may have another function, such as active site of pMMO or regulator of the activity. On the other hand, the iron signal (g=5.98) decreased by the addition of duroquinol, dioxygen and acetylene, showing an iron atom is contained in the active site of pMMO.     

Enzymes immobilized on montmorillonite: comparison of performance in batch and packed-bed reactors

Summary The enzymes a-amylase, invertase and glucoamylase were immobilized on acid activated montmorillonite using two techniques, viz. adsorption and covalent binding, and their activities were tested in a batch and packed-bed reactor and were compared. The packed-bed reactor showed an improved performance for all immobilized enzymes, which was attributed to lowering of diffusional restrictions to mass transfer. Lower activity in case of batch reactor for immobilized invertase was due to a combined effect of loss of native conformation of enzyme on account of immobilization and mass transfer resistances due to improper diffusion of substrate to the active site of enzyme. For immobilized glucoamylase, the packed-bed reactor demonstrated exceptionally high activity that was very close to the free enzyme. Covalently bound glucoamylase showed higher activity than the free enzyme.     

Oxygen economy of cytochrome P450: what is the origin of the mixed functionality as a dehydrogenase-oxidase enzyme compared with its normal function?

The economy of dioxygen consumption by enzymes constitutes a fundamental problem in enzymatic chemistry (ref 1). Sometimes, the enzyme converts ALL the oxygen into water, without affecting the organic substrate, thereby acting as an "oxidase" (ref 1). Other times, the enzyme converts all the oxygen into water and causes desaturation in the substrate, thus exhibiting a mixed function as both "oxidase" and "dehydrogenase" (refs 2-5). The present paper describes density functional calculations demonstrating that the oxidase-dehydrogenase mixed activity occurs from the cationic intermediate species and requires electro-steric inhibition of the rebound process. Furthermore, the calculations reveal that the carbocation is formally nascent from an excited state of the active species of the enzyme (2Cpd I), in which the Fe=O moiety is singlet coupled as in the 1Deltag state of dioxygen! Thus, our results resolve an important mechanism and reveal the factors that underlie its observability.     

CD and MCD studies of the non-heme ferrous active site in (4-hydroxyphenyl)pyruvate dioxygenase: correlation between oxygen activation in the extradiol and alpha-KG-dependent dioxygenases

(4-Hydroxyphenyl)pyruvate dioxygenase (HPPD) is an unusual alpha-keto acid-dependent non-heme iron dioxygenase as it incorporates both atoms of dioxygen into a single substrate, paralleling the extradiol dioxygenases. CD/MCD studies of the catalytically active ferrous site and its interaction with substrate reveal a geometic and electronic structure and mechanistic approach to oxygen activation which bridges those of the alpha-KG-dependent and the extradiol dioxygenases.     

Spectroscopic studies of 1-aminocyclopropane-1-carboxylic acid oxidase: molecular mechanism and CO(2) activation in the biosynthesis of ethylene

1-Aminocyclopropane-1-carboxylic acid (ACC) oxidase (ACCO) catalyzes the last step in the biosynthesis of the gaseous plant hormone ethylene, which is involved in development, including germination, fruit ripening, and senescence. ACCO is a mononuclear non-heme ferrous enzyme that couples the oxidation of the cosubstrate ascorbate to the oxidation of substrate ACC by dioxygen. In addition to substrate and cosubstrate, ACCO requires the activator CO(2) for continuous turnover. NIR circular dichroism and magnetic circular dichroism spectroscopies have been used to probe the geometric and electronic structure of the ferrous active site in ACCO to obtain molecular-level insight into its catalytic mechanism. Resting ACCO/Fe(II) is coordinatively saturated (six-coordinate). In the presence of CO(2), one ferrous ligand is displaced to yield a five-coordinate site only when both the substrate ACC and cosubstrate ascorbate are bound to the enzyme. The open coordination position allows rapid O(2) activation for the oxidation of both substrates. In the absence of CO(2), ACC binding alone converts the site to five-coordinate, which would react with O(2) in the absence of ascorbate and quickly deactivate the enzyme. These studies show that ACCO employs a general strategy similar to other non-heme iron enzymes in terms of opening iron coordination sites at the appropriate time in the reaction cycle and define the role of CO(2) as stabilizing the six-coordinate ACCO/Fe(II)/ACC complex, thus preventing the uncoupled reaction that inactivates the enzyme.     

Theoretical study of the mechanisms of substrate recognition by catalase

A variety of theoretical methods including classical molecular interaction potentials, classical molecular dynamics, and activated molecular dynamics have been used to analyze the substrate recognition mechanisms of peroxisomal catalase from Saccharomyces cerevisiae. Special attention is paid to the existence of channels connecting the heme group with the exterior of the protein. On the basis of these calculations a rationale is given for the unique catalytic properties of this enzyme, as well as for the change in enzyme efficiency related to key mutations. According to our calculations the water is expected to be a competitive inhibitor of the enzyme, blocking the access of hydrogen peroxide to the active site. The main channel is the preferred route for substrate access to the enzyme and shows a cooperative binding to hydrogen peroxide. However, the overall affinity of the main channel for H(2)O(2) is only slightly larger than that for H(2)O. Alternative channels connecting the heme group with the monomer interface and the NADP(H) binding site are detected. These secondary channels might be important for product release.     

A designed synthetic analogue of 4-OT is specific for a non-natural substrate

The substrate specificity of 4-oxalocrotonate tautomerase (4-OT) is characterized by electrostatic interactions between positively charged arginine (Arg) side chains on the enzyme and the dianionic substrate, 4-oxalocrotonate. To generate specific hydrogen-bonding interactions with a monoanionic substrate analogue, we have introduced a urea functional group into the active site by replacing arginine side chains with isosteric citrulline (Cit) residues. This design was based on the complementarity between the urea functionality of citrulline and the uncharged amide function of the substrate, as opposed to the guanidinium-carboxylate electrostatic interaction between the wild-type enzyme and the natural substrate. Indeed, the synthetic (Arg39Cit)4-OT analogue catalyzed the tautomerization of the non-natural monoamide-monoacid substrate while it was a poor catalyst for the natural diacid substrate. The specificity of (Arg39Cit)4-OT for the monoamide-monoacid substrate relative to that of the diacid substrate was found to be 740-fold greater than that of the wild-type enzyme for tautomerization of the non-natural substrate as compared with the natural one. The role of electrostatic interactions in the tautomerization of the monoamide-monoacid substrate was probed in detail with several other Arg to Cit analogues of this enzyme. This study has demonstrated that chemical manipulation of the functional groups within the active site of an enzyme can modify its catalytic activity and substrate specificity in a predictable way, suggesting that the incorporation of noncoded amino acids into proteins has great promise for the development of new enzymatic mechanisms and new binding interactions.     

Multiscale simulation reveals multiple pathways for H2 and O2 transport in a [NiFe]-hydrogenase

Hydrogenases are enzymes that catalyze the reversible conversion of hydrogen molecules to protons and electrons. The mechanism by which the gas molecules reach the active site is important for understanding the function of the enzyme and may play a role in the selectivity for hydrogen over inhibitor molecules. Here, we develop a general multiscale molecular simulation approach for the calculation of diffusion rates and determination of pathways by which substrate or inhibitor gases can reach the protein active site. Combining kinetic data from both equilibrium simulations and enhanced sampling, we construct a master equation describing the movement of gas molecules within the enzyme. We find that the time-dependent gas population of the active site can be fit to the same phenomenological rate law used to interpret experiments, with corresponding diffusion rates in very good agreement with experimental data. However, in contrast to the conventional picture, in which the gases follow a well-defined hydrophobic tunnel, we find that there is a diverse network of accessible pathways by which the gas molecules can reach the active site. The previously identified tunnel accounts for only about 60% of the total flux. Our results suggest that the dramatic decrease in the diffusion rate for mutations involving the residue Val74 could be in part due to the narrowing of the passage Val74-Arg476, immediately adjacent to the binding site, explaining why mutations of Leu122 had only a negligible effect in experiment. Our method is not specific to the [NiFe]-hydrogenase and should be generally applicable to the transport of small molecules in proteins.     

Co-C bond activation in methylmalonyl-CoA mutase by stabilization of the post-homolysis product Co2+ cobalamin

Despite decades of research, the mechanism by which coenzyme B12 (adenosylcobalamin, AdoCbl)-dependent enzymes promote homolytic cleavage of the cofactor's Co-C bond to initiate catalysis has continued to elude researchers. In this work, we utilized magnetic circular dichroism spectroscopy to explore how the electronic structure of the reduced B12 cofactor (i.e., the post-homolysis product Co2+ Cbl) is modulated by the enzyme methylmalonyl-CoA mutase. Our data reveal a fairly uniform stabilization of the Co 3d orbitals relative to the corrin pi/pi*-based molecular orbitals when Co2+ Cbl is bound to the enzyme active site, particularly in the presence of substrate. Contrastingly, our previous studies (Brooks, A. J.; Vlasie, M.; Banerjee, R.; Brunold, T. C. J. Am. Chem. Soc. 2004, 126, 8167-8180.) showed that when AdoCbl is bound to the MMCM active site, no enzymatic perturbation of the Co3+ Cbl electronic structure occurs, even in the presence of substrate (analogues). Collectively, these observations provide direct evidence that enzymatic Co-C bond activation involves stabilization of the post-homolysis product, Co2+ Cbl, rather than destabilization of the Co3+ Cbl "ground" state.