Evolution of an antibiotic resistance enzyme constrained by stability and activity trade-offs

Structural insights into the broadened substrate profile of the extended-spectrum β-lactamase OXY-1-1 fromKlebsiella oxytoca

Klebsiella oxytocais a pathogen that causes serious infections in hospital patients. It shows resistance to many clinically used β-lactam antibiotics by producing chromosomally encoded OXY-family β-lactamases. Here, the crystal structure of an OXY-family β-lactamase, OXY-1-1, determined at 1.93 Å resolution is reported. The structure shows that the OXY-1-1 β-lactamase has a typical class A β-lactamase fold and exhibits greater similarity to CTX-M-type β-lactamases than to TEM-family or SHV-family β-lactamases. It is also shown that the enzyme provides more space around the active cavity for theR1andR2substituents of β-lactam antibiotics. The half-positive/half-negative distribution of surface electrostatic potential in the substrate-binding pocket indicates the preferred properties of substrates or inhibitors of the enzyme. The results reported here provide a structural basis for the broadened substrate profile of the OXY-family β-lactamases.

A fundamental protein property, thermodynamic stability, revealed solely from large-scale measurements of protein function

The ability of a protein to carry out a given function results from fundamental physicochemical properties that include the protein's structure, mechanism of action, and thermodynamic stability. Traditional approaches to study these properties have typically required the direct measurement of the property of interest, oftentimes a laborious undertaking. Although protein properties can be probed by mutagenesis, this approach has been limited by its low throughput. Recent technological developments have enabled the rapid quantification of a protein's function, such as binding to a ligand, for numerous variants of that protein. Here, we measure the ability of 47,000 variants of a WW domain to bind to a peptide ligand and use these functional measurements to identify stabilizing mutations without directly assaying stability. Our approach is rooted in the well-established concept that protein function is closely related to stability. Protein function is generally reduced by destabilizing mutations, but this decrease can be rescued by stabilizing mutations. Based on this observation, we introduce partner potentiation, a metric that uses this rescue ability to identify stabilizing mutations, and identify 15 candidate stabilizing mutations in the WW domain. We tested six candidates by thermal denaturation and found two highly stabilizing mutations, one more stabilizing than any previously known mutation. Thus, physicochemical properties such as stability are latent within these large-scale protein functional data and can be revealed by systematic analysis. This approach should allow other protein properties to be discovered.

Deep sequencing of systematic combinatorial libraries reveals β-lactamase sequence constraints at high resolution

In this study, combinatorial libraries were used in conjunction with ultrahigh-throughput sequencing to comprehensively determine the impact of each of the 19 possible amino acid substitutions at each residue position in the TEM-1 β-lactamase enzyme. The libraries were introduced into Escherichiacoli, and mutants were selected for ampicillin resistance. The selected colonies were pooled and subjected to ultrahigh-throughput sequencing to reveal the sequence preferences at each position. The depth of sequencing provided a clear, statistically significant picture of what amino acids are favored for ampicillin hydrolysis for all 263 positions of the enzyme in one experiment. Although the enzyme is generally tolerant of amino acid substitutions, several surface positions far from the active site are sensitive to substitutions suggesting a role for these residues in enzyme stability, solubility, or catalysis. In addition, information on the frequency of substitutions was used to identify mutations that increase enzyme thermodynamic stability. Finally, a comparison of sequence requirements based on the mutagenesis results versus those inferred from sequence conservation in an alignment of 156 class A β-lactamases reveals significant differences in that several residues in TEM-1 do not tolerate substitutions and yet extensive variation is observed in the alignment and vice versa. An analysis of the TEM-1 and other class A structures suggests that residues that vary in the alignment may nevertheless make unique, but important, interactions within individual enzymes.

Engineering a model protein cavity to catalyze the Kemp elimination

Synthetic cavitands and protein cavities have been widely studied as models for ligand recognition. Here we investigate the Met102 → His substitution in the artificial L99A cavity in T4 lysozyme as a Kemp eliminase. The resulting enzyme had k(cat)/K(M) = 0.43 M(-1) s(-1) and a (k(cat)/K(M))/k(uncat) = 10(7) at pH 5.0. The crystal structure of this enzyme was determined at 1.30 Å, as were the structures of four complexes of substrate and product analogs. The absence of ordered waters or hydrogen bonding interactions, and the presence of a common catalytic base (His102) in an otherwise hydrophobic, buried cavity, facilitated detailed analysis of the reaction mechanism and its optimization. Subsequent substitutions increased eliminase activity by an additional four-fold. As activity-enhancing substitutions were engineered into the cavity, protein stability decreased, consistent with the stability-function trade-off hypothesis. This and related model cavities may provide templates for studying protein design principles in radically simplified environments.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design

Predicting which mutations proteins tolerate while maintaining their structure and function has important applications for modeling fundamental properties of proteins and their evolution; it also drives progress in protein design. Here we develop a computational model to predict the tolerated sequence space of HIV-1 protease reachable by single mutations. We assess the model by comparison to the observed variability in more than 50,000 HIV-1 protease sequences, one of the most comprehensive datasets on tolerated sequence space. We then extend the model to a second protein, reverse transcriptase. The model integrates multiple structural and functional constraints acting on a protein and uses ensembles of protein conformations. We find the model correctly captures a considerable fraction of protease and reverse-transcriptase mutational tolerance and shows comparable accuracy using either experimentally determined or computationally generated structural ensembles. Predictions of tolerated sequence space afforded by the model provide insights into stability-function tradeoffs in the emergence of resistance mutations and into strengths and limitations of the computational model.

Many related protein sequences can be consistent with the structure and function of a given protein, suggesting that proteins may be quite robust to mutations. This tolerance to mutations is frequently exploited by pathogens. In particular, pathogens can rapidly evolve mutated proteins that have a new function - resistance against a therapeutic inhibitor - without abandoning other functions essential for the pathogen. This principle may also hold more generally: Proteins tolerant to mutational changes can more easily acquire new functions while maintaining their existing properties. The ability to predict the tolerance of proteins to mutation could thus help both to analyze the emergence of resistance mutations in pathogens and to engineer proteins with new functions. Here we develop a computational model to predict protein mutational tolerance towards point mutations accessible by single nucleotide changes, and validate it using two important pathogenic proteins and therapeutic targets: the protease and reverse transcriptase from HIV-1. The model provides insights into how resistance emerges and makes testable predictions on mutations that have not been seen yet. Similar models of mutational tolerance should be useful for characterizing and reengineering the functions of other proteins for which a three-dimensional structure is available.

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.

Exploring the role of a conserved class A residue in the Ω-Loop of KPC-2 β-lactamase: a mechanism for ceftazidime hydrolysis

