Next Article in Journal
Exploring Immersion Coating as a Cost-Effective Method for Small-Scale Production of Enteric-Coated Gelatin Capsules
Next Article in Special Issue
In Vitro Antibiofilm Activity of Fosfomycin Alone and in Combination with Other Antibiotics against Multidrug-Resistant and Extensively Drug-Resistant Pseudomonas aeruginosa
Previous Article in Journal
Static Magnetic Field Reduces the Anticancer Effect of Hinokitiol on Melanoma Malignant Cells—Gene Expression and Redox Homeostasis Studies
Previous Article in Special Issue
Antimicrobial Activity of Ceragenins against Vancomycin-Susceptible and -Resistant Enterococcus spp.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications

by
Md. Kaderi Kibria
1,2,
Md. Ahad Ali
1,3,
Muhammad Yaseen
4,
Imran Ahmad Khan
5,
Mashooq Ahmad Bhat
6,
Md. Ariful Islam
1,
Rashidul Alam Mahumud
7 and
Md. Nurul Haque Mollah
1,*
1
Bioinformatics Laboratory, Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh
2
Department of Statistics, Hajee Mohammad Danesh Science and Technology University, Dinajpur 5200, Bangladesh
3
Department of Chemistry, University of Rajshahi, Rajshahi 6205, Bangladesh
4
Institute of Chemical Sciences, University of Swat, Main Campus, Charbagh 19130, Pakistan
5
Department of Chemistry, Government College University, Faisalabad 38000, Pakistan
6
Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh 11421, Saudi Arabia
7
NHMRC Clinical Trials Centre, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2024, 17(4), 432; https://doi.org/10.3390/ph17040432
Submission received: 26 February 2024 / Revised: 18 March 2024 / Accepted: 21 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Multidrug Resistance in Bacteria and New Therapeutic Options)

Abstract

:
SARS-CoV-2 infections, commonly referred to as COVID-19, remain a critical risk to both human life and global economies. Particularly, COVID-19 patients with weak immunity may suffer from different complications due to the bacterial co-infections/super-infections/secondary infections. Therefore, different variants of alternative antibacterial therapeutic agents are required to inhibit those infection-causing drug-resistant pathogenic bacteria. This study attempted to explore these bacterial pathogens and their inhibitors by using integrated statistical and bioinformatics approaches. By analyzing bacterial 16S rRNA sequence profiles, at first, we detected five bacterial genera and taxa (Bacteroides, Parabacteroides, Prevotella Clostridium, Atopobium, and Peptostreptococcus) based on differentially abundant bacteria between SARS-CoV-2 infection and control samples that are significantly enriched in 23 metabolic pathways. A total of 183 bacterial genes were found in the enriched pathways. Then, the top-ranked 10 bacterial genes (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) were selected as the pathogenic bacterial key genes (bKGs) by their protein–protein interaction (PPI) network analysis. Then, we detected bKG-guided top-ranked eight drug molecules (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Entreatinib, Digitoxin, and Elbasvir) by molecular docking. Finally, the binding stability of the top-ranked three drug molecules (Bemcentinib, Ledipasvir, and Velpatasvir) against three receptors (hldD, mlaA, and lptD) was investigated by computing their binding free energies with molecular dynamic (MD) simulation-based MM-PBSA techniques, respectively, and was found to be stable. Therefore, the findings of this study could be useful resources for developing a proper treatment plan against bacterial co-/super-/secondary-infection in SARS-CoV-2 infections.

Graphical Abstract

1. Introduction

The novel coronavirus disease 2019 (COVID-19), one of the most severe respiratory infections caused by the Severe Acute Respiratory Syndrome CoronaVirus-2 virus (SARS-CoV-2), remains a risk factor for human life. As of 7 January 2024, there have been over 774 million cases and over seven million deaths worldwide [1]. Though different vaccination programs have been reducing the rate of infection, a significant number of people are still being infected around the globe due to the random mutation in the RNA sequence of SARS-CoV-2 [2], which demands effective drugs as supplementary treatment for COVID-19. The spike (S) protein of SARS-CoV-2 interacts with the host ACE2 (angiotensin-converting enzyme 2) protein to initiate the infection [3,4]. Inhibition of this protein is essential to neutralize SARS-CoV-2 before the start of a cytokine storm, for which the immune system generates a number of inflammatory signals that can cause organ failure and patient death. The major protease (MPro/3ClPro), the papain-like protease (PLpro) [5,6] and the RNA-dependent RNA polymerase (RdRp) [7,8] of SARS-CoV-2 are associated with the infections by catalyzing the replication of RNA from an RNA template. Therefore, these proteins were used as targets/receptors to explore candidate drugs against SARS-CoV-2 infections [6,9]. However, COVID-19 patients with weak immunity may suffer from different complications due to microbial co-infections/super-infections/secondary-infections [10,11,12,13,14,15,16]. It should be noted here that co-infections refer to the simultaneous infection with two or more different microbial pathogens [14,17]. Super-infection occurs when already infected patients are re-infected with other microbial pathogens [15,16]. Secondary infections occur when already infected patients are re-infected, depending on the outcomes of the first infection [14,17,18]. COVID-19 patients often suffer from different symptoms, including fever [19], cough [20], sore throat [21], congestion or runny nose [22], difficulty breathing [23], loss of smell or taste [24], fatigue [21], muscle pain [25], headache [26], nausea or vomiting [27], diarrhea [28], meningitis [29], pneumonia [29] that are associated with different microbial pathogens [30,31,32,33]. Therefore, microbial pathogens may be considered the critical risk factors for different complications and the increasing mortality rates of COVID-19 patients [14,15,16,17,34]. In order to inhibit microbial pathogens, it is required to identify pathogenic genes that may lead to exploring potential inhibitors. It may be mentioned here that either disease-causing microbial proteins or host-proteins are considered the receptor proteins for exploring candidate drug molecules in order to inhabit microbial infections [6,9,35,36,37,38,39,40], since pathogenic proteins interact with host-proteins to develop infections [39,41,42].
There are some studies based on 16S rRNA-Seq profiles that have explored the association of microbial compositions with the complications of SARS-CoV-2 infections [43,44,45,46,47,48,49,50,51,52,53,54,55]. For instance, the bacterial composition of the respiratory tract and gut is associated with the occurrence and severity of disease in numerous respiratory viral infections (RVs) affecting the eventual respiratory health of COVID-19 [43]. The bacterial phyla groups Firmicutes and Bacteroidetes, as well as commensals of the phyla, are directly associated with the pathogenesis and severity of SARS-CoV-2 infections due to their regulatory function with the ACE2 host gene [53]. The Bacteroidetes phylum is associated with the downregulation of the ACE2 gene, while the Firmicutes phylum is associated with the upregulation, resulting in lower and higher infection rates, respectively [54]. A study has analyzed the oral and gut microbiome data and found that Rothia mucilaginosa and Granulicatella bacterial taxa are related to SARS-CoV-2 infections [55]. Thus, we observed that previous studies based on 16S rRNA-Seq profiles have explored only the bacterial taxa that are associated with co-/super-/secondary infections. They did not explore such infection-causing bacterial genes and their inhibitors. Two studies recommended antimicrobial agents in order to inhibit microbial pathogens. However, some microbial pathogens may progressively develop resistance to antimicrobial drugs due to environmental changes. This type of drug-resistant microbial pathogens is gradually increasing due to the availability of a small number of antimicrobial agents [56,57]. In this case, multi-targeted different variants of alternative antimicrobial drugs may be essential for effective treatment against co-/super-/secondary infections. Therefore, this study attempted to explore such infection-causing bacterial key genes (bKGs) based on 16S rRNA-Seq profiles and associated drug molecules in order to inhibit those infection-causing pathogenic bacteria, since 16S rRNA-Seq profile analysis performs better in order to identify novel pathogens and non-cultured bacteria [58,59].

2. Results

2.1. Diversity Analysis

All alpha diversity indices (richness measure: observed and Chao1, and evenness measure: Shannon and Simpson) of the oral and gut microbiota showed a significant difference (p < 0.01) between the COVID-19 patients and healthy group (Figure 1). The median observed index was 2941.5 (IQR: 2502, 3874) for the healthy group and 937.5 (IQR: 717.25, 1199.75) for the COVID-19 patients, indicating a significant decrease in observed species. Similarly, the median of Chao1 index was 4852.912 (IQR: 4250.33, 6318.36) for healthy patients and 1729.328 (IQR: 1303.84, 2200.93) for COVID-19 patients. The Shannon index was 4.875 (IQR: 4.455, 5.263) and 4.145 (IQR: 3.605, 4.556), and the Simpson index was 0.969 (IQR: 0.954, 0.977) for healthy and 0.953 (IQR: 0.919, 0.971) for COVID-19 patients, respectively (see Table S1 in Supplementary Materials). The findings indicate a significant decrease in four alpha diversity species among COVID-19 patients compared to healthy individuals, with a substantial effect size observed in Cliff’s Delta statistics for observed, Chao1, and Shannon indices, highlighting a significant distinction between the COVID-19 and healthy groups (see Table S1 in Supplementary Materials).
In addition to alpha diversity indices, beta diversity was employed to evaluate compositional changes in the microbial communities. Principle Coordinate Analysis (PCoA) based on Bray–Curtis distance revealed distinct clustering patterns among samples with varied distributions (see Figure 2a). Permutational Multivariate Analysis of Variance (PERMANOVA) on Bray–Curtis distance (F = 22.575, p-value: 0.001) indicated a significant difference between COVID-19 patients and healthy individuals.
A total of 2802 OTUs were found to be shared between COVID-19 patients and healthy controls, accounting for 37.71% and 32.60% of the total OTUs in healthy controls and COVID-19 patients, respectively (see Figure 2b). These results suggest that compositional differences were associated with the presence or absence of specific taxa and their relative abundances in the samples.

2.2. Taxonomy Analysis and Identification of Differentially Abundant Bacterial Compositions (DABCs)

The oral and gut microbiota profiles at the phylum and genus levels consisted of approximately 22 phyla and 210 genera, respectively. Among them, the top 10 most abundant phyla and the top 12 most abundant genera are shown in Figure 3. The results indicate that the abundance of the Bacteroidetes phylum was relatively higher in COVID-19 patients than in healthy subjects (26.9% vs. 19.6%), followed by Actinobacteria (13.941% vs. 7.703%), Fusobacteria (10.66% vs. 6.92%), TM7 (1.037% vs. 0.011%), and SR1 (0.382% vs. 0.122%). On the other hand, the relative abundances of Firmicutes, Proteobacteria, Cyanobacteria, and Verrucomicrobia phyla decreased in COVID-19 patients compared to healthy subjects. However, Thermi exhibited less significant changes in patients and healthy individuals at the phylum level (Figure 3a and Table S2 in Supplementary Materials). Also, among the 12 most abundant genera in the oral and gut microbiota in COVID-19 patients, Prevotella (37.21% vs. 27.85%), Rothia (5.49% vs. 1.84%), Actinomyces (5.17% vs. 2.95%), Clostridium (3.15% vs. 1.59%), Streptobacillus (3.06% vs. 0.68%), and Campylobacter (3.70% vs. 0.51%) showed increased relative abindances in COVID-19 patients compared to healthy subjects. Conversely, Bacteroides (25.73% vs. 34.08%), Corynebacterium (4.89% vs. 9.62%), Fusobacterium (4.02% vs. 8.82%), and Faecalibacterium (2.48% vs. 7.05%) were comparatively less abundant in COVID-19 patients than in healthy individuals (Figure 3b and Table S2 in Supplementary Materials). Alternatively, insignificant changes were observed in the cases of Porphyromonas and Bifidobacterium.
Based on their relative abundances, it is possible to identify which phyla and genera were more abundant, but it remains unclear which taxa were significantly differentially abundant in the SARS-CoV-2 infected group. LEfSe analysis was used to identify significantly abundant phyla and genera between groups (LDA score > 3.5, p < 0.05). The results reveal that five phyla were significantly different between the groups. Among them, Bacteroidetes, Actinobacteria, and Fusobacteria phyla were significantly enriched in COVID-19 patients, while Proteobacteria and Firmicutes were significantly decreased (see Figure 4). On the other hand, 25 genera showed different abundance between the groups, whereas the Bacteroides, Campylobacter, Actinomyces, and Rothia genera were highly enriched in COVID-19 patients, while Streptococcus, Neisseria Veillonella, Blautia, Megamonas Escherichia, and Megasphaera were less enriched compared to healthy patients. The LDA scores revealed a significant abundance of Bacteroides and Streptococcus as microbiological biomarkers in the treatment group and healthy group (see Figure 4). These results demonstrate a notable difference in gut microbiota composition between healthy and COVID-19 patients.
Further, the zero-inflated Gaussian (ZIG) model was employed using the metagenomeSeq R package to identify significantly differential bacterial genera based on the adjusted p-value and log2FC values. Almost 111 bacterial genera showed significant differences between the groups at p < 0.05, which were used for further analysis. All the genera identified as differential in the LEfSe analysis were included among the 111 differentially abundant genera identified in the ZIG model analysis. We represented the top 20 highly differential genera based on the adjusted p-value and log2FC in Table 1, noting that nine bacterial genera were upregulated while the remaining 11 were downregulated in COVID-19 patients. Particularly, genera such as, Neisseria, Streptococcus, Leptotrichia, and Veillonella were significantly decreased in COVID-19 patients, whereas Bacteroides, Parabacteroides, Prevotella, and Clostridium were significantly increased compared to healthy individuals. Consequently, in COVID-19 infections, the abundances of genera such as Bacteroides, Parabacteroides, Prevotella Clostridium, Atopobium, and Peptostreptococcus showed significant increases.

