A Theoretical Study on the Substrate Deacylation Mechanism of Class C β-Lactamase

Development of QM/MM (ABEEM polarizable force field) method to simulate the excision reaction mechanism of damaged cytosine

DNA damages are regarded as having harmful effects on cell. The base excision repair mechanism combats these effects by removing damaged bases. The deglycosylation mechanism of excising damaged bases by DNA glycosylase and the state of the leaving base have been controversial. The enzymatic reaction of DNA glycosylase to remove the damaged bases involves not only the formation and breaking of chemical bonds, but also complex polarization effect and charge transfer, which cannot be accurately simulated by the QM/MM method combined with the fixed charge force field. This work has developed the ABEEM fluctuating polarizable force field combining with the QM method, that is (QM/MM[ABEEM]), to accurately simulate the proton transfer, charge transfer and the charge distribution. The piecewise function is used as the valence-state electronegativity in the QM/MM (ABEEM) to realize the accurate fitting of the charge distribution in reaction. And the charge transfer is accurately simulated by the local charge conservation conditions. Four deglycosylation mechanisms including the monofunctional and difunctional mechanisms of four neutral and protonated cytosine derivatives are explored. It is confirmed that the monofunctional mechanism of Asp-activated nucleophile water is a better deglycosylation mechanism and the base is protonated before the reaction occurs. Protonization of the base reduced the activation energy by 10.00-17.00 kcal/mol. Asp provides the necessary charge for the reaction, and DNA glycosylase preferentially cleaves ɛC. This work provides a theoretical basis for the research of excising damaged bases by DNA glycosylase.

Secondary Structure and Conformational Stability of the Antigen Residues Making Contact with Antibodies

Antibodies are crucial biomolecules that bring high therapeutic efficacy in medicine and accurate molecular detection in diagnosis. Many studies have been devoted to analyzing the antigen-antibody interaction from the importance of understanding the antibody recognition mechanism. However, most of the previous studies examined the characteristic of the antibody for interaction. It is also informative to clarify the significant antigen residues contributing to the binding. To characterize the molecular interaction of antigens, we computationally analyzed 350 antigen-antibody complex structures by molecular mechanics (MM) calculations and molecular dynamics (MD) simulations. Based on the MM calculations, the antigen residues contributing to the binding were extracted from all the 350 complexes. The extracted residues are located at the antigen-antibody interface and are responsible for making contact with the antibody. The appearances of the charged polar residues, Asp, Glu, Arg, and Lys, were noticeably large. In contrast, the populations of the hydrophobic residues, Leu, Val, and Ala, were relatively low. The appearance frequencies of the other amino acid residues were almost close to the abundance of general proteins of eukaryotes. The binding score indicated that the hydrophilic interaction was dominant at the antigen-antibody contact instead of the hydrophobic one. The positively charged residues, Arg and Lys, remarkably contributed to the binding compared to the negatively charged ones, Asp and Glu. Considerable contributions were also observed for the noncharged polar residues, Asn and Gln. The analysis of the secondary structures of the extracted antigen residues suggested that there was no marked difference in recognition by antibodies among helix, sheet, turn, and coil. A long helix of the antigen sometimes made contact with antibody complementarity-determining regions, and a large sheet also frequently covered the antibody heavy and light chains. The turn structure was the most popularly observed at the contact with antibody among 350 complexes. Three typical complexes were picked up for each of the four secondary structures. MD simulations were performed to examine the stability of the interfacial structures of the antigens for these 12 complex models. The alterations of secondary structures were monitored through the simulations. The structural fluctuations of the contact residues were low compared with the other domains of antigen molecules. No drastic conversion was observed for every model during the 100 ns simulation. The motions of the interfacial antigen residues were small compared to the other residues on the protein surface. Therefore, diverse molecular conformations are possible for antibody recognition as long as the target areas are polar, nonflexible, and protruding on the protein surface.

Systematic Parameterization and Simulation of Boronic Acid-β-Lactamase Aqueous Solution in Developing the ABEEMσπ Polarizable Force Field

