NMR dynamics of PSE-4 beta-lactamase: an interplay of ps-ns order and mus-ms motions in the active site

Slow Protein Dynamics Elicits New Enzymatic Functions by Means of Epistatic Interactions

Protein evolution depends on the adaptation of these molecules to different functional challenges. This occurs by tuning their biochemical, biophysical, and structural traits through the accumulation of mutations. While the role of protein dynamics in biochemistry is well recognized, there are limited examples providing experimental evidence of the optimization of protein dynamics during evolution. Here we report an NMR study of four variants of the CTX-M β-lactamases, in which the interplay of two mutations outside the active site enhances the activity against a cephalosporin substrate, ceftazidime. The crystal structures of these enzymes do not account for this activity enhancement. By using NMR, here we show that the combination of these two mutations increases the backbone dynamics in a slow timescale and the exposure to the solvent of an otherwise buried β-sheet. The two mutations located in this β-sheet trigger conformational changes in loops located at the opposite side of the active site. We postulate that the most active variant explores alternative conformations that enable binding of the more challenging substrate ceftazidime. The impact of the mutations in the dynamics is context-dependent, in line with the epistatic effect observed in the catalytic activity of the different variants. These results reveal the existence of a dynamic network in CTX-M β-lactamases that has been exploited in evolution to provide a net gain-of-function, highlighting the role of alternative conformations in protein evolution.

Structural Comparisons of Cefotaximase (CTX-M-ase) Sub Family 1

The cefotaximase or CTX-M, family of serine-β-lactamases represents a significant clinical concern due to the ability for these enzymes to confer resistance to a broad array of β-lactam antibiotics an inhibitors. This behavior lends CTX-M-ases to be classified as extended spectrum β-lactamases (ESBL). Across the family of CTX-M-ases most closely related to CTX-M-1, the structures of CTX-M-15 with a library of different ligands have been solved and serve as the basis of comparison within this review. Herein we focus on the structural changes apparent in structures of CTX-M-15 in complex with diazabicyclooctane (DABCO) and boronic acid transition state analog inhibitors. Interactions between a positive surface patch near the active site and complementary functional groups of the bound inhibitor play key roles in the dictating the conformations of active site residues. The insights provided by analyzing structures of CTX-M-15 in complex with DABCO and boronic acid transition state analog inhibitors and analyzing existing structures of CTX-M-64 offer opportunities to move closer to making predictions as to how CTX-M-ases may interact with potential drug candidates, setting the stage for the further development of new antibiotics and β-lactamase inhibitors.

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.

The Structural Dynamics of Engineered β-Lactamases Vary Broadly on Three Timescales yet Sustain Native Function

Understanding the principles of protein dynamics will help guide engineering of protein function: altering protein motions may be a barrier to success or may be an enabling tool for protein engineering. The impact of dynamics on protein function is typically reported over a fraction of the full scope of motional timescales. If motional patterns vary significantly at different timescales, then only by monitoring motions broadly will we understand the impact of protein dynamics on engineering functional proteins. Using an integrative approach combining experimental and in silico methodologies, we elucidate protein dynamics over the entire span of fast to slow timescales (ps to ms) for a laboratory-engineered system composed of five interrelated β-lactamases: two natural homologs and three laboratory-recombined variants. Fast (ps-ns) and intermediate (ns-µs) dynamics were mostly conserved. However, slow motions (µs-ms) were few and conserved in the natural homologs yet were numerous and widely dispersed in their recombinants. Nonetheless, modified slow dynamics were functionally tolerated. Crystallographic B-factors from high-resolution X-ray structures were partly predictive of the conserved motions but not of the new slow motions captured in our solution studies. Our inspection of protein dynamics over a continuous range of timescales vividly illustrates the complexity of dynamic impacts of protein engineering as well as the functional tolerance of an engineered enzyme system to new slow motions.

Defining the architecture of KPC-2 Carbapenemase: identifying allosteric networks to fight antibiotics resistance

