Optimization of Conformational Dynamics in an Epistatic Evolutionary Trajectory.

The understanding of protein evolution depends on the ability to relate the impact of mutations on molecular traits to organismal fitness. Biological activity and robustness have been regarded as important features in shaping protein evolutionary landscapes. Conformational dynamics, which is essential for protein function, has received little attention in the context of evolutionary analyses. Here we employ NMR spectroscopy, the chief experimental tool to describe protein dynamics at atomic level in solution at room temperature, to study the intrinsic dynamic features of a metallo- Β: -lactamase enzyme and three variants identified during a directed evolution experiment that led to an expanded substrate profile. We show that conformational dynamics in the catalytically relevant microsecond to millisecond timescale is optimized along the favored evolutionary trajectory. In addition, we observe that the effects of mutations on dynamics are epistatic. Mutation Gly262Ser introduces slow dynamics on several residues that surround the active site when introduced in the wild-type enzyme. Mutation Asn70Ser removes the slow dynamics observed for few residues of the wild-type enzyme, but increases the number of residues that undergo slow dynamics when introduced in the Gly262Ser mutant. These effects on dynamics correlate with the epistatic interaction between these two mutations on the bacterial phenotype. These findings indicate that conformational dynamics is an evolvable trait, and that proteins endowed with more dynamic active sites also display a larger potential for promoting evolution.

Quantitative Description of a Protein Fitness Landscape Based on Molecular Features.

Understanding the driving forces behind protein evolution requires the ability to correlate the molecular impact of mutations with organismal fitness. To address this issue, we employ here metallo-ß-lactamases as a model system, which are Zn(II) dependent enzymes that mediate antibiotic resistance. We present a study of all the possible evolutionary pathways leading to a metallo-ß-lactamase variant optimized by directed evolution. By studying the activity, stability and Zn(II) binding capabilities of all mutants in the preferred evolutionary pathways, we show that this local fitness landscape is strongly conditioned by epistatic interactions arising from the pleiotropic effect of mutations in the different molecular features of the enzyme. Activity and stability assays in purified enzymes do not provide explanatory power. Instead, measurement of these molecular features in an environment resembling the native one provides an accurate description of the observed antibiotic resistance profile. We report that optimization of Zn(II) binding abilities of metallo-ß-lactamases during evolution is more critical than stabilization of the protein to enhance fitness. A global analysis of these parameters allows us to connect genotype with fitness based on quantitative biochemical and biophysical parameters.

Metallo-ß-lactamases withstand low Zn(II) conditions by tuning metal-ligand interactions.

A number of multiresistant bacterial pathogens inactivate antibiotics by producing Zn(II)-dependent ß-lactamases. We show that metal uptake leading to an active dinuclear enzyme in the periplasmic space of Gram-negative bacteria is ensured by a cysteine residue, an unusual metal ligand in oxidizing environments. Kinetic, structural and affinity data show that such Zn(II)-cysteine interaction is an adaptive trait that tunes the metal binding affinity, thus enabling antibiotic resistance at restrictive Zn(II) concentrations.

Longitudinal Evolution of the Pseudomonas-Derived Cephalosporinase (PDC) Structure and Activity in a Cystic Fibrosis Patient Treated with ß-Lactams.

