A secondary drug resistance mutation of TEM-1 beta-lactamase that suppresses misfolding and aggregation

Simultaneous enhancement of multiple functional properties using evolution-informed protein design

A major challenge in protein design is to augment existing functional proteins with multiple property enhancements. Altering several properties likely necessitates numerous primary sequence changes, and novel methods are needed to accurately predict combinations of mutations that maintain or enhance function. Models of sequence co-variation (e.g., EVcouplings), which leverage extensive information about various protein properties and activities from homologous protein sequences, have proven effective for many applications including structure determination and mutation effect prediction. We apply EVcouplings to computationally design variants of the model protein TEM-1 β-lactamase. Nearly all the 14 experimentally characterized designs were functional, including one with 84 mutations from the nearest natural homolog. The designs also had large increases in thermostability, increased activity on multiple substrates, and nearly identical structure to the wild type enzyme. This study highlights the efficacy of evolutionary models in guiding large sequence alterations to generate functional diversity for protein design applications.

Genetic diversity and prevalence of extended-spectrum beta-lactamase-producing Escherichia coli and Klebsiella pneumoniae in aquatic environments receiving untreated hospital effluents

The spread of extended-spectrum beta-lactamase (ESBL)-producing bacteria in the environment has been recognized as a challenge to public health. The aim of the present study was to assess the occurrence of ESBL-producing Escherichia coli and Klebsiella pneumoniae from selected water bodies receiving hospital effluents in Kerala, India. Nearly 69.8% of Enterobacteriaceae isolates were multi-drug resistant by the Kirby-Bauer disc diffusion method. The double disc synergy test was used to detect the ESBL production and the genes responsible for imparting resistance were detected by PCR. Conjugation experiments confirmed the mechanism of plasmid-mediated transfer of resistance. The prevalence of ESBL production in E. coli and K. pneumoniae was 49.2 and 46.8%, respectively. Among the ESBL-encoding genes, bla_CTX-M was the most prevalent group followed by bla_TEM, bla_OXA, bla_CMY, and bla_SHV. The results suggest that healthcare settings are one of the key contributors to the spread of ESBL-producing bacteria, not only through cross-transmission and ingestion of antibiotics but also through the discharge of waste without a proper treatment, leading to harmful effects on the aquatic environment. The high prevalence of ESBL-producing Enterobacteriaceae with resistance genes in public water bodies even post-treatment poses a serious threat.

Local and Global Protein Interactions Contribute to Residue Entrenchment in Beta-Lactamase TEM-1

Due to their rapid evolution and their impact on healthcare, beta-lactamases, protein degrading beta-lactam antibiotics, are used as generic models of protein evolution. Therefore, we investigated the mutation effects in two distant beta-lactamases, TEM-1 and CTX-M-15. Interestingly, we found a site with a complex pattern of genetic interactions. Mutation G251W in TEM-1 inactivates the protein’s function, just as the reciprocal mutation, W251G, does in CTX-M-15. The phylogenetic analysis revealed that mutation G has been entrenched in TEM-1’s background: while rarely observed throughout the phylogeny, it is essential in TEM-1. Using a rescue experiment, in the TEM-1 G251W mutant, we identified sites that alleviate the deviation from G to W. While few of these mutations could potentially involve local interactions, most of them were found on distant residues in the 3D structure. Many well-known mutations that have an impact on protein stability, such as M182T, were recovered. Our results therefore suggest that entrenchment of an amino acid may rely on diffuse interactions among multiple sites, with a major impact on protein stability.

Structure-Guided Engineering of a Family IV Cold-Adapted Esterase Expands Its Substrate Range

Cold active esterases have gained great interest in several industries. The recently determined structure of a family IV cold active esterase (EstN7) from Bacillus cohnii strain N1 was used to expand its substrate range and to probe its commercially valuable substrates. Database mining suggested that triacetin was a potential commercially valuable substrate for EstN7, which was subsequently proved experimentally with the final product being a single isomeric product, 1,2-glyceryl diacetate. Enzyme kinetics revealed that EstN7’s activity is restricted to C2 and C4 substrates due to a plug at the end of the acyl binding pocket that blocks access to a buried water-filled cavity. Residues M187, N211 and W206 were identified as key plug forming residues. N211A stabilised EstN7 allowing incorporation of the destabilising M187A mutation. The M187A-N211A double mutant had the broadest substrate range, capable of hydrolysing a C8 substrate. W206A did not appear to have any significant effect on substrate range either alone or when combined with the double mutant. Thus, the enzyme kinetics and engineering together with a recently determined structure of EstN7 provide new insights into substrate specificity and the role of acyl binding pocket plug residues in determining family IV esterase stability and substrate range.

Assessment of Phenotype Relevant Amino Acid Residues in TEM-β-Lactamases by Mathematical Modelling and Experimental Approval

Single substitutions or combinations of them alter the hydrolytic activity towards specific β-lactam-antibiotics and β-lactamase inhibitors of TEM-β-lactamases. The sequences and phenotypic classification of allelic TEM variants, as provided by the NCBI National Database of Antibiotic Resistant Organisms, does not attribute phenotypes to all variants. Some entries are doubtful as the data assessment differs strongly between the studies or no data on the methodology are provided at all. This complicates mathematical and bioinformatic predictions of phenotypes that rely on the database. The present work aimed to prove the role of specific substitutions on the resistance phenotype of TEM variants in, to our knowledge, the most extensive mutagenesis study. In parallel, the predictive power of extrapolation algorithms was assessed. Most well-known substitutions with direct impact on the phenotype could be reproduced, both mathematically and experimentally. Most discrepancies were found for supportive substitutions, where some resulted in antagonistic effects in contrast to previously described synergism. The mathematical modelling proved to predict the strongest phenotype-relevant substitutions accurately but showed difficulties in identifying less prevalent but still phenotype transforming ones. In general, mutations increasing cephalosporin resistance resulted in increased sensitivity to β-lactamase inhibitors and vice versa. Combining substitutions related to cephalosporin and β-lactamase inhibitor resistance in almost all cases increased BLI susceptibility, indicating the rarity of the combined phenotype.

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.

Detecting Protein-Protein Interaction Based on Protein Fragment Complementation Assay

Proteins are the most critical executive molecules by responding to the instructions stored in the genetic materials in any form of life. More frequently, proteins do their jobs by acting as a roleplayer that interacts with other protein(s), which is more evident when the function of a protein is examined in the real context of a cell. Identifying the interactions between (or amongst) proteins is very crucial for the biochemistry investigation of an individual protein and for the attempts aiming to draw a holo-picture for the interacting members at the scale of proteomics (or protein-protein interactions mapping). Here, we introduced the currently available reporting systems that can be used to probe the interaction between candidate protein pairs based on the fragment complementation of some particular proteins. Emphasis was put on the principles and details of experimental design. These systems are dihydrofolate reductase (DHFR), β-lactamase, tobacco etch virus (TEV) protease, luciferase, β- galactosidase, GAL4, horseradish peroxidase (HRP), focal adhesion kinase (FAK), green fluorescent protein (GFP), and ubiquitin.

