Molecular Analysis of β-Lactamase Genes to Understand their Differential Expression in Strains of Yersinia enterocolitica Biotype 1A

Beta-lactams are used as major therapeutic agents against a number of infectious agents. Due to widespread use of β-lactams, β-lactamases have evolved at a rapid pace leading to treatment failures. Yersinia enterocolitica causes many gastrointestinal problems. It is an extremely heterogeneous species comprising more than fifty serotypes and six biotypes which differ in their ecological niches, geographical distribution and pathogenic potential. Though biotype 1A strains have been associated with outbreaks of Yersiniosis, there has been a controversy regarding their pathogenicity. The strains of Y. enterocolitica isolated from India belonged to biotype 1A and possessed genes for two β-lactamases namely, blaA and blaB. An earlier study by us reported differential expression of blaA by strains of Y. enterocolitica biotype 1A. The present study has been carried out to understand the molecular bases which regulate the expression of blaA in Y. enterocolitica biotype 1A. We concluded that six types of blaA variants were present in strains of biotype 1A. Neither amino acid substitutions in blaA nor mutations in promoter regions of blaA contributed to differential expression of blaA in Y. enterocolitica biotype 1A. Rather, the secondary structures attained by mRNA of blaA might underlie the differential expression of blaA in Y. enterocolitica.

with high-level expression of chromosomal b-lactamases 17 . Studies on bla genes of Y. enterocolitica have reported differences in b-lactamase expression and b-lactam susceptibility among different biotypes, focusing primarily on 'pathogenic biotypes' [18][19][20][21] . The previous study on molecular characterization of bla genes and differences in b-lactamase expression in biotype 1A strains was reported from our laboratory 14,15 . Previous studies from our laboratory also reported that the expression of blaB in biotype 1A strains was quite similar among the strains, while the expression of blaA was quite variable among different strains 14,15 . To the best of our knowledge, till date no study has been carried out to understand the differential regulation of expression of blaA in Y. enterocolitica biotype 1A strains. The present study is an extension of our previous observations and was undertaken to understand the molecular basis of differential expression of b-lactamase genes of biotype 1A strains of Y. enterocolitica. This study also provides insights into the blaA variants in biotype 1A strains of Y. enterocolitica. Thus, the promoter regions and the structural regions of blaA genes in biotype 1A strains of Yersinia enterocolitica isolated from India and from other parts of the world were PCR-amplified, sequenced and compared. Three dimensional structures of blaA variants were modeled in-silico to co-relate the role of specific substitutions in differential expression of b-lactamases.