2.3. Identification of Pathway-Based Bacterial Key Genes (bKGs) from DABCs

To comprehend the functional role of COVID-19-related bacteria, PICRUSt2 analysis was conducted to predict gut microbiota metagenomes from the 111 differentially abundant genera identified in the ZIG model. Furthermore, the results of this functional analysis were compared to the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs to assess differences in predicted functional gene abundances. As a result, 23 MetaCyc signaling pathways differed significantly in the mean proportions between the COVID-19 and healthy groups (see Figure 5a). Of these, 14 MetaCyc pathways showed significant increases, including incomplete reductive TCA and pyruvate fermentation to acetate and lactate II, while 9 pathways, including palmitate biosynthesis II, palmitoleate biosynthesis I, and oleate biosynthesis IV, were significantly decreased in the COVID-19 patients group. Furthermore, we identified 183 KEGG genes corresponding to 23 differentially abundant pathways and hence explored the 10 most functional bKGs (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) through protein–protein interaction network analysis. These identified bKGs may represent possible functional gene contents altered due to secondary infection with SARS-CoV-2, and the findings were used for further molecular docking analyses (see Figure 5b).

2.4. bKGs-Guided Drug Repurposing by Molecular Docking

To select the top-ranked potential candidate drug molecules for COVID-19 treatment from the pools of 786 published molecules using molecular docking, we considered our proposed top-ranked 10 bKG-mediated proteins as receptors. Subsequently, we downloaded the 3D structures of our key proteins accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB from the Protein Data Bank (PDB), corresponding to the source codes 4hr7, 6h9o, 7qcf, 4ej0, 4isa, 4rhb, 5nuq, 7akc, 1fiy, and 5vtg, respectively. Heteroatoms, water molecules, and unused ligands were then removed using BIOVIA Discovery Studio and PyMOL. Subsequently, protein energy minimization was conducted using the GROMOS 43B1 force field. Finally, each protein was converted into an acceptable PDBQT format using AutoDock Tools. Molecular docking analyses were performed between the proposed 10 receptors and the 786 meta-drug agents to calculate the binding affinity scores (kcal/mol) for each receptor–ligand pair. These binding affinity scores were arranged in a matrix A = (Aij) in descending order, with receptor proteins in rows and ligands in columns (see Figure 6a). Consequently, we selected the top ranked 10 drugs (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Theaflavin digallate, Telcagepant, Entreatinib, Digitoxin, and Elbasvir) as candidate drug molecules, with an average binding affinity score of −9.618 kcal/mol against the proposed 10 receptors. To evaluate the binding performance of our proposed 10 candidate drug molecules against the top-ranked independent receptors, we reviewed 98 published articles identifying hub proteases/proteins with SARS-CoV-2 infections. This review yielded a total of 215 hub proteins (see Table S4 in Supplementary Materials).
Out of 215, only 9 hub proteins (ACE2, RdRp, 3CLpro, S, TMPRSS2, PLpro, IL6, TNF, and N) were found in at least six articles. These hub proteins were considered independent receptors to investigate their binding capacity with our proposed drug molecules. We downloaded those independent receptor proteins from the Protein Data Bank (PDB) with source codes 2ajf, 7bv2, 6lu7, 6vsb, 7meq, 6w9c, 1alu, 1a8m, and 6m3m, respectively. Then, molecular docking analysis was performed with these nine independent hub proteins against the same set of published 786 ligands. Then, their binding affinities were arranged in descending order and visualized in a matrix plot (see Figure 6b).
The top-ranked 10 drug molecules produced a binding affinity score lower than −7.0 kcal/mol against the top-ranked 9 independent receptors. Among the top 10 ligands (according to the binding affinity), 8 of them (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Telcagepant, Digitoxin, and Elbasvir) were common with our proposed 10 lead compounds (see Figure 6). Therefore, we considered these 8 drug molecules as the proposed candidates. Interestingly, three ligands (Bemcentinib, Ledipasvir, and Velpatasvir) were ranked highest against both proposed and independent receptors, suggesting their potential effectiveness against COVID-19. To assess the binding performance of our proposed drug molecules compared to the FDA approved two COVID-19 drugs (Molnupiravir and Nirmatrelvir) [60,61], we analyzed their overall docking scores against both proposed and independent receptors (see Table S5 in Supplementary Materials). Interestingly, neither Molnupiravir nor Nirmatrelvir ranked among the top 30 drug molecules based on their average binding affinity scores in both scenarios (see Figure 6). The interacting properties of these five target-ligand complexes are displayed in Figure 7 and also given in Table S6. The interacting complex (3D) surface view, pose view, and protein–ligand interactions are shown in Figure 7. The hldD-Bemcentinib complex was formed with one Pi-donor hydrogen bond at ARG213, one Pi-sigma bond at PHE200, two Pi-Pi T-shaped bonds at HIS186 and TYR292, an alkyl bond at VAL168, and a Pi-Alkyl bond at LYS198. Similarly, the other 4 target-ligand complexes were formed as given by their data in Table S6.
Additionally, the effectiveness and indemnity level of the top three lead compounds were assessed by evaluating their toxicological and pharmacokinetic properties (see Table S3 in Supplementary Materials). Central nervous system (CNS) permeability is less than −2, which indicates an inability to permeate the blood–brain barrier. Moreover, the number of H-bond acceptors and donors for the three compounds was less than ten and five, respectively, suggesting that the compounds are plausible drug candidates. Finally, no toxicity or carcinogenic profiles were observed for our three lead compounds (see Table S3 in Supplementary Materials).
A molecular dynamics (MD) simulation was conducted to verify the structural stability of the top three target-ligand complexes obtained from docking. The root mean square deviations (RMSD) of C-alpha atoms for the docked complexes are shown in Figure 8a. At the beginning of the simulation, all three complexes, hldD-Bemcentinib, mlaA-Ledipasvir, and lptD-Velpatasvir showed the same trend. However, the hldD-Bemcentinib complex exhibited a higher flexibility at 20–40 ns, while the mlaA-Ledipasivir complex showed a lower trend after 20 ns until 30 ns. Eventually, all three complexes reached a steady state after 60 ns throughout the simulation period. The RMSD value of the three complexes was less than 2.5 Å, indicating their stable and rigid structure. Furthermore, the stable and steady nature of the three complexes was assessed using the radius of gyration (Rg) and the solvent accessible surface area (SASA) (see Figure 8c,d). The lptD-velpatasvir complex exhibited a higher Rg, suggesting greater flexibility, while the mlaA-ledipasvir complex showed a lower Rg, indicating a more stable profile. On the other hand, the SASA analysis revealed changes in protein volume over time (see Figure 8d), with a lower SASA value indicating a condensed complex and a higher SASA value indicating an extended surface area.
The SASA values for the hldD-bemcentinib complex peaked at 25–50 ns, suggesting an expanded surface area, whereas the mlaA-ledipasvir complex exhibited lower SASA values at 15–45 ns, indicating a condensed surface. Conversely, the lptD-velpatasvir complex maintained a consistent profile throughout the simulation period. By 70 ns, all three complexes demonstrated stable and steady behavior. The binding energy of the three complexes was also calculated using the MM-PBSA method, with the more positive energy of the complexes indicating better binding (see Figure 8b). The mean binding energies of the hldD-bemcentinib, mlaA-ledipasvir, and lptD-velpatasvir complexes were 79.395, 148.595, and 12.078 kJ/mol, respectively. The mlaA-ledipasvir complex exhibited higher binding energy than the other two, suggesting better binding of the ligand molecules. The other two complexes had similar positive free energies, indicating better binding with these ligand molecules.

3. Discussion

