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) “
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.
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.
Not applicable.
Not applicable.
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.
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.
The authors declare no conflict of interest.
Footnotes
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.
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 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 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 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 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 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 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 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 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.
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.
Phylum | Family | Genus | Species | log2FC | Adj. p-Value |
---|---|---|---|---|---|
Proteobacteria | Neisseriaceae | Neisseria | oralis | −4.42 | 0.000 |
Firmicutes | Streptococcaceae | Streptococcus | infantis | −4.32 | 0.000 |
Firmicutes | Gemellaceae | Unclassified | Unclassified | −4.29 | 0.000 |
Firmicutes | Streptococcaceae | Streptococcus | Unclassified | −4.15 | 0.000 |
Fusobacteria | Leptotrichiaceae | Leptotrichia | Unclassified | −4 | 0.000 |
Proteobacteria | Oxalobacteraceae | Cupriavidus | Unclassified | −3.8 | 0.000 |
Firmicutes | Veillonellaceae | Veillonella | parvula | −3.69 | 0.000 |
Proteobacteria | Pasteurellaceae | Unclassified | Unclassified | −3.63 | 0.000 |
Firmicutes | Unclassified | Unclassified | Unclassified | −3.52 | 0.000 |
Proteobacteria | Burkholderiaceae | Lautropia | Unclassified | −3.49 | 0.000 |
Proteobacteria | Pasteurellaceae | Haemophilus | influenzae | −3.01 | 0.000 |
Firmicutes | Unclassified | Unclassified | Unclassified | 3.08 | 0.000 |
Bacteroidetes | Porphyromonadaceae | Parabacteroides | gordonii | 3.18 | 0.000 |
Actinobacteria | Actinomycetaceae | Actinomyces | hyovaginalis | 3.24 | 0.000 |
Proteobacteria | Campylobacteraceae | Campylobacter | fetus | 3.26 | 0.000 |
Firmicutes | Lachnospiraceae | Clostridium | difficile | 3.49 | 0.000 |
Bacteroidetes | Prevotellaceae | Prevotella | melaninogenica | 3.62 | 0.000 |
Actinobacteria | Coriobacteriaceae | Atopobium | Rimae | 3.87 | 0.000 |
Firmicutes | Peptostreptococcaceae | Peptostreptococcus | Anaerobius | 3.95 | 0.000 |
Bacteroidetes | Bacteroidaceae | Bacteroides | Acidifaciens | 5.1 | 0.000 |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. World Health Organization. COVID-19 Epidemiological Update—15 March 2024; Epidemiology WHO: Geneva, Switzerland, 2024; pp. 6-8.
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, pp. 361-379. [DOI: https://dx.doi.org/10.1038/s41579-023-00878-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37020110]
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, pp. 4225-4233. [DOI: https://dx.doi.org/10.1080/07391102.2020.1775129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32462996]
4. Chen, Y.; Liu, Q.; Guo, D. Emerging Coronaviruses: Genome Structure, Replication, and Pathogenesis. J. Med. Virol.; 2020; 92, pp. 418-423. [DOI: https://dx.doi.org/10.1002/jmv.25681] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31967327]
5. Mengist, H.M.; Dilnessa, T.; Jin, T. Structural Basis of Potential Inhibitors Targeting SARS-CoV-2 Main Protease. Front. Chem.; 2021; 9, 622898. [DOI: https://dx.doi.org/10.3389/fchem.2021.622898] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33889562]
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, pp. 2866-2879. [DOI: https://dx.doi.org/10.1021/acs.jmedchem.1c00566] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34570513]
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, pp. 325-337. [DOI: https://dx.doi.org/10.1080/13543776.2021.1880568] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33475441]
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, pp. 6039-6051. [DOI: https://dx.doi.org/10.1080/07391102.2021.1875886] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33491569]
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, pp. 784-790. [DOI: https://dx.doi.org/10.1016/j.csbj.2020.03.025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32280433]
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. [DOI: https://dx.doi.org/10.1186/s12941-023-00625-8]
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. [DOI: https://dx.doi.org/10.3390/pathogens12010117]
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, pp. 985-992. [DOI: https://dx.doi.org/10.1016/j.jmii.2022.09.006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36243668]
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, pp. 1201-1208. [DOI: https://dx.doi.org/10.1177/08850666211032175]
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, pp. 1622-1629. [DOI: https://dx.doi.org/10.1016/j.cmi.2020.07.016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32711058]
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. [DOI: https://dx.doi.org/10.3389/fmolb.2023.1113969] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36994427]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0254671]
17. Feldman, C.; Anderson, R. The Role of Co-Infections and Secondary Infections in Patients with COVID-19. Pneumonia; 2021; 13, 5. [DOI: https://dx.doi.org/10.1186/s41479-021-00083-w]
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, pp. 591-605. [DOI: https://dx.doi.org/10.1007/s15010-021-01602-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33709380]
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. [DOI: https://dx.doi.org/10.3389/fpubh.2021.712190] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34513787]
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, pp. 533-544. [DOI: https://dx.doi.org/10.1016/S2213-2600(21)00125-9]
21. Mahase, E. COVID-19: Sore Throat, Fatigue, and Myalgia Are More Common with New UK Variant. BMJ; 2021; 372, n288. [DOI: https://dx.doi.org/10.1136/bmj.n288]
22. Zhang, Z. Diagnosis and Treatment Protocol for COVID-19 Patients (Trial Version 9). Health Care Sci.; 2022; 1, pp. 14-28. [DOI: https://dx.doi.org/10.1002/hcs2.1]
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. [DOI: https://dx.doi.org/10.3389/fpubh.2022.945988] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36249235]
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, pp. 1602-1608. [DOI: https://dx.doi.org/10.1590/1806-9282.66.11.1602] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33295417]
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. [DOI: https://dx.doi.org/10.1186/s41983-020-00258-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33390724]
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. [DOI: https://dx.doi.org/10.1016/j.biopha.2022.112757] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35231696]
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. [DOI: https://dx.doi.org/10.3390/medicina58101500] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36295660]
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, pp. 1663-1672. [DOI: https://dx.doi.org/10.1016/j.cgh.2020.04.001]
29. Hurn, E.; Dickinson, L.; Abraham, J.A. Bacterial Meningitis and COVID-19: A Complex Patient Journey. BMJ Case Rep.; 2021; 14, e239533. [DOI: https://dx.doi.org/10.1136/bcr-2020-239533]
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. [DOI: https://dx.doi.org/10.1186/s13256-022-03386-8]
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. [DOI: https://dx.doi.org/10.1136/bmj.n436]
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. [DOI: https://dx.doi.org/10.1186/s12879-021-07007-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34961472]
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. [DOI: https://dx.doi.org/10.1080/19490976.2022.2027853] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35129072]
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. [DOI: https://dx.doi.org/10.1093/ijpp/riac019.028]
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. [DOI: https://dx.doi.org/10.1038/s41598-022-08073-8]
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. [DOI: https://dx.doi.org/10.3390/vaccines10081248] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36016137]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0266124]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0281981]
39. Weigert Muñoz, A.; Zhao, W.; Sieber, S.A. Monitoring Host-Pathogen Interactions Using Chemical Proteomics. RSC Chem. Biol.; 2023; 5, pp. 73-89. [DOI: https://dx.doi.org/10.1039/D3CB00135K]
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. [DOI: https://dx.doi.org/10.3389/fbinf.2021.717141]
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. [DOI: https://dx.doi.org/10.1038/ncomms14092]
42. Balint, D.; Brito, I.L. Human–Gut Bacterial Protein–Protein Interactions: Understudied but Impactful to Human Health. Trends Microbiol.; 2023; [DOI: https://dx.doi.org/10.1016/j.tim.2023.09.009]
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. [DOI: https://dx.doi.org/10.1038/s42003-021-01796-w]
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. [DOI: https://dx.doi.org/10.3389/fmed.2022.821777] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35237631]
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, pp. 1324-1337. [DOI: https://dx.doi.org/10.1093/bib/bbaa376] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33333559]
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, pp. 941-948. [DOI: https://dx.doi.org/10.1007/s12275-021-1206-5]
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. [DOI: https://dx.doi.org/10.3389/fcimb.2021.670424]
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, pp. 6443-6455. [DOI: https://dx.doi.org/10.1038/s41388-021-02023-w]
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. [DOI: https://dx.doi.org/10.1186/s12903-022-02093-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35241064]
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. [DOI: https://dx.doi.org/10.3390/ijerph191610107]
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. [DOI: https://dx.doi.org/10.3389/fcimb.2022.824578]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0247041]
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, pp. 944-955.e948. [DOI: https://dx.doi.org/10.1053/j.gastro.2020.05.048] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32442562]
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. [DOI: https://dx.doi.org/10.3389/fmicb.2021.761067] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34803986]
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. [DOI: https://dx.doi.org/10.1038/s41522-021-00232-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34294722]
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, pp. 97-118. [DOI: https://dx.doi.org/10.5772/intechopen.93635]
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, pp. 320-331. [DOI: https://dx.doi.org/10.1016/j.jiph.2022.12.022] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36657243]
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, pp. 840-862. [DOI: https://dx.doi.org/10.1128/CMR.17.4.840-862.2004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15489351]
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, pp. 2769-2773.
60. Syed, Y.Y. Molnupiravir: First Approval. Drugs; 2022; 82, pp. 455-460. [DOI: https://dx.doi.org/10.1007/s40265-022-01684-5]
61. Lamb, Y.N. Nirmatrelvir Plus Ritonavir: First Approval. Drugs; 2022; 82, pp. 585-591. [DOI: https://dx.doi.org/10.1007/s40265-022-01692-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35305258]
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. [DOI: https://dx.doi.org/10.3390/nu12051474]
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, pp. 797-803. [DOI: https://dx.doi.org/10.1086/591130]
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. [DOI: https://dx.doi.org/10.1038/s41598-018-35877-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30518941]
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, pp. e300-e310. [DOI: https://dx.doi.org/10.1016/S2666-5247(21)00035-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35544166]
66. Garcia, M.; Hirata, M.; Hossain, N.; Kunz, A. Prevotella Melaninogenica Meningitis and Abscess: Normal Flora in an Abnormal Location. Pediatrics; 2018; 142, 602. [DOI: https://dx.doi.org/10.1542/peds.142.1MA7.602]
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, pp. 54-63. [DOI: https://dx.doi.org/10.1136/gutjnl-2013-305178] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23812323]
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, pp. 530-540. [DOI: https://dx.doi.org/10.17235/reed.2016.4261/2016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27604361]
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, pp. 769-777. [DOI: https://dx.doi.org/10.21037/jgo.2018.04.07] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30151274]
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. [DOI: https://dx.doi.org/10.1080/19490976.2021.1926840] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34100340]
71. Newman, A.; Lambert, J.R. Campylobacter Jejuni Causing Flare-Up in Inflammatory Bowel Disease. Lancet; 1980; 316, 919. [DOI: https://dx.doi.org/10.1016/S0140-6736(80)92078-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/6107569]
72. Starr, J. Clostridium Difficile Associated Diarrhoea: Diagnosis and Treatment. Bmj; 2005; 331, 498. [DOI: https://dx.doi.org/10.1136/bmj.331.7515.498]
73. Kahler, C.M. Neisseria Species and Their Complicated Relationships with Human Health. Microbiol. Aust.; 2021; 42, pp. 79-83. [DOI: https://dx.doi.org/10.1071/MA21024]
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, pp. 5002-5017. [DOI: https://dx.doi.org/10.1128/JB.00542-10] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20656903]
75. Signat, B.; Roques, C.; Poulet, P.; Duffaut, D. Role of Fusobacterium Nucleatum in Periodontal Health and Disease. Curr. Issues Mol. Biol.; 2011; 13, pp. 25-36.
76. Knoop, F.C. Haemophilus Infections. xPharm: The Comprehensive Pharmacology Reference; Elsevier Inc.: Oxford, UK, 2011; pp. 1-4. [DOI: https://dx.doi.org/10.1016/B978-008055232-3.60874-8]
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, pp. 415-425. [DOI: https://dx.doi.org/10.1136/gutjnl-2014-307649] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26045134]
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, pp. 10-128. [DOI: https://dx.doi.org/10.1128/mBio.01012-14] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24692635]
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. [DOI: https://dx.doi.org/10.3389/fmicb.2016.01945] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27994578]
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, pp. 140-152. [DOI: https://dx.doi.org/10.1016/j.chom.2011.12.006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22341463]
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, pp. 1754-1762. [DOI: https://dx.doi.org/10.3945/jn.115.211441] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26108543]
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, pp. 109-112. [DOI: https://dx.doi.org/10.1016/j.mgene.2017.02.003]
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. [DOI: https://dx.doi.org/10.3390/biology10020112]
84. Ali, A.M.; Kunugi, H. Hypoproteinemia Predicts Disease Severity and Mortality in COVID-19: A Call for Action. Diagn. Pathol.; 2021; 16, 31. [DOI: https://dx.doi.org/10.1186/s13000-021-01092-5]
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. [DOI: https://dx.doi.org/10.3390/cells10051052] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33946869]
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, pp. 785-789. [DOI: https://dx.doi.org/10.4254/wjh.v8.i19.785] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27429714]
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. [DOI: https://dx.doi.org/10.1016/j.neuropharm.2020.108023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32142792]
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, pp. 970-977. [DOI: https://dx.doi.org/10.1212/WNL.0b013e3181b87942] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19770473]
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. [DOI: https://dx.doi.org/10.1038/s41598-020-66182-8]
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. [DOI: https://dx.doi.org/10.3389/fphys.2020.605792] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33551833]
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, pp. 349-368. [DOI: https://dx.doi.org/10.1016/j.jcmgh.2019.10.010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31689560]
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. [DOI: https://dx.doi.org/10.1371/journal.ppat.1009743]
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, pp. 137-138.
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, pp. 350-358. [DOI: https://dx.doi.org/10.25122/jml-2021-0175]
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, pp. 97-101. [DOI: https://dx.doi.org/10.1080/17425255.2020.1721459] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32003256]
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, pp. 269-271. [DOI: https://dx.doi.org/10.1038/s41422-020-0282-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32020029]
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, pp. 577-581. [DOI: https://dx.doi.org/10.1016/j.arcmed.2020.04.018] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32387040]
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, pp. S109-S121. [DOI: https://dx.doi.org/10.34172/PS.2021.3]
100. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics; 2014; 30, pp. 2114-2120. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu170] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24695404]
101. Gaspar, J.M. NGmerge: Merging Paired-End Reads via Novel Empirically-Derived Models of Sequencing Errors. BMC Bioinformatics; 2018; 19, 536. [DOI: https://dx.doi.org/10.1186/s12859-018-2579-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30572828]
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. [DOI: https://dx.doi.org/10.3389/fendo.2020.00272] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32508748]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0061217] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23630581]
104. Roswell, M.; Dushoff, J.; Winfree, R. A Conceptual Guide to Measuring Species Diversity. Oikos; 2021; 130, pp. 321-338. [DOI: https://dx.doi.org/10.1111/oik.07202]
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, pp. 11213-11226. [DOI: https://dx.doi.org/10.1002/ece3.3580]
106. Legendre, P.; Borcard, D.; Peres-Neto, P.R. Analyzing Beta Diversity: Partitioning the Spatial Variation of Community Composition Data. Ecol. Monogr.; 2005; 75, pp. 435-450. [DOI: https://dx.doi.org/10.1890/05-0549]
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, pp. 2116-2125. [DOI: https://dx.doi.org/10.1038/ismej.2013.106] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23823492]
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, pp. 2106-2113. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts342] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22711789]
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. [DOI: https://dx.doi.org/10.1371/journal.pone.0226372] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31891582]
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. [DOI: https://dx.doi.org/10.1172/jci.insight.98019] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29515036]
111. Zhang, X.; Guo, B.; Yi, N. Zero-Inflated Gaussian Mixed Models for Analyzing Longitudinal Microbiome Data. PLoS ONE; 2020; 15, e0242073. [DOI: https://dx.doi.org/10.1371/journal.pone.0242073] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33166356]
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, pp. 669-673. [DOI: https://dx.doi.org/10.1038/s41587-020-0548-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32483366]
113. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical Analysis of Taxonomic and Functional Profiles. Bioinformatics; 2014; 30, pp. 3123-3124. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu494]
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, pp. D362-D368. [DOI: https://dx.doi.org/10.1093/nar/gkw937]
115. Killcoyne, S.; Carter, G.W.; Smith, J.; Boyle, J. Cytoscape: A Community-Based Framework for Network Modeling. Methods Mol. Biol.; 2009; 563, pp. 219-239. [DOI: https://dx.doi.org/10.1007/978-1-60761-175-2_12]
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, pp. 213-218. [DOI: https://dx.doi.org/10.1007/s10822-016-0006-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28063067]
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, pp. 195-197. [DOI: https://dx.doi.org/10.1093/bib/2.2.195]
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. [DOI: https://dx.doi.org/10.1186/1758-2946-4-17] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22889332]
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, pp. 1426-1439. [DOI: https://dx.doi.org/10.1021/ct5001044]
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. [DOI: https://dx.doi.org/10.1186/1472-6807-11-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21269479]
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, pp. 3099-3105. [DOI: https://dx.doi.org/10.1021/ci300367a] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23092397]
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. [DOI: https://dx.doi.org/10.1038/srep42717] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28256516]
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, pp. 4066-4072. [DOI: https://dx.doi.org/10.1021/acs.jmedchem.5b00104]
126. Land, H.; Humble, M.S. YASARA: A Tool to Obtain Structural Guidance in Biocatalytic Investigations. Methods Mol. Biol.; 2018; 1685, pp. 43-67. [DOI: https://dx.doi.org/10.1007/978-1-4939-7366-8_4]
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. [DOI: https://dx.doi.org/10.1002/jcc.20290]
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. [DOI: https://dx.doi.org/10.3390/ijms20246256] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31835852]
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, pp. 102744-102753. [DOI: https://dx.doi.org/10.1039/C6RA21927F]
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, pp. 2067-2081. [DOI: https://dx.doi.org/10.1080/07391102.2020.1835729] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33089730]
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, pp. 4221-4231. [DOI: https://dx.doi.org/10.1080/07391102.2019.1680440] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31607222]
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, pp. 1139-1154. [DOI: https://dx.doi.org/10.1080/07391102.2020.1726815] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32037968]
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. [DOI: https://dx.doi.org/10.3390/biology10010002]
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, pp. 6317-6323. [DOI: https://dx.doi.org/10.1080/07391102.2020.1796813]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details





1 Bioinformatics Laboratory, Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh;
2 Bioinformatics Laboratory, Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh;
3 Institute of Chemical Sciences, University of Swat, Main Campus, Charbagh 19130, Pakistan;
4 Department of Chemistry, Government College University, Faisalabad 38000, Pakistan;
5 Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh 11421, Saudi Arabia;
6 Bioinformatics Laboratory, Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh;
7 NHMRC Clinical Trials Centre, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia;