Molecular dynamics of class A β-lactamases-effects of substrate binding

The Evolution of the Acylation Mechanism in β-Lactamase and Rapid Protein Dynamics

β-Lactamases are a class of well-studied enzymes that are known to have existed since billions of years ago, starting as a defense mechanism to stave off competitors and are now enzymes responsible for antibiotic resistance. Using ancestral sequence reconstruction, it is possible to study the crystal structure of a laboratory resurrected 2-3 billion year-old β-lactamase. Comparing the ancestral enzyme to its modern counterpart, a TEM-1 β-lactamase, the structural changes are minor, and it is probable that dynamic effects play an important role in the evolution of function. We used molecular dynamics simulations and employed transition path sampling methods to identify the presence of rate-enhancing dynamics at the femtosecond level in both systems, found that these fast motions are more efficiently coordinated in the modern enzyme, and examined how specific dynamics can pinpoint evolutionary effects that are essential for improving enzymatic catalysis.

Molecular Evolutionary Analyses of the Pseudomonas-Derived Cephalosporinase Gene

Despite the increasing evidence of the clinical impact of Pseudomonas-derived cephalosporinase (PDC) sequence polymorphisms, the molecular evolution of its encoding gene, blaPDC, remains elusive. To elucidate this, we performed a comprehensive evolutionary analysis of blaPDC. A Bayesian Markov Chain Monte Carlo phylogenetic tree revealed that a common ancestor of blaPDC diverged approximately 4660 years ago, leading to the formation of eight clonal variants (clusters A-H). The phylogenetic distances within clusters A to G were short, whereas those within cluster H were relatively long. Two positive selection sites and many negative selection sites were estimated. Two PDC active sites overlapped with negative selection sites. In docking simulation models based on samples selected from clusters A and H, piperacillin was bound to the serine and the threonine residues of the PDC active sites, with the same binding mode for both models. These results suggest that, in P. aeruginosa, blaPDC is highly conserved, and PDC exhibits similar antibiotic resistance functionality regardless of its genotype.

Detailed investigation of catalytically important residues of class A β-lactamase

An increasing global health challenge is antimicrobial resistance. Bacterial infections are often treated by using β-lactam antibiotics. But several resistance mechanisms have evolved in clinically mutated bacteria, which results in resistance against such antibiotics. Among which production of novel β-lactamase is the major one. This results in bacterial resistance against penicillin, cephalosporin, and carbapenems, which are considered to be the last resort of antibacterial treatment. Hence, β-lactamase enzymes produced by such bacteria are called extended-spectrum β-lactamase and carbapenemase enzymes. Further, these bacteria have developed resistance against many β-lactamase inhibitors as well. So, investigation of important residues that play an important role in altering and expanding the spectrum activity of these β-lactamase enzymes becomes necessary. This review aims to gather knowledge about the role of residues and their mutations in class A β-lactamase, which could be responsible for β-lactamase mediated resistance. Class A β-lactamase enzymes contain most of the clinically significant and expanded spectrum of β-lactamase enzymes. Ser70, Lys73, Ser130, Glu166, and Asn170 residues are mostly conserved and have a role in the enzyme's catalytic activity. In-depth investigation of 69, 130, 131, 132, 164, 165, 166, 170, 171, 173, 176, 178, 179, 182, 237, 244, 275 and 276 residues were done along with its kinetic analysis for knowing its significance. Further, detailed information from many previous studies was gathered to know the effect of mutations on the kinetic activity of class A β-lactamase enzymes with β-lactam antibiotics.Communicated by Ramaswamy H. Sarma.

On the force field optimisation of β -lactam cores using the force field Toolkit

When employing molecular dynamics (MD) simulations for computer-aided drug design, the quality of the used force fields is highly important. Here we present reparametrisations of the force fields for the core molecules from 9 differentβ -lactam classes, for which we utilized the force field Toolkit and Gaussian calculations. We focus on the parametrisation of the dihedral angles, with the goal of reproducing the optimised quantum geometry in MD simulations. Parameters taken from CGenFF turn out to be a good initial guess for the multiplicity of each dihedral angle, but the key to a successful parametrisation is found to lie in the phase shifts. Based on the optimised quantum geometry, we come up with a strategy for predicting the phase shifts prior to the dihedral potential fitting. This allows us to successfully parameterise 8 out of the 11 molecules studied here, while the remaining 3 molecules can also be parameterised with small adjustments. Our work highlights the importance of predicting the dihedral phase shifts in the ligand parametrisation protocol, and provides a simple yet valuable strategy for improving the process of parameterising force fields of drug-like molecules.