Results
Minimal Inhibitory Concentrations of b-lactam antibiotics. Y.enterocolitica strains showed variable susceptibilities to b-lactam antibiotics of different generations (Table 1). All biotype 1A strains of Indian and non-Indian origin except Y.enterocolitica strain E3 were resistant to amoxicillin and amoxicilln-clavulanate, antibiotic-inhibitor combination. The MIC of all biotype 1A strains to amoxicillin exceeded the upper range of E-test i.e . 256 mg/L. For six Indian strains, the MIC of amoxicilln-clavulanate combination also exceeded the upper range of E-test. The susceptibility of Y.enterocolitica to the third generation cephalosporins, cefpodoxime and cefotaxime was also tested by E-test. Antibiotic susceptibility of Y.enterocolitica strains to the third generation cephalosporin cefpodoxime varied and ranged from 0.125-8 mg/L. While most of the biotype 1A strains were either sensitive or intermediate to cefpodoxime, four strains were resistant to this antibiotic. Further, all Y. strains were sensitive to third generation cephalosporin, cefotaxime. For biotype 1A strains the MIC was lowest for cefotaxime and ranged between 0.032-0.64 mg/L. Strain E3 showed lowest MIC to all antibiotics including amoxicillin and amoxicillin-clavulanate.
Specific activities of b-lactamases. The mean specific activities of the blaA variants produced by Y. enterocolitica biotype 1A strains expressed as amount of nitrocefin hydrolyzed (mmol/min/mg of protein) are summarized in Table 2. The specific activity of the variant blaA2 was quite high than the mean specific activities of other variants, while that of blaA4 was the least. The mean specific activities of the variants blaA1 and blaA3 were similar to each other.
Promoter sequences and CCDS of blaA gene. The amplification of blaA gene using primers blaA11 and blaA12 resulted in the expected amplicon of 1.49 kb in 68 out of 81 strains of Y. enterocolitica. The band was excised from the agarose gel and sequenced after purification using QIAGEN purification kit. The BLAST analysis of the sequences confirmed these to be part of the blaA gene. Multiple sequence alignment of the promoter regions revealed that promoter sequences of all biotype 1A strains were identical (Fig. 1). The amplification of blaA gene using published primers 15 resulted in the amplification of desired product of 896 bp in all the strains. The band was excised from gel, purified using QIAGEN purification kit and sequenced. CCDSs of blaA were translated into their corresponding amino acid sequences and analyzed. Comparison of amino acid sequences of blaA using multiple sequence alignment revealed that though the sequence was conserved, some substitutions were identified. Since the initial 30 amino acids constitute the signal sequence and are cleaved before the mature enzyme is released in the periplasmic space, these were excluded from comparative studies. P32S and G60C was observed in 12 of  Table 3 and summarized in Fig. 2.
Phylogenetic analysis of blaA of Y. enterocolitica strains. The phylogenetic analysis of blaA of 35 strains of Y. enterocolitica biotype 1A was carried using the Maximum Likelihood method of MEGA5.5. The tree with the highest log likelihood (2939.3021) is shown (Fig. 3). The tree is drawn to scale, with branch lengths   represented by the number of substitutions per site. The MEGA analysis clustered isolates into three major groups.
Modeling analysis. Based on the identity percentage and E-values (expectation values), b-lactamase of Burkholderia multivorans (PDB code 3W4Q_A; Uniprot ID A9ANW2) was selected as the template protein due its high sequence similarity, identity and low E-value. Of the 20 models built using MODELLER, the 3D model with the lowest modeler objective function was selected as the final model.
Protein structure validation by PROCHECK and VERIFY3D. The results of PROCHECK as represented in Ramachandran contour plot showed that ca. 90% of the residues lay in the most favoured regions and only 1% of residues in disallowed regions. Verify3D profiles which were used to identify unreliable regions that had been modeled improperly showed that 90% region of the protein model of each blaA variant scored . 0.2 which was highly significant.
Molecular docking. Evaluation of binding affinity of b-lactam antibiotics with class A b-lactamases was performed using AutoDock Vina. The binding poses for each ligand molecule (amoxicillin, clavulanic acid) into the structures were determined and different poses were generated based on the total score (Dock score). The negative low free energy of binding of docked complexes and low inhibition constants indicated high affinity of b-lactamase for these antibiotics (Table 4 & 5, Fig. 4).
mRNA secondary structure analysis. The secondary structures of mRNA of the six blaA variants were analyzed. These formed branched secondary structures. The structure of blaA2 showed very high level of branching, with a number of bulge loops, interior loops and a smaller region of double stranded RNA helix. On the other hand the mRNA secondary structure of blaA4 was mostly composed of a double-stranded RNA helix of stacked base pairs with little branching (Fig. 5).

