Next Article in Journal
Syntheses, Structures, and Corrosion Inhibition of Various Alkali Metal Carboxylate Complexes
Previous Article in Journal
Structural, Spectroscopic, and Dynamic Properties of Li2+(X2∑g+) in Interaction with Krypton Atom
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 14
10.3390/molecules28145514
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans

by
Shakti Chandra Vadhana Marimuthu
1,
Jayaprabhakaran Murugesan
1,
Ewa Babkiewicz
2,3,
Piotr Maszczyk
2,
Murugesan Sankaranarayanan
4,
Esakkimuthu Thangamariappan
1,
Joseph Christina Rosy
1,
Sureshbabu Ram Kumar Pandian
1,
Selvaraj Kunjiappan
1,
Vanavil Balakrishnan
1 and
Krishnan Sundar
1,*
1
Department of Biotechnology, Kalasalingam Academy of Research and Education, Krishnankoil 626126, India
2
Department of Hydrobiology, Faculty of Biology, University of Warsaw, 02-089 Warsaw, Poland
3
Biological and Chemical Research Centre, University of Warsaw, 02-089 Warsaw, Poland
4
Medicinal Chemistry Research Laboratory, Department of Pharmacy, Birla Institute of Technology and Science-Pilani, Pilani 333031, India
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(14), 5514; https://doi.org/10.3390/molecules28145514
Submission received: 10 June 2023 / Revised: 15 July 2023 / Accepted: 17 July 2023 / Published: 19 July 2023
Browse Figures
Versions Notes

Abstract

:
Streptococcus mutans, a gram-positive oral pathogen, is the primary causative agent of dental caries. Biofilm formation, a critical characteristic of S. mutans, is regulated by quorum sensing (QS). This study aimed to utilize pharmacoinformatics techniques to screen and identify effective phytochemicals that can target specific proteins involved in the quorum sensing pathway of S. mutans. A computational approach involving homology modeling, model validation, molecular docking, and molecular dynamics (MD) simulation was employed. The 3D structures of the quorum sensing target proteins, namely SecA, SMU1784c, OppC, YidC2, CiaR, SpaR, and LepC, were modeled using SWISS-MODEL and validated using a Ramachandran plot. Metabolites from Azadirachta indica (Neem), Morinda citrifolia (Noni), and Salvadora persica (Miswak) were docked against these proteins using AutoDockTools. MD simulations were conducted to assess stable interactions between the highest-scoring ligands and the target proteins. Additionally, the ADMET properties of the ligands were evaluated using SwissADME and pkCSM tools. The results demonstrated that campesterol, meliantrol, stigmasterol, isofucosterol, and ursolic acid exhibited the strongest binding affinity for CiaR, LepC, OppC, SpaR, and Yidc2, respectively. Furthermore, citrostadienol showed the highest binding affinity for both SMU1784c and SecA. Notably, specific amino acid residues, including ASP86, ARG182, ILE179, GLU143, ASP237, PRO101, and VAL84 from CiaR, LepC, OppC, SecA, SMU1784c, SpaR, and YidC2, respectively, exhibited significant interactions with their respective ligands. While the docking study indicated favorable binding energies, the MD simulations and ADMET studies underscored the substantial binding affinity and stability of the ligands with the target proteins. However, further in vitro studies are necessary to validate the efficacy of these top hits against S. mutans.

1. Introduction

Streptococcus mutans is a gram-positive bacterium commonly found in the human oral cavity. While it is considered part of the normal microbial flora in the oral cavity, it is also the primary causative agent of dental caries [1,2]. Despite numerous studies that have been reported and implemented in the past to mitigate dental caries, none of them have succeeded in completely eradicating the disease. Additionally, S. mutans is known to be associated with extra-oral pathological conditions such as infective endocarditis, IgA nephropathy, cerebral microbleeds, and atherosclerosis [3,4]. Therefore, there is a pressing need to develop a practical approach for effectively eliminating S. mutans. The Darwinian theory postulates that organisms develop resistance to factors that threaten their survival. This phenomenon is a primary driver behind the development of antibiotic resistance in several pathogens [5]. Hence, it is imperative to develop a treatment that can selectively inhibit the pathogen’s virulence without endangering its overall survival.
One such strategy is controlling the quorum sensing of bacteria. Quorum sensing is a density-dependent bacterial communication mechanism that regulates biofilm formation and the expression of virulent characteristics. Therefore, targeting this mechanism may control their virulence without affecting their survival. Several virulence factors of S. mutans are involved in quorum sensing mechanisms, and targeting these factors may impair cell communication [6,7].
In the sequenced oral pathogen, S. mutans UA159, the CiaRH two-component signal transduction system (TCS) acts as a global regulator for multiple stress responses, including biofilm formation, acid tolerance, bacteriocin production, genetic competence, stress resistance, and pathogenesis [8]. Bacterial oligopeptide transport systems (Opp) play a role in regulating ATP hydrolysis and energy generation required for solute transport [9]. The Sec translocase and yidC are translocase systems that assist in the transport of proteins across the cell membrane [10]. Another virulent protein of S. mutans is encoded by the gene smu1784c, which is involved in stress management and biofilm formation [11].
Regularly used chemical substances such as fluoride, quaternary ammonium salts, and other commonly used antimicrobial agents can lead to undesirable side effects [12,13]. Prolonged use of antimicrobial medications has also contributed to the emergence of antimicrobial resistance, which is now recognized as a significant health problem [14,15]. Consequently, there has been a growing interest in exploring the potential of natural herbal plants as quorum-quenching agents to combat pathogens [16,17,18]. These natural products offer many advantages, such as reduced costs, fewer side effects, and increased efficacy [19]. For example, green tea extract has been found to suppress the growth of Porphyromonas gingivalis, and its potential in preventing and treating periodontitis is being investigated [15]. Furthermore, aqueous extracts of Azadirachta indica (Neem) sticks have been shown to inhibit bacterial aggregation, proliferation, adherence to hydroxyapatite, and synthesis of insoluble glucan, thus altering in vitro plaque development [20]. To expedite the process of identifying molecules with desired antibacterial activity from plant extracts, a time-saving approach is needed. Therefore, computational methods are utilized to accelerate the drug development process. Pharmacoinformatics tools play a crucial role in the rapid discovery of drugs for various diseases, offering a cost-effective approach to drug discovery. Several computer-based tools can be employed to search for potential drug targets in pathogens such as S. mutans and others.
Furthermore, these tools also aid in the study of drug–target interactions. In the current study, drug targets involved in the quorum sensing of S. mutans have been selected based on our previous study [21]. An in silico approach has been undertaken to screen potential phytochemicals from selected plants and to analyze their binding potential with the target proteins.

2. Results

2.1. Homology Modeling and Model Validation of the Target Proteins

SWISS-MODEL performs BLAST for each target protein against PDB and provides a list of templates. In the present study, for all the seven target proteins, model structures were built using SWISS-MODEL since none of them had a crystal structure available. For each of the selected target proteins, five templates were selected based on the identity and sequence coverage, and the 3D models were built.
The modeled structures that were built were assessed and validated using a Ramachandran plot using the SWISS-MODEL structure assessment tool, and one out of five models was selected for further studies (Supplementary Figure S1). The predicted structures of all the target proteins are depicted in Supplementary Figure S2.

2.2. Molecular Docking

Among the 110 ligands tested, 17 compounds were found to bind effectively with all the target proteins (Table 1). Among the 17 compounds that exhibited higher binding towards the selected targets of S. mutans, 15 compounds were from A. indica, and 1 each from M. citrifolia and S. persica. Among the 17 lead compounds, the top hits that exhibited lower binding energy and, hence, high binding affinity against each target were then selected for further simulation analysis (Table 2).
The 3D and 2D interactions between amino acid residues of the target proteins and their respective top hit compounds were also plotted using Discovery Studio. The list of amino acids present in each of the targets that interact with the ligands is presented in Table 3. The 3D and 2D plots depicting the interactions between target proteins and their ligands are presented in Figure 1 and Figure 2, respectively.

2.3. Molecular Dynamics (MD) Simulation

MD simulation assessed the inter-molecular protein–ligand interactions under an artificial environment with specified thermodynamical conditions such as temperature, volume, density, and pressure for 100 ns time duration. Further, a final production step aided in the exploration of the structural modification of the complex. The analysis of unique parameters such as root mean square deviation (RMSD) (Figure 3), simulation interactions diagram (Figure 4), protein–ligand contact (Figure 5), and timeline representation of interaction and contacts (Figure 6) between the target protein and their respective ligands, aided in the analysis of structural changes at each level.
The hydrogen bonding interaction between amino acid residues of the target protein and their respective ligand analyzed in MD simulation was found to correlate with that of AutoDock results (Table 3). The simulation study provided insight into the stability of the interaction between the top hits and their respective target proteins.
In order to perform the post-MM-GBSA (molecular mechanics with generalized Born and surface area solvation) analysis of the free binding energy calculation, 0-2002 frames with a sampling size of roughly 10 steps were produced. During the MM-GBSA calculation of the 100 ns MD data of the target proteins with their respective ligands, a total of 201 frames were processed and analyzed. All the complexes showed good binding affinity, thereby validating the docking and MD results (Table 4).

2.4. ADMET Analysis