The impact of long-distance mutations on the Ω-loop conformation in TEM type β-lactamases

β-lactamases are hydrolytic enzymes primarily responsible for occurrence and abundance of bacteria resistant to β-lactam antibiotics. TEM type β-lactamases are formed by the parent enzyme TEM-1 and more than two hundred of its mutants. Positions for the known amino acid substitutions cover ∼30% of TEM type enzyme's sequence. These substitutions are divided into the key mutations that lead to changes in catalytic properties of β-lactamases, and the secondary ones, which role is poorly understood. In this study, Residue Interaction Networks were constructed from molecular dynamic trajectories of β-lactamase TEM-1 and its variants with two key substitutions, G238S and E240K, and their combinations with secondary ones (M182T and Q39K). Particular attention was paid to a detailed analysis of the interactions that affect conformation and mobility of the Ω-loop, representing a part of the β-lactamase active site. It was shown that key mutations weakened the stability of contact inside the Ω-loop thus increasing its mobility. Combination of three amino acid substitutions, including the 182 residue, leads to the release of R65 promoting its new contacts with N175 and D176. As a result, Ω-loop is fixed on the protein globule. The second distal mutation Q39K prevents changes in spatial position of R65, which lead to the weakening of the effect of M182T substitution and the recovery of the Ω-loop mobility. Thus, the distal secondary mutations are directed for recovering the mobility of enzyme disturbed by the key mutations responsible for expansion of substrate specificity. AbbreviationsESBLextended spectrum beta-lactamasesIRinhibitor resistant beta-lactamasesMDmolecular dynamicsRINresidue interaction networksRMSDroot mean square deviationRMSFroot mean square fluctuations.Communicated by Ramaswamy H. Sarma.

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.

The kinetics of TEM1 antibiotic degrading enzymes that are displayed on Ure2 protein nanofibrils in a flow reactor

Enzymatic functionalization of cross-β structured protein nanofibrils has hitherto resulted in a severe reduction of the catalytic efficiency of high turnover biocatalysts. It has been speculated that steric restrictions and mass transport pose limits on the attached enzymes, but detailed kinetics analyzing this have not yet been reported. For a more comprehensive understanding, we studied protein nanofibrils endowed with TEM1, a β-lactamase from Escherichia coli. The packing density of TEM1 along the fibrils was controlled by co-fibrillation; in other words, the N-terminal ureidosuccinate transporter Ure2(1–80) from Saccharomyces cerevisiae was simultaneously aggregated with the chimeric proteins TEM1-Ure2(1–80). The mature fibrils were trapped in a column, and the rate of ampicillin hydrolysis was recorded using a continuous substrate flow. The turnover rate was plotted as a function of substrate molecules available per enzyme per second, which demonstrated that an elevated substrate availability counteracts mass transport limitations. To analyze this data set, we derived a kinetic model, which makes it possible to easily characterize and compare enzymes packed in columns. The functional TEM1 nanofibrils possess 80% of the catalytic turnover rate compared to free TEM1 in solution. Altogether, we have created protein nanofibrils that can effectively hydrolyze β-lactam antibiotic contaminations and provided a groundwork strategy for other highly functional enzymatic nanofibrils.

Structural and Mechanistic Basis for Extended-Spectrum Drug-Resistance Mutations in Altering the Specificity of TEM, CTX-M, and KPC β-lactamases

The most common mechanism of resistance to β-lactam antibiotics in Gram-negative bacteria is the production of β-lactamases that hydrolyze the drugs. Class A β-lactamases are serine active-site hydrolases that include the common TEM, CTX-M, and KPC enzymes. The TEM enzymes readily hydrolyze penicillins and older cephalosporins. Oxyimino-cephalosporins, such as cefotaxime and ceftazidime, however, are poor substrates for TEM-1 and were introduced, in part, to circumvent β-lactamase-mediated resistance. Nevertheless, the use of these antibiotics has lead to evolution of numerous variants of TEM with mutations that significantly increase the hydrolysis of the newer cephalosporins. The CTX-M enzymes emerged in the late 1980s and hydrolyze penicillins and older cephalosporins and derive their name from the ability to also hydrolyze cefotaxime. The CTX-M enzymes, however, do not efficiently hydrolyze ceftazidime. Variants of CTX-M enzymes, however, have evolved that exhibit increased hydrolysis of ceftazidime. Finally, the KPC enzyme emerged in the 1990s and is characterized by its broad specificity that includes penicillins, most cephalosporins, and carbapenems. The KPC enzyme, however, does not efficiently hydrolyze ceftazidime. As with the TEM and CTX-M enzymes, variants have recently evolved that extend the spectrum of KPC β-lactamase to include ceftazidime. This review discusses the structural and mechanistic basis for the expanded substrate specificity of each of these enzymes that result from natural mutations that confer oxyimino-cephalosporin resistance. For the TEM enzyme, extended-spectrum mutations act by establishing new interactions with the cephalosporin. These mutations increase the conformational heterogeneity of the active site to create sub-states that better accommodate the larger drugs. The mutations expanding the spectrum of CTX-M enzymes also affect the flexibility and conformation of the active site to accommodate ceftazidime. Although structural data are limited, extended-spectrum mutations in KPC may act by mediating new, direct interactions with substrate and/or altering conformations of the active site. In many cases, mutations that expand the substrate profile of these enzymes simultaneously decrease the thermodynamic stability. This leads to the emergence of additional global suppressor mutations that help correct the stability defects leading to increased protein expression and increased antibiotic resistance.

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.

Injection of T3SS effectors not resulting in invasion is the main targeting mechanism of Shigella toward human lymphocytes

The enteroinvasive bacterium Shigella is a facultative intracellular bacterium known, in vitro, to invade a large diversity of cells through the delivery of virulence effectors into the cell cytoplasm via a type III secretion system (T3SS). Here, we provide evidence that the injection of T3SS effectors does not necessarily result in cell invasion. Indeed, we demonstrate through optimization of a T3SS injection reporter that effector injection without subsequent cell invasion, termed the injection-only mechanism, is the main strategy used by Shigella to target human immune cells. We show that in vitro-activated human peripheral blood B, CD4⁺ T, and CD8⁺ T lymphocytes as well as switched memory B cells are mostly targeted by the injection-only mechanism. B and T lymphocytes residing in the human colonic lamina propria, encountered by Shigella upon its crossing of the mucosal barrier, are also mainly targeted by injection-only. These findings reveal that cells refractory to invasion can still be injected, thus extending the panel of host cells manipulated to the benefit of the pathogen. Future analysis of the functional consequences of the injection-only mechanism toward immune cells will contribute to the understanding of the priming of adaptive immunity, which is known to be altered during the course of natural Shigella infection.

Mistranslation can enhance fitness through purging of deleterious mutations

