Computational study of the aminolysis of esters. The reaction of methylformate with ammonia

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.

Precision Synthesis of Polysarcosine via Controlled Ring-Opening Polymerization of N-Carboxyanhydride: Fast Kinetics, Ultrahigh Molecular Weight, and Mechanistic Insights

The rapid and controlled synthesis of high-molecular-weight (HMW) polysarcosine (pSar), a potential polyethylene glycol (PEG) alternative, via the ring-opening polymerization (ROP) of N-carboxyanhydride (NCA) is rare and challenging. Here, we report the well-controlled ROP of sarcosine NCA (Sar-NCA) that is catalyzed by various carboxylic acids, which accelerate the polymerization rate up to 50 times, and enables the robust synthesis of pSar with an unprecedented ultrahigh molecular weight (UHMW) up to 586 kDa (DP ∼ 8200) and exceptionally narrow dispersity (D̵) below 1.05. Mechanistic experiments and density functional theory calculations together elucidate the role of carboxylic acid as a bifunctional catalyst that significantly facilitates proton transfer processes and avoids charge separation and suggest the ring opening of NCA, rather than decarboxylation, as the rate-determining step. UHMW pSar demonstrates improved thermal and mechanical properties over the low-molecular-weight counterparts. This work provides a simple yet highly efficient approach to UHMW pSar and generates a new fundamental understanding useful not only for the ROP of Sar-NCA but also for other NCAs.

An Exploration of the Hydrogen Bond Donor Ability of Ammonia

Ammonia is an important molecule due to its wide use in the fertiliser industry. It is also used in aminolysis reactions. Theoretical studies of the reaction mechanism predict that in reactive complexes and transition states, ammonia acts as a hydrogen bond donor forming N-H⋅⋅⋅O hydrogen bond. Experimental reports of N-H⋅⋅⋅O hydrogen bond, where ammonia acts as a hydrogen bond donor are scarce. Herein, the hydrogen bond donor ability of ammonia is investigated with three chalcogen atoms i. e. O, S, and Se using matrix isolation infrared spectroscopy and electronic structure calculations. In addition, the chalcogen bond acceptor ability of ammonia has also been investigated. The hydrogen bond acceptor molecules used here are O(CH3 )2 , S(CH3 )2 , and Se(CH3 )2 . The formation of the 1 : 1 complex has been monitored in the N-H symmetric and anti-symmetric stretching modes of ammonia. The nature of the complex has been delineated using Atoms in Molecules analysis, Natural Bond Orbital analysis, and Energy Decomposition Analysis. This work presents the first comparison of the hydrogen bond donor ability of ammonia with O, S, and Se.

Acylative kinetic resolution of racemic methyl-substituted cyclic alkylamines with 2,5-dioxopyrrolidin-1-yl (R)-2-phenoxypropanoate

The diastereoselective acylation of a number of racemic methyl-substituted cyclic alkylamines with active esters of 2-phenoxypropanoic acid was studied in detail. The ester of (R)-2-phenoxypropanoic acid and N-hydroxysuccinimide was found to be the most selective agent. The highest stereoselectivity was observed in the kinetic resolution of racemic 2-methylpiperidine in toluene at -40 °C (selectivity factor s = 73) with the predominant formation of (R,R)-amide (93.7% de). To explain the observed stereoselectivity, DFT modelling of the transition states in the reactions of the title acylating agent with 2-methylpiperidine and 2-methylpyrrolidine was performed. The calculated values were in good agreement with experimental data. It has been demonstrated that the acylation proceeds via a concerted mechanism, in which the addition of an amine occurs simultaneously with the elimination of the hydroxysuccinimide fragment. The high stereoselectivity of the (R,R)-amide formation is largely ensured by the lower steric hindrances in the transition states as compared to the formation of (R,S)-amide.

Mechanistic insights into lysine-targeting covalent inhibition through a theoretical study of ester aminolysis