Although SARS-CoV-2 is the main cause of COVID-19, there are some microbial compositions that are associated with the complications of COVID-19 by their co-infections, super-infections and secondary infections [14,15,16,17,34]. This study aims to explore co-/super-/secondary infections causing bacterial key genes (bKGs) and associated drug molecules through bioinformatics analysis. Alpha and beta diversity analyses showed that there are some bacterial compositions that are able to differentiate COVID-19 patients from healthy individuals (see Figure 3 and Figure 4). Bacterial compositions were observed at the phylum and genus level, revealing that at the phylum level, Bacteroidetes had a higher abundance in COVID-19 patients compared to healthy individuals. On the other hand, Firmicutes, Proteobacteria, and Cyanobacteria showed lower abundances in the gut microbiota of COVID patients than healthy patients. The ratio of Firmicutes to Bacteroidetes is known to influence health status [62], and an abundance of Bacteroidetes is also associated with acute diarrheal disease [63]. Also, the most abundant genera in COVID-19 patients compared to healthy individuals were Prevotella, othia, Bacteroides, and Campylobacter indicating higher abundances of these genera in infected individuals. While Prevotella has been associated with HIV in children and chronic obstructive pulmonary disease [64,65], Prevotella melaninogenica pathogen infection is also associated with meningitis diseases [66]. Bacteroidetes is also associated with gastric carcinogenesis [67,68] as well as colorectal cancer [69], whereas Bacteroides Acidifaciens, another differentially abundant bacterium, is implicated in liver damage [33]. Campylobacter enrichment in the gut microbiome of the patients is associated with diarrhea [28], gastrointestinal issues [70], and inflammatory bowel disease [71], while Campylobacter fetus pathogens can cause liver disease and diabetes mellitus [32]. Another upregulated genera identified in this study is Clostridium difficile, which is associated with diarrheal disease, another symptom of COVID-19 [72]. Alternatively, the genera Neisseria, Streptococcus, Fusobacterium, Faecalibacterium, Corynebacterium, and Staphylococcus were comparatively less abundant in COVID patients, suggesting that their reduced abundance may indicate COVID-19 infection. Among them, Neisseria contains two pathogens, Neisseria gonorrhoeae and Neisseria meningitides, which are associated with the meningococcal disease [73]. In another study, Streptococcus and Fusobacterium are often highly abundant in the oral cavity [74], and Fusobacterium is also highly associated with periodontal disease [75]. Another downregulated genus, Haemophilus influenzae, is implicated in meningitis and pneumonia diseases, which are common co-infections in COVID-19 patients [76]. On the other hand, the abundance of Faecalibacterium in the gut microbiome produces short-chain fatty acids, which have been associated with reduced intestinal inflammation [77]. Differentially abundant analysis revealed highly differential genera Neisseria, Streptococcus, Leptotrichia, and Bacteroides among the groups, belonging to the phyla Proteobacteria, Firmicutes, Fusobacteria, and Bacteroidetes, respectively.
We identified 23 significantly enriched metabolic pathways, including the incomplete reductive TCA and pyruvate fermentation to acetate and lactate II pathways, which were highly increased in COVID-19 patients, consistent with another study [78]. Although the pyruvate fermentation pathway does not release any ammonia, it does release butyrate, which can be harmful in conditions like chronic periodontitis [79]. On the other hand, palmitate biosynthesis II, palmitoleate biosynthesis I, and oleate biosynthesis IV significantly decreased COVID-19 infections. Some studies have shown that palmitate is essential in maintaining gut barrier integrity by controlling MUC2 secretion and function [80]. Another recent study found that palmitate increases the production of MUC2 in the goblet cells of the gut, leading to the formation of a thick mucus gel and maintaining the integrity of the gut barrier [81]. In addition to identifying metabolic pathways, we have identified the top 10 bacterial genes (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) as bacterial key genes (bKGs) capable of facilitating secondary infection with SARS-CoV-2 through protein–protein interaction network analysis. Previous studies have highlighted the importance of certain genes in lipid metabolism, such as accB, which is significantly associated with diabetes and obesity in the Pakistani population [82]. On the other hand, another study explored that ppc genes are significantly associated with the metabolic disease hyperproteinemia, which is also associated with COVID-19 and increases the severity and mortality rate [83,84].
To explore the potential drug molecules, we conducted molecular docking corresponding to the identified bKGs. We selected the eight top-ranked ligands (bemcentinib, ledipasvir, velpatasvir, tirilazad, acetyldigitoxin, entreatinib, digitoxin, and elbasvir) based on their binding affinity scores, which were also supported by the independent receptors. Of the eight proposed drug molecules, only Ledipasvir is FDA-approved for COVID-19 [85]. Bemcentinib is FDA-approved for acute myeloid leukemia and passed clinical phase 2 trials for COVID-19 [86]. Velpatasvir is FDA-approved for the Hepatitis C virus (HCV) [87] and is under clinical trial for COVID-19. Tirilazad is an FDA-approved drug for brain cancer and spinal cord injury [88] and is under clinical trial for COVID-19. Acetyldigitoxin is used for the treatment of heart failure (HF) and is passing clinical trial phases for COVID-19. Digitoxin is passing a clinical phase 1 trial for COVID-19 [9]. Telcagepant is in a clinical phase 3 trial for acute migraine [89], but it is not yet under the clinical phase for COVID-19. Elbasvir is an FDA-approved drug for the treatment of hepatitis C virus genotype 1b (HCV GT1b) patients [90], and it is supported for COVID-19 by other studies [91], but is not yet considered for a clinical trial for COVID-19. To compare our suggested eight drug molecules against the two FDA-approved drugs (Molnupiravir and Nirmatrelvir) for COVID-19, we investigated their docking scores (kcal/mol) with our proposed and top-ranked published target proteins and found that our proposed drug molecules produced higher binding affinity scores in both cases. Interestingly, the top three drug molecules were the same against our proposed and published targets, and hence, we can conclude that these top three drug molecules, bemcentinib, ledipasvir, and velpatasvir, might be the most effective ligands for the treatment of COVID-19 infections. A study demonstrated that bemcentinib is efficacious in decreasing viral infection in lung cells, has a significant impact on SARS-CoV2 infection [92], and is also more effective in cancer treatment [93]. This drug molecule has no toxicity in the human body and is safe for treatment [92]. It prevents SARS-CoV-2 infections, as shown by viral transcripts in RNAseq studies and viral load in qRT-PCR studies of human lung epithelium, Vero-E6, and A549-hACE2 cells [94]. It reduces viral internalization but has no impact on viral binding. Ledipasvir is an approved drug for both COVID-19 and hepatitis C virus treatment with a good safety profile [95,96]. The antiviral activity of ledipasvir was identified in a cell-based screening assay against SARS-CoV-2 [97,98]. A study showed that the velpatasvir molecule produced the best binding affinity score against COVID-19 protease, which is also more effective against the hepatitis C virus treatment [99]. Furthermore, toxicity testing using the ADMET approach showed that the three complexes are the least toxic and the safest possible drugs for the treatment of SARS-CoV-2 infection. In addition, molecular dynamics (MD) simulation studies confirmed the structural stability of drug–target complexes by the RMSD, binding free energy, radius of gyration, and SASA parameters. So, our proposed drug molecules may have potential for the treatment of bacterial co-/super-/secondary infections with SARS-CoV-2 infections.

4. Limitations and Commercial Applications

This study was entirely computational-based. Therefore, the outputs of this study, including co-/secondary infection-causing bacterial key genes (bKGs) and candidate therapeutic agents, require experimental validation in a wet lab before going into the commercial production of effective drugs in order to inhibit bacterial pathogens for COVID-19.

5. Materials and Methods

5.1. Data Source and Description

In this study, we considered bacterial 16S rRNA sequence profiles and metadata on drug molecules that are associated with COVID-19, as described below.

5.1.1. Collection of 16S rRNA Sequence Data

Bacterial 16S rRNA sequence profiles and associated metadata were downloaded from the NCBI online database with the bio-project number PRJNA684070. These samples were collected from the oral and gut microbiome signatures of hospitalized patients with COVID-19 and healthy individuals in Guangdong, China [55]. The dataset consisted of 297 samples, including 76 healthy and 221 COVID-19 patients. Previously, this dataset was analyzed to explore the bacterial taxa associated with SARS-CoV-2 infections for oral and gut samples separately [55]. In this study, we analyzed oral and gut samples jointly to explore common bacterial taxa/groups and their key genes that are associated with the severity of SARS-CoV-2 infections.

5.1.2. Collection of Metadata on Drug Molecules

To explore bacterial key genes (bKGs) guided candidate drug molecules by molecular docking simulation for the treatment against bacterial co-/super-/secondary infections with SARS-CoV-2 infections, we collected 786 meta-drug agents from different published articles on COVID-19 (see Table S4 in Supplementary Materials).

5.2. Statistics and Bioinformatics Analysis

5.2.1. Preprocessing of 16S rRNA Sequence Profiles

At first, we checked the quality of the raw sequence reads using FASTQC and then filtered the poor-quality reads using Trimmomatic-0.39 [100] with default parameters. After trimming the poor-quality reads, on average, almost 97% of the sequence reads survived for further study. Then, NGmerge (v0.2) was used to combine the overlapping areas between the paired-end trimmed reads with a minimum 5 bp overlap and a maximum of 10% mismatches [101]. Finally, the high-quality integrated reads were clustered into operational taxonomic units (OTUs), defined at 97% sequence similarity compared to the Greengenes reference database via Qiime2 software (version: 2023.9). A total of 12,671,011 OTUs with 42,663.33 mean numbers of OTUs per sample were obtained. For a large number of OTUs, filter the OTUs table with the minimum number of 5 samples that a feature observes. The filtered table contained 12,288,422 OTUs with a mean number of 41,797.35 (range: 80937-1778) OTUs per sample. Then, the OTU table was reduced to the lowest value of the total number of individuals observed in all samples.

5.2.2. Diversity Analysis

The bacterial diversities within samples were investigated by using alpha diversity, measured by using four alpha diversity indices: observed, Chao1, Shannon, and Simpson [102], determined using phyloseq [103] R-package after trimming the minimum counts of all samples. The Chao1 and observed species are used to measure bacterial richness, while Shannon and Simpson quantitative indices measure species richness and evenness [104]. The significant differences across samples were tested by the nonparametric Wilcoxon sign rank test and the results were visualized using the ggbox R package [105]. To assess differences in species complexity between samples, beta diversity was used and measured using the Bray–Curtis distance based on normalization and log10 transformed rarefied OTU abundances [106]. In the R package, MicrobiotaProcess was used to perform principal coordinate analysis (PCoA) to obtain principal coordinates [107]. The variations in bacterial composition between the COVID-19 and healthy groups were investigated using PERMANOVA with Bray–Curtis distance through the vegan R-package with permutations set to 999 [108].

5.2.3. Taxonomy Analysis and Identification of Differentially Abundant Bacterial Compositions (DABCs)

The Ribosomal Database Project (RDP Version 2.10) classifier was used to assign each representative sequence of OTU-clusters to the bacterial taxa in the Greengenes database (version: 13_5) by using Qiime2 (version: 2023.9) [109]. While the RDP classifier is based on the naïve Bayes model, which assigns a sequence to the closest match using a posterior score, we used linear discriminant analysis (LDA) effect size (LEfSe) on the Galaxy platform to identify the intergroup discriminant bacterial features at phylum and genus level across the patient groups (p < 0.05, LDA score > 3.5) “http://huttenhowersph.harva rd.edu/galaxy (accessed on 5 September 2023)” [110]. LEfSe uses the two-tailed nonparametric Kruskal–Wallis test to evaluate the significance of differences in OTUs in case–control groups. Additionally, a zero-inflated Gaussian mixture (ZIGM) model of mean group frequencies was applied to normalized count data to detect differentially abundant features or genera between case–control groups at a fixed p-value (p < 0.05) [111]. Bacterial count data were aggregated to the genus level, and abundant genera (>30 normalized counts per sample) were used for the ZIGM model. This analysis was performed using the R software package, version 4.3.3. The explored differential genera’s were employed as input for the functional analysis.

5.2.4. Identification of Pathway-Based Bacterial Key Genes (bKGs) from DABCs

The content of functional genes in the gut microbiota was predicted using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) software based on the identified features (differentially abundant bacterial compositions) [112]. PICRUSt2 predictions were based on the enzyme classification numbers and the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs (KO) (v77.1). This functional analysis was utilized to investigate all of the samples’ metabolic pathways and genes. After that, we used Statistical Analysis of Metagenomic Profiles (STAMP) software [113] to identify differentially abundant metabolic pathways between SARS-CoV-2 and control patients and then collected all genes from these identified pathways. To identify differential metabolic pathways between the two groups, we utilized Welch’s t statistic to produce q-values for the two-tailed test, and Welch’s inverted confidence interval method was applied to calculate the confidence interval. The Benjamini–Hochberg false discovery rate method was then applied to compute the adjusted p-value (p < 0.05). Then, microbial protein–protein interaction networking of those pathway genes was performed by using the STRING database [114], and the result was visualized by using the cytoscape software [115] to identify the top-ranked bacterial key genes (bKGs) that were used as drug targets for further molecular docking analysis (see Figure S1).

5.2.5. Bacterial Key Genes Guided Drug Repurposing by Molecular Docking

Molecular docking analysis was performed to repurpose potential drug molecules by using AutoDock Vina [116] for the treatment against bacterial co-/super-/secondary infections with SARS-CoV-2 infections. Initially, all receptor protein structures mediated by bKGs were downloaded from the Protein Data Bank (PDB), and heteroatoms, water molecules, and attached ligands were removed using BIOVIA Discovery Studio [117] and PyMOL [118]. Subsequently, protein energies were minimized with SwissPdb Viewer [119] and finally converted to an acceptable PDBQT format by fixing the grid box center in AutoDock tools. The ligand structures were energy-minimized using the MMFF94 force field in the Avogadro software [120] with specific parameters (total steps: 200, update steps: 1, energy difference: 0.1) [121]. Minimized ligands were processed and then converted to the PDBQT format using AutoDock tools, and protein–ligand complexes had their ionizable group pKa values estimated with PROPKA at physiological pH 7.0 [122]. Docking analysis was performed using AutoDock Vina (version 1.1.2) with an exhaustiveness parameter set to 10. We visualized the top-ranked drug–target complexes by arranging ligands and targets based on descending average binding affinity scores (kcal/mol) using BIOVIA Discovery Studio (version 3.0) and PyMOL software (version 2.3) to verify non-bond interactions. Furthermore, we conducted pharmacokinetic and toxicological assessments of the top-ranked drug molecules using their canonical SMILES as input on three online web servers: admetSAR [123], SwissADME [124], and pKCSM [125].
To evaluate the overall stability of the top three complexes, molecular dynamic simulations were conducted using YASARA software [126] and the AMBER14 force field [127], neutralized at 298 K, pH 7.4, and 0.9 percent NaCl. Energy minimization and annealing methods were applied before simulating for 100 ns, with trajectory snapshots saved every 100 ps. MM-Poisson–Boltzmann Surface Area (MM-PBSA) binding free energy for the complexes was calculated using captured trajectory snapshots [128].
Binding Energy = EpotRecept + EsolvRecept + EpotLigand + EsolvLigand − EpotComplex − EsolvComplex
In this study, the MM-PBSA binding free energy was calculated for the top-ranked 3 complexes using the YASARA macro, where positive energy signifies better binding [129]. Simulation snapshots were used to calculate the root-mean-square deviation (RMSD), the radius of gyration, and the surface-accessible surface area (SASA) [130,131,132,133,134]. The entire workflow of this study is displayed in Figure 9.

6. Conclusions

