Modeling the reaction mechanisms of the amide hydrolysis in an N-(o-carboxybenzoyl)-L-amino acid

Competing Reaction Mechanisms of Peptide Bond Formation in Water Revealed by Deep Potential Molecular Dynamics and Path Sampling

The formation of an amide bond is an essential step in the synthesis of materials and drugs, and in the assembly of amino acids to form peptides. The mechanism of this reaction has been studied extensively, in particular to understand how it can be catalyzed, but a representation capable of explaining all the experimental data is still lacking. Numerical simulation should provide the necessary molecular description, but the solvent involvement poses a number of challenges. Here, we combine the efficiency and accuracy of neural network potential-based reactive molecular dynamics with the extensive and unbiased exploration of reaction pathways provided by transition path sampling. Using microsecond-scale simulations at the density functional theory level, we show that this method reveals the presence of two competing distinct mechanisms for peptide bond formation between alanine esters in aqueous solution. We describe how both reaction pathways, via a general base catalysis mechanism and via direct cleavage of the tetrahedral intermediate respectively, change with pH. This result contrasts with the conventional mechanism involving a single pathway in which only the barrier heights are affected by pH. We show that this new proposal involving two competing mechanisms is consistent with the experimental data, and we discuss the implications for peptide bond formation under prebiotic conditions and in the ribosome. Our work shows that integrating deep potential molecular dynamics with path sampling provides a powerful approach for exploring complex chemical mechanisms.

Investigating Betrixaban Maleate drug degradation profiles, isolation and characterization of unknown degradation products by mass-triggered preparative HPLC, HRMS, and NMR

Betrixaban Maleate, a novel oral, once-daily factor Xa inhibitor drug substance, was subjected to stress testing under a wide range of degradation conditions, including acidic hydrolysis, alkaline hydrolysis, oxidative, thermal, and photolytic, to determine its inherent stability. The drug was biodegradable in acidic and alkaline environments, and three new degradation products were identified. Two degraded products are formed in an acidic environment, while the third is in alkaline conditions. The three degradants were identified using UPLC-ESI/MS and isolated using mass-triggered preparative HPLC, and their structures were unambiguously elucidated using HRMS and 2D NMR techniques. Based on spectral and chromatographic data, it was firmly proven that these distinct degradation products were the betrixaban chemical's hydrolysis components. The formation of the degradants has been hypothesized through several possible mechanisms.

Computed ammonia affinity for evaluating Lewis acidity of organoboronates and organoboronamides

Lewis acidity of organoboronates [B(pin), B(neop), B(cat), B(eg), B(nad)] and organoboronamides [B(dan), B(aam), B(mdan)] has been found to be unifiedly evaluated by computed ammonia affinity (AA), while other methods [LUMO energies, global electrophilicity index (GEI), fluoride ion affinity (FIA)] were only partially applicable. The relationships between the AA values and such structural characters including the B-X bond lengths, the X-B-X angles, and the changes in the B-X bond lengths in the formation of the ammonia adducts were also described.

A computational study of the reaction mechanism involved in the fast cleavage of an unconstrained amide bond assisted by an amine intramolecular nucleophilic attack

In the present work, the fast amide bond cleavage of [3-((1R,5S,7s)-3-azabicyclo[3.3.1]nonane-7-carbonyl)-3-azabicyclo[3.3.1]nonane-7-carboxylic acid (bi-ATDO)], through an intramolecular nucleophilic attack of an amine group is evaluated. First, six possible peptide bond cleavage mechanisms, two of them including a water molecule, are described at the ωB97XD/6-311 + G(d,p)//MP2/6-311 + G(d,p) level of theory. The reaction consisting of an intramolecular nitrogen nucleophilic attack followed by a proton transfer and the amide bond cleavage is determined as the most favorable mechanism. The activation free energy computed for the latter is 20.5 kcal mol^-1 , which agrees with the reported experimental result of 24.8 kcal mol^-1 . Inclusion of a water molecule to assist the first step of the reaction results in an activation free energy increase of about 17 kcal mol^-1 . All the steps in the most favorable mechanism are studied more in detail employing intrinsic reaction coordinate as well as the reaction force and reaction electronic flux analysis.

Mechanisms of peptide hydrolysis by aspartyl and metalloproteases

Peptide hydrolysis has been involved in a wide range of biological, biotechnological, and industrial applications.

Modeling the chelation of As(III) in lewisite by dithiols using density functional theory and solvent-assisted proton exchange