Gram-negative bacteria harboring KPC-2, a class A β-lactamase, are resistant to all β-lactam antibiotics and pose a major public health threat. Arg-164 is a conserved residue in all class A β-lactamases and is located in the solvent-exposed Ω-loop of KPC-2. To probe the role of this amino acid in KPC-2, we performed site-saturation mutagenesis. When compared with wild type, 11 of 19 variants at position Arg-164 in KPC-2 conferred increased resistance to the oxyimino-cephalosporin, ceftazidime (minimum inhibitory concentration; 32→128 mg/liter) when expressed in Escherichia coli. Using the R164S variant of KPC-2 as a representative β-lactamase for more detailed analysis, we observed only a modest 25% increase in k(cat)/K(m) for ceftazidime (0.015→0.019 μm(-1) s(-1)). Employing pre-steady-state kinetics and mass spectrometry, we determined that acylation is rate-limiting for ceftazidime hydrolysis by KPC-2, whereas deacylation is rate-limiting in the R164S variant, leading to accumulation of acyl-enzyme at steady-state. CD spectroscopy revealed that a conformational change occurred in the turnover of ceftazidime by KPC-2, but not the R164S variant, providing evidence for a different form of the enzyme at steady state. Molecular models constructed to explain these findings suggest that ceftazidime adopts a unique conformation, despite preservation of Ω-loop structure. We propose that the R164S substitution in KPC-2 enhances ceftazidime resistance by proceeding through "covalent trapping" of the substrate by a deacylation impaired enzyme with a lower K(m). Future antibiotic design must consider the distinctive behavior of the Ω-loop of KPC-2.

Enhancing the promiscuous phosphotriesterase activity of a thermostable lactonase (GkaP) for the efficient degradation of organophosphate pesticides

The phosphotriesterase-like lactonase (PLL) enzymes in the amidohydrolase superfamily hydrolyze various lactones and exhibit latent phosphotriesterase activities. These enzymes serve as attractive templates for in vitro evolution of neurotoxic organophosphates (OPs) with hydrolytic capabilities that can be used as bioremediation tools. Here, a thermostable PLL from Geobacillus kaustophilus HTA426 (GkaP) was targeted for joint laboratory evolution with the aim of enhancing its catalytic efficiency against OP pesticides. By a combination of site saturation mutagenesis and whole-gene error-prone PCR approaches, several improved variants were isolated. The most active variant, 26A8C, accumulated eight amino acid substitutions and demonstrated a 232-fold improvement over the wild-type enzyme in reactivity (k(cat)/K(m)) for the OP pesticide ethyl-paraoxon. Concomitantly, this variant showed a 767-fold decrease in lactonase activity with δ-decanolactone, imparting a specificity switch of 1.8 × 10(5)-fold. 26A8C also exhibited high hydrolytic activities (19- to 497-fold) for several OP pesticides, including parathion, diazinon, and chlorpyrifos. Analysis of the mutagenesis sites on the GkaP structure revealed that most mutations are located in loop 8, which determines substrate specificity in the amidohydrolase superfamily. Molecular dynamics simulation shed light on why 26A8C lost its native lactonase activity and improved the promiscuous phosphotriesterase activity. These results permit us to obtain further insights into the divergent evolution of promiscuous enzymes and suggest that laboratory evolution of GkaP may lead to potential biological solutions for the efficient decontamination of neurotoxic OP compounds.

Equilibrium fluctuations of a single folded protein reveal a multitude of potential cryptic allosteric sites

Cryptic allosteric sites--transient pockets in a folded protein that are invisible to conventional experiments but can alter enzymatic activity via allosteric communication with the active site--are a promising opportunity for facilitating drug design by greatly expanding the repertoire of available drug targets. Unfortunately, identifying these sites is difficult, typically requiring resource-intensive screening of large libraries of small molecules. Here, we demonstrate that Markov state models built from extensive computer simulations (totaling hundreds of microseconds of dynamics) can identify prospective cryptic sites from the equilibrium fluctuations of three medically relevant proteins--β-lactamase, interleukin-2, and RNase H--even in the absence of any ligand. As in previous studies, our methods reveal a surprising variety of conformations--including bound-like configurations--that implies a role for conformational selection in ligand binding. Moreover, our analyses lead to a number of unique insights. First, direct comparison of simulations with and without the ligand reveals that there is still an important role for an induced fit during ligand binding to cryptic sites and suggests new conformations for docking. Second, correlations between amino acid sidechains can convey allosteric signals even in the absence of substantial backbone motions. Most importantly, our extensive sampling reveals a multitude of potential cryptic sites--consisting of transient pockets coupled to the active site--even in a single protein. Based on these observations, we propose that cryptic allosteric sites may be even more ubiquitous than previously thought and that our methods should be a valuable means of guiding the search for such sites.

Quantifying the adaptive potential of an antibiotic resistance enzyme

For a quantitative understanding of the process of adaptation, we need to understand its "raw material," that is, the frequency and fitness effects of beneficial mutations. At present, most empirical evidence suggests an exponential distribution of fitness effects of beneficial mutations, as predicted for Gumbel-domain distributions by extreme value theory. Here, we study the distribution of mutation effects on cefotaxime (Ctx) resistance and fitness of 48 unique beneficial mutations in the bacterial enzyme TEM-1 β-lactamase, which were obtained by screening the products of random mutagenesis for increased Ctx resistance. Our contributions are threefold. First, based on the frequency of unique mutations among more than 300 sequenced isolates and correcting for mutation bias, we conservatively estimate that the total number of first-step mutations that increase Ctx resistance in this enzyme is 87 [95% CI 75-189], or 3.4% of all 2,583 possible base-pair substitutions. Of the 48 mutations, 10 are synonymous and the majority of the 38 non-synonymous mutations occur in the pocket surrounding the catalytic site. Second, we estimate the effects of the mutations on Ctx resistance by determining survival at various Ctx concentrations, and we derive their fitness effects by modeling reproduction and survival as a branching process. Third, we find that the distribution of both measures follows a Fréchet-type distribution characterized by a broad tail of a few exceptionally fit mutants. Such distributions have fundamental evolutionary implications, including an increased predictability of evolution, and may provide a partial explanation for recent observations of striking parallel evolution of antibiotic resistance.

Bridging the gaps in design methodologies by evolutionary optimization of the stability and proficiency of designed Kemp eliminase KE59

