Ab initio QM/MM dynamics simulation of the tetrahedral intermediate of serine proteases: insights into the active site hydrogen-bonding network

Ab initio QM/MM dynamics simulation is employed to examine the stability of the tetrahedral intermediate during the deacylation step in elastase-catalyzed hydrolysis of a simple peptide. An extended quantum region includes the catalytic triad, the tetrahedral structure, and the oxyanion hole. The calculations indicate that the tetrahedral intermediate of serine proteases is a stable species on the picosecond time scale. On the basis of geometrical and dynamical properties, and in agreement with many experimental and theoretical studies, it is suggested that the crucial hydrogen bonds involved in stabilizing this intermediate are between Asp-102 and His-57 and between the charged oxygen of the intermediate and the backbone N-H group of Gly-193 in the oxyanion hole. The mobility of the imidazolium ring between O(w) and O(gamma), two of the oxygens of the tetrahedral structure, shows how the intermediate could proceed toward the product state without a "ring-flip mechanism", proposed earlier on the basis of NMR data. In addition to the proposed C(epsilon)(1)-H.O hydrogen bond between the imidazolium ring and the backbone carbonyl of Ser-214, we observe an alternative C(epsilon)(1)-H.O hydrogen bond with the backbone carbonyl of Thr-213, that can stabilize the intermediate during the imidazolium movement. Proton hopping occurs between Asp-102 and His-57 during the simulation. The proton is, however, largely localized on the nitrogen, and hence it does not participate in a low-barrier hydrogen bond. The study also suggests factors that may be implicated in product release: breaking the hydrogen bond of the charged oxygen with the backbone of Ser-195 in the oxyanion hole and a loop opening between residues 216-225 that enables the breaking of a hydrogen bond in subsite S(3).     

Ab initio and density functional theory studies of the catalytic mechanism for ester hydrolysis in serine hydrolases

We present results from ab initio and density functional theory studies of the mechanism for serine hydrolase catalyzed ester hydrolysis. A model system containing both the catalytic triad and the oxyanion hole was studied. The catalytic triad was represented by formate anion, imidazole, and methanol. The oxyanion hole was represented by two water molecules. Methyl formate was used as the substrate. In the acylation step, our computations show that the cooperation of the Asp group and oxyanion hydrogen bonds is capable of lowering the activation barrier by about 15 kcal/mol. The transition state leading to the first tetrahedral intermediate in the acylation step is rate limiting with an activation barrier (ΔE₀) of 13.4 kcal/mol. The activation barrier in the deacylation step is smaller. The double-proton-transfer mechanism is energetically unfavorable by about 2 kcal/mol. The bonds between the Asp group and the His group, and the hydrogen bonds in the oxyanion hole, increase in strength going from the Michaelis complex toward the transition state and the tetrahedral intermediate. In the acylation step, the tetrahedral intermediate is a very shallow minimum on the energy surface and is not viable when molecular vibrations are included. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 89–103, 1998     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

Analysis of catalytic mechanism of serine proteases. Viability of the ring-flip hypothesis

Quantum calculations are applied to the active site of serine proteases, including four specific residues and a water molecule, as well as a substrate and proton donors in the oxyanion hole. Residues are tethered to the protein backbone of an X-ray structure but otherwise allowed to move freely to their lowest energy positions. The viability of the ring-flip hypothesis, which proposes that a 180 degrees rotation of the His-57 imidazole ring facilitates the catalysis, is assessed by comparison of energies of configurations both before and after such a flip. Specifically considered is the contribution to catalysis of the Ser-214 residue and a water molecule that is observed in the active site. The calculations provide detailed information concerning the nature, geometry, and strength of hydrogen bonds that are formed within the active site at each stage of the enzymatic mechanism.     

Theoretical studies on the deacylation step of serine protease catalysis in the gas phase, in solution, and in elastase