The Role of Hydrophobic Nodes in the Dynamics of Class A β-Lactamases

Class A β-lactamases are known for being able to rapidly gain broad spectrum catalytic efficiency against most β-lactamase inhibitor combinations as a result of elusively minor point mutations. The evolution in class A β-lactamases occurs through optimisation of their dynamic phenotypes at different timescales. At long-timescales, certain conformations are more catalytically permissive than others while at the short timescales, fine-grained optimisation of free energy barriers can improve efficiency in ligand processing by the active site. Free energy barriers, which define all coordinated movements, depend on the flexibility of the secondary structural elements. The most highly conserved residues in class A β-lactamases are hydrophobic nodes that stabilize the core. To assess how the stable hydrophobic core is linked to the structural dynamics of the active site, we carried out adaptively sampled molecular dynamics (MD) simulations in four representative class A β-lactamases (KPC-2, SME-1, TEM-1, and SHV-1). Using Markov State Models (MSM) and unsupervised deep learning, we show that the dynamics of the hydrophobic nodes is used as a metastable relay of kinetic information within the core and is coupled with the catalytically permissive conformation of the active site environment. Our results collectively demonstrate that the class A enzymes described here, share several important dynamic similarities and the hydrophobic nodes comprise of an informative set of dynamic variables in representative class A β-lactamases.

Toho-1 β-lactamase: backbone chemical shift assignments and changes in dynamics upon binding with avibactam

Backbone chemical shift assignments for the Toho-1 β-lactamase (263 amino acids, 28.9 kDa) are reported based on triple resonance solution-state NMR experiments performed on a uniformly 2H,13C,15N-labeled sample. These assignments allow for subsequent site-specific characterization at the chemical, structural, and dynamical levels. At the chemical level, titration with the non-β-lactam β-lactamase inhibitor avibactam is found to give chemical shift perturbations indicative of tight covalent binding that allow for mapping of the inhibitor binding site. At the structural level, protein secondary structure is predicted based on the backbone chemical shifts and protein residue sequence using TALOS-N and found to agree well with structural characterization from X-ray crystallography. At the dynamical level, model-free analysis of 15N relaxation data at a single field of 16.4 T reveals well-ordered structures for the ligand-free and avibactam-bound enzymes with generalized order parameters of ~ 0.85. Complementary relaxation dispersion experiments indicate that there is an escalation in motions on the millisecond timescale in the vicinity of the active site upon substrate binding. The combination of high rigidity on short timescales and active site flexibility on longer timescales is consistent with hypotheses for achieving both high catalytic efficiency and broad substrate specificity: the induced active site dynamics allows variously sized substrates to be accommodated and increases the probability that the optimal conformation for catalysis will be sampled.

Allosteric communication in class A β-lactamases occurs via cooperative coupling of loop dynamics

Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A β-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.

β-Lactamase of Mycobacterium tuberculosis Shows Dynamics in the Active Site That Increase upon Inhibitor Binding

The Mycobacterium tuberculosis β-lactamase BlaC is a broad-spectrum β-lactamase that can convert a range of β-lactam antibiotics. Enzymes with low specificity are expected to exhibit active-site flexibility. To probe the motions in BlaC, we studied the dynamic behavior in solution using nuclear magnetic resonance (NMR) spectroscopy. ¹⁵N relaxation experiments show that BlaC is mostly rigid on the pico- to nanosecond timescale. Saturation transfer experiments indicate that also on the high-millisecond timescale BlaC is not dynamic. Using relaxation dispersion experiments, clear evidence was obtained for dynamics in the low-millisecond range, with an exchange rate of ca. 860 s^-1 The dynamic amide groups are localized in the active site. Upon formation of an adduct with the inhibitor avibactam, extensive line broadening occurs, indicating an increase in magnitude of the active-site dynamics. Furthermore, the rate of the motions increases significantly. Upon reaction with the inhibitor clavulanic acid, similar line broadening is accompanied by duplication of NMR signals, indicative of at least one additional, slower exchange process (exchange rate, k _ex, of <100 s^-1), while for this inhibitor also loss of pico- to nanosecond timescale rigidity is observed for some amides in the α domain. Possible sources of the observed dynamics, such as motions in the omega loop and rearrangements of active-site residues, are discussed. The increase in dynamics upon ligand binding argues against a model of inhibitor binding through conformational selection. Rather, the induced dynamics may serve to maximize the likelihood of sampling the optimal conformation for hydrolysis of the bound ligand.

The Role of the Ω-Loop in Regulation of the Catalytic Activity of TEM-Type β-Lactamases

