Mechanism of mycolic acid cyclopropane synthase: a theoretical study

The His23 and Lys79 pair determines the high catalytic efficiency of the inorganic pyrophosphatase of the haloacid dehalogenase superfamily

Haloacid dehalogenase (HAD) superfamily members are mainly phosphomonoesterases, while BT2127 from Bacteroides thetaiotaomicron of the HAD superfamily is identified as an inorganic pyrophosphatase. In this study, to explore the roles of the Lys79 and His23 pair in the hydrolysis reaction of inorganic pyrophosphate (PP_i) catalyzed by BT2127, a series of models were designed. Calculations were performed by using the density functional theory (DFT) method with the dispersion energy D3-B3LYP. The His23 and Lys79 pair plays a key role in the high catalytic efficiency of BT2127 with PPi. First, the His23 and Lys79 pair prompts Asp13 to easily provide a proton to the leaving group, which remarkably reduces the energy barrier of the phospho-transfer step; then, Lys79 provides a proton to the first leaving phosphate group via His23, produces a more electrically stabilized phosphate (H₃PO₄), makes this step exothermal, and further promotes the subsequent phospho-enzyme intermediate hydrolysis. The results suggest that the Lys79-His23 pair helps BT2127 reach high catalytic efficiency by strengthening the acid catalysis. Our study provides detailed chemical insights into the evolution of the inorganic pyrophosphatase function of BT2127 from the phosphomonoesterase of the HAD superfamily and the biomimetic enzyme design.

Diversity of the reaction mechanisms of SAM-dependent enzymes

S-adenosylmethionine (SAM) is ubiquitous in living organisms and is of great significance in metabolism as a cofactor of various enzymes. Methyltransferases (MTases), a major group of SAM-dependent enzymes, catalyze methyl transfer from SAM to C, O, N, and S atoms in small-molecule secondary metabolites and macromolecules, including proteins and nucleic acids. MTases have long been a hot topic in biomedical research because of their crucial role in epigenetic regulation of macromolecules and biosynthesis of natural products with prolific pharmacological moieties. However, another group of SAM-dependent enzymes, sharing similar core domains with MTases, can catalyze nonmethylation reactions and have multiple functions. Herein, we mainly describe the nonmethylation reactions of SAM-dependent enzymes in biosynthesis. First, we compare the structural and mechanistic similarities and distinctions between SAM-dependent MTases and the non-methylating SAM-dependent enzymes. Second, we summarize the reactions catalyzed by these enzymes and explore the mechanisms. Finally, we discuss the structural conservation and catalytical diversity of class I-like non-methylating SAM-dependent enzymes and propose a possibility in enzymes evolution, suggesting future perspectives for enzyme-mediated chemistry and biotechnology, which will help the development of new methods for drug synthesis.

Sensory Perception of Non-Deuterated and Deuterated Organic Compounds

The chemical background of olfactory perception has been subject of intensive research, but no available model can fully explain the sense of smell. There are also inconsistent results on the role of the isotopology of molecules. In experiments with human subjects it was found that the isotope effect is weak with acetone and D6 -acetone. In contrast, clear differences were observed in the perception of octanoic acid and D15 -octanoic acid. Furthermore, a trained sniffer dog was initially able to distinguish between these isotopologues of octanoic acid. In chromatographic measurements, the respective deuterated molecule showed weaker interaction with a non-polar liquid phase. Quantum chemical calculations give evidence that deuterated octanoic acid binds more strongly to a model receptor than non-deuterated. In contrast, the binding of the non-deuterated molecule is stronger with acetone. The isotope effect is calculated in the framework of statistical mechanics. It results from a complicated interplay between various thermostatistical contributions to the non-covalent free binding energies and it turns out to be very molecule-specific. The vibrational terms including non-classical zero-point energies play about the same role as rotational/translational contributions and are larger than bond length effects for the differential isotope perception of odor for which general rules cannot be derived.

The Effect of Deuteration on the H2 Receptor Histamine Binding Profile: A Computational Insight into Modified Hydrogen Bonding Interactions

