The Molecular Mechanisms Underlying Hidden Phenotypic Variation among Metallo-β-Lactamases

Understanding epistatic networks in the B1 β-lactamases through coevolutionary statistical modeling and deep mutational scanning

Throughout evolution, protein families undergo substantial sequence divergence while preserving structure and function. Although most mutations are deleterious, evolution can explore sequence space via epistatic networks of intramolecular interactions that alleviate the harmful mutations. However, comprehensive analysis of such epistatic networks across protein families remains limited. Thus, we conduct a family wide analysis of the B1 metallo-β-lactamases, combining experiments (deep mutational scanning, DMS) on two distant homologs (NDM-1 and VIM-2) and computational analyses (in silico DMS based on Direct Coupling Analysis, DCA) of 100 homologs. The methods jointly reveal and quantify prevalent epistasis, as ~1/3rd of equivalent mutations are epistatic in DMS. From DCA, half of the positions have a >6.5 fold difference in effective number of tolerated mutations across the entire family. Notably, both methods locate residues with the strongest epistasis in regions of intermediate residue burial, suggesting a balance of residue packing and mutational freedom in forming epistatic networks. We identify entrenched WT residues between NDM-1 and VIM-2 in DMS, which display statistically distinct behaviors in DCA from non-entrenched residues. Entrenched residues are not easily compensated by changes in single nearby interactions, reinforcing existing findings where a complex epistatic network compounds smaller effects from many interacting residues.

β-lactamase expression induces collateral sensitivity in Escherichia coli

Major antibiotic groups are losing effectiveness due to the uncontrollable spread of antimicrobial resistance (AMR) genes. Among these, β-lactam resistance genes -encoding β-lactamases- stand as the most common resistance mechanism in Enterobacterales due to their frequent association with mobile genetic elements. In this context, novel approaches that counter mobile AMR are urgently needed. Collateral sensitivity (CS) occurs when the acquisition of resistance to one antibiotic increases susceptibility to another antibiotic and can be exploited to eliminate AMR selectively. However, most CS networks described so far emerge as a consequence of chromosomal mutations and cannot be leveraged to tackle mobile AMR. Here, we dissect the CS response elicited by the acquisition of a prevalent antibiotic resistance plasmid to reveal that the expression of the β-lactamase gene blaOXA-48 induces CS to colistin and azithromycin. We next show that other clinically relevant mobile β-lactamases produce similar CS responses in multiple, phylogenetically unrelated E. coli strains. Finally, by combining experiments with surveillance data comprising thousands of antibiotic susceptibility tests, we show that β-lactamase-induced CS is pervasive within Enterobacterales. These results highlight that the physiological side-effects of β-lactamases can be leveraged therapeutically, paving the way for the rational design of specific therapies to block mobile AMR or at least counteract their effects.

Identification and characterization of CIM-1, a carbapenemase that adds to the family of resistance factors against last resort antibiotics

The increasing rate of carbapenem-resistant bacteria within healthcare environments is an issue of great concern that needs urgent attention. This resistance is driven by metallo-β-lactamases (MBLs), which can catalyse the hydrolysis of almost all clinically available β-lactams and are resistant to all the clinically utilized β-lactamase inhibitors. In this study, an uncharacterized MBL is identified in a multidrug resistant isolate of the opportunistic pathogen, Chryseobacterium indologenes. Sequence analysis predicts this MBL (CIM-1) to be a lipoprotein with an atypical lipobox. Characterization of CIM-1 reveals it to be a high-affinity carbapenemase with a broad spectrum of activity that includes all cephalosporins and carbapenems. Results also shown that CIM-1 is potentially a membrane-associated MBL with an uncharacterized lipobox. Using prediction tools, we also identify more potentially lipidated MBLs with non-canonical lipoboxes highlighting the necessity of further investigation of lipidated MBLs.

Measuring differential fitness costs and interactions between genetic cassettes using fluorescent spectrophotometry