Phenotypic mutations are amino acid changes caused by mistranslation. How phenotypic mutations affect the adaptive evolution of new protein functions is unknown. Here we evolve the antibiotic resistance protein TEM-1 towards resistance on the antibiotic cefotaxime in an Escherichia coli strain with a high mistranslation rate. TEM-1 populations evolved in such strains endow host cells with a general growth advantage, not only on cefotaxime but also on several other antibiotics that ancestral TEM-1 had been unable to deactivate. High-throughput sequencing of TEM-1 populations shows that this advantage is associated with a lower incidence of weakly deleterious genotypic mutations. Our observations show that mistranslation is not just a source of noise that delays adaptive evolution. It could even facilitate adaptive evolution by exacerbating the effects of deleterious mutations and leading to their more efficient purging. The ubiquity of mistranslation and its effects render mistranslation an important factor in adaptive protein evolution.

Mistranslation results in amino acid changes in proteins known as phenotypic mutations and these occur at a much higher rate than DNA mutations. Here, the authors show that mistranslation can increase the response to directional selection by exacerbating the fitness effects of deleterious DNA mutations.

Identifying ionic interactions within a membrane using BLaTM, a genetic tool to measure homo- and heterotypic transmembrane helix-helix interactions

The assembly of integral membrane protein complexes is frequently supported by transmembrane domain (TMD) interactions. Here, we present the BLaTM assay that measures homotypic as well as heterotypic TMD-TMD interactions in a bacterial membrane. The system is based on complementation of β-lactamase fragments genetically fused to interacting TMDs, which confers ampicillin resistance to expressing cells. We validated BLaTM by showing that the assay faithfully reports known sequence-specific interactions of both types. In a practical application, we used BLaTM to screen a focussed combinatorial library for heterotypic interactions driven by electrostatic forces. The results reveal novel patterns of ionizable amino acids within the isolated TMD pairs. Those patterns indicate that formation of heterotypic TMD pairs is most efficiently supported by closely spaced ionizable residues of opposite charge. In addition, TMD heteromerization can apparently be driven by hydrogen bonding between basic or between acidic residues.

Modelling proteins' hidden conformations to predict antibiotic resistance

TEM β-lactamase confers bacteria with resistance to many antibiotics and rapidly evolves activity against new drugs. However, functional changes are not easily explained by differences in crystal structures. We employ Markov state models to identify hidden conformations and explore their role in determining TEM's specificity. We integrate these models with existing drug-design tools to create a new technique, called Boltzmann docking, which better predicts TEM specificity by accounting for conformational heterogeneity. Using our MSMs, we identify hidden states whose populations correlate with activity against cefotaxime. To experimentally detect our predicted hidden states, we use rapid mass spectrometric footprinting and confirm our models' prediction that increased cefotaxime activity correlates with reduced Ω-loop flexibility. Finally, we design novel variants to stabilize the hidden cefotaximase states, and find their populations predict activity against cefotaxime in vitro and in vivo. Therefore, we expect this framework to have numerous applications in drug and protein design.

Epistasis in protein evolution

The structure, function, and evolution of proteins depend on physical and genetic interactions among amino acids. Recent studies have used new strategies to explore the prevalence, biochemical mechanisms, and evolutionary implications of these interactions-called epistasis-within proteins. Here we describe an emerging picture of pervasive epistasis in which the physical and biological effects of mutations change over the course of evolution in a lineage-specific fashion. Epistasis can restrict the trajectories available to an evolving protein or open new paths to sequences and functions that would otherwise have been inaccessible. We describe two broad classes of epistatic interactions, which arise from different physical mechanisms and have different effects on evolutionary processes. Specific epistasis-in which one mutation influences the phenotypic effect of few other mutations-is caused by direct and indirect physical interactions between mutations, which nonadditively change the protein's physical properties, such as conformation, stability, or affinity for ligands. In contrast, nonspecific epistasis describes mutations that modify the effect of many others; these typically behave additively with respect to the physical properties of a protein but exhibit epistasis because of a nonlinear relationship between the physical properties and their biological effects, such as function or fitness. Both types of interaction are rampant, but specific epistasis has stronger effects on the rate and outcomes of evolution, because it imposes stricter constraints and modulates evolutionary potential more dramatically; it therefore makes evolution more contingent on low-probability historical events and leaves stronger marks on the sequences, structures, and functions of protein families.

Removal of the Side Chain at the Active-Site Serine by a Glycine Substitution Increases the Stability of a Wide Range of Serine β-Lactamases by Relieving Steric Strain

Serine β-lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. They utilize an active-site serine residue as a nucleophile, forming an acyl-enzyme intermediate during hydrolysis. In this study, thermal denaturation experiments as well as X-ray crystallography were performed to test the effect of substitution of the catalytic serine with glycine on protein stability in serine β-lactamases. Six different enzymes comprising representatives from each of the three classes of serine β-lactamases were examined, including TEM-1, CTX-M-14, and KPC-2 of class A, P99 of class C, and OXA-48 and OXA-163 of class D. For each enzyme, the wild type and a serine-to-glycine mutant were evaluated for stability. The glycine mutants all exhibited enhanced thermostability compared to that of the wild type. In contrast, alanine substitutions of the catalytic serine in TEM-1, OXA-48, and OXA-163 did not alter stability, suggesting removal of the Cβ atom is key to the stability increase associated with the glycine mutants. The X-ray crystal structures of P99 S64G, OXA-48 S70G and S70A, and OXA-163 S70G suggest that removal of the side chain of the catalytic serine releases steric strain to improve enzyme stability. Additionally, analysis of the torsion angles at the nucleophile position indicates that the glycine mutants exhibit improved distance and angular parameters of the intrahelical hydrogen bond network compared to those of the wild-type enzymes, which is also consistent with increased stability. The increased stability of the mutants indicates that the enzyme pays a price in stability for the presence of a side chain at the catalytic serine position but that the cost is necessary in that removal of the serine drastically impairs function. These findings support the stability-function hypothesis, which states that active-site residues are optimized for substrate binding and catalysis but that the requirements for catalysis are often not consistent with the requirements for optimal stability.

Residue proximity information and protein model discrimination using saturation-suppressor mutagenesis

Identification of residue-residue contacts from primary sequence can be used to guide protein structure prediction. Using Escherichia coli CcdB as the test case, we describe an experimental method termed saturation-suppressor mutagenesis to acquire residue contact information. In this methodology, for each of five inactive CcdB mutants, exhaustive screens for suppressors were performed. Proximal suppressors were accurately discriminated from distal suppressors based on their phenotypes when present as single mutants. Experimentally identified putative proximal pairs formed spatial constraints to recover >98% of native-like models of CcdB from a decoy dataset. Suppressor methodology was also applied to the integral membrane protein, diacylglycerol kinase A where the structures determined by X-ray crystallography and NMR were significantly different. Suppressor as well as sequence co-variation data clearly point to the X-ray structure being the functional one adopted in vivo. The methodology is applicable to any macromolecular system for which a convenient phenotypic assay exists.