Traditional studies on the evolution of antibiotic resistance development use approaches that can range from laboratory-based experimental studies, to epidemiological surveillance, to sequencing of clinical isolates. However, evolutionary trajectories also depend on the environment in which selection takes place, compelling the need to more deeply investigate the impact of environmental complexities and their dynamics over time. Herein, we explored the within-patient adaptive long-term evolution of a Pseudomonas aeruginosa hypermutator lineage in the airways of a cystic fibrosis (CF) patient by performing a chronological tracking of mutations that occurred in different subpopulations; our results demonstrated parallel evolution events in the chromosomally encoded class C ß-lactamase (blaPDC). These multiple mutations within blaPDC shaped diverse coexisting alleles, whose frequency dynamics responded to the changing antibiotic selective pressures for more than 26 years of chronic infection. Importantly, the combination of the cumulative mutations in blaPDC provided structural and functional protein changes that resulted in a continuous enhancement of its catalytic efficiency and high level of cephalosporin resistance. This evolution was linked to the persistent treatment with ceftazidime, which we demonstrated selected for variants with robust catalytic activity against this expanded-spectrum cephalosporin. A "gain of function" of collateral resistance toward ceftolozane, a more recently introduced cephalosporin that was not prescribed to this patient, was also observed, and the biochemical basis of this cross-resistance phenomenon was elucidated. This work unveils the evolutionary trajectories paved by bacteria toward a multidrug-resistant phenotype, driven by decades of antibiotic treatment in the natural CF environmental setting. IMPORTANCE Antibiotics are becoming increasingly ineffective to treat bacterial infections. It has been consequently predicted that infectious diseases will become the biggest challenge to human health in the near future. Pseudomonas aeruginosa is considered a paradigm in antimicrobial resistance as it exploits intrinsic and acquired resistance mechanisms to resist virtually all antibiotics known. AmpC ß-lactamase is the main mechanism driving resistance in this notorious pathogen to ß-lactams, one of the most widely used classes of antibiotics for cystic fibrosis infections. Here, we focus on the ß-lactamase gene as a model resistance determinant and unveil the trajectory P. aeruginosa undertakes on the path toward a multidrug-resistant phenotype during the course of two and a half decades of chronic infection in the airways of a cystic fibrosis patient. Integrating genetic and biochemical studies in the natural environment where evolution occurs, we provide a unique perspective on this challenging landscape, addressing fundamental molecular mechanisms of resistance.

Adaptive protein evolution grants organismal fitness by improving catalysis and flexibility.

Protein evolution is crucial for organismal adaptation and fitness. This process takes place by shaping a given 3-dimensional fold for its particular biochemical function within the metabolic requirements and constraints of the environment. The complex interplay between sequence, structure, functionality, and stability that gives rise to a particular phenotype has limited the identification of traits acquired through evolution. This is further complicated by the fact that mutations are pleiotropic, and interactions between mutations are not always understood. Antibiotic resistance mediated by beta-lactamases represents an evolutionary paradigm in which organismal fitness depends on the catalytic efficiency of a single enzyme. Based on this, we have dissected the structural and mechanistic features acquired by an optimized metallo-beta-lactamase (MbetaL) obtained by directed evolution. We show that antibiotic resistance mediated by this enzyme is driven by 2 mutations with sign epistasis. One mutation stabilizes a catalytically relevant intermediate by fine tuning the position of 1 metal ion; whereas the other acts by augmenting the protein flexibility. We found that enzyme evolution (and the associated antibiotic resistance) occurred at the expense of the protein stability, revealing that MbetaLs have not exhausted their stability threshold. Our results demonstrate that flexibility is an essential trait that can be acquired during evolution on stable protein scaffolds. Directed evolution aided by a thorough characterization of the selected proteins can be successfully used to predict future evolutionary events and design inhibitors with an evolutionary perspective.

Engineered mononuclear variants in Bacillus cereus metallo-beta-lactamase BcII are inactive.

Metallo-beta-lactamases (MbetaLs) are zinc enzymes able to hydrolyze almost all beta-lactam antibiotics, rendering them inactive, at the same time endowing bacteria high levels of resistance. The design of inhibitors active against all classes of MbetaLs has been hampered by their structural diversity and by the heterogeneity in metal content in enzymes from different sources. BcII is the metallo-beta-lactamase from Bacillus cereus, which is found in both the mononuclear and dinuclear forms. Despite extensive studies, there is still controversy about the nature of the active BcII species. Here we have designed two mutant enzymes in which each one of the metal binding sites was selectively removed. Both mutants were almost inactive, despite preserving most of the structural features of each metal site. These results reveal that neither site isolated in the MbetaL scaffold is sufficient to render a fully active enzyme. This suggests that only the dinuclear species is active or that the mononuclear variants can be active only if aided by other residues that would be metal ligands in the dinuclear species.

Biochemical characterization of metallo-beta-lactamase VIM-11 from a Pseudomonas aeruginosa clinical strain.

A detailed biochemical characterization of the Pseudomonas aeruginosa VIM-11 metallo-beta-lactamase (MbetaL) is reported. The only substitution differentiating VIM-11 from VIM-2 (N165S) promoted a slightly improved catalytic efficiency of the former on 3 out of 12 substrates, notably the bulky cephalosporins. Thus, MbetaL-mediated resistance also may be modulated by remote mutations.