The deacylation step of serine protease catalysis is studied using DFT and ab initio QM/MM calculations combined with MD/umbrella sampling calculations. Free energies of the entire reaction are calculated in the gas phase, in a continuum solvent, and in the enzyme elastase. The calculations show that a concerted mechanism in the gas phase is replaced by a stepwise mechanism when solvent effects or an acetate ion are added to the reference system, with the tetrahedral intermediate being a shallow minimum on the free energy surface. In the enzyme, the tetrahedral intermediate is a relatively stable species ( approximately 7 kcal/mol lower in energy than the transition state), mainly due to the electrostatic effects of the oxyanion hole and Asp102. It is formed in the first step of the reaction, as a result of a proton transfer from the nucleophilic water to His57 and of an attack of the remaining hydroxyl on the ester carbonyl. This is the rate-determining step of the reaction, which requires approximately 22 kcal/mol for activation, approximately 5 kcal/mol less than the reference reaction in water. In the second stage of the reaction, only small energy barriers are detected to facilitate the proton transfer from His57 to Ser195 and the breakdown of the tetrahedral intermediate. Those are attributed mainly to a movement of Ser195 and to a rotation of the His57 side chain. During the rotation, the imidazolium ion is stabilized by a strong H-bond with Asp102, and the C(epsilon)(1)-H...O H-bond with Ser214 is replaced by one with Thr213, suggesting that a "ring-flip mechanism" is not necessary as a driving force for the reaction. The movements of His57 and Ser195 are highly correlated with rearrangements of the binding site, suggesting that product release may be implicated in the deacylation process.     

Investigation of the mechanism of the cell wall DD-carboxypeptidase reaction of penicillin-binding protein 5 of Escherichia coli by quantum mechanics/molecular mechanics calculations

Penicillin-binding protein 5 (PBP 5) of Escherichia coli hydrolyzes the terminal D-Ala-D-Ala peptide bond of the stem peptides of the cell wall peptidoglycan. The mechanism of PBP 5 catalysis of amide bond hydrolysis is initial acylation of an active site serine by the peptide substrate, followed by hydrolytic deacylation of this acyl-enzyme intermediate to complete the turnover. The microscopic events of both the acylation and deacylation half-reactions have not been studied. This absence is addressed here by the use of explicit-solvent molecular dynamics simulations and ONIOM quantum mechanics/molecular mechanics (QM/MM) calculations. The potential-energy surface for the acylation reaction, based on MP2/6-31+G(d) calculations, reveals that Lys47 acts as the general base for proton abstraction from Ser44 in the serine acylation step. A discrete potential-energy minimum for the tetrahedral species is not found. The absence of such a minimum implies a conformational change in the transition state, concomitant with serine addition to the amide carbonyl, so as to enable the nitrogen atom of the scissile bond to accept the proton that is necessary for progression to the acyl-enzyme intermediate. Molecular dynamics simulations indicate that transiently protonated Lys47 is the proton donor in tetrahedral intermediate collapse to the acyl-enzyme species. Two pathways for this proton transfer are observed. One is the direct migration of a proton from Lys47. The second pathway is proton transfer via an intermediary water molecule. Although the energy barriers for the two pathways are similar, more conformers sample the latter pathway. The same water molecule that mediates the Lys47 proton transfer to the nitrogen of the departing D-Ala is well positioned, with respect to the Lys47 amine, to act as the hydrolytic water in the deacylation step. Deacylation occurs with the formation of a tetrahedral intermediate over a 24 kcal x mol(-1) barrier. This barrier is approximately 2 kcal x mol(-1) greater than the barrier (22 kcal x mol(-1)) for the formation of the tetrahedral species in acylation. The potential-energy surface for the collapse of the deacylation tetrahedral species gives a 24 kcal x mol(-1) higher energy species for the product, signifying that the complex would readily reorganize and pave the way for the expulsion of the product of the reaction from the active site and the regeneration of the catalyst. These computational data dovetail with the knowledge on the reaction from experimental approaches.     

Molecular orbital studies of enzyme mechanisms. II. Catalytic oxidation of alcohols by liver alcohol dehydrogenase

Semiempirical (AM1) molecular orbital theory has been used to investigate the oxidation of alcohols at the active site of liver alcohol dehydrogenase (LADH). The model active site consists of a zinc dication coordinated to two methyl-mercaptans (Cys-46, Cys-176), an imidazole (His-67), and a water. An imidazole (His-51) hydrogen bonded to a hydroxy-acetate (Ser-48) forms the remote base. AM1 calculations that address the two distinct steps in the catalytic mechanism of ethanol oxidation by LADH are reported. These two steps are: (1) the deprotonation of ethanol by imidazole (His-51) via hydrogen-bonded hydroxy-acetate (Ser-48), creating a proton relay system; and (2) the rate-limiting hydride transfer step from ethanol C1 to nicotinamide adenine dinucleotide (NAD⁺), leading to product formation. Detailed calculations have been used to resolve the unsolved problems of mechanisms that have been suggested on the basis of kinetic data and crystal structures of several LADH complexes. We investigated two possible mechanisms for the deprotonation of ethanol, by zinc-bound OH^? and by direct deprotonation of zinc-bound ethanol by imidazole via hydroxyacetate (Ser-48). Our calculations show that there is no need for LADH to activate a water molecule at the active site as in many other zinc enzymes. This result agrees with experimental evidence. Our calculations also indicate that substrates are bound in an inner-sphere-pentacoordinated complex to the active site zincion. In this case, spectroscopic investigations agree with our results but crystallographic data do not. The highest activation energy is found for the hydride transfer, in agreement with the experiment. Finally, we proposed an alternative mechanism for the mode of action of LADH based upon our results. © 1993 John Wiley & Sons, Inc.     

Semiempirical Molecular Orbital Studies of the Acylation Step in the Lipase‐Catalyzed Ester Hydrolysis

In this study, we present the results from the semiempirical molecular orbital calculations for the acylation step in the lipase‐catalyzed ester hydrolysis. The results reveal that the lowest energy path for the formation of the tetrahedral intermediate is for the serine residue of the catalytic triad to attack the substrate, followed by coupling heavy atom movement and proton transfer. The calculations of four active site models show that the cooperation of the aspartate group and the oxyanion hole is capable of lowering the activation energy by about 16 kcalmol^?1. Our results further suggest that the lipase‐catalyzed ester hydrolysis adopts the single proton transfer mechanism.     

Binding of active site directed ligands to phospholipase A2: Implications on the molecular constraints and catalytic mechanism

Molecular constraints for the localization of active site directed ligands (competitive inhibitors and substrates) in the active site of phospholipase A2 (PLA2) are characterized. Structure activity relationships with known inhibitors suggest that the head group interactions dominate the selectivity as well as a substantial part of the affinity. Theab initio fitting of the amide ligands in the active site was carried out to characterize the head group interactions. Based on a systematic coordinate space search, formamide is docked with known experimental constraints such as coordination of the carbonyl group to Ca²⁺ and hydrogen bond between amide nitrogen and ND1 of His48. An optimal position for a bound water molecule is identified and its significance for the catalytic mechanism is postulated. Unlike the traditional “pseudo-triad” mechanism, the “Ca-coordinated-oxyanion” mechanism proposed here invokes activation of the catalytic water to form the oxyanion in the coordination sphere of calcium. As it attacks the carbonyl carbon of the ester, a near-tetrahedral intermediate is formed. As the second proton of the catalytic water is abstracted by the ester oxygen, its reorientation and simultaneous cleavage form hydrogen bond with ND1 of His48. In this mechanism of esterolysis, a catalytic role for the water co-ordinated to Ca²⁺ is recognised.     

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.     

Divergent function in the crotonase superfamily: an anhydride intermediate in the reaction catalyzed by 3-hydroxyisobutyryl-CoA hydrolase

