The effects of NDM-5 on Escherichia coli and the screening of interacting proteins

Carbapenem-resistant Escherichia coli (E. coli) strains are widely distributed and spreading rapidly, creating significant challenges for clinical therapeutics. NDM-5, a novel mutant of New Delhi Metallo-β-Lactamase-1 (NDM-1), exhibits high hydrolase activity toward carbapenems. Since the genetic backgrounds of clinically isolated carbapenem-resistant E. coli are heterogeneous, it is difficult to accurately evaluate the impact of blaNDM–5 on antibiotic resistance. Herein, E. coli BL21 was transformed with a plasmid harboring blaNDM–5, and the resultant strain was named BL21 (pET-28a-blaNDM–5). Consistent with the findings of previous studies, the introduction of exogenous blaNDM–5 resulted in markedly greater resistance of E. coli to multiple β-lactam antibiotics. Compared with BL21 (pET-28a), BL21 (pET-28a-blaNDM–5) exhibited reduced motility but a significant increase in biofilm formation capacity. Furthermore, transcriptome sequencing was conducted to compare the transcriptional differences between BL21 (pET-28a) and BL21 (pET-28a-blaNDM–5). A total of 461 differentially expressed genes were identified, including those related to antibiotic resistance, such as genes associated with the active efflux system (yddA, mcbR and emrY), pili (csgC, csgF and fimD), biofilm formation (csgD, csgB and ecpR) and antioxidant processes (nuoG). Finally, the pGS21a plasmid harboring blaNDM–5 was transformed into E. coli Rosetta2, after which the expression of the NDM-5 protein was induced using isopropyl-β-D-thiogalactoside (IPTG). Using glutathione-S-transferase (GST) pull-down assays, total proteins from E. coli were scanned to screen out 82 proteins that potentially interacted with NDM-5. Our findings provide new insight into the identified proteins to identify potential antibiotic targets and design novel inhibitors of carbapenem-resistant bacteria.


