Optimization of Conformational Dynamics in an Epistatic Evolutionary Trajectory

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus

Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.

Epistasis arises from shifting the rate-limiting step during enzyme evolution of a β-lactamase

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a β-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity.

Pervasive epistasis exposes intramolecular networks in adaptive enzyme evolution

Enzyme evolution is characterized by constant alterations of the intramolecular residue networks supporting their functions. The rewiring of these network interactions can give rise to epistasis. As mutations accumulate, the epistasis observed across diverse genotypes may appear idiosyncratic, that is, exhibit unique effects in different genetic backgrounds. Here, we unveil a quantitative picture of the prevalence and patterns of epistasis in enzyme evolution by analyzing 41 fitness landscapes generated from seven enzymes. We show that >94% of all mutational and epistatic effects appear highly idiosyncratic, which greatly distorted the functional prediction of the evolved enzymes. By examining seemingly idiosyncratic changes in epistasis along adaptive trajectories, we expose several instances of higher-order, intramolecular rewiring. Using complementary structural data, we outline putative molecular mechanisms explaining higher-order epistasis along two enzyme trajectories. Our work emphasizes the prevalence of epistasis and provides an approach to exploring this phenomenon through a molecular lens.

Metallo-β-lactamase SMB-1 evolves into a more efficient hydrolase under the selective pressure of meropenem

Metallo-β-lactamases (MβLs) are the primary mechanism of resistance to carbapenem antibiotics. To elucidate how MβLs have evolved with the introduction and use of antibiotics, the mutation and evolution of SMB-1 from Serratia marcescens were investigated in microbial evolution plates containing discontinuous meropenem (MEM) concentration gradients. The results revealed 2-point mutations, A242G and S257R; 1 double-site mutation, C240G/E258G; and 3 frameshift mutations, M5, M12, and M13, which are all missense mutations situated at the C-terminus. Compared with that of the wild-type (WT), the minimum inhibitory concentrations (MICs) of MEM for A242G, C240G/E258G, M5, M12, and M13 increased at least 120-fold, and that of S257R increased 8-fold. The catalytic efficiency kcat/Km increased by 365% and 647%, respectively. Concerning the structural changes, the structure at the active site changed from an ordered structure to an unordered conformation. Simultaneously, the flexibility of loop 1 was enhanced. These changes increased the volume of the active site cavity; thus, this was more conducive to exposing the Zn2+ site, facilitating substrate binding and conversion to products. In A242G, structural changes in Gly-242 can be transmitted to the active region via a network of interactions between the side chains of Gly-242 and the amino acid side chains near the active pocket. Together, these results pointed to the process of persistent drug tolerance and resistance, the SMB-1 enzyme evolved into a more exquisite structure with increased flexibility and stability, and stronger hydrolysis activity via genetic mutations and structural changes.

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

β-Lactam antibiotics are the most important and widely used antibacterial agents across the world. However, the widespread dissemination of β-lactamases among pathogenic bacteria limits the efficacy of β-lactam antibiotics. This has created a major public health crisis. The use of β-lactamase inhibitors has proven useful in restoring the activity of β-lactam antibiotics, yet, effective clinically approved inhibitors against class B metallo-β-lactamases are not available. L1, a class B3 enzyme expressed by Stenotrophomonas maltophilia, is a significant contributor to the β-lactam resistance displayed by this opportunistic pathogen. Structurally, L1 is a tetramer with two elongated loops, α3-β7 and β12-α5, present around the active site of each monomer. Residues in these two loops influence substrate/inhibitor binding. To study how the conformational changes of the elongated loops affect the active site in each monomer, enhanced sampling molecular dynamics simulations were performed, Markov State Models were built, and convolutional variational autoencoder-based deep learning was applied. The key identified residues (D150a, H151, P225, Y227, and R236) were mutated and the activity of the generated L1 variants was evaluated in cell-based experiments. The results demonstrate that there are extremely significant gating interactions between α3-β7 and β12-α5 loops. Taken together, the gating interactions with the conformational changes of the key residues play an important role in the structural remodeling of the active site. These observations offer insights into the potential for novel drug development exploiting these gating interactions.

Bypassing evolutionary dead ends and switching the rate-limiting step of a human immunotherapeutic enzyme

The Trp metabolite kynurenine (KYN) accumulates in numerous solid tumours and mediates potent immunosuppression. Bacterial kynureninases (KYNases), which preferentially degrade kynurenine, can relieve immunosuppression in multiple cancer models, but immunogenicity concerns preclude their clinical use, while the human enzyme (HsKYNase) has very low activity for kynurenine and shows no therapeutic effect. Using fitness selections, we evolved a HsKYNase variant with 27-fold higher activity, beyond which exploration of >30 evolutionary trajectories involving the interrogation of >109 variants led to no further improvements. Introduction of two amino acid substitutions conserved in bacterial KYNases reduced enzyme fitness but potentiated rapid evolution of variants with ~500-fold improved activity and reversed substrate specificity, resulting in an enzyme capable of mediating strong anti-tumour effects in mice. Pre-steady-state kinetics revealed a switch in rate-determining step attributable to changes in both enzyme structure and conformational dynamics. Apart from its clinical significance, our work highlights how rationally designed substitutions can potentiate trajectories that overcome barriers in protein evolution.