3-Hydroxyisobutyryl-CoA hydrolase (HICH), a member of the enoyl-CoA (crotonase) superfamily, catalyzes the hydrolysis of 3-hydroxyisobutyryl-CoA to 3-hydroxyisobutyrate. Like other members of the superfamily, the sequence of HICH contains conserved sequences for an oxyanion hole that stabilizes anionic intermediates. In contrast to most members of the superfamily, the reaction catalyzed by HICH does not proceed via formation of a thioester enolate anion; instead, evidence based on substrate deuterium isotope effects, the reactivity of substrate analogues that cannot form thioester enolate anions, single-turnover experiments in H218O, and the kinetic phenotypes of site-directed mutants provide evidence for a mechanism involving the formation of an anhydride intermediate involving Glu143 in the active site. In the reactions catalyzed by many members of the superfamily, homologues of Glu143 abstract the alpha proton of the thioester substrate to generate the thioester enolate anion intermediate. Presumably, one or more of the anionic tetrahedral intermediates on the HICH reaction coordinate are stabilized by the oxyanion hole. Thus, we conclude that the conserved oxyanion hole in this superfamily can be used to stabilize a variety of anionic intermediates.     

Theoretical study of the role of low-barrier hydrogen bonds in enzyme catalysis: a model of proton transfer in serine protease

 A model of low-barrier hydrogen bonds (LBHBs) in enzymes has been studied by ab initio quantum mechanical calculations including the self-consistent reaction field solvent model. The hydrogen-bond strengths and the deprotonation energies for the hydrogen-bonded and non-hydrogen-bonded cis-urocanic acid were calculated at the HF/6-31 + G(d,p) level at various dielectric constants. The same calculations were performed for the α,β-dihydrourocanic acid to model the catalytic dyad of serine protease. The deprotonation energy of N^ɛ2 in α,β-dihydrourocanic acid is increased by formation of LBHBs and depends very much on the dielectric constant. This study suggests that the formation of LBHBs increases the basicity of the dyad, and the polarity change near the reaction center in the active site could help the proton abstraction from Ser 195 and the donation to the leaving group. Both the LBHBs and the environment can play crucial roles in the enzyme catalysis. Received: 8 March 2000 / Accepted: 3 January 2001 / Published online: 3 May 2001     

Discussion of the catalytic pathway of cysteine proteases based on AM1 calculations

The deacylation reaction of the cysteine protease papain was examined by AM1 reaction-coordinate calculations. The transition-state (TS) structure was extracted from the reaction pathway as corresponding to the maximum point along this minimum-energy pathway. Consistent with experimental kinetic data revealing that deacylation is about 60 times faster for thioester (---C(O)---S---) than dithioester (---C(S)---S---) intermediates, calculated E_a values are about 10 kcal mol⁻¹ lower for the former than the latter. The calculated partial atomic charges indicate that the C=O carbon in the thioester is a good site for nucleophilic attack whereas the corresponding C=S carbon in the dithioester is a poor site. The present calculations reveal that the enzyme's oxyanion hole contributes about 9 kcal mol⁻¹ toward reducing E_a for the anionic tetrahedral intermediate and TS structure. On the other hand, the presence of Asn in the putative Asn-His-Cys catalytic triad contributes only about l kcal mol⁻¹ toward reducing their E_a value. The presence of this Asn, however, did appear to stabilize His in its protonated form (ImH⁺) over its unprotonated form (Im). Two novel mechanisms are introduced to explain the unusual effect of a remote X substituent on the deacylation kinetics of the substrate family under consideration. The first mechanism invokes a “field effect” while the second mechanism embodies the concepts of induction and homoconjugation. The unique feature of these two mechanisms is that, unlike other proposed models, they circumvent the requirement for a close N…S interaction which has stimulated controversy.     

The Catalytic Mechanism of Fluoroacetate Dehalogenase: A Computational Exploration of Biological Dehalogenation