DOI:http://dx.doi.org/10.7554/eLife.09532.001

Common techniques to determine the three-dimensional structures of proteins can help researchers to understand these molecules’ activities, but are often time-consuming and do not work for all proteins. Proteins are made of chains of amino acids. When a protein chain folds, some of these amino acids interact with other amino acids and these contacts dictate the overall shape of the protein. This means that identifying the pairs of contacting amino acids could make it possible to predict the protein’s structure.

Interactions between pairs of contacting amino acids tend to remain conserved throughout evolution, and if a mutation alters one of the amino acids in a pair then a 'compensatory' change often occurs to alter the second amino acid as well. Compensatory mutations can suggest that two amino acids are close to each other in the three-dimensional shape of a protein, but the computational methods used to identify such amino acid pairs can sometimes be inaccurate.

In 2012, researchers generated mutants of a bacterial protein called CcdB with changes to single amino acids that caused the protein to fail to fold correctly. Now, Sahoo et al. – who include two of the researchers involved in the 2012 work – have developed an experimental method to identify contacting amino acids and use the CcdB protein as a test case. The approach involved searching for additional mutations that could restore the activity of five of the original mutant proteins when the proteins were produced in yeast cells. The rationale was that any secondary mutations that restored the activity must have corrected the folding defect caused by the original mutation. Sahoo et al. then predicted how close the amino acids affected by the secondary mutations were to the amino acids altered by the original mutations. This information was used to select reliable three-dimensional models of CcdB from a large set of possible structures that had been generated previously using computer models.

Next, the technique was applied to a protein called diacylglycerol kinase A. The structure of this protein had previously been inferred using techniques such as X-ray crystallography and nuclear magnetic resonance, but there was a mismatch between the two methods. Sahoo et al. found that the amino acid contacts derived from their experimental method matched those found in the crystal structure, suggesting that the functional protein structure in living cells is similar to the crystal structure. In the future, the experimental approach developed in this work could be combined with existing methods to reliably guide protein structure prediction.

DOI:http://dx.doi.org/10.7554/eLife.09532.002

Prevalence of bla TEM-220 gene in Penicillinase-producing Neisseria gonorrhoeae strains carrying Toronto/Rio plasmid in Argentina, 2002 - 2011

Background

Penicillinase-producing Neisseria gonorroheae (PPNG) was first isolated in 1976. PPNG strains carrying bla_TEM-1 and bla_TEM-135 gene have been described in different countries. Recently, a novel bla_TEM-220 allele was detected in PPNG isolates carrying Toronto/Rio plasmid. The prevalence and characteristics of TEM-220 strains worldwide are unknown, and therefore, it needs to be studied. The purpose of this study was to detect bla_TEM-220 gene in PPNG strains possessing Toronto/Rio plasmid over a period of ten years in Argentina, and to evaluate the proportion of isolates producing non-TEM-220 containing the T539C substitution in the bla_TEM allele.

Methods

One hundred and fifty one PPNG isolates carrying Toronto/Rio plasmid were studied between 2002 and 2011. A mismatch amplification mutation assay (MAMA) PCR was used to identify the T539C substitution in the bla_TEM allele and a MAMA-PCR protocol was developed to detect the G547A substitution in the bla_TEM-220. The reference agar dilution method of the Clinical and Laboratory Standard Institute (CLSI) was used for susceptibility testing to five β-lactams antibiotics, ciprofloxacin, tetracycline and azithromycin. In all TEM-220-producing isolates, the whole bla_TEM gene was sequenced and the isolates were typed using N. gonorroheae multiantigen sequence typing (NG-MAST).

Results

MAMA PCR successfully identified the G547A substitution in the bla_TEM-220 allele. The proportion of isolates that possessed the bla_TEM-220 allele was 2.6 %, and 93.2 % MAMA TEM-220 PCR-negative isolates showed the T539C substitution in the bla_TEM gene. No differences in the susceptibility to five beta-lactam antibiotics tested were observed in PPNG isolates TEM-220-producing and PPNG isolates carrying the T539C substitution in the bla_TEM gene. All TEM-220 isolates were indistinguishable by NG-MAST.

Conclusion

This is the first study which shows the prevalence of bla_TEM-220 in N. gonorrhoeae isolates carrying Toronto/Rio plasmid in Argentina. Although the bla_TEM-220 allele does not appear to be associated with an extended spectrum beta-lactamase (ESBL) phenotype of resistance, a single nucleotide polymorphism added to the bla_TEM-220 or bla_TEM containing the T539C substitution could lead to the emergence of ESBL. Thus, it is imperative to investigate in surveillance programs, not only the plasmid type in PPNG isolates and the bla_TEM allele associated, but phenotypical characteristics and geographical distribution of isolates.

SHV-129: A Gateway to Global Suppressors in the SHV β-Lactamase Family?

Enzymes are continually evolving in response to environmental pressures. In order to increase enzyme fitness, amino acid substitutions can occur leading to a changing function or an increased stability. These evolutionary drivers determine the activity of an enzyme and its success in future generations in response to changing conditions such as environmental stressors or to improve physiological function allowing continual persistence of the enzyme. With recent warning reports on antibiotic resistance and multidrug resistant bacterial infections, understanding the evolution of β-lactamase enzymes, which are a large contributor to antibiotic resistance, is increasingly important. Here, we investigated a variant of the SHV β-lactamase identified from a clinical isolate of Escherichia coli in 2011 (SHV-129, G238S-E240K-R275L-N276D) to identify the first instance of a global suppressor substitution in the SHV β-lactamase family. We have used this enzyme to show that several evolutionary principles are conserved in different class A β-lactamases, such as active site mutations reducing stability and requiring compensating suppressor substitutions in order to ensure evolutionary persistence of a given β-lactamase. However, the pathway taken by a given β-lactamase in order to reach its evolutionary peak under a given set of conditions is likely different. We also provide further evidence for a conserved stabilizing substitution among class A β-lactamases, the back to consensus M182T substitution. In addition to expanding the spectrum of β-lactamase activity to include the hydrolysis of cefepime, the amino acid substitutions found in SHV-129 provide the enzyme with an excess of stability, which expands the evolutionary landscape of this enzyme and may result in further evolution to potentially include resistance to carbapenems or β-lactamase inhibitors.

Characterization of the global stabilizing substitution A77V and its role in the evolution of CTX-M β-lactamases