Energy landscape reshaped by strain-specific mutations underlies epistasis in NS1 evolution of influenza A virus

Elucidating how individual mutations affect the protein energy landscape is crucial for understanding how proteins evolve. However, predicting mutational effects remains challenging because of epistasis—the nonadditive interactions between mutations. Here, we investigate the biophysical mechanism of strain-specific epistasis in the nonstructural protein 1 (NS1) of influenza A viruses (IAVs). We integrate structural, kinetic, thermodynamic, and conformational dynamics analyses of four NS1s of influenza strains that emerged between 1918 and 2004. Although functionally near-neutral, strain-specific NS1 mutations exhibit long-range epistatic interactions with residues at the p85β-binding interface. We reveal that strain-specific mutations reshaped the NS1 energy landscape during evolution. Using NMR spin dynamics, we find that the strain-specific mutations altered the conformational dynamics of the hidden network of tightly packed residues, underlying the evolution of long-range epistasis. This work shows how near-neutral mutations silently alter the biophysical energy landscapes, resulting in diverse background effects during molecular evolution.

Influenza A virus (IAV) nonstructural protein 1 (NS1) is a multifunctional virulence factor that interacts with several host factors such as phosphatidylinositol-3-kinase (PI3K). NS1 binds specifically to the p85β regulatory subunit of PI3K and subsequently activates PI3K signaling. Here, Kim et al. show that functionally near-neutral, strain-specific NS1 mutations lead to variations in binding kinetics to p85β exhibit long-range epistatic interactions. Applying NMR they provide evidence that the structural dynamics of the NS1 hydrophobic core have evolved over time and contributed to epistasis.

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.

Different Conformations Revealed by NMR Underlie Resistance to Ceftazidime/Avibactam and Susceptibility to Meropenem and Imipenem among D179Y Variants of KPC β-Lactamase

β-Lactamase-mediated resistance to ceftazidime-avibactam (CZA) is a serious limitation in the treatment of Gram-negative bacteria harboring Klebsiella pneumoniae carbapenemase (KPC). Herein, the basis of susceptibility to carbapenems and resistance to ceftazidime (CAZ) and CZA of the D179Y variant of KPC-2 and -3 was explored. First, we determined that resistance to CZA in a laboratory strain of Escherichia coli DH10B was not due to increased expression levels of the variant enzymes, as demonstrated by reverse transcription PCR (RT-PCR). Using timed mass spectrometry, the D179Y variant formed prolonged acyl-enzyme complexes with imipenem (IMI) and meropenem (MEM) in KPC-2 and KPC-3, which could be detected up to 24 h, suggesting that IMI and MEM act as covalent β-lactamase inhibitors more than as substrates for D179Y KPC-2 and -3. This prolonged acyl-enzyme complex of IMI and MEM by D179Y variants was not observed with wild-type (WT) KPCs. CAZ was studied and the D179Y variants also formed acyl-enzyme complexes (1 to 2 h). Thermal denaturation and differential scanning fluorimetry showed that the tyrosine substitution at position 179 destabilized the KPC β-lactamases (KPC-2/3 melting temperature [T_m] of 54 to 55°C versus D179Y T_m of 47.5 to 51°C), and the D179Y protein was 3% disordered compared to KPC-2 at 318 K. Heteronuclear ¹H/¹⁵N-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectroscopy also revealed that the D179Y variant, compared to KPC-2, is partially disordered. Based upon these observations, we discuss the impact of disordering of the Ω loop as a consequence of the D179Y substitution. These conformational changes and disorder in the overall structure as a result of D179Y contribute to this unanticipated phenotype.

Localising individual atoms of tryptophan side chains in the metallo-β-lactamase IMP-1 by pseudocontact shifts from paramagnetic lanthanoid tags at multiple sites

The metallo-β -lactamase IMP-1 features a flexible loop near the active site that assumes different conformations in single crystal structures, which may assist in substrate binding and enzymatic activity. To probe the position of this loop, we labelled the tryptophan residues of IMP-1 with 7-13 C-indole and the protein with lanthanoid tags at three different sites. The magnetic susceptibility anisotropy (Δ χ ) tensors were determined by measuring pseudocontact shifts (PCSs) of backbone amide protons. The Δ χ tensors were subsequently used to identify the atomic coordinates of the tryptophan side chains in the protein. The PCSs were sufficient to determine the location of Trp28, which is in the active site loop targeted by our experiments, with high accuracy. Its average atomic coordinates showed barely significant changes in response to the inhibitor captopril. It was found that localisation spaces could be defined with better accuracy by including only the PCSs of a single paramagnetic lanthanoid ion for each tag and tagging site. The effect was attributed to the shallow angle with which PCS isosurfaces tend to intersect if generated by tags and tagging sites that are identical except for the paramagnetic lanthanoid ion.

Second-Shell Amino Acid R266 Helps Determine N-Succinylamino Acid Racemase Reaction Specificity in Promiscuous N-Succinylamino Acid Racemase/o-Succinylbenzoate Synthase Enzymes