We used a range of computational techniques to reveal an increased histamine affinity for its H₂ receptor upon deuteration, which was interpreted through altered hydrogen bonding interactions within the receptor and the aqueous environment preceding the binding. Molecular docking identified the area between third and fifth transmembrane α-helices as the likely binding pocket for several histamine poses, with the most favorable binding energy of −7.4 kcal mol⁻¹ closely matching the experimental value of −5.9 kcal mol⁻¹. The subsequent molecular dynamics simulation and MM-GBSA analysis recognized Asp98 as the most dominant residue, accounting for 40% of the total binding energy, established through a persistent hydrogen bonding with the histamine −NH₃⁺ group, the latter further held in place through the N–H∙∙∙O hydrogen bonding with Tyr250. Unlike earlier literature proposals, the important role of Thr190 is not evident in hydrogen bonds through its −OH group, but rather in the C–H∙∙∙π contacts with the imidazole ring, while its former moiety is constantly engaged in the hydrogen bonding with Asp186. Lastly, quantum-chemical calculations within the receptor cluster model and utilizing the empirical quantization of the ionizable X–H bonds (X = N, O, S), supported the deuteration-induced affinity increase, with the calculated difference in the binding free energy of −0.85 kcal mol⁻¹, being in excellent agreement with an experimental value of −0.75 kcal mol⁻¹, thus confirming the relevance of hydrogen bonding for the H₂ receptor activation.

Designing novel inhibitors against cyclopropane mycolic acid synthase 3 (PcaA): targeting dormant state of Mycobacterium tuberculosis

Mycobacterium tuberculosis can sustain inside the host in dormant (non-replicating) state for years. It suppresses the host immune system by residing in the host alveolar macrophage, resulting in the development of latent tuberculosis. Despite many antibiotics available for the treatment of tuberculosis, the major hurdle in complete elimination is the ability of the bacilli to undergo dormancy and develop resistance against the existing drugs. Cyclopropanation of mycolic acids present in the cell wall of mycobacteria is required for its persistence and virulence. Cyclopropane synthases such as PcaA, CmaA1 and CmaA2, introduce site-specific modifications in mycolic acids. PcaA expression levels are high during dormancy and the gene mutants fails to persist, showing reduced survival in host macrophage. Hence, PcaA appears as a potential target to develop inhibitors against the dormant bacilli. In this study, we have identified compounds with maximum binding affinity against PcaA by in-silico virtual screening of anti-tuberculosis compounds and their structural analogues. In-silico docking followed molecular dynamic simulations and free energy calculations of the compounds with highest docking score in their respective libraries. This study reports novel inhibitors that can act as better anti-tuberculosis compounds targeting PcaA.Communicated by Ramaswamy H. Sarma.

Relevance of Hydrogen Bonds for the Histamine H2 Receptor-Ligand Interactions: A Lesson from Deuteration

We used a combination of density functional theory (DFT) calculations and the implicit quantization of the acidic N–H and O–H bonds to assess the effect of deuteration on the binding of agonists (2-methylhistamine and 4-methylhistamine) and antagonists (cimetidine and famotidine) to the histamine H2 receptor. The results show that deuteration significantly increases the affinity for 4-methylhistamine and reduces it for 2-methylhistamine, while leaving it unchanged for both antagonists, which is found in excellent agreement with experiments. The revealed trends are interpreted in the light of the altered strength of the hydrogen bonding upon deuteration, known as the Ubbelohde effect, which affects ligand interactions with both active sites residues and solvent molecules preceding the binding, thus providing strong evidence for the relevance of hydrogen bonding for this process. In addition, computations further underline an important role of the Tyr250 residue for the binding. The obtained insight is relevant for the therapy in the context of (per)deuterated drugs that are expected to enter therapeutic practice in the near future, while this approach may contribute towards understanding receptor activation and its discrimination between agonists and antagonists.

A radical S-adenosyl-L-methionine enzyme and a methyltransferase catalyze cyclopropane formation in natural product biosynthesis