Development of targeted covalent inhibitors in drug design has a broad and important interest and many efforts are currently being made in this direction. Targeted covalent inhibitors have special relevance in oncology due to the possibilities they offer to overcome the problems of acquired resistance. In recent experiments, lysine-targeting has been envisaged for the irreversible inhibition of the heterodimeric lipid kinase phosphoinositide 3-kinase delta (PI3Kδ). Activated esters have been evaluated and shown to be promising inhibitors of this enzyme, but the reaction mechanisms display specificities that are not yet fully understood. In the present work, we have carried out a theoretical study of the aminolysis reaction of model esters in aqueous solution to gain insights into the corresponding biological processes. We have found that phenolic esters bearing electron-withdrawing groups are particularly reactive. The predicted mechanism involves the formation of a tetrahedral zwitterionic intermediate, which dissociates into an alkoxide and a protonated amide, this charge separation being the driving force for the subsequent proton transfer and final product formation. Structure-reactivity relationships are reported and shown to be a useful tool for evaluating potential inhibitor candidates.

Theoretical Study of the Mechanism of Ribosomal Peptide Bond Formation Using the ONIOM Method

Peptide bond formation in living cells occurs at the peptidyl transferase center (PTC) of the large ribosomal subunit and involves the transfer of the peptidyl group from peptidyl-tRNA to aminoacyl-tRNA. Despite numerous kinetic and theoretical studies, many details of this reaction -such as whether it proceeds via a stepwise or concerted mechanism- remain unclear. In this study, we calculated the geometry and energy of the transition states and intermediates in peptide bond formation in the PTC environment using the ONIOM (our own n-layered integrated molecular orbital and molecular mechanics) method. The calculations indicated that the energy of the transition states of stepwise mechanisms are lower than those of concerted mechanisms and suggested that the reaction involves a neutral tetrahedral intermediate that is stabilized through the hydrogen-bonding network in the PTC environment. The results will lead to a better understanding of the mechanism of peptidyl transfer reaction, and resolve fundamental questions of the steps and molecular intermediates involved in peptide bond formation in the ribosome.

Strain-Induced Nucleophilic Ring Opening of Donor-Acceptor Cyclopropenes for Synthesis of Monosubstituted Succinic Acid Derivatives

1,2,3-Trisubstituted donor-acceptor cyclopropenes (DACPs) generated in situ from enoldiazo compounds react with nucleophiles to form α-substituted succinic acid derivatives in high yields. Initial dirhodium(II) carboxylate catalysis rapidly converts enoldiazo-acetates or -acetamides to DACPs that undergo catalyst-free Favorskii ring opening with amines, and also with anilines, alcohols, and thiols, when facilitated by catalytic amounts of 4-dimethylaminopyridine (DMAP). This methodology provides easy access to mixed esters and amides of monosubstituted succinic acids, including derivatives of naturally occurring compounds. It also affords dihydrazide, dihydroxamic acid, and diamide derivatives, as well as α-substituted tetrahydropyridazine-3,6-diones in high yields. Attempts to generate optically enriched DACPs were not successful because their populations having the R and S configurations formed with a chiral dirhodium catalyst are quite similar, and the loss of enantiocontrol likely originates from the DACP ring forming step which is reversible with its intermediate metal carbene.

Nucleofugality hierarchy, in the aminolysis reaction of 4-cyanophenyl 4-nitrophenyl carbonate and thionocarbonate. Experimental and theoretical study

Brønsted plot of the reaction of 9 with secondary alicyclic amines. The concave upward non-linear plot is in accordance with two parallel mechanistic pathways.

Experimental and theoretical study of the role of CH/π interactions in the aminolysis reaction of acetyl galactoside

Molecular recognition of saccharides is a growing field, which has many implications in cancer therapy, drug discovery, and cellular communication among others. The participation of CH/π interactions in this event is well known. Nevertheless, the intrinsic role of CH/π for modulating chemical reactions is still far from being applicable. In this experimental and computational work we have evaluated the participation of CH/π interactions in the aminolysis reaction of acetyl galactoside promoted with different 6-substituted 2(1H)-pyridones. Two features have been incorporated to the promoter molecular structure, on one end the promoting pyridone group and on the other end the recognition moiety, joined together by an alkyne spacer. The small increment in the observed pseudo-first-order rate constant values (k_obs) was related to the stability of the transition state provided by noncovalent interactions, including CH/π interactions. A longer alkyne spacer was necessary to improve the molecular recognition of the galactoside substrate. The trend of the calculated activation energy values (ΔE_RTS) was in good accordance with the experimental rate constant values.