Introduction
As an opportunistic pathogenic bacterium, Escherichia coli (E.coli) can cause severe infections in multiple local tissues and organs, such as the gastrointestinal tract and urinary tract, under certain conditions.The rapid dissemination of drug-resistant E. coli in the environment has placed extreme pressure on both livestock and human health Li et al. 10.3389/fmicb.2024.1328572(Steffan et al., 2020).Due to the widespread abuse of carbapenems, the emergence and transmission of carbapenem-resistant E. coli have become increasingly frequent.For instance, 34 of 100 urinary pathogenic E. coli (UPEC) strains isolated from clinical urine samples in Iran were found to be resistant to carbapenems (Khonsari et al., 2021).Carbapenem hydrolases, especially class B metalloenzymes, are mainly responsible for β-lactam resistance in Enterobacteriaceae.The different enzyme structures of these hydrolases can be divided into three categories, namely, B1, B2 and B3, with NDM, VIM, and IMP being the most important enzyme types (Behzadi et al., 2020;Bonnin et al., 2021).
Transcriptomics is a widely used high-throughput sequencing and analysis method that can be applied to extract comprehensive information about all transcripts in cells or tissues (Jovic et al., 2022).RNA-seq can help identify the critical factors involved in the emergence and dissemination of antibiotic resistance, facilitating the exploration of corresponding mechanisms and the screening of novel drug targets.For instance, RNA sequencing could be utilized to simultaneously evaluate the expression levels of thousands of genes under various circumstances, illustrating the mechanisms underlying the effects of multiple antibiotic resistance genes in bacteria and establishing a theoretical basis for designing novel antibiotics (Wang et al., 2009).Furthermore, Sinel et al. (2017) conducted RNA-seq analysis of Enterococcus plants treated with a subinhibitory concentration of daptomycin and detected alterations in the expression of 260 sRNA-encoding genes, 51% of which were correlated with carbohydrate metabolism and transport.
Interactions between bioactive proteins provide the foundation for cellular activity.Glutathione-S-transferase (GST) pull-down, which is a visual and fast method for detecting direct interactions between proteins, can be applied for the validation of interactions between known proteins and screening of unknown interacting candidates (Kim and Hakoshima, 2019).Mass spectrometry can be used to identify proteins via GST pull-down assays, playing a crucial role in the development of proteomics research.By combining GST pull-down with mass spectrometry, Lin et al. (2019) identified five proteins that interact with the porcine circovirus type 2 (PCV2) ORF4.Proteins interacting with Helicobacter pylori actin homologs (MreB) and influenza virus PB1-F2 were also successfully identified through a similar method (Zamarin et al., 2005;Zepeda Gurrola et al., 2017).
Carbapenem-resistant bacteria isolated from clinical environments are generally resistant to both carbapenems and other antibiotics.Therefore, the bla NDM−5 resistance gene was transformed into wild-type E. coli to eliminate interference by environmental factors.Transcriptome sequencing was performed to analyze alterations in the gene expression levels and functional pathways related to carbapenem resistance caused by bla NDM−5 .GST pull-down assays and mass spectrometry were used to identify proteins that interact with NDM-5 expressed in vitro, thus identifying potential drug targets for the treatment of carbapenem-resistant E. coli (Scheme 1).
Construction of E. coli BL21 (pET-28a-bla NDM−5 ).The bla NDM−5 DNA fragment was amplified from E. coli 120 by PCR.The primers used were designed with Primer Premier 5.0 according to the sequence of E. coli bla NDM−5 from GenBank (Gene ID: 2827929) (Supplementary Table 1; Zhao et al., 2023).After the amplified fragments and plasmid pET-28a were digested by HindIII/EcoRI and extracted by gel electrophoresis, the bla NDM−5 fragment and pET-28a were linked by Ligation Mix ligase.The recombinant plasmid was subsequently transformed into E. coli DH5α, after which the bacteria were cultivated on medium containing kanamycin (50 µg/mL) to construct positive recombinant plasmids.The selected plasmids were further validated by PCR and sequencing.Subsequently, the successfully constructed plasmid was transformed into BL21 competent cells by electroporation.Positive clones were identified on medium containing kanamycin (50 µg/mL) followed by PCR validation.The verified strain was named BL21 (pET-28a-bla NDM−5 ).Moreover, the plasmid pET-28a was electroporated into BL21 to obtain the control strain BL21 (pET-28a).
Antibiotic susceptibility assays.Single colonies of BL21, BL21 (pET-28a) and BL21 (pET-28a-bla NDM−5 ) were inoculated into 5 mL of cation-adjusted Mueller-Hinton broth (CAMHB) until the logarithmic growth period.The bacterial solution was adjusted to the 0.5 McFarland standard and diluted 100-fold.The minimum inhibitory concentrations (MICs) were determined by the micro broth dilution method according to the American Clinical and Laboratory Standards ("Clinical Lab Standards Institute.Performance standards for antimicrobial susceptibility testing."2016).The antibiotics were diluted twofold by CAMHB and mixed with an equal volume of bacterial suspension containing approximately 1.5 × 10 6 colony-forming units (CFUs) mL −1 in a clear ultraviolet (UV)-sterilized 96-well plate.After incubation at 37 • C for 18 h, the lowest concentrations of antibiotics with no visible growth of bacteria were determined as the MICs.Flowchart of the experiment.The biological characteristics of BL21 (PET-28a-bla NDM−5 ) were determined.Transcriptomic analysis and RT-qPCR were performed to identify differentially expressed genes in BL21 (PET-28a-bla NDM−5 ).The expression of the NDM-5 protein was induced, after which the NDM-5-interacting proteins were identified through a GST pull-down assay and mass spectrometry.
Growth curve and motility assays.One milliliter of bacterial liquid (optical density (OD) 600 = 0.2) was diluted 100-fold and inoculated in Luria-Bertani (LB) broth medium.The absorbance of each strain at a wavelength of 600 nm was monitored by a UV spectrophotometer every 2 h to construct a growth curve.
The motility of each strain was evaluated through the following method.Bacteria were cultivated in LB broth until the OD 600 reached 1.0.Then, 1 µL of the bacterial suspension was inoculated on the center of an LB semisolid culture plate (0.3% agar) and cultured at 37 • C for 16 h (Li et al., 2021).The diameter of the growth zone was used as an index to evaluate the movement ability of the bacteria.All the experiments were repeated three times.
Biofilm formation assays.200 µL of bacterial suspension (OD 600 = 0.2) was added to one well of a UV-sterilized 96-well plate.Each strain was plated in triplicate, and an equal amount of LB broth was added to the remaining wells as the control group.Then, the 96-well plates were placed in a wet box and incubated statically at 37 • C for 48 h.After the bacterial solutions had been discarded, each well was mixed with 200 µL of methanol for 15 min, followed by 3 washes with phosphate-buffered saline (PBS).Subsequently, each well was stained with 0.1% crystal violet, and then acetic acid was added to dissolve the precipitate.The OD 570 was measured by a multifunctional microplate reader (Tecan, Switzerland).All the experiments were repeated three times.
Scanning electron microscopy (SEM) analysis.The bacterial solution was adjusted to the 0.5 McFarland standard and diluted 100-fold.Then, each strain was cultured in LB broth until the logarithmic growth phase.The samples were washed with PBS and fixed in glutaraldehyde overnight at 4 • C. Finally, SEM was performed to observe the morphological changes in the bacteria.
Real-time fluorescent quantitative PCR (RT-qPCR).Seven differentially expressed genes (DEGs) related to antibiotic resistance were randomly selected from the transcriptome sequencing results, and the gene expression levels were detected via RT-qPCR.Briefly, total RNA was extracted with TRIzol reagent (Thermo Fisher Science, Carlsbad, Canada) and reverse transcribed to cDNA with a TransGen Reverse Transcription Kit (TransGen Biotech, Beijing).All the specific primers (Supplementary Table 2) used were designed according to the genome sequence of E. coli BL21 (NCBI reference sequence: CP020368.1).RT-qPCR was performed with the SYBR Green I method to analyze mRNA expression levels according to the instructions for the full gold quantitative PCR kit.16S rRNA was used as a housekeeping gene to standardize the expression levels, and the relative mRNA expression data were analyzed by the 2-Ct method (Li et al., 2023).Three replicates were performed for the RT-qPCR assay.
Expression of the NDM-5 protein and Western blot assays.The pGS21a plasmid was digested by KpnI/EcoR III and subsequently linked with the bla NDM−5 fragment by Ligation Mix ligase.The recombinant plasmid was transfected into E. coli Rosetta2.Positive colonies selected on LB medium containing 50 µg/mL ampicillin were identified by PCR and sequencing.E. coli Rosetta2 (pGS21a-bla NDM−5 ) was cultured until the OD 600 reached 0.6.After treatment with 0.5 mmol/L isopropyl-β-D-thiogalactoside (IPTG) for 4 h, 20 µL of each bacterial culture was centrifuged, collected, mixed with 5 µL of 5 × loading buffer, and then boiled at 100 • C for 10 min.The induced protein was purified following the instructions of the GST-Labeled Protein Purification Kit (Beyotime Biotechnology, Shanghai).Subsequently, the purified samples were separated via sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto a polyvinylidene difluoride (PVDF) membrane.Then, the membrane was blocked with 5% skim milk for 2 h at room temperature and further incubated with mouse GST monoclonal primary antibodies against GST-NDM-5 and horseradish peroxidase (HRP)-conjugated goat anti-mouse secondary antibodies.Finally, the membrane was visualized and recorded by a chemiluminescence imager.
GST pull-down assays.Approximately 25 µg of GST and GST-NDM-5 fusion proteins were immobilized in 50 µl of 50% glutathione resin suspension and equilibrated.Total protein from E. coli (25 µg) was added to the immobilized GST and GST-NDM-5.After being incubated at 4 • C for 4 h with gentle shaking, the mixture was centrifuged at 2000 rpm for 5 min.The precipitate was washed with PBS 3 times, after which the fusion proteins were eluted from the resin and analyzed via SDS-PAGE.