In this article, we present a method for designing, executing, and analyzing data from a microbial competition experiment. We use fluorescent reporters to label different competing strains and resolve individual growth curves using a fluorescent spectrophotometer. Our comprehensive data analysis pipeline integrates multiple experiments to simultaneously infer sources of variation, extract selection coefficients, and estimate the genetic contributions to fitness for various synthetic genetic cassettes (SGCs). To demonstrate the method, we employ a synthetic biological system based on Escherichia coli. Strains carry 1 of 10 different plasmids and one of three genomically integrated fluorescent markers. All strains are co-cultured to obtain real-time measurements of optical density (total population density) and fluorescence (sub-population densities). We identify challenges in calibrating between fluorescence and density and of fluorescent proteins maturing at different rates. To resolve these issues, we compare two methods of fluorescence calibration and correct for maturation by measuring in vivo maturation times. We provide evidence of genetic interactions occurring between our SGCs and further show how to use our statistical model to test some hypotheses about microbial growth and the costs of protein expression.IMPORTANCEFluorescently labeled co-cultures are becoming increasingly popular. The approach proposed here offers a high standard for experimental design and data analysis to measure selection coefficients and growth rates in competition. Measuring competitive differences is useful in many laboratory studies, allowing for fitness cost-correction of growth rates and ecological interactions and testing hypotheses in synthetic biology. Using time-resolved growth curves, rather than endpoint measurements, for competition assays allows us to construct a detailed scientific model that can be used to ask questions about fine-grained phenomena, such as bacterial growth dynamics, as well as higher-level phenomena, such as the interactions between synthetic cassette expression.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability

Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.

Deciphering the evolution of metallo-β-lactamases: a journey from the test tube to the bacterial periplasm

Understanding the evolution of metallo-β-lactamases (MBLs) is fundamental to deciphering the mechanistic basis of resistance to carbapenems in pathogenic and opportunistic bacteria. Presently, these MBL producing pathogens are linked to high rates of morbidity and mortality worldwide. However, the study of the biochemical and biophysical features of MBLs in vitro provides an incomplete picture of their evolutionary potential, since this limited and artificial environment disregards the physiological context where evolution and selection take place. Herein, we describe recent efforts aimed to address the evolutionary traits acquired by different clinical variants of MBLs in conditions mimicking their native environment (the bacterial periplasm) and considering whether they are soluble or membrane-bound proteins. This includes addressing the metal content of MBLs within the cell under zinc starvation conditions, and the context provided by different bacterial hosts that result in particular resistance phenotypes. Our analysis highlights recent progress bridging the gap between in vitro and in-cell studies.

Gene sharing among plasmids and chromosomes reveals barriers for antibiotic resistance gene transfer

The emergence of antibiotic resistant bacteria is a major threat to modern medicine. Rapid adaptation to antibiotics is often mediated by the acquisition of plasmids carrying antibiotic resistance (ABR) genes. Nonetheless, the determinants of plasmid-mediated ABR gene transfer remain debated. Here, we show that the propensity of ABR gene transfer via plasmids is higher for accessory chromosomal ABR genes in comparison with core chromosomal ABR genes, regardless of the resistance mechanism. Analysing the pattern of ABR gene occurrence in the genomes of 2635 Enterobacteriaceae isolates, we find that 33% of the 416 ABR genes are shared between chromosomes and plasmids. Phylogenetic reconstruction of ABR genes occurring on both plasmids and chromosomes supports their evolution by lateral gene transfer. Furthermore, accessory ABR genes (encoded in less than 10% of the chromosomes) occur more abundantly in plasmids in comparison with core ABR genes (encoded in greater than or equal to 90% of the chromosomes). The pattern of ABR gene occurrence in plasmids and chromosomes is similar to that in the total Escherichia genome. Our results thus indicate that the previously recognized barriers for gene acquisition by lateral gene transfer apply also to ABR genes. We propose that the functional complexity of the underlying ABR mechanism is an important determinant of ABR gene transferability. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Analyses of HIV proteases variants at the threshold of viability reveals relationships between processing efficiency and fitness

Investigating the relationships between protein function and fitness provides keys for understanding biochemical mechanisms that underly evolution. Mutations with partial fitness defects can delineate the threshold of biochemical function required for viability. We utilized a previous deep mutational scan of HIV-1 protease (PR) to identify variants with 15–45 per cent defects in replication and analysed the biochemical function of eight variants (L10M, L10S, V32C, V32I, A71V, A71S, Q92I, Q92N). We purified each variant and assessed the efficiency of peptide cleavage for three cut sites (MA-CA, TF-PR, and PR-RT) as well as gel-based analyses of processing of purified Gag. The cutting activity of at least one site was perturbed relative to WT protease for all variants, consistent with cutting activity being a primary determinant of fitness effects. We examined the correlation of fitness defects with cutting activity of different sites. MA-CA showed the weakest correlation (R² = 0.02) with fitness, suggesting relatively weak coupling with viral replication. In contrast, cutting of the TF-PR site showed the strongest correlation with fitness (R² = 0.53). Cutting at the TF-PR site creates a new PR protein with a free N-terminus that is critical for activity. Our findings indicate that increasing the pool of active PR is rate limiting for viral replication, making this an ideal step to target with inhibitors.