This study attempted to explore co-/super-/secondary infection-causing bacterial pathogens and their inhibitors for SARS-CoV-2 infection by using integrated statistical and bioinformatics approaches. At first, we analyzed the 16S rRNA bacterial sequence profiles of 297 individuals, of which 221 were COVID-19 patients and 76 were healthy, using Qiime2. It produced an OTU table of the count data between bacteria and patient samples. The diversity analysis based on the OTU table showed that SARS-CoV-2-infected samples are significantly separated from the control samples by the changes in the bacterial community. Several differentially abundant genus-taxa (Bacteroides, Parabacteroides, Prevotella Clostridium, Atopobium, and Peptostreptococcus) between SARS-CoV-2 infection and control samples were identified by using the OUT-table with the taxon-data table in the metagenomeSeq r-package. Then, we detected 23 significantly differentially abundant metabolic pathways based on the differentially abundant bacterial genus by using STAMP software. Among them, ‘incomplete reductive TCA’ and ‘pyruvate fermentation to acetate and lactate II’ pathways were significantly upregulated: palmitate biosynthesis II, palmitoleate biosynthesis I, and oleate biosynthesis IV pathways, were significantly downregulated. We obtained 183 bacterial genes from the identified differentially abundant pathways. Then, we selected the top-ranked 10 bKGs (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) by the protein–protein interaction (PPI) network analysis of those 183 pathway genes. Then, we detected the bKGs-guided top-ranked 10 anti-SARS-CoV-2 drug molecules by molecular docking. Similarly, the SARS-CoV-2 infection-causing top-ranked nine published key gene (ACE2, RdRp, 3CLpro, S, TMPRSS2, PLpro, IL6, TNF, and N)-guided the top-ranked 10 anti-SARS-CoV-2 drug molecules were also identified. Then, we considered the top-ranked eight common molecules (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Entreatinib, Digitoxin, and Elbasvir) out of ten as the candidate anti-SARS-CoV-2 drug molecules. Finally, the binding stability of the top-ranked three drug–target complexes was investigated by molecular dynamic (MD) simulation studies and found to be stable. Our in silico pipeline might be a useful guideline to explore potential targets and ligands for taking a proper treatment plan against microbial co-/secondary infections with other diseases as well.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph17040432/s1. Table S1: Comparison of different alpha diversity indices of bacterial communities with a significant difference between COVID-19 patients (n = 221) and healthy subjects (n = 76) based on Wilcoxon-Mann–Whitney test p values; Table S2: Comparative relative abundances of the oral and gut microbiota at the phylum and genus levels in the COVID-19 and healthy groups. Firmicutes, Bacteroidetes, Proteobacteria, and Fusobacteria were the most abundant phyla between the two groups, whereas Bacteroides, Prevotella, Actinomyces, Rothia, and Fusobacterium were the most dominant genera in the microbiome of the sample; Table S3: Pharmacological assessment of the top three potential ligand molecules derived from the admetSAR, SwissADME, and pKCSM web servers; Table S4: Metadata of SARS-CoV-2 infection obtained by reviewing published articles were used in this study; Table S5: Binding affinity score of two FDA-approved drugs corresponding to our proposed and published reviewed targets; Table S6: Non-bond interactions between top-ordered three receptors and ligand compounds based on their binding affinity; Figure S1: A workflow for exploring potential functional genes and their metabolic pathways from differentially abundant bacterial genera between SARS-CoV-2 infection and control samples. All the references in the supplementary materials have been mentioned in this paper.

Author Contributions

M.K.K. and M.N.H.M. conceived the idea of the study. M.K.K. analyzed 16S rRNA sequence profiles using statistics and other upstream and downstream analyses using different bioinformatics tools and M.K.K., M.A.A. and M.A.I. performed the drug screening steps by molecular docking and M.K.K. writing the whole manuscript. M.N.H.M., M.Y., I.A.K., M.A.B. and R.A.M. edited the manuscript. M.N.H.M. supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the research supporting project number (RSPD2024R740), King Saud University, Riyadh, Saudi Arabia, and BANBEIS research project number (MS20191106), Ministry of Education, Bangladesh.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw 16S rRNA sequence profile data set analyzed in this study is publicly available. It can be freely downloaded from the online NCBI database with bio-project number PRJNA684070. https://www.ncbi.nlm.nih.gov/bioproject/PRJNA684070 (access on 25 May 2023).

Acknowledgments