Dry lab section
Transcriptomic analysis.Single bacterial colonies were picked and cultured in LB broth overnight.A 1 mL bacterial suspension was inoculated into fresh LB broth at a ratio of 1:100 and cultured until the OD 600 reached 0.6.Then, the cell pellets were collected by centrifugation, and total RNA was extracted with TRIzol reagent (Thermo Fisher Scientific, Carlsbad, Canada).Subsequently, the integrity, purity and concentration of the extracted RNA were evaluated by an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., California, USA).The transcriptome library was constructed according to the instructions of the Illumina TruSeq RNA Sample Prep Kit (Thermo Fisher Science, Shanghai).After the library had been qualified and optimized to achieve standards, paired-end sequencing was performed on the Illumina NovaSeq 6000 platform.The raw paired-end reads were trimmed and quality controlled by Trimmomatic.Then, the clean reads were separately aligned to the reference genome in orientation mode using Rockhopper 1 software.If the proportion of clean reads was greater than 70%, EdgeR software was used to analyze the difference in the expression level of each transcript between samples.The DEGs were analyzed according to the fragments per kilobase of exon per million fragments mapped (FPKM) value (Behzadi and Ranjbar, 2019).Goatools 2 was used for Gene Ontology (GO) enrichment analysis, while KOBAS 3 was utilized for Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis.