The biological dehalogenation of fluoroacetate carried out by fluoroacetate dehalogenase is discussed by using quantum mechanical/molecular mechanical (QM/MM) calculations for a whole‐enzyme model of 10 800 atoms. Substrate fluoroacetate is anchored by a hydrogen‐bonding network with water molecules and the surrounding amino acid residues of Arg105, Arg108, His149, Trp150, and Tyr212 in the active site in a similar way to haloalkane dehalogenase. Asp104 is likely to act as a nucleophile to attack the α‐carbon of fluoroacetate, resulting in the formation of an ester intermediate, which is subsequently hydrolyzed by the nucleophilic attack of a water molecule to the carbonyl carbon atom. The cleavage of the strong C? F bond is greatly facilitated by the hydrogen‐bonding interactions between the leaving fluorine atom and the three amino acid residues of His149, Trp150, and Tyr212. The hydrolysis of the ester intermediate is initiated by a proton transfer from the water molecule to His271 and by the simultaneous nucleophilic attack of the water molecule. The transition state and produced tetrahedral intermediate are stabilized by Asp128 and the oxyanion hole composed of Phe34 and Arg105.     

Inhibition of class A beta-lactamases by carbapenems: crystallographic observation of two conformations of meropenem in SHV-1

Carbapenem antibiotics are often the "last resort" in the treatment of infections caused by bacteria resistant to penicillins and cephalosporins. To understand why meropenem is resistant to hydrolysis by the SHV-1 class A beta-lactamase, the atomic structure of meropenem inactivated SHV-1 was solved to 1.05 A resolution. Two conformations of the Ser70 acylated intermediate are observed in the SHV-1-meropenem complex; the meropenem carbonyl oxygen atom of the acyl-enzyme is in the oxyanion hole in one conformation, while in the other conformation it is not. Although the structures of the SHV-1 apoenzyme and the SHV-1-meropenem complex are very similar (0.29 A rmsd for Calpha atoms), the orientation of the conserved Ser130 is different. Notably, the Ser130-OH group of the SHV-1-meropenem complex is directed toward Lys234Nz, while the Ser130-OH of the apo enzyme is oriented toward the Lys73 amino group. This altered position may affect proton transfer via Ser130 and the rate of hydrolysis. A most intriguing finding is the crystallographic detection of protonation of the Glu166 known to be involved in the deacylation mechanism. The critical deacylation water molecule has an additional hydrogen-bonding interaction with the OH group of meropenem's 6alpha-1 R-hydroxyethyl substituent. This interaction may weaken the nucleophilicity and/or change the direction of the lone pair of electrons of the water molecule and result in poor turnover of meropenem by the SHV-1 beta-lactamase. Using timed mass spectrometry, we further show that meropenem is covalently attached to SHV-1 beta-lactamase for at least 60 min. These observations explain key properties of meropenem's ability to resist hydrolysis by SHV-1 and lead to important insights regarding future carbapenem and beta-lactamase inhibitor design.     

Quinoprotein methanol dehydrogenase: a molecular dynamics study and comparison with crystal structure

One of the mechanisms proposed for methanol oxidation by methanol dehydrogenase (MDH) involves a hydride transfer to the quinone carbonyl carbon C5 of 2,7,9-tricarboxy-1H-pyrrolo[2,3-f]-quinoline-4,5-dione (PQQ), initiated by abstraction of a proton from the substrate methanol by a general base. Molecular dynamics studies are performed on MDH-bound to the C5 reduced intermediate (C5RI) of PQQ, for 3 ns. The structural features of the MD and X-ray structures are compared. An interesting feature observed during simulations is the strong hydrogen bond between oxyanion O5 of C5RI and Asp297-CO₂H in the active site. Asp297-CO₂⁻ is suggested to be the general-base catalyst for removing the hydroxyl proton of methanol in concert with the hydride ion transfer from the putative methoxide to C5 carbonyl of PQQ. The formed Asp297-CO₂H acts as the required proton source for the immediate product. Anticorreleated motions observed in the MD structure are not across the active site to influence the reaction mechanism of MDH.     

Solvent effects on the electronic absorption spectrum of formamide studied by a sequential Monte Carlo/quantum mechanical approach

