- Research
- Open access
- Published:
Identification of new potent NLRP3 inhibitors by multi-level in-silico approaches
BMC Chemistry volume 18, Article number: 76 (2024)
Abstract
Nod-like receptor protein 3 (NLRP-3), is an intracellular sensor that is involved in inflammasome activation, and the aberrant expression of NLRP3 is responsible for diabetes mellitus, its complications, and many other inflammatory diseases. NLRP3 is considered a promising drug target for novel drug design. Here, a pharmacophore model was generated from the most potent inhibitor, and its validation was performed by the Gunner-Henry scoring method. The validated pharmacophore was used to screen selected compounds databases. As a result, 646 compounds were mapped on the pharmacophore model. After applying Lipinski's rule of five, 391 hits were obtained. All the hits were docked into the binding pocket of target protein. Based on docking scores and interactions with binding site residues, six compounds were selected potential hits. To check the stability of these compounds, 100Â ns molecular dynamic (MD) simulations were performed. The RMSD, RMSF, DCCM and hydrogen bond analysis showed that all the six compounds formed stable complex with NLRP3. The binding free energy with the MM-PBSA approach suggested that electrostatic force, and van der Waals interactions, played a significant role in the binding pattern of these compounds. Thus, the outcomes of the current study could provide insights into the identification of new potential NLRP3 inflammasome inhibitors against diabetes and its related disorders.
Introduction
The immune system defends against external and internal pathogens by using innate and adaptive immunity [1]. The innate immune system is the body’s first line of defense against foreign invaders. Various pattern recognition receptors (PRRs) are used by the innate immune system to determine consistent microbial patterns [2]. NOD-like receptors (NLRs) are described as intracellular PRRs which play an essential role in sensing the molecules linked with intracellular receptors and stress conditions. Thus, they sense different stimuli which present microbial infection and damage [3]. The NLRs consist of various members, including NLRP2, NLRP3, NLRP4, NLRP6, NLRP7, NLRP10, and NLRP12 [4]. Among them, NLRP3 is the most characterized member of the NLRs family. It is organized into three domains: an amino-terminal pyrin domain (PYD), a central nucleotide-binding domain (NACHT), which exhibits ATPase activity, also promotes self-oligomerization and a carboxy-terminal leucine-rich repeat (LRR) domain that acts as auto-inhibitory domain and plays an active role in signaling [5]. The PYD, NACHT, and LRR domains containing protein-3 were discovered as an important part of inflammasome [6, 7]. The inflammasomes are multi-protein complexes, composed of sensor protein (NLRP3), adaptor protein apoptosis-associated speck-like protein (ASC), and pro caspase-1 [8]. Under favorable conditions, the NLRP3 inflammasome localized at cytosol. Meanwhile, the presence of exogenous invaders or activators such as pathogens-associated molecular patterns, damage-associated molecular patterns, and environmental stress allows NLRP3 to associate with ASC by homotypic PYD interactions. Then, ASC binds with pro caspase-1 via homotypic caspase recruitment domain (CARD), and forming NLRP3 inflammasome complex. The caspase-1 facilitates the maturation of cytokines, (IL)-1β and interleukin-18, which starts the inflammation process, termed as pyroptosis [9,10,11]. Aberrant expression of NLRP3 contributes to multiple inflammatory disorders, including obesity, diabetes, dyslipidemia, hypertension, traumatic brain injury, and cerebrovascular disease [12,13,14,15,16,17]. Accumulating evidence shows the relationship between diabetes and NLRP3 inflammasome. Tschopp et al. first reported that NLRP3 act as sensor for metabolic danger that might facilitate the diabetes progression [18]. NLRP3 inflammasome activation affects insulin sensitivity and glucose tolerance, which promotes (IL)-1β and interleukin-18. The over expression of (IL)-1β causes endoplasmic reticulum stress and oxidative stress, which can lead to pancreatic cell death, affects T-cell activation, disrupt the function of islets β cells, and lead to diabetes [19]. A recent study reported that NLRP3 over expression contributes to diabetes related disorders, including atherosclerosis, diabetic nephropathy, and cardiomyopathy. Despite the advancement in therapeutic options, the management of diabetes and its complicated disorders become a significant challenge. The development of novel inhibitors by targeting the NLRP3 inflammasome, is a key player in pathogenesis of these disorders. Recently, computer aided drug discovery changed the way in which most effective drugs are designed by allowing the researchers to screen, and analyze the large database of compounds for therapeutic candidates against targeted proteins [20]. In this study, we used different computational approaches, including Pharmacophore-based virtual screening, molecular docking, MD simulations, and MM-PBSA calculations to predict the new and potential inhibitors against NLRP3 protein.
Computational methodology
Ligand-base pharmacophore generation
A pharmacophore is a chemical framework which contains the structure features of known active compounds. A ligand based pharmacophore model is developed based on known active, and specific compound, Tranilast [21] by using the pharmacophore query editor wizard of molecular operating environment (MOE). The MOE used an in-built set of pharmacophore features containing an aromatic center, H-bond donor, H-bond acceptor, cationic or anionic atom, hydrophobic atom, and Pi ring center [22]. In current study, the essential chemical features i.e. hydrogen bond donors, hydrogen bond acceptors, atom Q, and hydrophobicity were utilized to develop pharmacophore model.
Pharmacophore validation
Analysis of validated pharmacophore is essential to ensure the accuracy of molecular model and verify its reliability [23]. The developed pharmacophore model was validated by the Guner-Henry (GH) method [24], in which the internal database comprised of in-active (1319) (decoys) and experimentally active (150) compounds. The active inhibitors for NLRP3 were retrieved from literature and decoys were obtained in smile format from active inhibitors by LiDEB server, and builder in MOE was used for the generation of three dimensional structures [25], and the internal database screening was carried out through pharmacophore search protocol in MOE.
GH score was calculated by a mathematical formula.
The statistical parameter of Ht presents a number of compounds considered as hits; Ha indicates a number of actual actives in hits; A express number of active molecules in internal database; D indicates total number of molecules in internal database. The higher score of GH represents the performance of pharmacophore model is good. The Pharmacophore model was further used for virtual screening process.
Evaluation of drug like features
The pharmacophore-based virtual screening of databases are considered to be the most crucial tool for drug discovery and brings information about electronic and geometric features that are involved in binding interaction with receptor [26]. Prior to conducting the virtual screening, ZINC (12,000 molecules), In-house (1700 molecules), and Phytochemical (5000 molecules) databases were used. These ZINC, and Phytochemical databases are the largest and openly available databases, usually employed to identify the most effective inhibitors against different diseases. To generate more refined and precise databases, all the databases were cleaned, subjected to 3D protonation by using the MMFF94 force field and energy minimized with an RMS gradient of 0.05 by using the energy minimize application in MOE [27]. In addition, the force field partial charges were computed and hydrogen was added. Because of screening 349, 186, and 111 hits were extracted from ZINC, In-house, and Phytochemical libraries, respectively. Identification of drug-likeness predicts new hits by virtual screening [28]. Therefore, the screened compounds were filtered by Lipinski rule of five (RO5), and ADMET [29, 30], for searching the drug-likeness properties and pharmacokinetics, respectively. The screened compounds were short-listed by RO5 in MOE. According to RO5, a drug like molecule has a hydrogen bond donors (< 5), hydrogen bond acceptors (< 10), molecular weight (< 500 Da), polar surface area not greater than 140 Å and lipophilicity (logP < 5). Pan assay interference compounds (PAINS) (http://biosig.unimelb.edu.au/pkcsm/predication) is an electronic filter used to study the compounds quality in the database. PAINS analyze the compounds which are chemically reactive, more assay interfering, predict pharmacokinetics properties of screened compounds. Thus, it is important to study the combination of filtered compounds in order to obtain the desired pharmacokinetics properties. For the analysis, the smile format of the compounds were passed through filter in PAINS server and ADMET properties were predicted. All the compounds that obeyed the RO5 and passed the PAIN filters were subjected to molecular docking and MS simulation studies.
Molecular docking
Molecular docking analysis predicts the binding interaction between macromolecules or small molecules with protein receptors at the atomic level. It can identify the binding mode of small molecules into the protein-binding site [31,32,33,34].
Structure preparation
The three-dimensional structure of NLRP3 with (PDB ID: 7ALV) was retrieved in PDB format from RCSB Protein Data Bank. Before molecular docking, protein structure was prepared, correcting the protein structural issues, structure protonation, and structure minimized to a specific gradients, by using the quick prep option of MOE. Water molecules and other co-factors from protein structure were eliminated and H-atoms were added [35].
Active site residues prediction
For molecular docking study, active sites was selected from the literature. The active site residues of NLRP3 such as Ala226, Ala227, Arg351, Pro352, Arg578, and lys232, reported by ishania and dekker were used for molecular docking study [35, 36]. In the docking parameter of MOE, we set force field to MMFF94x, obtaining the gradient setting of 0.05 kcal/mol and the triangle matcher placement algorithm was applied. London dG method was used to score the poses [37]. After examining the docking results, total 6 hits from these databases were selected for MD simulations based on binding interaction, and S score. The S score observed the interactions, the lower S score with inhibitors interacting strongly with NLRP3 protein [38].
Molecular dynamic (MD) simulation
Combining the docking results with MD simulations allows the validation of docking results by confirming the conformational flexibility and structure stability of protein–ligand complexes. Therefore, the protein–ligand complexes were used to interact in a simulated environment for a specific period and trajectories were computed, affording the overall data by molecular motions as a function of time. MD simulation was employed via Amber2022 [39]. Force field plays significant role in MD simulation as they compute the potential energy of protein–ligand complexes [40]. In this work, generic AMBER force field (GAFF) was used for ligands and the FF14SB AMBER force field was applied for protein. Both topology and coordinates files for each system were built by tleap module of AMBER 2022 software, however, the atomic charges and topology file for ligands were built by antechamber suite in AMBER 2022 [41]. Solvation is essential as it enables studying the internal motion of protein at different temperature. A truncated octahedral box of TIP3P molecules of water was used for solvating the system and water molecules were added with LEap module of Amber2022. An appropriate amount of Na+ ions were added to neutralize the system. The Particle Mesh Ewald (PME) step was then employed to compute long range electrostatic and a cut-off distance was adjusted to 10 Å for the non-bonded interactions [42]. The SHAKE algorithm was monitored to constrain hydrogen containing bonds, after that, each system was heated from 0 to 300 K, and then and then equilibration was carried out at constant pressure and 300 K temperature. Finally, 100 ns production was completed of top six complexes along with reference compound [43, 44].
Post dynamic analysis
The generated trajectories from MD simulation of NLRP3 and hits were saved. Post MD analysis was performed by implementing the CPPTRAJ module in AMBER2022 package.
MMPBSA binding free energy
The binding free energy of protein–ligand complexes was calculated to signify their thermodynamic stability and binding affinity, which are closely associated with compound’s potency [45]. In present study, the Poisson–Boltzman (MM-PBSA) technique was used to compute the binding energy of all complexes.
The binding free energy was calculated by following equations.
where ΔG bind represents binding free energy; ΔG complex denotes free energy of complex; ΔG ligand and ΔG receptor are free energy of ligand and receptor in complex system, respectively.
Dynamical cross-correlation matrix (DCCM)
DCCM were generated to better understand protein dynamics by analyzing cross correlation shift of backbone atoms. DCCM helps to provide protein dynamic simulation, representing how atomic displacement is coupled [46,47,48], denote amino acids correlation. In order to determine DCCM, the following equation was used.
Both, Δri and Δrj were displacements from mean position of ith and jth atoms respectively. The resulting values computed were from (− 1 to 1). A negative value represents negatively correlated movement and a positive value implies positive correlated motion [49].
Result
Ligand base pharmacophore model generation and validation
The pharmacophore was generated by using known active compound, Tranilast with (Pub Chem ID 5282230). Then, internal database was generated, which is composed of 150 active compounds and 1318 inactive compounds to properly assess the discriminative ability of the generated pharmacophore model by GH method. The generated database underwent screening against the pharmacophore model to verify its accuracy. The model identified 8 compounds, of which six were active and two were inactive, which indicates that model can differentiate between active and inactive compounds. Several essential parameters were calculated including active hits (Ha), total hits (Ht), % ratio of actives, % yield of actives, Enrichment factor (E), and GH score which were illustrated in Table 1. GH score between 0.7 and 0.8 presents a good model [50]. The goodness score of our validated pharmacophore was 0.76, which indicated that our resultant pharmacophore was efficient to be further used against different databases for virtual screening purpose.
The generated pharmacophore model has seven pharmacophore features, including two Hydrophobic Atom (F1 and F3), one H-bond donor (F2), one H-bond acceptor (F4), one H-bond donor and acceptor (F5), Atom Q (F6) and one Hydrophobic (F7) as illustrated in (Fig. 1).
Virtual screening
The validated pharmacophore model was used to screen ZINC, In-house, and Phytochemical databases to identify new hits. Subsequently, because of screening, 349, 186, and 111 structurally diverse hits from ZINC, In-house, and Phytochemical databases were retrieved, respectively. Compounds that obey Lipinski’s rule of five are considered to be active in the human body [51]. Therefore, our study examined various properties including H-bond donor, molecular weight, logP value, and H-bond acceptor of compounds. As a result of Lipinski's rule of five filtering, 205, 101, and 85 hits that meet the criteria for drug-like characteristics were obtained for ZINC, In-house, and Phytochemical databases, respectively.
Molecular docking analysis
In the MOE (2016) software, molecular docking was carried out to explore the binding interaction of ligands against target protein. We docked 205, 101 and 85 diverse hits from ZINC, In-house and Phytochemical databases, respectively, into active sites of NLRP3 NACHT domain. The docking results from MOE revealed that hits could be strongly accommodated into the binding pocket of NLRP3. From each database, the top 50 ranked conformations of protein–ligand complexes were saved based on docking score. The resultant binding interaction between protein and hits were visualized. We identified 30 best hits after filtering those inhibitors, which possess interactions with active site residues.
Hits optimization on the basis of binding affinity and binding energy
The binding interaction and binding energy of selected inhibitors (30) from each library along with reference compound were studied. Finally, the 06 best hits (02 from each database) were finalized based on docking score and interaction (Table 2). As demonstrated in (Fig. 2a), the binding mode of reference compound in the active site of NLRP3 showed two hydrogen bonds with active site residues of NLRP3 (Ala227, Glu369) and one pi-cation (Arg578) along with binding score − 8.456 kcal/mol (Table 2). In ZINC database, compound ZINC12359085 exhibited one pi-H bond with Ile141 and four H-bonds with Pro352, Arg578, and Glu629 (Fig. 2b), with − 8.435 kcal/mol docking score. Similarly, compound ZINC72288245 demonstrated five hydrogen bonds (Thr439, Thr439, Arg578, Arg578 and Met661) in the binding site of NLRP3 (Fig. 2c). The compound ZINC72288245 revealed − 9.897 kcal/mol docking score.
In case of In-house database, compound BA-II-51 formed two H-bonds (Ala227, Ala228) and two pi-H bonds (Thr662) with binding pocket residues along with a docking score of − 12.654 kcal/mol (Fig. 2d). In addition, compound BA-II-45 exhibited four hydrogen bonds with various residues (Ala227, Ala228, Arg578, and Glu629) and two pi-H bonds with Thr662. For this compound, the docking score of − 10.987 kcal/mol was observed (Fig. 2e).
As depicted in Fig. 2f, compound 5,280,448, formed three hydrogen bonds with Ala227, Met408, Arg578 residues and − 7.987 kcal/mol docking score. Similarly, compound 115,089, formed three H-bonds (Ala227, Arg578, and Glu629) and pi-H bond interactions (Pro352) with − 9.478 kcal/mol docking score Fig. 2g.
The detail of binding interactions of the top 6 hits along with reference compound are presented in Fig. 2.
PAINS filter assay
During drug designing, it is important to pass the compounds into various filtration for desirable pharmacokinetics properties. Therefore, all the six compounds along with reference were successfully passed through electronic filter. The results indicate that these compounds have desirable pharmacokinetics properties. Table 3 presented the PAINS filter results, compound structures and their IUPAC names.
Molecular dynamic (MD) simulation and its analysis
The MD simulation was carried out to check the stability of selected compounds along with reference compound. The trajectories analysis was performed by using the CPPTRAJ module of AMBER software and post simulation analysis such as RMSD, RMSF, DCCM, PCA, hydrogen bond analysis, and binding free energy calculation were performed.
Root mean square deviation
Root mean square deviation evaluates the differences in the backbone of protein complexes, from its initial structural state to final conformation state [52].
The larger deviation curve represents lower stability and smaller deviation curve represents higher stability. The MD simulation results revealed that reference compound showed high deviation as compared with other selected compounds (Fig. 3). The RMSD curve for Tranilast-NLRP3 complex presenting 5–5.5 Å deviation throughout MD simulation period, presents maximum instability.
In ZINC database, the ZINC12359085-NLRP3 complex showed low RMSD value of 1.5–3.0 Å during 70 ns simulation period. Moreover, it oscillates at 3.5 Å around 75–85 ns, afterward gradually decreased (Fig. 3a). The ZINC72288245-NLRP3 complex revealed low deviations during 30 ns, after that it drastically jumped to 3.2 Å from 45 to 90 ns during MD simulation period (Fig. 3b). However, the complexes were more stable than reference compound.
The dynamic stability of BA-II-51 and BA-II-45 were computed to identify the binding stability. Interestingly, BA-II-51-NLRP3 complex presented local deviation, (1.2–2.5 Å), showed a dynamically stable behavior than other selected compounds (Fig. 3c). Therefore, our results showed that the BA-II-51-NLRP3 complex was more stable than Tranilast-NLRP3 complex, as it presented little deviation and more stable behavior than reference compound. As shown in (Fig. 3d), the RMSD curve for BA-II-45-NLRP3 complex showed 1–2.5 Å around 30 ns, but later jumped to 3.2–3.5 Å, during 30-80 ns simulation period. On the other hand, in phytochemical database, the RMSD graph for 5,280,448-NLRP3 complex exhibited 3.5-5 Å during MD simulation, showed unstable behavior like reference compound (Fig. 3e). In addition, the 115,089-NLRP3 complex exhibited high fluctuations 10–30 ns at 5 Å RMSD, which presented comparatively unstable behavior unlike other compounds. (Fig. 3f). We also calculated an average RMSD curve for complex systems. The average RMDS values for ZINC12359085-NLRP3, ZINC72288245-NLRP3, BA-II-51-NLRP3, BA-II-45-NLRP3, 5,280,448-NLRP3, 115,089-NLRP3, and Tranilast-NLRP3 were found to be 3.18 ± 0.006, 3.34 ± 0.062, 3.050 ± 004, 3.19 ± 0.005, 4.099 ± 0.008, 3.84 ± 0.006, and 4.35 ± 0.008 respectively. As minor fluctuations with less RMSD value indicate good system stability [53]. Overall the RMSD analysis indicates that the hits predicted as active against NLRP3 drug target from the ZINC, In-house and Phytochemical database, revealed low average RMSD values than reference compound (Table 4). The dynamic stability assessment of these compounds showed stable pharmacological behavior and may have better pharmacological efficiency in experimental set up.
Residues flexibility indexing
Root mean square fluctuations (RMSF) measure the residues flexibility for each complex. The region with low RMSF lead to rigidity, while the high RMSF value indicate more flexibility. In case of ZINC database, ZINC12359085-NLRP3 complex showed 1–2 Å low fluctuations for residues 50–78 and 0.1–0.4 Å for residues 120–130, as compared with reference compound. However, some regions show local fluctuations (Fig. 4a). The ZINC72288245-NLRP3 complex revealed higher fluctuations at the region of 20–40 and 50–70 while the region between 100 and 300 presented local fluctuations by comparing with Tranilast-NLRP3 complex (Fig. 4b). Interestingly, for BA-II-51-NLRP3 complex low fluctuations in 50–70 residues (2–4 Å) and 250–300 (1–1.5 Å) were observed. The BA-II-51-NLRP3 complex acts as a better candidate as compared with Tranilast-NLRP3 complex (Fig. 4c). In case of BA-II-45-NLRP3 complex same fluctuation was observed as ZINC12359085-NLRP3 and BA-II-51-NLRP3 complexes (Fig. 4d). In case of phytochemical database, 5,280,448-NLRP3 and 115,089-NLRP3 complexes, similar RMSF and local fluctuations were observed (Fig. 4e and f).
In addition, the average RMSF value of all complexes, including ZINC1239085-NLRP3, ZINC72288245-NLRP3, BA-II-51-NLRP3, BA-II-45-NLRP3, 5,280,448-NLRP3, 115,089-NLRP3, and Tranilast-NLRP3 were found to be 1.54 ± 0.05, 1.41 ± 0.04, 1.19 ± 0.03, 2.05 ± 0.06, 1.63 ± 0.04, 1.29 ± 0.02, and 1.43 ± 0.06, respectively.
The average RMSF values for all the compounds ZINC72288245, BA-II-51, and 115,089 indicating potential interaction with NLRP3 protein (Table 4).
Dynamic cross correlation matrix (DCCM)
DCCM was performed to understand the anti-correlated and correlated motion in NLRP3 and its docked complexes with ZINC12359085, ZINC72288245, BA-II-51, B-II-45, 5,280,448, 115,089. Seven plots for DCCM were generated, in which patterns of negative and positive correlated motions in protein–ligand complexes were presented in Fig. 5a–g. The deep green color indicates a positive correlation, whereas, brown color presents negative correlation. The positive correlated residues move in same direction and the negative correlated residues move in opposite direction. The DCCM analysis revealed that the binding site residues, where compound, BA-II-51 strongly bound, showed positive correlations as compared to reference compound (Fig. 5c). However, the BA-II-45, ZINC12359085, and ZINC72288245 in complex with NLRP3 showed significant positive correlation motion, with minor fluctuations were observed. In case of Phytochemical database, the compounds 5,280,448 and 115,089 in complex with NLRP3 protein, showed slightly weak negative correlation motion at 200–300 residues (Fig. 5e and f). The overall DCCM results showed that the selected compound might play significant role in the stability of these complexes.
Clustering of protein’s motion
Motion mode analysis was performed which explored the favorable conformational changes in the chemistry of bonded hits and target protein. The PCA analysis of top hits is given in Fig. 6. The ZINC12359805-NLRP3 and ZINC72288245-NLRP3 complexes consist of similar patterns of phase motion along with magnitude. The ZINC12359085-NLRP3 complex starts with red dot and end with blue dot, covering an area from − 100 to + 150 along PC1 and − 100 to + 100 along PC2 (Fig. 6a). The ZINC72288245-NLRP3 covering area of from − 175 to + 100 along PC1 and − 80 to + 80 along PC2 (Fig. 6b). The BA-II-51-NLRP3 complex was arranged and compact as compared with Tranilast-NLRP3 complex, beginning from red dot, ending to blue dot covering an area from − 120 to + 75 along PC1 and − 80 to + 80 along PC2. The Tranilast-NLRP3 complex was found to be so assembled and more dispersed, presenting that BA-II-51-NLRP3 complex showed good magnitude along with protein motion during 100 ns MD simulation period. The BA-II-51-NLRP3 complex residues are almost in same energy phase (Fig. 6c), showing stable behavior. As shown in (Fig. 6d), BA-II-45-NLRP3 complex, showed little disassembly but in arranged than Tranilast-NLRP3 complex, starting from red and ending in blue, covering space of − 75 to + 120 along PC1 and − 75 to + 90 along PC2. The PCA graph for 5,280,448-NLRP3 complex covering the space of − 150 to + 150 along PC1 and − 100 to + 130 along PC2 and 115,089 in complex with NLRP3 covering space of − 75 to + 150 along PC1 and − 75 to + 100 for PC2 (Fig. 6e and f). However, it was noted that Tranilast-NLRP3 complex was clustered and mixed (Fig. 6g).
Our results described that each compound showed different motions from other complexes. The change in motion of complexes indicated that binding of compounds induced less conformational impact on NLRP3 conformational dynamics, however, quite stabilization during MD simulation period.
Hydrogen bond analysis
Hydrogen bond analysis plays a significant role in observing the lifetime interactions between protein-ligands and helps to precise the atomic level analysis [54]. It’s possible to observe the hydrogen bond formation between NLRP3 and hits/reference at long production steps in MD in Fig. 7. The hydrogen bond analysis presented that both ZINC12359085 and ZINC72288245 compounds in complex with NLRP3 favored hydrogen bond formation as compared to Tranilast-NLRP3 complex (Fig. 7a and b). Binding of BA-II-51 with NLRP3 protein significantly increased the hydrogen bonding networks as comparison with Tranilast-NLRP3 complex (Fig. 7c). Interestingly, BA-II-45 has lesser effect in establishing hydrogen bonds (Fig. 7d). In case of 5,280,448 and 115,089 in complex with NLRP3 have less effect on hydrogen bond formation as compared with other compounds (Fig. 7e and f). Our hydrogen bond analysis indicated that most of our final compounds strongly bond with NLRP3, which may serve as a potent inhibitors.
MMPBSA analysis
The binding free energy of top hits with NLRP3 protein were computed by using the MMPBSA approach. The binding free energy and its components, van der Waals energy, surface area energy, electrostatic contribution to the solvation free energy, and electrostatic energy are enlisted in Table 5. The binding free energy and its components calculated by MMPBSA approach are enlisted in Table 5. Our result showed that ZINC12359085-NLRP3 complex has the lowest binding free energy − 27.6039 kJ/mol, followed by ZINC72288245-NLRP3 complex, with value of − 25.4018 kJ/mol. The BA-II-51-NLRP3 complex obtained binding free energy of − 29.9161 kJ/mol, lower than the Tranilast-NLRP3 complex (− 25.0211 kJ/mol). Also, the BA-II-45-NLRP3 complex has a relative binding free energy of − 29.5539 kJ/mol. In Phytochemical database, the 5,280,448 and 115,089 compounds with NLRP3 obtained an affinity of − 22.7076 and − 25.4018 kJ/mol values, respectively. Overall results revealed that selected compounds have stronger binding energy than reference compound and might have a role in complex formation with NLRP3. It was also noted that relative binding free energy of all complexes were in agreement with RMSF, RMSD, PCA, hydrogen bond analysis, and DCCM analysis.
Discussion
NLRP3 mediate caspase-1 activation and secrete cytokines IL-1β during microbial infection and cellular damage. However, its abnormal activation is accountable for diabetes, its complications, dementia and other inflammatory diseases [55,56,57]. A number of small molecule inhibitors for NLRP3 inflammasome have been identified, and some of them have demonstrated significant therapeutic potential. Some direct inhibitors are available for NLRP3 inflammasome including MCC950, CY-09, Tranilast, OLT1177, and Oridonin. However, to date, none of them have FDA approved [58]. As no FDA-approved drugs are available for NLRP3 inflammasome, so there is a need to develop safe and effective drugs for the NLRP3 inflammasome. The current study focuses on the development of new inhibitors for NLRP3 inflammasome.
Virtual screening is an essential computational approach, which is widely used in the process of drug discovery and development. The ligand-based pharmacophore model is commonly used to discover new and potent ligands or inhibitors by comparing their molecular similarity to known inhibitors for a specific drug target. This approach does not require information about the protein structure [59]. In previous study, we used pharmacophore based virtual screening for the identification of new inhibitors against STAT3 drug target [60]. In current study, we used pharmacophore based virtual screening for the identification of new inhibitors against NLRP3. Pharmacophore model was generated from reference compound and validated by GH score (0.76). The validated pharmacophore was used for the virtual screening of three databases. After pharmacophore based virtual screening, total number of 646 active hits were selected from these databases. During drug development process, it is important to analyze the pharmacological efficiency of a drug candidate. During drug-likeness analysis, total of 391 compounds were obtained from these databases and docked with NACHT domain of NLRP3. Based on docking score and binding interactions, the top six compounds, ZINC12359085, ZINC72288245, BA-II-51, BA-II-45, 5,280,448, and 115,089 were selected for MD simulation. During the docking simulation, it was observed that the interaction of six compounds exhibited good docking score and strong binding interactions compared to reference compound. In a recently reported study, molecular docking of berberine against NLRP3 was carried out [61]. The docking score of berberine was predicted as − 7.27 kcal/mol. The docking scores of our compounds ZINC12359085 (− 8.43), ZINC72288245 (− 9.89), BA-II-51 (− 12.65), BA-II-45 (− 10.98), 5,280,448 (7.98), and 115,089 (− 9.47) were good than berberine. The selected six compounds exhibited binding interaction with same active sites, as previously identified by Dekker et al. [62], Dos Santos et al. [36]. The critical residues of the NACHT domain consist of Ala227, Ala228, Arg351, Arg578, Glu629, and Thr662. By binding to these residues, it is possible to hinder the conformational change which is necessary for the initial stage of inflammasome activation [62]. By comparing this, the docking scores and binding interaction of six compounds were found to be superior compared to reference compound.
After docking, MD simulation was carried out for six compounds. These compounds represented good stability in MD simulation. The RMSD analysis indicated that the six diverse inhibitors, having better stability than reference compound. The RMSD analysis was further supported by RMSF, hydrogen bond analysis, and DCCM analysis, which revealed that all the six inhibitors revealed more stability as compared to reference compound. In addition, the MM/PBSA calculated from the last 500 frames found more negative ∆G binding value, than reference compound. Recent work, done by Patli et al. reported a similar pattern of RMSD and RMSF [63]. One of the active compounds predicted in this study with PubChem CID 5,280,448 (calycosin) is reported to be a primary active ingredient in Astragalus mungoholicus Bunge [64]. This compound possesses a number of pharmacological effects, such as anti-inflammatory [65, 66], anti-oxidant [67], and anti-osteoporotic [68] properties. The compound 5,280,448 (calycosin) can also be used to treat Osteoarthritis. It was reported that calycosin inhibits NF-κB activation and decreases the synthesis of pro-inflammatory cytokines, such as IL-1β, IL-6, and TNF-α. The active compound PubChem CID 115,089 (Enterodiol) was categorized as phytoestrogens because of its plant origin. Studies in the fields of epidemiology and pharmacology have demonstrated the protective properties of END, and in particular, its oxidation product enterolactone (ENL), against osteoporosis, cardiovascular disorders, hyperlipemia, breast, colon, prostate, and menopausal syndrome [69]. One of the identified active compounds ZINC12359085 contains carbonyl, amino functional groups, and thiadiazole ring. Thiadiazole is a five-membered scaffold, which is a common and crucial component of many drugs and it is reported that this compound may reveal antibacterial, antitubercular, analgesic, antiepileptic, antiviral, and anticancer properties [70]. The identified compound ZINC72288245 contains ethoxy, amide, sulphonyl functional group, and thiophen ring. Sulfonyl or sulfonamide functional groups are widely employed as pharmaceutical agents, and their uses in medicinal chemistry cannot be ignored [71]. Similarly, thiophene moiety is present in commercially available drugs, including Olanzapine, Benzocyclidine, Sertaconazole, Tioconazol, Dorzolamide, Tipepidine, and Tiquizium Bromide [72]. Overall, our study found that the newly identified inhibitors would be able to inhibit NLRP3 protein, which ultimately cures the disease.
Conclusion
The current study aimed to design novel NLRP3 inhibitors by employing multi-level in-silico approaches. The molecular docking analysis resulted the top six compounds, demonstrated excellent binding interactions and docking score. The hydrogen and hydrophobic interactions between NLRP3 protein and these compounds revealed the involvement of key residues, namely, Ala227, Ala228, Arg351, Arg578, Glu629, and Thr662. They also exhibited drug likeness properties. MD simulations brought valuable insights into stability and behavior of six compounds, exhibiting notable affinity and stability than reference compound. These findings provide an excellent example of pharmacophore based virtual screening as practical approach to deign novel NLRP3 inhibitors against diabetes. The outlook of the current work is that some experimental approaches are require to confirm the activity of these compounds against NLRP3.
Data availability
All the data and its links are available in the manuscript.
References
Xiao Y, Zhang L. Mechanistic and therapeutic insights into the function of NLRP3 inflammasome in sterile arthritis. Front Immunol. 2023;14:1273174.
Huang Y, Xu W, Zhou R. NLRP3 inflammasome activation and cell death. Cell Mol Immunol. 2021;18(9):2114–27.
Almeida-da-Silva CLC, et al. The role of NOD-like receptors in innate immunity. Front Immunol. 2023;14:1122586.
Gupta L, et al. Novel NLRP12 variant presenting with familial cold autoimmunity syndrome phenotype. Ann Rheum Dis. 2019. https://doi.org/10.1136/annrheumdis-2019-216158.
He Y, Hara H, Núñez G. Mechanism and regulation of NLRP3 inflammasome activation. Trends Biochem Sci. 2016;41(12):1012–21.
Jin Y, Fu J. Novel insights into the NLRP 3 inflammasome in atherosclerosis. J Am Heart Assoc. 2019;8(12): e012219.
Martinon F, Burns K, Tschopp J. The inflammasome: a molecular platform triggering activation of inflammatory caspases and processing of proIL-β. Mol Cell. 2002;10(2):417–26.
Riaz M, et al. Predicting multi-interfacial binding mechanisms of NLRP3 and ASC pyrin domains in inflammasome activation. ACS Chem Neurosci. 2021;12(4):603–12.
Mahendra J, et al. The expression of allele changes in NLRP3 (rs35829419) and IL-1β (+3954) gene polymorphisms in periodontitis and coronary artery disease. Materials. 2021;14:5103.
Li H, et al. Colchicine ameliorates myocardial injury induced by coronary microembolization through suppressing pyroptosis via the AMPK/SIRT1/NLRP3 signaling pathway. BMC Cardiovasc Disord. 2024;24(1):23.
Jorgensen I, Miao EA. Pyroptotic cell death defends against intracellular pathogens. Immunol Rev. 2015;265(1):130–42.
Whitmer RA, et al. Obesity in middle age and future risk of dementia: a 27 year longitudinal population based study. BMJ. 2005;330(7504):1360.
Amin AH, et al. Immunomodulatory effect of papaya (carica papaya) pulp and seed extracts as a potential natural treatment for bacterial stress. J Food Biochem. 2019;43;e13050.
Eglit GM, et al. Hypertension and Alzheimer’s disease: Indirect effects through circle of Willis atherosclerosis. Brain Commun. 2020;2(2):fcaa114.
Ramos-Cejudo J, et al. Traumatic brain injury and Alzheimer’s disease: the cerebrovascular link. EBioMedicine. 2018;28:21–30.
Love S, Miners JS. Cerebrovascular disease in ageing and Alzheimer’s disease. Acta Neuropathol. 2016;131:645–58.
Lin T-M, et al. Autoimmune rheumatic diseases increase dementia risk in middle-aged patients: a nationwide cohort study. PLoS ONE. 2018;13(1): e0186475.
Schroder K, Zhou R, Tschopp J. The NLRP3 inflammasome: a sensor for metabolic danger? Science. 2010;327(5963):296–300.
Verma G, Datta M. IL-1β induces ER stress in a JNK dependent manner that determines cell death in human pancreatic epithelial MIA PaCa-2 cells. Apoptosis. 2010;15:864–76.
dos Nascimento IJS, de Aquino TM, da Silva-Júnior EF. The new era of drug discovery: the power of computer-aided drug design (CADD). Lett Drug Des Discov. 2022;19(11):951–5.
Campoccia D, et al. Orthopedic implant infections: incompetence of Staphylococcus epidermidis, Staphylococcus lugdunensis, and Enterococcus faecalis to invade osteoblasts. J Biomed Mater Res Part A. 2016;104(3):788–801.
Ejaz SA, et al. Pharmacophore-based virtual screening and in-silico explorations of biomolecules (curcumin derivatives) of curcuma longa as potential lead inhibitors of ERBB and VEGFR-2 for the treatment of colorectal cancer. Molecules. 2023;28(10):4044.
Molla MHR, et al. Integrative ligand-based pharmacophore modeling, virtual screening, and molecular docking simulation approaches identified potential lead compounds against pancreatic cancer by targeting FAK1. Pharmaceuticals. 2023;16(1):120.
Kumar SP. Receptor pharmacophore ensemble (REPHARMBLE): a probabilistic pharmacophore modeling approach using multiple protein-ligand complexes. J Mol Model. 2018;24(10):282.
Gori DNP, et al. LIDeB tools: a Latin American resource of freely available, open-source cheminformatics apps. Artif Intell Life Sci. 2022;2:100049.
Dhanjal JK, et al. Use of ligand-based pharmacophore modeling and docking approach to find novel acetylcholinesterase inhibitors for treating Alzheimer’s. Biomed Pharmacother. 2015;71:146–52.
Rauf A, et al. Enzyme inhibitory activities of extracts and carpachromene from the stem of Ficus benghalensis. BioMed Res Int. 2022. https://doi.org/10.1155/2022/7053655.
Muegge I. Selection criteria for drug-like compounds. Med Res Rev. 2003;23(3):302–21.
Lipinski CA. Lead-and drug-like compounds: the rule-of-five revolution. Drug Discov Today Technol. 2004;1(4):337–41.
Wadood A, et al. In silico drug designing for ala438 deleted ribosomal protein S1 (RpsA) on the basis of the active compound Zrl 15. ACS Omega. 2021;7(1):397–408.
Singh S, Baker QB, Singh DB. Molecular docking and molecular dynamics simulation. In: Singh DB, Pathak RK, editors. Bioinformatics. Amsterdam: Elsevier; 2022. p. 291–304.
Noreen S, et al. Synthesis, characterization, molecular docking and network pharmacology of bioactive metallic sulfonamide-isatin ligands against promising drug targets. J Mol Struct. 2023;1277:134780.
Mustafa G, et al. Synthesis, structure elucidation, SC-XRD/DFT, molecular modelling simulations and DNA binding studies of 3, 5-diphenyl-4, 5-dihydro-1 H-pyrazole chalcones. J Biomol Struct Dyn. 2023. https://doi.org/10.1080/07391102.2023.2293260.
Noreen S, Sumrra SH. Correlating the charge transfer efficiency of metallic sulfa-isatins to design efficient NLO materials with better drug designs. Biometals. 2022;35(3):519–48.
Al-Karmalawy AA, et al. Molecular docking and dynamics simulation revealed the potential inhibitory activity of ACEIs against SARS-CoV-2 targeting the h ACE2 receptor. Front Chem. 2021;9:661230.
dos Santos Nascimento IJ, de Aquino TM, da Silva-Júnior EF. Molecular docking and dynamics simulation studies of a dataset of NLRP3 inflammasome inhibitors. Recent Adv Inflamm Allergy Drug Discov. 2021;15(2):80–6.
Zhang Y, et al. Identification of novel Src inhibitors: Pharmacophore-based virtual screening, molecular docking and molecular dynamics simulations. Molecules. 2020;25(18):4094.
Ajmal A, et al. Computer-assisted drug repurposing for thymidylate kinase drug target in monkeypox virus. Front Cell Infect Microbiol. 2023;13:618.
Salomon-Ferrer R, Case DA, Walker RC. An overview of the Amber biomolecular simulation package. Wiley Interdiscip Rev: Comput Mol Sci. 2013;3(2):198–210.
Grinter SZ, Zou X. Challenges, applications, and recent advances of protein-ligand docking in structure-based drug design. Molecules. 2014;19(7):10150–76.
Mahmood A, et al. A novel biallelic variant in the Popeye domain-containing protein 1 (POPDC1) underlies limb girdle muscle dystrophy type 25. Clin Genet. 2023;103(2):219–25.
Darden T, York D, Pedersen L. Particle mesh Ewald: an N⋅log (N) method for Ewald sums in large systems. J Chem Phys. 1993;98(12):10089–92.
Ryckaert J-P, Ciccotti G, Berendsen HJ. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977;23(3):327–41.
Salomon-Ferrer R, et al. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh Ewald. J Chem Theory Comput. 2013;9(9):3878–88.
Hata H, et al. Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model. Biophys Physicobiol. 2021;18:305–16.
Li X, et al. How Parkinson’s disease-related mutations disrupt the dimerization of WD40 domain in LRRK2: a comparative molecular dynamics simulation study. Phys Chem Chem Phys. 2020;22(36):20421–33.
Ni D, et al. Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy. Chem Sci. 2021;12(1):464–76.
Lu S, et al. Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators. Acta Pharm Sin B. 2021;11(5):1355–61.
Amadei A, Linssen AB, Berendsen HJ. Essential dynamics of proteins. Prot Struct Funct Bioinf. 1993;17(4):412–25.
Zhou Y, Di B, Niu M-M. Structure-based pharmacophore design and virtual screening for novel tubulin inhibitors with potential anticancer activity. Molecules. 2019;24(17):3181.
Babu S, et al. Identification of potent and selective JAK1 lead compounds through ligand-based drug design approaches. Front Pharmacol. 2022;13:837369.
Ali A, et al. Comparative binding analysis of WGX50 and Alpha-M with APP family proteins APLP1 and APLP2 using structural-dynamics and free energy calculation approaches. Phys Chem Chem Phys. 2023;25(21):14887–97.
Pundir H, et al. Using Chou’s 5-steps rule to study pharmacophore-based virtual screening of SARS-CoV-2 Mpro inhibitors. Mol Diversity. 2021;25:1731–44.
Shahab M, et al. Machine learning-based drug design for identification of thymidylate kinase inhibitors as a potential anti-Mycobacterium tuberculosis. J Biomol Struct Dyn. 2023. https://doi.org/10.1080/07391102.2023.2216278.
Rheinheimer J, et al. Current role of the NLRP3 inflammasome on obesity and insulin resistance: a systematic review. Metabolism. 2017;74:1–9.
Yang J, Wise L, Fukuchi K-I. TLR4 cross-talk with NLRP3 inflammasome and complement signaling pathways in Alzheimer’s disease. Front Immunol. 2020;11:724.
Zhou B, et al. FNDC5 attenuates oxidative stress and NLRP3 inflammasome activation in vascular smooth muscle cells via activating the AMPK-SIRT1 signal pathway. Oxide Med Cell Longev. 2020;2020:1–15.
Zahid A, et al. Pharmacological inhibitors of the NLRP3 inflammasome. Front Immunol. 2019;10:2538.
Ghufran M, et al. New lead compounds identification against KRas mediated cancers through pharmacophore-based virtual screening and in vitro assays. J Biomol Struct Dyn. 2023;41(16):8053–67.
Rafiq H, et al. Identification of novel STAT3 inhibitors for liver fibrosis, using pharmacophore-based virtual screening, molecular docking, and biomolecular dynamics simulations. Sci Rep. 2023;13(1):20147.
Chen J, et al. Berberine ameliorates inflammation in acute lung injury via NF-κB/Nlrp3 signaling pathway. Front Nutr. 2022;9:851255.
Dekker C, et al. Crystal structure of NLRP3 NACHT domain with an inhibitor defines mechanism of inflammasome inhibition. J Mol Biol. 2021;433(24):167309.
Patil SM, et al. Computational screening of benzophenone integrated derivatives (BIDs) targeting the NACHT domain of the potential target NLRP3 inflammasome. Adv Cancer Biol-Metastasis. 2022;5:100056.
Deng M, et al. Calycosin: a review of its pharmacological effects and application prospects. Expert Rev Anti Infect Ther. 2021;19(7):911–25.
Dong L, et al. Anti-inflammatory effect of Calycosin glycoside on lipopolysaccharide-induced inflammatory responses in RAW 264.7 cells. Gene. 2018;675:94–101.
Hu T, et al. Identification of bioactives from Astragalus chinensis Lf and their antioxidant, anti-inflammatory and anti-proliferative effects. J Food Sci Technol. 2017;54:4315–23.
Elsherbiny NM, et al. Renoprotective effect of calycosin in high fat diet-fed/STZ injected rats: effect on IL-33/ST2 signaling, oxidative stress and fibrosis suppression. Chem Biol Interact. 2020;315:108897.
Fang Y, et al. Calycosin stimulates the osteogenic differentiation of rat calvarial osteoblasts by activating the IGF1R/PI3K/Akt signaling pathway. Cell Biol Int. 2019;43(3):323–32.
Wang C-Z, et al. Production of enterodiol from defatted flaxseeds through biotransformation by human intestinal bacteria. BMC Microbiol. 2010;10(1):1–9.
Atmaram UA, Roopan SM. Biological activity of oxadiazole and thiadiazole derivatives. Appl Microbiol Biotechnol. 2022;106(9–10):3489–505.
Zhao C, et al. Pharmaceutical and medicinal significance of sulfur (SVI)-containing motifs for drug discovery: a critical review. Eur J Med Chem. 2019;162:679–734.
da Cruz RMD, et al. Thiophene-based compounds with potential anti-inflammatory activity. Pharmaceuticals. 2021;14(7):692.
Acknowledgements
The authors extent their appreciation to the Deputyship of Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project with number: ISP23-81.
Funding
The study was funded by the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia (Project number: ISP23-81).
Author information
Authors and Affiliations
Contributions
Conceptualization: Abdul Wadood, Asaad Khalid, Vetriselvan Subramaniyan and Mubarak A. Alamri Methodology: Chandni Hayat, Vetriselvan Subramaniyan and Ling Shing Wong Software: Mubarak A. Alamri, Ashraf N. Abdalla and Sahib Gul Afridi Data curation: Chandni Hayat, Vetriselvan Subramaniyan, Ling Shing Wong, Mubarak A. Alamri and Ashraf N. Abdalla Writing—Original draft preparation: Chandni Hayat and Abdul Wadood Visualization: Asaad Khalid, Vetriselvan Subramaniyan, Ling Shing Wong and Sahib Gul Afridi Investigation: Vetriselvan Subramaniyan, Ling Shing Wong, Abdul Wadood, Asaad Khalid and Mubarak A. Alamri Supervision: Abdul Wadood, Asaad Khalid and Vinoth Kumarasamy, Validation: Vetriselvan Subramaniyan, Sahib Gul Afridi, Mubarak A. Alamri, Ashraf N. Abdalla and Vinoth Kumarasamy Writing—Reviewing and Editing: Abdul Wadood, Vinoth Kumarasamy, Vetriselvan Subramaniyan, Ling Shing Wong, Mubarak A. Alamri, and Asaad Khalid Funding: Vetriselvan Subramaniyan and Asaad Khalid.
Corresponding authors
Ethics declarations
Ethical approval consent to participate
This article does not contain any studies with human participants or animals performed by any of the authors. We have not performed work on human/animal, therefore, there is no need of informed consent.
Consent for publication
Not applicable.
Competing interests
The authors have no competing of interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Hayat, C., Subramaniyan, V., Alamri, M.A. et al. Identification of new potent NLRP3 inhibitors by multi-level in-silico approaches. BMC Chemistry 18, 76 (2024). https://doi.org/10.1186/s13065-024-01178-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13065-024-01178-3