Boronic acid, an inhibitor of β-lactamase, has begun to be applied to the treatment of biological infections and tumors. Scientists are working to develop new and more effective boronic acid. Molecular dynamics (MD) simulation provides a powerful auxiliary tool for drug design. However, the current force fields have no boron-related parameters. In this work, an atom-bond electronegativity equalization method at the σπ level (ABEEMσπ) polarizable force field (ABEEMσπ PFF) of boronic acid and β-lactamase has been developed to determine the potential functions and parameters. The interaction between boron and serine in β-lactamase is regarded as a bonded mode. The interaction between them is simulated by the Morse potential energy function, which is close to the experimental change of the stretching potential energy in a large range. The potential energy surfaces of the bond length, bond angle, and dihedral angle of boronic acid-β-lactamase have the same stability point and change trend as M06-2X/6-311G**. For 47 boronic acid-β-lactamase training molecules, the linear correlation coefficient (R) of the charge distribution between the ABEEMσπ PFF and HF/STO-3G is greater than 0.96. Attributed to the fact that the charge distribution of the ABEEMσπ PFF can fluctuate with the change of geometry and environment, the polarization effect and charge-transfer effect are well reflected. The binding ability of different boronic acids with the same β-lactamase is different. A total of 10 boronic acid-β-lactamase model molecules and 10 boronic acid-β-lactamase and water complexes are simulated. The order of binding energy of five large model molecules calculated by the ABEEMσπ PFF is consistent with that of the MP2 method. The binding energies of boronic acid-β-lactamase and water complexes are close to those of the MP2 method. The results of MD simulation of five aqueous boronic acid-β-lactamase complexes in the NVT ensemble verify the rationality of boron-related parameters of the ABEEMσπ PFF, which have a good application prospect. This study lays a solid theoretical foundation for further study of the inhibition of boronic acid on β-lactamase.

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.

Deacylation Mechanism and Kinetics of Acyl-Enzyme Complex of Class C β-Lactamase and Cephalothin

Understanding the molecular details of antibiotic resistance by the bacterial enzymes β-lactamases is vital for the development of novel antibiotics and inhibitors. In this spirit, the detailed mechanism of deacylation of the acyl-enzyme complex formed by cephalothin and class C β-lactamase is investigated here using hybrid quantum-mechanical/molecular-mechanical molecular dynamics methods. The roles of various active-site residues and substrate in the deacylation reaction are elucidated. We identify the base that activates the hydrolyzing water molecule and the residue that protonates the catalytic serine (Ser64). Conformational changes in the active sites and proton transfers that potentiate the efficiency of the deacylation reaction are presented. We have also characterized the oxyanion holes and other H-bonding interactions that stabilize the reaction intermediates. Together with the kinetic and mechanistic details of the acylation reaction, we analyze the complete mechanism and the overall kinetics of the drug hydrolysis. Finally, the apparent rate-determining step in the drug hydrolysis is scrutinized.

Structure of a highly acidic β-lactamase from the moderate halophile Chromohalobacter sp. 560 and the discovery of a Cs(+)-selective binding site

Environmentally friendly absorbents are needed for Sr(2+) and Cs(+), as the removal of the radioactive Sr(2+) and Cs(+) that has leaked from the Fukushima Nuclear Power Plant is one of the most important problems in Japan. Halophilic proteins are known to have many acidic residues on their surface that can provide specific binding sites for metal ions such as Cs(+) or Sr(2+). The crystal structure of a halophilic β-lactamase from Chromohalobacter sp. 560 (HaBLA) was determined to resolutions of between 1.8 and 2.9 Å in space group P31 using X-ray crystallography. Moreover, the locations of bound Sr(2+) and Cs(+) ions were identified by anomalous X-ray diffraction. The location of one Cs(+)-specific binding site was identified in HaBLA even in the presence of a ninefold molar excess of Na(+) (90 mM Na(+)/10 mM Cs(+)). From an activity assay using isothermal titration calorimetry, the bound Sr(2+) and Cs(+) ions do not significantly affect the enzymatic function of HaBLA. The observation of a selective and high-affinity Cs(+)-binding site provides important information that is useful for the design of artificial Cs(+)-binding sites that may be useful in the bioremediation of radioactive isotopes.

A triple mutant in the Ω-loop of TEM-1 β-lactamase changes the substrate profile via a large conformational change and an altered general base for catalysis

β-Lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. TEM-1 is a prevalent plasmid-encoded β-lactamase in Gram-negative bacteria that efficiently catalyzes the hydrolysis of penicillins and early cephalosporins but not oxyimino-cephalosporins. A previous random mutagenesis study identified a W165Y/E166Y/P167G triple mutant that displays greatly altered substrate specificity with increased activity for the oxyimino-cephalosporin, ceftazidime, and decreased activity toward all other β-lactams tested. Surprisingly, this mutant lacks the conserved Glu-166 residue critical for enzyme function. Ceftazidime contains a large, bulky side chain that does not fit optimally in the wild-type TEM-1 active site. Therefore, it was hypothesized that the substitutions in the mutant expand the binding site in the enzyme. To investigate structural changes and address whether there is an enlargement in the active site, the crystal structure of the triple mutant was solved to 1.44 Å. The structure reveals a large conformational change of the active site Ω-loop structure to create additional space for the ceftazidime side chain. The position of the hydroxyl group of Tyr-166 and an observed shift in the pH profile of the triple mutant suggests that Tyr-166 participates in the hydrolytic mechanism of the enzyme. These findings indicate that the highly conserved Glu-166 residue can be substituted in the mechanism of serine β-lactamases. The results reveal that the robustness of the overall β-lactamase fold coupled with the plasticity of an active site loop facilitates the evolution of enzyme specificity and mechanism.

Mechanism of acyl-enzyme complex formation from the Henry-Michaelis complex of class C β-lactamases with β-lactam antibiotics