Heterologous production of cellulose- and starch-degrading hydrolases to expand Saccharomyces cerevisiae substrate utilization: Lessons learnt

Selected strains of Saccharomyces cerevisiae are used for commercial bioethanol production from cellulose and starch, but the high cost of exogenous enzymes for substrate hydrolysis remains a challenge. This can be addressed through consolidated bioprocessing (CBP) where S. cerevisiae strains are engineered to express recombinant glycoside hydrolases during fermentation. Looking back at numerous strategies undertaken over the past four decades to improve recombinant protein production in S. cerevisiae, it is evident that various steps in the protein production "pipeline" can be manipulated depending on the protein of interest and its anticipated application. In this review, we briefly introduce some of the strategies and highlight lessons learned with regards to improved transcription, translation, post-translational modification and protein secretion of heterologous hydrolases. We examine how host strain selection and modification, as well as enzyme compatibility, are crucial determinants for overall success. Finally, we discuss how lessons from heterologous hydrolase expression can inform modern synthetic biology and genome editing tools to provide process-ready yeast strains in future. However, it is clear that the successful expression of any particular enzyme is still unpredictable and requires a trial-and-error approach.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Structural and biochemical characterization of the environmental MBLs MYO-1, ECV-1 and SHD-1

BACKGROUND

MBLs form a large and heterogeneous group of bacterial enzymes conferring resistance to β-lactam antibiotics, including carbapenems. A large environmental reservoir of MBLs has been identified, which can act as a source for transfer into human pathogens. Therefore, structural investigation of environmental and clinically rare MBLs can give new insights into structure-activity relationships to explore the role of catalytic and second shell residues, which are under selective pressure.

OBJECTIVES

To investigate the structure and activity of the environmental subclass B1 MBLs MYO-1, SHD-1 and ECV-1.

METHODS

The respective genes of these MBLs were cloned into vectors and expressed in Escherichia coli. Purified enzymes were characterized with respect to their catalytic efficiency (kcat/Km). The enzymatic activities and MICs were determined for a panel of different β-lactams, including penicillins, cephalosporins and carbapenems. Thermostability was measured and structures were solved using X-ray crystallography (MYO-1 and ECV-1) or generated by homology modelling (SHD-1).

RESULTS

Expression of the environmental MBLs in E. coli resulted in the characteristic MBL profile, not affecting aztreonam susceptibility and decreasing susceptibility to carbapenems, cephalosporins and penicillins. The purified enzymes showed variable catalytic activity in the order of <5% to ∼70% compared with the clinically widespread NDM-1. The thermostability of ECV-1 and SHD-1 was up to 8°C higher than that of MYO-1 and NDM-1. Using solved structures and molecular modelling, we identified differences in their second shell composition, possibly responsible for their relatively low hydrolytic activity.

CONCLUSIONS

These results show the importance of environmental species acting as reservoirs for MBL-encoding genes.

Avoidance of protein unfolding constrains protein stability in long-term evolution

Every amino acid residue can influence a protein's overall stability, making stability highly susceptible to change throughout evolution. We consider the distribution of protein stabilities evolutionarily permittable under two previously reported protein fitness functions: flux dynamics and misfolding avoidance. We develop an evolutionary dynamics theory and find that it agrees better with an extensive protein stability data set for dihydrofolate reductase orthologs under the misfolding avoidance fitness function rather than the flux dynamics fitness function. Further investigation with ribonuclease H data demonstrates that not any misfolded state is avoided; rather, it is only the unfolded state. At the end, we discuss how our work pertains to the universal protein abundance-evolutionary rate correlation seen across organisms' proteomes. We derive a closed-form expression relating protein abundance to evolutionary rate that captures Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens experimental trends without fitted parameters.

Cryptic β-Lactamase Evolution Is Driven by Low β-Lactam Concentrations

Very low antibiotic concentrations have been shown to drive the evolution of antimicrobial resistance. While substantial progress has been made to understand the driving role of low concentrations during resistance development for different antimicrobial classes, the importance of β-lactams, the most commonly used antibiotics, is still poorly studied.