Computational design is a test of our understanding of enzyme catalysis and a means of engineering novel, tailor-made enzymes. While the de novo computational design of catalytically efficient enzymes remains a challenge, designed enzymes may comprise unique starting points for further optimization by directed evolution. Directed evolution of two computationally designed Kemp eliminases, KE07 and KE70, led to low to moderately efficient enzymes (k(cat)/K(m) values of ≤ 5 10(4) M(-1)s(-1)). Here we describe the optimization of a third design, KE59. Although KE59 was the most catalytically efficient Kemp eliminase from this design series (by k(cat)/K(m), and by catalyzing the elimination of nonactivated benzisoxazoles), its impaired stability prevented its evolutionary optimization. To boost KE59's evolvability, stabilizing consensus mutations were included in the libraries throughout the directed evolution process. The libraries were also screened with less activated substrates. Sixteen rounds of mutation and selection led to > 2,000-fold increase in catalytic efficiency, mainly via higher k(cat) values. The best KE59 variants exhibited k(cat)/K(m) values up to 0.6 10(6) M(-1)s(-1), and k(cat)/k(uncat) values of ≤ 10(7) almost regardless of substrate reactivity. Biochemical, structural, and molecular dynamics (MD) simulation studies provided insights regarding the optimization of KE59. Overall, the directed evolution of three different designed Kemp eliminases, KE07, KE70, and KE59, demonstrates that computational designs are highly evolvable and can be optimized to high catalytic efficiencies.

Exploring the molecular linkage of protein stability traits for enzyme optimization by iterative truncation and evolution

The stability of proteins is paramount for their therapeutic and industrial use and, thus, is a major task for protein engineering. Several types of chemical and physical stabilities are desired, and discussion revolves around whether each stability trait needs to be addressed separately and how specific and compatible stabilizing mutations act. We demonstrate a stepwise perturbation-compensation strategy, which identifies mutations rescuing the activity of a truncated TEM β-lactamase. Analyses relating structural stress with the external stresses of heat, denaturants, and proteases reveal our second-site suppressors as general stability centers that also improve the full-length enzyme. A library of lactamase variants truncated by 15 N-terminal and three C-terminal residues (Bla-NΔ15CΔ3) was subjected to activity selection and DNA shuffling. The resulting clone with the best in vivo performance harbored eight mutations, surpassed the full-length wild-type protein by 5.3 °C in T(m), displayed significantly higher catalytic activity at elevated temperatures, and showed delayed guanidine-induced denaturation. The crystal structure of this mutant was determined and provided insights into its stability determinants. Stepwise reconstitution of the N- and C-termini increased its thermal, denaturant, and proteolytic resistance successively, leading to a full-length enzyme with a T(m) increased by 15.3 °C and a half-denaturation concentration shifted from 0.53 to 1.75 M guanidinium relative to that of the wild type. These improvements demonstrate that iterative truncation-optimization cycles can exploit stability-trait linkages in proteins and are exceptionally suited for the creation of progressively stabilized variants and/or downsized proteins without the need for detailed structural or mechanistic information.

Twelve positions in a β-lactamase that can expand its substrate spectrum with a single amino acid substitution

The continuous evolution of β-lactamases resulting in bacterial resistance to β-lactam antibiotics is a major concern in public health, and yet the underlying molecular basis or the pattern of such evolution is largely unknown. We investigated the mechanics of the substrate fspectrum expansion of the class A β-lactamase using PenA of Burkholderia thailandensis as a model. By analyzing 516 mutated enzymes that acquired the ceftazidime-hydrolyzing activity, we found twelve positions with single amino acid substitutions (altogether twenty-nine different substitutions), co-localized at the active-site pocket area. The ceftazidime MIC (minimum inhibitory concentration) levels and the relative frequency in the occurrence of substitutions did not correlate well with each other, and the latter appeared be largely influenced by the intrinsic mutational biases present in bacteria. Simulation studies suggested that all substitutions caused a congruent effect, expanding the space in a conserved structure called the omega loop, which in turn increased flexibility at the active site. A second phase of selection, in which the mutants were placed under increased antibiotic pressure, did not result in a second mutation in the coding region, but a mutation that increased gene expression arose in the promoter. This result suggests that the twelve amino acid positions and their specific substitutions in PenA may represent a comprehensive repertoire of the enzyme's adaptability to a new substrate. These mapped substitutions represent a comprehensive set of general mechanical paths to substrate spectrum expansion in class A β-lactamases that all share a functional evolutionary mechanism using common conserved residues.

Substrate spectrum extension of PenA in Burkholderia thailandensis with a single amino acid deletion, Glu168del

We describe a deletion mutation in a class A β-lactamase, PenA, of Burkholderia thailandensis that extended the substrate spectrum of the enzyme to include ceftazidime. Glu168del was located in a functional domain called the omega loop causing expansion of the space in the loop, which in turn increased flexibility at the active site. This deletion mutation represents a rare but significant alternative mechanical path to substrate spectrum extension in PenA besides more common substitution mutations.

Prospective screening of novel antibacterial inhibitors of dihydrofolate reductase for mutational resistance

Resistance to trimethoprim (TMP) resulting from point mutations in the enzyme drug target dihydrofolate reductase (DHFR) drives the development of new antifolate inhibitors effective against methicillin-resistant Staphylococcus aureus (MRSA). For the past several years we have used structure-based design to create propargyl-linked antifolates that are highly potent antibacterial agents. In order to focus priority on the development of lead compounds with a low propensity to induce resistance, we prospectively evaluated resistance profiles for two of these inhibitors in an MRSA strain. By selection with the lead inhibitors, we generated resistant strains that contain single point mutations F98Y and H30N associated with TMP resistance and one novel mutation, F98I, in DHFR. Encouragingly, the pyridyl propargyl-linked inhibitor selects mutants at low frequency (6.85 × 10(-10) to 1.65 × 10(-9)) and maintains a low MIC (2.5 μg/ml) and a low mutant prevention concentration (1.25 μg/ml), strongly supporting its position as a lead compound. Results from this prospective screening method inform the continued design of antifolates effective against mutations at the Phe 98 position. Furthermore, the method can be used broadly to incorporate ideas for overcoming resistance early in the development process.

Winning the arms race by improving drug discovery against mutating targets

Enzymes are often excellent drug targets. Yet drug pressure on an enzyme target often fosters the rise of cells with resistance-conferring mutations, some of which may compromise fitness and others that compensate to restore fitness. This review presents, first, a structural analysis of a diverse group of wild-type and mutant enzyme targets and, second, an in-depth analysis of five diverse targets to elucidate a broader perspective of the effects of resistance-conferring mutations on protein or organismal fitness. The structural analysis reveals that resistance-conferring mutations may introduce steric hindrance or eliminate critical interactions, as expected, but that they may also have indirect effects such as altering protein dynamics and enzyme kinetics. The structure-based development of the latest generation of inhibitors targeting HIV reverse transcriptase, P. falciparum and S. aureus dihydrofolate reductase, neuraminidase, and epithelial growth factor receptor (EGFR) tyrosine kinase, is highlighted to emphasize lessons that may be applied to future drug discovery to overcome mutation-induced resistance. Successful next-generation drugs tend to be more flexible and exploit a greater number of interactions mimicking those of the substrate with conserved residues.

Chemical chaperones assist intracellular folding to buffer mutational variations