Results
3.1 Successful construction of E. coli BL21 (pET-28a-bla NDM−5 ) Due to the complex genetic background of clinically isolated bla NDM−5 -positive E. coli, the plasmid pET-28a containing the bla NDM−5 fragment was transformed into wild-type E. coli BL21 to exclude potential interference.After confirming the existence of bla NDM−5 in E. coli strain 120 by PCR (Figure 1A) and sequencing (Supplementary Table 3), the bla NDM−5 fragment was amplified and linked with plasmid pET-28a by double enzyme digestion.As shown in Figures 1B, C, the recombinant plasmid pET-28abla NDM−5 was successfully constructed.The sequencing data for   1D), confirming the successful construction of E. coli BL21 (pET-28a-bla NDM−5 ).
3.2 Increased resistance to β-lactam antibiotics in BL21 (pET-28a-bla NDM−5 ) The MICs of meropenem and imipenem against all the bacterial strains are available in Table 1.Transformation with the vector pET-28a did not change the susceptibility of BL21 to meropenem or imipenem.In contrast, the MIC of meropenem increased 520fold (MIC = 16 µg/mL) and that of imipenem increased 128-fold (MIC = 32 µg/mL) for bla NDM−5 -positive BL21.Subsequently, the susceptibility of BL21 (pET-28a) and BL21 (pET-28a-bla NDM−5 ) to 14 commonly used antibiotics was evaluated (Table 2).Compared with that for BL21 (pET-28a), for BL21 (pET-28a-bla NDM−5 ), the MIC of ceftazidime increased 128-fold and that of amoxicillin increased 32-fold, but those of colistin and polymyxin B decreased 2-fold.In addition, the susceptibility to non-β-lactam drugs such as ofloxacin and gentamicin did not significantly change after the introduction of the bla NDM−5 resistance gene.The fold changes in MICs are shown in Figure 2A.
As shown in Figure 2B, there was no significant difference in growth rate among the three bacterial strains (p > 0.05), demonstrating that transformation with the bla NDM−5 recombinant plasmid did not affect the growth rate of E. coli.The biofilm formation ability of the strains was subsequently evaluated via the crystal violet staining method.Compared with that of BL21 and BL21 (pET-28a), the biofilm formation capacity of BL21 (pET-28a-bla NDM−5 ) increased markedly (p < 0.01) (Figure 2C).In contrast, a significant reduction in the growth zone diameter was observed for BL21 (pET-28a-bla NDM−5 ), indicating that the acquisition of carbapenem resistance diminishes the motility of E. coli BL21 (p < 0.05) (Figure 2D).As shown in Figure 2E, bla NDM−5 had little influence on the morphology of the bacteria.SEM revealed that only an extremely tiny fraction of the BL21 cells (pET-28a-bla NDM−5 ) exhibited slight deformation of cell morphology.

Results of transcriptome analysis
DEGs were screened out as genes with | log2 (fold change)| > 1 and Q-value < 0.005.A total of 665 DEGs were detected in the BL21 (pET-28a)/BL21 group, including 414 upregulated genes and 251 downregulated genes.Moreover, 461 DEGs, including 194 upregulated genes and 267 downregulated genes, were detected in the BL21 (pET-28a)/BL21 (pET-28a-bla NDM−5 ) group (Figure 3A).The distribution of DEGs in the BL21 (pET-28a)/BL21 (pET-28abla NDM−5 ) group is shown in Figure 3B.The genes related to antibiotic resistance were selected from the identified DEGs, and some of them are listed in Supplementary Table 5.These genes included active efflux system genes (yddA, mcbR and emrY), pilus genes (csgC, csgF and fimD), biofilm formation genes (csgD, csgB and ecpR) and antioxidant process-related genes (nuoG).GO enrichment classification was used to annotate the functional associations of the genes.A total of 307 DEGs from the BL21 (pET-28a-bla NDM−5 )/BL21 (pET-28a) group were annotated into three categories: biological process (BP), cellular component (CC), and molecular function (MF).As shown in Figure 3C, DEGs in the biological process category were mainly enriched in the cellular process, metabolic process, response to stimulus, biological regulation, and cellular localization terms.In the cellular component category, DEGs were enriched in terms such as intracellular, protein-containing complexes and cellular antioxidant entities.In addition, in the molecular function category, terms related to catalytic activity, transport activity, structural molecule activity and binding reactions were foremost highly enriched by the DEGs.Subsequently, the KEGG database was used to annotate the DEGs into 99 pathways, and they were chiefly enriched in ribosome, oxidative phosphorylation, RNA polymerase, butyrate metabolism, lysine degradation, the tricarboxylic acid cycle, alanine metabolism and other functional pathways (Figure 3D).