Catalytic promiscuity is the coincidental ability to catalyze nonbiological reactions in the same active site as the native biological reaction. Several lines of evidence show that catalytic promiscuity plays a role in the evolution of new enzyme functions. Thus, studying catalytic promiscuity can help identify structural features that predispose an enzyme to evolve new functions. This study identifies a potentially preadaptive residue in a promiscuous N-succinylamino acid racemase/o-succinylbenzoate synthase (NSAR/OSBS) enzyme from Amycolatopsis sp. T-1-60. This enzyme belongs to a branch of the OSBS family which includes many catalytically promiscuous NSAR/OSBS enzymes. R266 is conserved in all members of the NSAR/OSBS subfamily. However, the homologous position is usually hydrophobic in other OSBS subfamilies, whose enzymes lack NSAR activity. The second-shell amino acid R266 is close to the catalytic acid/base K263, but it does not contact the substrate, suggesting that R266 could affect the catalytic mechanism. Mutating R266 to glutamine in Amycolatopsis NSAR/OSBS profoundly reduces NSAR activity but moderately reduces OSBS activity. This is due to a 1000-fold decrease in the rate of proton exchange between the substrate and the general acid/base catalyst K263. This mutation is less deleterious for the OSBS reaction because K263 forms a cation-π interaction with the OSBS substrate and/or the intermediate, rather than acting as a general acid/base catalyst. Together, the data explain how R266 contributes to NSAR reaction specificity and was likely an essential preadaptation for the evolution of NSAR activity.

Two β-Lactamase Variants with Reduced Clavulanic Acid Inhibition Display Different Millisecond Dynamics

The β-lactamase of Mycobacterium tuberculosis, BlaC, is susceptible to inhibition by clavulanic acid. The ability of this enzyme to escape inhibition through mutation was probed using error-prone PCR combined with functional screening in Escherichia coli. The variant that was found to confer the most inhibitor resistance, K234R, as well as variant G132N that was found previously were characterized using X-ray crystallography and nuclear magnetic resonance (NMR) relaxation experiments to probe structural and dynamic properties. The G132N mutant exists in solution in two almost equally populated conformations that exchange with a rate of ca. 88 s-1. The conformational change affects a broad region of the enzyme. The crystal structure reveals that the Asn132 side chain forces the peptide bond between Ser104 and Ile105 in a cis-conformation. The crystal structure suggests multiple conformations for several side chains (e.g., Ser104 and Ser130) and a short loop (positions 214 to 216). In the K234R mutant, the active-site dynamics are significantly diminished with respect to the wild-type enzyme. These results show that multiple evolutionary routes are available to increase inhibitor resistance in BlaC and that active-site dynamics on the millisecond time scale are not required for catalytic function.

Metallo-β-lactamases in the Age of Multidrug Resistance: From Structure and Mechanism to Evolution, Dissemination, and Inhibitor Design

Antimicrobial resistance is one of the major problems in current practical medicine. The spread of genes coding for resistance determinants among bacteria challenges the use of approved antibiotics, narrowing the options for treatment. Resistance to carbapenems, last resort antibiotics, is a major concern. Metallo-β-lactamases (MBLs) hydrolyze carbapenems, penicillins, and cephalosporins, becoming central to this problem. These enzymes diverge with respect to serine-β-lactamases by exhibiting a different fold, active site, and catalytic features. Elucidating their catalytic mechanism has been a big challenge in the field that has limited the development of useful inhibitors. This review covers exhaustively the details of the active-site chemistries, the diversity of MBL alleles, the catalytic mechanism against different substrates, and how this information has helped developing inhibitors. We also discuss here different aspects critical to understand the success of MBLs in conferring resistance: the molecular determinants of their dissemination, their cell physiology, from the biogenesis to the processing involved in the transit to the periplasm, and the uptake of the Zn(II) ions upon metal starvation conditions, such as those encountered during an infection. In this regard, the chemical, biochemical and microbiological aspects provide an integrative view of the current knowledge of MBLs.

Reviewing Challenges of Predicting Protein Melting Temperature Change Upon Mutation Through the Full Analysis of a Highly Detailed Dataset with High-Resolution Structures

Predicting the effects of mutations on protein stability is a key problem in fundamental and applied biology, still unsolved even for the relatively simple case of small, soluble, globular, monomeric, two-state-folder proteins. Many articles discuss the limitations of prediction methods and of the datasets used to train them, which result in low reliability for actual applications despite globally capturing trends. Here, we review these and other issues by analyzing one of the most detailed, carefully curated datasets of melting temperature change (ΔTm) upon mutation for proteins with high-resolution structures. After examining the composition of this dataset to discuss imbalances and biases, we inspect several of its entries assisted by an online app for data navigation and structure display and aided by a neural network that predicts ΔTm with accuracy close to that of programs available to this end. We pose that the ΔTm predictions of our network, and also likely those of other programs, account only for a baseline-like general effect of each type of amino acid substitution which then requires substantial corrections to reproduce the actual stability changes. The corrections are very different for each specific case and arise from fine structural details which are not well represented in the dataset and which, despite appearing reasonable upon visual inspection of the structures, are hard to encode and parametrize. Based on these observations, additional analyses, and a review of recent literature, we propose recommendations for developers of stability prediction methods and for efforts aimed at improving the datasets used for training. We leave our interactive interface for analysis available online at http://lucianoabriata.altervista.org/papersdata/proteinstability2021/s1626navigation.html so that users can further explore the dataset and baseline predictions, possibly serving as a tool useful in the context of structural biology and protein biotechnology research and as material for education in protein biophysics.

