Energetic and stereochemical effects of the protein environment on substrate: a theoretical study of methylmalonyl-CoA mutase

Functional characterization and categorization of missense mutations that cause methylmalonyl-CoA mutase (MUT) deficiency

Methylmalonyl-CoA mutase (MUT) is an essential enzyme in propionate catabolism that requires adenosylcobalamin as a cofactor. Almost 250 inherited mutations in the MUT gene are known to cause the devastating disorder methylmalonic aciduria; however, the mechanism of dysfunction of these mutations, more than half of which are missense changes, has not been thoroughly investigated. Here, we examined 23 patient missense mutations covering a spectrum of exonic/structural regions, clinical phenotypes, and ethnic populations in order to determine their influence on protein stability, using two recombinant expression systems and a thermostability assay, and enzymatic function by measuring MUT activity and affinity for its cofactor and substrate. Our data stratify MUT missense mutations into categories of biochemical defects, including (1) reduced protein level due to misfolding, (2) increased thermolability, (3) impaired enzyme activity, and (4) reduced cofactor response in substrate turnover. We further demonstrate the stabilization of wild-type and thermolabile mutants by chemical chaperones in vitro and in bacterial cells. This in-depth mutation study illustrates the tools available for MUT enzyme characterization, guides future categorization of further missense mutations, and supports the development of alternative, chaperone-based therapy for patients not responding to current treatment.

Role of vitamin B12 on methylmalonyl-CoA mutase activity

Vitamin B(12) is an organometallic compound with important metabolic derivatives that act as cofactors of certain enzymes, which have been grouped into three subfamilies depending on their cofactors. Among them, methylmalonyl-CoA mutase (MCM) has been extensively studied. This enzyme catalyzes the reversible isomerization of L-methylmalonyl-CoA to succinyl-CoA using adenosylcobalamin (AdoCbl) as a cofactor participating in the generation of radicals that allow isomerization of the substrate. The crystal structure of MCM determined in Propionibacterium freudenreichii var. shermanii has helped to elucidate the role of this cofactor AdoCbl in the reaction to specify the mechanism by which radicals are generated from the coenzyme and to clarify the interactions between the enzyme, coenzyme, and substrate. The existence of human methylmalonic acidemia (MMA) due to the presence of mutations in MCM shows the importance of its role in metabolism. The recent crystallization of the human MCM has shown that despite being similar to the bacterial protein, there are significant differences in the structural organization of the two proteins. Recent studies have identified the involvement of an accessory protein called MMAA, which interacts with MCM to prevent MCM's inactivation or acts as a chaperone to promote regeneration of inactivated enzyme. The interdisciplinary studies using this protein as a model in different organisms have helped to elucidate the mechanism of action of this isomerase, the impact of mutations at a functional level and their repercussion in the development and progression of MMA in humans. It is still necessary to study the mechanisms involved in more detail using new methods.

Combined spectroscopic/computational studies of vitamin B12 precursors: geometric and electronic structures of cobinamides

Vitamin B(12) (cyanocobalamin) and its biologically active derivatives, methylcobalamin and adenosylcobalamin, are members of the family of corrinoids, which also includes cobinamides. As biological precursors to cobalamins, cobinamides possess the same structural core, consisting of a low-spin Co(3+) ion that is ligated equatorially by the four nitrogens of a highly substituted tetrapyrrole macrocycle (the corrin ring), but differ with respect to the lower axial ligation. Specifically, cobinamides possess a water molecule instead of the nucleotide loop that coordinates axially to Co(3+)cobalamins via its dimethylbenzimidazole (DMB) base. Compared to the cobalamin species, cobinamides have proven much more difficult to study experimentally, thus far eluding characterization by X-ray crystallography. In this study, we have utilized combined quantum mechanics/molecular mechanics (QM/MM) computations to generate complete structural models of a representative set of cobinamide species with varying upper axial ligands. To validate the use of this approach, analogous QM/MM geometry optimizations were carried out on entire models of the cobalamin counterparts for which high-resolution X-ray structural data are available. The accuracy of the cobinamide structures was assessed further by comparing electronic absorption spectra computed using time-dependent density functional theory to those obtained experimentally. Collectively, the results obtained in this study indicate that the DMB → H(2)O lower axial ligand switch primarily affects the energies of the Co 3d(z(2))-based molecular orbital (MO) and, to a lesser extent, the other Co 3d-based MOs as well as the corrin π-based highest energy MO. Thus, while the energy of the lowest-energy electronic transition of cobalamins changes considerably as a function of the upper axial ligand, it is nearly invariant for the cobinamides.