The widespread use of oxyimino-cephalosporin antibiotics drives the evolution of the CTX-M family of β-lactamases that hydrolyze these drugs and confer antibiotic resistance. Clinically isolated CTX-M enzymes carrying the P167S or D240G active site-associated adaptive mutation have a broadened substrate profile that includes the oxyimino-cephalosporin antibiotic ceftazidime. The D240G substitution is known to reduce the stability of CTX-M-14 β-lactamase, and the P167S substitution is shown here to also destabilize the enzyme. Proteins are marginally stable entities, and second-site mutations that stabilize the enzyme can offset a loss in stability caused by mutations that enhance enzyme activity. Therefore, the evolution of antibiotic resistance enzymes can be dependent on the acquisition of stabilizing mutations. The A77V substitution is present in CTX-M extended-spectrum β-lactamases (ESBLs) from a number of clinical isolates, suggesting that it may be important in the evolution of antibiotic resistance in this family of β-lactamases. In this study, the effects of the A77V substitution in the CTX-M-14 model enzyme were characterized with regard to the kinetic parameters for antibiotic hydrolysis as well as enzyme expression levels in vivo and protein stability in vitro. The A77V substitution has little effect on the kinetics of oxyimino-cephalosporin hydrolysis, but it stabilizes the CTX-M enzyme and compensates for the loss of stability resulting from the P167S and D240G mutations. The acquisition of global stabilizing mutations, such as A77V, is an important feature in β-lactamase evolution and a common mechanism in protein evolution.

Quantitative Description of a Protein Fitness Landscape Based on Molecular Features

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

A triple mutant in the Ω-loop of TEM-1 β-lactamase changes the substrate profile via a large conformational change and an altered general base for catalysis

β-Lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. TEM-1 is a prevalent plasmid-encoded β-lactamase in Gram-negative bacteria that efficiently catalyzes the hydrolysis of penicillins and early cephalosporins but not oxyimino-cephalosporins. A previous random mutagenesis study identified a W165Y/E166Y/P167G triple mutant that displays greatly altered substrate specificity with increased activity for the oxyimino-cephalosporin, ceftazidime, and decreased activity toward all other β-lactams tested. Surprisingly, this mutant lacks the conserved Glu-166 residue critical for enzyme function. Ceftazidime contains a large, bulky side chain that does not fit optimally in the wild-type TEM-1 active site. Therefore, it was hypothesized that the substitutions in the mutant expand the binding site in the enzyme. To investigate structural changes and address whether there is an enlargement in the active site, the crystal structure of the triple mutant was solved to 1.44 Å. The structure reveals a large conformational change of the active site Ω-loop structure to create additional space for the ceftazidime side chain. The position of the hydroxyl group of Tyr-166 and an observed shift in the pH profile of the triple mutant suggests that Tyr-166 participates in the hydrolytic mechanism of the enzyme. These findings indicate that the highly conserved Glu-166 residue can be substituted in the mechanism of serine β-lactamases. The results reveal that the robustness of the overall β-lactamase fold coupled with the plasticity of an active site loop facilitates the evolution of enzyme specificity and mechanism.

Identification of TEM-135 β-lactamase in Neisseria gonorrhoeae strains carrying African and Toronto plasmids in Argentina

One hundred forty-three penicillinase-producing Neisseria gonorrhoeae (PPNG) isolates obtained in Argentina from 2008 and 2012 were examined to detect blaTEM-135 genes and to investigate plasmid profiles and multiantigen sequence types. Forty-two PPNG isolates were found to carry TEM-135, and two contained a new TEM derivative characterized as TEM-220. The blaTEM-135 allele was carried by the Toronto/Rio and African plasmids. Molecular epidemiology revealed that two blaTEM-135 isolates were related to previously described isolates from Thailand and China, indicating a common evolutionary origin.

Epistatically interacting substitutions are enriched during adaptive protein evolution

Most experimental studies of epistasis in evolution have focused on adaptive changes—but adaptation accounts for only a portion of total evolutionary change. Are the patterns of epistasis during adaptation representative of evolution more broadly? We address this question by examining a pair of protein homologs, of which only one is subject to a well-defined pressure for adaptive change. Specifically, we compare the nucleoproteins from human and swine influenza. Human influenza is under continual selection to evade recognition by acquired immune memory, while swine influenza experiences less such selection due to the fact that pigs are less likely to be infected with influenza repeatedly in a lifetime. Mutations in some types of immune epitopes are therefore much more strongly adaptive to human than swine influenza—here we focus on epitopes targeted by human cytotoxic T lymphocytes. The nucleoproteins of human and swine influenza possess nearly identical numbers of such epitopes. However, mutations in these epitopes are fixed significantly more frequently in human than in swine influenza, presumably because these epitope mutations are adaptive only to human influenza. Experimentally, we find that epistatically constrained mutations are fixed only in the adaptively evolving human influenza lineage, where they occur at sites that are enriched in epitopes. Overall, our results demonstrate that epistatically interacting substitutions are enriched during adaptation, suggesting that the prevalence of epistasis is dependent on the underlying evolutionary forces at play.

Mutations can fix during evolution for two reasons: they can be beneficial and fix for adaptive reasons, or they can be neutral or deleterious and fix solely by chance. Most studies focus on adaptation, where the evolving population is increasing in fitness due to a new selection pressure. Such studies have found an important evolutionary role for epistasis, the phenomenon where the effect of one mutation depends on another mutation. But adaptation only accounts for a fraction of overall evolutionary change. Here we investigate whether epistasis is as common during non-adaptive as adaptive evolution. We do this by comparing the same protein from human and swine influenza. Human influenza is constantly adapting to escape from the immunity that people acquire from previous influenza infections. But swine influenza is under less pressure to escape from acquired immunity since pigs have shorter lifetimes and are less likely to be infected with influenza multiple times. We find that epistasis is less common during the evolution of the swine influenza protein than its human influenza counterpart. Overall, our results suggest that mutations that interact via epistasis are more likely to fix during adaptive evolution.

Experimental Simulation of the Effects of an Initial Antibiotic Treatment on a Subsequent Treatment after Initial Therapy Failure

Therapy failure of empirical antibiotic treatments prescribed by primary care physicians occurs commonly. The effect of such a treatment on the susceptibility to second line antimicrobial drugs is unknown. Resistance to amoxicillin was rapidly induced or selected in E. coli at concentrations expected in the patient's body. Strains with reduced susceptibility outcompeted the wild-type whenever antibiotics were present, even in low concentrations that did not affect the growth rates of both strains. Exposure of E. coli to amoxicillin caused moderate resistance to cefotaxime. The combined evidence suggests that initial treatment by amoxicillin has a negative effect on subsequent therapy with beta-lactam antibiotics.

Capturing the mutational landscape of the beta-lactamase TEM-1

Adaptation proceeds through the selection of mutations. The distribution of mutant fitness effect and the forces shaping this distribution are therefore keys to predict the evolutionary fate of organisms and their constituents such as enzymes. Here, by producing and sequencing a comprehensive collection of 10,000 mutants, we explore the mutational landscape of one enzyme involved in the spread of antibiotic resistance, the beta-lactamase TEM-1. We measured mutation impact on the enzyme activity through the estimation of amoxicillin minimum inhibitory concentration on a subset of 990 mutants carrying a unique missense mutation, representing 64% of possible amino acid changes in that protein reachable by point mutation. We established that mutation type, solvent accessibility of residues, and the predicted effect of mutations on protein stability primarily determined alone or in combination changes in minimum inhibitory concentration of mutants. Moreover, we were able to capture the drastic modification of the mutational landscape induced by a single stabilizing point mutation (M182T) by a simple model of protein stability. This work thereby provides an integrated framework to study mutation effects and a tool to understand/define better the epistatic interactions.