Dithiols such as British anti-lewisite (BAL, rac-2,3-dimercaptopropanol) are an important class of antidotes for the blister agent lewisite (trans-2-chlorovinyldichloroarsine) and, more generally, are chelating agents for arsenic and other toxic metals. The reaction of the vicinal thiols of BAL with lewisite through the chelation of the As(III) center has been modeled using density functional theory (DFT) and solvent-assisted proton exchange (SAPE), a microsolvation method that uses a network of water molecules to mimic the role of bulk solvent in models of aqueous phase chemical reactions. The small activation barriers for the stepwise SN2-type nucleophilic attack of BAL on lewisite (0.7-4.9kcal/mol) are consistent with the favorable leaving group properties of the chloride and the affinity of As(III) for soft sulfur nucleophiles. Small, but insignificant, differences in activation barriers were found for the initial attack of the primary versus secondary thiol of BAL and the R vs S enantiomer. An examination of the relative stability of various dithiol-lewisite complexes shows that ethanedithiols like BAL form the most favorable chelation complexes because the angles formed in five-membered ring are most consistent with the hybridization of As(III). More obtuse S-As-S angles are required for larger chelate rings, but internal As⋯N or As⋯O interactions can enhance the stability of moderate-sized rings. The low barriers for lewisite detoxification by BAL and the greater stability of the chelation complexes of small dithiols are consistent with the rapid reversal of toxicity demonstrated in previously reported animal models.

On the possible biological relevance of HSNO isomers: a computational investigation

Thionitrous acid (HSNO), the smallest S-nitrosothiol, has been identified as a potential biologically active molecule that connects the biochemistries of two important gasotransmitters, nitric oxide (NO) and hydrogen sulfide (H2S). Here, we computationally explore possible isomerization reactions of HSNO that may occur under physiological conditions using high-level coupled-cluster as well as density functional theory and composite CBS-QB3 methodology calculations. Gas-phase calculations show that the formation of the tautomeric form HONS and the Y-isomer SN(H)O is thermodynamically feasible, as they are energetically close, within ∼6 kcal mol(-1), to HSNO, while the recently proposed three-membered ring isomer is not thermodynamically or kinetically accessible. The gas-phase intramolecular proton-transfer reactions required for HSNO isomerization into HONS and SN(H)O are predicted to have prohibitively high reaction barriers, 30-50 kcal mol(-1). However, the polar aqueous environment and water-assisted proton shuttle should decrease these barriers to ∼9 kcal mol(-1), which makes these two isomers kinetically accessible under physiological conditions. Our calculations also support the possibility of an aqueous reaction between the Y-isomer SN(H)O and H2S leading to biologically active nitroxyl HNO. These results suggest that the formation of HSNO in biological milieu can lead to various derivative species with their own, possibly biologically relevant, activity.

Experimental and theoretical studies on the rearrangement of 2-oxoazepane α,α-amino acids into 2'-oxopiperidine β(2,3,3) -amino acids: an example of intramolecular catalysis

Enantiopure β-amino acids represent interesting scaffolds for peptidomimetics, foldamers and bioactive compounds. However, the synthesis of highly substituted analogues is still a major challenge. Herein, we describe the spontaneous rearrangement of 4-carboxy-2-oxoazepane α,α-amino acids to lead to 2'-oxopiperidine-containing β(2,3,3) -amino acids, upon basic or acid hydrolysis of the 2-oxoazepane α,α-amino acid ester. Under acidic conditions, a totally stereoselective synthetic route has been developed. The reordering process involved the spontaneous breakdown of an amide bond, which typically requires strong conditions, and the formation of a new bond leading to the six-membered heterocycle. A quantum mechanical study was carried out to obtain insight into the remarkable ease of this rearrangement, which occurs at room temperature, either in solution or upon storage of the 4-carboxylic acid substituted 2-oxoazepane derivatives. This theoretical study suggests that the rearrangement process occurs through a concerted mechanism, in which the energy of the transition states can be lowered by the participation of a catalytic water molecule. Interestingly, it also suggested a role for the carboxylic acid at position 4 of the 2-oxoazepane ring, which facilitates this rearrangement, participating directly in the intramolecular catalysis.

Elucidating the catalytic mechanism of β-secretase (BACE1): a quantum mechanics/molecular mechanics (QM/MM) approach