Our current understanding of how low antibiotic concentrations shape the evolution of contemporary β-lactamases is limited. Using the widespread carbapenemase OXA-48, we tested the long-standing hypothesis that selective compartments with low antibiotic concentrations cause standing genetic diversity that could act as a gateway to developing clinical resistance. Here, we subjected Escherichia coli expressing bla_OXA-48, on a clinical plasmid, to experimental evolution at sub-MICs of ceftazidime. We identified and characterized seven single variants of OXA-48. Susceptibility profiles and dose-response curves showed that they increased resistance only marginally. However, in competition experiments at sub-MICs of ceftazidime, they demonstrated strong selectable fitness benefits. Increased resistance was also reflected in elevated catalytic efficiencies toward ceftazidime. These changes are likely caused by enhanced flexibility of the Ω- and β5-β6 loops and fine-tuning of preexisting active site residues. In conclusion, low-level concentrations of β-lactams can drive the evolution of β-lactamases through cryptic phenotypes which may act as stepping-stones toward clinical resistance.

IMPORTANCE Very low antibiotic concentrations have been shown to drive the evolution of antimicrobial resistance. While substantial progress has been made to understand the driving role of low concentrations during resistance development for different antimicrobial classes, the importance of β-lactams, the most commonly used antibiotics, is still poorly studied. Here, we shed light on the evolutionary impact of low β-lactam concentrations on the widespread β-lactamase OXA-48. Our data indicate that the exposure to β-lactams at very low concentrations enhances β-lactamase diversity and drives the evolution of β-lactamases by significantly influencing their substrate specificity. Thus, in contrast to high concentrations, low levels of these drugs may substantially contribute to the diversification and divergent evolution of these enzymes, providing a standing genetic diversity that can be selected and mobilized when antibiotic pressure increases.

Antibiotic resistance: turning evolutionary principles into clinical reality

Antibiotic resistance is one of the major challenges facing modern medicine worldwide. The past few decades have witnessed rapid progress in our understanding of the multiple factors that affect the emergence and spread of antibiotic resistance at the population level and the level of the individual patient. However, the process of translating this progress into health policy and clinical practice has been slow. Here, we attempt to consolidate current knowledge about the evolution and ecology of antibiotic resistance into a roadmap for future research as well as clinical and environmental control of antibiotic resistance. At the population level, we examine emergence, transmission and dissemination of antibiotic resistance, and at the patient level, we examine adaptation involving bacterial physiology and host resilience. Finally, we describe new approaches and technologies for improving diagnosis and treatment and minimizing the spread of resistance.

Comprehensive exploration of the translocation, stability and substrate recognition requirements in VIM-2 lactamase

Metallo-β-lactamases (MBLs) degrade a broad spectrum of β-lactam antibiotics, and are a major disseminating source for multidrug resistant bacteria. Despite many biochemical studies in diverse MBLs, molecular understanding of the roles of residues in the enzyme's stability and function, and especially substrate specificity, is lacking. Here, we employ deep mutational scanning (DMS) to generate comprehensive single amino acid variant data on a major clinical MBL, VIM-2, by measuring the effect of thousands of VIM-2 mutants on the degradation of three representative classes of β-lactams (ampicillin, cefotaxime, and meropenem) and at two different temperatures (25°C and 37°C). We revealed residues responsible for expression and translocation, and mutations that increase resistance and/or alter substrate specificity. The distribution of specificity-altering mutations unveiled distinct molecular recognition of the three substrates. Moreover, these function-altering mutations are frequently observed among naturally occurring variants, suggesting that the enzymes have continuously evolved to become more potent resistance genes.

Protein determinants of dissemination and host specificity of metallo-β-lactamases

The worldwide dissemination of metallo-β-lactamases (MBLs), mediating resistance to carbapenem antibiotics, is a major public health problem. The extent of dissemination of MBLs such as VIM-2, SPM-1 and NDM among Gram-negative pathogens cannot be explained solely based on the associated mobile genetic elements or the resistance phenotype. Here, we report that MBL host range is determined by the impact of MBL expression on bacterial fitness. The signal peptide sequence of MBLs dictates their adaptability to each host. In uncommon hosts, inefficient processing of MBLs leads to accumulation of toxic intermediates that compromises bacterial growth. This fitness cost explains the exclusion of VIM-2 and SPM-1 from Escherichia coli and Acinetobacter baumannii, and their confinement to Pseudomonas aeruginosa. By contrast, NDMs are expressed without any apparent fitness cost in different bacteria, and are secreted into outer membrane vesicles. We propose that the successful dissemination and adaptation of MBLs to different bacterial hosts depend on protein determinants that enable host adaptability and carbapenem resistance.

Metallo-β-lactamases (MBLs) confer resistance to carbapenem antibiotics. Here, López et al. show that the host range of MBLs depends on the efficiency of MBL signal peptide processing and secretion into outer membrane vesicles, which affects bacterial fitness.