Discussion
The MIC of different b-lactam antibiotics for Y. enterocolitica biotype 1A strains using E-test revealed that all strains were resistant to amoxicillin and amoxicillin-clavulanate. Y. enterocolitica strains have been reported to be 'naturally resistant' to ampicillin 30,31 .
Although the MIC of amoxicillin-clavulanate was perceptibly less than that of amoxicillin, it did not change the resistance phenotype, implying the existence of inhibitor resistant blaAs in Y. enterocolitica biotype 1A. Earlier studies have shown that Y. enterocolitica strains of serotypes O: 3, O: 5 and O: 9 31 were all resistant to ampicillin and a few were resistant to amoxicillin-clavulanate. Szych et al. 32 reported that though all strains of bioserotypes 4/O: 3 and 1B/O: 8 were resistant to ampicillin, yet 90% of these were sensitive to amoxicillin-clavulanate. They also reported that strains of bioserotype 2/O: 9 were resistant to both ampicillin and amoxicillin-clavulanate. Recent studies 8,33 also showed that strains of Y.enterocolitica bioserotype 4/ O: 3 were sensitive to amoxicillin-clavulanate where as that of bioserotypes 2/O: 5, 27 and 2/O: 9 were resistant. From the available information it is apparent that both the biotype and serotype play role in affecting the susceptibility of Y. enterocolitica to amoxicillinclavulanate. Since all strains in the present study belonged to biotype 1A, it is not possible to comment on the influence of biotype on antibiotic susceptibility; but based on our findings it may be concluded that serotype had no effect on the susceptibility of Y. enterocoloitica biotype 1A strains to amoxicillin-clavulanate. The biotype 1A strains were either sensitive or showed intermediate susceptibility to the third generation cephalosporin, cefpodoxime and only a few strains were resistant. This is the same as observed in an earlier study 8 . Y. enterocolitica strain E3 was sensitive to all the b-lactam antibiotics tested and showed lowest MIC to all the b-lactam antibiotics tested including amoxicillin and amoxicillin-clavulanate combination.
The mean specific activities of variants blaA1 and blaA3 were almost similar and the mutations/substitutions which lead to these variants had no effect on the enzyme activity. Similarly mutations that generated variants blaA5 and blaA6 did not affect the enzyme activity significantly. However, the mutations leading to variants blaA2 and blaA4 significantly increased or decreased the specific activities of the enzyme respectively. Bacteria have been reported to hyper produce b-lactamases by different mechanisms viz. mutation in the promoters, integration of insertion sequences containing efficient promoters or point mutations in the structural gene 17 . In the family Enterobacteriaceae, mutations and insertions, mostly in the promoter regions of the b-lactamase genes, have been associated with hyper production of b-lactamases [34][35][36] . In the present study, we attempted to understand the molecular basis of differential expression of chromosomal blaA of biotype1A strains of Y. enterocolitica. For this, the nucleotide sequences of promoter regions of blaA of biotype 1A strains were PCR-amplified and analyzed. Genetic events that affect the transcription efficacy of bacterial promoters usually occur in the 235 or 210 regions or sometimes in the space between these regions. It is well known that promoter strength is linked to the sequence of the 235 and 210 boxes, creating the 70 bp contact region for RNA polymerase. The 210 sequence TATAAT and the 235 sequence GTGGTG of blaA of Y. enteocolitica were similar to the canonical E.coli sequence indicating these were strong promoters. Alignment and comparison of blaA of biotype 1A strains by Clustal revealed that the ribosomal binding site (RBS), the 210 region and 235 regions were identical in all strains indicating that the promoter region had no role in differential expression of b-lactamases. Thus it was argued that the differential expression of class A beta-lactamase enzyme in Y. enterocolitica biotype 1A might be due to the changes in gene sequence of blaA.
To investigate the role of changes in gene sequence of blaA in differential expression of enzyme, blaA genes of biotype 1A strains were PCR-amplified and their CCDSs translated into corresponding amino acid sequences. The four significant motifs conserved across the Ambler class A b-lactamase viz. SXXK, SDN, EXXLN and KTG, (where X represents any amino acid) were present in blaA of all Y. enterocolitica strains, their respective sequence and position being 75 STFK 78 , 136 SDN 138 , 171 EPDLN 175 and 239 KTG 241 . One of the four conserved domains 171 EPDLN 175 was located in the omega loop, a structural domain constituting a part of active-site pocket. Clustal alignment and comparison of amino acid sequences of structural blactamase blaA of biotype 1A strains revealed that amino acid substitutions were present at some places. Since the initial 30 amino acids constitute the signal sequence and are cleaved before the      uncharged amino acid, while aspartate is a negatively charged amino acid; substitution of N101D in Y. enterocolitica strain 19 may influence enzyme activity due to change in the properties of the side chain. Substitution of I111V and T225S in Y. enterocolitica strain E4 and L144V in Y. enterocolitica strains E2 and E4 would also affect only the bulkiness of the side chain of the protein without any profound effect on enzyme activity because of the similarity in the chemical nature of the substituent and the substitute. For example isoleucine and valine are nonaliphatic amino acids, threonine and serine are polar amino acids and leucine and valine both nonpolar amino acids. Phylogenetic analysis based on amino acid sequences of blaA clustered the biotype 1A strains into 3 major groups. In our earlier study also the Y. enterocolitica biotype 1A strains clustered into three groups based on the profiles generated by the restriction analysis of blaA 15 .
The analysis of amino acid sequences of blaA of Y. enterocolitica biotype 1A strains showed that although none of the substitution/ mutation occurred in the four major motifs conserved across class A b-lactamases, or in the omega-loop region of blaA, yet there were differences in the specific activities of blaA. Amino acid substitutions which do not occur in the active site of the enzyme may create diverse local changes in the 3D structure of the enzyme and increase its conformational flexibility which may in turn affect the binding affinity with the substrate. To evaluate the role of such amino acid substitutions on binding affinity we classified the blaA of biotype 1A strains of Y. enterocolitica into six variants based on the mutation/ substitution in the amino acid sequence. Three dimensional modeling of the blaA variants was carried out to co-relate structural changes with the higher enzyme specific activities and MIC of different strains. In-silico docking of the blaA variants was carried out with amoxicillin and clavulanic acid to analyze if amino acid substitutions affected the residues involved in hydrogen bonding and hydrophobic interactions in the enzyme-ligand complex. It is well known that hydrogen bonds and hydrophobic interactions not only play a crucial role in molecular recognition and the overall stability of the protein structure; these also stabilize the protein-ligand complexes. Amoxicillin and clavulanic acid were chosen for docking experiments as both are penicillins and blaA is a penicillinase. Intermolecular hydrogen bonds formed between the modeled structures of blaA variants and amoxicillin/clavulanic acid were determined. Pymol analysis revealed that each blaA variant had its own drug-interacting active-site amino acids. In all blaA variants S-135, A-137 and S-242 participated in hydrogen bonding with amoxicillin, and S-135 in hydrogen bonding with clavulanic acid. Kumar et al. 37 reported that the number of the residues involved in hydrogen bond interactions was an indicator of binding affinity of the enzyme with the ligand which might be decreased in mutant forms of class A blactamase. However, in our study no co-relation was observed between the number of residues involved in hydrogen bond interactions in blaA variants and their MIC for amoxicillin or clavulanic acid. Yi et al. 38 also observed that the cetadazidime MIC and the relative frequency of substitutions in PenA, a class A b-lactamase of Burkholderia thailandensis did not co-relate well. It may be worthwhile to mention that the Y. enterocolitica isolates used in the present work were isolated from environmental as well as clinical settings. Variations/mutations associated with their respective blaAs might have evolved over time either in nature or in clinical settings in patients treated with antibiotics. Salverda et al. 39 have also reported that most of the mutations that evolved in b-lactamases in clinical settings do not clearly co-relate with the stability or enhanced hydrolytic capability of the enzyme.
Variation in amino acid sequence is not the only factor that might contribute to differences in specific activities of b-lactamases. Enzyme activity may also be affected at the level of translation. In this regard, the secondary structure of mRNA plays an important role. When the secondary structures of mRNA of blaA variants were predicted using the web server mfold it was found that except for blaA4, all variants formed branched secondary structures. The mRNA secondary structure of blaA2 showed extensive branching and a lesser region of double stranded RNA helix compared to other blaA variants. On the other hand the mRNA secondary structure of blaA4 comprised mainly of double-stranded RNA helix of stacked base pairs with little branching. Several studies have shown that longer and stronger-paired zones in mRNA secondary structure tend to be more stable, preventing the ribosome from breaking the pairing and helping translation to proceed 40,41 . A lower translation rate means lesser number of copies of protein and hence a lower enzyme specific activity. Thus the secondary structures attained by mRNA of blaA2 and blaA4 might be the probable reasons for the higher and lower specific activities of their respective b-lactamases.
Overall the present work revealed that variations in promoter regions were not responsible for differential expression of blaA in biotype 1A strains of Y. enterocolitica. Analysis of amino acid sequences revealed that six types of blaA variants were present among biotype 1A strains of Y. enterocolitica. In-silico analysis of the blaA variants showed that amino acid substitutions did not contribute to the differential expression of the enzyme activities either. Rather the secondary structures attained by mRNA of blaA might underlie the differential expression of blaA in Y. enterocolitica. Further studies on secondary structures of mRNA of blaA gene will be helpful in better understanding the mechanisms that might underlie the differential expression of b-lactamases in Y. enterocolitica.