82 NDM-5-interacting proteins were screened out and analyzed via mass spectrometry
Double enzyme-digested fragments of pGS21a-bla NDM−5 were verified by gel electrophoresis (Figure 4A).The sequencing results are available in Supplementary Table 6.After pGS21a-bla NDM−5 was transformed into E. coli Rosetta2, IPTG was used to induce protein expression.As expected, SDS-PAGE analysis detected a specific band at 54 kDa, indicating successful expression of the NDM-5 protein (Figure 4B).Subsequently, the NDM-5 protein was purified and linked with GST labels.As shown in Figure 4C, a single and distinct band could be observed at 54.5 kDa by SDS-PAGE, which is consistent with the results of Western blot analysis.
GST pull-down assays were performed for GST-NDM-5 and total E. coli proteins to investigate the proteins potentially interacting with NDM-5.The purified proteins were verified by SDS-PAGE.In addition to the NDM-5 band at 54.5 kDa, other protein bands were also observed, demonstrating the presence of NDM-5-interacting proteins (Figure 4D).The interacting proteins were extracted and purified, followed by mass spectrometry identification.A total of 665 interacting proteins were identified in Frontiers in Microbiology 06 frontiersin.orgthe GST + /E. coli + control group, and 466 potential interacting proteins were identified in the GST-NDM-5 + /E. coli + interacting group.Since 384 proteins were detected in both groups, these proteins were removed from the GST-NDM-5 + /E. coli + group, leaving the remaining 82 proteins as candidates that potentially interact with NDM-5 (Supplementary Figure 1).Some of the identified NDM-5-interacting proteins are listed in Supplementary Table 7. GO and KEGG analyzes were performed to elucidate the biological function of the 82 interacting proteins.A total of 307 GO-enriched functional pathways were annotated.Among them, 88 pathways were involved in molecular functions, mainly binding, catalytic and transport activities.In addition, 100 pathways were correlated with cellular components, including cells, cellular components and macromolecular complexes.The remaining 119 pathways were related to biological processes, such as metabolic processes, cellular processes, localization processes, biological regulation processes and responses to stimuli (Figure 5A).Moreover, 90 enriched functional pathways were annotated by KEGG analysis, including 65 metabolic pathways, 16 environmental information processing pathways, and 3 genetic information processing pathways (Figure 5B).