The rise of multi-drug resistance in bacterial pathogens is one of the grand challenges facing medical science. A major concern is the speed of development of β-lactamase-mediated resistance in Gram-negative species, thus putting at risk the efficacy of the most recently approved antibiotics and inhibitors, including carbapenems and avibactam, respectively. New strategies to overcome resistance are urgently required, which will ultimately be facilitated by a deeper understanding of the mechanisms that regulate the function of β-lactamases such as the Klebsiella Pneumoniae carbapenemases (KPCs). Using enhanced sampling computational methods together with site-directed mutagenesis, we report the identification of two “hydrophobic networks” in the KPC-2 enzyme, the integrity of which has been found to be essential for protein stability and corresponding resistance. Present throughout the structure, these networks are responsible for the structural integrity and allosteric signaling. Disruption of the networks leads to a loss of the KPC-2 mediated resistance phenotype, resulting in restored susceptibility to different classes of β-lactam antibiotics including carbapenems and cephalosporins. The ”hydrophobic networks” were found to be highly conserved among class-A β-lactamases, which implies their suitability for exploitation as a potential target for therapeutic intervention.

A Distal Disulfide Bridge in OXA-1 β-Lactamase Stabilizes the Catalytic Center and Alters the Dynamics of the Specificity Determining Ω Loop

Widespread antibiotic resistance, particularly when mediated by broad-spectrum β-lactamases, has major implications for public health. Substitutions in the active site often allow broad-spectrum enzymes to accommodate diverse types of β-lactams. Substitutions observed outside the active site are thought to compensate for the loss of thermal stability. The OXA-1 clade of class D β-lactamases contains a pair of conserved cysteines located outside the active site that forms a disulfide bond in the periplasm. Here, the effect of the distal disulfide bond on the structure and dynamics of OXA-1 was investigated via 4 μs molecular dynamics simulations. The results reveal that the disulfide promotes the preorganized orientation of the catalytic residues and affects the conformation of the functionally important Ω loop. Furthermore, principal component analysis reveals differences in the global dynamics between the oxidized and reduced forms, especially in the motions involving the Ω loop. A dynamical network analysis indicates that, in the oxidized form, in addition to its role in ligand binding, the KTG family motif is a central hub of the global dynamics. As activity of OXA-1 has been measured only in the reduced form, we suggest that accurate assessment of its functional profile would require oxidative conditions mimicking periplasm.

15N, 13C and 1H backbone resonance assignments of an artificially engineered TEM-1/PSE-4 class A β-lactamase chimera and its deconvoluted mutant

The widespread use of β-lactam antibiotics has given rise to a dramatic increase in clinically-relevant β-lactamases. Understanding the structure/function relation in these variants is essential to better address the ever-growing incidence of antibiotic resistance. We previously reported the backbone resonance assignments of a chimeric protein constituted of segments of the class A β-lactamases TEM-1 and PSE-4 (Morin et al. in Biomol NMR Assign 4:127-130, 2010. doi: 10.1007/s12104-010-9227-8 ). That chimera, cTEM17m, held 17 amino acid substitutions relative to TEM-1 β-lactamase, resulting in a well-folded and fully functional protein with increased dynamics. Here we report the (1)H, (13)C and (15)N backbone resonance assignments of chimera cTEM-19m, which includes 19 substitutions and exhibits increased active-site perturbation, as well as one of its deconvoluted variants, as the first step in the analysis of their dynamic behaviours.

Maintenance of native-like protein dynamics may not be required for engineering functional proteins

Proteins are dynamic systems, and understanding dynamics is critical for fully understanding protein function. Therefore, the question of whether laboratory engineering has an impact on protein dynamics is of general interest. Here, we demonstrate that two homologous, naturally evolved enzymes with high degrees of structural and functional conservation also exhibit conserved dynamics. Their similar set of slow timescale dynamics is highly restricted, consistent with evolutionary conservation of a functionally important feature. However, we also show that dynamics of a laboratory-engineered chimeric enzyme obtained by recombination of the two homologs exhibits striking difference on the millisecond timescale, despite function and high-resolution crystal structure (1.05 Å) being conserved. The laboratory-engineered chimera is thus functionally tolerant to modified dynamics on the timescale of catalytic turnover. Tolerance to dynamic variation implies that maintenance of native-like protein dynamics may not be required when engineering functional proteins.

Chimeric Avidin--NMR structure and dynamics of a 56 kDa homotetrameric thermostable protein