Exploiting Ring-Opening Aminolysis-Condensation as a Polymerization Pathway to Structurally Diverse Biobased Polyamides

A pathway to biobased polyamides (PAs) via ring-opening aminolysis-condensation (ROAC) under benign conditions with diverse structure was designed. Ethylene brassylate (EB), a plant oil-derived cyclic dilactone, was used in combination with an array of diamines of diverse chemical structure, and ring-opening of the cyclic dilactone EB was revealed as a driving force for the reaction. The ROAC reactions were adjusted, and reaction conditions of 100 °C under atmospheric pressure using 1,5,7-triazabicyclo[4.4.0]dec-5-ene (TBD) as a catalyst for 24 h were optimal. The structures of the polyamides were confirmed by mass spectroscopy, FTIR, and NMR, and the PAs had viscosity average molecular weights ( M_η) of ∼5-8 kDa. Glassy or semicrystalline PAs with glass transition temperatures between 48 and 55 °C, melting temperatures of 120-200 °C for the semicrystalline PAs, and thermal stabilities above 400 °C were obtained and were comparable to the existing PAs with similar structures. As a proof-of-concept of their usage, one of the PAs was shown to form fibers by electrospinning and films by melt pressing. Compared to conventional methods for PA synthesis, the ROAC route portrayed a reaction temperature at least 60-80 °C lower, could be readily carried out without a low-pressure environment, and eliminated the use of solvents and toxic chemicals. Together with the plant oil-derived monomer (EB), the ROAC route provided a sustainable alternative to design biobased PAs.

Mechanistic investigation of the uncatalyzed esterification reaction of acetic acid and acid halides with methanol: a DFT study

Implementation of catalysts to drive reactions from reactants to products remains a burden to synthetic and organic chemists. In spite of investigations into the kinetics and mechanism of catalyzed esterification reactions, less effort has been made to explore the possibility of an uncatalyzed esterification process. Therefore, a comprehensive mechanistic perspective for the uncatalyzed mechanism at the molecular level is presented. Herein, we describe the non-catalyzed esterification reaction of acetic acid and its halide derivatives (XAc, where X= OH, F, Cl, Br, I) with methanol (MeOH) through a concerted process. The reaction in vacuum and methanol was performed using the density functional theory (DFT) method at M06-2X level with def2-TZVP basis set after a careful literature survey and computations. Esterification through cyclic 4- or 6-membered transition state structures in one- or two-step concerted mechanisms were investigated. The present study outlines the possible cyclic geometry conformations that may occur during experiments at simple ratio of reactants. The free energy of activation for acetic acid and acetyl chloride are 36 kcal mol(-1) and 21 kcal mol(-1), respectively. These are in good agreement with available experimental results from the literature. The selected quantum chemical descriptors proved to be useful tools in chemical reactivity prediction for the reaction mechanism. This quantum mechanics study can serve as a necessary step towards revisiting uncatalyzed reaction mechanisms in some classical organic reactions.

Theoretical estimation of kinetic parameters for nucleophilic substitution reactions in solution: an application of a solution translational entropy model

The kinetic parameters, such as activation entropy, activation enthalpy, activation free-energy, and reaction rate constant, for a series of nucleophilic substitution (SN) reactions in solution, are investigated using both a solution-phase translational entropy model and an ideal gas-phase translational entropy model. The results obtained from the solution translational entropy model are in excellent agreement with the experimental values, while the overestimation of activation free-energy from the ideal gas-phase translational entropy model is as large as 6.9 kcal mol(-1). For some of the reactions studied, such as and in methanol, and and in aqueous solution, the explicit + implicit model, namely, a cluster-continuum type model, should be employed to account for the strong solvent-solute interactions. In addition, the explicit + implicit models have also been applied to the DMSO-H2O mixtures, which would open up a door to investigate the reactions in a mixed solvent using density functional theory (DFT) methods.