Analysis of Mesorhizobium loti glycogen operon: effect of phosphoglucomutase (pgm) and glycogen synthase (g/gA) null mutants on nodulation of Lotus tenuis.

The phosphoglucomutase (pgm) gene codes for a key enzyme required for the formation of UDP-glucose and ADP-glucose, the sugar donors for the biosynthesis of glucose containing polysaccharides. A Mesorhizobium loti pgm null mutant obtained in this study contains an altered form of lipopolysaccharide (LPS), lacks exopolysaccharide (EPS), beta cyclic glucan, and glycogen and is unable to nodulate Lotus tenuis. The nonnodulating phenotype of the pgm mutant was not due to the absence of glycogen, since a glycogen synthase (glgA) null mutant effectively nodulates this legume. In M. loti, pgm is part of the glycogen metabolism gene cluster formed by GlgP (glycogen phosphorylase), glgB (glycogen branching), glgC (ADP-glucose pyrophosphorylase), glgA, pgm, and glgX (glycogen debranching). The genes are transcribed as a single transcript from glgP to at least pgm under the control of a strong promoter (promoter I) upstream of glgP. An alternative promoter (promoter II), mapping in a 154-bp DNA fragment spanning 85 bp upstream of the glgA start codon and the first 69 bp of the glgA coding region, controls the expression of glgA and pgm, independently of the rest of the upstream genes. Primer extension experiments showed that transcription starts 19 bp upstream of the glgA start codon.

Sequence-function-stability relationships in proteins from datasets of functionally annotated variants: the case of TEM ß-lactamases.

A dataset of TEM lactamase variants with different substrate and inhibition profiles was compiled and analyzed. Trends show that loops are the main evolvable regions in these enzymes, gradually accumulating mutations to generate increasingly complex functions. Notably, many mutations present in evolved enzymes are also found in simpler variants, probably originating functional promiscuity. Following a function-stability tradeoff, the increase in functional complexity driven by accumulation of mutations fosters the incorporation of other stability-restoring substitutions, although our analysis suggests they might not be as "global" as generally accepted and seem instead specific to different networks of protein sites. Finally, we show how this dataset can be used to model functional changes in TEMs based on the physicochemical properties of the amino acids.

A minimalistic approach to identify substrate binding features in B1 Metallo-beta-lactamases.

The 2-oxoazetidinylacetate sodium salt was synthesized as a model of a minimal beta-lactam drug. This compound and the monobactam aztreonam were assayed as substrates of the Metallo-beta-lactamase BcII. None of them was hydrolyzed by the enzyme. While the azetidinone was not able to bind BcII, aztreonam was shown to bind in a nonproductive mode. These results provide an explanation for the unability of Metallo-beta-lactamases to inactive monobactams and give some clues for inhibitor design.

Mimicking natural evolution in metallo-beta-lactamases through second-shell ligand mutations.

Metallo-beta-lactamases (MBLs) represent the latest generation of beta-lactamases. The structural diversity and broad substrate profile of MBLs allow them to confer resistance to most beta-lactam antibiotics. To explore the evolutionary potential of these enzymes, we have subjected the Bacillus cereus MBL (BcII) to a directed evolution scheme, which resulted in an increased hydrolytic efficiency toward cephalexin. A systematic study of the hydrolytic profile, substrate binding, and active-site features of the evolved lactamase reveal that directed evolution has shaped the active site by means of remote mutations to better hydrolyze cephalosporins with small, uncharged C-3 substituents. One of these mutations is found in related enzymes from pathogenic bacteria and is responsible for the increase in that enzyme's hydrolytic profile. The mutations lowered the activation energy of the rate-limiting step rather than improved the affinity of the enzyme toward these substrates. The following conclusions can be made: (i) MBLs are able to expand their substrate spectrum without sacrificing their inherent hydrolytic capabilities; (ii) directed evolution is able to mimic mutations that occur in nature; (iii) the metal-ligand strength is tuned by second-shell mutations, thereby influencing the catalytic efficiency; and (iv) changes in the position of the second Zn(II) ion in MBLs affect the substrate positioning in the active site. Overall, these results show that the evolution of enzymatic catalysis can take place by remote mutations controlling reactivity.