An official website of the United States government
Pathogen resistance to β-lactam antibiotics compromises effective treatments of superbug infections. One major source of β-lactam resistance is the bacterial production of β-lactamases, which could effectively hydrolyze β-lactam drugs. In this thesis, the hydrolysis of various β-lactam antibiotics by class A serine-based β-lactamases (ASβLs) were investigated using hybrid Quantum Mechanical / Molecular Mechanical (QM/MM) minimum energy pathway (MEP) calculations and explainable machine learning (ML) approaches. The TEM-1/benzylpenicillin acylation reaction with QM/MM chain-of-states reaction pathways was firstly revisited. I proposed two decomposition methods for energy contribution analysis based on perturbing ML regression models. Both methods were shown to be model implementation invariant and successfully bridged the discrepancies between two pioneering mechanistic studies. The Toho-1 ASβL acylations of ampicillin and cefalexin were then investigated. I reported that the acylation pathway selection can be ligand dependent: ampicillin could undergo acylation via Lys73 or Glu166 acting as the general base while cefalexin acylation is limited to Lys73 as the general base. An explainable artificial intelligence (XAI) method, the Boltzmann-weighted Cumulative Integrated Gradients (BCIG), was developed to explain the different acylation pathway viability found for ampicillin and cefalexin. Lastly, conformational factors determining the GES-5/imipenem deacylation activity was investigated using edge-conditioned convolutional graph-learning (GL) methods. Critical mechanistic insights were derived from perturbative response of the GL latent representations, which explained the different deacylation reactivity between the two imipenem pyrroline tautomer states and identified the orientation of the carbapenem 6α-hydroxyethyl as the key factor that impacts the deacylation barrier heights. In summary, my thesis focuses on bridging QM/MM chain-of-states reaction pathway calculations and explainable ML to derive essential mechanistic insights into β-lactam resistance driven by ASβLs.
more »
« less
Unraveling the energetic significance of chemical events in enzyme catalysis via machine-learning based regression approach
Abstract The bacterial enzyme class of β-lactamases are involved in benzylpenicillin acylation reactions, which are currently being revisited using hybrid quantum mechanical molecular mechanical (QM/MM) chain-of-states pathway optimizations. Minimum energy pathways are sampled by reoptimizing pathway geometry under different representative protein environments obtained through constrained molecular dynamics simulations. Predictive potential energy surface models in the reaction space are trained with machine-learning regression techniques. Herein, using TEM-1/benzylpenicillin acylation reaction as the model system, we introduce two model-independent criteria for delineating the energetic contributions and correlations in the predicted reaction space. Both methods are demonstrated to effectively quantify the energetic contribution of each chemical process and identify the rate limiting step of enzymatic reaction with high degrees of freedom. The consistency of the current workflow is tested under seven levels of quantum chemistry theory and three non-linear machine-learning regression models. The proposed approaches are validated to provide qualitative compliance with experimental mutagenesis studies.
more »
« less
- Award ID(s):
- 1753167
- PAR ID:
- 10197328
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Communications Chemistry
- Volume:
- 3
- Issue:
- 1
- ISSN:
- 2399-3669
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Chemical energy ferroelectrics are generally solid macromolecules showing spontaneous polarization and chemical bonding energy. These materials still suffer drawbacks, including the limited control of energy release rate, and thermal decomposition energy well below total chemical energy. To overcome these drawbacks, we report the integrated molecular ferroelectric and energetic material from machine learning-directed additive manufacturing coupled with the ice-templating assembly. The resultant aligned porous architecture shows a low density of 0.35 g cm−3, polarization-controlled energy release, and an anisotropic thermal conductivity ratio of 15. Thermal analysis suggests that the chlorine radicals react with macromolecules enabling a large exothermic enthalpy of reaction (6180 kJ kg−1). In addition, the estimated detonation velocity of molecular ferroelectrics can be tuned from 6.69 ± 0.21 to 7.79 ± 0.25 km s−1by switching the polarization state. These results provide a pathway toward spatially programmed energetic ferroelectrics for controlled energy release rates.more » « less
-
Lian, T (Ed.)The fast and accurate simulation of chemical reactions is a major goal of computational chemistry. Recently, the pursuit of this goal has been aided by machine learning interatomic potentials (MLIPs), which provide energies and forces at quantum mechanical accuracy but at a fraction of the cost of the reference quantum mechanical calculations. Assembling the training set of relevant configurations is key to building the MLIP. Here, we demonstrate two approaches to training reactive MLIPs based on reaction pathway information. One approach exploits reaction datasets containing reactant, product, and transition state structures. Using an SN2 reaction dataset, we accurately locate reaction pathways and transition state geometries of up to 170 unseen reactions. In another approach, which does not depend on data availability, we present an efficient active learning procedure that yields an accurate MLIP and converged minimum energy path given only the reaction end point structures, avoiding quantum mechanics driven reaction pathway search at any stage of training set construction. We demonstrate this procedure on an SN2 reaction in the gas phase and with a small number of solvating water molecules, predicting reaction barriers within 20 meV of the reference quantum chemistry method. We then apply the active learning procedure on a more complex reaction involving a nucleophilic aromatic substitution and proton transfer, comparing the results against the reactive ReaxFF force field. Our active learning procedure, in addition to rapidly finding reaction paths for individual reactions, provides an approach to building large reaction path databases for training transferable reactive machine learning potentials.more » « less
-
Efficient mechanism-based design of antibiotics that are not susceptible to β-lactamases is hindered by the lack of comprehensive knowledge on the energetic landscapes for the hydrolysis of various β-lactams. Herein, we adopted efficient quantum mechanics/molecular mechanics simulations to explore the acylation reaction catalyzed by CTX-M-44 (Toho-1) β-lactamase. We show that the catalytic pathways for β-lactam hydrolysis are correlated to substrate scaffolds: using Glu166 as the only general base for acylation is viable for ampicillin but prohibitive for cefalexin. The present computational workflow provides quantitative insights to facilitate the optimization of future β-lactam antibiotics.more » « less
-
We present a thermodynamically consistent energetic variational model for active nematics driven by ATP hydrolysis. Extending the classical Toner–Tu framework, we introduce a chemo-mechanical coupling mechanism in which the self-advection and polarization dynamics are modulated by the ATP hydrolysis rate. The model is derived using an energetic variational approach that integrates both chemical free energy and mechanical energy into a unified energy dissipation law. The reaction rate equation explicitly incorporates mechanical feedback, revealing how active transport and alignment interactions influence chemical fluxes and vice versa. This formulation not only preserves consistency with non-equilibrium thermodynamics but also provides a transparent pathway for modeling energy transduction in active systems. We also present numerical simulations demonstrating the positive energy transduction under a specific choice of model parameters. The new modeling framework offers new insights into energy transduction and regulation mechanisms in biologically related active systems.more » « less
