A reconstructed genome-scale metabolic model of Helicobacter pylori for predicting putative drug targets in clarithromycin and rifampicin resistance conditions
Sepideh Mofidifar
Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
Search for more papers by this authorCorresponding Author
Abbas Yadegar
Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Gastroenterology and Liver Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Correspondence
Abbas Yadegar, Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Shahid Arabi Ave., Yemen St., Velenjak, Tehran, Iran.
Email: a.yadegar@sbmu.ac.ir; babak_y1983@yahoo.com
Mohammad Hossein Karimi-Jafari, Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
Email: mhkarimijafari@ut.ac.ir
Search for more papers by this authorCorresponding Author
Mohammad Hossein Karimi-Jafari
Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
Correspondence
Abbas Yadegar, Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Shahid Arabi Ave., Yemen St., Velenjak, Tehran, Iran.
Email: a.yadegar@sbmu.ac.ir; babak_y1983@yahoo.com
Mohammad Hossein Karimi-Jafari, Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
Email: mhkarimijafari@ut.ac.ir
Search for more papers by this authorSepideh Mofidifar
Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
Search for more papers by this authorCorresponding Author
Abbas Yadegar
Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Gastroenterology and Liver Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Correspondence
Abbas Yadegar, Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Shahid Arabi Ave., Yemen St., Velenjak, Tehran, Iran.
Email: a.yadegar@sbmu.ac.ir; babak_y1983@yahoo.com
Mohammad Hossein Karimi-Jafari, Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
Email: mhkarimijafari@ut.ac.ir
Search for more papers by this authorCorresponding Author
Mohammad Hossein Karimi-Jafari
Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
Correspondence
Abbas Yadegar, Foodborne and Waterborne Diseases Research Center, Research Institute for Gastroenterology and Liver Diseases, Shahid Beheshti University of Medical Sciences, Shahid Arabi Ave., Yemen St., Velenjak, Tehran, Iran.
Email: a.yadegar@sbmu.ac.ir; babak_y1983@yahoo.com
Mohammad Hossein Karimi-Jafari, Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
Email: mhkarimijafari@ut.ac.ir
Search for more papers by this authorAbstract
Background
Helicobacter pylori is considered a true human pathogen for which rising drug resistance constitutes a drastic concern globally. The present study aimed to reconstruct a genome-scale metabolic model (GSMM) to decipher the metabolic capability of H. pylori strains in response to clarithromycin and rifampicin along with identification of novel drug targets.
Materials and Methods
The iIT341 model of H. pylori was updated based on genome annotation data, and biochemical knowledge from literature and databases. Context-specific models were generated by integrating the transcriptomic data of clarithromycin and rifampicin resistance into the model. Flux balance analysis was employed for identifying essential genes in each strain, which were further prioritized upon being nonhomologs to humans, virulence factor analysis, druggability, and broad-spectrum analysis. Additionally, metabolic differences between sensitive and resistant strains were also investigated based on flux variability analysis and pathway enrichment analysis of transcriptomic data.
Results
The reconstructed GSMM was named as HpM485 model. Pathway enrichment and flux variability analyses demonstrated reduced activity in the ribosomal pathway in both clarithromycin- and rifampicin-resistant strains. Also, a significant decrease was detected in the activity of metabolic pathways of clarithromycin-resistant strain. Moreover, 23 and 16 essential genes were exclusively detected in clarithromycin- and rifampicin-resistant strains, respectively. Based on prioritization analysis, cyclopropane fatty acid synthase and phosphoenolpyruvate synthase were identified as putative drug targets in clarithromycin- and rifampicin-resistant strains, respectively.
Conclusions
We present a robust and reliable metabolic model of H. pylori. This model can predict novel drug targets to combat drug resistance and explore the metabolic capability of H. pylori in various conditions.
CONFLICT OF INTEREST STATEMENT
The authors have no conflicts of interest to declare.
Open Research
DATA AVAILABILITY STATEMENT
Data sharing is not applicable to this article as no new data were created or analyzed in this study.
Supporting Information
Filename | Description |
---|---|
hel13074-sup-0001-FigureS1.docxWord 2007 document , 279.7 KB |
Figure S1. |
hel13074-sup-0002-FigureS2.docxWord 2007 document , 313.2 KB |
Figure S2. |
hel13074-sup-0003-TableS1.docxWord 2007 document , 29.4 KB |
Table S1. |
hel13074-sup-0004-TableS2.docxWord 2007 document , 37 KB |
Table S2. |
hel13074-sup-0005-FileS1.xlsxExcel 2007 spreadsheet , 28.9 KB |
Supplementary file 1. |
hel13074-sup-0006-FileS2.xlsxExcel 2007 spreadsheet , 66 KB |
Supplementary file 2. |
hel13074-sup-0007-FileS3.xlsxExcel 2007 spreadsheet , 63.9 KB |
Supplementary file 3. |
Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.
REFERENCES
- 1Malfertheiner P, Camargo MC, El-Omar E, et al. Helicobacter pylori infection. Nat Rev Dis Primers. 2023; 9(1): 19.
- 2Alavifard H, Nabavi-Rad A, Baghaei K, Sadeghi A, Yadegar A, Zali MR. Pyrosequencing analysis for rapid and accurate detection of clarithromycin resistance-associated mutations in Iranian Helicobacter pylori isolates. BMC Res Notes. 2023; 16(1): 136.
- 3Shu X, Ye D, Hu C, et al. Alarming antibiotics resistance of Helicobacter pylori from children in Southeast China over 6 years. Sci Rep. 2022; 12(1): 17754.
- 4Chauhan N, Tay ACY, Marshall BJ, Jain U. Helicobacter pylori VacA, a distinct toxin exerts diverse functionalities in numerous cells: an overview. Helicobacter. 2019; 24(1):e12544.
- 5Roszczenko-Jasińska P, Wojtyś MI, Jagusztyn-Krynicka EK. Helicobacter pylori treatment in the post-antibiotics era—searching for new drug targets. Appl Microbiol Biotechnol. 2020; 104: 9891-9905.
- 6Godavarthy PK, Puli C. From antibiotic resistance to antibiotic renaissance: a new era in Helicobacter pylori treatment. Cureus. 2023; 15(3):e36041.
- 7Pruitt KD, Tatusova T, Maglott DR. NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005; 33(suppl_1): D501-D504.
- 8Munoz-Ramirez ZY, Mendez-Tenorio A, Kato I, et al. Whole genome sequence and phylogenetic analysis show Helicobacter pylori strains from Latin America have followed a unique evolution pathway. Front Cell Infect Microbiol. 2017; 7: 50.
- 9Cheok YY, Lee CYQ, Cheong HC, et al. An overview of Helicobacter pylori survival tactics in the hostile human stomach environment. Microorganisms. 2021; 9(12): 2502.
- 10Thiele I, Palsson BØ. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010; 5(1): 93-121.
- 11Gu C, Kim GB, Kim WJ, Kim HU, Lee SY. Current status and applications of genome-scale metabolic models. Genome Biol. 2019; 20(1): 1-18.
- 12Kim WJ, Kim HU, Lee SY. Current state and applications of microbial genome-scale metabolic models. Curr Opin Syst Biol. 2017; 2: 10-18.
10.1016/j.coisb.2017.03.001 Google Scholar
- 13Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010; 28(3): 245-248.
- 14Paul A, Anand R, Karmakar SP, Rawat S, Bairagi N, Chatterjee S. Exploring gene knockout strategies to identify potential drug targets using genome-scale metabolic models. Sci Rep. 2021; 11(1): 213.
- 15Sertbas M, Ulgen KO. Genome-scale metabolic modeling for unraveling molecular mechanisms of high threat pathogens. Front Cell Dev Biol. 2020; 8:566702.
- 16Alonso-Vásquez T, Fondi M, Perrin E. Understanding antimicrobial resistance using genome-scale metabolic modeling. Antibiotics. 2023; 12(5): 896.
- 17Schilling CH, Covert MW, Famili I, Church GM, Edwards JS, Palsson BO. Genome-scale metabolic model of Helicobacter pylori 26695. J Bacteriol. 2002; 184(16): 4582-4593.
- 18Thiele I, Vo TD, Price ND, Palsson BØ. Expanded metabolic reconstruction of Helicobacter pylori (iIT341 GSM/GPR): an in silico genome-scale characterization of single-and double-deletion mutants. J Bacteriol. 2005; 187(16): 5818-5830.
- 19Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016; 44(D1): D457-D462.
- 20Caspi R, Billington R, Keseler IM, et al. The MetaCyc database of metabolic pathways and enzymes—a 2019 update. Nucleic Acids Res. 2020; 48(D1): D445-D453.
- 21King ZA, Lu J, Dräger A, et al. BiGG models: a platform for integrating, standardizing and sharing genome-scale models. Nucleic Acids Res. 2016; 44(D1): D515-D522.
- 22Wang H, Marcišauskas S, Sánchez BJ, et al. RAVEN 2.0: a versatile toolbox for metabolic network reconstruction and a case study on Streptomyces coelicolor. PLoS Comput Biol. 2018; 14(10):e1006541.
- 23Monk JM, Lloyd CJ, Brunk E, et al. i ML1515, a knowledgebase that computes Escherichia coli traits. Nat Biotechnol. 2017; 35(10): 904-908.
- 24Elbourne LD, Tetu SG, Hassan KA, Paulsen IT. TransportDB 2.0: a database for exploring membrane transporters in sequenced genomes from all domains of life. Nucleic Acids Res. 2017; 45(D1): D320-D324.
- 25Lee WC, Goh KL, Loke MF, Vadivelu J. Elucidation of the metabolic network of Helicobacter pylori J99 and Malaysian clinical strains by phenotype microarray. Helicobacter. 2017; 22(1):e12321.
- 26UniProt: the universal protein knowledgebase in 2023. Nucleic Acids Res. 2023; 51(D1): D523-D531.
- 27Alm RA, Ling L-SL, Moir DT, et al. Genomic-sequence comparison of two unrelated isolates of the human gastric pathogen Helicobacter pylori. Nature. 1999; 397(6715): 176-180.
- 28Lachance J-C, Lloyd CJ, Monk JM, et al. BOFdat: generating biomass objective functions for genome-scale metabolic models from experimental data. PLoS Comput Biol. 2019; 15(4):e1006971.
- 29Brooks MJ, Rajasimha HK, Roger JE, Swaroop A. Next-generation sequencing facilitates quantitative analysis of wild-type and Nrl−/− retinal transcriptomes. Mol Vis. 2011; 17: 3034.
- 30Ischenko D, Alexeev D, Shitikov E, et al. Large scale analysis of amino acid substitutions in bacterial proteomics. BMC Bioinformatics. 2016; 17(1): 1-10.
- 31Hirai Y, Haque M, Yoshida T, Yokota K, Yasuda T, Oguma K. Unique cholesteryl glucosides in Helicobacter pylori: composition and structural analysis. J Bacteriol. 1995; 177(18): 5327-5333.
- 32Testerman TL, Conn PB, Mobley HL, McGee DJ. Nutritional requirements and antibiotic resistance patterns of Helicobacter species in chemically defined media. J Clin Microbiol. 2006; 44(5): 1650-1658.
- 33Salama NR, Shepherd B, Falkow S. Global transposon mutagenesis and essential gene analysis of Helicobacter pylori. J Bacteriol. 2004; 186(23): 7926-7935.
- 34Chalker AF, Minehart HW, Hughes NJ, et al. Systematic identification of selective essential genes in Helicobacter pylori by genome prioritization and allelic replacement mutagenesis. J Bacteriol. 2001; 183(4): 1259-1268.
- 35Ebrahim A, Lerman JA, Palsson BO, Hyduke DR. COBRApy: constraints-based reconstruction and analysis for python. BMC Syst Biol. 2013; 7: 1-6.
- 36Jensen PA, Papin JA. Functional integration of a metabolic network model and expression data without arbitrary thresholding. Bioinformatics. 2011; 27(4): 541-547.
- 37Jensen PA, Lutz KA, Papin JA. TIGER: toolbox for integrating genome-scale metabolic models, expression data, and transcriptional regulatory networks. BMC Syst Biol. 2011; 5: 1-12.
- 38Gudmundsson S, Thiele I. Computationally efficient flux variability analysis. BMC Bioinformatics. 2010; 11(1): 1-3.
- 39Wu T, Hu E, Xu S, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. The Innovations. 2021; 2(3):100141.
- 40Liu B, Zheng D, Zhou S, Chen L, Yang J. VFDB 2022: a general classification scheme for bacterial virulence factors. Nucleic Acids Res. 2022; 50(D1): D912-D917.
- 41Wishart DS, Feunang YD, Guo AC, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018; 46(D1): D1074-D1082.
- 42Shende G, Haldankar H, Barai RS, Bharmal MH, Shetty V, Idicula-Thomas S. PBIT: pipeline builder for identification of drug targets for infectious diseases. Bioinformatics. 2017; 33(6): 929-931.
- 43Liechti G, Goldberg JB. Helicobacter pylori relies primarily on the purine salvage pathway for purine nucleotide biosynthesis. J Bacteriol. 2012; 194(4): 839-854.
- 44Iwatani S, Nagashima H, Reddy R, Shiota S, Graham DY, Yamaoka Y. Identification of the genes that contribute to lactate utilization in Helicobacter pylori. PLoS One. 2014; 9(7):e103506.
- 45Caspi R, Billington R, Fulcher CA, et al. BioCyc: a genomic and metabolic web portal with multiple omics analytical tools. FASEB J. 2019; 33(S1):4732.
10.1096/fasebj.2019.33.1_supplement.473.2 Google Scholar
- 46Barrett T, Wilhite SE, Ledoux P, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2012; 41(D1): D991-D995.
- 47Geng X, Li W, Chen Z, et al. The bifunctional enzyme SpoT is involved in the clarithromycin tolerance of Helicobacter pylori by upregulating the transporters HP0939, HP1017, HP0497, and HP0471. Antimicrob Agents Chemother. 2017; 61(5):10-1128. doi:10.1128/aac.02011-02016
- 48Mosaei H, Zenkin N. Inhibition of RNA polymerase by rifampicin and rifamycin-like molecules. EcoSal Plus. 2020; 9(1). doi:10.1128/ecosalplus.ESP-0017-2019
- 49Dascălu RI, Bolocan A, Păduaru DN, et al. Multidrug resistance in Helicobacter pylori infection. Front Microbiol. 2023; 14:1128497.
- 50Naorem RS, Pangabam BD, Bora SS, et al. Identification of putative vaccine and drug targets against the methicillin-resistant Staphylococcus aureus by reverse vaccinology and subtractive genomics approaches. Molecules. 2022; 27(7): 2083.
- 51Beceiro A, Tomás M, Bou G. Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? Clin Microbiol Rev. 2013; 26(2): 185-230.
- 52Geisinger E, Isberg RR. Interplay between antibiotic resistance and virulence during disease promoted by multidrug-resistant bacteria. J Infect Dis. 2017; 215(suppl_1): S9-S17.
- 53Hussain A, Mazumder R, Asadulghani M, Clark TG, Mondal D. Combination of virulence and antibiotic resistance: a successful bacterial strategy to survive under hostile environments. In: A Kumar, S Tenguria, eds. Bacterial Survival in the Hostile Environment. Academic Press; 2023: 101-117.
10.1016/B978-0-323-91806-0.00004-7 Google Scholar
- 54Agoni C, Olotu FA, Ramharack P, Soliman ME. Druggability and drug-likeness concepts in drug design: are biomodelling and predictive tools having their say? J Mol Model. 2020; 26: 1-11.
- 55Nguyen M, Olson R, Shukla M, VanOeffelen M, Davis JJ. Predicting antimicrobial resistance using conserved genes. PLoS Comput Biol. 2020; 16(10):e1008319.
- 56Zhu Y, Zhao J, Li J. Genome-scale metabolic modeling in antimicrobial pharmacology. Eng Microbiol. 2022; 2(2):100021.
- 57Sholeh M, Khoshnood S, Azimi T, et al. The prevalence of clarithromycin-resistant Helicobacter pylori isolates: a systematic review and meta-analysis. PeerJ. 2023; 11:e15121.
- 58Schubert JP, Warner MS, Rayner CK, et al. Increasing Helicobacter pylori clarithromycin resistance in Australia over 20 years. Intern Med J. 2022; 52(9): 1554-1560.
- 59Boyanova L, Hadzhiyski P, Gergova R, Markovska R. Evolution of Helicobacter pylori resistance to antibiotics: a topic of increasing concern. Antibiotics. 2023; 12(2): 332.
- 60Regnath T, Raecke O, Enninger A, Ignatius R. Increasing metronidazole and rifampicin resistance of Helicobacter pylori isolates obtained from children and adolescents between 2002 and 2015 in southwest Germany. Helicobacter. 2017; 22(1):e12327.
- 61Marques AT, Vítor JM, Santos A, Oleastro M, Vale FF. Trends in Helicobacter pylori resistance to clarithromycin: from phenotypic to genomic approaches. Microbial Genomics. 2020; 6(3):e000344.
- 62Hussein RA, Al-Ouqaili MT, Majeed YH. Detection of clarithromycin resistance and 23SrRNA point mutations in clinical isolates of Helicobacter pylori isolates: phenotypic and molecular methods. Saudi J Biol Sci. 2022; 29(1): 513-520.
- 63Hamouche L, Poljak L, Carpousis AJ. Ribosomal RNA degradation induced by the bacterial RNA polymerase inhibitor rifampicin. RNA (New York, NY). 2021; 27(8): 946-958.
- 64Bhargava P, Collins JJ. Boosting bacterial metabolism to combat antibiotic resistance. Cell Metab. 2015; 21(2): 154-155.
- 65Lopatkin AJ, Bening SC, Manson AL, et al. Clinically relevant mutations in core metabolic genes confer antibiotic resistance. Science. 2021; 371(6531):eaba0862.
- 66Schrader SM, Botella H, Jansen R, et al. Multiform antimicrobial resistance from a metabolic mutation. Sci Adv. 2021; 7(35):eabh2037.
- 67Wareth G, Neubauer H, Sprague LD. A silent network's resounding success: how mutations of core metabolic genes confer antibiotic resistance. Signal Transduct Target Ther. 2021; 6(1): 301.
- 68Zafar M, Jahan H, Shafeeq S, et al. Clarithromycin exerts an antibiofilm effect against Salmonella enterica Serovar Typhimurium rdar biofilm formation and transforms the physiology towards an apparent oxygen-depleted energy and carbon metabolism. Infect Immun. 2020; 88(11):e00510-20. doi:10.1128/iai.00510-00520
- 69Raj DS, Kumar Kesavan D, Muthusamy N, Umamaheswari S. Efflux pumps potential drug targets to circumvent drug resistance—multi drug efflux pumps of Helicobacter pylori. Mater Today Proc. 2021; 45: 2976-2981.
- 70Kontoghiorghes GJ, Kontoghiorghe CN. Iron and chelation in biochemistry and medicine: new approaches to controlling iron metabolism and treating related diseases. Cells. 2020; 9(6): 1456.
- 71Mobarra N, Shanaki M, Ehteram H, et al. A review on iron chelators in treatment of iron overload syndromes. Int J Hematol Oncol Stem Cell Res. 2016; 10(4): 239-247.
- 72Chhabra R, Saha A, Chamani A, Schneider N, Shah R, Nanjundan M. Iron pathways and iron chelation approaches in viral, microbial, and fungal infections. Pharmaceuticals. 2020; 13(10): 275.
- 73Ge X, Cai Y, Chen Z, et al. Bifunctional enzyme SpoT is involved in biofilm formation of Helicobacter pylori with multidrug resistance by upregulating efflux pump Hp1174 (gluP). Antimicrob Agents Chemother. 2018; 62(11):e00957-18. doi:10.1128/aac.00957-00918
- 74Jiang X, Duan Y, Zhou B, et al. The cyclopropane fatty acid synthase mediates antibiotic resistance and gastric colonization of Helicobacter pylori. J Bacteriol. 2019; 201(20):e00374-19. doi:10.1128/jb.00374-00319
- 75Gao G, Liu X, Pavlovsky A, Viola RE. Identification of selective enzyme inhibitors by fragment library screening. J Biomol Screen. 2010; 15(9): 1042-1050.
- 76Fong P, Hao CH, Io CC, Sin PI, Meng LR. In silico and in vitro anti-Helicobacter pylori effects of combinations of phytochemicals and antibiotics. Molecules. 2019; 24(19):3608.