Protein dynamics at slow timescales in engineered ß-lactamases does not limit evolvability

Understanding the underlying mechanisms in the evolution of new protein functions is key to better directing enzyme engineering efforts. Intragenic epistasis (the non-additive interaction of mutations affecting function) is a key feature of protein evolution. For example, in TEM-1 ß-lactamase, the mutations E104K and G238S show positive epistasis in that their combination show a greater than expected increase in antibiotic resistance1. Here, we aim to understand the impact of protein dynamics at slow timescales on epistasis. Large conformational rearrangements associated with ligand-binding, turnover of substrate or allostery occur at this timescale2. The readily evolvable antibiotic-resistant TEM-1 ß-lactamase and two of its variants that are more dynamic at slow timescales serve as models for this study. Our models show similar catalytic activity and substrate recognition, thermal stability, as well as conserved motions in fast (ps-ns) and intermediate (ns-ms) timescales but different motions at slow timescales (ms-ms)123. Please click Additional Files below to see the full abstract

The Bacterial Genomic Context of Highly Trimethoprim-Resistant DfrB Dihydrofolate Reductases Highlights an Emerging Threat to Public Health

Type B dihydrofolate reductase (dfrb) genes were identified following the introduction of trimethoprim in the 1960s. Although they intrinsically confer resistance to trimethoprim (TMP) that is orders of magnitude greater than through other mechanisms, the distribution and prevalence of these short (237 bp) genes is unknown. Indeed, this knowledge has been hampered by systematic biases in search methodologies. Here, we investigate the genomic context of dfrbs to gain information on their current distribution in bacterial genomes. Upon searching publicly available databases, we identified 61 sequences containing dfrbs within an analyzable genomic context. The majority (70%) of those sequences also harbor virulence genes and 97% of the dfrbs are found near a mobile genetic element, representing a potential risk for antibiotic resistance genes. We further identified and confirmed the TMP-resistant phenotype of two new members of the family, dfrb10 and dfrb11. Dfrbs are found both in Betaproteobacteria and Gammaproteobacteria, a majority (59%) being in Pseudomonas aeruginosa. Previously labelled as strictly plasmid-borne, we found 69% of dfrbs in the chromosome of pathogenic bacteria. Our results demonstrate that the intrinsically TMP-resistant dfrbs are a potential emerging threat to public health and justify closer surveillance of these genes

Known evolutionary paths are accessible to engineered ß-Lactamases having altered protein motions at the timescale of catalytic turnover

The evolution of new protein functions is dependent upon inherent biophysical features of proteins. Whereas, it has been shown that changes in protein dynamics can occur in the course of directed molecular evolution trajectories and contribute to new function, it is not known whether varying protein dynamics modify the course of evolution. We investigate this question using three related ß-lactamases displaying dynamics that differ broadly at the slow timescale that corresponds to catalytic turnover yet have similar fast dynamics, thermal stability, catalytic, and substrate recognition profiles. Introduction of substitutions E104K and G238S, that are known to have a synergistic effect on function in the parent ß-lactamase, showed similar increases in catalytic efficiency toward cefotaxime in the related ß-lactamases. Molecular simulations using Protein Energy Landscape Exploration reveal that this results from stabilizing the catalytically-productive conformations, demonstrating the dominance of the synergistic effect of the E014K and G238S substitutions in vitro in contexts that vary in terms of sequence and dynamics. Furthermore, three rounds of directed molecular evolution demonstrated that known cefotaximase-enhancing mutations were accessible regardless of the differences in dynamics. Interestingly, specific sequence differences between the related ß-lactamases were shown to have a higher effect in evolutionary outcomes than did differences in dynamics. Overall, these ß-lactamase models show tolerance to protein dynamics at the timescale of catalytic turnover in the evolution of a new function.This work was supported by operating grant RGPIN-2018-04686 to JP from the Natural Science and Engineering Research Council of Canada (NSERC), infrastructure grant 11510 to JP from the Canada Foundation for Innovation and operating grant PID2019-106370RB-I00 to VG from the Spanish Ministry of Science and Innovation. LA and CL-S-D are grateful to FQRNT, Université de Montréal, and NSERC for scholarships.Peer ReviewedPostprint (published version

Discovery of Highly Trimethoprim-Resistant DfrB Dihydrofolate Reductases in Diverse Environmental Settings Suggests an Evolutionary Advantage Unrelated to Antibiotic Resistance

Type B dihydrofolate reductases (DfrB) are intrinsically highly resistant to the widely used antibiotic trimethoprim, posing a threat to global public health. The ten known DfrB family members have been strongly associated with genetic material related to the application of antibiotics. Several dfrB genes were associated with multidrug resistance contexts and mobile genetic elements, integrated both in chromosomes and plasmids. However, little is known regarding their presence in other environments. Here, we investigated the presence of dfrB beyond the traditional areas of enquiry by conducting metagenomic database searches from environmental settings where antibiotics are not prevalent. Thirty putative DfrB homologues that share 62 to 95% identity with characterized DfrB were identified. Expression of ten representative homologues verified trimethoprim resistance in all and dihydrofolate reductase activity in most. Contrary to samples associated with the use of antibiotics, the newly identified dfrB were rarely associated with mobile genetic elements or antibiotic resistance genes. Instead, association with metabolic enzymes was observed, suggesting an evolutionary advantage unrelated to antibiotic resistance. Our results are consistent with the hypothesis that multiple dfrB exist in diverse environments from which dfrB were mobilized into the clinically relevant resistome. Our observations reinforce the need to closely monitor their progression