Hidden genetic variations harbor potential for the evolution of new traits. Molecular chaperones, that assist protein folding, may conceal genetic variations in protein coding regions. Here, we investigate if the chemical milieu of cells has the potential to alleviate intracellular protein folding; potentially implicating a role of osmolytes in concealing genetic variations. Using the model osmolyte TMAO, we uncover that it can buffer mutations that impose kinetic traps in the folding pathways of two model proteins. Using this information, we rationally designed TMAO-dependent mutants in vivo, starting from a TMAO-independent protein. Strikingly, we delineate different osmolytes to have a unique spectrum of buffered-mutations. Consequently, the chemical milieu of cells may alter the folding pathways of unique mutant variants in polymorphic populations and lead to unanticipated spectra of genetic buffering.

Experimental approaches to evaluate the contributions of candidate protein-coding mutations to phenotypic evolution

Identifying mechanisms of molecular adaptation can provide important insights into the process of phenotypic evolution, but it can be exceedingly difficult to quantify the phenotypic effects of specific mutational changes. To verify the adaptive significance of genetically based changes in protein function, it is necessary to document functional differences between the products of derived and wild-type alleles and to demonstrate that such differences impinge on higher-level physiological processes (and ultimately, fitness). In the case of metabolic enzymes, this requires documenting in vivo differences in reaction rate that give rise to differences in flux through the pathway in which the enzymes function. These measured differences in pathway flux should then give rise to differences in cellular or systemic physiology that affect fitness-related variation in whole-organism performance. Efforts to establish these causal connections between genotype, phenotype, and fitness require experiments that carefully control for environmental variation and background genetic variation. Here, we discuss experimental approaches to evaluate the contributions of amino-acid mutations to adaptive phenotypic change. We discuss conceptual and methodological issues associated with in vitro and in vivo studies of protein function, and the evolutionary insights that can be gleaned from such studies. We also discuss the importance of isolating the effects of individual mutations to distinguish between positively selected substitutions that directly contribute to improvements in protein function versus positively selected, compensatory substitutions that mitigate negative pleiotropic effects of antecedent changes.

The transition of human estrogen sulfotransferase from generalist to specialist using directed enzyme evolution

Broad specificity is believed to be a property of primordial enzymes that diverged during natural protein evolution to produce highly specific and efficient enzymes. Human estrogen sulfotransferase (SULT1E1) is a broad-specificity enzyme that detoxifies a variety of chemicals, including estrogens, by the transfer of sulfate. To study the molecular basis for the broad specificity of this enzyme and to investigate the process of SULT1E1 specialization, we have adopted a directed enzyme evolution approach. Using two iterative rounds of evolution, we generated SULT1E1 mutants with increased thermostability and narrower specificity from the broadly specific wild-type enzyme. To identify mutants with enhanced specificity, we developed an unbiased screening assay to assess sulfate transfer to three different acceptors in parallel. Such an assay enabled the isolation of SULT1E1 mutants with enhanced or wild-type activity toward an estrogen acceptor and significantly reduced activity for phenol or coumarin type of acceptors, leading to up to 3 orders of magnitude increase in specificity. We found that mutations conferring novel specificity are located in the vicinity of the active site and thus may play a direct role in reshaping the acceptor-binding site. Finally, such mutations resulted in reduced SULT1E1 thermostability, revealing a trade-off between SULT1E1 thermostability and acquisition of novel function.

Molecular analyses of TEM genes and their corresponding penicillinase-producing Neisseria gonorrhoeae isolates in Bangkok, Thailand

Neisseria gonorrhoeae is a major public health problem globally, especially because the bacterium has developed resistance to most antimicrobials introduced for first-line treatment of gonorrhea. In the present study, 96 N. gonorrhoeae isolates with high-level resistance to penicillin from 121 clinical isolates in Thailand were examined to investigate changes related to their plasmid-mediated penicillin resistance and their molecular epidemiological relationships. A β-lactamase (TEM) gene variant, bla(TEM-135), that may be a precursor in the transitional stage of a traditional bla(TEM-1) gene into an extended-spectrum β-lactamase (ESBL), possibly causing high resistance to all extended-spectrum cephalosporins in N. gonorrhoeae, was identified. Clonal analysis using multilocus sequence typing (MLST) and N. gonorrhoeae multiantigen sequence typing (NG-MAST) revealed the existence of a sexual network among patients from Japan and Thailand. Molecular analysis of the bla(TEM-135) gene showed that the emergence of this allele might not be a rare genetic event and that the allele has evolved in different plasmid backgrounds, which results possibly indicate that it is selected due to antimicrobial pressure. The presence of the bla(TEM-135) allele in the penicillinase-producing N. gonorrhoeae population may call for monitoring for the possible emergence of ESBL-producing N. gonorrhoeae in the future. This study identified a bla(TEM) variant (bla(TEM-135)) that is a possible intermediate precursor for an ESBL, which warrants international awareness.

Directed evolution to re-adapt a co-evolved network within an enzyme

We have previously used targeted active-site saturation mutagenesis to identify a number of transketolase single mutants that improved activity towards either glycolaldehyde (GA), or the non-natural substrate propionaldehyde (PA). Here, all attempts to recombine the singles into double mutants led to unexpected losses of specific activity towards both substrates. A typical trade-off occurred between soluble expression levels and specific activity for all single mutants, but many double mutants decreased both properties more severely suggesting a critical loss of protein stability or native folding. Statistical coupling analysis (SCA) of a large multiple sequence alignment revealed a network of nine co-evolved residues that affected all but one double mutant. Such networks maintain important functional properties such as activity, specificity, folding, stability, and solubility and may be rapidly disrupted by introducing one or more non-naturally occurring mutations. To identify variants of this network that would accept and improve upon our best D469 mutants for activity towards PA, we created a library of random single, double and triple mutants across seven of the co-evolved residues, combining our D469 variants with only naturally occurring mutations at the remaining sites. A triple mutant cluster at D469, E498 and R520 was found to behave synergistically for the specific activity towards PA. Protein expression was severely reduced by E498D and improved by R520Q, yet variants containing both mutations led to improved specific activity and enzyme expression, but with loss of solubility and the formation of inclusion bodies. D469S and R520Q combined synergistically to improve k_cat 20-fold for PA, more than for any previous transketolase mutant. R520Q also doubled the specific activity of the previously identified D469T to create our most active transketolase mutant to date. Our results show that recombining active-site mutants obtained by saturation mutagenesis can rapidly destabilise critical networks of co-evolved residues, whereas beneficial single mutants can be retained and improved upon by randomly recombining them with natural variants at other positions in the network.

The molecular basis for the broad substrate specificity of human sulfotransferase 1A1