Epistasis and intramolecular networks in protein evolution

Proteins are molecular machines composed of complex, highly connected amino acid networks. Their functional optimization requires the reorganization of these intramolecular networks by evolution. In this review, we discuss the mechanisms by which epistasis, that is, the dependence of the effect of a mutation on the genetic background, rewires intramolecular interactions to alter protein function. Deciphering the biophysical basis of epistasis is crucial to our understanding of evolutionary dynamics and the elucidation of sequence-structure-function relationships. We featured recent studies that provide insights into the molecular mechanisms giving rise to epistasis, particularly at the structural level. These studies illustrate the convoluted and fascinating nature of the intramolecular networks co-opted by epistasis during the evolution of protein function.

Evolutionary Trajectory of the Tet(X) Family: Critical Residue Changes towards High-Level Tigecycline Resistance

The emergence of the plasmid-mediated high-level tigecycline resistance mechanism Tet(X) threatens the role of tigecycline as the "last-resort" antibiotic in the treatment of infections caused by carbapenem-resistant Gram-negative bacteria. Compared with that of the prototypical Tet(X), the enzymatic activities of Tet(X3) and Tet(X4) were significantly enhanced, correlating with high-level tigecycline resistance, but the underlying mechanisms remain unclear. In this study, we probed the key amino acid changes leading to the enhancement of Tet(X) function and clarified the structural characteristics and evolutionary path of Tet(X) based upon the key residue changes. Through domain exchange and site-directed mutagenesis experiments, we successfully identified five candidate residues mutations (L282S, A339T, D340N, V350I, and K351E), involved in Tet(X2) activity enhancement. Importantly, these 5 residue changes were 100% conserved among all reported high-activity Tet(X) orthologs, Tet(X3) to Tet(X7), suggesting the important role of these residue changes in the molecular evolution of Tet(X). Structural analysis suggested that the mutant residues did not directly participate in the substrate and flavin adenine dinucleotide (FAD) recognition or binding, but indirectly altered the conformational dynamics of the enzyme through the interaction with adjacent residues. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) and UV full-wavelength scanning experiments confirmed that each mutation led to an increase in activity without changing the biochemical properties of the Tet(X) enzyme. Further phylogenetic analysis suggested that Riemerella anatipestifer served as an important incubator and a main bridge vector for the resistance enhancement and spread of Tet(X). This study expands the knowledge of the structure and function of Tet(X) and provides insights into the evolutionary relationship between Tet(X) orthologs.IMPORTANCE The newly emerged tigecycline-inactivating enzymes Tet(X3) and Tet(X4), which are associated with high-level tigecycline resistance, demonstrated significantly higher activities in comparison to that of the prototypical Tet(X) enzyme, threatening the clinical efficacy of tigecycline as a last-resort antibiotic to treat multidrug-resistant (MDR) Gram-negative bacterial infections. However, the molecular mechanisms leading to high-level tigecycline resistance remain elusive. Here, we identified 5 key residue changes that lead to enhanced Tet(X) activity through domain swapping and site-directed mutagenesis. Instead of direct involvement with substrate binding or catalysis, these residue changes indirectly alter the conformational dynamics and allosterically affect enzyme activities. These findings further broaden the understanding of the structural characteristics and functional evolution of Tet(X) and provide a basis for the subsequent screening of specific inhibitors and the development of novel tetracycline antibiotics.

Molecular Evolution of Transition Metal Bioavailability at the Host-Pathogen Interface

The molecular evolution of the adaptive response at the host-pathogen interface has been frequently referred to as an 'arms race' between the host and bacterial pathogens. The innate immune system employs multiple strategies to starve microbes of metals. Pathogens, in turn, develop successful strategies to maintain access to bioavailable metal ions under conditions of extreme restriction of transition metals, or nutritional immunity. However, the processes by which evolution repurposes or re-engineers host and pathogen proteins to perform or refine new functions have been explored only recently. Here we review the molecular evolution of several human metalloproteins charged with restricting bacterial access to transition metals. These include the transition metal-chelating S100 proteins, natural resistance-associated macrophage protein-1 (NRAMP-1), transferrin, lactoferrin, and heme-binding proteins. We examine their coevolution with bacterial transition metal acquisition systems, involving siderophores and membrane-spanning metal importers, and the biological specificity of allosteric transcriptional regulatory proteins tasked with maintaining bacterial metallostasis. We also discuss the evolution of metallo-β-lactamases; this illustrates how rapid antibiotic-mediated evolution of a zinc metalloenzyme obligatorily occurs in the context of host-imposed nutritional immunity.

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.

The evolution and engineering of enzyme activity through tuning conformational landscapes