Bacteria that cause most of the hospital-acquired infections make use of class C β-lactamase (CBL) among other enzymes to resist a wide spectrum of modern antibiotics and pose a major public health concern. Other than the general features, details of the defensive mechanism by CBL, leading to the hydrolysis of drug molecules, remain a matter of debate, in particular the identification of the general base and role of the active site residues and substrate. In an attempt to unravel the detailed molecular mechanism, we carried out extensive hybrid quantum mechanical/molecular mechanical Car-Parrinello molecular dynamics simulation of the reaction with the aid of the metadynamics technique. On this basis, we report here the mechanism of the formation of the acyl-enzyme complex from the Henry-Michaelis complex formed by β-lactam antibiotics and CBL. We considered two β-lactam antibiotics, namely, cephalothin and aztreonam, belonging to two different subfamilies. A general mechanism for the formation of a β-lactam antibiotic-CBL acyl-enzyme complex is elicited, and the individual roles of the active site residues and substrate are probed. The general base in the acylation step has been identified as Lys67, while Tyr150 aids the protonation of the β-lactam nitrogen through either the substrate carboxylate group or a water molecule.

Investigation of the acylation mechanism of class C beta-lactamase: pKa calculation, molecular dynamics simulation and quantum mechanical calculation

β-Lactamases are bacterial enzymes that act as a bacterial defense system against β-lactam antibiotics. β-Lactamase cleaves the β-lactam ring of the antibiotic by a two step mechanism involving acylation and deacylation steps. Although class C β-lactamases have been investigated extensively, the details of their mechanism of action are not well understood at the molecular level. In this study, we investigated the mechanism of the acylation step of class C β-lactamase using pKa calculations, molecular dynamics (MD) simulations and quantum mechanical (QM) calculations. Serine64 (Ser64) is an active site residue that attacks the β-lactam ring. In this study, we considered three possible scenarios for activation of the nucleophile Ser64, where the activation base is (1) Tyrosine150 (Tyr150), (2) Lysine67 (Lys67), or (3) substrate. From the pKa calculation, we found that Tyr150 and Lys67 are likely to remain in their protonated states in the pre-covalent complex between the enzyme and substrate, although their role as activator would require them to be in the deprotonated state. It was found that the carboxylate group of the substrate remained close to Ser64 for most of the simulation. The energy barrier for hydrogen abstraction from Ser64 by the substrate was calculated quantum mechanically using a large truncated model of the enzyme active site and found to be close to the experimental energy barrier, which suggests that the substrate can initiate the acylation mechanism in class C β-lactamase.

Ab initio protein structure prediction with force field parameters derived from water-phase quantum chemical calculation

Molecular dynamics (MD) simulations are extensively used in the study of the structures and functions of proteins. Ab initio protein structure prediction is one of the most important subjects in computational biology, and many trials have been performed using MD simulation so far. Since the results of MD simulations largely depend on the force field, reliable force field parameters are indispensable for the success of MD simulation. In this work, we have modified atom charges in a standard force field on the basis of water-phase quantum chemical calculations. The modified force field turned out appropriate for ab initio protein structure prediction by the MD simulation with the generalized Born method. Detailed analysis was performed in terms of the conformational stability of amino acid residues, the stability of secondary structure of proteins, and the accuracy for prediction of protein tertiary structure, comparing the modified force field with a standard one. The energy balance between alpha-helix and beta-sheet structures was significantly improved by the modification of charge parameters.

Substrate deacylation mechanisms of serine-beta-lactamases

The substrate deacylation mechanisms of serine-beta-lactamases (classes A, C and D) were investigated by theoretical calculations. The deacylation of class A proceeds via four elementary reactions. The rate-determining process is the tetrahedral intermediate (TI) formation and the activation energy is 24.6 kcal/mol at the DFT level. The deacylation does not proceed only by Glu166, which acts as a general base, but Lys73 also participates in the reaction. The C3-carboxyl group of the substrate reduces the barrier height at the TI formation (substrate-assisted catalysis). In the case of class C, the deacylation consists of two elementary processes. The activation energy of the TI formation has been estimated to be 30.5 kcal/mol. Tyr150Oeta is stabilized in the deprotonated state in the acyl-enzyme complex and works as a general base. This situation can exist due to the interaction with two positively charged side chains of lysine (Lys67 and Lys315). The deacylation of class D also consists of two elementary reaction processes. The activation energy of the TI formation is ca. 30 kcal/mol. It is thought that the side chain of Lys70 is deprotonated and acts as a general base. When Lys70 is carbamylated, the activation energy is reduced to less than 20 kcal/mol. This suggests that the high hydrolysis activity of class D with carbamylated Lys70 is due to the reduction of activation energy for deacylation. From these results, it is concluded that the contribution of the lysine residue adjacent to the serine residue is indispensable for the enzymatic reactions by serine-beta-lactamases.