Cytosolic sulfotransferases (SULTs) are mammalian enzymes that detoxify a wide variety of chemicals through the addition of a sulfate group. Despite extensive research, the molecular basis for the broad specificity of SULTs is still not understood. Here, structural, protein engineering and kinetic approaches were employed to obtain deep understanding of the molecular basis for the broad specificity, catalytic activity and substrate inhibition of SULT1A1. We have determined five new structures of SULT1A1 in complex with different acceptors, and utilized a directed evolution approach to generate SULT1A1 mutants with enhanced thermostability and increased catalytic activity. We found that active site plasticity enables binding of different acceptors and identified dramatic structural changes in the SULT1A1 active site leading to the binding of a second acceptor molecule in a conserved yet non-productive manner. Our combined approach highlights the dominant role of SULT1A1 structural flexibility in controlling the specificity and activity of this enzyme.

The linear interaction energy method for the prediction of protein stability changes upon mutation

The coupling of protein energetics and sequence changes is a critical aspect of computational protein design, as well as for the understanding of protein evolution, human disease, and drug resistance. To study the molecular basis for this coupling, computational tools must be sufficiently accurate and computationally inexpensive enough to handle large amounts of sequence data. We have developed a computational approach based on the linear interaction energy (LIE) approximation to predict the changes in the free-energy of the native state induced by a single mutation. This approach was applied to a set of 822 mutations in 10 proteins which resulted in an average unsigned error of 0.82 kcal/mol and a correlation coefficient of 0.72 between the calculated and experimental ΔΔG values. The method is able to accurately identify destabilizing hot spot mutations; however, it has difficulty in distinguishing between stabilizing and destabilizing mutations because of the distribution of stability changes for the set of mutations used to parameterize the model. In addition, the model also performs quite well in initial tests on a small set of double mutations. On the basis of these promising results, we can begin to examine the relationship between protein stability and fitness, correlated mutations, and drug resistance.

The causes of epistasis

Since Bateson's discovery that genes can suppress the phenotypic effects of other genes, gene interactions-called epistasis-have been the topic of a vast research effort. Systems and developmental biologists study epistasis to understand the genotype-phenotype map, whereas evolutionary biologists recognize the fundamental importance of epistasis for evolution. Depending on its form, epistasis may lead to divergence and speciation, provide evolutionary benefits to sex and affect the robustness and evolvability of organisms. That epistasis can itself be shaped by evolution has only recently been realized. Here, we review the empirical pattern of epistasis, and some of the factors that may affect the form and extent of epistasis. Based on their divergent consequences, we distinguish between interactions with or without mean effect, and those affecting the magnitude of fitness effects or their sign. Empirical work has begun to quantify epistasis in multiple dimensions in the context of metabolic and fitness landscape models. We discuss possible proximate causes (such as protein function and metabolic networks) and ultimate factors (including mutation, recombination, and the importance of natural selection and genetic drift). We conclude that, in general, pleiotropy is an important prerequisite for epistasis, and that epistasis may evolve as an adaptive or intrinsic consequence of changes in genetic robustness and evolvability.

First-step mutations for adaptation at elevated temperature increase capsid stability in a virus

The relationship between mutation, protein stability and protein function plays a central role in molecular evolution. Mutations tend to be destabilizing, including those that would confer novel functions such as host-switching or antibiotic resistance. Elevated temperature may play an important role in preadapting a protein for such novel functions by selecting for stabilizing mutations. In this study, we test the stability change conferred by single mutations that arise in a G4-like bacteriophage adapting to elevated temperature. The vast majority of these mutations map to interfaces between viral coat proteins, suggesting they affect protein-protein interactions. We assess their effects by estimating thermodynamic stability using molecular dynamic simulations and measuring kinetic stability using experimental decay assays. The results indicate that most, though not all, of the observed mutations are stabilizing.

Network models of TEM β-lactamase mutations coevolving under antibiotic selection show modular structure and anticipate evolutionary trajectories

Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of β-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, β-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.

Understanding how new biological activities evolve on the molecular level has critical implications for biotechnology and for human health. Here we collect a database of mutations that contribute to the evolution of β-lactamase resistance to inhibitors and to new β-lactam antibiotics in bacterial pathogens, such as Escherichia coli. We compiled a database of TEM β-lactamase sequences evolved under antibiotic pressure and identified functional interactions between individual residue positions. We visualized these complex molecular interactions as a network and used network theory to derive information regarding the origin of individual mutations and their contribution to the observed resistance. Our approach should help interpret sequence databases for clinically relevant proteins undergoing high mutation rates and under selective (drug, immune) pressure, such as surface proteins of pathogens (particularly of RNA viruses such as HIV) or targets for chemotherapy in microbial pathogen or tumor cells. Notably, our approach only requires sequence data; detailed phylogenetic or tertiary structure information for the target gene is not necessary. Our analysis of how individual mutations work together to produce new biological activities should help anticipate evolution driven by a variety of clinically-relevant selections such as drug resistance, virulence, and immunity.

Accessory mutations maintain stability in drug-resistant HIV-1 protease

The underlying mechanisms driving the evolution of drug resistance in human immunodeficiency virus (HIV) are only partially understood. We investigated the evolutionary cost of the major resistance mutations in HIV-1 protease in terms of protein stability. The accumulation of resistance mutations destabilizes the protease, limiting further adaptation. From an analysis of clinical isolates, we identified specific accessory mutations that were able to restore the stability of the protease or even increase it beyond the wild-type baseline. Resistance mutations were also found to decrease the activity of HIV protease near neutral pH values, while incorporating stabilizing mutations improved the enzyme's pH tolerance. These findings help us to explain the prevalence of mutations located far from the active site and underscore the importance of protein stability during the evolution of drug resistance in HIV.

Oseltamivir-resistant influenza viruses get by with a little help from permissive mutations

Influenza A viruses (IAVs) encode two critical glycoproteins, hemagglutinin and neuraminidase (NA). Hemagglutinin promotes viral docking onto cells via interactions with IAV's receptor, sialic acid and NA facilitates release of newly synthesized virions by cleaving cellular and viral sialic acid. NA inhibitors, such as oseltamivir, are widely used drugs that work by binding to the active site of NA. Although oseltamivir-resistant viruses were easily generated years ago in laboratory experiments, it was widely believed that these viruses would not be able to circulate in the human population as they did not replicate efficiently. However, oseltamivir-resistant H1N1 viruses rapidly spread during the 2007-2008 IAV season and these viruses contained precisely the same exact drug-resistance mutation identified years prior, a histidine to tyrosine substitution at NA residue 274 (H274Y). Unlike the experimentally derived NA inhibitor-resistant viruses, 2007-2008 H1N1 viruses containing H274Y replicated efficiently. Bloom et al. have solved this riddle by identifying permissive NA mutations that allow viruses to tolerate H274Y. Here, we discuss these important findings and speculate how these studies may facilitate early detection of drug-resistant strains in the future.

Determinants of homodimerization specificity in histidine kinases