Proteins are dynamic molecules whose structures consist of an ensemble of conformational states. Dynamics contribute to protein function and a link to protein evolution has begun to emerge. This increased appreciation for the evolutionary impact of conformational sampling has grown from our developing structural biology capabilities and the exploration of directed evolution approaches, which have allowed evolutionary trajectories to be mapped. Recent studies have provided empirical examples of how proteins can evolve via conformational landscape alterations. Moreover, minor conformational substates have been shown to be involved in the emergence of new enzyme functions as they can become enriched through evolution. The role of remote mutations in stabilizing new active site geometries has also granted insight into the molecular basis underpinning poorly understood epistatic effects that guide protein evolution. Finally, we discuss how the growth of our understanding of remote mutations is beginning to refine our approach to engineering enzymes.

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.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

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.

Metallo-β-Lactamase Inhibitors Inspired on Snapshots from the Catalytic Mechanism

β-Lactam antibiotics are the most widely prescribed antibacterial drugs due to their low toxicity and broad spectrum. Their action is counteracted by different resistance mechanisms developed by bacteria. Among them, the most common strategy is the expression of β-lactamases, enzymes that hydrolyze the amide bond present in all β-lactam compounds. There are several inhibitors against serine-β-lactamases (SBLs). Metallo-β-lactamases (MBLs) are Zn(II)-dependent enzymes able to hydrolyze most β-lactam antibiotics, and no clinically useful inhibitors against them have yet been approved. Despite their large structural diversity, MBLs have a common catalytic mechanism with similar reaction species. Here, we describe a number of MBL inhibitors that mimic different species formed during the hydrolysis process: substrate, transition state, intermediate, or product. Recent advances in the development of boron-based and thiol-based inhibitors are discussed in the light of the mechanism of MBLs. We also discuss the use of chelators as a possible strategy, since Zn(II) ions are essential for substrate binding and catalysis.

β-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.

Active Site-Induced Evolutionary Constraints Follow Fold Polarity Principles in Soluble Globular Enzymes

A recent analysis of evolutionary rates in >500 globular soluble enzymes revealed pervasive conservation gradients toward catalytic residues. By looking at amino acid preference profiles rather than evolutionary rates in the same data set, we quantified the effects of active sites on site-specific constraints for physicochemical traits. We found that conservation gradients respond to constraints for polarity, hydrophobicity, flexibility, rigidity and structure in ways consistent with fold polarity principles; while sites far from active sites seem to experience no physicochemical constraint, rather being highly variable and favoring amino acids of low metabolic cost. Globally, our results highlight that amino acid variation contains finer information about protein structure than usually regarded in evolutionary models, and that this information is retrievable automatically with simple fits. We propose that analyses of the kind presented here incorporated into models of protein evolution should allow for better description of the physical chemistry that underlies molecular evolution.

Population Structure, Molecular Epidemiology, and β-Lactamase Diversity among Stenotrophomonas maltophilia Isolates in the United States

Stenotrophomonas maltophilia is a Gram-negative, nonfermenting, environmental bacillus that is an important cause of nosocomial infections, primarily associated with the respiratory tract in the immunocompromised population. Aiming to understand the population structure, microbiological characteristics and impact of allelic variation on β-lactamase structure and function, we collected 130 clinical isolates from across the United States. Identification of 90 different sequence types (STs), of which 63 are new allelic combinations, demonstrates the high diversity of this species. The majority of the isolates (45%) belong to genomic group 6. We also report excellent activity of the ceftazidime-avibactam and aztreonam combination, especially against strains recovered from blood and respiratory infections for which the susceptibility is higher than the susceptibility to trimethoprim-sulfamethoxazole, considered the "first-line" antibiotic to treat S. maltophilia Analysis of 73 blaL1 and 116 blaL2 genes identified 35 and 43 novel variants of L1 and L2 β-lactamases, respectively. Investigation of the derived amino acid sequences showed that substitutions are mostly conservative and scattered throughout the protein, preferentially affecting positions that do not compromise enzyme function but that may have an impact on substrate and inhibitor binding. Interestingly, we detected a probable association between a specific type of L1 and L2 and genomic group 6. Taken together, our results provide an overview of the molecular epidemiology of S. maltophilia clinical strains from the United States. In particular, the discovery of new L1 and L2 variants warrants further study to fully understand the relationship between them and the β-lactam resistance phenotype in this pathogen.IMPORTANCE Multiple antibiotic resistance mechanisms, including two β-lactamases, L1, a metallo-β-lactamase, and L2, a class A cephalosporinase, make S. maltophilia naturally multidrug resistant. Thus, infections caused by S. maltophilia pose a big therapeutic challenge. Our study aims to understand the microbiological and molecular characteristics of S. maltophilia isolates recovered from human sources. A highlight of the resistance profile of this collection is the excellent activity of the ceftazidime-avibactam and aztreonam combination. We hope this result prompts controlled and observational studies to add clinical data on the utility and safety of this therapy. We also identify 35 and 43 novel variants of L1 and L2, respectively, some of which harbor novel substitutions that could potentially affect substrate and/or inhibitor binding. We believe our results provide valuable knowledge to understand the epidemiology of this species and to advance mechanism-based inhibitor design to add to the limited arsenal of antibiotics active against this pathogen.

Unintended specificity of an engineered ligand-binding protein facilitated by unpredicted plasticity of the protein fold