Fisher's geometric model of adaptation meets the functional synthesis: data on pairwise epistasis for fitness yields insights into the shape and size of phenotype space

The functional synthesis uses experimental methods from molecular biology, biochemistry and structural biology to decompose evolutionarily important mutations into their more proximal mechanistic determinants. However these methods are technically challenging and expensive. Noting strong formal parallels between R.A. Fisher's geometric model of adaptation and a recent model for the phenotypic basis of protein evolution, we sought to use the former to make inferences into the latter using data on pairwise fitness epistasis between mutations. We present an analytic framework for classifying pairs of mutations with respect to similarity of underlying mechanism on this basis, and also show that these data can yield an estimate of the number of mutationally labile phenotypes underlying fitness effects. We use computer simulations to explore the robustness of our approach to violations of analytic assumptions and analyze several recently published datasets. This work provides a theoretical complement to the functional synthesis as well as a novel test of Fisher's geometric model.

Evolution of broad spectrum β-lactam resistance in an engineered metallo-β-lactamase

The extensive use and misuse of antibiotics during the last seven decades has led to the evolution and global spread of a variety of resistance mechanisms in bacteria. Of high medical importance are β-lactamases, a group of enzymes inactivating β-lactam antibiotics. Metallo-β-lactamases (MBLs) are particularly problematic because of their ability to act on virtually all classes of β-lactam antibiotics. An engineered MBL (evMBL9) characterized by low level activity with several β-lactam antibiotics was constructed and employed as a parental MBL in an experiment to examine how an enzyme can evolve toward increased activity with a variety of β-lactam antibiotics. We designed and synthesized a mutant library in which the substrate activity profile was varied by randomizing six active site amino acid residues. The library was expressed in Salmonella typhimurium, clones with increased resistance against seven different β-lactam antibiotics (penicillin G, ampicillin, cephalothin, cefaclor, cefuroxime, cefoperazone, and cefotaxime) were isolated, and the MBL variants were characterized. For the majority of the mutants, bacterial resistance was significantly increased despite marked reductions in both mRNA and protein levels relative to those of parental evMBL9, indicating that the catalytic activities of these mutant MBLs were highly increased. Multivariate analysis showed that the majority of the mutant enzymes were generalists, conferring increased resistance against most of the examined β-lactams.

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.

Mutagenesis of zinc ligand residue Cys221 reveals plasticity in the IMP-1 metallo-β-lactamase active site

Metallo-β-lactamases catalyze the hydrolysis of a broad range of β-lactam antibiotics and are a concern for the spread of drug resistance. To analyze the determinants of enzyme structure and function, the sequence requirements for the subclass B1 IMP-1 β-lactamase zinc binding residue Cys221 were tested by saturation mutagenesis and evaluated for protein expression, as well as hydrolysis of β-lactam substrates. The results indicated that most substitutions at position 221 destabilized the enzyme. Only the enzymes containing C221D and C221G substitutions were expressed well in Escherichia coli and exhibited catalytic activity toward β-lactam antibiotics. Despite the lack of a metal-chelating group at position 221, the C221G enzyme exhibited high levels of catalytic activity in the presence of exogenous zinc. Molecular modeling suggests the glycine substitution is unique among substitutions in that the complete removal of the cysteine side chain allows space for a water molecule to replace the thiol and coordinate zinc at the Zn2 zinc binding site to restore function. Multiple methods were used to estimate the C221G Zn2 binding constant to be 17 to 43 μM. Studies of enzyme function in vivo in E. coli grown on minimal medium showed that both IMP-1 and the C221G mutant exhibited compromised activity when zinc availability was low. Finally, substitutions at residue 121, which is the IMP-1 equivalent of the subclass B3 zinc-chelating position, failed to rescue C221G function, suggesting the coordination schemes of subclasses B1 and B3 are not interchangeable.

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.

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.

Molecular characterization of resistance to extended-spectrum cephalosporins in clinical Escherichia coli isolates from companion animals in the United States

Resistance to extended-spectrum cephalosporins (ESC) among members of the family Enterobacteriaceae occurs worldwide; however, little is known about ESC resistance in Escherichia coli strains from companion animals. Clinical isolates of E. coli were collected from veterinary diagnostic laboratories throughout the United States from 2008 to 2009. E. coli isolates (n = 54) with reduced susceptibility to ceftazidime or cefotaxime (MIC ≥ 16 μg/ml) and extended-spectrum-β-lactamase (ESBL) phenotypes were analyzed. PCR and sequencing were used to detect mutations in ESBL-encoding genes and the regulatory region of the chromosomal gene ampC. Conjugation experiments and plasmid identification were conducted to examine the transferability of resistance to ESCs. All isolates carried the bla(CTX-M-1)-group β-lactamase genes in addition to one or more of the following β-lactamase genes: bla(TEM), bla(SHV-3), bla(CMY-2), bla(CTX-M-14-like), and bla(OXA-1.) Different bla(TEM) sequence variants were detected in some isolates (n = 40). Three isolates harbored a bla(TEM-181) gene with a novel mutation resulting in an Ala184Val substitution. Approximately 78% of the isolates had mutations in promoter/attenuator regions of the chromosomal gene ampC, one of which was a novel insertion of adenine between bases -28 and -29. Plasmids ranging in size from 11 to 233 kbp were detected in the isolates, with a common plasmid size of 93 kbp identified in 60% of isolates. Plasmid-mediated transfer of β-lactamase genes increased the MICs (≥ 16-fold) of ESCs for transconjugants. Replicon typing among isolates revealed the predominance of IncI and IncFIA plasmids, followed by IncFIB plasmids. This study shows the emergence of conjugative plasmid-borne ESBLs among E. coli strains from companion animals in the United States, which may compromise the effective therapeutic use of ESCs in veterinary medicine.

Analysis of the functional contributions of Asn233 in metallo-β-lactamase IMP-1

Metallo-β-lactamases, such as IMP-1, are a major global health threat, as they provide for bacterial resistance to a wide range of β-lactam antibiotics, including carbapenems. Understanding the molecular details of the enzymatic process and the sequence requirements for function are essential aids in overcoming β-lactamase-mediated resistance. An asparagine residue is conserved at position 233 in approximately 67% of all metallo-β-lactamases. Despite its conservation, the molecular basis of Asn233 function is poorly understood and remains controversial. It has previously been shown that mutations at this site exhibit context-dependent sequence requirements in that the importance of a given amino acid depends on the antibiotic being tested. To provide a more thorough examination as to the function and sequence requirements at this position, a collection of IMP-1 mutants encoding each of the 19 possible amino acid substitutions was generated. The resistance levels toward four β-lactam antibiotics were measured for Escherichia coli containing each of these mutants. The sequence requirements at position 233 for wild-type levels of resistance toward two cephalosporins were the most relaxed, while there were more stringent sequence requirements for resistance to ampicillin or imipenem. Enzyme kinetic analysis and determinations of steady-state protein levels indicated that the effects of the substitutions on resistance are due to changes in the kinetic parameters of the enzyme. Taken together, the results indicate that substitutions at position 233 significantly alter the kinetic parameters of the enzyme, but most substituted enzymes are able to provide for a high level of resistance to a broad range of β-lactams.

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.