Two-component signal transduction pathways consisting of a histidine kinase and a response regulator are used by prokaryotes to respond to diverse environmental and intracellular stimuli. Most species encode numerous paralogous histidine kinases that exhibit significant structural similarity. Yet in almost all known examples, histidine kinases are thought to function as homodimers. We investigated the molecular basis of dimerization specificity, focusing on the model histidine kinase EnvZ and RstB, its closest paralog in Escherichia coli. Direct binding studies showed that the cytoplasmic domains of these proteins each form specific homodimers in vitro. Using a series of chimeric proteins, we identified specificity determinants at the base of the four-helix bundle in the dimerization and histidine phosphotransfer domain. Guided by molecular coevolution predictions and EnvZ structural information, we identified sets of residues in this region that are sufficient to establish homospecificity. Mutating these residues in EnvZ to the corresponding residues in RstB produced a functional kinase that preferentially homodimerized over interacting with EnvZ. EnvZ and RstB likely diverged following gene duplication to yield two homodimers that cannot heterodimerize, and the mutants we identified represent possible evolutionary intermediates in this process.

A computational-experimental approach identifies mutations that enhance surface expression of an oseltamivir-resistant influenza neuraminidase

The His274→Tyr (H274Y) oseltamivir (Tamiflu) resistance mutation causes a substantial decrease in the total levels of surface-expressed neuraminidase protein and activity in early isolates of human seasonal H1N1 influenza, and in the swine-origin pandemic H1N1. In seasonal H1N1, H274Y only became widespread after the occurrence of secondary mutations that counteracted this decrease. H274Y is currently rare in pandemic H1N1, and it remains unclear whether secondary mutations exist that might similarly counteract the decreased neuraminidase surface expression associated with this resistance mutation in pandemic H1N1. Here we investigate the possibility of predicting such secondary mutations. We first test the ability of several computational approaches to retrospectively identify the secondary mutations that enhanced levels of surface-expressed neuraminidase protein and activity in seasonal H1N1 shortly before the emergence of oseltamivir resistance. We then use the most successful computational approach to predict a set of candidate secondary mutations to the pandemic H1N1 neuraminidase. We experimentally screen these mutations, and find that several of them do indeed partially counteract the decrease in neuraminidase surface expression caused by H274Y. Two of the secondary mutations together restore surface-expressed neuraminidase activity to wildtype levels, and also eliminate the very slight decrease in viral growth in tissue-culture caused by H274Y. Our work therefore demonstrates a combined computational-experimental approach for identifying mutations that enhance neuraminidase surface expression, and describes several specific mutations with the potential to be of relevance to the spread of oseltamivir resistance in pandemic H1N1.

Resistance to the third-generation cephalosporin ceftazidime by a deacylation-deficient mutant of the TEM β-lactamase by the uncommon covalent-trapping mechanism

The Glu166Arg/Met182Thr mutant of Escherichia coli TEM(pTZ19-3) β-lactamase produces a 128-fold increase in the level of resistance to the antibiotic ceftazidime in comparison to that of the parental wild-type enzyme. The single Glu166Arg mutation resulted in a dramatic decrease in both the level of enzyme expression in bacteria and the resistance to penicillins, with a concomitant 4-fold increase in the resistance to ceftazidime, a third-generation cephalosporin. Introduction of the second amino acid substitution, Met182Thr, restored enzyme expression to a level comparable to that of the wild-type enzyme and resulted in an additional 32-fold increase in the minimal inhibitory concentration of ceftazidime to 64 μg/mL. The double mutant formed a stable covalent complex with ceftazidime that remained intact for the entire duration of the monitoring, which exceeded a time period of 40 bacterial generations. Compared to those of the wild-type enzyme, the affinity of the TEM(pTZ19-3) Glu166Arg/Met182Thr mutant for ceftazidime increased by at least 110-fold and the acylation rate constant was augmented by at least 16-fold. The collective experimental data and computer modeling indicate that the deacylation-deficient Glu166Arg/Met182Thr mutant of TEM(pTZ19-3) produces resistance to the third-generation cephalosporin ceftazidime by an uncommon covalent-trapping mechanism. This is the first documentation of such a mechanism by a class A β-lactamase in a manifestation of resistance.

Molecular mechanisms of epistasis within and between genes

'Disease-causing' mutations do not cause disease in all individuals. One possible important reason for this is that the outcome of a mutation can depend upon other genetic variants in a genome. These epistatic interactions between mutations occur both within and between molecules, and studies in model organisms show that they are extremely prevalent. However, epistatic interactions are still poorly understood at the molecular level, and consequently difficult to predict de novo. Here I provide an overview of our current understanding of the molecular mechanisms that can cause epistasis, and areas where more research is needed. A more complete understanding of epistasis will be vital for making accurate predictions about the phenotypes of individuals.

Exploring sequence requirements for C₃/C₄ carboxylate recognition in the Pseudomonas aeruginosa cephalosporinase: Insights into plasticity of the AmpC β-lactamase

In Pseudomonas aeruginosa, the chromosomally encoded class C cephalosporinase (AmpC β-lactamase) is often responsible for high-level resistance to β-lactam antibiotics. Despite years of study of these important β-lactamases, knowledge regarding how amino acid sequence dictates function of the AmpC Pseudomonas-derived cephalosporinase (PDC) remains scarce. Insights into structure-function relationships are crucial to the design of both β-lactams and high-affinity inhibitors. In order to understand how PDC recognizes the C₃/C₄ carboxylate of β-lactams, we first examined a molecular model of a P. aeruginosa AmpC β-lactamase, PDC-3, in complex with a boronate inhibitor that possesses a side chain that mimics the thiazolidine/dihydrothiazine ring and the C₃/C₄ carboxylate characteristic of β-lactam substrates. We next tested the hypothesis generated by our model, i.e. that more than one amino acid residue is involved in recognition of the C₃/C₄ β-lactam carboxylate, and engineered alanine variants at three putative carboxylate binding amino acids. Antimicrobial susceptibility testing showed that the PDC-3 β-lactamase maintains a high level of activity despite the substitution of C₃/C₄ β-lactam carboxylate recognition residues. Enzyme kinetics were determined for a panel of nine penicillin and cephalosporin analog boronates synthesized as active site probes of the PDC-3 enzyme and the Arg349Ala variant. Our examination of the PDC-3 active site revealed that more than one residue could serve to interact with the C₃/C₄ carboxylate of the β-lactam. This functional versatility has implications for novel drug design, protein evolution, and resistance profile of this enzyme.

Linking genotype and phenotype of Saccharomyces cerevisiae strains reveals metabolic engineering targets and leads to triterpene hyper-producers

BACKGROUND

Metabolic engineering is an attractive approach in order to improve the microbial production of drugs. Triterpenes is a chemically diverse class of compounds and many among them are of interest from a human health perspective. A systematic experimental or computational survey of all feasible gene modifications to determine the genotype yielding the optimal triterpene production phenotype is a laborious and time-consuming process.