Discussion
Carbapenem resistant E. coli strains usually contain other resistance genes.Therefore, herein, we used the pET-28a plasmid to transform bla NDM−5 into wild-type E. coli to exclude interference by other resistance genes.Compared with the wild type, the BL21 (pET-28a-bla NDM−5 ) strain exhibited markedly increased resistance to meropenem (520-fold) and imipenem (128-fold) coupled with markedly increased resistance to other β-lactam antibiotics, such as amoxicillin and ceftazidime.Meanwhile, the MICs of non-β-lactam antibiotics did not change, indicating that NDM-5 cannot hydrolyze non-β-lactam drugs.However, twofold Frontiers in Microbiology 08 frontiersin.orgdecreases in the MICs of colistin and polymyxin B were observed after the introduction of bla NDM−5 .This increase in susceptibility could be attributed to the downregulation of eptB expression since the chromosomally encoded PEA (phosphatidyl-ethanolamine) transferases (e.g., EptA, EptB and EptC) can mediate polymyxin resistance in gram-negative bacteria through the modification of Lipid A (Huang et al., 2018; Supplementary Table 5).
Our results indicated that the transformation of pET-28a-bla NDM−5 marginally reduced the growth rate of BL21, significantly restricted bacterial motility and strongly enhanced biofilm formation capacity.It is speculated that the acquisition of carbapenem resistance is probably caused by these phenotypic changes.For instance, compared with those of susceptible strains, the motility and growth rate of polymyxin-resistant bacteria are significantly lower (Tietgen et al., 2019;Yang et al., 2020).In addition, biofilms can prevent antibiotics from penetrating cell membranes and entering the interior of cells, thus enhancing bacterial resistance to multiple antibiotics (Romeu et al., 2020).Lee et al. (2020) transformed plasmids harboring bla NDM−1 and bla OXA−232 into E. coli DH5α and reported that both transformants had a significantly enhanced biofilm formation capacity.
The transcriptome analysis results indicated that the expression of biofilm-related genes (ycgZ, csgD and mcbR) was significantly upregulated, which is consistent with the observed bacterial phenotypes.Therefore, it is speculated that csgD and mcbR increase bacterial resistance through the modulation of motility and biofilm formation.The CSG region encodes pilus-related genes, while CsgD controls cell aggregation and inhibits flagellar synthesis through direct regulation of Curli pili (Herrington and Sitaras, 2013;Behzadi and Behzadi, 2017;Ogasawara et al., 2020).In addition, McbR can regulate the formation of biofilms and mucus by inhibiting the expression of the periplasmic protein YbiM and binding to the yciG promoter to activate the transcription of the yciGFE operon in E. coli K-12 (Yu et al., 2019).We also identified several efflux pump-related genes, including yddA, mcbR, emrY and mdtl.The efflux pump can expel intracellular antibiotics to reduce the effective concentration of drugs, thus conferring antibiotic resistance to bacteria (Lee et al., 2017).As a prokaryotic homolog of P-gp, YddA is a representative member of the superfamily of ATP-binding cassette (ABC) transporters that can recognize and transport a variety of substrates (Feng et al., 2020).In avian pathogenic E. coli, an electrophoretic mobility shift assay (EMSA) was used to confirm that McbR could specifically bind to the promoter regions of acrAB, acrD, acrR, emrD and mdtD, revealing that McbR modulates bacterial resistance by directly activating the transcription of efflux pump genes (Yu et al., 2020).Furthermore, the amino acid sequence of EmrY shares 63.3% homology with that of EmrB, and these proteins together constitute a multidrug efflux pump (Tanabe et al., 1997;Wei et al., 2022).The RNA-seq results verified the upregulated expression of such efflux pump genes, partially explaining the correlation between bla NDM−5 and carbapenem resistance.pET-28a is a plasmid encoding an RNA polymerase, an ori gene, a lacI gene, and an antimicrobial gene.Accordingly, protein expression of BL21 (pET-28a) and BL21 alone should be similar.However, a total of 665 DEGs were detected in the BL21 (pET-28a)/BL21 group, which is quite different from the expected results and is one of the limitations of this study.
Due to the high solubility, distinct affinity and high-efficiency expression of GST tags, GST fusion proteins are commonly utilized for the identification of interacting proteins (Pan et al., 2011).E. coli grows best at 37 • C, but at that temperature, the expression of specific proteins may not reach the optimal level due to "leakage" of some expression vectors and metabolic pressure on bacterial cells (Sørensen and Mortensen, 2005).Low temperature can slow the synthesis rate of proteins and alter the dynamics of peptide folding to increase the possibility of correct folding, thus facilitating the soluble expression of fusion proteins (Donovan et al., 2000;Sadraeian et al., 2013).In this paper, we compared the concentration ratios of soluble proteins at 15 • C and 37 • C and determined that 15 • C was the optimal incubation temperature for the expression of NDM-5.Subsequently, 82 proteins that interact with NDM-5 were identified by GST pulldown assays.Among them, the transcriptional regulator AraC, the multidrug efflux RND transporter permease AcrD, and the multidrug efflux transporter outer membrane subunit MdtP are closely associated with the function of efflux pumps.In gramnegative bacteria, the overexpression of the AcrAB, AcrD and MdtABC efflux pumps can cause antibiotic resistance (Yu and Schneiders, 2012).Therefore, the upregulated expression of efflux pump-related genes detected by transcriptome sequencing is presumably caused by the interaction between the NDM-5 protein and efflux pump-related proteins.The detected histidine kinase precursor ZraS and histidine kinase YpdA are two-component system sensors that participate in various cellular processes, including antibiotic resistance, quorum sensing and osmotic pressure (Zschiedrich et al., 2016).It is speculated that NDM-5 could interact with a two-component system to regulate the expression of antibiotic resistance genes, facilitating the recognition of the host environment and adaptation to adverse environments.ABC transporters are a family of proteins with highly conserved transmembrane structures that transport substances through the energy generated by ATP hydrolysis (Thomas et al., 2020).The screened proteins related to ABC transporters, such as A0A140N5T4, A0A140NBB8 and A0A140NFJ3, may interact with NDM-5 to participate in the regulation of bacterial resistance.The major limitations of the present study are that the proteins potentially interacting with NDM-5 were screened, and the mechanism of interaction could not be accurately investigated (Figure 5C).Considerably more work is needed to identify proteins that interact with NDM-5, with the aim of providing a theoretical basis for exploring new targets for antimicrobial drugs, identifying new inhibitors against carbapenem-resistant E. coli, and preventing and controlling superresistant bacteria.Although the potential mechanisms through which NDM-5 interacts with proteins need further investigation, the identified interacting proteins provide an outstanding foundation for exploring the function of NDM-5 in bacteria.