Attempts to create novel ligand-binding proteins often focus on formation of a binding pocket with shape complementarity against the desired ligand (particularly for compounds that lack distinct polar moieties). Although designed proteins often exhibit binding of the desired ligand, in some cases they display unintended recognition behavior. One such designed protein, that was originally intended to bind tetrahydrocannabinol (THC), was found instead to display binding of 25-hydroxy-cholecalciferol (25-D3) and was subjected to biochemical characterization, further selections for enhanced 25-D3 binding affinity and crystallographic analyses. The deviation in specificity is due in part to unexpected altertion of its conformation, corresponding to a significant change of the orientation of an α-helix and an equally large movement of a loop, both of which flank the designed ligand-binding pocket. Those changes led to engineered protein constructs that exhibit significantly more contacts and complementarity towards the 25-D3 ligand than the initial designed protein had been predicted to form towards its intended THC ligand. Molecular dynamics simulations imply that the initial computationally designed mutations may contribute to the movement of the helix. These analyses collectively indicate that accurate prediction and control of backbone dynamics conformation, through a combination of improved conformational sampling and/or de novo structure design, represents a key area of further development for the design and optimization of engineered ligand-binding proteins.

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.

A DNA aptamer reveals an allosteric site for inhibition in metallo-β-lactamases

The hydrolysis of β-lactam antibiotics by β-lactamase enzymes is the most prominent antibiotic resistance mechanism for many pathogenic bacteria. Out of this broad class of enzymes, metallo-β-lactamases are of special clinical interest because of their broad substrate specificities. Several in vitro inhibitors for various metallo-β-lactamases have been reported with no clinical efficacy. Previously, we described a 10-nucleotide single stranded DNA aptamer (10-mer) that inhibits Bacillus cereus 5/B/6 metallo-β-lactamase very effectively. Here, we find that the aptamer shows uncompetitive inhibition of Bacillus cereus 5/B/6 metallo-β-lactamase during cefuroxime hydrolysis. To understand the mechanism of inhibition, we report a 2.5 Å resolution X-ray crystal structure and solution-state NMR analysis of the free enzyme. Chemical shift perturbations were observed in the HSQC spectra for several residues upon titrating with increasing concentrations of the 10-mer. In the X-ray crystal structure, these residues are distal to the active site, suggesting an allosteric mechanism for the aptamer inhibition of the enzyme. HADDOCK molecular docking simulations suggest that the 10-mer docks 26 Å from the active site. We then mutated the three lysine residues in the basic binding patch to glutamine and measured the catalytic activity and inhibition by the 10-mer. No significant inhibition of these mutants was observed by the 10-mer as compared to wild type. Interestingly, mutation of Lys50 (Lys78; according to standard MBL numbering system) resulted in reduced enzymatic activity relative to wild type in the absence of inhibitor, further highlighting an allosteric mechanism for inhibition.

Kinetic Profile and Molecular Dynamic Studies Show that Y229W Substitution in an NDM-1/L209F Variant Restores the Hydrolytic Activity of the Enzyme toward Penicillins, Cephalosporins, and Carbapenems

The New Delhi metallo-β-lactamase-1 (NDM-1) enzyme is the most common metallo-β-lactamase identified in many Gram-negative bacteria causing severe nosocomial infections. The aim of this study was to focus the attention on non-active-site residues L209 and Y229 of NDM-1 and to investigate their role in the catalytic mechanism. Specifically, the effect of the Y229W substitution in the L209F variant was evaluated by antimicrobial susceptibility testing, kinetic, and molecular dynamic (MD) studies. The Y229W single mutant and L209F-Y229W double mutant were generated by site-directed mutagenesis. The K_m , k _cat, and k _cat/K_m kinetic constants, calculated for the two mutants, were compared with those of (wild-type) NDM-1 and the L209F variant. Compared to the L209F single mutant, the L209F-Y229W double mutant showed a remarkable increase in k _cat values of 100-, 240-, 250-, and 420-fold for imipenem, meropenem, benzylpenicillin, and cefepime, respectively. In the L209F-Y229W enzyme, we observed a remarkable increase in k _cat/K_m of 370-, 140-, and 80-fold for cefepime, meropenem, and cefazolin, respectively. The same behavior was noted using the antimicrobial susceptibility test. MD simulations were carried out on both L209F and L209F-Y229W enzymes complexed with benzylpenicillin, focusing attention on the overall mechanical features and on the differences between the two systems. With respect to the L209F variant, the L209F-Y229W double mutant showed mechanical stabilization of loop 10 and the N-terminal region. In addition, Y229W substitution destabilized both the C-terminal region and the region from residues 149 to 154. The epistatic effect of the Y229W mutation jointly with the stabilization of loop 10 led to a better catalytic efficiency of β-lactams. NDM numbering is used in order to facilitate the comparison with other NDM-1 studies.

The Reaction Mechanism of Metallo-β-Lactamases Is Tuned by the Conformation of an Active-Site Mobile Loop