Cyclopropanation of unactivated olefinic bonds via addition of a reactive one-carbon species is well developed in synthetic chemistry, whereas natural cyclopropane biosynthesis employing this strategy is very limited. Here, we identify a two-component cyclopropanase system, composed of a HemN-like radical S-adenosyl-L-methionine (SAM) enzyme C10P and a methyltransferase C10Q, catalyzes chemically challenging cyclopropanation in the antitumor antibiotic CC-1065 biosynthesis. C10P uses its [4Fe-4S] cluster for reductive cleavage of the first SAM to yield a highly reactive 5'-deoxyadenosyl radical, which abstracts a hydrogen from the second SAM to produce a SAM methylene radical that adds to an sp²-hybridized carbon of substrate to form a SAM-substrate adduct. C10Q converts this adduct to CC-1065 via an intramolecular S_N2 cyclization mechanism with elimination of S-adenosylhomocysteine. This cyclopropanation strategy not only expands the enzymatic reactions catalyzed by the radical SAM enzymes and methyltransferases, but also sheds light on previously unnoticed aspects of the versatile SAM-based biochemistry.

Structural and mechanistic comparison of the Cyclopropane Mycolic Acid Synthases (CMAS) protein family of Mycobacterium tuberculosis

Tuberculosis (TB) is a chronic disease caused by the bacillus Mycobacterium tuberculosis(Mtb) and remains a leading cause of mortality worldwide. The bacteria has an external wall which protects it from being killed, and the enzymes involved in the biosynthesis of the cell wall components have been proposed as promising targets for future drug development efforts. Cyclopropane Mycolic Acid Synthases (CMAS) constitute a group of ten homologous enzymes which belong to the mycolic acid biosynthesis pathway. These enzymes have S-adenosyl-l-methionine (SAM) dependent methyltransferase activity with a peculiarity, each one of them has strong substrate selectivity and reaction specificity, being able to produce among other things cyclopropanes or methyl-alcohol groups from the lipid olefin group. How each CMAS processes its substrate and how the specificity and selectivity are encoded in the protein sequence and structure, is still unclear. In this work, by using a combination of modeling tools, including comparative modeling, docking, all-atom MD and QM/MM methodologies we studied in detail the reaction mechanism of cmaA2, mmaA4, and mmaA1 CMAS and described the molecular determinants that lead to different products. We have modeled the protein-substrate complex structure and determined the free energy pathway for the reaction. The combination of modeling tools at different levels of complexity allows having a complete picture of the CMAS structure-activity relationship.

The Quantum Nature of Drug-Receptor Interactions: Deuteration Changes Binding Affinities for Histamine Receptor Ligands

In this article we report a combined experimental and computational study concerning the effects of deuteration on the binding of histamine and two other histaminergic agonists to 3H-tiotidine-labeled histamine H2 receptor in neonatal rat astrocytes. Binding affinities were measured by displacing radiolabeled tiotidine from H2 receptor binding sites present on cultured neonatal rat astrocytes. Quantum-chemical calculations were performed by employing the empirical quantization of nuclear motion within a cluster model of the receptor binding site extracted from the homology model of the entire H2 receptor. Structure of H2 receptor built by homology modelling is attached in the supporting information (S1 Table) Experiments clearly demonstrate that deuteration affects the binding by increasing the affinity for histamine and reducing it for 2-methylhistamine, while basically leaving it unchanged for 4-methylhistamine. Ab initio quantum-chemical calculations on the cluster system extracted from the homology H2 model along with the implicit quantization of the acidic N-H and O-H bonds demonstrate that these changes in the binding can be rationalized by the altered strength of the hydrogen bonding upon deuteration known as the Ubbelohde effect. Our computational analysis also reveals a new mechanism of histamine binding, which underlines an important role of Tyr250 residue. The present work is, to our best knowledge, the first study of nuclear quantum effects on ligand receptor binding. The ligand H/D substitution is relevant for therapy in the context of perdeuterated and thus more stable drugs that are expected to enter therapeutic practice in the near future. Moreover, presented approach may contribute towards understanding receptor activation, while a distant goal remains in silico discrimination between agonists and antagonists based on the receptor structure.

Theoretically exploring the key role of the Lys412 residue in the conversion of N2O to N2by nitrous oxide reductase from Achromobacter cycloclastes

The active Cu_Zcluster of NOR provides strong back-donation to coordinated N₂O and activates the O atom of the N₂O group facilitating H-bonding and protonationviathe Lys412 residue.