Sequential Monte Carlo/quantum mechanical calculations are performed to study the solvent effects on the electronic absorption spectrum of formamide (FMA) in aqueous solution, varying from hydrogen bonds to the outer solvation shells. Full quantum-mechanical intermediate neglect of differential overlap/singly excited configuration interaction calculations are performed in the supermolecular structures generated by the Monte Carlo simulation. The largest calculation involves the ensemble average of 75 statistically uncorrelated quantum mechanical results obtained with the FMA solute surrounded by 150 water solvent molecules. We find that the n → π* transition suffers a blueshift of 1,600 cm⁻¹ upon solvation and the π → π* transition undergoes a redshift of 800 cm⁻¹. On average, 1.5 hydrogen bonds are formed between FMA and water and these contribute with about 20% and about 30% of the total solvation shifts of the n → π* and π → π* transitions, respectively. The autocorrelation function of the energy is used to sample configurations from the Monte Carlo simulation, and the solvation shifts are shown to be converged values. Received: 14 March 2002 / Accepted: 3 April 2002 / Published online: 24 June 2002     

A theoretical approach for the acylation/deacylation mechanisms of avibactam in the reversible inhibition of KPC-2

Klebsiella pneumoniae carbapenemase (KPC-2) is the most commonly encountered class A β-lactamase variant worldwide, which confer high-level resistance to most available antibiotics. In this article we address the issue by a combined approach involving molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations. The study contributes to improve the understanding, at molecular level, of the acylation and deacylation stages of avibactam involved in the inhibition of KPC-2. The results show that both mechanisms, acylation and deacylation, the reaction occur via the formation of a tetrahedral intermediate. The formation of this intermediate corresponds to the rate limiting stage. The activation barriers are 19.5 kcal/mol and 23.0 kcal/mol for the acylation and deacylation stages, respectively. The associated rate constants calculated, using the Eyring equation, are 1.2?×?10^?1 and 3.9?×?10^?4 (s^?1). These values allow estimating a value of 3.3?×?10^?3 for the inhibition constant, in good agreement with the experimental value.

     

Acylation and deacylation mechanism and kinetics of penicillin G reaction with Streptomyces R61 DD-peptidase

Two quantum mechanical (QM)-cluster models are built for studying the acylation and deacylation mechanism and kinetics of Streptomyces R61 DD-peptidase with the penicillin G at atomic level detail. DD-peptidases are bacterial enzymes involved in the cross-linking of peptidoglycan to form the cell wall, necessary for bacterial survival. The cross-linking can be inhibited by antibiotic beta-lactam derivatives through acylation, preventing the acyl-enzyme complex from undergoing further deacylation. The deacylation step was predicted to be rate-limiting. Transition state and intermediate structures are found using density functional theory in this study, and thermodynamic and kinetic properties of the proposed mechanism are evaluated. The acyl-enzyme complex is found lying in a deep thermodynamic sink, and deacylation is indeed the severely rate-limiting step, leading to suicide inhibition of the peptidoglycan cross-linking. The usage of QM-cluster models is a promising technique to understand, improve, and design antibiotics to disrupt function of the Streptomyces R61 DD-peptidase.     

Enzymatic Ring-Opening Polymerization of Lactones by Lipase Catalyst: Mechanistic Aspects

Mechanistic aspects of lipase-catalyzed ring-opening polymerization (ROP) of lactones to give polyesters are discussed from accumulated experimental data and new insight. Comparison of the ROP reactivity by lipase catalyst with the anionic ROP reactivity by a metal-catalyst clearly demonstrates the characteristics of lipase catalysis; the larger ring-sized monomers with lower ring strain showed higher polymerizability than medium ring-sized ones, in contrast to the anionic ROP showing the reverse direction where the ring strain of monomer is operative. The enzyme-catalysis involves an acyl-enzyme intermediate formation as a key-step. From the copolymerization results a new mechanism is proposed, that involves the formation of the acyl-enzyme intermediate (acylation step) and/or the nucleophilic attack of the propagationg alcohol end to the carbonyl carbon of the intermediate to open the monomer ring (deacylation step) as the rate-determining step. The structure of the propagating alcohol end (primary or secondary) affects much on which step is more operative.