Formation Of Tetrahedral Intermediate Research Articles

New approach to the equilibrium of D-glucaric acid in water: DFT study of tetrahedral intermediates

The present study aims at investigating, for the first time, the formation of tetrahedral intermediates during the equilibrium mechanism of D-glucaric acid in water. In this work, carried out at B3LYP/6-31+G(d,p) computational level in vacuum and using the polarizable continuum model, two systems are proposed and simulated: the first describes an intramolecular proton transfer, and the second requires the intervention of solvent molecules for the proton transfer. The second system, simulated in vacuum, indicates greater ease in establishing this equilibrium due to a significant reduction in angular tensions. Also, the idea of a tetrahedral intermediate carrying two geminal diol groups constitutes a favorable hypothesis due to the low energy barriers calculated as well as a better structural flexibility generated following the equilibrium of the tetrahedral intermediates between them. Experimental confirmation of the presence of such a powerful tetrahedral intermediate would be a considerable advance in resolving the equilibrium mechanisms of D-glucaric acid.

Decarboxylation Mechanism of iso-Orotate Decarboxylase Revisited.

iso-Orotate decarboxylase (IDCase), which is involved in the thymidine salvage pathway, has attracted considerable interest owing to its chemical similarity to a hypothetical DNA decarboxylase in mammals. Although valuable insights into the active DNA demethylation of 5-methyl-cytosine can be obtained from the decarboxylation mechanism of 5-carboxyl-uracil (5caU) catalyzed by IDCase, this mechanism remains under debate. In this study, the catalytic mechanism of 5caU decarboxylation by IDCase was studied using hybrid quantum mechanics/molecular mechanics (QM/MM) methodologies and density functional theory (DFT) calculations with a truncated model. The calculations supported a mechanism involving three sequential stages: activation of the 5caU substrate via proton transfer from an arginine (R262') to the carboxyl group of 5caU, formation of a tetrahedral intermediate, and decarboxylation of the tetrahedral intermediate to generate uracil as the product. The reaction pathways and structures obtained using the QM/MM and DFT methods coincided with each other. These simulations provided detailed insights into the unique mechanism of IDCase, clarifying various unresolved issues, such as the critical role of R262'. In addition, aspartate D323 was found to act as a general base in the tetrahedral intermediate formation step and a general acid in the later C-C bond cleavage step.

One-Pot, Telescoped Alkenylation of Amides via Stable Tetrahedral Intermediates as Lithium Enolate Precursors.

A mild and efficient telescoped procedure for the stereoselective alkenylation of simple, non-activated amides using LiCH2SiMe3 and carbonyl compounds as surrogates of alkenyllithium reagents is reported. Our methodology relies on the formation of stable tetrahedral intermediates, which, upon collapse into highly reactive lithium enolates in a solvent-dependent fashion, allows for the assembly of α,β-unsaturated ketones in a single synthetic operation with high stereoselectivity.

Synthesis and Application of Constrained Amidoboronic Acids Using Amphoteric Boron-Containing Building Blocks

Amidoboronic acid-containing peptidomimetics are an important class of scaffolds in chemistry and drug discovery. Despite increasing interest in boron-based enzyme inhibitors, constrained amidoboronic acids have received little attention due to the limited options available for their synthesis. We describe a new methodology to prepare both α- and β-amidoboronic acids that impose restrictions on backbone angles. Lewis acid-promoted Boyer-Schmidt-Aube lactam ring expansions using an azidoalkylboronate enabled generation of constrained α-amidoboronic acid derivatives, whereas assembly of the homologous β-amidoboronic acids was achieved through a novel boronic acid-mediated lactamization process stemming from an α-boryl aldehyde. The results of quantum chemical calculations suggest carboxylate-boron coordination to be rate-limiting for small ring sizes, whereas the tetrahedral intermediate formation is rate limiting in the case of larger rings. As part of this study, an application of β-amidoboronic acid derivatives as novel VIM-2 metallo-β-lactamase inhibitors has been demonstrated.

Structure-based virtual screening suggests inhibitors of 3-Chymotrypsin-Like Protease of SARS-CoV-2 from Vernonia amygdalina and Occinum gratissimum.

Antiviral culinary plants are potential bioresources for preventive nutraceuticals and/or antiviral drugs in COVID-19. Structure-based virtual screening was undertaken to screen 173 compounds previously reported from Vernonia amygdalina and Occinum gratissimum for direct interaction with the active site of the 3-Chymotrypsin-Like Protease (3CLpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Based on docking scores and comparison with reference inhibitors, a hit-list of 10 top phytocompounds was defined, which also had strong interactions with the catalytic centre of 3CLpro from three related strains of coronavirus (SARS-CoV, MERS-CoV, HKU4). Among these, six compounds (neoandrographolide, vernolide, isorhamnetin, chicoric acid, luteolin, and myricetin) exhibited the highest binding tendencies to the equilibrated conformers of SARS-CoV-2 3CLpro in an in-depth docking analysis to 5 different representative conformations from the cluster analysis of the molecular dynamics simulation (MDS) trajectories of the protein. In silico drug-likeness analyses revealed two drug-like terpenoids viz: neoandrographolide and vernolide as promising inhibitors of SARS-CoV-2 3CLpro. These structures were accommodated within the substrate-binding pocket; and interacted with the catalytic dyad (Cys145 and His41), the oxyanion loop (residues 138–145), and the S1/S2 sub-sites of the enzyme active site through the formation of an array of hydrogen bonds and hydrophobic interactions. Molecular dynamics simulation and binding free energy calculation revealed that the terpenoid-enzyme complexes exhibit strong interactions and structural stability. Therefore, these compounds may stabilize the conformation of the flexible oxyanion loop; and thereby interfere with the tetrahedral oxyanion intermediate formation during the proteolytic activity of the enzyme.

Hydrolytic deamination mechanisms of guanosine monophosphate: A computational study

The computational importance of guanosine monophosphate (GMP) comes from its importance in the synthesis of nucleic acids. Deamination reaction mechanisms, kinetics, and thermodynamics of GMP with 3H2O and 2H2O/OH− to produce xanthosine monophosphate (XMP) have been investigated. The proposed reactions were optimized in the gas-phase and in an aqueous medium. Two key steps were found for the proposed pathways: a tetrahedral intermediate formation and a single water-mediated six-membered ring transition state, yielding the desired product. The deamination of GMP with 2H2O/OH− yield deprotonated XMP− with only one pathway, which is highly exothermic, and the activation energy is significantly lower (33 kJ mol−1 at M11/6-31G(d)) compared to the other pathways. Our results suggest that pathway C is the most kinetically and thermodynamically favorable mechanism for the deamination of GMP with 3H2O with the lowest energy barrier of 123 kJ mol−1, due to multiple functions of water that enhance the reaction profile.

Transition-Metal-Free Activation of Amides by N-C Bond Cleavage.

The amide bond N-C activation represents a powerful strategy in organic synthesis to functionalize the historically inert amide linkage. This personal account highlights recent remarkable advances in transition-metal-free activation of amides by N-C bond cleavage, focusing on both (1) mechanistic aspects of ground-state-destabilization of the amide bond enabling formation of tetrahedral intermediates directly from amides with unprecedented selectivity, and (2) synthetic utility of the developed transformations. Direct nucleophilic addition to amides enables a myriad of powerful methods for the formation of C-C, C-N, C-O and C-S bonds, providing a straightforward and more synthetically useful alternative to acyl-metals.

First Step of the Transglutaminase Reaction Catalyzed by Activated Factor XIII Subunit A, Hybrid Quantum Chemistry/Molecular Mechanics Calculations.

The A subunit of blood coagulation factor XIII belongs to the family of transglutaminase enzymes. Its active form (FXIIIa) catalyzes isopeptide bond formation between Glu and Lys residues of specific substrates. Little data are available on the mechanism of this reaction. In this work, the first step of the proposed two-step process was investigated using two different protocols of hybrid QM/molecular mechanics (MM) calculations: an ONIOM-based model as well as QM/MM/molecular dynamics (MD) metadynamics simulations in explicit TIP3P solvent with Gromacs, PLUMED, and a DFTB3 package. Based on calculations involving a truncated system derived from docking of a peptide substrate, our study confirms the higher stability of a zwitterionic form of the catalytic Cys and His residues in the Michaelis complex as well as the "resting" state of the enzyme. Potential energy surfaces, obtained by geometry optimizations with Gaussian, show a two-step reaction mechanism with a zwitterionic tetrahedral intermediate formation in the first and NH3 dissociation in the second step in the case of our ONIOM system. In contrast, in QM/MM MD metadynamics simulations, all three steps occurred in a concerted manner. As a conclusion, our model is able to provide insights into the reaction mechanism of this enzyme.

Quantum Chemical Study of Addition–Elimination Reactions of Dimethyl Carbonate with Methylamine

The mechanisms of noncatalytic and autocatalytic addition–elimination reactions of dimethyl carbonate with methylamine have been studied in terms of wB97XD/6-311++G(df,p), M06/6-311++G(df,p), and PBE0/6-311++G(df,p) quantum chemical methods, and thermodynamic and activation parameters of these reactions have been calculated. The rate-determining stage is the formation of tetrahedral carbon intermediate. The addition–elimination mechanism is preferred over the autocatalytic SN2 mechanism determined by enhanced donor–acceptor and acid–base properties of amine complexes with alcohol associates in comparison to amines.

Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent.

Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.

Understanding the hydrolysis mechanism of ethyl acetate catalyzed by an aqueous molybdocene: a computational chemistry investigation.

A thoroughly mechanistic investigation on the [Cp2Mo(OH)(OH2)](+)-catalyzed hydrolysis of ethyl acetate has been performed using density functional theory methodology together with continuum and discrete-continuum solvation models. The use of explicit water molecules in the PCM-B3LYP/aug-cc-pVTZ (aug-cc-pVTZ-PP for Mo)//PCM-B3LYP/aug-cc-pVDZ (aug-cc-pVDZ-PP for Mo) computations is crucial to show that the intramolecular hydroxo ligand attack is the preferred mechanism in agreement with experimental suggestions. Besides, the most stable intermediate located along this mechanism is analogous to that experimentally reported for the norbornenyl acetate hydrolysis catalyzed by molybdocenes. The three most relevant steps are the formation and cleavage of the tetrahedral intermediate immediately formed after the hydroxo ligand attack and the acetic acid formation, with the second one being the rate-determining step with a Gibbs energy barrier of 36.7 kcal/mol. Among several functionals checked, B3LYP-D3 and M06 give the best agreement with experiment as the rate-determining Gibbs energy barrier obtained only differs 0.2 and 0.7 kcal/mol, respectively, from that derived from the experimental kinetic constant measured at 296.15 K. In both cases, the acetic acid elimination becomes now the rate-determining step of the overall process as it is 0.4 kcal/mol less stable than the tetrahedral intermediate cleavage. Apart from clarifying the identity of the cyclic intermediate and discarding the tetrahedral intermediate formation as the rate-determining step for the mechanism of the acetyl acetate hydrolysis catalyzed by molybdocenes, the small difference in the Gibbs energy barrier found between the acetic acid formation and the tetrahedral intermediate cleavage also uncovers that the rate-determining step could change when studying the reactivity of carboxylic esters other than ethyl acetate substrate specific toward molybdocenes or other transition metal complexes. Therefore, in general, the information reported here could be of interest in designing new catalysts and understanding the reaction mechanism of these and other metal-catalyzed hydrolysis reactions.

Theoretical studies on the deacylation step of acylated Candida Antarctica lipase B: Structural and reaction pathway analysis

The deacylation step of acylated Candida Antarctica lipase B, which was acylated with methylcaprylate (MEC) and acetylcholine (ACh), has been studied by using density functional theory method. Free energies of the entire reaction were calculated for enzyme deacylation by water and hydrogen peroxide that represented hydrolysis and perhydrolysis reactions, respectively. The calculations displayed that a stepwise mechanism there was with the enzyme–product complex being a deep minimum on the free energy surfaces of both of two reactions. The tetrahedral intermediate formation was the rate-determining step of all reactions, which needed 8.1 to 10.5kcalmol−1 for activation in different reactions. In the second stage of the reaction, fewer free energy barriers, between 4.7 and 5.9kcalmol−1, were identified to enable the proton transfer from His224 to Ser105 and the breakdown of the tetrahedral intermediate. These calculated activation free energies approved theoretical possibility of both of two reactions for two substrates. Finally, an applied tool examined the interactions role in the stability and energy levels of different chemical species.

Further theoretical insight into the reaction mechanism of the hepatitis C NS3/NS4A serine protease

The main reactions of the hepatitis C virus NS3/NS4A serine protease are studied using the second-order Møller–Plesset ab initio method and rather large basis sets to correct the previously reported AM1/CHARMM22 potential energy surfaces. The reaction efficiencies measured for the different substrates are explained in terms of the tetrahedral intermediate formation step (the rate-limiting process). The energies of the barrier and the corresponding intermediate are so close that the possibility of a concerted mechanism is open (especially for the NS5A/5B substrate). This is in contrast to the suggested general reaction mechanism of serine proteases, where a two-step mechanism is postulated.

Hemiacetal stabilization in a chymotrypsin inhibitor complex and the reactivity of the hydroxyl group of the catalytic serine residue of chymotrypsin

The aldehyde inhibitor Z-Ala-Ala-Phe-CHO has been synthesized and shown by 13C-NMR to react with the active site serine hydroxyl group of alpha-chymotrypsin to form two diastereomeric hemiacetals. For both hemiacetals oxyanion formation occurs with a pKa value of ~7 showing that chymotrypsin reduces the oxyanion pKa values by ~5.6 pKa units and stabilizes the oxyanions of both diastereoisomers by ~32kJmol−1. As pH has only a small effect on binding we conclude that oxyanion formation does not have a significant effect on binding the aldehyde inhibitor. By comparing the binding of Z-Ala-Ala-Phe-CHO with that of Z-Ala-Ala-Phe-H we estimate that the aldehyde group increases binding ~100 fold. At pH7.2 the effective molarity of the active site serine hydroxy group is ~6000 which is ~7× less effective than with the corresponding glyoxal inhibitor. Using 1H-NMR we have shown that at both 4 and 25°C the histidine pKa is ~7.3 in free chymotrypsin and it is raised to ~8 when Z-Ala-Ala-Phe-CHO is bound. We conclude that oxyanion formation only has a minor role in raising the histidine pKa and that the aldehyde hydrogen must be replaced by a larger group to raise the histidine pKa>10 and give stereospecific formation of tetrahedral intermediates. The results show that a large increase in the pKa of the active site histidine is not needed for the active site serine hydroxyl group to have an effective molarity of 6000.

Design, synthesis and in vitro kinetic study of tranexamic acid prodrugs for the treatment of bleeding conditions

Based on density functional theory (DFT) calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives four tranexamic acid prodrugs were designed. The DFT results on the acid catalyzed hydrolysis revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. When the amine leaving group was a primary amine or tranexamic acid moiety, the tetrahedral intermediate collapse was the rate-limiting step, whereas in the cases by which the amine leaving group was aciclovir or cefuroxime the rate-limiting step was the tetrahedral intermediate formation. The linear correlation between the calculated DFT and experimental rates for N-methylmaleamic acids 1-7 provided a credible basis for designing tranexamic acid prodrugs that have the potential to release the parent drug in a sustained release fashion. For example, based on the calculated B3LYP/6-31G(d,p) rates the predicted t1/2 (a time needed for 50% of the prodrug to be converted into drug) values for tranexamic acid prodrugs ProD 1-ProD 4 at pH 2 were 556h [50.5h as calculated by B3LYP/311+G(d,p)] and 6.2h as calculated by GGA: MPW1K), 253h, 70s and 1.7h, respectively. Kinetic study on the interconversion of the newly synthesized tranexamic acid prodrug ProD 1 revealed that the t1/2 for its conversion to the parent drug was largely affected by the pH of the medium. The experimental t1/2 values in 1N HCl, buffer pH 2 and buffer pH 5 were 54min, 23.9 and 270h, respectively.

Prodrugs for masking bitter taste of antibacterial drugs—a computational approach

DFT calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. Further, it was established that when the amine leaving group was a secondary amine, acyclovir or cefuroxime moiety the tetrahedral intermediate formation was the rate-limiting step such as in the cases of acyclovir ProD 1- ProD 4 and cefuroxime ProD 1- ProD 4. In addition, the linear correlation between the calculated and experimental rates provided a credible basis for designing prodrugs for masking bitter taste of the corresponding parental drugs which have the potential to release the parent drug in a sustained release fashion. For example, based on the DFT calculated rates the predicted t₁/₂ (a time needed for 50 % of the reactant to be hydrolyzed to products) for cefuroxime prodrugs, cefuroxime ProD 1- ProD 4, were 12 min, 18 min, 200 min and 123 min, respectively.

Importance of Tetrahedral Intermediate Formation in the Catalytic Mechanism of the Serine Proteases Chymotrypsin and Subtilisin

Two new inhibitors in which the terminal α-carboxyl groups of Z-Ala-Ala-Phe-COOH and Z-Ala-Pro-Phe-COOH have been replaced with a proton to give Z-Ala-Ala-Phe-H and Z-Ala-Pro-Phe-H, respectively, have been synthesized. Using these inhibitors, we estimate that for α-chymotrypsin and subtilisin Carlsberg the terminal carboxylate group decreases the level of inhibitor binding 3-4-fold while a glyoxal group increases the level of binding by 500-2000-fold. We show that at pH 7.2 the effective molarities of the catalytic hydroxyl group of the active site serine are 41000-229000 and 101000-159000 for α-chymotrypsin and subtilisin Carlsberg, respectively. It is estimated that oxyanion stabilization and the increased effective molarity of the catalytic serine hydroxyl group can account for the catalytic efficiency of the reaction. We argue that substrate binding induces the formation of a strong hydrogen bond or low-barrier hydrogen bond between histidine-57 and aspartate-102 that increases the pK(a) of the active site histidine, allowing it to be an effective general base catalyst for the formation of the tetrahedral intermediate and increasing the effective molarity of the catalytic hydroxyl group of serine-195. A catalytic mechanism for acyl intermediate formation in the serine proteases is proposed.

Effect of Structural Variation on Aldehyde Oxidase-Catalyzed Oxidation of Zoniporide

Current studies explored the effect of structural changes on the aldehyde oxidase (AO)-mediated metabolism of zoniporide (1). Zoniporide analogs with modifications of the acylguanidine moiety, the cyclopropyl group on the pyrazole ring, and the quinoline ring were studied for their AO-catalyzed metabolism using the human S9 fraction. Analysis of the half-lives suggested that subtle changes in the structure of 1 influenced its metabolism and that the guanidine and the quinoline moieties were prerequisites for AO-catalyzed oxidation to 2-oxozoniporide (M1). In contrast, replacement of the cyclopropyl group with other alkyl groups was tolerated. The effect of structural variation on AO properties was rationalized by docking 1 and its analogs into the human AO homology model. These studies indicated the importance of electrostatic, π-π stacking and hydrophobic interactions of the three motifs with residues in the active site. Differences in substrate properties were also rationalized by comparing their half-lives with cLogD, electrophilicity parameters [electrostatic potential (ESP) charges and energy of lowest unoccupied molecular orbitals (E(LUMO))], and the energies of formation of tetrahedral intermediates (J Med Chem 50:4642-4647, 2007). Whereas the success of energetics in predicting the AO substrate properties of analogs was 87%, the predictive ability of other descriptors was none (cLogD) to 60% (ESP charges and E(LUMO)). Overall, the structure-metabolism relationship could be rationalized using a combination of both the energy calculations and docking studies. This combination method can be incorporated into a strategy for mitigating AO liabilities observed in the lead candidate or studying structure-metabolism relationships of other AO substrates.

Prodrugs of Acyclovir – A Computational Approach

Density functional theory calculation results demonstrated that the efficiency of the acid-catalyzed hydrolysis of Kirby's acid amides 1-15 is strongly dependent on the substitution on the C-C double bond and the nature of the amide N-alkyl group. Further, the results established that while in the gas phase the hydrolysis rate-limiting step is the tetrahedral intermediate formation in polar solvents such as water, the rate-limiting step could be either the formation or the collapse of the tetrahedral intermediate depending on the substitution on the C-C double bond and on the amide nitrogen substituent. Based on a linear correlation between the calculated and experimental effective molarities, the study on the systems reported herein could provide a good basis for designing prodrug systems that are less hydrophilic than their parental drugs and can be used, in different dosage forms, to release the parent drug in a controlled manner. For example, based on the calculated log effective molarities values, the predicted t(1/2) (a time needed for 50% of the reactant to be hydrolyzed to products) for acyclovir prodrugs, ProD 1-4, was 29.2 h, 6097 days, 4.6 min, and 8.34 h, respectively. Hence, the rate by which acyclovir prodrug releases acyclovir can be determined according to the structural features of the linker (Kirby's acid amide moiety).

Capture hydrolysis signals in the microsomal stability assay: Molecular mechanisms of the alkyl ester drug and prodrug metabolism

The hydrolysis of carboxylic acid esters is often catalyzed by carboxylesterases in human liver microsomes, which is also a common 'noise' in the microsomal stability assay, a widely used screening protocol in drug discovery to monitor the activity of cytochrome P450 enzymes. Herein, we captured this 'noise', the hydrolysis signal of small alkyl ester drugs and prodrugs with a unique pairwise analysis of Pfizer's microsomal clearance database. The hydrolysis mechanisms were further elucidated with density functional theory and molecular docking approaches. The results suggested that the electronic properties of ester moieties, tetrahedral intermediate formation energies, and specific drug-enzyme molecular interactions are key factors for the determination of the metabolic fate of the studied alkyl esters, but individually these factors failed to correlate with the observed rate of hydrolysis.