Analysis of ADMET properties of the top hits was performed using the SwissADME online tool and the pkCSM tool, and the results are presented in Table 5 and Table 6, respectively. All six compounds had molecular weights greater than 400 g/mol. Except for meliantrol, which was predicted to be moderately soluble in water, all other selected bioactive compounds were predicted to be poorly soluble in water. Solubility of ligands aid in the proper solvation and absorption into the host body. This also aids in the formulation of a solvent or carrier for delivery. Due to their poor solubility, all the hits violated one out of five of Lipinski’s rules. The RADAR plot (Figure 7) shows the overall drug-likeness of the top hit ligands. Meliantrol passes all the criteria for drug-likeness; however, other ligands violate the drug-likeness due to their poor solubility in water. The BOILED-Egg plot (Figure 8) shows all the compounds, except citrostadienol, in the area of intestinal absorption. However, only meliantrol was predicted to exhibit good GI absorption. None of the hits cross the blood–brain barrier. The results were correlated with that of pkCSM.

3. Discussion

Streptococcus mutans is one of the oral commensals that opportunistically cause dental caries. One of the pathogenic mechanisms adopted by S. mutans to cause dental caries is the formation of biofilm on the tooth surface which results in the erosion of tooth enamel through acid production. Biofilm formation requires quorum sensing, a density-dependent communication mechanism. The development of antimicrobial resistance in these pathogens can be avoided by using newer approaches rather than conventional antibiotics. These approaches should control the virulence of the pathogen instead of killing them. Hence, the objective of this study was to develop a compound that can control the quorum sensing in S. mutans rather than killing it. Pharmacoinformatics-based approaches such as drug target identification, molecular docking, and molecular dynamics simulations have accelerated the process of drug discovery. In the present study, we used bioinformatics tools and software to screen phytochemicals from A. indica, M. citrifolia, and S. persica for their potential as anti-quorum sensing agents against the selected target proteins of S. mutans.
The target proteins involved in quorum sensing were selected from our previous study [21]. These seven targets included CiaR (putative response regulator CiaR), LepC (signal peptidase I), OppC (putative transmembrane protein, permease OppC), SecA (protein translocase subunit SecA), SMU1784c (putative Eep protein-like protein), SMU_659 (putative response regulator SpaR), and YidC2 (membrane protein insertase YidC2). The target protein, CiaR, is a response regulator in a two-component signal transport system that controls several virulent characteristics of S. mutans. These virulent characteristics include mutacin I activity, oxidative stress tolerance, acid tolerance, and biofilm formation [22,23]. Studies on S. sanguinis showed the development of a fragile biofilm as a result of the CiaR mutation, which resulted in decreased polysaccharide synthesis [24]. The product of the lepC target gene is a signal peptidase that aids in the export of several virulent proteins. It has also been used as a housekeeping gene for several studies [10,25]. The target protein, OppC, is an oligopeptide permease of the ABC transporter family. It helps the bacteria to regulate XIP (sig X inducing peptide) production and in competence development [9,26]. SecA is a membrane protein translocase that helps the bacteria export proteins, leading to an increase in the virulence of the organism [10,27,28]. Another target, SMU1784c, plays an important role in the management of oxidative and acid stress, EPS production, and biofilm formation [11]. The SpaR protein is a response regulator of the spa (surface protein antigen) family, which is one of the virulence factors of S. mutans [28,29]. The membrane protein insertase of S. mutans, Yidc2, helps in EPS production and biofilm formation [10,30,31].
The resolved structures of the target proteins were not available in PDB or any other structural databases. Hence, the SWISS-MODEL online tool was used to predict the 3D structure of target molecules. Once the protein sequence is submitted as a query, SWISS-MODEL performs BLAST against PDB and gives a list of templates. Based on high identity and sequence coverage, five templates were selected for each protein, and the 3D models were built. A similar approach has been followed by researchers who have used SWISS-MODEL to model NOX 2 of S. mutans by using the crystal structure of NADH oxidase from S. pyogenes as a template [32,33]. Various proteins of S. mutans that were modeled using SWISS-MODEL include domain V of glucosyltransferase (GTF-SI) [34]; SMU.63, an amyloid-like secretory protein of S. mutans [35]; the Spase I protein and β-sheet-rich N-terminal collagen-binding domain (CBD) of Cnm, a collagen- and laminin-binding surface adhesin protein of S. mutans [36,37]. Since the models are predicted in silico, it requires validation before further processing. Hence, the Ramachandran plot from the structure assessment tool of SWISS-MODEL was used for validation of the predicted structures. All five models had more than 90% residues in the allowed region. The models from other templates with less than 90% residues in the allowed region were rejected. The remaining selected models were taken for further docking studies. The modeled structures of fibronectin/fibrinogen binding protein (FBP) from S. mutans, phospholipase D (F13) protein of monkeypox virus, 3-chymotrypsin and papain-like proteases of SARS-CoV2, and U box domain-containing protein gene (GsPUB8) from Glycine soja were all validated using the Ramachandran plot [38,39,40,41].
Computational drug discovery studies evaluate the binding of ligands with the target protein, but the agonist or antagonistic effect of ligands on the target protein is only validated through in vitro and in vivo experiments [42,43]. In the current study, based on an initial screening of 110 compounds using AutoDock, 17 high binders were selected that were found to bind with all the target proteins efficiently. Among the 17 compounds, 15 were from A. indica, and 1 each from M. citrifolia and S. persica. Campesterol, meliantrol, stigmasterol, isofucosterol, and ursolic acid were the top binders specific for the target proteins, CiaR, LepC, OppC, SpaR, and Yidc2, with a binding energy of −8.76, −10.16, −6.75, −9.1, and −9.73 kcal/mol, respectively. Citrostadienol was the high binder against two of the selected targets, SMU1784c and SecA, with a binding score of −8.45 and −9.88 kcal/mol, respectively. The source of ursolic acid is M. citrifolia, whereas all other leads are constituents of A. indica. Molecular docking is a versatile tool that is very useful in screening hundreds of compounds before testing the effective ones using in vitro studies. Molecular docking using AutoDock tools has been used previously to study both agonist and antagonist activity of various natural and synthetic ligands. In silico and in vitro agonistic activity of ligands have been studied for the treatment of diabetes [44,45], Parkinson’s disease [46], and cardiac diseases [47]. The antagonistic activity of ligands against S. mutans [48,49], Leishmania donovani [50], Helicobacter pylori [51], and SARS-CoV-2 [52] has also been studied. All these research works corroborate the necessity of in vitro experiments in the validation of computational analysis. However, they also demonstrate that phytocompounds and synthetic compounds can cause competitive or non-competitive inhibition of target proteins involved in microbial diseases.
The RMSD of the protein–ligand complex is plotted to evaluate the stability of the interaction between the protein and the ligand. The RMSD plot of target protein CiaR in complex with its top binding ligand, campesterol, displayed a fluctuation in RMSD up to 7 Å to 13 Å (Figure 5a). Ligand RMSD was stable, and fluctuations were between 9 and 10.5 Å. The complex of LepC and its top hit ligand meliantrol shows RMSD fluctuation up to 8 and 12 Å (Figure 5b). Ligand RMSD was stable, and fluctuations were between 4 and 10.6 Å. The simulation of target protein OppC in complex with stigmasterol shows the stability of protein at 9 to 16 Å (Figure 5c). Ligand RMSD was stable, and fluctuations were between 9 and 1.6 Å. In the SecA–citrostadienol complex, the protein remains stable at 3 to 4 Å (Figure 5d), whereas the ligand fluctuation is at 14–16 Å and, hence, shows stable interaction between the protein and the ligand. The SMU1784c–citrostadienol complex shows fluctuation up to 2 to 4 Å (Figure 5e), whereas the ligand fluctuation is at 3–6 Å and, hence, shows stable interaction between the protein and the ligand. The Yidc2 and ursolic acid complex shows good RMSD results. The protein fluctuation is up to 6 Å and the ligand up to 9 Å. The RMSD plot converges from 20 ns till the end of the simulation at 100 ns, thus showing a stable interaction between the protein and the ligand. The hydrogen bond and other interactions plotted in the simulation interactions diagram correlate with the docking results. The hydrogen bond interaction between residues ASP86, LEU155, and PRO101 in CiaR, OppC, and SpaR with their respective ligands was observed both in AutoDock and MD simulations. Similarly, in the LepC protein, amino acid residues ILE69 and GLU178 were observed. In SecA, SMU1784c, and Yidc2 though there were no common residues binding through a hydrogen bond, the same residues bind with ligands through other types of bonds. The highly interacting amino acids were ASP86 of CiaR, ARG182 of LepC, ILE179 of OppC, GLU143 of SecA, ASP237 of SMU_1784c, PRO101 of SpaR, and VAL84 of YidC2. Similarly, molecular dynamic simulations and an energy calculation method have been used by researchers to study the LPXTG sequence in the C-terminus of surface proteins, the substrate of the cysteine transpeptidase sortase A (SrtA) enzyme, to better understand how leucine residue affects the dynamics of the enzyme-substrate complex structure. According to the findings, the substrate’s ‘Leu’ residue appears to be essential for anchoring and guiding the conformational shift of the enzyme SrtA [53]. Molecular docking and dynamics simulation studies have been exploited to study the inhibition of glucan sucrase-mediated biofilm formation of S. mutans by thiosemicarbazide derivatives [54]. Similar techniques have also been used to investigate the stability of phosphodiesterase type 5 (PDE5) in complexes with bioactive compounds from Mimosa pudica to understand their aphrodisiac performance [55]. Similarly, a pharmacoinformatics-based molecular docking and dynamics simulation analysis of bioactive components from Indian cuisine, rasam, was conducted against MAPK6 (mitogen-activated protein kinase 6), a family of serine/threonine protein kinases that is crucial in regulating extracellular signaling into a variety of cellular functions, including ROS production [56].
MM-GBSA also validates the molecular docking and MD simulation studies as it shows binding energy ranging from −49 to −79 kcal/mol in all the protein–ligand complexes studied. MM-GBSA has previously been used to study and validate in silico interaction of ligands with SARS-CoV-2 protease [57]. A similar MD simulation approach has also been utilized for screening substrate analog inhibitors of L-Ornithine-N5-monooxygenase (PvdA) to control Pseudomonas aeruginosa infections [58].
Analysis of ligands for ADMET properties using SwissADME and pkCSM shows that the molecular weight of all six hits was larger than 400g/mol. Only meliantrol was expected to have moderate water solubility. All other hits had low water solubility. The right solvation and absorption into the host body are made possible by the ligands’ solubility. Additionally, it helps with the formulation of the solvent or delivery vehicle. Except for meliantrol, none of the selected compounds exhibited good GI absorption. All the hits broke one out of five of Lipinski’s rules because of their poor solubility. The Bioavailability RADAR plot depicts the overall drug-likeness of the top binding ligands. Meliantrol fulfills all criteria for drug-likeness, while other ligands fail due to their low solubility in water. The BOILED-Egg plot displays all of the hits in the area of intestinal absorption except citrostadienol. None of the ligands were predicted to penetrate the blood–brain barrier. Analysis of ADMET properties of an array of ligands using Swiss ADME and pkCSM tools has been previously followed by many researchers to control S. mutans biofilm [59,60].