Mechanistic studies of the inactivation of TEM-1 and P99 by NXL104, a novel non-beta-lactam beta-lactamase inhibitor

NXL104 is a potent inhibitor of class A and C serine β-lactamases, including KPC carbapenemases. Native and NXL104-inhibited TEM-1 and P99 β-lactamases analyzed by liquid chromatography-electrospray ionization-time of flight mass spectrometry revealed that the inactivated enzymes formed a covalent adduct with NXL104. The principal inhibitory characteristics of NXL104 against TEM-1 and P99 β-lactamases were determined, including partition ratios, dissociation constants (K), rate constants for deactivation (k(2)), and reactivation rates. NXL104 is a potent inhibitor of TEM-1 and P99, characterized by high carbamylation efficiencies (k(2)/K of 3.7 × 10(5) M(-1) s(-1) for TEM-1 and 1 × 10(4) M(-1) s(-1) for P99) and slow decarbamylation. Complete loss of β-lactamase activity was obtained at a 1/1 enzyme/NXL104 ratio, with a k(3) value (rate constant for formation of product and free enzyme) close to zero for TEM-1 and P99. Fifty percent inhibitory concentrations (IC(50)s) were evaluated on selected β-lactamases, and NXL104 was shown to be a very potent inhibitor of class A and C β-lactamases. IC(50)s obtained with NXL104 (from 3 nM to 170 nM) were globally comparable on the β-lactamases CTX-M-15 and SHV-4 with those obtained with the comparators (clavulanate, tazobactam, and sulbactam) but were far lower on TEM-1, KPC-2, P99, and AmpC than those of the comparators. In-depth studies on TEM-1 and P99 demonstrated that NXL104 had a comparable or better affinity and inactivation rate than clavulanate and tazobactam and in all cases an improved stability of the covalent enzyme/inhibitor complex.

Identification and characterization of beta-lactamase inhibitor protein-II (BLIP-II) interactions with beta-lactamases using phage display

Protein-protein interactions are critical to cellular processes yet the ability to predict and rationally design interactions is limited because of incomplete knowledge of the principles governing these interactions. The beta-lactamase inhibitory protein (BLIP)/beta-lactamase interaction has become a model system to investigate protein-protein interactions and has been the focus of several structural, thermodynamic and binding specificity studies. BLIP-II also inhibits beta-lactamase but has no sequence homology with BLIP. The structure of BLIP-II in complex with TEM-1 beta-lactamase revealed that BLIP-II has a completely different structure than BLIP but it interacts with the same protruding loop-helix region of TEM-1 as does BLIP. The significance of the individual interacting residues in molecular recognition by BLIP-II is currently unknown. Therefore, a phage display vector was developed with the purpose of expressing BLIP-II onto the surface of the M13 filamentous bacteriophage. The BLIP-II displayed phage bound to TEM-1 with picomolar affinity indicating that BLIP-II is properly folded while on the surface of the phage. The phage system, as well as enzyme inhibition assays with purified proteins, revealed that BLIP-II is a more potent inhibitor than BLIP for several class A beta-lactamases with K(i) values in the low picomolar range.

Structural bases for stability-function tradeoffs in antibiotic resistance

Preorganization of enzyme active sites for substrate recognition typically comes at a cost to the stability of the folded form of the protein; consequently, enzymes can be dramatically stabilized by substitutions that attenuate the size and preorganization "strain" of the active site. How this stability-activity tradeoff constrains enzyme evolution has remained less certain, and it is unclear whether one should expect major stability insults as enzymes mutate towards new activities or how these new activities manifest structurally. These questions are both germane and easy to study in beta-lactamases, which are evolving on the timescale of years to confer resistance to an ever-broader spectrum of beta-lactam antibiotics. To explore whether stability is a substantial constraint on this antibiotic resistance evolution, we investigated extended-spectrum mutants of class C beta-lactamases, which had evolved new activity versus third-generation cephalosporins. Five mutant enzymes had between 100-fold and 200-fold increased activity against the antibiotic cefotaxime in enzyme assays, and the mutant enzymes all lost thermodynamic stability (from 1.7 kcal mol(-)(1) to 4.1 kcal mol(-)(1)), consistent with the stability-function hypothesis. Intriguingly, several of the substitutions were 10-20 A from the catalytic serine; the question of how they conferred extended-spectrum activity arose. Eight structures, including complexes with inhibitors and extended-spectrum antibiotics, were determined by X-ray crystallography. Distinct mechanisms of action, including changes in the flexibility and ground-state structures of the enzyme, are revealed for each mutant. These results explain the structural bases for the antibiotic resistance conferred by these substitutions and their corresponding decrease in protein stability, which will constrain the evolution of new antibiotic resistance.

Potential role of phenotypic mutations in the evolution of protein expression and stability

Phenotypic mutations (errors occurring during protein synthesis) are orders of magnitude more frequent than genetic mutations. Consequently, the sequences of individual protein molecules transcribed and translated from the same gene can differ. To test the effects of such mutations, we established a bacterial system in which an antibiotic resistance gene (TEM-1 beta-lactamase) was transcribed by either a high-fidelity RNA polymerase or its error-prone mutant. This setup enabled the analysis of individual mRNA transcripts that were synthesized under normal or error-prone conditions. We found that an increase of approximately 20-fold in the frequency of transcription errors promoted the evolution of higher TEM-1 expression levels and of more stable enzyme variants. The stabilized variants exhibited a distinct advantage under error-prone transcription, although under normal transcription they conferred resistance similar to wild-type TEM-1. They did so, primarily, by increasing TEM-1's tolerance to destabilizing deleterious mutations that arise from transcriptional errors. The stabilized TEM-1 variants also showed increased tolerance to genetic mutations. Thus, although phenotypic mutations are not individually subjected to inheritance and natural selection, as are genetic mutations, they collectively exert a direct and immediate effect on protein fitness. They may therefore play a role in shaping protein traits such as expression levels, stability, and tolerance to genetic mutations.

Adaptive protein evolution grants organismal fitness by improving catalysis and flexibility

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

Genetic and structural characterization of an L201P global suppressor substitution in TEM-1 beta-lactamase