METHODOLOGY/PRINCIPAL FINDINGS

Based on the recent genome-wide sequencing of Saccharomyces cerevisiae CEN.PK 113-7D and its phenotypic differences with the S288C strain, we implemented a strategy for the construction of a β-amyrin production platform. The genes Erg8, Erg9 and HFA1 contained non-silent SNPs that were computationally analyzed to evaluate the changes that cause in the respective protein structures. Subsequently, Erg8, Erg9 and HFA1 were correlated with the increased levels of ergosterol and fatty acids in CEN.PK 113-7D and single, double, and triple gene over-expression strains were constructed.

CONCLUSIONS

The six out of seven gene over-expression constructs had a considerable impact on both ergosterol and β-amyrin production. In the case of β-amyrin formation the triple over-expression construct exhibited a nearly 500% increase over the control strain making our metabolic engineering strategy the most successful design of triterpene microbial producers.

The contribution of coevolving residues to the stability of KDO8P synthase

Background

The evolutionary tree of 3-deoxy-D-manno-octulosonate 8-phosphate (KDO8P) synthase (KDO8PS), a bacterial enzyme that catalyzes a key step in the biosynthesis of bacterial endotoxin, is evenly divided between metal and non-metal forms, both having similar structures, but diverging in various degrees in amino acid sequence. Mutagenesis, crystallographic and computational studies have established that only a few residues determine whether or not KDO8PS requires a metal for function. The remaining divergence in the amino acid sequence of KDO8PSs is apparently unrelated to the underlying catalytic mechanism.

Methodology/Principal Findings

The multiple alignment of all known KDO8PS sequences reveals that several residue pairs coevolved, an indication of their possible linkage to a structural constraint. In this study we investigated by computational means the contribution of coevolving residues to the stability of KDO8PS. We found that about 1/4 of all strongly coevolving pairs probably originated from cycles of mutation (decreasing stability) and suppression (restoring it), while the remaining pairs are best explained by a succession of neutral or nearly neutral covarions.

Conclusions/Significance

Both sequence conservation and coevolution are involved in the preservation of the core structure of KDO8PS, but the contribution of coevolving residues is, in proportion, smaller. This is because small stability gains or losses associated with selection of certain residues in some regions of the stability landscape of KDO8PS are easily offset by a large number of possible changes in other regions. While this effect increases the tolerance of KDO8PS to deleterious mutations, it also decreases the probability that specific pairs of residues could have a strong contribution to the thermodynamic stability of the protein.

Structure of the extended-spectrum β-lactamase TEM-72 inhibited by citrate

TEM-72, a class A β-lactamase identified in isolates of Enterobacteriaceae, is a quadruple mutant of TEM-1 (Q39K, M182T, G238S and E240K) and shows extended-spectrum β-lactamase (ESBL) properties arising from the G238S and E240K substitutions. Although many structures of TEM variants have been published, they do not include an enzyme with the simultaneous presence of both of the ESBL-conferring G238S and E240K substitutions. Furthermore, the structure shows the presence of a citrate anion bound to the TEM-72 active site, where it interacts with all of the conserved residues of class A β-lactamases. The present structure supports the use of polycarboxylates as a scaffold for the design of broad-spectrum inhibitors of serine β-lactamases.

Ligand-dependent disorder of the Omega loop observed in extended-spectrum SHV-type beta-lactamase

Among Gram-negative bacteria, resistance to β-lactams is mediated primarily by β-lactamases (EC 3.2.6.5), periplasmic enzymes that inactivate β-lactam antibiotics. Substitutions at critical amino acid positions in the class A β-lactamase families result in enzymes that can hydrolyze extended-spectrum cephalosporins, thus demonstrating an "extended-spectrum" β-lactamase (ESBL) phenotype. Using SHV ESBLs with substitutions in the Ω loop (R164H and R164S) as target enzymes to understand this enhanced biochemical capability and to serve as a basis for novel β-lactamase inhibitor development, we determined the spectra of activity and crystal structures of these variants. We also studied the inactivation of the R164H and R164S mutants with tazobactam and SA2-13, a unique β-lactamase inhibitor that undergoes a distinctive reaction chemistry in the active site. We noted that the reduced Ki values for the R164H and R164S mutants with SA2-13 are comparable to those with tazobactam (submicromolar). The apo enzyme crystal structures of the R164H and R164S SHV variants revealed an ordered Ω loop architecture that became disordered when SA2-13 was bound. Important structural alterations that result from the binding of SA2-13 explain the enhanced susceptibility of these ESBL enzymes to this inhibitor and highlight ligand-dependent Ω loop flexibility as a mechanism for accommodating and hydrolyzing β-lactam substrates.

Natural evolution of TEM-1 β-lactamase: experimental reconstruction and clinical relevance

TEM-1 β-lactamase is one of the most well-known antibiotic resistance determinants around. It confers resistance to penicillins and early cephalosporins and has shown an astonishing functional plasticity in response to the introduction of novel drugs derived from these antibiotics. Since its discovery in the 1960s, over 170 variants of TEM-1 - with different amino acid sequences and often resistance phenotypes - have been isolated in hospitals and clinics worldwide. Next to this well-documented 'natural' evolution, the in vitro evolution of TEM-1 has been the focus of attention of many experimental studies. In this review, we compare the natural and laboratory evolution of TEM-1 in order to address the question to what extent the evolution of antibiotic resistance can be repeated, and hence might have been predicted, under laboratory conditions. We also use the comparison to gain an insight into the adaptive relevance of hitherto uncharacterized substitutions present in clinical isolates and to predict substitutions not yet observed in nature. Based on new structural insights, we review what is known about substitutions in TEM-1 that contribute to the extension of its resistance phenotype. Finally, we address the clinical relevance of TEM alleles during the past decade, which has been dominated by the emergence of another β-lactamase, CTX-M.

Multiple global suppressors of protein stability defects facilitate the evolution of extended-spectrum TEM β-lactamases

The introduction of extended-spectrum cephalosporins and β-lactamase inhibitors has driven the evolution of extended-spectrum β-lactamases (ESBLs) that possess the ability to hydrolyze these drugs. The evolved TEM ESBLs from clinical isolates of bacteria often contain substitutions that occur in the active site and alter the catalytic properties of the enzyme to provide an increased hydrolysis of extended-spectrum cephalosporins or an increased resistance to inhibitors. These active-site substitutions often result in a cost in the form of reduced enzyme stability. The evolution of TEM ESBLs is facilitated by mutations that act as global suppressors of protein stability defects in that they allow the enzyme to absorb multiple amino acid changes despite incremental losses in stability associated with the substitutions. The best-studied example is the M182T substitution, which corrects protein stability defects and is commonly found in TEM ESBLs or inhibitor-resistant variants from clinical isolates. In this study, a genetic selection for second-site mutations that could partially restore function to a severely destabilized primary mutant enabled the identification of A184V, T265M, R275Q, and N276D, which are known to occur in TEM ESBLs from clinical isolates, as suppressors of TEM-1 protein stability defects. Further characterization demonstrated that these substitutions increased the thermal stability of TEM-1 and were able to correct the stability defects of two different sets of destabilizing mutations. The acquisition of compensatory global suppressors of stability costs associated with active-site mutations may be a common mechanism for the evolution of novel protein function.