We are thankful to the researchers supporting project number (RSPD2024R740), King Saud University, Riyadh, Saudi Arabia and BANBEIS research project number (MS20191106), Ministry of Education, Bangladesh. We also would like to acknowledge the editor and all reviewers for their valuable comments and suggestions that help us to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. COVID-19 Epidemiological Update—15 March 2024; Epidemiology; WHO: Geneva, Switzerland, 2024; pp. 6–8. [Google Scholar]
  2. Markov, P.V.; Ghafari, M.; Beer, M.; Lythgoe, K.; Simmonds, P.; Stilianakis, N.I.; Katzourakis, A. The Evolution of SARS-CoV-2. Nat. Rev. Microbiol. 2023, 21, 361–379. [Google Scholar] [CrossRef] [PubMed]
  3. Ahmad, S.; Abbasi, H.W.; Shahid, S.; Gul, S.; Abbasi, S.W. Molecular Docking, Simulation and MM-PBSA Studies of Nigella Sativa Compounds: A Computational Quest to Identify Potential Natural Antiviral for COVID-19 Treatment. J. Biomol. Struct. Dyn. 2021, 39, 4225–4233. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, Y.; Liu, Q.; Guo, D. Emerging Coronaviruses: Genome Structure, Replication, and Pathogenesis. J. Med. Virol. 2020, 92, 418–423. [Google Scholar] [CrossRef] [PubMed]
  5. Mengist, H.M.; Dilnessa, T.; Jin, T. Structural Basis of Potential Inhibitors Targeting SARS-CoV-2 Main Protease. Front. Chem. 2021, 9, 622898. [Google Scholar] [CrossRef] [PubMed]
  6. Huff, S.; Kummetha, I.R.; Tiwari, S.K.; Huante, M.B.; Clark, A.E.; Wang, S.; Bray, W.; Smith, D.; Carlin, A.F.; Endsley, M.; et al. Discovery and Mechanism of SARS-CoV-2 Main Protease Inhibitors. J. Med. Chem. 2022, 65, 2866–2879. [Google Scholar] [CrossRef] [PubMed]
  7. Vicenti, I.; Zazzi, M.; Saladini, F. SARS-CoV-2 RNA-Dependent RNA Polymerase as a Therapeutic Target for COVID-19. Expert Opin. Ther. Pat. 2021, 31, 325–337. [Google Scholar] [CrossRef] [PubMed]
  8. Mishra, A.; Rathore, A.S. RNA Dependent RNA Polymerase (RdRp) as a Drug Target for SARS-CoV-2. J. Biomol. Struct. Dyn. 2022, 40, 6039–6051. [Google Scholar] [CrossRef] [PubMed]
  9. Beck, B.R.; Shin, B.; Choi, Y.; Park, S.; Kang, K. Predicting Commercially Available Antiviral Drugs That May Act on the Novel Coronavirus (SARS-CoV-2) through a Drug-Target Interaction Deep Learning Model. Comput. Struct. Biotechnol. J. 2020, 18, 784–790. [Google Scholar] [CrossRef] [PubMed]
  10. Pałka, A.; Kujawska, A.; Hareza, D.; Gajda, M.; Wordliczek, J.; Jachowicz-Matczak, E.; Owsianka, I.; Żółtowska, B.; Chmielarczyk, A.; Romaniszyn, D.; et al. Secondary Bacterial Infections & Extensively Drug-Resistant Bacteria among COVID-19 Hospitalized Patients at the University Hospital in Kraków. Ann. Clin. Microbiol. Antimicrob. 2023, 22, 77. [Google Scholar] [CrossRef]
  11. Bedenić, B.; Bratić, V.; Mihaljević, S.; Lukić, A.; Vidović, K.; Reiner, K.; Schöenthaler, S.; Barišić, I.; Zarfel, G.; Grisold, A. Multidrug-Resistant Bacteria in a COVID-19 Hospital in Zagreb. Pathogens 2023, 12, 117. [Google Scholar] [CrossRef]
  12. Wu, H.Y.; Chang, P.H.; Chen, K.Y.; Lin, I.F.; Hsih, W.H.; Tsai, W.L.; Chen, J.A.; Lee, S.S.J. Coronavirus Disease 2019 (COVID-19) Associated Bacterial Coinfection: Incidence, Diagnosis and Treatment. J. Microbiol. Immunol. Infect. 2022, 55, 985–992. [Google Scholar] [CrossRef] [PubMed]
  13. Ritter, L.A.; Britton, N.; Heil, E.L.; Teeter, W.A.; Murthi, S.B.; Chow, J.H.; Ricotta, E.; Chertow, D.S.; Grazioli, A.; Levine, A.R. The Impact of Corticosteroids on Secondary Infection and Mortality in Critically Ill COVID-19 Patients. J. Intensive Care Med. 2021, 36, 1201–1208. [Google Scholar] [CrossRef]
  14. Langford, B.J.; So, M.; Raybardhan, S.; Leung, V.; Westwood, D.; MacFadden, D.R.; Soucy, J.P.R.; Daneman, N. Bacterial Co-Infection and Secondary Infection in Patients with COVID-19: A Living Rapid Review and Meta-Analysis. Clin. Microbiol. Infect. 2020, 26, 1622–1629. [Google Scholar] [CrossRef] [PubMed]
  15. Gomez, A.C.; Ortiz, T.; Valenzuela, A.; Egoávil-Espejo, R.; Huerto-Huanuco, R.; Pinto, J.A.; Lagos, J.; Ruiz, J. Super-Infection by Multiple Microorganisms in COVID-19 Patients. Front. Mol. Biosci. 2023, 10, 1113969. [Google Scholar] [CrossRef] [PubMed]
  16. Cataño-Correa, J.C.; Cardona-Arias, J.A.; Mancilla, J.P.P.; García, M.T. Bacterial Superinfection in Adults with COVID-19 Hospitalized in Two Clinics in Medellín-Colombia, 2020. PLoS ONE 2021, 16, e0254671. [Google Scholar] [CrossRef]
  17. Feldman, C.; Anderson, R. The Role of Co-Infections and Secondary Infections in Patients with COVID-19. Pneumonia 2021, 13, 5. [Google Scholar] [CrossRef]
  18. Chong, W.H.; Saha, B.K.; Ramani, A.; Chopra, A. State-of-the-Art Review of Secondary Pulmonary Infections in Patients with COVID-19 Pneumonia. Infection 2021, 49, 591–605. [Google Scholar] [CrossRef] [PubMed]
  19. Ding, F.M.; Feng, Y.; Han, L.; Zhou, Y.; Ji, Y.; Hao, H.J.; Xue, Y.S.; Yin, D.N.; Xu, Z.C.; Luo, S.; et al. Early Fever Is Associated with Clinical Outcomes in Patients with Coronavirus Disease. Front. Public Health 2021, 9, 712190. [Google Scholar] [CrossRef] [PubMed]
  20. Song, W.J.; Hui, C.K.; Hull, J.H.; Birring, S.S.; McGarvey, L.; Mazzone, S.B.; Chung, K.F. Confronting COVID-19-Associated Cough and the Post-COVID Syndrome: Role of Viral Neurotropism, Neuroinflammation, and Neuroimmune Responses. Lancet Respir. Med. 2021, 9, 533–544. [Google Scholar] [CrossRef]
  21. Mahase, E. COVID-19: Sore Throat, Fatigue, and Myalgia Are More Common with New UK Variant. BMJ 2021, 372, n288. [Google Scholar] [CrossRef]
  22. Zhang, Z. Diagnosis and Treatment Protocol for COVID-19 Patients (Trial Version 9). Health Care Sci. 2022, 1, 14–28. [Google Scholar] [CrossRef]
  23. Rain, M.; Puri, G.D.; Bhalla, A.; Avti, P.; Subramaniam, B.; Kaushal, V.; Srivastava, V.; Mahajan, P.; Singh, M.; Pandey, N.; et al. Effect of Breathing Intervention in Patients with COVID and Healthcare Workers. Front. Public Health 2022, 10, 945988. [Google Scholar] [CrossRef] [PubMed]
  24. Gonçalves, L.F.; Gonzáles, A.I.; Paiva, K.M.; Patatt, F.S.A.; Stolz, J.V.; Haas, P. Smell and Taste Alterations in COVID-19 Patients: A Systematic Review. Rev. Assoc. Med. Bras. 2021, 66, 1602–1608. [Google Scholar] [CrossRef] [PubMed]
  25. Widyadharma, I.P.E.; Sari, N.N.S.P.; Pradnyaswari, K.E.; Yuwana, K.T.; Adikarya, I.P.G.D.; Tertia, C.; Wijayanti, I.A.S.; Indrayani, I.A.S.; Utami, D.K.I. Pain as Clinical Manifestations of COVID-19 Infection and Its Management in the Pandemic Era: A Literature Review. Egypt. J. Neurol. Psychiatry Neurosurg. 2020, 56, 121. [Google Scholar] [CrossRef] [PubMed]
  26. Daria, S.; Katarzyna, S.A.; Jerzy, F.; Ivan, K. Muscle Pain and Muscle Weakness in COVID19 Patients: Cross-Talk with Statins—Preliminary Results. Biomed. Pharmacother. 2022, 148, 112757. [Google Scholar] [CrossRef] [PubMed]
  27. Mutiawati, E.; Kusuma, H.I.; Fahriani, M.; Harapan, H.; Syahrul, S.; Musadir, N. Headache in Post-COVID-19 Patients: Its Characteristics and Relationship with the Quality of Life. Medicina 2022, 58, 1500. [Google Scholar] [CrossRef] [PubMed]
  28. D’Amico, F.; Baumgart, D.C.; Danese, S.; Peyrin-Biroulet, L. Diarrhea during COVID-19 Infection: Pathogenesis, Epidemiology, Prevention, and Management. Clin. Gastroenterol. Hepatol. 2020, 18, 1663–1672. [Google Scholar] [CrossRef]
  29. Hurn, E.; Dickinson, L.; Abraham, J.A. Bacterial Meningitis and COVID-19: A Complex Patient Journey. BMJ Case Rep. 2021, 14, e239533. [Google Scholar] [CrossRef]
  30. Richards, T.; Stephen, J.; Lui, C.L. Severe Disseminated Veillonella parvula Infection Including Endocarditis, Bilateral Psoas Abscess, Discitis, and Osteomyelitis but Sparing Spinal and Hip Prostheses: A Case Report. J. Med. Case Rep. 2022, 16, 157. [Google Scholar] [CrossRef]
  31. Attaway, A.H.; Scheraga, R.G.; Bhimraj, A.; Biehl, M.; Hatipoğlu, U. Severe COVID-19 Pneumonia: Pathogenesis and Clinical Management. BMJ 2021, 372, n436. [Google Scholar] [CrossRef]
  32. Seong, Y.J.; Lee, S.H.; Kim, E.J.; Choi, Y.H.; Kim, T.J.; Lee, W.G.; Heo, J.Y. Campylobacter fetus Subspecies venerealis Meningitis Associated with a Companion Dog in a Young Adult: A Case Report. BMC Infect. Dis. 2021, 21, 1280. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, H.; Wang, Q.; Yang, C.; Guo, M.; Cui, X.; Jing, Z.; Liu, Y.; Qiao, W.; Qi, H.; Zhang, H.; et al. Bacteroides Acidifaciens in the Gut Plays a Protective Role against CD95-Mediated Liver Injury. Gut Microbes 2022, 14, 2027853. [Google Scholar] [CrossRef] [PubMed]
  34. Alshaikh, F.S.; Baker, A.I. Prevalence of Bacterial Coinfection and Patterns of Antibiotics Prescribing in Patients with COVID-19: A Systematic Review and Meta-Analysis. Int. J. Pharm. Pract. 2022, 30, i21. [Google Scholar] [CrossRef]
  35. Mosharaf, M.P.; Reza, M.S.; Kibria, M.K.; Ahmed, F.F.; Kabir, M.H.; Hasan, S.; Mollah, M.N.H. Computational Identification of Host Genomic Biomarkers Highlighting Their Functions, Pathways and Regulators That Influence SARS-CoV-2 Infections and Drug Repurposing. Sci. Rep. 2022, 12, 4279. [Google Scholar] [CrossRef]
  36. Mosharaf, M.P.; Kibria, M.K.; Hossen, M.B.; Islam, M.A.; Reza, M.S.; Mahumud, R.A.; Alam, K.; Gow, J.; Mollah, M.N.H. Meta-Data Analysis to Explore the Hub of the Hub-Genes That Influence SARS-CoV-2 Infections Highlighting Their Pathogenetic Processes and Drugs Repurposing. Vaccines 2022, 10, 1248. [Google Scholar] [CrossRef] [PubMed]
  37. Ahmed, F.F.; Reza, M.S.; Sarker, M.S.; Islam, M.S.; Mosharaf, M.P.; Hasan, S.; Mollah, M.N.H. Identification of Host Transcriptome-Guided Repurposable Drugs for SARS-CoV-1 Infections and Their Validation with SARS-CoV-2 Infections by Using the Integrated Bioinformatics Approaches. PLoS ONE 2022, 17, e0266124. [Google Scholar] [CrossRef]
  38. Sarker, B.; Rahaman, M.M.; Islam, M.A.; Alamin, M.H.; Husain, M.M.; Ferdousi, F.; Ahsan, M.A.; Mollah, M.N.H. Identification of Host Genomic Biomarkers from Multiple Transcriptomics Datasets for Diagnosis and Therapies of SARS-CoV-2 Infections. PLoS ONE 2023, 18, e0281981. [Google Scholar] [CrossRef]
  39. Weigert Muñoz, A.; Zhao, W.; Sieber, S.A. Monitoring Host-Pathogen Interactions Using Chemical Proteomics. RSC Chem. Biol. 2023, 5, 73–89. [Google Scholar] [CrossRef]
  40. Mahmud, S.; Hasan, M.R.; Biswas, S.; Paul, G.K.; Afrose, S.; Mita, M.A.; Sultana Shimu, M.S.; Promi, M.M.; Hani, U.; Rahamathulla, M.; et al. Screening of Potent Phytochemical Inhibitors Against SARS-CoV-2 Main Protease: An Integrative Computational Approach. Front. Bioinforma. 2021, 1, 717141. [Google Scholar] [CrossRef]
  41. Crua Asensio, N.; Munõz Giner, E.; De Groot, N.S.; Torrent Burgas, M. Centrality in the Host-Pathogen Interactome Is Associated with Pathogen Fitness during Infection. Nat. Commun. 2017, 8, 14092. [Google Scholar] [CrossRef]
  42. Balint, D.; Brito, I.L. Human–Gut Bacterial Protein–Protein Interactions: Understudied but Impactful to Human Health. Trends Microbiol. 2023. [Google Scholar] [CrossRef]
  43. Xu, R.; Lu, R.; Zhang, T.; Wu, Q.; Cai, W.; Han, X.; Wan, Z.; Jin, X.; Zhang, Z.; Zhang, C. Temporal Association between Human Upper Respiratory and Gut Bacterial Microbiomes during the Course of COVID-19 in Adults. Commun. Biol. 2021, 4, 240. [Google Scholar] [CrossRef]
  44. Rafiqul Islam, S.M.; Foysal, M.J.; Hoque, M.N.; Mehedi, H.M.H.; Rob, M.A.; Salauddin, A.; Tanzina, A.Y.; Biswas, S.; Noyon, S.H.; Siddiki, A.M.A.M.Z.; et al. Dysbiosis of Oral and Gut Microbiomes in SARS-CoV-2 Infected Patients in Bangladesh: Elucidating the Role of Opportunistic Gut Microbes. Front. Med. 2022, 9, 821777. [Google Scholar] [CrossRef] [PubMed]
  45. Moni, M.A.; Quinn, J.M.W.; Sinmaz, N.; Summers, M.A. Gene Expression Profiling of SARS-CoV-2 Infections Reveal Distinct Primary Lung Cell and Systemic Immune Infection Responses That Identify Pathways Relevant in COVID-19 Disease. Brief. Bioinform. 2021, 22, 1324–1337. [Google Scholar] [CrossRef] [PubMed]
  46. Zhou, Y.; Zhang, J.; Zhang, D.; Ma, W.L.; Wang, X. Linking the Gut Microbiota to Persistent Symptoms in Survivors of COVID-19 after Discharge. J. Microbiol. 2021, 59, 941–948. [Google Scholar] [CrossRef]
  47. Gaibani, P.; D’Amico, F.; Bartoletti, M.; Lombardo, D.; Rampelli, S.; Fornaro, G.; Coladonato, S.; Siniscalchi, A.; Re, M.C.; Viale, P.; et al. The Gut Microbiota of Critically Ill Patients with COVID-19. Front. Cell. Infect. Microbiol. 2021, 11, 670424. [Google Scholar] [CrossRef]
  48. Hua, Q.; Zhang, B.; Xu, G.; Wang, L.; Wang, H.; Lin, Z.; Yu, D.; Ren, J.; Zhang, D.; Zhao, L.; et al. CEMIP, a Novel Adaptor Protein of OGT, Promotes Colorectal Cancer Metastasis through Glutamine Metabolic Reprogramming via Reciprocal Regulation of β-Catenin. Oncogene 2021, 40, 6443–6455. [Google Scholar] [CrossRef]
  49. Uehara, O.; Abiko, Y.; Nagasawa, T.; Morikawa, T.; Hiraki, D.; Harada, F.; Kawano, Y.; Toraya, S.; Matsuoka, H.; Paudel, D.; et al. Alterations in the Oral Microbiome of Individuals with a Healthy Oral Environment Following COVID-19 Vaccination. BMC Oral Health 2022, 22, 50. [Google Scholar] [CrossRef] [PubMed]
  50. Perrone, M.R.; Romano, S.; De Maria, G.; Tundo, P.; Bruno, A.R.; Tagliaferro, L.; Maffia, M.; Fragola, M. Compositional Data Analysis of 16S RRNA Gene Sequencing Results from Hospital Airborne Microbiome Samples. Int. J. Environ. Res. Public Health 2022, 19, 10107. [Google Scholar] [CrossRef]
  51. Shi, Y.L.; He, M.Z.; Han, M.Z.; Gui, H.Y.; Wang, P.; Yu, J.L.; Ge, Y.L.; Sun, Y.; Huang, S.H. Characterization of Altered Oropharyngeal Microbiota in Hospitalized Patients with Mild SARS-CoV-2 Infection. Front. Cell. Infect. Microbiol. 2022, 12, 824578. [Google Scholar] [CrossRef]
  52. Mazzarelli, A.; Giancola, M.L.; Farina, A.; Marchioni, L.; Rueca, M.; Gruber, C.E.M.; Bartolini, B.; Bartoli, T.A.; Maffongelli, G.; Capobianchi, M.R.; et al. 16S RRNA Gene Sequencing of Rectal Swab in Patients Affected by COVID-19. PLoS ONE 2021, 16, e0247041. [Google Scholar] [CrossRef]
  53. Zuo, T.; Zhang, F.; Lui, G.C.Y.; Yeoh, Y.K.; Li, A.Y.L.; Zhan, H.; Wan, Y.; Chung, A.C.K.; Cheung, C.P.; Chen, N.; et al. Alterations in Gut Microbiota of Patients with COVID-19 during Time of Hospitalization. Gastroenterology 2020, 159, 944–955.e948. [Google Scholar] [CrossRef] [PubMed]
  54. Al Bataineh, M.T.; Henschel, A.; Mousa, M.; Daou, M.; Waasia, F.; Kannout, H.; Khalili, M.; Kayasseh, M.A.; Alkhajeh, A.; Uddin, M.; et al. Gut Microbiota Interplay with COVID-19 Reveals Links to Host Lipid Metabolism among Middle Eastern Populations. Front. Microbiol. 2021, 12, 761067. [Google Scholar] [CrossRef] [PubMed]
  55. Wu, Y.; Cheng, X.; Jiang, G.; Tang, H.; Ming, S.; Tang, L.; Lu, J.; Guo, C.; Shan, H.; Huang, X. Altered Oral and Gut Microbiota and Its Association with SARS-CoV-2 Viral Load in COVID-19 Patients during Hospitalization. npj Biofilms Microbiomes 2021, 7, 61. [Google Scholar] [CrossRef] [PubMed]
  56. Vila Domínguez, A.; Enrique Jiménez Mejías, M.; Smani, Y. Drugs Repurposing for Multi-Drug Resistant Bacterial Infections. Drug Repurposing—Hypothesis Mol. Asp. Ther. Appl. 2020, 6, 97–118. [Google Scholar] [CrossRef]
  57. Abubakar, U.; Al-Anazi, M.; Alanazi, Z.; Rodríguez-Baño, J. Impact of COVID-19 Pandemic on Multidrug Resistant Gram Positive and Gram Negative Pathogens: A Systematic Review. J. Infect. Public Health 2023, 16, 320–331. [Google Scholar] [CrossRef] [PubMed]
  58. Clarridge, J.E., 3rd. Impact of 16S RRNA Gene Sequence Analysis for Identification of Bacteria on Clinical Microbiology and Infectious Diseases. Clin. Microbiol. Rev. 2004, 17, 840–862. [Google Scholar] [CrossRef] [PubMed]
  59. Sharma Manjul, A.; Shirkot, P.; Anshul Sharma Manjul, C. 6S RRNA Gene Sequencing for Bacterial Identification of Pullulanase Synthesizing Thermophilic Bacteria Contributing to Big Data. Int. J. Chem. Stud. 2018, 6, 2769–2773. [Google Scholar]
  60. Syed, Y.Y. Molnupiravir: First Approval. Drugs 2022, 82, 455–460. [Google Scholar] [CrossRef]
  61. Lamb, Y.N. Nirmatrelvir Plus Ritonavir: First Approval. Drugs 2022, 82, 585–591. [Google Scholar] [CrossRef] [PubMed]
  62. Magne, F.; Gotteland, M.; Gauthier, L.; Zazueta, A.; Pesoa, S.; Navarrete, P.; Balamurugan, R. The Firmicutes/Bacteroidetes Ratio: A Relevant Marker of Gut Dysbiosis in Obese Patients? Nutrients 2020, 12, 1474. [Google Scholar] [CrossRef]
  63. Sears, C.L.; Islam, S.; Saha, A.; Arjumand, M.; Alam, N.H.; Faruque, A.S.G.; Salam, M.A.; Shin, J.; Hecht, D.; Weintraub, A.; et al. Association of Enterotoxigenic Bacteroides Fragilis Infection with Inflammatory Diarrhea. Clin. Infect. Dis. 2008, 47, 797–803. [Google Scholar] [CrossRef]
  64. Kaur, U.S.; Shet, A.; Rajnala, N.; Gopalan, B.P.; Moar, P.; Himanshu, D.; Singh, B.P.; Chaturvedi, R.; Tandon, R. High Abundance of Genus Prevotella in the Gut of Perinatally HIV-Infected Children Is Associated with IP-10 Levels despite Therapy. Sci. Rep. 2018, 8, 17679. [Google Scholar] [CrossRef] [PubMed]
  65. Ramsheh, M.Y.; Haldar, K.; Esteve-Codina, A.; Purser, L.F.; Richardson, M.; Müller-Quernheim, J.; Greulich, T.; Nowinski, A.; Barta, I.; Stendardo, M.; et al. Lung Microbiome Composition and Bronchial Epithelial Gene Expression in Patients with COPD versus Healthy Individuals: A Bacterial 16S RRNA Gene Sequencing and Host Transcriptomic Analysis. Lancet Microbe 2021, 2, e300–e310. [Google Scholar] [CrossRef] [PubMed]
  66. Garcia, M.; Hirata, M.; Hossain, N.; Kunz, A. Prevotella Melaninogenica Meningitis and Abscess: Normal Flora in an Abnormal Location. Pediatrics 2018, 142, 602. [Google Scholar] [CrossRef]
  67. Lertpiriyapong, K.; Whary, M.T.; Muthupalani, S.; Lofgren, J.L.; Gamazon, E.R.; Feng, Y.; Ge, Z.; Wang, T.C.; Fox, J.G. Gastric Colonisation with a Restricted Commensal Microbiota Replicates the Promotion of Neoplastic Lesions by Diverse Intestinal Microbiota in the Helicobacter Pylori INS-GAS Mouse Model of Gastric Carcinogenesis. Gut 2014, 63, 54–63. [Google Scholar] [CrossRef] [PubMed]
  68. Dias-Jácome, E.; Libânio, D.; Borges-Canha, M.; Galaghar, A.; Pimentel-Nunes, P. Gastric Microbiota and Carcinogenesis: The Role of Non-Helicobacter Pylori Bacteria—A Systematic Review. Rev. Esp. Enfermedades Dig. 2016, 108, 530–540. [Google Scholar] [CrossRef] [PubMed]
  69. Dahmus, J.D.; Kotler, D.L.; Kastenberg, D.M.; Kistler, C.A. The Gut Microbiome and Colorectal Cancer: A Review of Bacterial Pathogenesis. J. Gastrointest. Oncol. 2018, 9, 769–777. [Google Scholar] [CrossRef] [PubMed]
  70. Newsome, R.C.; Gauthier, J.; Hernandez, M.C.; Abraham, G.E.; Robinson, T.O.; Williams, H.B.; Sloan, M.; Owings, A.; Laird, H.; Christian, T.; et al. The Gut Microbiome of COVID-19 Recovered Patients Returns to Uninfected Status in a Minority-Dominated United States Cohort. Gut Microbes 2021, 13, 1926840. [Google Scholar] [CrossRef] [PubMed]
  71. Newman, A.; Lambert, J.R. Campylobacter Jejuni Causing Flare-Up in Inflammatory Bowel Disease. Lancet 1980, 316, 919. [Google Scholar] [CrossRef] [PubMed]
  72. Starr, J. Clostridium Difficile Associated Diarrhoea: Diagnosis and Treatment. Bmj 2005, 331, 498. [Google Scholar] [CrossRef]
  73. Kahler, C.M. Neisseria Species and Their Complicated Relationships with Human Health. Microbiol. Aust. 2021, 42, 79–83. [Google Scholar] [CrossRef]
  74. Dewhirst, F.E.; Chen, T.; Izard, J.; Paster, B.J.; Tanner, A.C.R.; Yu, W.H.; Lakshmanan, A.; Wade, W.G. The Human Oral Microbiome. J. Bacteriol. 2010, 192, 5002–5017. [Google Scholar] [CrossRef] [PubMed]
  75. Signat, B.; Roques, C.; Poulet, P.; Duffaut, D. Role of Fusobacterium Nucleatum in Periodontal Health and Disease. Curr. Issues Mol. Biol. 2011, 13, 25–36. [Google Scholar]
  76. Knoop, F.C. Haemophilus Infections. In xPharm: The Comprehensive Pharmacology Reference; Elsevier Inc.: Oxford, UK, 2011; pp. 1–4. [Google Scholar] [CrossRef]
  77. Quévrain, E.; Maubert, M.A.; Michon, C.; Chain, F.; Marquant, R.; Tailhades, J.; Miquel, S.; Carlier, L.; Bermúdez-Humarán, L.G.; Pigneur, B.; et al. Identification of an Anti-Inflammatory Protein from Faecalibacterium Prausnitzii, a Commensal Bacterium Deficient in Crohn’s Disease. Gut 2016, 65, 415–425. [Google Scholar] [CrossRef] [PubMed]
  78. Jorth, P.; Turner, K.H.; Gumus, P.; Nizam, N.; Buduneli, N.; Whiteley, M. Metatranscriptomics of the Human Oral Microbiome during Health and Disease. MBio 2014, 5, 10–128. [Google Scholar] [CrossRef] [PubMed]
  79. Anand, S.; Kaur, H.; Mande, S.S. Comparative in Silico Analysis of Butyrate Production Pathways in Gut Commensals and Pathogens. Front. Microbiol. 2016, 7, 1945. [Google Scholar] [CrossRef] [PubMed]
  80. Wei, X.; Yang, Z.; Rey, F.E.; Ridaura, V.K.; Davidson, N.O.; Gordon, J.I.; Semenkovich, C.F. Fatty Acid Synthase Modulates Intestinal Barrier Function through Palmitoylation of Mucin 2. Cell Host Microbe 2012, 11, 140–152. [Google Scholar] [CrossRef] [PubMed]
  81. Benoit, B.; Bruno, J.; Kayal, F.; Estienne, M.; Debard, C.; Ducroc, R.; Plaisancié, P. Saturated and Unsaturated Fatty Acids Differently Modulate Colonic Goblet Cells In Vitro and in Rat Pups. J. Nutr. 2015, 145, 1754–1762. [Google Scholar] [CrossRef] [PubMed]
  82. Zain, M.; Awan, F.R.; Najam, S.S.; Islam, M.; Khan, A.R.; Bilal, A.; Bellili, N.; Marre, M.; Roussel, R.; Fumeron, F. Association of ACACB Gene Polymorphism (Rs2268388, G > A) with Type 2 Diabetes and End Stage Renal Disease in Pakistani Punjabi Population. Meta Gene 2017, 12, 109–112. [Google Scholar] [CrossRef]
  83. Wang, Y.F.; Wang, G.; Li, J.L.; Qu, Y.X.; Liang, X.Y.; Chen, X.D.; Sima, Y.H.; Xu, S.Q. Influence of Hyperproteinemia on Insect Innate Immune Function of the Circulatory System in Bombyx mori. Biology 2021, 10, 112. [Google Scholar] [CrossRef]
  84. Ali, A.M.; Kunugi, H. Hypoproteinemia Predicts Disease Severity and Mortality in COVID-19: A Call for Action. Diagn. Pathol. 2021, 16, 31. [Google Scholar] [CrossRef]
  85. Pirzada, R.H.; Haseeb, M.; Batool, M.; Kim, M.S.; Choi, S. Remdesivir and Ledipasvir among the FDA-Approved Antiviral Drugs Have Potential to Inhibit SARS-CoV-2 Replication. Cells 2021, 10, 1052. [Google Scholar] [CrossRef] [PubMed]
  86. Gabra, H.; BerGenBio ASA. NCT04890509 A Study of Bemcentinib for the Treatment of COVID-19 in Hospitalised Patients. 2021. Available online: https://clinicaltrials.gov/show/NCT04890509 (accessed on 14 October 2021).
  87. Bonaventura, A.; Montecucco, F. Sofosbuvir/Velpatasvir: A Promising Combination. World J. Hepatol. 2016, 8, 785–789. [Google Scholar] [CrossRef] [PubMed]
  88. Hill, R.L.; Singh, I.N.; Brelsfoard, J.; Hall, E.D. Pharmacological Inhibition of Lipid Peroxidative Damage by the 21-Aminosteroid U-74389G Improves Cortical Mitochondrial Function Following Traumatic Brain Injury in Young Adult Male Rats. Neuropharmacology 2020, 170, 108023. [Google Scholar] [CrossRef] [PubMed]
  89. Connor, K.M.; Shapiro, R.E.; Diener, H.C.; Lucas, S.; Kost, J.; Fan, X.; Fei, K.; Assaid, C.; Lines, C.; Ho, T.W. Randomized, Controlled Trial of Telcagepant for the Acute Treatment of Migraine. Neurology 2009, 73, 970–977. [Google Scholar] [CrossRef] [PubMed]
  90. Liu, C.H.; Peng, C.Y.; Fang, Y.J.; Kao, W.Y.; Yang, S.S.; Lin, C.K.; Lai, H.C.; Su, W.P.; Fang, S.U.; Chang, C.C.; et al. Elbasvir/Grazoprevir for Hepatitis C Virus Genotype 1b East-Asian Patients Receiving Hemodialysis. Sci. Rep. 2020, 10, 9180. [Google Scholar] [CrossRef]
  91. Fang, C.; Mei, J.; Tian, H.; Liou, Y.L.; Rong, D.; Zhang, W.; Liao, Q.; Wu, N. CSF3 Is a Potential Drug Target for the Treatment of COVID-19. Front. Physiol. 2021, 11, 605792. [Google Scholar] [CrossRef] [PubMed]
  92. Tutusaus, A.; de Gregorio, E.; Cucarull, B.; Cristóbal, H.; Aresté, C.; Graupera, I.; Coll, M.; Colell, A.; Gausdal, G.; Lorens, J.B.; et al. A Functional Role of GAS6/TAM in Nonalcoholic Steatohepatitis Progression Implicates AXL as Therapeutic Target. Cmgh 2020, 9, 349–368. [Google Scholar] [CrossRef] [PubMed]
  93. Bohan, D.; Van Ert, H.; Ruggio, N.; Rogers, K.J.; Badreddine, M.; Briseno, J.A.A.; Elliff, J.M.; Chavez, R.A.R.; Gao, B.; Stokowy, T.; et al. Phosphatidylserine Receptors Enhance SARS-CoV-2 Infection. PLoS Pathog. 2021, 17, e1009743. [Google Scholar] [CrossRef]
  94. Bohan, D.; Van Ert, H.; Dhakal, S.; Ruggio, N.; Rogers, K.; Stokowy, T.; Gausdal, G.; Jackson, A.; Gabra, H.; Micklem, D.R.; et al. Targeting the Receptor Axl by Bemcentinib Prevents SARS-CoV-2 Infection. Top. Antivir. Med. 2021, 29, 137–138. [Google Scholar]
  95. Elgohary, M.A.S.; Hasan, E.M.; Ibrahim, A.A.; Abdelsalam, M.F.A.; Abdel-Rahman, R.Z.; Zaki, A.I.; Elaatar, M.B.; Elnagar, M.T.; Emam, M.E.; Hamada, M.M.; et al. Efficacy of Sofosbuvir plus Ledipasvir in Egyptian Patients with COVID-19 Compared to Standard Treatment: A Randomized Controlled Trial. J. Med. Life 2022, 15, 350–358. [Google Scholar] [CrossRef]
  96. Zuccaro, V.; Lombardi, A.; Asperges, E.; Sacchi, P.; Bruno, R. PK/PD and Antiviral Activity of Anti-HCV Therapy: Is There Still a Role in the Choice of Treatment? Expert Opin. Drug Metab. Toxicol. 2020, 16, 97–101. [Google Scholar] [CrossRef] [PubMed]
  97. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and Chloroquine Effectively Inhibit the Recently Emerged Novel Coronavirus (2019-NCoV) In Vitro. Cell Res. 2020, 30, 269–271. [Google Scholar] [CrossRef] [PubMed]
  98. Sayad, B.; Sobhani, M.; Khodarahmi, R. Sofosbuvir as Repurposed Antiviral Drug Against COVID-19: Why Were We Convinced to Evaluate the Drug in a Registered/Approved Clinical Trial? Arch. Med. Res. 2020, 51, 577–581. [Google Scholar] [CrossRef] [PubMed]
  99. Al-Karmalawy, A.A.; Alnajjar, R.; Dahabd, M.M.; Metwaly, A.M.; Eissa, I.H. Molecular Docking and Dynamics Simulations Reveal the Potential of Anti-HCV Drugs to Inhibit COVID-19 Main Protease. Pharm. Sci. 2021, 27, S109–S121. [Google Scholar] [CrossRef]
  100. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  101. Gaspar, J.M. NGmerge: Merging Paired-End Reads via Novel Empirically-Derived Models of Sequencing Errors. BMC Bioinformatics 2018, 19, 536. [Google Scholar] [CrossRef] [PubMed]
  102. Ning, Y.; Yang, G.; Chen, Y.; Zhao, X.; Qian, H.; Liu, Y.; Chen, S.; Shi, G. Characteristics of the Urinary Microbiome from Patients with Gout: A Prospective Study. Front. Endocrinol. 2020, 11, 272. [Google Scholar] [CrossRef] [PubMed]
  103. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed]
  104. Roswell, M.; Dushoff, J.; Winfree, R. A Conceptual Guide to Measuring Species Diversity. Oikos 2021, 130, 321–338. [Google Scholar] [CrossRef]
  105. Cox, K.D.; Black, M.J.; Filip, N.; Miller, M.R.; Mohns, K.; Mortimor, J.; Freitas, T.R.; Greiter Loerzer, R.; Gerwing, T.G.; Juanes, F.; et al. Community Assessment Techniques and the Implications for Rarefaction and Extrapolation with Hill Numbers. Ecol. Evol. 2017, 7, 11213–11226. [Google Scholar] [CrossRef]
  106. Legendre, P.; Borcard, D.; Peres-Neto, P.R. Analyzing Beta Diversity: Partitioning the Spatial Variation of Community Composition Data. Ecol. Monogr. 2005, 75, 435–450. [Google Scholar] [CrossRef]
  107. McCafferty, J.; Mühlbauer, M.; Gharaibeh, R.Z.; Arthur, J.C.; Perez-Chanona, E.; Sha, W.; Jobin, C.; Fodor, A.A. Stochastic Changes over Time and Not Founder Effects Drive Cage Effects in Microbial Community Assembly in a Mouse Model. ISME J. 2013, 7, 2116–2125. [Google Scholar] [CrossRef] [PubMed]
  108. Chen, J.; Bittinger, K.; Charlson, E.S.; Hoffmann, C.; Lewis, J.; Wu, G.D.; Collman, R.G.; Bushman, F.D.; Li, H. Associating Microbiome Composition with Environmental Covariates Using Generalized UniFrac Distances. Bioinformatics 2012, 28, 2106–2113. [Google Scholar] [CrossRef] [PubMed]
  109. Ahmad, A.; Yang, W.; Chen, G.; Shafiq, M.; Javed, S.; Zaidi, S.S.A.; Shahid, R.; Liu, C.; Bokhari, H. Analysis of Gut Microbiota of Obese Individuals with Type 2 Diabetes and Healthy Individuals. PLoS ONE 2019, 14, e0226372. [Google Scholar] [CrossRef] [PubMed]
  110. Bajaj, J.S.; Sikaroodi, M.; Gillevet, P.M.; Bajaj, J.S.; Thacker, L.R.; Fagan, A.; White, M.B.; Gavis, E.A.; Hylemon, P.B.; Brown, R.; et al. Gut Microbial RNA and DNA Analysis Predicts Hospitalizations in Cirrhosis. JCI Insight 2018, 3, e98019. [Google Scholar] [CrossRef] [PubMed]
  111. Zhang, X.; Guo, B.; Yi, N. Zero-Inflated Gaussian Mixed Models for Analyzing Longitudinal Microbiome Data. PLoS ONE 2020, 15, e0242073. [Google Scholar] [CrossRef] [PubMed]
  112. Douglas, G.M.; Maffei, V.J.; Zaneveld, J.R.; Yurgel, S.N.; Brown, J.R.; Taylor, C.M.; Huttenhower, C.; Langille, M.G. PICRUSt2 for Prediction of Metagenome Functions. Nat. Biotechnol. 2020, 38, 669–673. [Google Scholar] [CrossRef] [PubMed]
  113. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical Analysis of Taxonomic and Functional Profiles. Bioinformatics 2014, 30, 3123–3124. [Google Scholar] [CrossRef]
  114. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P.; et al. The STRING Database in 2017: Quality-Controlled Protein-Protein Association Networks, Made Broadly Accessible. Nucleic Acids Res. 2017, 45, D362–D368. [Google Scholar] [CrossRef]
  115. Killcoyne, S.; Carter, G.W.; Smith, J.; Boyle, J. Cytoscape: A Community-Based Framework for Network Modeling. Methods Mol. Biol. 2009, 563, 219–239. [Google Scholar] [CrossRef]
  116. Di Muzio, E.; Toti, D.; Polticelli, F. DockingApp: A User Friendly Interface for Facilitated Docking Simulations with AutoDock Vina. J. Comput. Aided. Mol. Des. 2017, 31, 213–218. [Google Scholar] [CrossRef] [PubMed]
  117. Studio, D. Dassault Systemes BIOVIA, Discovery Studio Modelling Environment, Release 4.5; Accelrys Software Inc.: San Diego, CA, USA, 2015.
  118. DeLano, W.L. The PyMOL Molecular Graphics System, Version 2.3; Schrödinger LLC: New York, NY, USA, 2020.
  119. Kaplan, W.; Littlejohn, T.G. Swiss-PDB Viewer (Deep View). Brief. Bioinform. 2001, 2, 195–197. [Google Scholar] [CrossRef]
  120. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeerschd, T.; Zurek, E.; Hutchison, G.R. Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed]
  121. Jaramillo-Botero, A.; Naserifar, S.; Goddard, W.A. General Multiobjective Force Field Optimization Framework, with Application to Reactive Force Fields for Silicon Carbide. J. Chem. Theory Comput. 2014, 10, 1426–1439. [Google Scholar] [CrossRef]
  122. Rostkowski, M.; Olsson, M.H.; Søndergaard, C.R.; Jensen, J.H. Graphical Analysis of PH-Dependent Properties of Proteins Predicted Using PROPKA. BMC Struct. Biol. 2011, 11, 6. [Google Scholar] [CrossRef] [PubMed]
  123. Cheng, F.; Li, W.; Zhou, Y.; Shen, J.; Wu, Z.; Liu, G.; Lee, P.W.; Tang, Y. AdmetSAR: A Comprehensive Source and Free Tool for Assessment of Chemical ADMET Properties. J. Chem. Inf. Model. 2012, 52, 3099–3105. [Google Scholar] [CrossRef] [PubMed]
  124. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed]
  125. Pires, D.E.V.; Blundell, T.L.; Ascher, D.B. PkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures. J. Med. Chem. 2015, 58, 4066–4072. [Google Scholar] [CrossRef]
  126. Land, H.; Humble, M.S. YASARA: A Tool to Obtain Structural Guidance in Biocatalytic Investigations. Methods Mol. Biol. 2018, 1685, 43–67. [Google Scholar] [CrossRef]
  127. Case, D.A.; Cheatham, T.E.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber Biomolecular Simulation Programs. J. Comput. Chem. 2005, 26, 1668. [Google Scholar] [CrossRef]
  128. Dash, R.; Ali, M.C.; Dash, N.; Azad, M.A.K.; Zahid Hosen, S.M.; Hannan, M.A.; Moon, I.S. Structural and Dynamic Characterizations Highlight the Deleterious Role of SULT1A1 R213H Polymorphism in Substrate Binding. Int. J. Mol. Sci. 2019, 20, 6256. [Google Scholar] [CrossRef] [PubMed]
  129. Srinivasan, E.; Rajasekaran, R. Computational Investigation of Curcumin, a Natural Polyphenol That Inhibits the Destabilization and the Aggregation of Human SOD1 Mutant (Ala4Val). RSC Adv. 2016, 6, 102744–102753. [Google Scholar] [CrossRef]
  130. Swargiary, A.; Mahmud, S.; Saleh, M.A. Screening of Phytochemicals as Potent Inhibitor of 3-Chymotrypsin and Papain-like Proteases of SARS-CoV-2: An in Silico Approach to Combat COVID-19. J. Biomol. Struct. Dyn. 2020, 40, 2067–2081. [Google Scholar] [CrossRef] [PubMed]
  131. Mahmud, S.; Parves, M.R.; Riza, Y.M.; Sujon, K.M.; Ray, S.; Tithi, F.A.; Zaoti, Z.F.; Alam, S.; Absar, N. Exploring the Potent Inhibitors and Binding Modes of Phospholipase A2 through in Silico Investigation. J. Biomol. Struct. Dyn. 2020, 38, 4221–4231. [Google Scholar] [CrossRef] [PubMed]
  132. Bappy, S.S.; Sultana, S.; Adhikari, J.; Mahmud, S.; Khan, M.A.; Kibria, K.M.K.; Rahman, M.M.; Shibly, A.Z. Extensive Immunoinformatics Study for the Prediction of Novel Peptide-Based Epitope Vaccine with Docking Confirmation against Envelope Protein of Chikungunya Virus: A Computational Biology Approach. J. Biomol. Struct. Dyn. 2021, 39, 1139–1154. [Google Scholar] [CrossRef] [PubMed]
  133. Chowdhury, K.H.; Chowdhury, M.R.; Mahmud, S.; Tareq, A.M.; Hanif, N.B.; Banu, N.; Ali Reza, A.S.M.; Bin Emran, T.; Simal-Gandara, J. Drug Repurposing Approach against Novel Coronavirus Disease (COVID-19) through Virtual Screening Targeting SARS-CoV-2 Main Protease. Biology 2021, 10, 2. [Google Scholar] [CrossRef]
  134. Khan, M.A.; Mahmud, S.; Alam, A.S.M.R.U.; Rahman, M.E.; Ahmed, F.; Rahmatullah, M. Comparative Molecular Investigation of the Potential Inhibitors against SARS-CoV-2 Main Protease: A Molecular Docking Study. J. Biomol. Struct. Dyn. 2020, 39, 6317–6323. [Google Scholar] [CrossRef]