N-heterocyclic carbene copper(I) catalysed N-methylation of amines using CO2

The N-methylation of amines using CO2 and PhSiH3 as source of CH3 was efficiently performed using a N-heterocyclic carbene copper(i) complex. The methodology was found compatible with aromatic and aliphatic primary and secondary amines. Synthetic and computational studies have been carried out to support the proposed reaction mechanism for this transformation.

Concerted amidation of activated esters: reaction path and origins of selectivity in the kinetic resolution of cyclic amines via N-heterocyclic carbenes and hydroxamic acid cocatalyzed acyl transfer

The N-heterocyclic carbene and hydroxamic acid cocatalyzed kinetic resolution of cyclic amines generates enantioenriched amines and amides with selectivity factors up to 127. In this report, a quantum mechanical study of the reaction mechanism indicates that the selectivity-determining aminolysis step occurs via a novel concerted pathway in which the hydroxamic acid plays a key role in directing proton transfer from the incoming amine. This modality was found to be general in amide bond formation from a number of activated esters including those generated from HOBt and HOAt, reagents that are broadly used in peptide coupling. For the kinetic resolution, the proposed model accurately predicts the faster reacting enantiomer. A breakdown of the steric and electronic control elements shows that a gearing effect in the transition state is responsible for the observed selectivity.

Microsolvated transition state models for improved insight into chemical properties and reaction mechanisms

Over the years, several methods have been developed to effectively represent the chemical behavior of solutes in solvents. The environmental effects arising due to solvation can generally be achieved either through inclusion of discrete solvent molecules or by inscribing into a cavity in a homogeneous and continuum dielectric medium. In both these approaches of computational origin, the perturbations on the solute induced by the surrounding solvent are at the focus of the problem. While the rigor and method of inclusion of solvent effects vary, such solvation models have found widespread applications, as evident from modern chemical literature. A hybrid method, commonly referred to as cluster-continuum model (CCM), brings together the key advantages of discrete and continuum models. In this perspective, we intend to highlight the latent potential of CCM toward obtaining accurate estimates on a number of properties as well as reactions of contemporary significance. The objective has generally been achieved by choosing illustrative examples from the literature, besides expending efforts to bring out the complementary advantages of CCM as compared to continuum or discrete solvation models. The majority of examples emanate from the prevalent applications of CCM to organic reactions, although a handful of interesting organometallic reactions have also been discussed. In addition, increasingly accurate computations of properties like pK(a) and solvation of ions obtained using the CCM protocol are also presented.

Mechanism of the aminolysis of Fischer alkoxy and thiocarbene complexes: a DFT study

B3LYP calculations have been carried out to study the reaction mechanism of the aminolysis of Fischer carbene complexes of the type (CO)(5)Cr=C(XMe)R (X = O and S; R = Me and Ph). We have explored different possible reaction mechanisms either through neutral or zwitterionic intermediates as well as a general base catalysis assisted by an ammonia molecule. Our results show that the most favorable pathway for the aminolysis of Fischer carbene complexes is through a stepwise reaction via a zwitterionic intermediate generated by the initial nucleophilic attack. We have found that the ammonia-catalyzed mechanism entails a significantly lower barrier for the rate-determining step than the uncatalyzed one. At lower pressure gas-phase conditions, the rate-determining step corresponds to the concerted proton transfer and MeXH elimination. Thiocarbene complexes show a higher energy barrier for this rate-determining step due to the lower basicity of the MeS(-) substituent. At higher pressure or in solution, the rate-determining step corresponds to the initial nucleophilic attack. Our results indicate that the transition state of the nucleophilic attack is more advanced and has a higher barrier for alkoxycarbene than thiocarbene complexes due to the stronger pi-donor character of the alkoxy group that reduces the electrophilicity of the attacked carbene atom making the nucleophilic attack more difficult.