Theoretical study on the chemical mechanism of enoyl-CoA hydratase and the form of inhibitor binding

Enoyl-CoA hydratase (ECH) catalyzes the second step in the vital β-oxidation pathway of fatty acid metabolism. This enzyme catalyzes the syn-addition of a water molecule across the double bond of 4-(N,N-dimethylamino) cinnamoyl-CoA (DAC-CoA). In this work, the reaction mechanisms of ECH were investigated using the density functional theory (DFT) methods. The different protonation states in which the important residues Glu164 and Glu144 are either neutral or ionized were considered. Four models of the active site were designed based on the X-ray crystal structure of the enzyme. The calculations gave strong support to the proposed mechanism and confirmed that both Glu164 and Glu144 are in a deprotonated state in the reaction mechanism of ECH. In addition, we constructed a model of the active site with the inhibitor acetoacetyl-CoA based on the crystal structure. Caomparison of the calculated energy barriers showed that binding of the keto-enol form of the inhibitor is more reasonable than that of the di-keto form in the inhibition process. Moreover, acetoacetyl-CoA was found to exhibit a keto-enol tautomerism when it acts as an inhibitor in the reaction. The present theoretical results indicated that both residues Glu164 and Glu144 are unprotonated in ECH with the substrate bound, while only Glu164 is unprotonated when the inhibitor binds ECH.

Characterization of yatakemycin gene cluster revealing a radical S-adenosylmethionine dependent methyltransferase and highlighting spirocyclopropane biosynthesis

Yatakemycin (YTM), an antitumor natural product, represents the most potent member of a class of potent anticancer natural products including CC-1065 and duocarmycins. Herein we describe the biosynthetic gene cluster of YTM, which was identified by genome scanning of Streptomyces sp. TP-A0356. This cluster consists of 31 open reading frames (ORFs) and was localized to a 36 kb DNA segment. Moreover, its involvement in YTM biosynthesis was confirmed by cluster deletion, gene replacement, and complementation. Inactivation of ytkT, which encodes a radical S-adenosylmethionine (SAM) protein, created a mutant strain that failed to produce YTM but accumulated a new metabolite, which was structurally elucidated as a precursor that was related to the formation of the cyclopropane ring. More importantly, biochemical characterization of the radical SAM-dependent enzyme YtkT revealed that it is a novel C-methyltransferase and contributes to an advanced intermediate during formation of the cyclopropane ring through a radical mechanism in the YTM biosynthetic pathway. On the basis of in silico analysis, genetic experiments, structure elucidation of the novel intermediate, and biochemical characterization, a biosynthetic pathway for yatakemycin was proposed, which sets the stage to further investigate the novel enzymatic mechanisms and engineer the biosynthetic machinery for the production of novel analogues.

Structure of geranyl diphosphate C-methyltransferase from Streptomyces coelicolor and implications for the mechanism of isoprenoid modification

Geranyl diphosphate C-methyltransferase (GPPMT) from Streptomyces coelicolor A3(2) is the first methyltransferase discovered that modifies an acyclic isoprenoid diphosphate, geranyl diphosphate (GPP), to yield a noncanonical acyclic allylic diphosphate product, 2-methylgeranyl diphosphate, which serves as the substrate for a subsequent cyclization reaction catalyzed by a terpenoid cyclase, methylisoborneol synthase. Here, we report the crystal structures of GPPMT in complex with GPP or the substrate analogue geranyl S-thiolodiphosphate (GSPP) along with S-adenosyl-L-homocysteine in the cofactor binding site, resulting from in situ demethylation of S-adenosyl-L-methionine, at 2.05 or 1.82 Å resolution, respectively. These structures suggest that both GPP and GSPP can undergo catalytic methylation in crystalline GPPMT, followed by dissociation of the isoprenoid product. S-Adenosyl-L-homocysteine remains bound in the active site, however, and does not exchange with a fresh molecule of cofactor S-adenosyl-L-methionine. These structures provide important clues about the molecular mechanism of the reaction, especially with regard to the face of the 2,3 double bond of GPP that is methylated as well as the stabilization of the resulting carbocation intermediate through cation-π interactions.