In this quantum mechanics/molecular mechanics (QM/MM) study, the mechanisms of the hydrolytic cleavage of the Met2-Asp3 and Leu2-Asp3 peptide bonds of the amyloid precursor protein (WT-substrate) and its Swedish mutant (SW) respectively catalyzed by β-secretase (BACE1) have been investigated by explicitly including the electrostatic and steric effects of the protein environment in the calculations. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely acknowledged as a promising therapeutic target. The general acid-base mechanism followed by the enzyme proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. The formation of the gem-diol intermediate occurs with the barriers of 19.6 and 16.1 kcal/mol for the WT- and SW-substrate respectively. The QM/MM energetics predict that with the barriers of 21.9 and 17.2 kcal/mol for the WT- and SW-substrate respectively the cleavage of the peptide bond occurs in the rate-determining step. The computed barriers are in excellent agreement with the measured barrier of ∼18.0 kcal/mol for the SW-substrate and in line with the experimental observation that the cleavage of this substrate is sixty times more efficient than the WT-substrate.

SOLVENT AND SUBSTITUENT EFFECTS ON THE INTRAMOLECULAR AMIDE HYDROLYSIS OF N-METHYLMALEAMIC ACID

Intramolecular amide hydrolysis of N-methylmaleamic acid is revisited at the B3LYP/6-311G(2df, p)//B3LYP/6-31G(d, p) + ZVPE level, including solvent effects at the CPCM-B3LYP/6-311G(2df, p)//Onsager-B3LYP/6-31G(d, p) + ZPVE level. The concerted reaction mechanism is energetically favorable over stepwise reaction mechanisms in both the gas phase and solution. The calculated reaction barriers are significantly lower in solution than in the gas phase. In addition, it is concluded that the substituents of the four N-methylmaleamic acid derivatives considered herein have a significant effect on the gas-phase reaction barriers but a smaller, or little, effect on the barriers in solution.

Acid-catalyzed reactions of twisted amides in water solution: competition between hydration and hydrolysis

The acid-catalyzed reactions of twisted amides in water solution were investigated by using cluster-continuum model calculations. In contrast to the previous widely suggested concerted hydration of the C=O group, our calculations show that the reaction proceeds in a practically stepwise manner, and that the hydration and hydrolysis channels of the C-N bond compete. The Eigen ion (H(3)O(+)) is the key species involved in the reaction, and it modulates the hydration and hydrolysis reaction pathways. The phenyl substitution in the twisted amide not only activates the N-CO bond, but also stabilizes the hydrolysis product through n(N)→π(phenyl) delocalization, leading exclusively to the hydrolysis product of the ring-opened carboxylic acid. Generally, the twisted amides are more active than the planar amides, and such a rate acceleration results mainly from the increase in exothermicity in the first N-protonation step; the second step of the nucleophilic attack is less affected by the twisting of the amide bond. The present results show good agreement with the available experimental observations.

Computational modeling of substrate specificity and catalysis of the β-secretase (BACE1) enzyme

In this combined MD simulation and DFT study, interactions of the wild-type (WT) amyloid precursor protein (APP) and its Swedish variant (SW), Lys670 → Asn and Met671 → Leu, with the beta-secretase (BACE1) enzyme and their cleavage mechanisms have been investigated. BACE1 catalyzes the rate-limiting step in the generation of 40-42 amino acid long Alzheimer amyloid beta (Aβ) peptides. All key structural parameters such as position of the flap, volume of the active site, electrostatic binding energy, structures, and positions of the inserts A, D, and F and 10s loop obtained from the MD simulations show that, in comparison to the WT-substrate, BACE1 exhibits greater affinity for the SW-substrate and orients it in a more reactive conformation. The enzyme-substrate models derived from the MD simulations were further utilized to investigate the general acid/base mechanism used by BACE1 to hydrolytically cleave these substrates. This mechanism proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. For the WT-substrate, the overall barrier of 22.4 kcal/mol for formation of the gem-diol intermediate is 3.3 kcal/mol higher than for the SW-substrate (19.1 kcal/mol). This process is found to be the rate-limiting in the entire mechanism. The computed barrier is in agreement with the measured barrier of ca. 18.00 kcal/mol for the WT-substrate and supports the experimental observation that the cleavage of the SW-substrate is 60 times more efficient than the WT-substrate.

Effects of Mixed H2O—CH3CN and H2O—HCONMe2 Solvents on the Rate of Cleavage of N-Benzylphthalamic Acid

The values of pseudo first-order rate constants, k1, for the O-cyclization and k3 for N-cyclization of N-benzylphthalamic acid 1, obtained at 0.02 M HCl and 35°C, show a respective nonlinear increase from 17.6 × 10−6 to 29.7 × 10−6 s−1 and decrease from 12.5 × 10−7 to 3.75 × 10−7 s−1 with increase in the CH3CN content from 2 to 60% v/v in mixed aqueous solvent. Increase in the CH3CN content from 60 to 90% v/v decreases k1 from 29.7 × 10−6 to 24.4 × 10−6 s−1 and increases k3 from 3.75 × 10−7 to 14.7 × 10−7 s−1. The respective values of k1 and k3 decrease nonlinearly from 16.0 × 10−6 to 3.0 × 10−6 s−1 and 14.0 × 10−7 to ~4.0 × 10−7 s−1 with increase in the content of N,N-dimethylformamide (HCONMe2) from 2 to 80% v/v in mixed aqueous solvents. These results indicate rather mild rate-retarding effects of mixed H2O—CH3CN and H2O—HCONMe2 solvents on the rates of the intramolecular reactions.

Modeling Competitive Reaction Mechanisms of Peroxynitrite Oxidation of Guanine

5-Guanidino-4-nitroimidazole is a stable product from the peroxynitrite induced one-electron oxidation of guanine. Reaction mechanisms to form the 5-guanidino-4-nitroimidazole as well as 8-nitroguanine, through the combination of the guanine radical cation and nitrogen dioxide radical and through the combination of the deprotonated neutral guanine radical and nitrogen dioxide radical, have been investigated by the use of the B3LYP method of density functional theory. Our calculations suggest that the guanine radical cation mechanism is preferred over the neutral guanine radical mechanism and that a water molecule is involved in the reaction as a catalyst or as a reactant.

A self-stabilized model of the chymotrypsin catalytic pocket. The energy profile of the overall catalytic cycle

A model of the catalytic triad of chymotrypsin is built assuring the arrangement and properties as they are within the complete enzyme. The model contains 18 amino acid residues of chymotrypsin and its substrate. A total of 135 atoms (including 70 heavy atoms) were subjected to full ab initio geometry optimizations through 127 individual steps along the reaction coordinate of the complete catalytic mechanism. It was shown that the described model of the catalytic apparatus forms a self-stabilized molecule ensemble without the rest of the enzyme and substrate. According to the calculations, the formations of the first and second tetrahedral intermediates in the model have 20.3 and 15.7 kcal/mol activation energy barriers, respectively. Removing elements of the catalytic apparatus such as the (1) catalytic aspartate or (2) the anion hole, as well as (3) inserting a water molecule "early" in the catalytic process, or (4) introducing conformational rigidity of the substrate, results in an increase of the above energy barrier of the first catalytic step in the model by 6.4, 13.7, 3.7, and 4.1 kcal/mol, respectively. Based on the calculated process one can conclude that the catalytic reaction in this model is much more similar to the reaction in the enzyme than to the reference reaction. To our knowledge, this is the first model system that mimics the complete catalytic mechanism.

The Mechanism of Formamide Hydrolysis in Water from Ab Initio Calculations and Simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).

Influence of the number of water molecules on the mechanism of N-sulfinylaniline hydrolysis

The mechanism of the uncatalyzed hydrolysis of N-sulfinylaniline (Ph-N=S=O) has been studied with B3LYP/6-31+G(2d,2p) in the gas phase, with explicit treatment of water molecules. Hydrolysis involves water attack on sulfur, with a close to perpendicular alignment of a water molecule and the NSO plane in both prereaction complexes and transition states for the rate-determining step. Consequently, the distance of the weak S···O interaction, together with the efficiency of protonation of either nitrogen (attack across the N=S bond) or oxygen (attack across the S=O bond) atoms of the NSO group, determines the height of the activation barrier for hydrolysis. While the reaction with one water molecule is characterized by an unreasonably high enthalpy of activation, a cooperative effect from the weak interactions appears with the inclusion of a second water molecule, where both participate in the reaction, and the activation enthalpy drops significantly. The preference for attack across the S=O bond that is found in the reaction with one water molecule switches to a dominance of attack across the N=S bond in the reaction with three water molecules.Key words: N-sulfinylaniline, hydrolysis, mechanism, density functional theory (DFT).

Theoretical Studies on the Mechanism of Acid-Promoted Hydrolysis of N-Formylaziridine in Comparison with Formamide

[reaction: see text] We present an ab initio study of the acid-promoted hydrolysis reaction mechanism of N-formylaziridine in comparison with formamide. Since the rate of amide hydrolysis reactions depends on the formation of the tetrahedral intermediate, we focused our attention mainly on the reactant complex, the tetrahedral intermediate, and the transition state connecting these two stationary points. Geometries were optimized using the density functional theory, and the energetics were refined using ab initio theory including electron correlation. Solvent effects were investigated by using polarizable continuum method calculations. The proton-transfer reaction between the O-protonated and N-protonated amides was investigated. In acidic media, despite that the N-protonated species is more stable than the O-protonated one, it is predicted that both N-protonated and O-protonated pathways compete in the hydrolysis reaction of N-formylaziridine.