Figure 1. Violin plot for alpha-diversity results. These results demonstrate the differences in bacterial community between COVID-19 patients and healthy samples based on the (a) Chao1, (b) observed, (c) Shannon, and (d) Simpson index. The Wilcoxon signed-rank test was used to test the significance (p-value < 0.05) difference between patients and healthy individuals.
Figure 1. Violin plot for alpha-diversity results. These results demonstrate the differences in bacterial community between COVID-19 patients and healthy samples based on the (a) Chao1, (b) observed, (c) Shannon, and (d) Simpson index. The Wilcoxon signed-rank test was used to test the significance (p-value < 0.05) difference between patients and healthy individuals.
Pharmaceuticals 17 00432 g001
Figure 2. (a) Bacterial genetic diversity between COVID-19 and control samples was investigated by principal coordinates analysis (PCoA) with the beta diversity indices, where healthy and patient samples were represented by dot and a triangle symbols, respectively. (b) Venn diagrams showed the percentage of shared OTUs between COVID-19 patients and healthy controls.
Figure 2. (a) Bacterial genetic diversity between COVID-19 and control samples was investigated by principal coordinates analysis (PCoA) with the beta diversity indices, where healthy and patient samples were represented by dot and a triangle symbols, respectively. (b) Venn diagrams showed the percentage of shared OTUs between COVID-19 patients and healthy controls.
Pharmaceuticals 17 00432 g002
Figure 3. Composition of bacterial abundances at the phylum and genus levels across COVID-19 and control samples. (a) Histograms of relative abundances of bacteria for the 10 most abundant bacterial phyla across COVID-19 and control samples. (b) Histogram of relative abundance of bacteria for the 12 most abundant bacterial genera across COVID-19 and control samples.
Figure 3. Composition of bacterial abundances at the phylum and genus levels across COVID-19 and control samples. (a) Histograms of relative abundances of bacteria for the 10 most abundant bacterial phyla across COVID-19 and control samples. (b) Histogram of relative abundance of bacteria for the 12 most abundant bacterial genera across COVID-19 and control samples.
Pharmaceuticals 17 00432 g003
Figure 4. Identify differentially abundant bacteria between COVID-19 and control samples at the phylum and genus levels. (a) Histogram of LDA scores to display the significantly differentially abundant phyla between two groups (LDA score > 3.5 and p < 0.05). (b) Histogram of LDA scores to display the significantly differentially abundant genus between two groups (LDA score > 3.5 and p < 0.05). (c,e) The relative bacterial abundance for two groups of two significant phyla levels at Bacteroidetes and Bacteroides was displayed by using a Box plot. (d,f) The relative bacterial abundance of two significant genus levels at Proteobacteria and Streptococcus for two groups was displayed by using a box plot.
Figure 4. Identify differentially abundant bacteria between COVID-19 and control samples at the phylum and genus levels. (a) Histogram of LDA scores to display the significantly differentially abundant phyla between two groups (LDA score > 3.5 and p < 0.05). (b) Histogram of LDA scores to display the significantly differentially abundant genus between two groups (LDA score > 3.5 and p < 0.05). (c,e) The relative bacterial abundance for two groups of two significant phyla levels at Bacteroidetes and Bacteroides was displayed by using a Box plot. (d,f) The relative bacterial abundance of two significant genus levels at Proteobacteria and Streptococcus for two groups was displayed by using a box plot.
Pharmaceuticals 17 00432 g004
Figure 5. Identifying metabolic pathways and associated functional key genes that were associated with SARS-CoV-2 infection. (a) Prediction of metabolic pathway contents correlated with COVID-19 (patients) and control (healthy) groups. Extended error bar-plot representing the differences in mean proportions for each pair of groups with each pathway. Welch’s t statistic produced q-values for the two-tailed test, which were adjusted by controlling FDR at 0.05 with the Benjamini–Hochberg method. (b) Protein–protein interaction networking to select bKGs, where the orange color indicates bKGs.
Figure 5. Identifying metabolic pathways and associated functional key genes that were associated with SARS-CoV-2 infection. (a) Prediction of metabolic pathway contents correlated with COVID-19 (patients) and control (healthy) groups. Extended error bar-plot representing the differences in mean proportions for each pair of groups with each pathway. Welch’s t statistic produced q-values for the two-tailed test, which were adjusted by controlling FDR at 0.05 with the Benjamini–Hochberg method. (b) Protein–protein interaction networking to select bKGs, where the orange color indicates bKGs.
Pharmaceuticals 17 00432 g005
Figure 6. The protein–ligand binding affinity score matrices were displayed, where the X-axis represents the top-ranked 30 drug molecules and the y-axis represents the corresponding ordered bKGs. (a) The image of the score matrix indicates the binding affinities between the proposed 10 bKGs/ receptors and the 30 top-ranked drug molecules. The top 10 drug molecules, determined by their binding affinity scores, are highlighted in purple, and (b) the score matrix indicates the binding affinities between the COVID-19-causing top-ranked 9 published key receptors and the 30 top-ranked drug molecules with 8 drugs are being common with the proposed top 10 drug indicated by the purple color.
Figure 6. The protein–ligand binding affinity score matrices were displayed, where the X-axis represents the top-ranked 30 drug molecules and the y-axis represents the corresponding ordered bKGs. (a) The image of the score matrix indicates the binding affinities between the proposed 10 bKGs/ receptors and the 30 top-ranked drug molecules. The top 10 drug molecules, determined by their binding affinity scores, are highlighted in purple, and (b) the score matrix indicates the binding affinities between the COVID-19-causing top-ranked 9 published key receptors and the 30 top-ranked drug molecules with 8 drugs are being common with the proposed top 10 drug indicated by the purple color.
Pharmaceuticals 17 00432 g006
Figure 7. Visualization of top-scorer target-ligand complexes, with lightseagreen color representing the protein and gold indicating the ligand. The interacting complex (3D) surface view, pose view, and interactions were shown in the 1st, 2nd, and 3rd columns, respectively. (a) Overall top-ranked complex between Bemcentinib and hldD. (b) Overall 2nd top-ranked complex between Ledipasvir and mlaA. (c) Overall 3rd top-ranked complex Velpatasvir and lptD. (d) Complex between lptD and Molnupiravir, which is the highest scorer complex with Molnupiravir only. (e) Complex between mlaA and the Nirmatrelvir drug, which is the highest scorer complex with Nirmatrelvir only. Here, Ledipasvir, Molnupiravir, and Nirmatrelvir are three FDA-approved COVID-19 drugs. The last two complexes are displayed here to compare them with the overall top-scoring three proposed complexes (the first three).
Figure 7. Visualization of top-scorer target-ligand complexes, with lightseagreen color representing the protein and gold indicating the ligand. The interacting complex (3D) surface view, pose view, and interactions were shown in the 1st, 2nd, and 3rd columns, respectively. (a) Overall top-ranked complex between Bemcentinib and hldD. (b) Overall 2nd top-ranked complex between Ledipasvir and mlaA. (c) Overall 3rd top-ranked complex Velpatasvir and lptD. (d) Complex between lptD and Molnupiravir, which is the highest scorer complex with Molnupiravir only. (e) Complex between mlaA and the Nirmatrelvir drug, which is the highest scorer complex with Nirmatrelvir only. Here, Ledipasvir, Molnupiravir, and Nirmatrelvir are three FDA-approved COVID-19 drugs. The last two complexes are displayed here to compare them with the overall top-scoring three proposed complexes (the first three).
Pharmaceuticals 17 00432 g007
Figure 8. The molecular dynamic simulation of the proposed three docked complexes. (a) root mean square deviation of the alpha carbon atoms, (b) binding free energy of the complexes, where a more positive value indicates better binding, (c) degree of rigidity and compactness analysis of the complexes, and (d) protein volume with expansion analysis.
Figure 8. The molecular dynamic simulation of the proposed three docked complexes. (a) root mean square deviation of the alpha carbon atoms, (b) binding free energy of the complexes, where a more positive value indicates better binding, (c) degree of rigidity and compactness analysis of the complexes, and (d) protein volume with expansion analysis.
Pharmaceuticals 17 00432 g008
Figure 9. The workflow of the entire study involves different stages: upstream analysis is represented by green color, downstream analysis by burlywood color, and molecular docking and molecular dynamic simulation study by gold color.
Figure 9. The workflow of the entire study involves different stages: upstream analysis is represented by green color, downstream analysis by burlywood color, and molecular docking and molecular dynamic simulation study by gold color.
Pharmaceuticals 17 00432 g009
Table 1. The top 20 differentially abundant gut microbiomes between COVID-19 and healthy patients were taken based on the adjusted p-value and log2FC value. The MetagenomeSeq R package was used to perform the analysis based on the zero-inflated Gaussian method. Among them, nine bacterial genera were upregulated, and the remaining eleven genera were downregulated.
Table 1. The top 20 differentially abundant gut microbiomes between COVID-19 and healthy patients were taken based on the adjusted p-value and log2FC value. The MetagenomeSeq R package was used to perform the analysis based on the zero-inflated Gaussian method. Among them, nine bacterial genera were upregulated, and the remaining eleven genera were downregulated.
PhylumFamilyGenusSpecieslog2FCAdj. p-Value
ProteobacteriaNeisseriaceaeNeisseriaoralis−4.420.000
FirmicutesStreptococcaceaeStreptococcusinfantis−4.320.000
FirmicutesGemellaceaeUnclassifiedUnclassified−4.290.000
FirmicutesStreptococcaceaeStreptococcusUnclassified−4.150.000
FusobacteriaLeptotrichiaceaeLeptotrichiaUnclassified−40.000
ProteobacteriaOxalobacteraceaeCupriavidusUnclassified−3.80.000
FirmicutesVeillonellaceaeVeillonellaparvula−3.690.000
ProteobacteriaPasteurellaceaeUnclassifiedUnclassified−3.630.000
FirmicutesUnclassifiedUnclassifiedUnclassified−3.520.000
ProteobacteriaBurkholderiaceaeLautropiaUnclassified−3.490.000
ProteobacteriaPasteurellaceaeHaemophilusinfluenzae−3.010.000
FirmicutesUnclassifiedUnclassifiedUnclassified3.080.000
BacteroidetesPorphyromonadaceaeParabacteroidesgordonii3.180.000
ActinobacteriaActinomycetaceaeActinomyceshyovaginalis3.240.000
ProteobacteriaCampylobacteraceaeCampylobacterfetus3.260.000
FirmicutesLachnospiraceaeClostridiumdifficile3.490.000
BacteroidetesPrevotellaceaePrevotellamelaninogenica3.620.000
ActinobacteriaCoriobacteriaceaeAtopobiumRimae3.870.000
FirmicutesPeptostreptococcaceaePeptostreptococcusAnaerobius3.950.000
BacteroidetesBacteroidaceaeBacteroidesAcidifaciens5.10.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kibria, M.K.; Ali, M.A.; Yaseen, M.; Khan, I.A.; Bhat, M.A.; Islam, M.A.; Mahumud, R.A.; Mollah, M.N.H. Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications. Pharmaceuticals 2024, 17, 432. https://doi.org/10.3390/ph17040432

AMA Style

Kibria MK, Ali MA, Yaseen M, Khan IA, Bhat MA, Islam MA, Mahumud RA, Mollah MNH. Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications. Pharmaceuticals. 2024; 17(4):432. https://doi.org/10.3390/ph17040432

Chicago/Turabian Style

Kibria, Md. Kaderi, Md. Ahad Ali, Muhammad Yaseen, Imran Ahmad Khan, Mashooq Ahmad Bhat, Md. Ariful Islam, Rashidul Alam Mahumud, and Md. Nurul Haque Mollah. 2024. "Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications" Pharmaceuticals 17, no. 4: 432. https://doi.org/10.3390/ph17040432

APA Style

Kibria, M. K., Ali, M. A., Yaseen, M., Khan, I. A., Bhat, M. A., Islam, M. A., Mahumud, R. A., & Mollah, M. N. H. (2024). Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications. Pharmaceuticals, 17(4), 432. https://doi.org/10.3390/ph17040432

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop