Gating interactions steer loop conformational changes in the active site of the L1 metallo-β-lactamase

  1. Zhuoran Zhao
  2. Xiayu Shen
  3. Shuang Chen
  4. Jing Gu
  5. Haun Wang
  6. Maria F Mojica
  7. Moumita Samanta
  8. Debsindhu Bhowmik
  9. Alejandro J Vila
  10. Robert A Bonomo
  11. Shozeb Haider  Is a corresponding author
  1. Department of Pharmaceutical and Biological Chemistry, School of Pharmacy, University College London, United Kingdom
  2. Department of Molecular Biology and Microbiology, Case Western Reserve University School of Medicine, United States
  3. Louis Stokes Cleveland Department of Veterans Affairs Medical Center, United States
  4. CWRU-Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES), United States
  5. College of Computing, Georgia Institute of Technology, United States
  6. Computer Science and Engineering Division, Oak Ridge National Laboratories, United States
  7. Laboratorio de Metaloproteínas, Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), Argentina
  8. Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Argentina
  9. Departments of Medicine, Biochemistry, Pharmacology, and Proteomics and Bioinformatics, Case Western Reserve University School of Medicine, United States
  10. UCL Centre for Advanced Research Computing, University College London, United Kingdom
12 figures, 2 tables and 1 additional file

Figures

Figure 1 with 2 supplements
Overall structural features of the L1 metallo-β-lactamase (MBL).

(A) The homo-tetrameric conformation of L1 MBL (PDB entry 1SML). Three monomers are illustrated in surface representation; the fourth monomer is shown in green cartoon. (B) A close up of the …

Figure 1—figure supplement 1
Catalytic cycle of L1 metallo-β-lactamase (MBL).

Proposed mechanism of β-lactam hydrolysis by S. maltophilia L1 MBL.

Figure 1—figure supplement 2
Sequence alignment of seven class B metallo-β-lactamase with the BBL numbering.

(A) Conserved secondary structure elements for the three subclasses are indicated above the sequences: α-helix, β-sheets, and L-loop. Amino acid insertions in the newly sequences enzymes are …

Figure 2 with 1 supplement
L1 metallo-β-lactamase Markov State Model.

(A) Implied timescales (ITS) plot. The ITS plot describes the convergence behavior of the implied timescales associated with the processes at the lag time of 5 ns. (B) Net flux pathways plot …

Figure 2—figure supplement 1
Chapman–Kolmogorov (CK) test plot.

The CK test plot, along with the ITS plot (Figure 2A) was used to assess the Markovian behavior of the constructed model.

Figure 3 with 6 supplements
The L1 metallo-β-lactamase active site remodeling.

Structures exist in the closed, intermediate, and open states defined by a salt bridge between D150c-R236, π–π stacking interaction between H151-Y227 and the conformation of the cyclic side chain of …

Figure 3—figure supplement 1
Probability weighted distance distributions in the open, intermediate, and the closed states.

The distances between the D150c-R236 salt bridge (blue), H151-Y227 π–π stacking (orange), A224-Q310 (red), and P225-Zn (green) have been illustrated.

Figure 3—figure supplement 2
Salt bridge interaction between D150c-R236.

Violin plots highlighting the calculated D150c-R236 salt bridge interaction in each of the four subunits in different metastable states.

Figure 3—figure supplement 3
π–π stacking interaction between H151-Y227.

Violin plots highlighting the calculated H151-Y227 π–π stacking interaction in each of the four subunits in different metastable states.

Figure 3—figure supplement 4
Distance between P225 and the centroid of zinc ions.

Violin plots highlighting the calculated P225-Zn distance in each of the four subunits in different metastable states.

Figure 3—figure supplement 5
Interaction between A224-Q310.

Violin plots highlighting the calculated A224(bb)-Q310(sc) distance in each of the four subunits in different metastable states.

Figure 3—figure supplement 6
Structural comparison of the L1 metallo-β-lactamases.

A Cα-RMSD superimposition of state 7 (green) with a monomer from apo state crystal structure (PDB id 1SML; orange; 0.79 Å) and with inhibitor complex (PDB id 5DPX; purple; 0.77 Å) structure of L1 …

Figure 4 with 2 supplements
Convolutional variational autoencoder (CVAE)-based deep learning analysis.

CVAE-learned features of high-dimensional data represented in 2D following t-distributed stochastic neighbor embedding (tSNE) treatment. (A) The results show that A and C subunits show similar …

Figure 4—figure supplement 1
Convolutional variational autoencoder (CVAE)-based deep learning analysis.

(A) The training and validation loss plotted as assessed simultaneously over consecutive epochs at various reduced dimensions. Dimension 14 was selected for further analysis as it displayed the …

Figure 4—figure supplement 2
Clustering metastable states via t-distributed stochastic neighbor embedding (tSNE).

Individual subunits are plotted on the tSNE space after clustering of states 1–7. A comparison with Figure 4 highlights the open (state 1), intermediate open (state 4), intermediate closed (state …

Trans conformation of P225.

The cis/trans conformation of a peptide bond is determined by the torsion of the ω angle. The ω torsion angle is close to 0° in the cis configuration, or ca. 180° in the trans configuration. The ω …

Author response image 1
CVAE based deep learning analysis from a 28x28 matrix with 80 (training):20 (validation) data split.

(a) The training and validation loss plot as assessed over consecutive epochs at various reduced dimensions. (b) Validation loss during CVAE implementation is plotted at different latent dimensions …

Author response image 2
CVAE based deep learning analysis from a 28x28 matrix with 70 (training):30 (validation) data split.

(a) The training and validation loss plot as assessed over consecutive epochs at various reduced dimensions. (b) Validation loss during CVAE implementation is plotted at different latent dimensions …

Author response image 3
CVAE based deep learning analysis from a 17x17 matrix with 80 (training):20 (validation) data split.

(a) The training and validation loss plot as assessed over consecutive epochs at various reduced dimensions. (b) Validation loss during CVAE implementation is plotted at different latent dimensions …

Author response image 4
CVAE based deep learning analysis from a 17x17 matrix with 70 (training): 30 (validation) data split.

(a) The training and validation loss plot as assessed over consecutive epochs at various reduced dimensions. (b) Validation loss during CVAE implementation is plotted at different latent dimensions …

Author response image 5
Author response image 6
Author response image 7

Tables

Table 1
Elongated loops.
α3-β7GGSDDLHFGDGITYPPANAD Residues 149–171
β12-α5ADSLSAPGYQLQGNPRYPH Residues 219–239
Table 2
Susceptibility results of various β-lactams against L1 variants expressed in the pBC SK (+) vector in E. coli DH10B cells.

The reported values are the mode of at least three biological replicates.

StrainMIC value (mg/l)*
CAZFEPIMIMERTEBTEB + CS319
DH10B pBC SK blaL1a25621616322
DH10B pBC SK blaL1a_H151V20.030.1250.03≤0.03≤0.03
DH10B pBC SK blaL1a_D150cV25618840.25
DH10B pBC SK blaL1a_P225A640.1250.2520.125≤0.03
DH10B pBC SK blaL1a_R236A25628840.25
DH10B pBC SK blaL1a_Y227A640.250.2520.5≤0.03
  1. *

    CAZ, ceftazidime; FEP, cefepime; IMI, imipenem; MER, meropenem; TEB, tebipenem.

  2. CS319 was added to a final concentration of 100 mg/l.

Additional files

Download links