A deep learning training set is established for novel anti-H. pylori agents exploration
First, the dataset was curated from reputable sources and ensured the diversity in chemical structures and activity levels. A sizable collection of 938 compounds with known anti-H. pylori properties was established. This dataset included 801 reported anti-H. pylori compounds with structural diversity from ChEMBL database,29 as well as 137 self-established BBR derivatives.30,31,32 An MIC value of 16g/mL was set as the critical value. The compounds with MICs 16g/mL were defined as active (label 1) and MICs > 16g/mL as inactive (label 0). The proposed deep learning framework firstly represented compounds with molecular graph, and extracted the molecular extended-connectivity fingerprints (ECFP)33 which preserve rich functional group information. Feature engineering was performed to extract the ECFP that captured essential functional group information, and leveraged message passing deep neural network to extract properties directly from molecular structure.
Since the significant interactions between atomic pairs with topologically distant could also affect the overall molecular properties (Fig. 1a), a deep graph neural network (Attentive FP)34 was applied to learn the embeddings of molecular graph, including both local and nonlocal features of the molecular structures. More specifically, every compound was represented with molecular graph, where nodes denoted atoms, and edges denoted bonds (Fig. 1a). By leveraging RDKit and DGL-LifeSci packages, vectors with a length of 39 for nodes and 11 for edges were obtained to represent the chemical properties of atoms and bonds, respectively. Attentive FP was used to translate the molecular graph with node and edge features into a continuous vector, which was the compound representation. Attentive FP iteratively aggregated the features of atoms and bonds with graph attention network (GAT)35 in the messaging phases, which allowed an atom to focus on most relevant neighborhoods. Then, it retained and filtered information with a gated recurrent network unit (GRU)36 in the readout phases, which allowed the model to capture the implicit effects among distant atoms. After obtaining the molecular graph representation, an attention mechanism to self-adaptively integrate molecular graph representation and ECFP fingerprints was introduced.
Establishment of the deep learning model. a Deep learning-based anti-H. pylori compound discovery. SMILES simplified molecular input line entry system. b A pie chart for data distributions, including three pre-train sets, a fine-tune set and a test set. c ROC-AUC plot evaluating model performance under the ten-fold cross-validation. d t-Distributed stochastic neighbor embedding (t-SNE) of all molecules from the pre-training, fine tune, and test set, revealing chemical relationships between these compounds
Considering that the 938 compounds with known anti-H. pylori properties were insufficient for training a successful deep learning model, we utilized the pre-train-then-fine-tuning paradigm,37 which pre-trained the deep learning model on large-scale bioassays related to H. pylori from PubChem databas,38 and fine-tuned the pre-trained deep learning model on the collection of 938 compounds. The pre-train database included 8999, 892, and 2809 compounds (Fig. 1b), respectively. All the above-mentioned training set information was provided as supplementary data sets. In the fine-tune phase, the parameters of the nonlinear multilayer perceptron network (MLP) in the pre-trained deep learning model were initialized and the model was further optimized on the collection of 938 compounds for capturing task-specific patterns. Finally, the molecular fingerprint features and molecular graph embeddings were self-adaptively integrated to form the compound feature vectors and then an MLP layer39 was leveraged to predict their activity against H. pylori.
The predictive accuracy of the model was assessed through ten-fold cross-validation on the training dataset and external validation on the independent dataset. Cross-validation techniques were applied to validate the robustness and reliability of the model.40 The performance of our final model was quantified as follows: the area under the receiver operating characteristic curve (ROC-AUC) attained a value of 0.9033, signifying commendable discriminative capacity; the area under precision-recall curve (AUPR) registered at 0.9615, indicating a robust precision-recall balance. Moreover, the F1-score, a composite metric denoting the harmonious interplay of precision and recall, manifested at 0.8797, attesting to a noteworthy equilibrium between these facets. The model also attained an accuracy rate of 0.8326, representing the proportion of accurately classified instances. Furthermore, the recall, an indicator of the models ability to correctly identify actual positives, attained a value of 0.8454, while the precision, signifying the proportion of predicted positives correctly classified, was recorded at 0.9169. These metrics collectively corroborated the models effectiveness in addressing specific classification tasks within the ambit of H. pylori inhibition.
Thus, this established deep learning model enabled the establishment of a correlation between the structural characteristics of these compounds and their antibacterial activity against H. pylori. To validate the effectiveness of this deep learning model, a series of novel BBR derivatives were strategically designed for prediction.
It is reported that modifications on the D-ring of BBR/PMT (Fig. 2a), such as 9-position mono-substitution, have limited enhancements of anti-H. pylori activity.32 While modifications were conducted on 13-position of ring C (Fig. 2a), the corresponding derivatives only exhibited moderate anti-H. pylori potencies.41 Meanwhile, there is scarce literature reporting on the anti-H. pylori activity of A-ring modified derivatives, making them highly attractive for novel anti-H. pylori drug discovery utilizing deep learning models. Considering the synthetic accessibility, we selectively chose 3-position of the A-ring for modifications with various types of substituents, including chain alkanes, cycloalkanes and substituted phenyls. Thus, a set of 3-substituted novel BBR/PMT derivatives was virtually designed for prediction. Two of them (5 and 6) were positively predicted and the rest nine were predicted to be negative (14, 913). To verify the accuracy and reliability of the deep learning model employed, all designed compounds were synthesized through an easy-to-operate one-step synthetic procedure as shown in Supplementary Scheme 1, and subsequently subjected to the antibacterial activity evaluation. Simultaneously, two 3,13-disubstituted derivatives (7 and 8) were accidentally obtained and identified during the synthesis of 5 and 6, respectively, with the existence of an excessive -C containing electrophilic reagent. Compared to previously reported procedures involving more than three steps,42 the disubstituted derivatives could be obtained with satisfactory yields ranging from 6167%. These two compounds were also predicted to be positive (78).
In vivo antibacterial evaluations for compound 8. a Chemical structures of BBR and 8. b Serum biochemical indices of liver and kidney functions for mice in different treatment groups (n=6). c Plasma and stomach concentrationtime profiles of 8 following a single oral dose of 30mg/kg (n=4). d The schematic diagram of H. pylori infection and treatment process in C57BL/6 mice. e, g The viable counts in the stomach of mice infected with H. pylori CCPM(A)-P-3722159 in each group (n=5) after different treatments. The administration dosage of each treatment component is as follows: OPZ (200g/kg); 8 (30mg/kg); AMX (15mg/kg); CLA (15mg/kg); CMC, carboxymethyl cellulose; AC, AMX+CLA; Bi, bismuth citrate (5mg/kg). f Hematoxylin and eosin (H&E) staining of stomach tissues
All constructed BBR/PMT derivatives were first evaluated for their activity against six different H. pylori strains, including two American Type Culture Collection (ATCC) reference strains of ATCC 43504 and ATCC 700392, and other four clinical isolates, taking BBR, PMT, CLA, AMX, LEV, and MTZ as positive controls. The tested strains included CLA-resistant strains (CCPM(A)-P-3716289 and CCPM(A)-P-3716370), MTZ-resistant strains (ATCC 43504 and CCPM(A)-P-3716289), LEV-resistant strains (CCPM(A)-P-3716289 and CCPM(A)-P-2316370), and an AMX-resistant strain (SS1). The chemical structures of the designed compounds, the deep learning prediction results, and their MIC values against the tested H. pylori strains are listed in Table 1. The results demonstrate a notable degree of predictive success, as evidenced by the MIC values. Specifically, the positively predicted compounds (58) exhibited substantially lower MIC values, ranging from 0.258g/mL. In contrast, for the negatively predicted compounds (14, 913), the MIC values went up to a range of 16 to >256g/mL. Therefore, compounds 5, 7, and 8 with the best antibacterial potencies were selected as representative compounds for further investigation. This approach exemplifies a judicious combination of computational prediction through deep learning models and experimental validation, constituting a powerful strategy for candidate exploration in future anti-H. pylori drug development.
The effects of predicted hits 5, 7, and 8 on cell viability were evaluated using the MTT assay in gastric epithelial cells (GES-1), hepatocellular carcinoma (HepG2), human non-small lung cancer (H460) and human embryonic kidney (293T) cells. The cell viability was determined after the exposure to varying concentrations of these compounds. As presented in Supplementary Table S1, compound 8 (Fig. 2a) exhibited lower cytotoxicity with the median toxic concentration (TC50) values ranging from 50.59 to 57.07M, compared to those of 5 (17.6824.96M) and 7 (8.8112.70M). Compound 8 exhibited the best anti-H. pylori activity and the lowest cytotoxicity, as well as the most favorable therapeutic index. Therefore, it was selected as a potential candidate for further studies.
The acute oral toxicity test of compound 8 was conducted in Kunming mice. The mice were closely monitored for 14 days, and the medium lethal dose (LD50) value of 8 was over 500mg/kg, which indicated a satisfactory safety profile of 8 for oral administration. Then, the blood samples collected from the above mice were assessed for the biochemical indices of liver and kidney functions. As illustrated in Fig. 2b, 8 did not lead to obviously elevation of glutamic oxalacetic transaminase (GOT), glutamic pyruvic transaminase (GPT), blood urea nitrogen (BUN) or creatine (CRE), indicating no detectable adverse effect of 8 on liver or kidney function.
To explore the pharmacokinetic profile of compound 8, the stomachs and plasma of C57BL/6 mice were collected and detected at different time points after a single oral dose of 30mg/kg. As illustrated in Fig. 2c, the gastric concentrations of 8 maintained above its MIC value (0.5g/mL) after 24h (3.251.51g/g, Supplementary Table S2), indicating an ideal gastric retention that could ensure its anti-H. pylori efficacy in vivo. Meanwhile, the maximum concentration (Cmax) of 8 in plasma was below 0.1g/mL (Supplementary Table S3), and it became undetectable (below the detection limit of 0.001g/mL) after 6h, suggesting a low possibility of systemic side effects. Besides, the acid stability of 8 was also assessed under the pH values of 1.0 and 3.0 (to simulate the acidic environment in gastric acid), at different time points (2, 8, and 24h). As shown in Supplementary Table S4, the content of 8 was still above 90% after 24h treatment in the acidic environment. Taken together, the favorable acid stability and pharmacokinetic properties of 8, including bare absorption to system circulation and long gastrointestinal retention, make it suitable for being developed as an anti-H. pylori agent.
Twenty-seven clinically isolated H. pylori strains were employed for further potency evaluation of 8. As shown in Table 2, compound 8 exhibited a robust activity with an MIC of 0.5g/mL against all tested strains (14 CLA-resistant strains, 11 MET-resistant strains, 10 LEV-resistant strains, 2 AMX-resistant strains, and 6 MDR strains, and all the resistant information is highlighted in dark color in Table 2).
Compound 8 was then challenged over a 36-day serial passage assay to determine the rate of potential resistance induction on H. pylori ATCC 43504, which is susceptible to CLA and AMX originally. As shown in Supplementary Fig. S1, repeated exposure to sub-MIC level of 8 or AMX did not develop resistance in the tested H. pylori strain by serial passage (12 passages). After 12 passages under permanent selective pressure of CLA, the bacteria showed resistance to CLA with the MIC reaching and stabilizing at 4g/mL (256-fold of initial MIC).
Checkerboard assay was performed to test the combined effects of 8 and AMX or CLA. As displayed in Supplementary Table S5, when combined with CLA, synergistic effects (fractional inhibitory concentration index, FICI0.5) could be observed in 10 out of 25 tested strains (5 out of 9 CLA-resistant strains) with the FICI values of 0.1880.50. Meanwhile, only additive effect (0.5 The in vivo antibacterial activity of compound 8 was evaluated in the C57BL/6 mouse gastric infection model (Fig. 2d). The mice were first randomly assigned into five groups: an uninfected control group and four infected groups with different treatments, which included a vehicle carboxymethyl cellulose (CMC) control group, dual therapy group (OPZ plus 8 [OPZ+8]), triple therapy group (OPZ plus AMX and CLA [OPZ+AC]), and quadruple therapy group (OPZ plus AMX, CLA, and 8 [OPZ+AC+8]), respectively. The mice in the infected groups were orally administrated via gavage with H. pylori CCPM(A)-P-3722159, a mouse-adapted MDR strain (resistant to AMX, CLA, and LEV), every other day for four times. After a two-week colonization period, the different treatments were performed as above for five consecutive days. The therapeutic efficacy was evaluated by comparing the viable bacteria counts in the mouse stomachs. As shown in Fig. 2e, treatment with OPZ+8 (30mg/kg) significantly decreased the gastric bacteria load of the infected mice from 1.3105 to 6.5102 CFU/g (2.2-log reduction in comparison to CMC group), which was similar to that of the triple-therapy group (OPZ+AC, 1.8-log reduction in bacterial burden). Remarkably, the quadruple-therapy treatment (OPZ+AC+8) further decreased the bacteria load to 2.0102 CFU/g (2.8-log reduction), representing a 99.8% inhibition of stomach colonization compared with CMC group. These results suggest that, with the pretreatment of OPZ, 8 could exert comparable eradicative efficacy to the combination of OPZ, AMX and CLA in vivo, and exhibited improved activity when combined with AMX and CLA, thereby increasing the clearance of the colonized multidrug-resistant H. pylori. Additionally, there was no significant body weight loss after the different treatment, as shown in Supplementary Fig. S2. Histopathological examination of fixed stomach sections revealed that H. pylori infection led to a more porous and bloated structure of the gastric gland, the obvious inflammatory infiltration, and the increase of pepsinogen (high pepsinogen usually related to H. pylori infection, peptic ulcer, and gastritis) compared with the uninfected tissue (Fig. 2f). The dual, triple, and quadruple-therapy treatments alleviated the gastric inflammation in some degree and decreased the level of pepsinogen, indicating the eradication of the pathogens. The long-term use of antibiotics often leads to a disturbance of the intestinal flora and a decrease in gut microbiota diversity. To investigate whether 8 affects the gut microbiota, stool samples were collected from each group, and 16S rRNA gene sequencing was employed to analyze the gut microbiota constitution. The Venn graph was used to analyze the characteristic sequence numbers of each group. As shown in Fig. 3a, the largest number of same specific characteristic sequences between 8 treatment group (T8: OPZ+8) and the uninfected group were observed, compared with other comparisons. Using alpha diversity (Pieloi_e) analysis (Fig. 3b), the microbiota diversity in the vehicle control group (CMC) and triple therapy group (OPZ+AC) was found to be significantly decreased compared with that in the uninfected group at the genus level. It is worth noting that, the box diagram showed that the diversity of intestinal flora of mice in group T8 was close to that in the healthy group (p>0.05), higher than that in CMC group and group OPZ+AC. Principal coordinate analysis (PCoA) showed that, in comparison to CMC and OPZ+AC groups, the composition of intestinal flora of the T8 group exhibited more similarity to the uninfected group (Fig. 3c). Gut microbiome analysis in different treatment groups (n=5). Uninfect, the uninfected group; CMC, vehicle control group; T8, dual therapy group (OPZ+8); AC, triple therapy group (OPZ+AC); AC8, quadruple therapy group (OPZ+AC+8). a The Venn diagram of microbial characteristic sequences of each treatment group. b Alpha diversity analysis on microbiota diversity of each treatment group. c Beta diversity of PCoA analysis. d A bar plot analysis at the genus level (ten bacterial genera with the highest abundance). e A heatmap analysis at the genus level (ten bacterial genera with the highest abundance). f LDA value distribution histogram revealed by LEfSe software. When species with LDA Score >4 are statistically different, the length of the histogram (LDA Score) represent the impact size of the different species. g Evolutionary branching trees from the inside out in a clade represent the level of phylum, class, order, family, genus Next, a bar plot and a heat map analysis at the genus level were performed to show the ten bacterial genera with the highest abundance of each treatment group (Fig. 3d, e). Through the relative abundance analysis at genus level, the intestinal flora disorder was observed in AC group, with the overgrowth of several genera, including Klebsiella, Escherichia-Shigella, and Bacteroides. In contrast to AC group, the microbiota constitution of the dual therapy group (T8) was sustained and the abundance of probiotics, including Lactobacillus and Dubosiella was partially restored. The bacterial genera with the highest abundance in each mouse was also displayed in Supplementary Fig. S3. In addition, Bifidobacterium, another kind of well-known probiotics (not belonging to ten highest abundance), was also significantly enriched in the dual therapy group compared with AC group (Supplementary Fig. S4), confirming that 8 has the tendency to avoid dysbiosis of intestinal flora. To further display the observed differences in the microbiome composition, linear discriminant analysis (LDA) effect size (LEfSe) analysis (Fig. 3f) was performed, and the cladogram was generated based on LEfSe analysis (Fig. 3g). Consistent with the above results, there was a significant increase in the abundance of Lactobacillus (LDA (log10)>4.0, p<0.05) in the dual therapy group. These results suggest that 8 might not exert an impact on the diversity of the intestinal flora, and increase the abundance of some probiotics while eradicating H. pylori. To figure out why compound 8 could exhibit anti-H. pylori activity without exerting an impact on intestinal microbiota, the antibacterial spectrum of 8 was evaluated. The antibacterial activities of 8 against common gram-positive and negative bacteria were shown in Supplementary Table S7. Compound 8 only exhibited a moderate antibacterial efficacy against Staphylococcus aureus ATCC 29213 (MIC value: 8g/mL), while being ineffective against all tested gram-negative bacteria. Therefore, the antibacterial spectrum indicates the specific inhibitory effect of compound 8 against H. pylori, while exerting minor impact on the intestinal microbiota. Proton pump inhibitor including OPZ is recommended to take before meals to avoid the over-production of gastric acid, so as to increase the stability of antibiotics. Considering that compound 8 possessed an ideal profile of acid stability, in vivo activity of compound 8 itself was evaluated, without the co-administration of OPZ. As shown in Fig. 2g, the mono-therapy of 8 showed a comparable potency compared with both the triple-therapy (OPZ+AMX+CLA) and the quadruple-therapy (OPZ+AMX+CLA+bismuth citrate). These results indicated that mono treatment of compound 8 may be applied as an alternative therapy of traditional triple or quadruple H. pylori eradication regimen. Bacterial cell morphologic changes can provide valuable clues on the antibacterial mode of action, and are often used for pilot mechanism investigation. Therefore, we performed scanning electron microscopy (SEM) and transmission electron microscopy (TEM) analysis on H. pylori ATCC 43504 after the treatment of compound 8. Bacterial cells were incubated with or without sub-MIC (1/2 MIC, 0.25g/mL) level of 8 for 2 days. The SEM and TEM analysis results showed that the integrity of the H. pylori outer membrane was compromised, and obvious perforations were observed compared to the untreated control group (Fig. 4a, b). This suggests that the mechanism of action of 8 might be related to its impact on the integrity of the bacterial outer membrane, which warrants further investigation. Mechanism of action and direct targets exploration on compound 8. a, b Images for morphology of H. pylori under electron microscope (a) SEM images of H. pylori treated without (upper) or with (lower) 8. b TEM images of H. pylori treated without (upper) or with (lower) 8. c The structure of the active photoaffinity probe 8-O. d Cy3-labeled target proteins were identified using fluorescent gel imaging. SecA (e) and BamD (f) were pulled down from H. pylori by using probe 😯 in immunoblot assay. SecA and BamD pulled down by 😯 were competitively inhibited by 8. The recombinant SecA (g) and BamD (h) proteins pulled down by 😯 were competitively inhibited by 8. Surface plasmon resonance (SPR) sensorgrams obtained on SecA (i)/BamD (j)-coated chips at different concentrations of 8. The thermal stability of SecA (k)/BamD (l) proteins with or without 8-treatment (n=3) The effectiveness of 8 against both drug-susceptible and resistant H. pylori strains suggests that it might possess a unique mechanism of action distinct from those of the first-line antibiotics used for the treatment of H. pylori infection. Hence, it is of great significance to identify the direct targets of 8 and further elucidate its specific mechanism of action. ABPP technique, a chemical biological tool for target protein exploration,43,44,45 was applied for the target fishing and identifying of 8 in this study, and the workflow of the specific process was described in Supplementary Fig. S5. Due to the lack of functional groups of 8 that can form covalent bonds with its target proteins, a photoaffinity probe (8-O, Fig. 4c) of 8 containing a diazirine photo cross-linking tag and an alkynyl functional group on position 3 was constructed. As mentioned above, mono substitution at position 3, and di-substitutions at positions 3 and 13 were beneficial for anti-H. pylori activity. Considering both structural similarity and synthetic feasibility, we opted for a probe design with a mono substitution at position 3. To make sure that probe 😯 possessed a similar mechanism as compound 8 and is suitable for target exploration, we assessed the effects of 😯 on the integrity of the H. pylori membrane through SEM and TEM analysis. As shown in the Supplementary Fig. S6, similar to compound 8, probe 😯 induced rupture and perforation of the H. pylori outer membrane. Subsequently, the probes activity against H. pylori was evaluated. As expected, 😯 exhibited comparable potency against the tested strains, with MICs ranging from 0.52g/mL as illustrated in Supplementary Table S8, indicating a similar mechanism with 8. Consequently, 😯 was deemed a viable functional probe for subsequent target exploration and verification. Following by the addition of probe 😯 (25M) to the lysate of H. pylori ATCC 43504, the mixture was incubated for 1h (Supplementary Fig. S5). Upon exposure to 365nm light, the diazirine photo cross-linking tag of 😯 could generate free radical fragments. These fragments could then form covalent bonds with adjacent hydroxyl groups of target proteins. Next, the alkyne reporter group of the 8-O/protein conjugate was coupled with an azide-modified fluorescent dye (Cy3) via a click reaction. The Cy3-labeled complex was separated using SDS-polyacrylamide gel electrophoresis (SDS-PAGE), with DMSO treatment serving as the blank control. Fluorescent bands with molecular weights (MW) ranging from 25150kDa were observed, and the addition of 8 competitively weakened several of these bands, as depicted in Fig. 4d. This result demonstrated that 😯 might partially occupy the binding sites of 8s targets, and was suitable for further verifications as a chemical tool. Similarly, a biotin-labeled complex was formed by coupling 8-O/protein conjugate with biotin-azide (Supplementary Fig. S5). After being purified and enriched, the complex was identified through liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis in three biological replicates. Totally, 24 proteins were identified twice in the analysis (Supplementary Table S9). Among these, two proteins belonging to the bacterial general secretory pathway (Sec pathway) and -barrel assembly machinery (BAM), namely protein translocase subunit SecA (SecA) and outer membrane protein assembly factor BamD (BamD), were selected for further verification, respectively. Since Sec pathway and BAM complex are known to be responsible for transporting and assembling the majority of OMPs to the outer membrane, targeting this system could potentially affect the integrity of the bacterial outer membrane, which is consistent with the findings in SEM and TEM analysis on 8-treated H. pylori cells. Thus, SecA and BamD were given priority for further investigation. Firstly, after pre-treatment of 8 in live H. pylori, SecA and BamD were successfully confirmed to be the potential direct targets of 8 through immunoblot assays using the 😯 probe in the pull-down experiments (Fig. 4e, f). Obvious competitive inhibition could be detected when 8 was pre-treated in situ, indicating possible specific interactions between 8 and these two proteins. Meanwhile, the recombinant H. pylori SecA and BamD proteins were also expressed and purified for further verification. In the presence of both UV (365nm) exposure and the active probe 😯 treatment, SecA/8-O conjugate with Cy3-labeling was successfully pulled down (Fig. 4g). Whereas, the fluorescent band was significantly weakened when either UV exposure or the active probe was absent, indicating the necessity of covalent bond formation between SecA and 😯 for successful pull-down. The fluorescence also faded when SecA was pre-treated with 8, indicating possible competitive inhibitions. Moreover, the fluorescent band of the 8-O/SecA complex almost vanished under the condition of 95C, suggesting that the active labeling of 😯 binding with SecA only occurred in the natively folded state rather than in the heat-treated unfolded state. Similar results were also observed in the BamD treatment group (Fig. 4h). It was found that 8 could dose-dependently bind to immobilized SecA and BamD with Kd values of 3.39 and 21.21M (Fig. 4i, j), respectively, in surface plasmon resonance (SPR) analysis. These results further confirmed the direct interactions between 8 and SecA or BamD. Besides, the cellular thermal shift assay (CESTA) was applied for further validation of their specific interactions, as displayed in Fig. 4k, l. Taking DMSO as the blank control, the thermal stability of the SecA protein decreased with a serial increase in temperatures ranging from 44 to 76C. However, with the addition of 8, the stability of SecA improved significantly, indicating the possible formation of an 8/SecA complex. The same trend was observed for BamD. These findings demonstrated that 8 might serve as a potential substrate of SecA as well as BamD and enhance the thermostability of these two proteins. To further figure out the specific binding sites and amino acid residues interacting with 8, protein mass spectrometry analysis was conducted. As shown in Fig. 5a, Escherichia coli (E. coli) strain Rosetta overexpressing H. pylori SecA or BamD was pretreated with or without 8 before probe 😯 was added. After proteome labeling and coupling with biotin, the specific peptide differences between the probe treatment and competitive inhibition group were analyzed through peptide fragment identification. Mass spectrometry analysis of the characteristic peaks was performed on the specific peptides of SecA/BamD, which might interact with 8. These characteristic peaks revealed that three different active cavities of SecA might serve as the potential binding sites of 8 (Supplementary Fig. S7). Then, the docking pattern analysis (Fig. 5b) was simulated in Discovery Studio 4.5 software (BIOVIA, San Diego, California, USA) for the prediction of the dominant contribution of each amino acid residue in these three cavities, and four potential residues forming hydrogen-bond interactions were selected for single-mutation verification. After being single mutated to alanine, the specific binding site was verified (KAENLFGVDNLYKIENAALSHHLDQALK), and 239-arginine inside this cavity was found to play a key role in SecA-8 interaction (the bright red ball, Fig. 5c). The two- and three-dimensional specific binding modes were displayed in Fig. 5c. Similarly, two adjacent peptide segments of BamD in space (one cavity), including YRPYVEYMQIKFILGQNELNRAIANVYK and IDETLEK, might contribute together to the interaction between BamD and 8 (Supplementary Fig. S8). Guided by the docking pattern and single mutation analysis, 171-glutamic acid and 209-serine were further confirmed to play the key roles among these residues. These findings provide solid evidences for the therapeutic targets verification of 8 and valuable insights for the exploration of novel candidates against H. pylori. The exploration of active binding sites between 8 and SecA/BamD. a Experimental workflow for binding site and interaction residues investigation and validation based on LC-MS/MS analysis. The predicted docking patterns between 8 and SecA (b)/BamD (d) were performed by Discovery Studio 4.5 software based on the peptide fragment difference identification results of LC-MS/MS analysis. Specific binding pattern between 8 and SecA (c)/BamD (e) Transcriptomic analysis was performed to gain comprehensive understanding of the antibacterial mechanism of 8 and verify its impact on OMPs (Fig. 6a, b). The inhibition of Sec pathway has been reported to impair the secretion of unfolded intracellular OMPs into the periplasmic space, leading to the over accumulation of OMPs within the intracellular space.46 As depicted in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, ribosome synthesis related genes were obviously down-regulated, which might due to the excessive accumulation of intracellular proteins. Specifically, after the treatment of 8, groEL and groES responsible for intracellular protein folding were significantly up-regulated, which might be used to deal with the excessive unfolded proteins (Fig. 6c). Lipopolysaccharide (LPS) transport highly dependents on Lpt machinery system, which consists of LptB located in cytoplasm and the other components in inner membrane (LptF, LptG), periplasmic space (LptA, LptC) or outer membrane (LptD, LptE). The impaired outer membrane transport will also result in a hampered LPS transport. It is worth noting that, as a cytoplasmic protein, the transcriptional level of LptB was significantly up-regulated after the treatment of 8 for the compensation of LPS deficiency in outer membrane. While as the Sec and Bam pathway was suppressed, the proteins located outside the inner membrane (LptA, LptD, LptE) could not be transported out and stacked in cytoplasm, which led to a negative regulation in the transcription of their coding genes (Fig. 6c). The transcription levels of H. pylori adhesion proteins in outer membrane, including BabA, SabA, and OipA were also significantly decreased in the transcriptome study (data not shown) and RT-qPCR validation (Fig. 6d). Collectively, these data suggest that the treatment of 8 arouses OMP aggregation in the cytoplasmic and periplasmic spaces and ineffective transportation, which is consistent with the Sec pathway and Bam machinery dysfunction. Compound 8 disturbs the OMPs related gene transcription and inhibits the protein function of SecA and BamD. a, b Transcriptome analysis of H. pylori with or without the treatment of 8 (n=3). a Volcano plot analysis (Red dots: 239 up-regulated genes; Green dots: 302 down-regulated genes), and (b) KEGG analysis. c The differential expression genes at transcriptional level related to the OMPs secretion and transport dysfunction. d RT-qPCR verifications on gene transcription of the key H. pylori OMPs after the treatment of 8 (n=3). e Inhibition of 8 on the ATPase activity of SecA (n=3). f The interaction of BamA and BamD was inhibited by 8 in Co-IP analysis. g The change of the total amount of H. pylori OMPs after the treatment of 8. h, i Confocal analysis on adhesive effect of 8-treated H. pylori to GES-1 cells. No treatment group (h); 8 treatment group (i). For cell nucleic acid staining: 4,6-diamidino2-phenylindole (DAPI); for cell membrane staining: 1,1-Dioctadecyl-3,3,3,3-tetramethylindodicarbocyanine, 4-chlorobenzenesulfonate salt (DiD); for bacteria staining: fluorescein isothiocyanate (FITC) SecA plays an indispensable role in the Sec complex as an ATPase.47 Therefore, the ATPase activity of SecA in the presence of 8 was measured. As depicted in Fig. 6e, 8 could dose dependently inhibit SecA, with an IC50 value of 11.53g/mL. Furthermore, to demonstrate the potential of SecA as an anti-H. pylori target, a previously reported SecA inhibitor CJ-21058 (IC50=7.0M)48 was evaluated for its anti-H. pylori potency. The MIC values of CJ-21058 against tested H. pylori strains were found to be in the range of 4-8g/mL (Supplementary Table S10), suggesting that SecA has the potential to be an attractive anti-H. pylori target and screening for SecA inhibitors could be an effective strategy for developing novel anti-H. pylori candidates. In gram-negative bacteria, the assembly of OMPs requires the Bam machinery complex, in which BamA is the central component. The -barrel domain of BamA interacts with four lipoproteins, including the essential lipoprotein BamD that directly interacts with BamA, and the other accessory lipoproteins BamB, BamC, and BamE.49 BamD facilitates the delivery of OMP substrates to BamA -barrel and the subsequent assembly. To investigate if the function of BamD was affected by 8, a Co-Immunoprecipitation (Co-IP) test was performed using GST-tagged BamD and His-tagged BamA. As depicted in Fig. 6f, BamD exhibited a strong interaction with BamA, and this effect was suppressed by 8, indicating that 8 might inhibit the function of the BAM machinery by affecting the BamA-BamD interaction.