Carbapenems are "last resort" β-lactam antibiotics used to treat serious and life-threatening health care-associated infections caused by multidrug-resistant Gram-negative bacteria. Unfortunately, the worldwide spread of genes coding for carbapenemases among these bacteria is threatening these life-saving drugs. Metallo-β-lactamases (MβLs) are the largest family of carbapenemases. These are Zn(II)-dependent hydrolases that are active against almost all β-lactam antibiotics. Their catalytic mechanism and the features driving substrate specificity have been matter of intense debate. The active sites of MβLs are flanked by two loops, one of which, loop L3, was shown to adopt different conformations upon substrate or inhibitor binding, and thus are expected to play a role in substrate recognition. However, the sequence heterogeneity observed in this loop in different MβLs has limited the generalizations about its role. Here, we report the engineering of different loops within the scaffold of the clinically relevant carbapenemase NDM-1. We found that the loop sequence dictates its conformation in the unbound form of the enzyme, eliciting different degrees of active-site exposure. However, these structural changes have a minor impact on the substrate profile. Instead, we report that the loop conformation determines the protonation rate of key reaction intermediates accumulated during the hydrolysis of different β-lactams in all MβLs. This study demonstrates the existence of a direct link between the conformation of this loop and the mechanistic features of the enzyme, bringing to light an unexplored function of active-site loops on MβLs.

Synergistic effects of functionally distinct substitutions in β-lactamase variants shed light on the evolution of bacterial drug resistance

The CTX-M β-lactamases have emerged as the most widespread extended-spectrum β-lactamases (ESBLs) in Gram-negative bacteria. These enzymes rapidly hydrolyze cefotaxime, but not the related cephalosporin, ceftazidime. ESBL variants have evolved, however, that provide enhanced ceftazidime resistance. We show here that a natural variant at a nonactive site, i.e. second-shell residue N106S, enhances enzyme stability but reduces catalytic efficiency for cefotaxime and ceftazidime and decreases resistance levels. However, when the N106S variant was combined with an active-site variant, D240G, that enhances enzyme catalytic efficiency, but decreases stability, the resultant double mutant exhibited higher resistance levels than predicted on the basis of the phenotypes of each variant. We found that this epistasis is due to compensatory effects, whereby increased stability provided by N106S overrides its cost of decreased catalytic activity. X-ray structures of the variant enzymes in complex with cefotaxime revealed conformational changes in the active-site loop spanning residues 103-106 that were caused by the N106S substitution and relieve steric strain to stabilize the enzyme, but also alter contacts with cefotaxime and thereby reduce catalytic activity. We noted that the 103-106 loop conformation in the N106S-containing variants is different from that of WT CTX-M but nearly identical to that of the non-ESBL, TEM-1 β-lactamase, having a serine at the 106 position. Therefore, residue 106 may serve as a "switch" that toggles the conformations of the 103-106 loop. When it is serine, the loop is in the non-ESBL, TEM-like conformation, and when it is asparagine, the loop is in a CTX-M-like, cefotaximase-favorable conformation.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

The evolution of substrate discrimination in macrolide antibiotic resistance enzymes

The production of antibiotics by microbes in the environment and their use in medicine and agriculture select for existing and emerging resistance. To address this inevitability, prudent development of antibiotic drugs requires careful consideration of resistance evolution. Here, we identify the molecular basis for expanded substrate specificity in MphI, a macrolide kinase (Mph) that does not confer resistance to erythromycin, in contrast to other known Mphs. Using a combination of phylogenetics, drug-resistance phenotypes, and in vitro enzyme assays, we find that MphI and MphK phosphorylate erythromycin poorly resulting in an antibiotic-sensitive phenotype. Using likelihood reconstruction of ancestral sequences and site-saturation combinatorial mutagenesis, supported by Mph crystal structures, we determine that two non-obvious mutations in combination expand the substrate range. This approach should be applicable for studying the functional evolution of any antibiotic resistance enzyme and for evaluating the evolvability of resistance enzymes to new generations of antibiotic scaffolds.

New antibiotics with reduced potential for resistance are urgently needed. Here, the authors use a multidisciplinary approach to characterize substrate discrimination in macrolide resistance kinases and present a strategy for the prediction of mutations that expand the substrate range of antibiotic-inactivating enzymes.

Enzyme Efficiency but Not Thermostability Drives Cefotaxime Resistance Evolution in TEM-1 β-Lactamase

A leading intellectual challenge in evolutionary genetics is to identify the specific phenotypes that drive adaptation. Enzymes offer a particularly promising opportunity to pursue this question, because many enzymes' contributions to organismal fitness depend on a comparatively small number of experimentally accessible properties. Moreover, on first principles the demands of enzyme thermostability stand in opposition to the demands of catalytic activity. This observation, coupled with the fact that enzymes are only marginally thermostable, motivates the widely held hypothesis that mutations conferring functional improvement require compensatory mutations to restore thermostability. Here, we explicitly test this hypothesis for the first time, using four missense mutations in TEM-1 β-lactamase that jointly increase cefotaxime Minimum Inhibitory Concentration (MIC) ∼1500-fold. First, we report enzymatic efficiency (kcat/KM) and thermostability (Tm, and thence ΔG of folding) for all combinations of these mutations. Next, we fit a quantitative model that predicts MIC as a function of kcat/KM and ΔG. While kcat/KM explains ∼54% of the variance in cefotaxime MIC (∼92% after log transformation), ΔG does not improve explanatory power of the model. We also find that cefotaxime MIC rises more slowly in kcat/KM than predicted. Several explanations for these discrepancies are suggested. Finally, we demonstrate substantial sign epistasis in MIC and kcat/KM, and antagonistic pleiotropy between phenotypes, in spite of near numerical additivity in the system. Thus constraints on selectively accessible trajectories, as well as limitations in our ability to explain such constraints in terms of underlying mechanisms are observed in a comparatively "well-behaved" system.