Integrating computational methods to retrofit enzymes to synthetic pathways

Microbial production of desired compounds provides an efficient framework for the development of renewable energy resources. To be competitive to traditional chemistry, one requirement is to utilize the full capacity of the microorganism to produce target compounds with high yields and turnover rates. We use integrated computational methods to generate and quantify the performance of novel biosynthetic routes that contain highly optimized catalysts. Engineering a novel reaction pathway entails addressing feasibility on multiple levels, which involves handling the complexity of large-scale biochemical networks while respecting the critical chemical phenomena at the atomistic scale. To pursue this multi-layer challenge, our strategy merges knowledge-based metabolic engineering methods with computational chemistry methods. By bridging multiple disciplines, we provide an integral computational framework that could accelerate the discovery and implementation of novel biosynthetic production routes. Using this approach, we have identified and optimized a novel biosynthetic route for the production of 3HP from pyruvate.

The fragmentation-recombination mechanism of the enzyme glutamate mutase studied by QM/MM simulations

The radical mechanism of the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is studied with quantum mechanical/molecular mechanical (QM/MM) simulations based on density functional theory (DFT/MM). The hydrogen transfer between the substrate and the cofactor is found to be rate limiting with a barrier of 101.1 kJ mol(-1). A careful comparison to the uncatalyzed reaction in water is performed. The protein influences the reaction predominantly electrostatically and to a lesser degree sterically. Our calculations shed light on the atomistic details of the reaction mechanism. The well-known arginine claw and Glu 171 ( Clostridium cochlearium notation) are found to have the strongest influence on the reaction. However, a catalytic role of Glu 214, Lys 322, Gln 147, Glu 330, Lys 326, and Met 294 is found as well. The arginine claw keeps the intermediates in place and is probably responsible for the enantioselectivity. Glu 171 temporarily accepts a proton from the glutamyl radical intermediate and donates it back at the end of the reaction. We relate our results to experimental data when available. Our simulations lead to further understanding of how glutamate mutase catalyzes the carbon skeleton rearrangement of glutamate.

DFT and ONIOM(DFT:MM) studies on Co-C bond cleavage and hydrogen transfer in B12-dependent methylmalonyl-CoA mutase. Stepwise or concerted mechanism?

The considerable protein effect on the homolytic Co-C bond cleavage to form the 5'-deoxyadenosyl (Ado) radical and cob(II)alamin and the subsequent hydrogen transfer from the methylmalonyl-CoA substrate to the Ado radical in the methylmalonyl-CoA mutase (MMCM) have been extensively studied by DFT and ONIOM(DFT/MM) methods. Several quantum models have been used to systematically study the protein effect. The calculations have shown that the Co-C bond dissociation energy is very much reduced in the protein, compared to that in the gas phase. The large protein effect can be decomposed into the cage effect, the effect of coenzyme geometrical distortion, and the protein MM effect. The largest contributor is the MM effect, which mainly consists of the interaction of the QM part of the coenzyme with the MM part of the coenzyme and the surrounding residues. In particular, Glu370 plays an important role in the Co-C bond cleavage process. These effects tremendously enhance the stability of the Co-C bond cleavage state in the protein. The initial Co-C bond cleavage and the subsequent hydrogen transfer were found to occur in a stepwise manner in the protein, although the concerted pathway for the Co-C bond cleavage coupled with the hydrogen transfer is more favored in the gas phase. The assumed concerted transition state in the protein has more deformation of the coenzyme and the substrate and has less interaction with the protein than the stepwise route. Key factors and residues in promoting the enzymatic reaction rate have been discussed in detail.