Bacterial resistance to β-lactams, the most commonly used class of antibiotics, poses a global challenge. This resistance is caused by the production of bacterial enzymes that are termed β-lactamases (βLs). The evolution of serine-class A β-lactamases from penicillin-binding proteins (PBPs) is related to the formation of the Ω-loop at the entrance to the enzyme’s active site. In this loop, the Glu166 residue plays a key role in the two-step catalytic cycle of hydrolysis. This residue in TEM–type β-lactamases, together with Asn170, is involved in the formation of a hydrogen bonding network with a water molecule, leading to the deacylation of the acyl–enzyme complex and the hydrolysis of the β-lactam ring of the antibiotic. The activity exhibited by the Ω-loop is attributed to the positioning of its N-terminal residues near the catalytically important residues of the active site. The structure of the Ω-loop of TEM-type β-lactamases is characterized by low mutability, a stable topology, and structural flexibility. All of the revealed features of the Ω-loop, as well as the mechanisms related to its involvement in catalysis, make it a potential target for novel allosteric inhibitors of β-lactamases.

Machine Learning Classification Model for Functional Binding Modes of TEM-1 β-Lactamase

TEM family of enzymes is one of the most commonly encountered β-lactamases groups with different catalytic capabilities against various antibiotics. Despite the studies investigating the catalytic mechanism of TEM β-lactamases, the binding modes of these enzymes against ligands in different functional catalytic states have been largely overlooked. But the binding modes may play a critical role in the function and even the evolution of these proteins. In this work, a newly developed machine learning analysis approach to the recognition of protein dynamics states was applied to compare the binding modes of TEM-1 β-lactamase with regard to penicillin in different catalytic states. While conventional analysis methods, including principal components analysis (PCA), could not differentiate TEM-1 in different binding modes, the application of a machine learning method led to excellent classification models differentiating these states. It was also revealed that both reactant/product states and apo/product states are more differentiable than the apo/reactant states. The feature importance generated by the training procedure of the machine learning model was utilized to evaluate the contribution from residues at active sites and in different secondary structures. Key active site residues, Ser70 and Ser130, play a critical role in differentiating reactant/product states, while other active site residues are more important for differentiating apo/product states. Overall, this study provides new insights into the different dynamical function states of TEM-1 and may open a new venue for β-lactamases functional and evolutional studies in general.

Influence of substrates and inhibitors on the structure of Klebsiella pneumoniae carbapenemase-2

The hydrolysis of last resort carbapenem antibiotics by Klebsiella pneumoniae carbapenemase-2 (KPC-2) presents a significant danger to global health. Combined with horizontal gene transfer, the emergence KPC-2 threatens to quickly expand carbapenemase activity to ever increasing numbers of pathogens. Our understanding of KPC-2 has greatly increased over the past decade thanks, in great part, to 20 crystal structures solved by groups around the world. These include apo KPC-2 structures, along with structures featuring a library of 10 different inhibitors representing diverse structural and functional classes. Herein we focus on cataloging the available KPC-2 structures and presenting a discussion of key aspects of each structure and important relationships between structures. Although the available structures do not provide information on dynamic motions with KPC-2, and the family of structures indicates small conformational changes across a wide array of bound inhibitors, substrates, and products, the structures provide a strong foundation for additional studies in the coming years to discover new KPC-2 inhibitors.

Impact statement

The work herein is important to the field as it provides a clear and succinct accounting of available KPC-2 structures. The work advances the field by collecting and analyzing differences and similarities across the available structures. This work features new analyses and interpretations of the existing structures which will impact the field in a positive way by making structural insights more widely available among the beta-lactamase community.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

[Investigation the role of mutations M182T and Q39K in structure of beta-lactamase TEM-72 by molecular dynamics method]

Synthesis of b-lactamases is one of the common mechanisms of bacterial resistance to b-lactam antibiotics including penicillins and cephalosporins. The widespread use of antibiotics results in appearance of numerous extended-spectrum b-lactamase variants or resistance to inhibitors. Mutations of 92 residues of TEM type were found. Several mutations are the key mutations that determine the extension of spectrum of substrates. However, roles of the most associated mutations, located far from active site, remain unknown. We have investigated the role of associated mutations in structure of b-lactamase TEM-72, which contain two key mutation (G238S, E240K) and two associated mutations (Q39K, M182T) by means of simulation of molecular dynamics. The key mutation lead to destabilization of the protein globule, characterized by increased mobility of amino acid residues at high temperature of modelling. Mutation M182T lead to stabilization protein, whereas mutation Q39K is destabilizing mutation. It seems that the last mutation serves for optimization of conformational mobility of b-lactamase and may influence on enzyme activity.