Mutational effects and the evolution of new protein functions

The divergence of new genes and proteins occurs through mutations that modulate protein function. However, mutations are pleiotropic and can have different effects on organismal fitness depending on the environment, as well as opposite effects on protein function and dosage. We review the pleiotropic effects of mutations. We discuss how they affect the evolution of gene and protein function, and how these complex mutational effects dictate the likelihood and mechanism of gene duplication and divergence. We propose several factors that can affect the divergence of new protein functions, including mutational trade-offs and hidden, or apparently neutral, variation.

The effect of chaperonin buffering on protein evolution

Molecular chaperones are highly conserved and ubiquitous proteins that help other proteins in the cell to fold. Pioneering work by Rutherford and Lindquist suggested that the chaperone Hsp90 could buffer (i.e., suppress) phenotypic variation in its client proteins and that alternate periods of buffering and expression of these variants might be important in adaptive evolution. More recently, Tokuriki and Tawfik presented an explicit mechanism for chaperone-dependent evolution, in which the Escherichia coli chaperonin GroEL facilitated the folding of clients that had accumulated structurally destabilizing but neofunctionalizing mutations in the protein core. But how important an evolutionary force is chaperonin-mediated buffering in nature? Here, we address this question by modeling the per-residue evolutionary rate of the crystallized E. coli proteome, evaluating the relative contributions of chaperonin buffering, functional importance, and structural features such as residue contact density. Previous findings suggest an interaction between codon bias and GroEL in limiting the effects of misfolding errors. Our results suggest that the buffering of deleterious mutations by GroEL increases the evolutionary rate of client proteins. We then examine the evolutionary fate of GroEL clients in the Mycoplasmas, a group of bacteria containing the only known organisms that lack chaperonins. We show that GroEL was lost once in the common ancestor of a monophyletic subgroup of Mycoplasmas, and we evaluate the effect of this loss on the subsequent evolution of client proteins, providing evidence that client homologs in 11 Mycoplasma species have lost their obligate dependency on GroEL for folding. Our analyses indicate that individual molecules such as chaperonins can have significant effects on proteome evolution through their modulation of protein folding.

Role of conservative mutations in protein multi-property adaptation

Protein physicochemical properties must undergo complex changes during evolution, as a response to modifications in the organism environment, the result of the proteins taking up new roles or because of the need to cope with the evolution of molecular interacting partners. Recent work has emphasized the role of stability and stability-function trade-offs in these protein adaptation processes. In the present study, on the other hand, we report that combinations of a few conservative, high-frequency-of-fixation mutations in the thioredoxin molecule lead to largely independent changes in both stability and the diversity of catalytic mechanisms, as revealed by single-molecule atomic force spectroscopy. Furthermore, the changes found are evolutionarily significant, as they combine typically hyperthermophilic stability enhancements with modulations in function that span the ranges defined by the quite different catalytic patterns of thioredoxins from bacterial and eukaryotic origin. These results suggest that evolutionary protein adaptation may use, in some cases at least, the potential of conservative mutations to originate a multiplicity of evolutionarily allowed mutational paths leading to a variety of protein modulation patterns. In addition the results support the feasibility of using evolutionary information to achieve protein multi-feature optimization, an important biotechnological goal.

Enzyme promiscuity: a mechanistic and evolutionary perspective

Many, if not most, enzymes can promiscuously catalyze reactions, or act on substrates, other than those for which they evolved. Here, we discuss the structural, mechanistic, and evolutionary implications of this manifestation of infidelity of molecular recognition. We define promiscuity and related phenomena and also address their generality and physiological implications. We discuss the mechanistic enzymology of promiscuity--how enzymes, which generally exert exquisite specificity, catalyze other, and sometimes barely related, reactions. Finally, we address the hypothesis that promiscuous enzymatic activities serve as evolutionary starting points and highlight the unique evolutionary features of promiscuous enzyme functions.

Permissive secondary mutations enable the evolution of influenza oseltamivir resistance

The His274-->Tyr274 (H274Y) mutation confers oseltamivir resistance on N1 influenza neuraminidase but had long been thought to compromise viral fitness. However, beginning in 2007-2008, viruses containing H274Y rapidly became predominant among human seasonal H1N1 isolates. We show that H274Y decreases the amount of neuraminidase that reaches the cell surface and that this defect can be counteracted by secondary mutations that also restore viral fitness. Two such mutations occurred in seasonal H1N1 shortly before the widespread appearance of H274Y. The evolution of oseltamivir resistance was therefore enabled by "permissive" mutations that allowed the virus to tolerate subsequent occurrences of H274Y. An understanding of this process may provide a basis for predicting the evolution of oseltamivir resistance in other influenza strains.

The population genetics of antibiotic resistance: integrating molecular mechanisms and treatment contexts

Despite efforts from a range of disciplines, our ability to predict and combat the evolution of antibiotic resistance in pathogenic bacteria is limited. This is because resistance evolution involves a complex interplay between the specific drug, bacterial genetics and both natural and treatment ecology. Incorporating details of the molecular mechanisms of drug resistance and ecology into evolutionary models has proved useful in predicting the dynamics of resistance evolution. However, putting these models to practical use will require extensive collaboration between mathematicians, molecular biologists, evolutionary ecologists and clinicians.

Structural insights into substrate recognition and product expulsion in CTX-M enzymes

beta-Lactamase-mediated resistance to beta-lactam antibiotics poses a major threat to our antibiotic armamentarium. Among beta-lactamases, a significant threat comes from enzymes that hydrolyze extended-spectrum cephalosporins such as cefotaxime. Among the enzymes that exhibit this phenotype, the CTX-M family is found worldwide. These enzymes have a small active site, which makes it difficult to explain how they hydrolyze the bulky extended-spectrum cephalosporins into the binding site. We investigated noncovalent substrate recognition and product release in CTX-M enzymes using steered molecular dynamics simulation and X-ray diffraction. An arginine residue located far from the binding site favors the capture and tracking of substrates during entrance into the catalytic pocket. We show that the accommodation of extended-spectrum cephalosporins by CTX-M enzymes induced subtle changes in the active site and established a high density of electrostatic interactions. Interestingly, the product of the catalytic reaction initiates its own release because of steric hindrances and electrostatic repulsions. This suggests that there exists a general mechanism for product release for all members of the beta-lactamase family and probably for most carboxypeptidases.