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

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Origin of pKa Shifts of Internal Lysine Residues in SNase Studied Via Equal-Molar VMMS Simulations in Explicit Water

Protein internal ionizable groups can exhibit large shifts in pK_a values. Although the environment and interaction changes have been extensively studied both experimentally and computationally, direct calculation of pK_a values of these internal ionizable groups in explicit water is challenging due to energy barriers in solvent interaction and in conformational transition. The virtual mixture of multiple states (VMMS) method is a new approach designed to study chemical state equilibrium. This method constructs a virtual mixture of multiple chemical states in order to sample the conformational space of all states simultaneously and to avoid crossing energy barriers related to state transition. By applying VMMS to 25 variants of staphylococcal nuclease with lysine residues at internal positions, we obtained the pK_a values of these lysine residues and investigated the physics underlining the pK_a shifts. Our calculation results agree reasonably well with experimental measurements, validating the VMMS method for pK_a calculation and providing molecular details of the protonation equilibrium for protein internal ionizable groups. Based on our analyses of protein conformation relaxation, lysine side chain flexibility, water penetration, and the microenvironment, we conclude that the hydrophobicity of the microenvironment around the lysine side chain (which affects water penetration differently for different protonation states) plays an important role in the pK_a shifts.

Free Energy Contribution Analysis Using Response Kernel Approximation: Insights into the Acylation Reaction of a Beta-Lactamase

A widely applicable free energy contribution analysis (FECA) method based on the quantum mechanical/molecular mechanical (QM/MM) approximation using response kernel approaches has been proposed to investigate the influences of environmental residues and/or atoms in the QM region on the free energy profile. This method can evaluate atomic contributions to the free energy along the reaction path including polarization effects on the QM region within a dramatically reduced computational time. The rate-limiting step in the deactivation of the β-lactam antibiotic cefalotin (CLS) by β-lactamase was studied using this method. The experimentally observed activation barrier was successfully reproduced by free energy perturbation calculations along the optimized reaction path that involved activation by the carboxylate moiety in CLS. It was found that the free energy profile in the QM region was slightly higher than the isolated energy and that two residues, Lys67 and Lys315, as well as water molecules deeply influenced the QM atoms associated with the bond alternation reaction in the acyl-enzyme intermediate. These facts suggested that the surrounding residues are favorable for the reactant complex and prevent the intermediate from being too stabilized to proceed to the following deacylation reaction. We have demonstrated that the free energy contribution analysis should be a useful method to investigate enzyme catalysis and to facilitate intelligent molecular design.

Nylon-Oligomer Hydrolase Promoting Cleavage Reactions in Unnatural Amide Compounds

The active site of 6-aminohexanoate-dimer hydrolase, a nylon-6 byproduct-degrading enzyme with a β-lactamase fold, possesses a Ser112/Lys115/Tyr215 catalytic triad similar to the one of penicillin-recognizing family of serine-reactive hydrolases but includes a unique Tyr170 residue. By using a reactive quantum mechanics/molecular mechanics (QM/MM) approach, we work out its catalytic mechanism and related functional/structural specificities. At variance with other peptidases, we show that the involvement of Tyr170 in the enzyme-substrate interactions is responsible for a structural variation in the substrate-binding state. The acylation via a tetrahedral intermediate is the rate-limiting step, with a free-energy barrier of ∼21 kcal/mol, driven by the catalytic triad Ser112, Lys115, and Tyr215, acting as a nucleophile, general base, and general acid, respectively. The functional interaction of Tyr170 with this triad leads to an efficient disruption of the tetrahedral intermediate, promoting a conformational change of the substrate favorable for proton donation from the general acid.

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.