Chimeric avidin (ChiAVD) is a product of rational protein engineering remarkably resistant to heat and harsh conditions. In quest of the fundamentals behind factors affecting stability we have elucidated the solution NMR spectroscopic structure of the biotin–bound form of ChiAVD and characterized the protein dynamics through ¹⁵N relaxation and hydrogen/deuterium (H/D) exchange of this and the biotin–free form. To surmount the challenges arising from the very large size of the protein for NMR spectroscopy, we took advantage of its high thermostability. Conventional triple resonance experiments for fully protonated proteins combined with methyl–detection optimized experiments acquired at 58°C were adequate for the structure determination of this 56 kDa protein. The model–free parameters derived from the ¹⁵N relaxation data reveal a remarkably rigid protein at 58°C in both the biotin–bound and the free forms. The H/D exchange experiments indicate a notable increase in hydrogen protection upon biotin binding.

Chimeric β-lactamases: global conservation of parental function and fast time-scale dynamics with increased slow motions

Enzyme engineering has been facilitated by recombination of close homologues, followed by functional screening. In one such effort, chimeras of two class-A β-lactamases – TEM-1 and PSE-4 – were created according to structure-guided protein recombination and selected for their capacity to promote bacterial proliferation in the presence of ampicillin (Voigt et al., Nat. Struct. Biol. 2002 9:553). To provide a more detailed assessment of the effects of protein recombination on the structure and function of the resulting chimeric enzymes, we characterized a series of functional TEM-1/PSE-4 chimeras possessing between 17 and 92 substitutions relative to TEM-1 β-lactamase. Circular dichroism and thermal scanning fluorimetry revealed that the chimeras were generally well folded. Despite harbouring important sequence variation relative to either of the two ‘parental’ β-lactamases, the chimeric β-lactamases displayed substrate recognition spectra and reactivity similar to their most closely-related parent. To gain further insight into the changes induced by chimerization, the chimera with 17 substitutions was investigated by NMR spin relaxation. While high order was conserved on the ps-ns timescale, a hallmark of class A β-lactamases, evidence of additional slow motions on the µs-ms timescale was extracted from model-free calculations. This is consistent with the greater number of resonances that could not be assigned in this chimera relative to the parental β-lactamases, and is consistent with this well-folded and functional chimeric β-lactamase displaying increased slow time-scale motions.

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

The effects of substrate binding on class A β-lactamase dynamics were studied using molecular dynamics simulations of two model enzymes; 40 100-ns trajectories of the free and substrate-bound forms of TEM-1 (with benzylpenicillin) and PSE-4 (with carbenicillin) were recorded (totaling 4.0 μs). Substrates were parameterized with the CHARMM General Force Field. In both enzymes, the Ω loop exhibits a marked flexibility increase upon substrate binding, supporting the hypothesis of substrate gating. However, specific interactions that are formed or broken in the Ω loop upon binding differ between the two enzymes: dynamics are conserved, but not specific interactions. Substrate binding also has a global structuring effect on TEM-1, but not on PSE-4. Changes in TEM-1's normal modes show long-range effects of substrate binding on enzyme dynamics. Hydrogen bonds observed in the active site are mostly preserved upon substrate binding, and new, transient interactions are also formed. Agreement between NMR relaxation parameters and our theoretical results highlights the dynamic duality of class A β-lactamases: enzymes that are highly structured on the ps-ns timescale, with important flexibility on the μs-ms timescale in regions such as the Ω loop.

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe²⁺-O₂ + (•)NO → trHbN-Fe³⁺-OH₂ + NO₃⁻) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (•)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (•)NO and O₂ to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.

relaxGUI: a new software for fast and simple NMR relaxation data analysis and calculation of ps-ns and μs motion of proteins

Investigation of protein dynamics on the ps-ns and μs-ms timeframes provides detailed insight into the mechanisms of enzymes and the binding properties of proteins. Nuclear magnetic resonance (NMR) is an excellent tool for studying protein dynamics at atomic resolution. Analysis of relaxation data using model-free analysis can be a tedious and time consuming process, which requires good knowledge of scripting procedures. The software relaxGUI was developed for fast and simple model-free analysis and is fully integrated into the software package relax. It is written in Python and uses wxPython to build the graphical user interface (GUI) for maximum performance and multi-platform use. This software allows the analysis of NMR relaxation data with ease and the generation of publication quality graphs as well as color coded images of molecular structures. The interface is designed for simple data analysis and management. The software was tested and validated against the command line version of relax.