Quantum catalysis in B12-dependent methylmalonyl-CoA mutase: experimental and computational insights

B12-dependent methylmalonyl-CoA mutase catalyses the interchange of a hydrogen atom and the carbonyl-CoA group on adjacent carbons of methylmalonyl-CoA to give the rearranged product, succinyl-CoA. The first step in this reaction involves the transient generation of cofactor radicals by homolytic rupture of the cobalt-carbon bond to generate the deoxyadenosyl radical and cob(II)alamin. This step exhibits a curious sensitivity to isotopic substitution in the substrate, methylmalonyl-CoA, which has been interpreted as evidence for kinetic coupling. The magnitude of the isotopic discrimination is large and a deuterium isotope effect ranging from 35.6 at 20 degrees C to 49.9 at 5 degrees C has been recorded. Arrhenius analysis of the temperature dependence of this isotope effect provides evidence for quantum tunnelling in this hydrogen transfer step. The mechanistic complexity of the observed rate constant for cobalt-carbon bond homolysis together with the spectroscopically silent nature of many of the component steps limits the insights that can be derived by experimental approaches alone. Computational studies using a newly developed geometry optimization scheme that allows determination of the transition state in the full quantum mechanical/molecular mechanical coordinate space have yielded novel insights into the strategy deployed for labilizing the cobalt-carbon bond and poising the resulting deoxyadenosyl radical for subsequent hydrogen atom abstraction.

Computational study on mechanistic details of the aminoethanol rearrangement catalyzed by the vitamin B12-dependent ethanolamine ammonia lyase: His and Asp/Glu acting simultaneously as catalytic auxiliaries

The rearrangement of aminoethanol catalyzed by ethanolamine ammonia lyase is investigated by computational means employing DFT (B3LYP/6-31G) and ab initio molecular orbital theory (QCISD/cc-pVDZ). The study aims at providing a detailed account on various crucial aspects, in particular a distinction between a direct intramolecular migration of the partially protonated NH(2) group vs elimination of NH(4)(+). Three mechanistic scenarios were explored: (i) According to the calculations, irrespective of the nature of the protonating species, intramolecular migration of the NH(3) group is energetically less demanding than elimination of NH(4)(+). However, all computed activation enthalpies exceed the experimentally derived activation enthalpy (15 kcal/mol) associated with the rate-determining step, i.e., the hydrogen abstraction from the 5'-deoxyadenosine by the product radical. For example, when imidazole is used as a model system for His interacting with the NH(3) group of the substrate, the activation enthalpy for the migration process amounts to 27.4 kcal/mol. If acetic acid is employed to mimic Asp or Glu, the activation enthalpy is somewhat lower, being equal to 24.2 kcal/mol. (ii) For a partial deprotonation of the substrate 2 at the OH group, the rearrangement mechanism consists of the dissociation of an NH(2) radical from C(2) and its association at C(1) atom. For all investigated proton acceptors (i.e., OH(-), HCOO(-), CH(3)COO(-), CH(2)NH, imidazole), the activation enthalpy for the dissociation step also exceeds 15 kcal/mol. Typical data are 20.2 kcal/mol for Ac(-) and 23.8 kcal/mol for imidazole. (iii) However, in a synergistic action of partial protonation of the NH(2) group and partial deprotonation of the OH group by the two conceivable catalytic auxiliaries Asp/Glu and His, the activation enthalpy computed is compatible with the experimental data. For imidazole and acetate as model systems, the activation enthalpy is equal to 13.7 kcal/mol. This synergistic action of the two catalytic groups is expected to take place in a physiologically realistic pH range of 6-9.5, and the present computational findings may help to further characterize the yet unknown structural details of the ethanolamine ammonia lyase's active site.