TEM-1 beta-lactamase confers bacterial resistance to penicillin antibiotics and has acquired mutations that permit the enzyme to hydrolyze extended-spectrum cephalosporins or to avoid inactivation by beta-lactamase inhibitors. However, many of these substitutions have been shown to reduce activity against penicillin antibiotics and/or result in loss of stability for the enzyme. In order to gain more information concerning the tradeoffs associated with active site substitutions, a genetic selection was used to find second site mutations that partially restore ampicillin resistance levels conferred by an R244A active site TEM-1 beta-lactamase mutant. An L201P substitution distant from the active site that enhanced ampicillin resistance levels and increased protein expression levels of the R244A TEM-1 mutant was identified. The L201P substitution also increases the ampicillin resistance levels and restores expression levels of a poorly expressed TEM-1 mutant with a core-disrupting substitution. In vitro thermal denaturation of purified protein indicated that the L201P mutation increases the T(m) value of the TEM-1 enzyme. The X-ray structure of the L201P TEM-1 mutant was determined to gain insight into the increase in enzyme stability. The proline substitution occurs at the N-terminus of an alpha-helix and may stabilize the enzyme by reducing the helix dipole, as well as by lowering the conformational entropy cost of folding due to the reduced number of conformations available in the unfolded state. Collectively, the data suggest that L201P promotes tolerance of some deleterious TEM-1 mutations by enhancing the protein stability of these mutants.

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 beta-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 beta-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of beta-lactamase by almost 20 degrees C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol(-1) relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring beta-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.

Characterization of extended-spectrum beta-lactamases in Salmonella spp. at a tertiary hospital in Durban, South Africa

Extended-spectrum beta-lactamases (ESBLs) were characterized in 41 Salmonella spp. isolates from patients admitted to a pediatric ward of a tertiary hospital in Durban, South Africa. The most common (17/ 41) serotype was Salmonella enterica serotype Typhimurium, followed by S. enterica serotype Isangi (16/41), S. enterica serotype Saint-paul (2/41), S. enterica serotype Kissi (2/41), S. enterica serotype Kivu (2/41), and S. enterica serotype Reading (1/41). All isolates were resistant to ampicillin, amoxicillin-clavulanate, piperacillin, ceftazidime, and aztreonam but susceptible to meropenem. SHV-12 found in 39% of the isolates was the most common ESBL. TEM-63 was produced in 29% and TEM-116 in 10% of the isolates, and TEM-131 was found in 1 isolate. Other ESBLs that were identified included SHV-2 (n = 2), CTX-M-3 (n = 1), CTX-M-15 (n = 2), and CTX-M-37 (n = 5). In addition, CMY-2 (n = 3) and the OXA-1(n = 1) beta-lactamase were also detected. The diversity of ESBLs suggests that its incidence in Salmonellae needs to be monitored.

Intense neutral drifts yield robust and evolvable consensus proteins

What changes occur when a natural protein that had been under low mutation rates is subjected to a neutral drift at high mutational loads, thus generating genetically diverse (polymorphic) gene ensembles that all maintain the protein's original function and structure? To address this question we subjected large populations of TEM-1 beta-lactamase to a prolonged neutral drift, applying high mutation rates and purifying selection to maintain TEM-1's existing penicillinase activity. Purging of deleterious mutations and enrichment of beneficial ones maintained the sequence of these ensembles closer to TEM-1's family consensus and inferred ancestor. In particular, back-to-consensus/ancestor mutations that increase TEM-1's kinetic and thermodynamic stability were enriched. These acted as global suppressors and enabled the tolerance of a broad range of deleterious mutations, thus further increasing the genetic diversity of the drifting populations. The probability of a new function emerging (cefotaxime degradation) was also substantially increased in these ensembles owing to the presence of many gene variants carrying the global suppressors. Our findings indicate the unique features of large, polymorphic neutral ensembles generated under high mutational loads and prompt the speculation that the progenitors of today's proteins may have evolved under high mutational loads. The results also suggest that predictable back-to-consensus/ancestor changes can be used in the laboratory to generate highly diverse and evolvable gene libraries.

E240V substitution increases catalytic efficiency toward ceftazidime in a new natural TEM-type extended-spectrum beta-lactamase, TEM-149, from Enterobacter aerogenes and Serratia marcescens clinical isolates

The aim of this study was to characterize a novel extended-spectrum beta-lactamase that belongs to the TEM family, the TEM-149 enzyme, and that was isolated from the urine of two hospitalized patients from different hospitals in southern Italy. The peculiarity of this enzyme was the finding of a valine residue at position 240. The array of amino acid substitutions found in TEM-149 was as follows: E104K, R164S, M182T, and E240V. A reversion of a threonine residue at position 182 was also performed to create a new mutant, TEM-149 T182M, in order to assess the contribution of this substitution on the kinetic profile and the stability of TEM-149. The bla TEM-149 and bla TEM-149/T182M genes were cloned into pBC-SK, and the corresponding enzymes were purified from recombinant Escherichia coli HB101 by the same procedure. Both enzymes hydrolyzed all beta-lactams tested, with a preference for ceftazidime, which was found to be the best substrate. By comparison of the kinetic parameters of the TEM-149 and the TEM-149 T182M enzymes, a reduction of the catalytic efficiency for the TEM-149 T182M mutant was observed against all substrates tested except benzylpenicillin, cefotaxime, and aztreonam. Tazobactam, clavulanic acid, and sulbactam were good inhibitors of the TEM-149 beta-lactamase.

The Bacillus licheniformis BlaP beta-lactamase as a model protein scaffold to study the insertion of protein fragments

Using genetic engineering technologies, the chitin-binding domain (ChBD) of the human macrophage chitotriosidase has been inserted into the host protein BlaP, a class A beta-lactamase produced by Bacillus licheniformis. The product of this construction behaved as a soluble chimeric protein that conserves both the capacity to bind chitin and to hydrolyze beta-lactam moiety. Here we describe the biochemical and biophysical properties of this protein (BlaPChBD). This work contributes to a better understanding of the reciprocal structural and functional effects of the insertion on the host protein scaffold and the heterologous structured protein fragments. The use of BlaP as a protein carrier represents an efficient approach to the functional study of heterologous protein fragments.

Polyhead formation in phage P22 pinpoints a region in coat protein required for conformational switching

Eighteen single amino acid substitutions in phage P22 coat protein cause temperature-sensitive folding defects (tsf). Three intragenic global suppressor (su) substitutions (D163G, T166I and F170L), localized to a flexible loop, rescue the folding of several tsf coat proteins. Here we investigate the su substitutions in the absence of the original tsf substitutions. None of the su variant coat proteins displayed protein folding defects. Individual su substitutions had little effect on phage production in vivo; yet double and triple combinations resulted in a cold-sensitive (cs) phenotype, consistent with a defect in assembly. During virus assembly and maturation, conformational switching of capsid subunits is required when chemically identical capsid subunits form an icosahedron. Analysis of double- and triple-su phage-infected cell lysates by negative-stain electron microscopy reveals an increase in aberrant structures at the cs temperature. In vitro assembly of F170L coat protein causes production of polyheads, never seen before in phage P22. Purified procapsids composed of all of the su coat proteins showed defects in expansion, which mimics maturation in vitro. Our results suggest that a previously identified surface-exposed loop in coat protein is critical in conformational switching of subunits during both procapsid assembly and maturation.