Backbone resonance assignments of an artificially engineered TEM-1/PSE-4 Class A β-lactamase chimera

The rapid evolution of Class A β-lactamases, which procure resistance to an increasingly broad panel of β-lactam antibiotics, underscores the urgency to better understand the relation between their sequence variation and their structural and functional features. To date, more than 300 clinically-relevant β-lactamase variants have been reported, and this number continues to increase. With the aim of obtaining insights into the evolutionary potential of β-lactamases, an artificially engineered, catalytically active chimera of the Class A TEM-1 and PSE-4 β-lactamases is under study by kinetics and NMR. Here we report the (1)H, (13)C and (15)N backbone resonance assignments for the 30 kDa chimera cTEM-17m. Despite its high molecular weight, the data provide evidence that this artificially-evolved chimeric enzyme is well folded. The hydrolytic activity of cTEM-17m was determined using the chromogenic substrate CENTA, with K (M) = 160 ± 35 μM and k (cat) = 20 ± 4 s(-1), which is in the same range as the values for TEM-1 and PSE-4 β-lactamases.

Interpretation of biomolecular NMR spin relaxation parameters

Biomolecular nuclear magnetic resonance (NMR) spin relaxation experiments provide exquisite information on the picosecond to nanosecond timescale motions of bond vectors. Spin-lattice (T1) and spin-spin (T2) relaxation times and the steady-state nuclear Overhauser effect (NOE) are the first set of parameters extracted from typical 15N or 13C NMR relaxation experiments. Therefore, verifying that T1, T2, and NOE are consistent with theoretical predictions is an important step before carrying out the more detailed model-free and reduced spectral density mapping analyses commonly employed. In this mini-review, we discuss the essential motional parameters used to describe biomolecular dynamics in the context of a variety of examples of folded and intrinsically disordered proteins and peptides in aqueous and membrane mimetic environments. Estimates of these parameters can be used as input for an online interface, introduced herein, allowing plotting of trends of T1, T2, and NOE with magnetic field strength. The plots may serve as a first-check to the spectroscopist preparing to embark on a detailed NMR relaxation analysis.

TEM-1 backbone dynamics-insights from combined molecular dynamics and nuclear magnetic resonance

Dynamic properties of class A beta-lactamase TEM-1 are investigated from molecular dynamics (MD) simulations. Comparison of MD-derived order parameters with those obtained from model-free analysis of nuclear magnetic resonance (NMR) relaxation data shows high agreement for N-H moieties within alpha- and beta-secondary structures, but significant deviation for those in loops. This was expected, because motions slower than the protein global tumbling often take place in loop regions. As previously shown using NMR, TEM-1 is a highly ordered protein. Motions are observed within the Omega loop that could, upon substrate binding, stabilize E166 in a catalytically efficient position as the cavity between the protein core and the Omega loop is partially filled. The rigidity of active site residues is consistent with the enzyme high turnover number. MD data are also shown to be useful during the model selection step of model-free analysis: local N-H motions observed over the course of the trajectories help assess whether a peptide plan undergoes low or high amplitude motions on one or more timescales. This joint use of MD and NMR provides a better description of protein dynamics than would be possible using either technique alone.

Simple tests for the validation of multiple field spin relaxation data

(15)N spin relaxation data is widely used to extract detailed dynamic information regarding bond vectors such as the amide N-H bond of the protein backbone. Analysis is typically carried using the Lipari-Szabo model-free approach. Even though the original model-free equation can be determined from single field R (1), R (2) and NOE, over-determination of more complex motional models is dependent on the recording of multiple field datasets. This is especially important for the characterization of conformational exchange which affects R (2) in a field dependent manner. However, severe artifacts can be introduced if inconsistencies arise between experimental setups with different magnets (or samples). Here, we propose the use of simple tests as validation tools for the assessment of consistency between different datasets recorded at multiple magnetic fields. Synthetic data are used to show the effects of inconsistencies on the proposed tests. Moreover, an analysis of data currently deposited in the BMRB is performed. Finally, two cases from our laboratory are presented. These tests are implemented in the open-source program relax, and we propose their use as a routine check-up for assessment of multiple field dataset consistency prior to any analysis such as model-free calculations. We believe this will aid in the extraction of higher quality dynamics information from (15)N spin relaxation data.