4. Materials and Methods

4.1. Target Proteins in S. mutans

The target proteins were selected based on our earlier study, which demonstrated an in silico subtractive proteomics approach for screening drug targets in S. mutans UA159. Among the 13 novel drug targets that were identified, 7 proteins were found to be involved in quorum sensing [21]. Hence, these target proteins have been subsequently used in the current study. The target proteins include two putative response regulators, a signal peptidase, a putative EEP protein-like protein, a putative transmembrane protein-permease, a protein translocase subunit, and a membrane protein insertase (Table 7).

4.2. Selection of Ligands

The bioactive compounds from three Indian medicinal plants, Azadirachta indica, Morinda citrifolia, and Salvodora persica, were selected as ligands from the literature, as well as various databases such as Dr. Duke’s Phytochemical and Ethnobotanical Databases (https://phytochem.nal.usda.gov/; accessed on 11 April 2022), IMPPAT (https://cb.imsc.res.in/imppat/; accessed on 13 April 2022), PubChem (https://pubchem.ncbi.nlm.nih.gov/; accessed on 15 April 2022), and Drug Bank (https://go.drugbank.com/; accessed on 18 April 2022). However, among the hundreds of compounds, the ones with a background of antimicrobial activity in the literature were selected. The screening process involved manually searching each compound in the PubMed database for its antimicrobial activity against S. mutans or other oral pathogens. Finally, 110 compounds were shortlisted for analysis. The 3D structures of these ligands were either collected from PubChem or drawn using ChemSketch Version 12.00 (Table 8).

4.3. Computational Analysis

Computational analysis was performed using Debian Linux Operating System running on a 3.10 GHz Intel® Core™ i5-4440 CPU (Acer, Bengaluru, India) with 8 GB RAM. For the analysis of docking results, the Accelrys® Discovery Studio Visualizer (Accelrys Software Inc., San Diego, CA, USA) that runs on the Windows operating system was used.

4.4. Homology Modeling and Model Validation

The 3D structures of target proteins were modeled using the SWISS-MODEL online tool (https://swissmodel.expasy.org/ accessed on 20 April 2022) [61]. The modeled structures were assessed and validated using the Ramachandran plot analysis.

4.5. Molecular Docking

The target proteins were prepared by removing the water molecules using Discovery Studio. The parameters for docking were set up by using the graphical user interface (GUI) of AutoDock Tools [62]. The preparation of the target protein structure involves the addition of hydrogen atoms and formal charges. The number of rotatable bonds and determination of root for the ligand were set to default. Later, molecular docking was performed using AutoDock 4.2. The AutoGrid4 program of AutoDock allowed the generation of grid maps for target proteins embedded in a three-dimensional grid of manually set parameters (Table 9).
The binding energies between the target protein and ligand were calculated by running the AutoDock4 program using pre-set grid maps. The result analysis was performed using the GUI of AutoDock Tools and Discovery Studio. The binding energies for various conformations of the ligand with the target proteins were determined, and the best conformation was chosen based on the binding energy and the number of hydrogen bonds that they formed with the protein.

4.6. Molecular Dynamics

The protein–ligand complexes (PLCs) were prepared using PyMOL [63], and these PLCs were taken as input for Molecular Dynamics (MD) simulations (Table 10). The MD simulation of PLCs aids in the visualization of target binding sites and also provides information concerning the binding stability of the PLCs under physiological conditions. The Desmond module of Schrödinger was used to perform MD simulations (academic license, Version 2020-1) [64]. Initially, an explicit water model was prepared as an orthorhombic simulation box with Simple Point-Charge (SPC). This was designed with the builder panel of the system in such a way that a minimum distance of 10 Å is maintained between the protein surface and the solvent surface. Then, the PLCs with receptors were solvated with the orthorhombic TIP3P water model [65]. The addition of counter ions and limitation of salt concentration of the physiological system to 0.15M was performed for the neutralization of the solvated system (Table 11).
Then, the PLC system was designated with the OPLS AA force field [66]. For the RESPA (reversible reference system propagator algorithms) integrator [67], Nose–Hoover chain thermostat [68], and Martyna–Tobias–Klein barostat, a relaxation time of two seconds was used. The equilibrated system was used to perform the final production process of MD simulations. In this study, the default parameters set for relaxation before MD simulations include a time duration of 100 ns, temperature of 310 K, a pressure of 1.0 bar with NPT (Isothermal–Isobaric ensemble, constant temperature, constant pressure, constant number of particles) ensemble [69] (Table 12).
The MD simulation tool was used to run the simulation, and the output file was used to retrieve trajectories and create the movie. The output file in .cms format was imported into the software, and the created simulation movie was exported at a higher resolution of 1280 × 1024 px for better quality. The trajectory of the overall MD simulation process was written in 1000 frames. The frames of the protein backbone were aligned to the backbone of the first frame. This provides a better understanding of the stability of the PLCs during MD simulations. Finally, the Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuation (RMSF) of PLCs and the simulation interaction diagrams were analyzed [57,70]. For each protein–ligand complex, three simulations studies were run, and the best run was taken for analysis.
The free binding energies of the protein and ligand complexes were also examined using Molecular Mechanics, the Generalized Born model, and Solvent Accessibility (MM-GBSA). The “thermal_MMGBSA.py” script and “MM-GBSA ΔG Bind:” parameter from the Prime/Desmond module of the Schrodinger suite were employed to carry out the post-MM-GBSA analysis [71].

4.7. ADMET Analysis

The ADMET properties of the high binding ligands were analyzed using the SwissADME online tool (http://www.swissadme.ch; accessed on 10 May 2022) [72] and pkCSM online tools (http://biosig.unimelb.edu.au/pkcsm/prediction; accessed on 10 May 2022) [73]. The input was given as SMILES notation of the ligands, and the output was retrieved as a “.csv” file. In SwissADME, the BOILED-Egg plot and radar plot for all the leads were also retrieved as a single image [74]. The ADMET analysis helps to predict the drug-likeness of the top hits.

5. Conclusions

S. mutans is a major contributor to dental caries. In this study, an in silico approach was employed to identify a plant metabolite with potential efficacy against specific target proteins of S. mutans involved in quorum sensing. The selected target proteins included: Membrane protein insertase YidC 2, Permease OppC, Putative Eep protein-like protein, Putative transmembrane protein, Putative response regulator CiaR, Putative response regulator SpaR, and Signal peptidase I LepC. The three-dimensional models of these target proteins were generated using SWISS-MODEL and validated using the Ramachandran plot. Plant metabolites derived from A. indica (Neem), M. citrifolia (Noni), and S. persica (Miswak) were evaluated for their potential binding affinity to the selected target proteins. Molecular docking studies were performed using AutoDock Tools. From a total of 110 ligands, 6 top hits were identified for each target protein. These ligands, namely campesterol, meliantrol, citrostadienol, stigmasterol, isofucosterol, and ursolic acid, were subjected to molecular simulation analysis to assess their stability and interaction patterns. The highly interacted amino acid residues identified were ASP86, ARG182, ILE179, GLU143, ASP237, PRO101, and VAL84, corresponding to the proteins CiaR, LepC, OppC, SecA, SMU1784c, SpaR, YidC2, respectively. Furthermore, the ADME characteristics of the identified ligands were evaluated using the SwissADME program. Collectively, the results suggest that these phytosterols have the potential to serve as quorum-quenching agents. However, further in vitro analysis is required to confirm their efficacy against S. mutans and to enable their application in the treatment of dental caries.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28145514/s1. Figure S1: Ramachandran plot of modelled target proteins using structure assessment tool of the SWISS-MODEL: (a) ciaR; (b) lepC; (c) oppC; (d) secA; (e) smu1784c; (f) spaR and (g) yidc2, (h) represents the transmembrane localization of target protein yidc2. Figure S2: 3D structure of target proteins predicted using SWISS-MODEL online tool: (a) ciaR; (b) lepC; (c) oppC; (d) secA; (e) smu1784c; (f) spaR and (g) yidc2.

Author Contributions

Conceptualization, S.C.V.M. and K.S.; methodology, S.C.V.M., J.M., E.B., P.M., M.S., E.T., J.C.R., S.R.K.P., S.K. and V.B.; software, M.S. and P.M.; validation, K.S., S.K. and S.R.K.P.; formal analysis, S.C.V.M., E.B., P.M. and E.T.; investigation, S.C.V.M. and J.C.R.; resources, K.S. and P.M.; data curation, S.C.V.M.; writing—original draft preparation, S.C.V.M. and E.T.; writing—review and editing, S.C.V.M., E.T., S.R.K.P., S.K. and K.S.; visualization, S.C.V.M. and K.S.; supervision, K.S.; project administration, K.S.; funding acquisition, P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Engineering Research Board of India (Grant number: EMR/2016/003035 to K.S.) and the Department of Biotechnology (Grant number: BT/PR36633/TRM/120/277/2020 to K.S. and S.K.). The article processing charges were funded by the National Science Centre, Poland, Grant number 2019/35/B/NZ8/04523 to P.M.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Additional data may be available upon personal request.

Acknowledgments

S.C.V.M. thanks the Kalasalingam Academy of Research and Education for University Research Fellowship. S.M. is thankful to BITS-Pilani, Pilani campus, for providing the Computational facilities.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Loesche, W.J. Role of Streptococcus mutans in human dental decay. Microbiol. Rev. 1986, 50, 353–380. [Google Scholar] [CrossRef] [PubMed]
  2. Ajdić, D.; McShan, W.M.; McLaughlin, R.E.; Savić, G.; Chang, J.; Carson, M.B.; Primeaux, C.; Tian, R.; Kenton, S.; Jia, H.; et al. Genome sequence of Streptococcus mutans UA159, a cariogenic dental pathogen. Proc. Natl. Acad. Sci. USA 2002, 99, 14434–14439. [Google Scholar] [CrossRef] [PubMed]
  3. Lemos, J.A.; Quivey, R.G.; Koo, H.; Abranches, J. Streptococcus mutans: A new gram-positive paradigm? Microbiology 2013, 159, 436–445. [Google Scholar] [CrossRef] [PubMed]
  4. Lemos, J.A.; Palmer, S.R.; Zeng, L.; Wen, Z.T.; Kajfasz, J.K.; Freires, I.A.; Abranches, J.; Brady, L.J. The Biology of Streptococcus mutans. Microbiol. Spectr. 2019, 7, 1–18. [Google Scholar] [CrossRef]
  5. Salazar, J.A.G.; Díaz, J.R.M.; Segura, E.R.; Ferrara, J.G.T. What are the origins of growing microbial resistance? Both Lamarck and Darwin were right. Expert. Rev. Anti. Infect. Ther. 2021, 19, 563–569. [Google Scholar] [CrossRef]
  6. Kaur, G.; Rajesh, S.; Princy, S.A. Plausible drug targets in the Streptococcus mutans quorum sensing pathways to combat dental biofilms and associated risks. Indian J. Microbiol. 2015, 55, 349–356. [Google Scholar] [CrossRef]
  7. Muras, A.; Mayer, C.; Romero, M.; Camino, T.; Ferrer, M.D.; Mira, A.; Otero, A. Inhibition of Streptococcus mutans biofilm formation by extracts of Tenacibaculum sp. 20J, a bacterium with wide-spectrum quorum quenching activity. J. Oral Microbiol. 2018, 10, 1429788. [Google Scholar] [CrossRef] [Green Version]
  8. Wu, C.; Ayala, E.A.; Downey, J.S.; Merritt, J.; Goodman, S.D.; Qi, F. Regulation of ciaXRH operon expression and identification of the CiaR regulon in Streptococcus mutans. J. Bacteriol. 2010, 192, 4669–4679. [Google Scholar] [CrossRef] [Green Version]
  9. Nepomuceno, R.S.L.; Tavares, M.B.; Lemos, J.A.; Griswold, A.R.; Ribeiro, J.L.; Balan, A.; Guimaraes, K.S.; Cai, S.; Burne, R.A.; Ferreira, L.C.D.S.; et al. The oligopeptide (opp) gene cluster of Streptococcus mutans: Identification, prevalence, and characterization. Oral Microbiol. Immunol. 2007, 22, 277–284. [Google Scholar] [CrossRef]
  10. Xie, K.; Dalbey, R.E. Inserting proteins into the bacterial cytoplasmic membrane using the Sec and YidC translocases. Nat. Rev. Microbiol. 2008, 6, 234–244. [Google Scholar] [CrossRef]
  11. Galvão, L.C.D. Construction and Phenotypic Characterization of Streptococcus mutans Mutant Strains of Genes Related to Virulence. P.h.D. Thesis, The State University of Campinas (UNICAMP), Sao Paulo, Brazil, 2016. [Google Scholar]
  12. Rolim, W.R.; Lamilla, C.; Pieretti, J.C.; Díaz, M.; Tortella, G.R.; Cristina Diez, M.; Barrientos, L.; Seabra, A.B.; Rubilar, O. Comparison of antibacterial and antibiofilm activities of biologically synthesized silver nanoparticles against several bacterial strains of medical interest. Energy Ecol. Environ. 2019, 4, 143–149. [Google Scholar] [CrossRef]
  13. Mehrishi, P.; Agarwal, P.; Broor, S.; Sharma, A. Antibacterial and antibiofilm properties of medicinal plant extracts against multi drug resistant Staphylococcus species and non-fermenter bacteria. J. Pure Appl. Microbiol. 2020, 14, 403–413. [Google Scholar] [CrossRef] [Green Version]
  14. Rezoqe, N.S.; Dhahir, S.A.; Ascar, E.F. Biological effect of aqueous & alcoholic extracts of green tea leaves of some pathogenic bacteria in vitro. J. Kufa Chem. Sci. 2012, 4, 1–5. [Google Scholar]
  15. Fajriani, S.S.; Handayani, H.; Putri, D. The role of green tea extract on inhibiting Porphyromonas gingivalis as a major periodontitis pathogen: In vitro study. Syst. Rev. Pharm. 2020, 11, 152–155. [Google Scholar]
  16. Kim, Y.G.; Baltabekova, A.Z.; Zhiyenbay, E.E.; Aksambayeva, A.S.; Shagyrova, Z.S.; Khannanov, R.; Ramanculov, E.M.; Shustov, A.V. Recombinant Vaccinia virus-coded interferon inhibitor B18R: Expression, refolding and a use in a mammalian expression system with a RNA-vector. PLoS ONE 2017, 12, e0189308. [Google Scholar] [CrossRef] [Green Version]
  17. Murugan, K.; Sekar, K.; Sangeetha, S.; Ranjitha, S.; Sohaibani, S.A. Antibiofilm and quorum sensing inhibitory activity of Achyranthes aspera on cariogenic Streptococcus mutans: An in vitro and in silico study. Pharm. Biol. 2013, 51, 728–736. [Google Scholar] [CrossRef]
  18. Li, X.; Yin, L.; Ramage, G.; Li, B.; Tao, Y.; Zhi, Q.; Lin, H.; Zhou, Y. Assessing the impact of curcumin on dual-species biofilms formed by Streptococcus mutans and Candida albicans. Microbiol. Open 2019, 8, e937. [Google Scholar] [CrossRef] [Green Version]
  19. Zayed, S.M.; Aboulwafa, M.M.; Hashem, A.M.; Saleh, S.E. Biofilm formation by Streptococcus mutans and its inhibition by green tea extracts. AMB Express 2021, 11, 73. [Google Scholar] [CrossRef]
  20. Wolinsky, L.E.; Mania, S.; Nachnani, S.; Ling, S. The inhibiting effect of aqueous Azadirachta indica (Neem) extract upon bacterial properties influencing in vitro plaque formation. J. Dent. Res. 1995, 75, 816–822. [Google Scholar] [CrossRef]
  21. Marimuthu, S.C.V.; Ravinarayanan, H.; Rosy, J.C.; Sundar, K. Mining the proteome of Streptococcus mutans for putative drug targets. Infect. Disord. Drug Targets 2021, 21, 429–438. [Google Scholar] [CrossRef]
  22. Qi, F.; Merritt, J.; Lux, R.; Shi, W. Inactivation of the ciaH Gene in Streptococcus mutans diminishes mutacin production and competence development, alters sucrose-dependent biofilm formation, and reduces stress tolerance. Infect. Immun. 2004, 72, 4895–4899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Shanker, E.; Federle, M.J. Quorum sensing regulation of competence and bacteriocins in Streptococcus pneumoniae and mutans. Genes 2017, 8, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhu, B.; Ge, X.; Stone, V.; Kong, X.; El-Rami, F.; Liu, Y.; Kitten, T.; Xu, P. CiaR impacts biofilm formation by regulating an arginine biosynthesis pathway in Streptococcus sanguinis SK36. Sci. Rep. 2017, 7, 17183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Nakano, K.; Lapirattanakul, J.; Nomura, R.; Nemoto, H.; Alaluusua, S.; Grönroos, L.; Vaara, M.; Hamada, S.; Ooshima, T.; Nakagawa, I. Streptococcus mutans clonal variation revealed by multilocus sequence typing. J. Clin. Microbiol. 2007, 45, 2616–2625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Guo, Q.; Ahn, S.J.; Kaspar, J.; Zhou, X.; Burne, R.A. Growth phase and pH influence peptide signaling for competence development in Streptococcus mutans. J. Bacteriol. 2014, 196, 227–236. [Google Scholar] [CrossRef] [Green Version]
  27. Lewis, N.E.; Brady, L.J. Breaking the bacterial protein targeting and translocation model: Oral organisms as a case in point. Mol. Oral Microbiol. 2015, 30, 186–197. [Google Scholar] [CrossRef] [Green Version]
  28. Petersen, F.C.; Assev, S.; Van Der Mei, H.C.; Busscher, H.J.; Scheie, A.A. Functional variation of the antigen I/II surface protein in Streptococcus mutans and Streptococcus intermedius. Infect. Immun. 2002, 70, 249–256. [Google Scholar] [CrossRef] [Green Version]
  29. Saito, M.; Otake, S.; Ohmura, M.; Hirasawa, M.; Takada, K.; Mega, J.; Takahashi, I.; Kiyono, H.; McGhee, J.R.; Takeda, Y.; et al. Protective immunity to Streptococcus mutans induced by nasal vaccination with surface protein antigen and mutant cholera toxin adjuvant. J. Infect. Dis. 2001, 183, 823–826. [Google Scholar] [CrossRef]
  30. Palmer, S.R.; Ren, Z.; Hwang, G.; Liu, Y.; Combs, A.; Söderström, B.; Lara Vasquez, P.; Khosravi, Y.; Brady, L.J.; Koo, H.; et al. Streptococcus mutans yidC1 and yidC2 impact cell envelope biogenesis, the biofilm matrix, and biofilm biophysical properties. J. Bacteriol. 2019, 201, e00396-18. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, Z.C.; De Keyzer, J.; Berrelkamp-lahpor, G.A.; Driessen, A.J.M. Interaction of Streptococcus mutans YidC1 and YidC2 with translating and non-translating ribosomes. J. Bacteriol. 2013, 195, 4545–4551. [Google Scholar] [CrossRef] [Green Version]
  32. Petschacher, B.; Staunig, N.; Müller, M.; Schürmann, M.; Mink, D.; De Wildeman, S.; Gruber, K.; Glieder, A. Cofactor specificity engineering of Streptococcus mutans NADH oxidase 2 for NAD(P)+ regeneration in biocatalytic oxidations. Comput. Struct. Biotechnol. J. 2014, 9, e201402005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lu, S.; Wu, J.; Gao, Y.; Han, G.; Ding, W.; Huang, X. Cloning, expression, characterization and homology modeling of a novel water-forming NADH oxidase from Streptococcus mutans ATCC 25175. Int. J. Biol. Macromol. 2018, 113, 1073–1079. [Google Scholar] [CrossRef]
  34. Osorio, M.I.; Zúñiga, M.A.; Mendoza, F.; Jaña, G.A.; Jiménez, V.A. Modulation of glucan-enzyme interactions by domain V in GTF-SI from Streptococcus mutans. Proteins Struct. Funct. Bioinform. 2019, 87, 74–80. [Google Scholar] [CrossRef] [PubMed]
  35. Pandey, S.D.; Biswas, I. Expression of an Extracellular Protein (SMU.63) Is Regulated by SprV in Streptococcus mutans. Appl. Environ. Microbiol. 2020, 86, e01647-20. [Google Scholar] [CrossRef] [PubMed]
  36. di Cologna, N.M.; Samaddar, S.; Valle, C.A.; Vargas, J.; Aviles-Reyes, A.; Morales, J.; Ganguly, T.; Pileggi, R.; Brady, L.J.; Lemos, J.A.; et al. Amyloid aggregation of Streptococcus mutans Cnm influences its collagen-binding activity. Appl. Environ. Microbiol. 2021, 87, e01149-21. [Google Scholar] [CrossRef]
  37. Al-Khafaji, Z.M.; Mahmood, S.B.; Mahmood, A.B. In silico identification of potential inhibitors for Streptococcus mutans signal peptidase I using structure-based drug design. Lat Am. J. Pharm. 2021, 40, 221–227. [Google Scholar]
  38. Hosseini, B.; Behbahani, M.; Dini, G.; Mohabatkar, H.; Keyhanfar, M. Investigating the anti-streptococcal biofilm effect of ssDNA aptamer-silver nanoparticles complex on a titanium-based substrate. RSC Adv. 2022, 12, 24876–24886. [Google Scholar] [CrossRef]
  39. Li, D.; Liu, Y.; Li, K.; Zhang, L. Targeting F13 from monkeypox virus and variola virus by tecovirimat: Molecular simulation analysis. J. Infect. 2022, 85, e99–e101. [Google Scholar] [CrossRef]
  40. Mallano, A.I.; Nisa, Z.U.; Khaliq, B.; Ali, N.; Ali, Q.; Chao, C.; Yanming, Z. Cloning and in silico characterization of an abiotic stress-inducible U-box domain-containing protein gene Gs PUB8 from Glycine soja. Sci. Rep. 2022, 12, 17146. [Google Scholar] [CrossRef]
  41. 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. 2022, 40, 2067–2081. [Google Scholar] [CrossRef]
  42. Julaton, T.; Taclendo, A.; Oyong, G.; Rempillo, O.; Galvez, M.C.; Vallar, E. In silico insights on the pro-inflammatory potential of polycyclic aromatic hydrocarbons and the prospective anti-inflammatory capacity of Andrographis paniculata phytocompounds. Int. J. Environ. Res. Public Health 2022, 19, 8588. [Google Scholar] [CrossRef] [PubMed]
  43. Deshpande, S.H.; Muhsinah, A.B.; Bagewadi, Z.K.; Ankad, G.M.; Mahnashi, M.H.; Yaraguppi, D.A.; Shaikh, I.A.; Khan, A.A.; Hegde, H.V.; Roy, S. In silico study on the interactions, molecular docking, dynamics and simulation of potential compounds from Withania somnifera (L.) Dunal root against cancer by targeting KAT6A. Molecules 2023, 28, 1117. [Google Scholar] [CrossRef] [PubMed]
  44. Lu, S.Y.; Jiang, Y.J.; Lv, J.; Wu, T.X.; Yu, Q.S.; Zhu, W.L. Molecular docking and molecular dynamics simulation studies of GPR40 receptor–agonist interactions. J. Mol. Graph. Modell. 2010, 28, 766–774. [Google Scholar] [CrossRef] [PubMed]
  45. Prabhu, D.S.; Rajeswari, V.D. In vitro and in silico analyses of Vicia faba L. on Peroxisome proliferator–activated receptor gamma. J. Cell. Biochem. 2018, 119, 7729–7737. [Google Scholar] [CrossRef]
  46. Mirza, M.U.; Mirza, A.H.; Ghori, N.U.; Ferdous, S. Glycyrrhetinic acid and E. resveratroloside act as potential plant derived compounds against dopamine receptor D3 for Parkinson’s disease: A pharmacoinformatics study. Drug Des. Dev. Ther. 2014, 18, 187–198. [Google Scholar] [CrossRef] [Green Version]
  47. Ugwu, D.I.; Okoro, U.C.; Mishra, N.K.; Okafor, S.N. Novel Phenoxazinones as potent agonist of PPAR-α: Design, synthesis, molecular docking and in vivo studies. Lipids Health Dis. 2018, 17, 120. [Google Scholar] [CrossRef] [Green Version]
  48. Ochoa, R.; Martínez-Pabón, M.C.; Arismendi-Echeverri, M.A.; Rendón-Osorio, W.L.; Muskus-López, C.E. In silico search of inhibitors of Streptococcus mutans for the control of dental plaque. Arch. Oral Biol. 2017, 83, 68–75. [Google Scholar] [CrossRef]
  49. Islam, M.I.; Arokiyaraj, S.; Kuralarasan, M.; Kumar, V.S.; Harikrishnan, P.; Saravanan, S.; Ashok, G.; Chellappandian, M.; Bharanidharan, R.; Muralidaran, S.; et al. Inhibitory potential of EGCG on Streptococcus mutans biofilm: A new approach to prevent cariogenesis. Microb. Pathog. 2020, 143, 104129. [Google Scholar] [CrossRef]
  50. Vadloori, B.; Sharath, A.K.; Prabhu, N.P.; Maurya, R. Homology modelling, molecular docking, and molecular dynamics simulations reveal the inhibition of Leishmania donovani dihydrofolate reductase-thymidylate synthase enzyme by Withaferin-A. BMC Res. Notes 2018, 11, 246. [Google Scholar] [CrossRef] [Green Version]
  51. Salehi, N.; Attaran, B.; Eskini, N.; Esmaeili, M.; Sharifirad, A.; Sadeghi, M.; Mohammadi, M. New insights into resistance of Helicobacter pylori against third-and fourth-generation fluoroquinolones: A molecular docking study of prevalent GyrA mutations. Helicobacter 2019, 24, e12628. [Google Scholar] [CrossRef]
  52. Naz, A.; Asif, S.; Alwutayd, K.M.; Sarfaraz, S.; Abbasi, S.W.; Abbasi, A.; Alenazi, A.M.; Hasan, M.E. Repurposing FIASMAs against Acid Sphingomyelinase for COVID-19: A computational molecular docking and dynamic simulation approach. Molecules 2023, 28, 2989. [Google Scholar] [CrossRef] [PubMed]
  53. Selvaraj, C.; Priya, R.B.; Lee, J.K.; Singh, S.K. Mechanistic insights of SrtA-LPXTG blockers targeting the transpeptidase mechanism in Streptococcus mutans. RSC Adv. 2015, 5, 100498–100510. [Google Scholar] [CrossRef]
  54. Miazga-Karska, M.; Wos, M.; Kaczor, A.A.; Pachuta-Stec, A.; Ginalska, G.; Pitucha, M. Inhibitory effects of 1, 4-disubstituted thiosemicarbazide derivatives on Streptococcus mutans and Streptococcus sanguinis mono-species biofilms. Lett. Drug Des. Discov. 2018, 15, 843–850. [Google Scholar] [CrossRef]
  55. Palanichamy, C.; Pavadai, P.; Panneerselvam, T.; Arunachalam, S.; Babkiewicz, E.; Ram Kumar Pandian, S.; Shanmugampillai Jeyarajaguru, K.; Nayak Ammunje, D.; Kannan, S.; Chandrasekaran, J.; et al. Aphrodisiac performance of bioactive compounds from Mimosa pudica Linn.: In silico molecular docking and dynamics simulation approach. Molecules 2022, 27, 3799. [Google Scholar] [CrossRef] [PubMed]
  56. Kalimuthu, A.K.; Panneerselvam, T.; Pavadai, P.; Pandian, S.R.K.; Sundar, K.; Murugesan, S.; Ammunje, D.N.; Kumar, S.; Arunachalam, S.; Kunjiappan, S. Pharmacoinformatics-based investigation of bioactive compounds of Rasam (South Indian recipe) against human cancer. Sci. Rep. 2021, 11, 21488. [Google Scholar] [CrossRef]
  57. Kumar, B.K.; Faheem, N.; Sekhar, K.V.G.C.; Ojha, R.; Prajapati, V.K.; Pai, A.; Murugesan, S. Pharmacophore based virtual screening, molecular docking, molecular dynamics and MM-GBSA approach for identification of prospective SARS-CoV-2 inhibitor from natural product databases. J. Biomol. Struct. Dyn. 2022, 40, 1363–1386. [Google Scholar] [CrossRef]
  58. Rosy, J.C.; Babkiewicz, E.; Maszczyk, P.; Mrówka, P.; Kumar, B.K.; Murugesan, S.; Kunjiappan, S.; Sundar, K. l-Ornithine-N5-monooxygenase (PvdA) Substrate Analogue Inhibitors for Pseudomonas aeruginosa Infections Treatment: Drug Repurposing Computational Studies. Biomolecules 2022, 12, 887. [Google Scholar] [CrossRef]
  59. Atta, L.; Khalil, R.; Khan, K.M.; Zehra, M.; Saleem, F.; Nur-e-Alam, M.; Ul-Haq, Z. Virtual screening, synthesis and biological evaluation of Streptococcus mutans mediated biofilm inhibitors. Molecules 2022, 27, 1455. [Google Scholar] [CrossRef]
  60. Azzam, K.A.L. SwissADME and pkCSM webservers predictors: An integrated online platform for accurate and comprehensive predictions for in silico ADME/T properties of artemisinin and its derivatives. Kompleks. Ispolz. Miner. Syra 2023, 325, 14–21. [Google Scholar] [CrossRef]
  61. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.; Clamp, M.; Barton, G.J. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [Green Version]
  62. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Schrödinger, LLC. The PyMOL Molecular Graphics System, Version 2.0; Schrödinger, LLC: New York, NY, USA, 2017. [Google Scholar]
  64. D.E. Shaw Research; Schrödinger. Schrödinger Release 2020-4: Desmond Molecular Dynamics System; Maestro-Desmond Interoperability Tools; D.E. Shaw Research: New York, NY, USA; Schrödinger: New York, NY, USA, 2020. [Google Scholar]
  65. Mark, P.; Nilsson, L. Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J. Phys. Chem. A 2001, 105, 9954–9960. [Google Scholar] [CrossRef]
  66. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  67. Tuckerman, M.; Berne, B.J.; Martyna, G.J. Reversible multiple time scale molecular dynamics. J. Chem. Phys. 1992, 97, 1990–2001. [Google Scholar] [CrossRef] [Green Version]
  68. Cheng, A.; Merz, K.M. Application of the Nosé-Hoover chain algorithm to the study of protein dynamics. J. Phys. Chem. 1996, 100, 1927–1937. [Google Scholar] [CrossRef]
  69. Kalibaeva, G.; Ferrario, M.; Ciccotti, G. Constant pressure-constant temperature molecular dynamics: A correct constrained NPT ensemble using the molecular virial. Mol. Phys. 2003, 101, 765–778. [Google Scholar] [CrossRef]
  70. Kumar, A.; Patel, S.; Bhatkar, D.; Sarode, S.C.; Sharma, N.K. Design, synthesis and evaluation of novel β-carboline ester analogues as potential anti-leishmanial agents. J. Biomol. Struct. Dyn. 2022, 40, 12592–12607. [Google Scholar] [CrossRef]
  71. Masetti, M.; Falchi, F.; Gioia, D.; Recanatini, M.; Ciurli, S.; Musiani, F. Targeting the protein tunnels of the urease accessory complex: A theoretical investigation. Molecules 2020, 25, 2911. [Google Scholar] [CrossRef]
  72. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  73. Pires, D.E.V.; Blundell, T.L.; Ascher, D.B. pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem. 2015, 58, 4066–4072. [Google Scholar] [CrossRef]
  74. Daina, A.; Zoete, V. A BOILED-Egg to predict gastrointestinal absorption and brain penetration of small molecules. ChemMedChem 2016, 11, 1117–1121. [Google Scholar] [CrossRef] [Green Version]
Molecules 28 05514 g001 550
Figure 1. 3D interaction of target proteins with their respective top hit ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 1. 3D interaction of target proteins with their respective top hit ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g001
Molecules 28 05514 g002 550
Figure 2. 2D interaction of target proteins with their respective top hit ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 2. 2D interaction of target proteins with their respective top hit ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g002
Molecules 28 05514 g003 550
Figure 3. RMSD plot of protein–ligand complexes: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 3. RMSD plot of protein–ligand complexes: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g003
Molecules 28 05514 g004 550
Figure 4. Hydrogen and hydrophobic interactions of protein–ligand complexes: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 4. Hydrogen and hydrophobic interactions of protein–ligand complexes: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g004
Molecules 28 05514 g005 550
Figure 5. Protein–ligand contacts between target proteins and their respective ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 5. Protein–ligand contacts between target proteins and their respective ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g005
Molecules 28 05514 g006 550
Figure 6. Timeline representation of interaction and contacts between target proteins and their respective ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Figure 6. Timeline representation of interaction and contacts between target proteins and their respective ligands: (a) CiaR–campesterol; (b) LepC–meliantrol; (c) OppC–stigmasterol; (d) SecA–citrostadienol; (e) SMU1784c–citrostadienol; (f) SpaR–isofucosterol; (g) Yidc2–ursolic acid.
Molecules 28 05514 g006
Molecules 28 05514 g007 550
Figure 7. Bioavailability RADAR Plot analysis of the selected top hit ligands: (a) campesterol; (b) citrostadienol; (c) isofucosterol; (d) meliantrol; (e) stigmasterol; (f) ursolic acid.
Figure 7. Bioavailability RADAR Plot analysis of the selected top hit ligands: (a) campesterol; (b) citrostadienol; (c) isofucosterol; (d) meliantrol; (e) stigmasterol; (f) ursolic acid.
Molecules 28 05514 g007
Molecules 28 05514 g008 550
Figure 8. BOILED-Egg plot of selected high-binding ligands as predicted using SwissADME.
Figure 8. BOILED-Egg plot of selected high-binding ligands as predicted using SwissADME.
Molecules 28 05514 g008
Table
Table 1. List of top hits from different plant sources against each target protein, along with their binding energy (kcal/mol).
Table 1. List of top hits from different plant sources against each target protein, along with their binding energy (kcal/mol).
LigandsCIARLEPCSMU1784COPPCSPARSECAYIDC2Plant Source
5α-Androstan-16-one, cyclic ethylene mercaptole−7.47−9.5−7.46−6.23−8.29−9.26−8.24S. persica
Azadiradione−7.91−9.18−8.32−6.22−8.1−8.45−8.23A. indica
Azadirone−7.81−8.85−8.34−5.83−8.34−8.16−8.85A. indica
Campesterol−8.76−9.91−8.08−6.41−8.56−9.29−8.86A. indica
Citrostadienol−7.55−9.81−8.45−6.43−8.9−9.88−9.4A. indica
Isofucosterol−8.44−9.31−8.14−6.68−9.1−9.58−8.94A. indica
Margocin−7.24−8.85−7.58−5.94−8.11−7.81−8.07A. indica
Meliantrol−8.44−10.16−7.43−5.95−8.14−8.66−8.61A. indica
Nimbinene−7.41−7.45−7.88−5.38−7.39−7.99−7.39A. indica
Nimbione−7.09−7.61−7.63−5.91−7.65−7.7−7.46A. indica
Nimbolide−7.62−8.14−7.57−5.8−8.45−7.56−8.3A. indica
Nimbolin b−7.9−9.09−7.58−6.09−9.07−9.77−7.41A. indica
Nimocinol−7.95−8.55−8.27−5.98−8.59−7.84−8.41A. indica
Stigmasterol−7.78−9.87−8.35−6.75−8.76−9.12−9.04A. indica
Ursolic acid−7.98−9.42−8.31−6.12−8.14−7.54−9.73M. citrifolia
Vepinin−7.81−10.06−7.69−5.51−7.85−9−8.38A. indica
Zafaral−7.71−9.12−7.3−5.88−7.8−8.09−8.1A. indica
Table
Table 2. List of top hits from different plant sources against each target protein, along with their binding energy.
Table 2. List of top hits from different plant sources against each target protein, along with their binding energy.
Top BindersBinding Energy (kcal/mol)Target
Campesterol−8.76CIAR
Meliantrol−10.16LEPC
Citrostadienol−8.45SMU1784C
Stigmasterol−6.75OPPC
Isofucosterol−9.1SPAR
Citrostadienol−9.88SECA
Ursolic acid−9.73YIDC2
Table
Table 3. List of amino acid residues interacting with the ligand. The residues which interact with ligands through conventional hydrogen bonds during the docking process are bold-faced.
Table 3. List of amino acid residues interacting with the ligand. The residues which interact with ligands through conventional hydrogen bonds during the docking process are bold-faced.
Amino Acid Residues Interacting through Hydrogen Bond
CiaRASP86; LYS87; ASN144; GLY145; GLU147
LepCASN66; ASN68; ILE69; VAL70; LYS71; ARG72; VAL73; ILE74; LYS93; THR98; TYR151; LYS174; VAL176; GLY177; GLU178; VAL179; LYS180; PHE181; ARG182
OppCVAL148; ALA149; THR151; LEU155; TYR17; GLN177; SER184; SER185; VAL186; PRO187; ASN188
SecAGLU143; ASN464; ALA469; GLN470
SMU1784cASN230; LYS231; ASP237; ASN240
SpaRLYS100; PRO101; GLY103; GLN105; GLU106
Yidc2ARG73; LEU77; LEU81; GLN90; GLY137; GLN204; THR208; HIS250
Table
Table 4. MMGBSA Score of Protein–Ligand Complexes.
Table 4. MMGBSA Score of Protein–Ligand Complexes.
Protein–Ligand ComplexMMGBSA Score (kcal/mol)
CiaR–campesterol−59.03099478
LepC–meliantrol−49.87694241
OppC–stigmasterol−65.96184631
SecA–citrostadienol−60.69835805
SMU1784c–citrostadienol−57.96969001
SpaR–isofucosterol−58.58904881
Yidc2–ursolic acid−79.57923832
Table
Table 5. ADME properties of top hits assessed in SWISS-ADME.
Table 5. ADME properties of top hits assessed in SWISS-ADME.
MoleculeCampesterolCitrostadienolIsofucosterolMeliantrolStigmasterolUrsolic Acid
Physicochemical Properties
FormulaC28H48OC30H50OC29H48OC28H46O5C29H48OC30H48O3
Molecular Weight400.68426.72412.69462.66412.69456.7
#Heavy atoms293130333033
#Aromatic heavy atoms000000
Fraction Csp30.930.870.860.930.860.9
#Rotatable bonds555351
#H-bond acceptors111513
#H-bond donors111412
MR128.42137.56132.75130.66132.75136.91
TPSA20.2320.2320.2390.1520.2357.53
Lipophilicity
iLOGP4.925.415.153.765.013.71
XLOGP38.89.038.854.188.567.34
WLOGP7.638.197.944.177.87.09
MLOGP6.546.826.623.456.625.82
Silicos-IT Log P6.6376.883.466.865.46
Consensus Log P6.97.297.093.86.975.88
Water Solubility
ESOL Log S−7.54−7.84−7.64−5.14−7.46−7.23
ESOL Solubility (mg/mL)0.00001160.00000610.000009360.003320.00001430.0000269
ESOL Solubility (mol/l)0.0000000290.00000001430.00000002270.000007180.00000003460.0000000589
ESOL ClassPoorly solublePoorly solublePoorly solubleModerately solublePoorly solublePoorly soluble
Ali Log S−9.11−9.35−9.16−5.78−8.86−8.38
Ali Solubility (mg/mL)0.0000003130.0000001920.0000002860.0007640.0000005710.00000192
Ali Solubility (mol/L)0.000000000780.000000000450.0000000006920.000001650.000000001380.00000000421
Ali ClassPoorly solublePoorly solublePoorly solubleModerately solublePoorly solublePoorly soluble
Silicos-IT LogSw−5.79−5.97−5.83−3.14−5.47−5.67
Silicos-IT Solubility (mg/mL)0.0006420.0004570.0006160.3320.00140.000972
Silicos-IT Solubility (mol/L)0.00000160.000001070.000001490.0007180.000003390.00000213
Silicos-IT classModerately solubleModerately solubleModerately solubleSolubleModerately solubleModerately soluble
Pharmacokinetics
GI absorptionLowLowLowHighLowLow
BBB permeantNoNoNoNoNoNo
Pgp substrateNoNoNoYesNoNo
CYP1A2 inhibitorNoNoNoNoNoNo
CYP2C19 inhibitorNoNoNoNoNoNo
CYP2C9 inhibitorNoNoNoNoYesNo
CYP2D6 inhibitorNoNoNoNoNoNo
CYP3A4 inhibitorNoNoNoNoNoNo
log Kp (cm/s)−2.5−2.49−2.53−6.15−2.74−3.87
Druglikeness
Lipinski #violations1; MLOGP > 4.151; MLOGP > 4.151; MLOGP > 4.1501; MLOGP > 4.151; MLOGP > 4.15
Ghose #violations2; WLOGP > 5.6, #atoms > 703; WLOGP > 5.6, MR > 130, #atoms > 703; WLOGP > 5.6, MR > 130, #atoms > 702; MR > 130, #atoms > 703; WLOGP > 5.6, MR > 130, #atoms > 703; WLOGP > 5.6, MR > 130, #atoms > 70
Veber #violations000000
Egan #violations1; WLOGP > 5.881; WLOGP > 5.881; WLOGP > 5.8801; WLOGP > 5.881; WLOGP > 5.88
Muegge #violations2; XLOGP3 > 5, Heteroatoms < 22; XLOGP3 > 5, Heteroatoms < 22; XLOGP3 > 5, Heteroatoms < 202; XLOGP3 > 5, Heteroatoms < 21; XLOGP3 > 5
Bioavailability Score0.550.550.550.550.550.85
Medicinal Chemistry
PAINS #alerts000000
Brenk #alerts1; isolated_alkene1; isolated_alkene1; isolated_alkene1; isolated_alkene1; isolated_alkene1; isolated_alkene
Leadlikeness #violations2; MW > 350, XLOGP3 > 3.52; MW > 350, XLOGP3 > 3.52; MW > 350, XLOGP3 > 3.52; MW > 350, XLOGP3 > 3.52; MW > 350, XLOGP3 > 3.52; MW > 350, XLOGP3 > 3.5
Synthetic Accessibility6.176.226.156.716.216.21
# denotes number.
Table
Table 6. ADMET properties of top hit ligands assessed in pkCSM.
Table 6. ADMET properties of top hit ligands assessed in pkCSM.
DescriptorPredicted ValueUnit
CampesterolCitrostadienolIsofucosterolMeliantrolStigmasterolUrsolic acid
Molecular Weight400.691426.729412.702462.671412.702456.711g/mol
LogP7.63478.19097.94494.17177.80087.0895
#Rotatable Bonds555351
#Acceptors111512
#Donors111412
Surface Area180.674192.714186.349199.164186.349201.354
PropertyModel Name
AbsorptionWater solubility−7.194−6.74−6.917−4.096−6.882−3.193Numeric (log mol/L)
AbsorptionCaCO2 permeability1.2841.2851.2791.5261.281.286Numeric (log Papp in 10−6 cm/s)
AbsorptionIntestinal absorption (human)95.74995.90796.06197.23896.39100Numeric (% Absorbed)
AbsorptionSkin Permeability−2.756−2.691−2.7−2.726−2.702−2.732Numeric (log Kp)
AbsorptionP-glycoprotein substrateNoYesNoYesNoNoCategorical (Yes/No)
AbsorptionP-glycoprotein I inhibitorYesYesYesYesYesNoCategorical (Yes/No)
AbsorptionP-glycoprotein II inhibitorYesYesYesYesYesNoCategorical (Yes/No)
DistributionVDss (human)0.351−0.0280.104−0.0740.102−1.104Numeric (log L/kg)
DistributionFraction unbound (human)000000Numeric (Fu)
DistributionBBB permeability0.8040.8160.807−0.8480.814−0.182Numeric (log BB)
DistributionCNS permeability−1.43−1.232−1.326−1.84−1.326−1.118Numeric (log PS)
MetabolismCYP2D6 substrateNoNoNoNoNoNoCategorical (Yes/No)
MetabolismCYP3A4 substrateYesYesYesYesYesYesCategorical (Yes/No)
MetabolismCYP1A2 inhibitorNoNoNoNoNoNoCategorical (Yes/No)
MetabolismCYP2C19 inhibitorNoNoNoNoNoNoCategorical (Yes/No)
MetabolismCYP2C9 inhibitorNoNoNoNoNoNoCategorical (Yes/No)
MetabolismCYP2D6 inhibitorNoNoNoNoNoNoCategorical (Yes/No)
MetabolismCYP3A4 inhibitorNoNoNoNoNoNoCategorical (Yes/No)
ExcretionTotal Clearance0.5720.5850.6190.4460.6180.083Numeric (log mL/min/kg)
ExcretionRenal OCT2 substrateNoNoNoNoNoNoCategorical (Yes/No)
ToxicityAMES toxicityNoNoNoNoNoNoCategorical (Yes/No)
ToxicityMax. tolerated dose (human)−0.193−0.394−0.374−0.726−0.3850.239Numeric (log mg/kg/day)
ToxicityhERG I inhibitorNoNoNoNoNoNoCategorical (Yes/No)
ToxicityhERG II inhibitorYesYesYesNoYesNoCategorical (Yes/No)
ToxicityOral Rat Acute Toxicity (LD50)2.3552.962.8474.4882.8362.81Numeric (mol/kg)
ToxicityOral Rat Chronic Toxicity (LOAEL)1.1251.0961.1192.2491.1022.128Numeric (log mg/kg_bw/day)
ToxicityHepatotoxicityNoNoNoNoNoYesCategorical (Yes/No)
ToxicitySkin SensitisationNoNoNoNoNoNoCategorical (Yes/No)
ToxicityT. pyriformis toxicity0.6760.4050.480.2880.4810.285Numeric (log ug/L)
ToxicityMinnow toxicity−2.071−2.047−1.9882.749−1.952−1.204Numeric (log mM)
# denotes number.
Table
Table 7. List of selected target proteins involved in quorum sensing.
Table 7. List of selected target proteins involved in quorum sensing.
UniProt Id.Gene NameProtein Name
Q8DU28CiaRPutative response regulator CiaR
Q8DSC7LepCSignal peptidase I
Q8DW23OppCPutative transmembrane protein, permease OppC
Q8DSF0SecAProtein translocase subunit SecA
Q8DSJ8SMU1784cPutative Eep protein-like protein
Q8DV49SMU659Putative response regulator SpaR
Q8DSP8YidC2Membrane protein insertase YidC 2
Table
Table 8. List of ligands selected from three plants, A. indica, M. citrifolia, and S. persica.
Table 8. List of ligands selected from three plants, A. indica, M. citrifolia, and S. persica.
Azadirachta indica
6-desacetyllnimbineneMargocinOleic acid
Arachidic acidMeliacinanhydridePalmitic acid
Azadirachtin AMeliantrolPalmitoleic acid
Azadirachtin BNimbandiolQuercetin
AzadiradioneNimbidiolSalanin
AzadironeNimbinStearic acid
Behenic acidNimbineneStigmasterol
Beta-sitosterolNimbinoneVepinin
CampesterolNimbioneVilasinin
CitrostadienolNimbolicinZafaral
FraxinelloneNimbolideα-tocopherol
Gadoleic acidNimbolin Aβ-tocopherol
IsomeldeninNimbolin Bγ-tocopherol
Linoleic acidNimocinolIsofucosterol
Linolenic acidNimocinolideδ-tocopherol
Morinda citrifolia
2-methoxy-1,3,6-trihydroxyanthraquinoneDeacetylasperulosideNordamnacanthal
6R-hydroxyadoxosideDehydromethoxygaertnerosideOctanoic acid (caprylic acid)
Americanin AEpi-dihydrocorninRetinol
Ascorbic acidL-asperulosideRubiadin
Asperulosidic acidMethyl alpha-d-fructofuranosideRubiadin-1-methyl ether
BorreriageninMethyl beta-d-fructofuranosideRutin
Caproic acidMorindacinScopoletin
Citrifolinin B epimer aMorindoneUrsolic acid
Citrifolinin B epimer bNarcissoside
CitrifolinosideNicotifloroside
Salvadora persica
1H-Pyrazole-1-carbothioamideCatechinPyrrolidine
1-triacantanolCineoleSalvadoside
3, 5-Dithiahexanol 5, 5-dioxideEpicatechinSalvadourea
5-O-caffeoylquinic acidFarnesolSyringin
5α-Androstan-16-one, cyclic ethylene mercaptoleGlycerinTheobromine
Benzyl ureaHumuleneTrigonelline
Beta-sitosteryl arabinosyl vanilloyl stearateKaempferolα-thujones
BorneolLimoneneβ-cymene
Bornyl acetateLinaloolβ-myrcene
ButanediamideNaringenineβ-santatol
CaffeineN-benzyl- 2-phenylacetamideβ-thujones
CamphorN-benzyl-benzamide
CaryophylleneOctacosanol
Table
Table 9. List of parameters used to generate grid on target proteins for molecular docking.
Table 9. List of parameters used to generate grid on target proteins for molecular docking.
Gene NameNo. of Grid Points in Axes (X, Y, Z)No. of Grid Points in ÅGrid Center Coordinates (X, Y, Z)
CiaR260, 160, 3000.37549.202, 21.493, −15.155
LepC150, 250, 2000.37532.025, 56.24, 14.572
OppC150, 100, 1200.37567.03, 139.688, 246.646
SecA300, 300, 3000.375−7.808, 13.046, 09.175
SMU1784c100, 100, 1200.375−30.648, 25.537, 8.584
SMU659 (SpaR)250, 260, 3000.375−56.159, 24.702, −77.258
Table
Table 10. Details of Protein–Ligand Complex (PLC).
Table 10. Details of Protein–Ligand Complex (PLC).
Protein NameCiaRLepCOppCSecASMU1784cSpaRYidc2
Total Residues4461304975883434219
Protein Chain(s)A, BACAAA, BA
Residue in Chain(s)223, 2231304975883217, 217219
No. of Atoms7188210474612072132471223509
No. of Heavy Atoms35721043356605963335201710
Atoms Charge−249−1−241449
Ligand nameCampesterolMeliantrolStigmasterolCitrostadienolCitrostadienolIsofucosterolUrsolic acid
No. of Atoms77 (total) 29 (heavy)79 (total) 33 (heavy)78 (total) 30 (heavy)81 (total) 31 (heavy)81 (total) 31 (heavy)78 (total) 30 (heavy)81 (total) 33 (heavy)
Atomic Mass [a.u.]400.694462.676412.705426.732426.732412.705456.715
Molecular FormulaC28H48OC28H46O5C29H48OC30H50OC30H50OC29H48OC30H48O3
No. of Fragments1122221
No. of Rotatable Bonds6766663
Table
Table 11. Details of counterions and salt concentration added for simulation.
Table 11. Details of counterions and salt concentration added for simulation.
Protein–Ligand ComplexTypeNo. of AtomsConcentration [mM]Total Charge
CiaR–campesterolNa9268.69892
Cl6850.777−68
LepC–meliantrolNa3250.49232
Cl4164.693−41
OppC–stigmasterolNa14.2831
Cl---
SecA–citrostadienolNa10166.292101
Cl7750.54−77
SMU1784c–citrostadienolNa1250.81112
Cl26110.09−26
SpaR–isofucosterolNa5550.91655
Cl5954.62−59
Yidc2–ursolic acidNa---
Cl915.519−9
Table
Table 12. Details of default parameters set before simulation.
Table 12. Details of default parameters set before simulation.
Protein–Ligand ComplexEnsembleTemperature [K]Simulation Time [ns]No. of AtomsNo. of Water MoleculesCharge
CiaR_campesterolNPT310.1100.10280,47224,3490
LepC_meliantrolNPT310.110136,82511,5230
OppC_stigmasterolNPT310.1100.10213,56042450
SecA_citrostadienolNPT310.110195,43427,7010
SMU1784c_citrostadienolNPT310.110114,32542940
SpaR_isofucosterolNPT310.110166,23419,6400
Yidc2_ursolic acidNPT310.1100.10235,23110,5440
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marimuthu, S.C.V.; Murugesan, J.; Babkiewicz, E.; Maszczyk, P.; Sankaranarayanan, M.; Thangamariappan, E.; Rosy, J.C.; Ram Kumar Pandian, S.; Kunjiappan, S.; Balakrishnan, V.; et al. Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans. Molecules 2023, 28, 5514. https://doi.org/10.3390/molecules28145514

AMA Style

Marimuthu SCV, Murugesan J, Babkiewicz E, Maszczyk P, Sankaranarayanan M, Thangamariappan E, Rosy JC, Ram Kumar Pandian S, Kunjiappan S, Balakrishnan V, et al. Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans. Molecules. 2023; 28(14):5514. https://doi.org/10.3390/molecules28145514

Chicago/Turabian Style

Marimuthu, Shakti Chandra Vadhana, Jayaprabhakaran Murugesan, Ewa Babkiewicz, Piotr Maszczyk, Murugesan Sankaranarayanan, Esakkimuthu Thangamariappan, Joseph Christina Rosy, Sureshbabu Ram Kumar Pandian, Selvaraj Kunjiappan, Vanavil Balakrishnan, and et al. 2023. "Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans" Molecules 28, no. 14: 5514. https://doi.org/10.3390/molecules28145514

APA Style

Marimuthu, S. C. V., Murugesan, J., Babkiewicz, E., Maszczyk, P., Sankaranarayanan, M., Thangamariappan, E., Rosy, J. C., Ram Kumar Pandian, S., Kunjiappan, S., Balakrishnan, V., & Sundar, K. (2023). Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans. Molecules, 28(14), 5514. https://doi.org/10.3390/molecules28145514

Article Metrics

Back to TopTop