- Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
Introduction: Ever since the outbreak of listeriosis and other related illnesses caused by the dreadful pathogen Listeria monocytogenes, the lives of immunocompromised individuals have been at risk.
Objectives and Methods: The main goal of this study is to comprehend the potential of terpenes, a major class of secondary metabolites in inhibiting one of the disease-causing protein Internalin A (InlA) of the pathogen via in silico approaches.
Results: The best binding affinity value of −9.5 kcal/mol was observed for Bipinnatin and Epispongiadiol according to the molecular docking studies. The compounds were further subjected to ADMET and biological activity estimation which confirmed their good pharmacokinetic properties and antibacterial activity.
Discussion: Molecular dynamic simulation for a timescale of 100 ns finally revealed Epispongiadiol to be a promising drug-like compound that could possibly pave the way to the treatment of this disease.
1 Introduction
Plants have been an ideal source of medicine since prehistoric times due to their beneficial effects on human health with minimum or no side effects. For this reason, medicinal plants and other plant-derived products are used extensively for treating various infectious diseases and ailments by the pharmaceutical industries. Medicinal plants are considered a ‘storehouse’ of an array of biologically active compounds with varied therapeutic applications (Pagare et al., 2015; Crozier et al., 2006; Tiwari and Rana, 2015; Yang et al., 2018). These plants possess anticancer, antibacterial, antiviral, analgesic and anti-inflammatory properties among their many other medicinal uses. More than 80% of the world’s developing population relies on traditional medical practices involving the usage of plants for their primary wellness as reported by the World Health Organization (WHO) (Bennett and Wallsgrove, 1994; Teoh and Teoh, 2016). The advancements in modern research has facilitated a significant surge in the study of plant metabolites throughout the last century. Secondary metabolites (SMs) are vital compounds that are typically produced by plants and have a distinct carbon skeleton framework. These metabolites are the end products of primary metabolites generated from the biosynthetic modifications such as methylation, glycosylation and hydroxylation. Although not directly involved in respiratory and photosynthetic metabolism, these metabolites assist plants in withstanding specific environmental conditions. In fact, the SMs are essential in defending the plant against potential threats like animals and various microbes. Moreover, certain plants utilize SMs to attract seed dispensers and pollinators, as well as signals to mediate the symbiotic relationship between plants and microorganisms (Aye et al., 2019; Twaij and Hasan, 2022). The side chains and structural composition certainly make SMs more complex. SMs can be categorized based on factors including their solubility in different solvents, chemical structure, biosynthetic pathways and composition. Terpenes (like sterols, carotenoids, plant volatiles, cardiac glycosides), phenolics (like coumarins, flavonoids, phenolic acids, lignin, tannins, stilbenes, lignans), nitrogen containing (like non-protein amino acids, alkaloids, cyanogenic glucosides) and sulfur containing compounds (like thionins, glutathione, phytoalexins, defensins) are some of the main groups belonging to SMs (Roba, 2020; Theis and Lerdau, 2003; Tetali, 2019).
Among the different phytochemicals, terpenes are the largest class of natural products with over 30,000 members that have been utilized for a wide range of applications such as flavoring, perfume, medicine, cosmetics, biofuels, food and beverages. They are the most diverse family that range in both structure and size. These volatile, unsaturated five-carbon cyclic molecules are composed of isoprene units, and the type of terpene that is formed varies depending on the number of isoprene units present (Costa et al., 2012). It is the hemiterpenes (1 isoprene unit) - dimethylallyl pyrophosphate (DMAPP) and isopentenyl pyrophosphate (IPP) that give rise to the many subclasses of terpenes. Terpene subclasses include mono- (2 isoprene units), sesqui- (3 isoprene units), Di- (4 isoprene units), sester- (5 isoprene units), tri- (6 isoprene units) and carotenoids (8 isoprene units). Numerous studies have elucidated the ability of terpenes to improve skin penetration, prevent inflammatory disorders, shield living organisms from biotic and abiotic stresses and also their antimicrobial property to combat infectious diseases (Croteau et al., 2000; Boncan et al., 2020; Cox-Georgian et al., 2019). The possible antibacterial potential of terpenes will be further explored in this work.
The food and economic sectors are being severely disrupted by the outbreak of foodborne illnesses and the subsequent death rates. Thus, public health initiatives have primarily focused on widely recognized foodborne diseases and pathogens in the food chain as a response to this growing concern (Altekruse et al., 1997; Motarjemi and Kiiferstein, 1997; Newell et al., 2010). Multiple factors contribute to the rise in foodborne diseases which include increase in immunocompromised individuals, lack of proper microbiological safety regulations before global food trading, transportation conditions that allow the pathogens to survive on the food and get to the consumers in a viable state, consumption of raw vegetables, meat and unpasteurized dairy products as well as temperature related changes favorable to the pathogen growth and production (Tauxe, 1997; Thakur et al., 2018). Listeriosis caused by Listeria monocytogenes (L. monocytogenes) is a systemic disease that arises from consuming contaminated food, notably ready-to-eat food and can be fatal to those with weak immune systems. Pregnant women are more vulnerable to this disease than the elderly people or infants, as several reports have stated 16%–17% of Listeria monocytogenes infections in this population (Matereke and Okoh, 2020; Ireton et al., 2021; Lecuit et al., 1997). Since this Gram-positive, psychrotolerant pathogen can pass through the blood brain and placental barriers, it often causes spontaneous miscarriages and meningoencephalitis in pregnant women. The pathogen’s tolerance to environmental stresses and sublethal doses of antibiotics or other antimicrobial drugs are significant factors to its antimicrobial resistance. Despite the advancements in the medical field to fight L. monocytogenes, the bacteria still pose a hazard to food safety as proved by the listeriosis outbreak that occurred in South Africa in 2017–2018. Of all the foodborne diseases, listeriosis has the highest death rates often ranging between 20% and 30% and sometimes even exceeds 50% (Bonazzi et al., 2009a; b; Lebrun et al., 1996). This pathogen expresses a wide range of virulence factors in order to proliferate and survive in the gastrointestinal system as well as to pass through the various host biological barriers. Internalins are one of the major virulence factors of L. monocytogenes that aids the pathogen to adhere and invade the host cell membrane. Among the internalins, the most researched ones are internalin A (InlA) and internalin B that play a vital role in its pathogenesis. InlA is an 88 kDa protein that comprises of a signal peptide at the N-terminus followed by 15 leucine rich repeats (LRR). An important structural component of InlA is an inter-repeat (IR) domain located downstream of the LRR domain that is essential for the binding of LRR region to its host cell receptor E-cadherin (Martin, 2003). This widely expressed and evenly distributed surface protein is covalently attached to the cell wall and upon binding to the E-cadherin receptor facilitates the pathogen’s internalization into the enterocytes and entry to the primary site of the host cell, i.e., the intestinal epithelial barrier. An LPTXG motif that anchors InlA to the bacterial cell wall is located at the C-terminal end of the protein followed by a sorting peptide (Lopes-Luz et al., 2021; Kathariou, 2002).
To date, no in silico investigations have been performed to predict the potential of terpenes in suppressing the activity of InlA protein. Thus, this study was intended to investigate the binding interactions of InlA with 80 terpenes through molecular docking and dynamic simulation approaches with an ultimate goal of identifying potential inhibitors of the virulence protein.
2 Materials and methods
2.1 Preparation of protein and ligands
The 3D structure of InlA protein from L. monocytogenes with resolution of 1.60 Å (PDB ID: 1O6T) (Figure 1) was obtained in PDB format from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (https://www.rcsb.org/, accessed 12 February 2024) (Berman et al., 2000). This protein was then subjected to preparation prior to docking using the AutoDock Tools 1.5.6 by eliminating water molecules, heteroatoms and adding polar hydrogens and Kollman charges. Finally, the prepared target was saved in pdbqt format (Morris et al., 2009).
In the case of ligands, 80 terpenes from different plant species with a range of medicinal properties were selected based on ‘Lipinski’s rule of five’ (Table 1) and the 3D structures of these compounds were obtained in SDF format from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/, accessed12 February 2024) (Benet et al., 2016; Kim, 2021). The compounds were energy minimized using the steepest descent algorithm of Avogadro software and Open Babel program was further utilized to generate their pdbqt files (Hanwell et al., 2012; O'Boyle et al., 2011).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-t001.jpg)
Table 1. Selected terpenes that are Lipinski compliant. The abbreviations HBD, HBA, Log P, RB and TPSA stands for Hydrogen bond donor, hydrogen bond acceptor, octanol-water partition coefficient, rotatable bonds and topological polar surface area.
2.2 Molecular docking of InlA with the selected terpenes
AutoDock Vina program was used to dock all the 80 ligands against InlA protein inorder to gain an insight into the various interactions involved during the binding process (Trott and Olson, 2010). This program employs an efficient optimization technique that depends on a special scoring function (combination of empirical and knowledge-based approaches) and a gradient-based local search genetic algorithm to predict the binding modes thereby producing reliable information concerning the receptor and ligand interactions (Jaghoori et al., 2016; Sarkar et al., 2024). The binding affinities were estimated for each InlA-ligand complex. The protein was kept rigid and the ligands were flexible throughout the docking study. With the grid box centered at x = −15.289, y = −1.522 and z = 26.283 and dimensions set at points 126 × 82 × 126 separated by 1 Å, all the input files needed for docking were created in ADT 1.5.6. Based on the best docking scores, the top ten compounds from among the 80 terpenes were chosen for further interaction analysis. A compound with greater negative binding affinity value revealed its stronger binding to the target receptor (Rolta et al., 2022).
2.3 Visualization of binding interactions
Different interactions and the amino acid residues involved during the binding of InlA protein with the top 10 compounds were examined using the BIOVIA discovery studio (DS) visualizer and Ligplot+ software (Studio, 2008; Laskowski and Swindells, 2011). In particular, Ligplot+ was used to study the two-dimensional (2D) protein-ligand interactions involving hydrogen and hydrophobic interactions, and DS visualizer was used to investigate the different forms of hydrogen and hydrophobic bonding that were present in the complexes. Additionally, PyMOL software was also utilized for the three-dimensional (3D) visualization of receptor-ligand complexes (Yuan et al., 2017).
2.4 ADMET and biological activity prediction
The pharmacokinetic properties and the efficacy of the top 10 compounds obtained upon docking were computed by the admetSAR server (http://lmmd.ecust.edu.cn/admetsar2/, accessed 15 February 2024). The server predicts the carcinogenicity, blood-brain barrier penetration (BBB), Human Intestinal Absorption (HIA), subcellular localization, human oral bioavailability, toxicity and various other parameters of drugs or drug-like compounds that play a crucial role in the drug development process (Cheng et al., 2012). Additionally, the PASS (Prediction of Activity Spectra for Substances) online server (https://www.way2drug.com/passonline/, accessed 18 February 2024) was used to further explore the possible therapeutic activities of these compounds (Filimonov et al., 2014). This server estimates a compound’s expected activity spectrum as ‘probability of activity (Pa)’ and ‘probability of inactivity (Pi)’. Pa and Pi have values ranging from 0.000 to 1.000 and the compound activity is considered possible only if Pa > Pi (Chakraborty et al., 2016).
2.5 Molecular dynamic (MD) simulation and MM/GBSA binding free energy estimation
MD simulation of protein-ligand complexes with the best docking scores and the protein alone (InlA) was carried out using the GROMACS software (version 2023) for a 100 nanosecond (ns) timescale to gain an insight into their stability as well as the conformational behaviour upon ligand binding (Bauer et al., 2022). CHARMM27 forcefield of GROMACS and SwissParam server were used to generate the protein and ligand topologies (Bjelkmar et al., 2010; Bugnon et al., 2023). Both the protein and complex systems were further solvated in cubic box and neutralized by adding Na+/Cl− ions. Post solvation and neutralization, using steepest descent algorithm the systems were energy minimized to eliminate the steric clashes and then equilibrated utilizing the NVT and NPT ensembles so as to maintain the temperature (300K) and pressure (1atm) of the systems (Berendsen et al., 1984; Andersen, 1980; Petersen, 1995). Ultimately, MD run was initiated for 100 ns and the root mean square deviation (RMSD), root mean square fluctuation (rmsf), radius of gyration (Rg), solvent accessible surface area (SASA) of the protein-ligand and protein systems were determined. Furthermore, the protein-ligand complexes were subjected to principal component analysis (PCA). Xmgrace software was used to obtain the 2D plots of all the aforementioned parameters (Turner, 2005).
Post MD simulation, MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) binding free energy of the protein-ligand complexes with the best docking scores were estimated for the last 20 ns MD trajectory using the gmx_MMPBSA program (v1.6.0) (Valdés-Tresanco et al., 2021). In MM/GBSA method, the total binding free energy is determined by the following equation:
where Gcomplex is the free energy of the protein-ligand complexes, while Gprotein and Gligand are the free energies of the protein and ligands. The conformational entropy upon ligand binding is represented by T∆S. ∆Ggas is the sum of internal energies (ΔEint - bond, angle, dihedral), van der Waals energy (ΔEvdW) and electrostatic component of internal energy (ΔEele). ∆Gsol is the combination of polar (∆EGB) and non-polar (∆ESURF) components of solvation energy.
3 Results
3.1 Lipinski rule evaluation of terpenes
All the selected compounds obeyed ‘Lipinski’s rule of 5’ without any violation suggesting their possibility of oral administration. According to this rule, a drug or a drug-like compound can be orally administered only if it satisfies the following criteria: molecular weight ≤500 g/mol, HBD ≤5, HBA ≤10 and LogP ≤ 5 (Walters, 2012; Saini et al., 2021). The molecular weight of all the terpenes varied from 136.23 to 398.8 g/mol. Alpha-Fenchene and Isoterpinolene had the lowest molecular weight of 136.23 g/mol while the highest value of 398.8 g/mol was shown by Chlorojanerin. Moreover, the HBDs and the HBAs of all the compounds ranged from 0 to 3 and 0 to 7, respectively. The lowest logP value of 0.1 was noted for Janerin whereas Cespitularin A had the maximum logP value of 4.4. These observations collectively demonstrate the potential of the chosen terpenes in oral delivery.
3.2 Molecular docking analysis of protein with the terpenes
All the 80 terpenes were blindly docked against the InlA protein with the grid box covering the entire protein surface. Additionally, a reference ligand Ampicillin, which is one of the most widely used antibiotic to treat illnesses caused by L. monocytogenes, was also docked against the virulence protein in order to determine whether the terpenes were more or less effective compared to the drug. We observed exceptional binding affinity values for each protein-ligand complexes and the top 10 compounds (Figure 2) with the highest scores were finalised for further analysis. The important amino acid residues and the different interactions involved in the binding of the best 10 ligands with the InlA protein from DS visualizer and Ligplot+ software are displayed in table 2. The binding affinity values obtained for all the 80 terpenes and Ampicillin are highlighted in table 3.
According to numerous studies, higher negative binding affinity values often imply a ligand’s increased capacity to inhibit the disease-causing target protein (Deepasree and Subhashree, 2023; Sahu and Pattanayak, 2019; Islam et al., 2021; Correa-Basurto et al., 2015). The top 10 terpenes had varied binding affinities ranging from −8.8 to −9.5 kcal/mol. Among all the selected terpenes, Bipinnatin (lig.13) and Epispongiadiol (lig.54) were found to have displayed the best interaction with InlA protein with a binding affinity value of −9.5 kcal/mol followed by Artemisinin and Bonducellpin G showing the next highest docking score of −9.1 kcal/mol. In DS visualizer, Bipinnatin displayed 2 hydrogen bonds with Ser173 and Asp213 amino acid residues of InlA along with 3 hydrophobic interactions (2 pi-alkyl and 1 pi-sigma bonds) with Phe150 (Figure 3A). Similarly, the ligand13-InlA interaction analysis with Ligplot+ revealed a maximum of 3 hydrogen bonds to Ser172, Ser173 and Asp213 residues formed at a distance of 2.95 Å, 3.02 Å and 3.23 Å. Here, Asn129, Phe150 and Asn151 had hydrophobic interactions with Bipinnatin (Figure 3B). The 3D representation of ligand13-InlA complex is shown in Figure 3C. In the case of Epispongiadiol, the ligand showed 1 conventional and 1 pi-donor hydrogen bond with Asn129, 1 carbon hydrogen bond with Ser172 and 3 pi-alkyl bonds with Phe150 in DS visualizer (Figure 4A) while 2 hydrogen bonds to Asn129 and Glu 170 (formed at a distance of 3.29 Å and 3.04 Å) were noted in Ligplot+ (Figure 4B). Moreover, the compound had hydrophobic interactions with Asn151, Phe150 and Ser172 residues of the virulence protein. Figure 4C depicts the 3D view of ligand54-InlA complex.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g003.jpg)
Figure 3. (A) and (B) 2D images of Bipinnatin interacting with InlA protein observed in discovery studio and Ligplot+. Green dotted lines indicate the hydrogen bonds involved while the red pockets containing amino acid residues reveal the residues involved in hydrophobic interactions (C) 3D view of Bipinnatin-InlA complex. Yellow dotted lines represent the hydrogen bonds formed during the binding process.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g004.jpg)
Figure 4. (A) and (B) 2D images of Epispongiadiol interacting with InlA protein observed in discovery studio and Ligplot+. (C) 3D view of Epispongiadiol-InlA complex.
Artemisinin (ligand 27) that exhibited the next best affinity value (−9.1 kcal/mol) had 1 hydrogen bond with Ser172, 1 pi-alkyl and 1 pi-sigma bond with Phe150 residue of the protein (DS visualizer- Figure 5A). No hydrogen bonds were revealed upon the visualization of ligand 27-InlA complex in Ligplot+ and only hydrophobic interactions to Ser173, Ser172, Phe150, Glu170 and Asp213 were identified (Figure 5B). Also, Bonducellpin G (ligand 75) sharing the same value as that of ligand 27 possessed 5 hydrogen bonds with residues including Asn151, Asn129, Ser173, Ser215 and Asp213 as well as 1 pi-alkyl bonding with Phe150. The ligand also had an additional pi-pi stacked bonding with Phe150 residue (DS visualizer- Figure 6A). Figure 6B (Ligplot+) shows that ligand 75 formed hydrogen bonds with Asn129 and Ser173 each at a distance of 3.15 Å and 3.23 Å. Moreover, this compound had hydrophobic interactions with residues, namely, Ser215, Ser216, Ser172, Phe150 and Glu170.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g005.jpg)
Figure 5. (A) and (B) 2D images of Artemisinin interacting with InlA protein observed in discovery studio and Ligplot+.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g006.jpg)
Figure 6. (A) and (B) 2D images of Bonducellpin G interacting with InlA protein observed in discovery studio and Ligplot+.
Ligand 9 (Inuchinenolide C), ligand 39 (Arbusculin C), and ligand 78 (Scabequinone) had a docking score of −9.0 kcal/mol respectively. Out of the best 10 compounds, Axinysone C (ligand 44) showed the lowest binding affinity value of −8.8 kcal/mol with InlA protein indicating its weak interaction with the target protein when compared to the other compounds in the top 10 list. In contrast to the top 10 terpenes with the best docking scores, Ampicillin displayed a remarkably low binding affinity value of −8.7 kcal/mol. This result reveals that the top 10 terpenes exhibited stronger interaction with InlA in comparison to Ampicillin, pointing to the terpenes’ potential to be more bactericidal than the drug. Figures 7A, B (2D images) and Figure 7C (3D images) depict the various interactions involved between Ampicillin and InlA. Therefore, based on the docking results, Bipinnatin and Epispongiadiol with the highest docking scores could possibly act as potential antibacterial agents that can suppress the virulence protein InlA of L. monocytogenes. Furthermore, the key residues that were observed predominantly in the docked complexes were Ser172, Ser173, Phe150 and Asp213. This may also highlight the vital role played by these residues in inhibiting the protein upon binding with the antibacterial drugs thereby turning the virulence protein inactive.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g007.jpg)
Figure 7. (A) and (B) 2D images of Ampicillin interacting with InlA protein observed in discovery studio and Ligplot+. (C) 3D view of Ampicillin-InlA complex.
3.3 ADMET and biological activity analysis
Post docking, the 10 terpenes were subjected to ADMET and biological activity analysis for a deeper understanding of their pharmacokinetic profiles and wide spectrum of therapeutic functions particularly the antibacterial activity. Some of the significant pharmacokinetic characteristics and therapeutic activities observed for the top 10 terpenes are shown in tables 4, 5.
One of the biggest challenges faced during oral drug development is assessing the drug’s capability to cross the intestinal epithelial barrier that estimates the rate and extent to which the drug is absorbed into the human body eventually impacting its bioavailability (Rudrapal et al., 2022; Aja et al., 2021; Qaddir et al., 2017). For this reason, both HIA and Caco-2 intestinal cell permeability of the compounds were estimated. Except for Tanacetin and Arbusculin C, the remaining terpenes exhibited positive response for HIA indicating that each of those terpenes can undergo intestinal absorption upon oral administration. Also, Inuchinenolide C and Bonducellpin G did not show caco-2 permeability while all other terpenes had positive response for the same. P-glycoprotein, an ABC transporter is a protein that actively engages in a number of biological processes including drug absorption, metabolism, distribution and excretion thereby detoxifying and protecting the body from harmful substances. It was discovered that none of the 10 compounds inhibited this protein. With the exception of Tanacetin and Arbusculin C, the remaining terpenes revealed BBB penetration. Moreover, all the terpenes were observed to be non-carcinogens. The subcellular localization of the terpenes were found to be in the mitochondria except for Furodysinin in plasma membrane. In the case of acute oral toxicity, the terpenes were categorized into four levels-mainly category I, II, III and IV based on their LD50 values. Category I consists of compounds with LD50 ≤ 50 mg/kg while category II has compounds with LD50 > 50 mg/kg but less than 500 mg/kg. Compounds with LD50 values higher than 500 mg/kg but less than 5,000 mg/kg are grouped into category III and those with LD50 value greater than 5,000 mg/kg comes under category IV. Generally, compounds that fall in category III and IV are often considered non-toxic while II is slightly toxic and category I compounds are highly toxic (Guan et al., 2019). In this study, none of the terpenes belonged to category I for acute oral toxicity. However, Inuchinenolide C and Bipinnatin displayed category II toxicity indicating their slightly toxic nature. Aside from Artemisinin and Tanacetin that belonged to category IV, all other terpenes were placed in category III. Thus, the top 10 compounds were predicted to have satisfactory ADMET findings for the major parameters.
The biological activities possessed by the 10 terpenes were further explored using the PASS web server. Antineoplastic, antiprotozoal, immunosuppressant, antiinflammatory, antifungal, antiviral and antileukemic activities were some of the predominantly expressed functions of the 10 terpene compounds obtained from the server and the Pa value was significantly higher than Pi in all the activity predictions. Furthermore, all the terpenes except for Inuchinenolide C, Furodysinin and Axinysone C exhibited antibacterial activity supporting the core aim of the study.
3.4 MD simulation analysis
Bipinnatin (ligand 13) and Epispongiadiol (ligand 54) showcased the highest binding affinities (−9.5 kcal/mol) based on the docking data obtained and hence MD simulation was performed for the complexes of these compounds with InlA protein as well as the apo form of InlA over a 100 ns timescale. Here, the RMSD, RMSF, Rg and SASA of the complexes and the apo protein were computed as a function of time inorder to monitor the variations in their stability, conformation and the amino acid residue interactions (Joshi et al., 2021; Surti et al., 2020; Nagamalla et al., 2022).
3.4.1 RMSD analysis
The RMSD was assessed to understand the stability of both the complex systems and apo protein that provides information on the degree of protein deviation from its native conformation upon ligand binding (Mani et al., 2023; Rampogu et al., 2022). The backbone RMSD plots of ligand13-InlA, ligand54-InlA and apo form are depicted in Figure 8. The apo protein started off with an RMSD value of 0.25 nm and then attained stability at 0.3 nm after a time period of 65 ns. This value was maintained till the end of 100 ns MD run. The Bipinnatin-InlA complex (complex13), displayed stability upto 20 ns at 0.3 nm RMSD but then significant deviations were seen until 70 ns at 0.45 nm. After 70 ns, the complex regained stability at the initial backbone RMSD value of 0.3 nm and remained stable with the same RMSD value until the 100 ns timescale. On the other hand, Epispongiadiol-InlA complex (complex54) was found to be stable after 50 ns maintaining a backbone RMSD value of 0.3 nm. A slight decrease in RMSD was noticed at 90 ns. Thus, ligand54- InlA complex was observed to be more stable than ligand13-InlA complex and apo protein according to the obtained RMSD results.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g008.jpg)
Figure 8. RMSD plots of Apo protein (black), Bipinnatin-InlA complex (red) and Epispongiadiol-InlA complex (green).
3.4.2 RMSF analysis
To gain a better understanding of the fluctuations in protein atoms or residues resulting from ligand interaction, RMSF analysis was carried out for all the 3 systems: apo protein, Bipinnatin-InlA and Epispongiadiol- InlA complexes. The presence of more flexible regions (turns, loops) are indicated by higher backbone RMSF values while lower value confirms the presence of secondary structures (helices or sheets) (Bhatt et al., 2021; Shahlaei et al., 2011; Krushna et al., 2017). The initial RMSF of apo protein was recorded to be 0.29 nm and as the 100 ns run proceeded, maximum fluctuations were observed at 0.33 nm between residues 420–440 and 450–480. Similarly, the Bipinnatin-InlA complex also began with an RMSF value of 0.29 nm and higher fluctuations were visible at 0.3 nm between residues 490 and 500. Minimum fluctuations with an average RMSF value of 0.2 nm were noted between the start and end of the 100 ns trajectory. The starting RMSF value of Epispongiadiol-InlA complex was noted to be 0.31 nm and higher variations were observed between residues 450–480 at 0.3 nm. Hence, from the above results it can be inferred that both ligand13-InlA and ligand54-InlA complexes did not exhibit much fluctuations from the apo protein (Figure 9) indicating the presence of more secondary structures than flexible regions in all the systems. This can further reveal that the protein has undergone less structural changes as a result of ligand binding.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g009.jpg)
Figure 9. RMSF plots of Apo protein (black), Bipinnatin-InlA complex (red) and Epispongiadiol-InlA complex (green).
3.4.3 Rg analysis
The radius of gyration (Rg) determined the structural compactness of all the systems under study by measuring the atomic distribution around the center of mass. In other words, Rg examines if the system steadily folds or unfolds during the specific simulation time period (Sharma et al., 2020; Thirumal Kumar et al., 2017; Apicella et al., 2017). Rg plots generated for apo protein, Bipinnatin-InlA and Epispongiadiol-InlA complexes are displayed in Figure 10. It was discovered that both the protein and protein-ligand complexes had relatively same and consistent Rg values over the course of 100 ns trajectory. The average Rg value for all the systems were noted to be 3.05 nm. Therefore, this finding implies that the protein-ligand complexes and apo protein were perfectly compact and in folded condition during the entire 100 ns timescale.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g010.jpg)
Figure 10. Rg plots of Apo protein (black), Bipinnatin-InlA complex (red) and Epispongiadiol-InlA complex (green).
3.4.4 SASA analysis
Solvent accessible surface area (SASA) is a parameter that estimates the area of protein’s surface exposed to solvent. Usually, higher SASA indicate greater exposure of protein residues to water which in turn lead to smaller binding interface area. Hence, a lower SASA value is expected during the MD run which indicates less protein expansion and better stability of the protein structure (Rahman et al., 2021; Saini et al., 2019). Initially, both apo protein and the protein-ligand complexes (ligand 13-InlA and ligand 54-InlA) displayed similar SASA values of 210 nm2. As the 100 ns simulation progressed, slight variations in the values were observed for each of the three systems. By the end of 100 ns, both the apo protein and Epispongiadiol-InlA complex had the same SASA value of 205 nm2 while the Bipinnatin-InlA complex exhibited a marginally higher value of 215 nm2. Thus, the above results suggest that upon the binding of Bipinnatin, the InlA protein had more exposure to solvent when compared to Epispongiadiol. Moreover, the binding of Epispongiadiol had least affected the stability of InlA protein as both the apo form and ligand 54-InlA complex showcased similar SASA values by the end of 100 ns trajectory (Figure 11).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g011.jpg)
Figure 11. SASA plots of Apo protein (black), Bipinnatin-InlA complex (red) and Epispongiadiol-InlA complex (green).
3.4.5 Principal component analysis
The collective motion of Bipinnatin-InlA and Epispongiadiol-InlA complexes were predicted by PCA. In other words, PCA determines the structural changes that occur in the protein upon ligand binding (Maisuradze et al., 2009). The dynamics of the 2 protein-ligand complexes were estimated with respect to the backbone. The initial principal components (PCs) help in the recognition of dominating motions and also in obtaining a significant portion of variance in the data. In order to understand the motive variations amongst the complex systems, the 2D projection of two eigenvectors (eigenvector 1 in X-axis and eigenvector 3 in the Y-axis) were assessed using the essential dynamics (ED) approach. A system is considered stable if it occupies a small phase space with stable clusters whereas systems that occupies more space and exhibits erratic clusters are categorized as less stable (Tong et al., 2021). Figure 12 indicates that both the complexes almost shared similar conformational spaces suggesting the stable nature of InlA protein upon the binding of Bipinnatin and Epispongiadiol.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1463750/fbinf-04-1463750-HTML/image_m/fbinf-04-1463750-g012.jpg)
Figure 12. 2D trajectory projection of Bipinnatin-InlA (red colour) and Epispongiadiol-InlA (green colour) complexes during MD simulation.
4 Binding free energy analysis
The total binding free energy of Bipinnatin-InlA and Epispongiadiol-InlA complexes were computed after 100 ns MD simulation. The overall binding free energy value of Bipinnatin-InlA complex was noted to be −12.12 kcal/mol which indicates the good stability and strong binding interactions between the protein and Bipinnatin (Figure 13). According to the data, the molecular mechanical energy changes in the gas phase (ΔGGas) displayed the highest negative energy value of −49.96 kcal/mol followed by the electrostatic component of internal energy (ΔEele) contributing to a value of −36.76 kcal/mol. In the case of Epispongiadiol-InlA complex, a total binding free energy value of −18.36 kcal/mol was observed which further suggests the strong binding and stable nature of Epispongiadiol-InlA complex (Figure 14). Here, the major contributor to the binding free energy was ΔGGas (−34.16 kcal/mol) and ΔEvdW (van der Waal’s energy) was the next best contributor with an energy value of −19.15 kcal/mol. Hence, based on the energy data, Epispongiadiol-InlA complex was more stable and had stronger interactions with the protein compared to Bipinnatin due their greater negative binding free energy value. Tables 6, 7 depicts the various energy contributions that resulted in the stable complexes of Bipinnatin and Epispongiadiol with InlA protein.
5 Discussion
Listeriosis caused by L. monocytogenes is not just a single disease rather a series of illnesses including meningitis, bacteremia, encephalitis or serious pregnancy related infections that can be fatal to those individuals with weakened immune system. According to an earlier study, the Centers for Disease Control and Prevention (CDC) estimated an annual occurrence of 1662 invasive listerial infections in the United States leading to 1520 hospitalizations and 266 deaths (Cartwright et al., 2013). Development of effective medications or antidotes against this disease-causing pathogen continue to be a significant challenge, despite enormous advancements in the medical sector. This failure can be attributed to the bacteria’s growing resistance to antibiotics. The bioactive compounds or phytochemicals found in plants have been explored by several scientists and researchers over the years in an effort to learn more about how these compounds might be able to treat some of the major illnesses caused by different microorganisms.
In the current work, we utilized in silico methods to investigate the possible antibacterial effect of terpenes on one of the virulence proteins of L. monocytogenes, InlA. Many research studies have reported a range of therapeutic properties of terpenes including anticancer, antiinflammatory, antiviral, antidiabetic, antiplasmodial and antimicrobial activities (Tetali, 2019). For example, researchers have discovered terpenes and its derivatives to be potential inhibitors of SARS-CoV-2 proteases as well as a promising natural antagonists of cancer through in silico studies (Diniz et al., 2021; Muhseen and Li, 2019). In our previously published work, the inhibitory effect of terpenoids (derived from terpenes) on another virulence protein of L. monocytogenes were predicted utilizing computational methods (Deepasree and Subhashree, 2023). Computational methodologies have proven to be very beneficial in pharmaceutical research because of their capacity to identify and generate novel, potent drugs, especially through the use of molecular docking and molecular dynamic simulation approaches. Every terpene chosen for this investigation complied with “Lipinski’s rule of five” and when these compounds were docked against the target protein InlA, it became evident that Bipinnatin and Epispongiadiol had the highest negative binding affinity values (−9.5 kcal/mol), which made them the most potent drug-like compounds. This work stands out for offering innovation as there are only very few reports linked to the aforementioned compounds. Selecting only the top ten compounds with the best docking scores, we proceeded to investigate the pharmacokinetic properties and also predicted the possible biological activities that these ten compounds might exhibit. The drugs chosen for analysis provided acceptable ADMET results for the compounds and also found that all the compounds with an exception of Inuchinenolide C, Furodysinin and Axinysone C displayed antibacterial activity which supports the goal of the study. Finally, the two compounds with the best docking scores were subjected to MD simulation for 100 ns and the results obtained suggests Epispongiadiol to be the best compound over Bipinnatin due to its greater structural stability. However, the aforementioned computational approaches have certain limitations as these methods don’t completely capture the cellular environments and other crucial interactions. In addition, the exact mechanism of action or resistance mechanisms may not be predicted by computational studies alone. Thus, preclinical and clinical investigations would be required to explore the in silico predictions in biological models for further therapeutic applications.
6 Conclusion
This study aimed to understand the potential of 80 terpenes in inhibiting one of the major virulence protein, Internalin A of the pathogen L. monocytogenes. Among the terpenes, Bipinnatin and Epispongiadiol were discovered to be the potential drug-like candidates that could be targeted against InlA due to their good binding affinity value of −9.5 kcal/mol. Although both the compounds possessed antibacterial activity upon biological activity evaluation, Bipinnatin belongs to a category of compounds that induces less toxicity only at lower doses according to ADMET analysis while Epispongiadiol is non-toxic even at higher doses. Results of the MD simulation for 100 ns revealed Epispongiadiol to be better than Bipinnatin due to their overall structural stability despite both the compounds having identical docking scores. The MM/GBSA binding free energy analysis confirmed the stable nature and stronger interaction of Epispongiadiol with the virulence protein. Hence, based on the in silico investigation, this work reports for the first time that Epispongiadiol could be a possible antibacterial agent that can be directed against InlA protein. Future studies are intended for the in vivo/in vitro validation with the best compound identified from this work.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
DK: Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing–original draft. SV: Conceptualization, Project administration, Supervision, Writing–review & editing.
Funding
The author(s) declare that financial support was received for the publication of this article from Vellore Institute of Technology.
Acknowledgments
The authors would like to place their sincere thanks to acknowledge the support offered by Vellore Institute of Technology towards the study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Aja, P. M., Agu, P. C., Ezeh, E. M., Awoke, J. N., Ogwoni, H. A., Deusdedit, T., et al. (2021). Prospect into therapeutic potentials of Moringa oleifera phytocompounds against cancer upsurge: de novo synthesis of test compounds, molecular docking, and ADMET studies. Bull. Natl. Res. Centre 45 (1), 99. doi:10.1186/s42269-021-00554-6
Altekruse, S. F., Cohen, M. L., and Swerdlow, D. L. (1997). Emerging foodborne diseases. Emerg. Infect. Dis. 3 (3), 285–293. doi:10.3201/eid0303.970304
Andersen, H. C. (1980). Molecular dynamics simulations at constant pressure and/or temperature. J. Chem. Phys. 72 (4), 2384–2393. doi:10.1063/1.439486
Apicella, A., Marascio, M., Colangelo, V., Soncini, M., Gautieri, A., and Plummer, C. J. (2017). Molecular dynamics simulations of the intrinsically disordered protein amelogenin. J. Biomol. Struct. Dyn. 35 (8), 1813–1823. doi:10.1080/07391102.2016.1196151
Aye, M. M., Aung, H. T., Sein, M. M., and Armijos, C. (2019). A review on the phytochemistry, medicinal properties and pharmacological activities of 15 selected Myanmar medicinal plants. Molecules 24 (2), 293. doi:10.3390/molecules24020293
Bauer, P., Hess, B., and Lindahl, E. (2022). GROMACS 2022.3. Manual, 673. doi:10.5281/zenodo.7037337
Benet, L. Z., Hosey, C. M., Ursu, O., and Oprea, T. I. (2016). BDDCS, the Rule of 5 and drugability. Adv. drug Deliv. Rev. 101, 89–98. doi:10.1016/j.addr.2016.05.007
Bennett, R. N., and Wallsgrove, R. M. (1994). Secondary metabolites in plant defence mechanisms. New phytol. 127 (4), 617–633. doi:10.1111/j.1469-8137.1994.tb02968.x
Berendsen, H. J., Postma, J. V., Van Gunsteren, W. F., DiNola, A. R. H. J., and Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81 (8), 3684–3690. doi:10.1063/1.448118
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein data bank. Nucleic acids Res. 28 (1), 235–242. doi:10.1093/nar/28.1.235
Bhatt, P., Joshi, T., Bhatt, K., Zhang, W., Huang, Y., and Chen, S. (2021). Binding interaction of glyphosate with glyphosate oxidoreductase and C–P lyase: molecular docking and molecular dynamics simulation studies. J. Hazard. Mater. 409, 124927. doi:10.1016/j.jhazmat.2020.124927
Bjelkmar, P., Larsson, P., Cuendet, M. A., Hess, B., and Lindahl, E. (2010). Implementation of the CHARMM force field in GROMACS: analysis of protein stability effects from correction maps, virtual interaction sites, and water models. J. Chem. theory Comput. 6 (2), 459–466. doi:10.1021/ct900549r
Bonazzi, M., Lecuit, M., and Cossart, P. (2009a). Listeria monocytogenes internalin and E-cadherin: from structure to pathogenesis. Cell. Microbiol. 11 (5), 693–702. doi:10.1111/j.1462-5822.2009.01293.x
Bonazzi, M., Lecuit, M., and Cossart, P. (2009b). Listeria monocytogenes internalin and E-cadherin: from bench to bedside. Cold Spring Harb. Perspect. Biol. 1 (4), a003087. doi:10.1101/cshperspect.a003087
Boncan, D. A. T., Tsang, S. S., Li, C., Lee, I. H., Lam, H. M., Chan, T. F., et al. (2020). Terpenes and terpenoids in plants: interactions with environment and insects. Int. J. Mol. Sci. 21 (19), 7382. doi:10.3390/ijms21197382
Bugnon, M., Goullieux, M., Rohrig, U. F., Perez, M. A., Daina, A., Michielin, O., et al. (2023). SwissParam 2023: a modern web-based tool for efficient small molecule parametrization. J. Chem. Inf. Model. 63 (21), 6469–6475. doi:10.1021/acs.jcim.3c01053
Cartwright, E. J., Jackson, K. A., Johnson, S. D., Graves, L. M., Silk, B. J., and Mahon, B. E. (2013). Listeriosis outbreaks and associated food vehicles, United States, 1998–2008. Emerg. Infect. Dis. 19 (1), 1–9. doi:10.3201/eid1901.120393
Chakraborty, J., Kabir, M. S. H., Banik, H., Majumder, S., Sayem, M. A., Karim, T., et al. (2016). In silico PASS prediction, molecular docking and ADME/T property analysis of isolated compounds from Tinospora cordifolia for new thrombolytic. drug Discov. doi:10.20959/wjpr20167-6632
Cheng, F., Li, W., Zhou, Y., Shen, J., Wu, Z., Liu, G., et al. (2012). admetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties. J. Chem. Inf. Model. 52, 3099–3105. doi:10.1021/ci300367a
Correa-Basurto, J., Cuevas-Hernández, R. I., Phillips-Farfán, B. V., Martínez-Archundia, M., Romo-Mancillas, A., Ramírez-Salinas, G. L., et al. (2015). Identification of the antiepileptic racetam binding site in the synaptic vesicle protein 2A by molecular dynamics and docking simulations. Front. Cell. Neurosci. 9, 125. doi:10.3389/fncel.2015.00125
Costa, T. D. S. A., Vieira, R. F., Bizzo, H. R., Silveira, D., and Gimenes, M. A., 2012. Secondary metabolites.
Cox-Georgian, D., Ramadoss, N., Dona, C., and Basu, C. (2019). Therapeutic and medicinal uses of terpenes. Med. plants farm Pharm., 333–359. doi:10.1007/978-3-030-31269-5_15
Croteau, R., Kutchan, T. M., and Lewis, N. G. (2000). Natural products (secondary metabolites). Biochem. Mol. Biol. plants 24, 1250–1319.
Crozier, A., Clifford, M. N., and Ashihara, H. (2006). Plant secondary metabolites. Occurrence, Structure and Role in the human diet. New Jersey, United States: Blackwell-Publishers. doi:10.1002/9780470988558
Deepasree, K., and Subhashree, V. (2023). Molecular docking and dynamic simulation studies of terpenoid compounds against phosphatidylinositol-specific phospholipase C from Listeria monocytogenes. Inf. Med. Unlocked 39, 101252. doi:10.1016/j.imu.2023.101252
Diniz, L. R. L., Perez-Castillo, Y., Elshabrawy, H. A., Filho, C. D. S. M. B., and de Sousa, D. P. (2021). Bioactive terpenes and their derivatives as potential SARS-CoV-2 proteases inhibitors from molecular modeling studies. Biomolecules 11 (1), 74. doi:10.3390/biom11010074
Filimonov, D. A., Lagunin, A. A., Gloriozova, T. A., Rudik, A. V., Druzhilovskii, D. S., Pogodin, P. V., et al. (2014). Prediction of the biological activity spectra of organic compounds using the PASS online web resource. Chem. Heterocycl. Compd. 50, 444–457. doi:10.1007/s10593-014-1496-1
Guan, L., Yang, H., Cai, Y., Sun, L., Di, P., Li, W., et al. (2019). ADMET-score–a comprehensive scoring function for evaluation of chemical drug-likeness. Medchemcomm 10 (1), 148–157. doi:10.1039/c8md00472b
Hanwell, M. D., Curtis, D. E., Lonie, D. C., Vandermeersch, T., Zurek, E., and Hutchison, G. R. (2012). Avogadro: an advanced semantic chemical editor, visualization, and analysis platform. J. cheminformatics 4, 17. doi:10.1186/1758-2946-4-17
Ireton, K., Mortuza, R., Gyanwali, G. C., Gianfelice, A., and Hussain, M. (2021). Role of internalin proteins in the pathogenesis of Listeria monocytogenes. Mol. Microbiol. 116 (6), 1407–1419. doi:10.1111/mmi.14836
Islam, R., Parves, M. R., Paul, A. S., Uddin, N., Rahman, M. S., Mamun, A. A., et al. (2021). A molecular modeling approach to identify effective antiviral phytochemicals against the main protease of SARS-CoV-2. J. Biomol. Struct. Dyn. 39 (9), 3213–3224. doi:10.1080/07391102.2020.1761883
Jaghoori, M. M., Bleijlevens, B., and Olabarriaga, S. D. (2016). 1001 Ways to run AutoDock Vina for virtual screening. J. computer-aided Mol. Des. 30, 237–249. doi:10.1007/s10822-016-9900-9
Joshi, T., Joshi, T., Sharma, P., Chandra, S., and Pande, V. (2021). Molecular docking and molecular dynamics simulation approach to screen natural compounds for inhibition of Xanthomonas oryzae pv. Oryzae by targeting peptide deformylase. J. Biomol. Struct. Dyn. 39 (3), 823–840. doi:10.1080/07391102.2020.1719200
Kathariou, S. (2002). Listeria monocytogenes virulence and pathogenicity, a food safety perspective. J. food Prot. 65 (11), 1811–1829. doi:10.4315/0362-028x-65.11.1811
Kim, S. (2021). Exploring chemical information in PubChem. Curr. Protoc. 1 (8), e217. doi:10.1002/cpz1.217
Krushna, G. S., Shivaranjani, V. L., Umamaheswari, J., Srinivasulu, C., Hussain, S. A., Kareem, M. A., et al. (2017). In vivo and molecular docking studies using whole extract and phytocompounds of Aegle marmelos fruit protective effects against isoproterenol-induced myocardial infarction in rats. Biomed. & Pharmacother. 91, 880–889. doi:10.1016/j.biopha.2017.04.115
Laskowski, R. A., and Swindells, M. B. (2011). LigPlot+: multiple ligand–protein interaction diagrams for drug discovery. J. Chem. Inf. Model. 51, 2778–2786. doi:10.1021/ci200227u
Lebrun, M., Mengaud, J., Ohayon, H., Nato, F., and Cossart, P. (1996). Internalin must be on the bacterial surface to mediate entry of Listeria monocytogenes into epithelial cells. Mol. Microbiol. 21 (3), 579–592. doi:10.1111/j.1365-2958.1996.tb02566.x
Lecuit, M., Ohayon, H., Braun, L., Mengaud, J., and Cossart, P. (1997). Internalin of Listeria monocytogenes with an intact leucine-rich repeat region is sufficient to promote internalization. Infect. Immun. 65 (12), 5309–5319. doi:10.1128/iai.65.12.5309-5319.1997
Lopes-Luz, L., Mendonça, M., Bernardes Fogaça, M., Kipnis, A., Bhunia, A. K., and Bührer-Sékula, S. (2021). Listeria monocytogenes: review of pathogenesis and virulence determinants-targeted immunological assays. Crit. Rev. Microbiol. 47 (5), 647–666. doi:10.1080/1040841X.2021.1911930
Maisuradze, G. G., Liwo, A., and Scheraga, H. A. (2009). Principal component analysis for protein folding dynamics. J. Mol. Biol. 385 (1), 312–329. doi:10.1016/j.jmb.2008.10.018
Mani, H., Chang, C. C., Hsu, H. J., Yang, C. H., Yen, J. H., and Liou, J. W. (2023). Comparison, analysis, and molecular dynamics simulations of structures of a viral protein modeled using various computational Tools. Bioengineering 10 (9), 1004. doi:10.3390/bioengineering10091004
Martin, S. E. (2003). Listeria monocytogenes. Microb. food Saf. animal Agric. Curr. Top., 257–266. doi:10.1002/9780470752616.ch26
Matereke, L. T., and Okoh, A. I. (2020). Listeria monocytogenes virulence, antimicrobial resistance and environmental persistence: a review. Pathogens 9 (7), 528. doi:10.3390/pathogens9070528
Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J. Comput. Chem. 30 (16), 2785–2791. doi:10.1002/jcc.21256
Motarjemi, Y., and Kiiferstein, F. K. (1997). Global estimation of food borne diseases. World health stat. Q. 50 (1-2), 5–11.
Muhseen, Z. T., and Li, G. (2019). Promising terpenes as natural antagonists of cancer: an in-silico approach. Molecules 25 (1), 155. doi:10.3390/molecules25010155
Nagamalla, L., Kumar, J. S., Shaik, M. R., Sanjay, C., Alsamhan, A. M., Kasim, M. A., et al. (2022). Identification of novel AXL kinase inhibitors using ligand-based pharmacophore screening and molecular dynamics simulations. Crystals 12 (8), 1158. doi:10.3390/cryst12081158
Newell, D. G., Koopmans, M., Verhoef, L., Duizer, E., Aidara-Kane, A., Sprong, H., et al. (2010). Food-borne diseases—the challenges of 20 years ago still persist while new ones continue to emerge. Int. J. food Microbiol. 139, S3–S15. doi:10.1016/j.ijfoodmicro.2010.01.021
O'Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., and Hutchison, G. R. (2011). Open Babel: an open chemical toolbox. J. cheminformatics 3, 1–14. Available at: http://www.jcheminf.com/content/3/1/33.
Pagare, S., Bhatia, M., Tripathi, N., Pagare, S., and Bansal, Y. K. (2015). Secondary metabolites of plants and their role: overview. Curr. trends Biotechnol. Pharm. 9 (3), 293–304.
Petersen, H. G. (1995). Accuracy and efficiency of the particle mesh Ewald method. J. Chem. Phys. 103 (9), 3668–3679. doi:10.1063/1.470043
Qaddir, I., Rasool, N., Hussain, W., and Mahmood, S. (2017). Computer-aided analysis of phytochemicals as potential dengue virus inhibitors based on molecular docking, ADMET and DFT studies. J. vector borne Dis. 54 (3), 255–262. doi:10.4103/0972-9062.217617
Rahman, M. M., Saha, T., Islam, K. J., Suman, R. H., Biswas, S., Rahat, E. U., et al. (2021). Virtual screening, molecular dynamics and structure–activity relationship studies to identify potent approved drugs for Covid-19 treatment. J. Biomol. Struct. Dyn. 39 (16), 6231–6241. doi:10.1080/07391102.2020.1794974
Rampogu, S., Lemuel, M. R., and Lee, K. W. (2022). Virtual screening, molecular docking, molecular dynamics simulations and free energy calculations to discover potential DDX3 inhibitors. Adv. Cancer Biology-Metastasis 4, 100022. doi:10.1016/j.adcanc.2021.100022
Rolta, R., Salaria, D., Kumar, V., Patel, C. N., Sourirajan, A., Baumler, D. J., et al. (2022). Molecular docking studies of phytocompounds of Rheum emodi Wall with proteins responsible for antibiotic resistance in bacterial and fungal pathogens: in silico approach to enhance the bio-availability of antibiotics. J. Biomol. Struct. Dyn. 40 (8), 3789–3803. doi:10.1080/07391102.2020.1850364
Rudrapal, M., Gogoi, N., Chetia, D., Khan, J., Banwas, S., Alshehri, B., et al. (2022). Repurposing of phytomedicine-derived bioactive compounds with promising anti-SARS-CoV-2 potential: molecular docking, MD simulation and drug-likeness/ADMET studies. Saudi J. Biol. Sci. 29 (4), 2432–2446. doi:10.1016/j.sjbs.2021.12.018
Sahu, S. N., and Pattanayak, S. K. (2019). Molecular docking and molecular dynamics simulation studies on PLCE1 encoded protein. J. Mol. Struct. 1198, 126936. doi:10.1016/j.molstruc.2019.126936
Saini, G., Dalal, V., Gupta, D. N., Sharma, N., Kumar, P., and Sharma, A. K. (2021). A molecular docking and dynamic approach to screen inhibitors against ZnuA1 of Candidatus Liberibacter asiaticus. Mol. Simul. 47 (6), 510–525. doi:10.1080/08927022.2021.1888948
Saini, G., Dalal, V., Savita, B. K., Sharma, N., Kumar, P., and Sharma, A. K. (2019). Molecular docking and dynamic approach to virtual screen inhibitors against Esbp of Candidatus Liberibacter asiaticus. J. Mol. Graph. Model. 92, 329–340. doi:10.1016/j.jmgm.2019.08.012
Sarkar, A., Concilio, S., Sessa, L., Marrafino, F., and Piotto, S. (2024). Advancements and novel approaches in modified autodock vina algorithms for enhanced molecular docking. Results Chem. 7, 101319. doi:10.1016/j.rechem.2024.101319
Shahlaei, M., Madadkar-Sobhani, A., Mahnam, K., Fassihi, A., Saghaie, L., and Mansourian, M. (2011). Homology modeling of human CCR5 and analysis of its binding properties through molecular docking and molecular dynamics simulation. Biochimica Biophysica Acta (BBA)-Biomembranes 1808 (3), 802–817. doi:10.1016/j.bbamem.2010.12.004
Sharma, V., Panwar, A., and Sharma, A. K. (2020). Molecular dynamic simulation study on chromones and flavonoids for the in silico designing of a potential ligand inhibiting mTOR pathway in breast cancer. Curr. Pharmacol. Rep. 6, 373–379. doi:10.1007/s40495-020-00246-1
Surti, M., Patel, M., Adnan, M., Moin, A., Ashraf, S. A., Siddiqui, A. J., et al. (2020). Ilimaquinone (marine sponge metabolite) as a novel inhibitor of SARS-CoV-2 key target proteins in comparison with suggested COVID-19 drugs: designing, docking and molecular dynamics simulation study. RSC Adv. 10 (62), 37707–37720. doi:10.1039/d0ra06379g
Tauxe, R. V. (1997). Emerging foodborne diseases: an evolving public health challenge. Emerg. Infect. Dis. 3 (4), 425–434. doi:10.3201/eid0304.970403
Teoh, E. S., and Teoh, E. S. (2016) “Secondary metabolites of plants,” in Medicinal orchids of asia, 59–73.
Tetali, S. D. (2019). Terpenes and isoprenoids: a wealth of compounds for global use. Planta 249, 1–8. doi:10.1007/s00425-018-3056-x
Thakur, M., Asrani, R. K., and Patial, V. (2018). “Listeria monocytogenes: a food-borne pathogen,” in Foodborne diseases (Academic Press), 157–192. doi:10.1016/B978-0-12-811444-5.00006-3
Theis, N., and Lerdau, M. (2003). The evolution of function in plant secondary metabolites. Int. J. plant Sci. 164 (S3), S93–S102. doi:10.1086/374190
Thirumal Kumar, D., Lavanya, P., George Priya Doss, C., Tayubi, I. A., Naveen Kumar, D. R., Francis Yesurajan, I., et al. (2017). A molecular docking and dynamics approach to screen potent inhibitors against fosfomycin resistant enzyme in clinical klebsiella pneumoniae. J. Cell. Biochem. 118 (11), 4088–4094. doi:10.1002/jcb.26064
Tiwari, R., and Rana, C. S. (2015). Plant secondary metabolites: a review. Int. J. Eng. Res. General Sci. 3 (5), 661–670.
Tong, J. B., Luo, D., Xu, H. Y., Bian, S., Zhang, X., Xiao, X. C., et al. (2021). A computational approach for designing novel SARS-CoV-2 M pro inhibitors: combined QSAR, molecular docking, and molecular dynamics simulation techniques. New J. Chem. 45 (26), 11512–11529. doi:10.1039/d1nj02127c
Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31 (2), 455–461. doi:10.1002/jcc.21334
Turner, P. J. (2005). XMGRACE, version 5.1. 19. Center for Coastal and land-margin research, OR, 2. Oregon Graduate Institute of Science and Technology, Beaverton.
Twaij, B. M., and Hasan, M. N. (2022). Bioactive secondary metabolites from plant sources: types, synthesis, and their therapeutic uses. Int. J. Plant Biol. 13 (1), 4–14. doi:10.3390/ijpb13010003
Valdés-Tresanco, M. S., Valdés-Tresanco, M. E., Valiente, P. A., and Moreno, E. (2021). gmx_MMPBSA: a new tool to perform end-state free energy calculations with GROMACS. J. Chem. theory Comput. 17 (10), 6281–6291. doi:10.1021/acs.jctc.1c00645
Walters, W. P. (2012). Going further than Lipinski's rule in drug design. Expert Opin. drug Discov. 7 (2), 99–107. doi:10.1517/17460441.2012.648612
Yang, L., Wen, K. S., Ruan, X., Zhao, Y. X., Wei, F., and Wang, Q. (2018). Response of plant secondary metabolites to environmental factors. Molecules 23 (4), 762. doi:10.3390/molecules23040762
Keywords: Internalin A, Lipinski’s rule, listeriosis, secondary metabolites, terpenes
Citation: K D and Venugopal S (2024) Molecular docking and molecular dynamic simulation studies to identify potential terpenes against Internalin A protein of Listeria monocytogenes. Front. Bioinform. 4:1463750. doi: 10.3389/fbinf.2024.1463750
Received: 12 July 2024; Accepted: 26 August 2024;
Published: 06 September 2024.
Edited by:
Vikram Dalal, Washington University in St. Louis, United StatesReviewed by:
Preeti Verma, University of Virginia, United StatesGunjan Saini, Purdue University, United States
Sanjaykumar Patel, University of Bridgeport, United States
Onkar Nerkar, Wieland (Germany), Germany
Copyright © 2024 K and Venugopal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Subhashree Venugopal, vsubhashree@vit.ac.in