Methods
Bacterial strains. A total of 81 strains of Y. enterocolitica belonging to biotype 1A maintained in our laboratory (University of Delhi South Campus) were examined. Of these, sixty-five strains were isolated from India, which comprised 35 clinical, 18 wastewater, 7 swine and 5 pork isolates. All isolates were authenticated by, and have been deposited with Yersinia National Reference Laboratory and WHO Collaborating Center, Pasteur Institute, Paris (France). All these strains have also been deposited with the National repository i.e. Microbial Type Culture Collection (MTCC) and Gene Bank located at Institute of Microbial Technology, Chandigarh. Ten strains of Y. enterocolitica biotype 1A belonging to European and American origin were obtained from E. Carniel (Yersinia National Reference Laboratory and WHO Collaborating Center, Pasteur Institute, Paris, France). Another six biotype 1A strains of European origin were provided by J. Heesemann (Max von Pattenkofer Institute, Munich, Germany). The serovars, source of isolation, country of origin and reference laboratory accession numbers of these strains have been reported previously 22 . All strains were maintained on trypticase soy agar at 4uC. Spectrophotometeric assay of b-lactamases. Thirty-five strains of Y. enterocolitica biotype 1A were used for assaying b-lactamase activity as reported by Stock and Weidmann 18 with slight modification. A loopful of bacterial growth was inoculated in 10 ml of LB broth and incubated overnight at 28uC. The bacteria were sedimented by centrifugation at 4uC, washed thrice with chilled phosphate buffer (0.05 M, pH 7.0) and finally resuspended (1 g wet weight) in 5 ml of sonication buffer. The cells were disrupted by sonication on ice with three pulses of 30 sec each. The clear cell lysate, prepared by centrifugation at 12,000 g for 30 min at 4uC was used to assay the blactamase activity. The protein content of the lysates was estimated by Bradford's method 24 . b-lactamase activity was assessed spectrophotometrically by hydrolysis of nitrocefin. The assay mixture contained 110 ml nitrocefin (0.5 mg/ml), 30 mg of crude cell lysates containing b-lactamase in a final volume of 1.0 ml of 50 mM phosphate buffer (pH 7.0). b-lactamase activity was monitored by measuring decrease in the optical density at 486 nm for 10 min at 28uC. The specific activity of the enzyme was expressed as mmol of nitrocefin hydrolyzed/min/mg of protein and was calculated based on the molar extinction coefficient of 20,500 M 21 cm 21 for nitrocefin 25 .
Preparation of genomic DNA. Each strain was grown overnight in trypticase soy broth at 28uC. One ml of the culture was pelleted at 8,000 3 g for 10 min. The genomic DNA was prepared from each strain by using DNeasy Tissue kit (Qiagen, Hilden, Germany) with modifications for Gram-negative bacteria according to the manufacturer's recommendations. Purified DNA was eluted in sterile water and quantitated spectrophotometrically at 260 nm.
Primer designing for PCR amplification of blaA gene including the promoter and downstream region. For amplification of partial coding sequence of blaA including its promoter, primers A11f 59CAAATCGCCAGCAGAACAGA 39 and A12r 59GGGCAGCGAGCAGGGGTAAA 39 were designed by retrieving the whole genome sequence of Y. enterocolitica biotype 1B and the corresponding sequence of blaA along with the 350 bp each upstream and downstream of the gene. The primers were synthesized from Microsynth GmbH (Balgach, Germany).
PCR amplification of the promoter region and partial coding sequence of blaA was carried out for 81 strains of Y. enterocolitica in a Thermal Cycler (My Cycler TM , Bio-Rad, CA, USA). The PCR reaction mixture comprised of 13 PCR buffer (10 mM Tris-HCL, 1.5 mM MgCl 2 , 1.5 mM KCl and 0.1% Triton X-100), 200 mM of each of the four dNTPs (MBI Fermentas GmbH, Germany), 10 pmol each of forward and reverse primers, 2 U of Taq DNA polymerse (NewEngland Biolabs, MA, USA) and 50-100 ng of genomic DNA in a total volume of 25 ml. For PCR, the initial denaturation was performed at 95uC for 5 min followed by 25 consecutive cycles of 30 sec at 95uC, 30 sec at 56uC and 90 sec at 72uC and a final extension for 10 min at 72uC. The PCR products were analyzed by electrophoresis on 1% (w/v) agarose gel. The gels were stained with ethidium bromide (0.5 mg/ml) and visualized under UV transilluminator.
PCR amplification of blaA gene. PCR amplification of complete coding sequence (CCDS) of blaA was performed for 35 strains using published primers 15 in a Thermal Cycler (My Cycler TM , Bio-Rad, CA, USA). The PCR conditions were same as that for partial coding region of blaA and its flanking regions.
Sequencing of CCDS of blaA and its upstream and downstream regions. PCR amplicons were excised from gel, purified using QIA Quick Gel Extraction Kit (Qiagen, Hilden, Germany). The purified products were sequenced using the Big Dye Terminator Cycle Sequencing Ready Reaction kit and analyzed by capillary electrophoresis on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems, CA, U.S.A). The nucleotide sequences amplified using published primers 15 were aligned with those sequences which were amplified using primer pair A11f and A12r. The sequences which were common to the amplicons obtained using the two sets of primers and the extra sequences which were not common were used to construct the full sequence of blaA including the upstream and downstream sequences. CCDSs of blaA of 35 strains of Y. enterocolitica biotype 1A were translated into their corresponding amino acid sequences using the software expasy (www.expasy.org/translate). The amino acid sequences of blaA were aligned by Clustal . (www.ebi.uc.ak/clustal ).
Phylogenetic analysis of blaA of Y. enterocolitica strains. The amino acid sequences of blaA of 35 strains of Y. enterocolitica strains biotype 1A were also analyzed from evolutionary perspectives using the software MEGA5.2 (Molecular Evolutionary Genetic Analysis). The evolutionary history was inferred by using the Maximum Likelihood method based on the JTT matrix-based model (www.megasoftware.net).
Selection of Y. enterocolitica strains for modeling of blaA. To investigate the role of point substitutions in amino acid sequence in differential expression of blaA, analysis of protein models of blaA variants was carried out. After analysis by Clustal, amino acid sequences of blaA of Y. enterocolitica strains representing type/(s) of substitution/(s) were selected for modeling of their b-lactamases. The blaA of Y. strains 8, 19, 20, 29, E2 and E4 were selected for modeling and accordingly named as blaA1, blaA2, blaA3, blaA4, blaA5 and blaA6.
Homology modeling of b-lactamases. Since the three dimensional (3D) structure of blaA of Y. enterocolitica is not known, the 3D structures of b-lactamase variants blaA1, blaA2, blaA3, blaA4, blaA5 and blaA6 were predicted by homology modeling. The pairwise alignment between the target and template sequences was performed with PDB-BLAST. Three-dimensional structures of blaA1, blaA2, blaA3, blaA4, blaA5 and blaA6 were built using the software MODELLER 9.12 (http://salilab.org/ modeller/). Twenty models satisfying the spatial restraints based on target-template alignment were built for each protein, of which the 3D model with the lowest modeler objective function was selected. Following the modeling process, the structure was validated by PROCHECK and VERIFY3D. The PROCHECK program performs a detailed analysis of the stereochemistry of the modeled protein structure in the form of Ramachandran plot by plotting the phi-psi torsion angles. Population of phi-psi torsion angles in the plot gives an assessment of the overall quality of the modeled structure and also highlights regions that may need further investigation 26 . The quality of the three-dimensional models was also assessed by VERIFY3D, which tests the compatibility of a protein structure with its amino acid sequence. For highresolution, experimentally determined structures, the VERIFY3D scores are positive and consistently high i.e . 0.2 27 .
Molecular docking. Molecular docking analysis of the modeled protein structures was carried out to identify if substitutions in the amino acid sequences of blaA resulted in variations in their binding affinity to b-lactam antibiotic amoxicillin and the blactamase inhibitor clavulanic acid. The structures of amoxicillin and clavulanic acid were retrieved from the PubChem compound database (http://pubchem.ncbi.nlm.nih. gov/) 28 . Subsequently the downloaded SDF file was converted to PDB format using iBabel. Protein structures were pre-processed for docking using latest AutoDockTools (ADT) version by adding polar hydrogen atoms, removing all non-protein molecules from the modeled structure, adding Kollman charges to the protein and converting it to PDBQT format. Ligand structures were also prepared using ADT, using the default method for preparing ligands for docking that added hydrogen and charges. The default rotatable bonds were accepted as well, and the structure was converted to PDBQT format. The search space for docking was determined visually by centering the Grid Box to the known binding site of the ligand and expanding the dimensions of the cubic search space to completely encompass the site. Docking was performed using AutoDock Vina (http://vina.scripps.edu/index.html) with default parameter settings other than the search space specification described above, and the mean predicted binding affinity from the set of predicted binding poses was accepted as the true binding affinity for each docking run. The maximum number of poses per ligand was set to 10. The predicted binding of the ligand was compared to the known binding sites of blactamase in order to make predictions about differential responses with respect to each of the drugs. Interaction between ligand and target protein was analyzed by PyMOL 29 .
Analysis of mRNA secondary structure. The secondary structure of mRNA of blaA variants was predicted using the mfold webserver at default parameters (http://mfold. rna.albany.edu/). The mfold predicts energetically optimal secondary structure of an RNA molecule taking into account many realistic physical parameters that affect the RNA folding, such as pH, temperature and local biases of RNA.