Molecular Effects of Concentrated Solutes on Protein Hydration, Dynamics, and Electrostatics

Most studies of protein structure and function are performed in dilute conditions, but proteins typically experience high solute concentrations in their physiological scenarios and biotechnological applications. High solute concentrations have well-known effects on coarse protein traits like stability, diffusion, and shape, but likely also perturb other traits through finer effects pertinent at the residue and atomic levels. Here, NMR and molecular dynamics investigations on ubiquitin disclose variable interactions with concentrated solutes that lead to localized perturbations of the protein's surface, hydration, electrostatics, and dynamics, all dependent on solute size and chemical properties. Most strikingly, small polar uncharged molecules are sticky on the protein surface, whereas charged small molecules are not, but the latter still perturb the internal protein electrostatics as they diffuse nearby. Meanwhile, interactions with macromolecular crowders are favored mainly through hydrophobic, but not through polar, surface patches. All the tested small solutes strongly slow down water exchange at the protein surface, whereas macromolecular crowders do not exert such strong perturbation. Finally, molecular dynamics simulations predict that unspecific interactions slow down microsecond- to millisecond-timescale protein dynamics despite having only mild effects on pico- to nanosecond fluctuations as corroborated by NMR. We discuss our results in the light of recent advances in understanding proteins inside living cells, focusing on the physical chemistry of quinary structure and cellular organization, and we reinforce the idea that proteins should be studied in native-like media to achieve a faithful description of their function.

Insights into an evolutionary strategy leading to antibiotic resistance

Metallo-β-lactamases (MBLs) with activity towards a broad-spectrum of β-lactam antibiotics have become a major threat to public health, not least due to their ability to rapidly adapt their substrate preference. In this study, the capability of the MBL AIM-1 to evade antibiotic pressure by introducing specific mutations was probed by two alternative methods, i.e. site-saturation mutagenesis (SSM) of active site residues and in vitro evolution. Both approaches demonstrated that a single mutation in AIM-1 can greatly enhance a pathogen's resistance towards broad spectrum antibiotics without significantly compromising the catalytic efficiency of the enzyme. Importantly, the evolution experiments demonstrated that relevant amino acids are not necessarily in close proximity to the catalytic centre of the enzyme. This observation is a powerful demonstration that MBLs have a diverse array of possibilities to adapt to new selection pressures, avenues that cannot easily be predicted from a crystal structure alone.

Distinct Roles for Conformational Dynamics in Protein-Ligand Interactions

Conformational dynamics has an established role in enzyme catalysis, but its contribution to ligand binding and specificity is largely unexplored. Here we used the Tiam1 PDZ domain and an engineered variant (QM PDZ) with broadened specificity to investigate the role of structure and conformational dynamics in molecular recognition. Crystal structures of the QM PDZ domain both free and bound to ligands showed structural features central to binding (enthalpy), while nuclear-magnetic-resonance-based methyl relaxation experiments and isothermal titration calorimetry revealed that conformational entropy contributes to affinity. In addition to motions relevant to thermodynamics, slower microsecond to millisecond switching was prevalent in the QM PDZ ligand-binding site consistent with a role in ligand specificity. Our data indicate that conformational dynamics plays distinct and fundamental roles in tuning the affinity (conformational entropy) and specificity (excited-state conformations) of molecular interactions. More broadly, our results have important implications for the evolution, regulation, and design of protein-ligand interactions.

Detection and sequence/structure mapping of biophysical constraints to protein variation in saturated mutational libraries and protein sequence alignments with a dedicated server

Background

Protein variability can now be studied by measuring high-resolution tolerance-to-substitution maps and fitness landscapes in saturated mutational libraries. But these rich and expensive datasets are typically interpreted coarsely, restricting detailed analyses to positions of extremely high or low variability or dubbed important beforehand based on existing knowledge about active sites, interaction surfaces, (de)stabilizing mutations, etc.

Results

Our new webserver PsychoProt (freely available without registration at http://psychoprot.epfl.ch or at http://lucianoabriata.altervista.org/psychoprot/index.html) helps to detect, quantify, and sequence/structure map the biophysical and biochemical traits that shape amino acid preferences throughout a protein as determined by deep-sequencing of saturated mutational libraries or from large alignments of naturally occurring variants.

Discussion

We exemplify how PsychoProt helps to (i) unveil protein structure-function relationships from experiments and from alignments that are consistent with structures according to coevolution analysis, (ii) recall global information about structural and functional features and identify hitherto unknown constraints to variation in alignments, and (iii) point at different sources of variation among related experimental datasets or between experimental and alignment-based data. Remarkably, metabolic costs of the amino acids pose strong constraints to variability at protein surfaces in nature but not in the laboratory. This and other differences call for caution when extrapolating results from in vitro experiments to natural scenarios in, for example, studies of protein evolution.

Conclusion

We show through examples how PsychoProt can be a useful tool for the broad communities of structural biology and molecular evolution, particularly for studies about protein modeling, evolution and design.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1124-4) contains supplementary material, which is available to authorized users.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.