Conclusion
In this study, we used the wild-type strain of E. coli and the plasmid pET-28a to construct a bla NDM−5− positive strain.The antibiotic susceptibility and biological characteristics of E. coli BL21 (pET-28a-bla NDM−5 ) were evaluated, revealing that bla NDM−5 has a significant impact on the bacterial phenotype.Transcriptome analysis and RT-qPCR demonstrated that bla NDM−5 could cooperate with other antibiotic resistance genes (yddA, oppD, mcbR, emrY, etc.) to increase bacterial resistance to various βlactam antibiotics, especially imipenem and meropenem.We successfully expressed the soluble NDM-5 protein in vitro and screened 82 proteins that interact with NDM-5.In conclusion, our study could provide a theoretical basis for the exploration of novel drug targets, but further research focused on the interaction mechanism between specific proteins and NDM-5 is needed.

FIGURE 1
FIGURE 1 Identification of the constructed strains.(A) PCR amplification results of the NDM-5 gene.(B) PCR amplification results of the recombinant plasmid.(C) Construction and detection results of the pET-28a-bla NDM−5 recombinant plasmid.(D) PCR verification of the BL21 (pET-28a-bla NDM−5 ) strain.

FIGURE 2
FIGURE 2Evaluation of biological characteristics.(A) Antimicrobial susceptibility of the bacterial strains.(B) Growth curves of the strains.The growth rate was evaluated based on the OD 600 value of the bacteria after2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22  and 24h of incubation.(C) The biofilm formation ability of the strains.The OD 570 value of the bacteria in 96-well cell culture plates measured by an enzymometer was used as an index to evaluate the biofilm formation ability of the bacteria.**p < 0.01.(D) Mobility of the strains.The bacterial colony diameter after 16h of incubation was used as an index to evaluate the motility of the strain.*p < 0.05.(E) Observation of bacterial morphology by scanning electron microscopy.Red arrows: abnormal changes.
and 24h of incubation.(C) The biofilm formation ability of the strains.The OD 570 value of the bacteria in 96-well cell culture plates measured by an enzymometer was used as an index to evaluate the biofilm formation ability of the bacteria.**p < 0.01.(D) Mobility of the strains.The bacterial colony diameter after 16h of incubation was used as an index to evaluate the motility of the strain.*p < 0.05.(E) Observation of bacterial morphology by scanning electron microscopy.Red arrows: abnormal changes.

FIGURE 3
FIGURE 3 Results of transcriptome analysis.(A) Numbers of differentially expressed genes between groups.(B) Volcano map of differentially expressed genes in BL21 (pET-28a)/BL21 (pET-28a-bla NDM−5 ) cells.(C) GO classification of differentially expressed genes.(D) KEGG functional enrichment of the differentially expressed genes.The number of DEGs is proportional to the size of the dots.(E) RT-qPCR results.

FIGURE 5
FIGURE 5 Mass spectrometry analysis results.(A) GO annotation of potential NDM-5-interacting proteins.(B) KEGG annotation of potential NDM-5-interacting proteins.(C) Summary of the effects of NDM-5 on E. coli and the screening of interacting proteins.

TABLE 1
Results of the sensitivity of the test strains to carbapenem drugs.