Next Article in Journal
A Turn-On Lipid Droplet-Targeted Near-Infrared Fluorescent Probe with a Large Stokes Shift for Detection of Intracellular Carboxylesterases and Cell Viability Imaging
Next Article in Special Issue
Identification of Potential Inhibitors for the Treatment of Alkaptonuria Using an Integrated In Silico Computational Strategy
Previous Article in Journal
Subcritical Water Extraction of Onosma mutabilis: Process Optimization and Chemical Profile of the Extracts
Previous Article in Special Issue
Virtual Screening of Hepatitis B Virus Pre-Genomic RNA as a Novel Therapeutic Target
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 5
10.3390/molecules28052315
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Rational Computational Approaches in Drug Discovery: Potential Inhibitors for Allosteric Regulation of Mutant Isocitrate Dehydrogenase-1 Enzyme in Cancers

by
Masthan Thamim
1,
Ashish Kumar Agrahari
2,
Pawan Gupta
3 and
Krishnan Thirumoorthy
1,*
1
Department of Chemistry, School of Advanced Sciences, Vellore Institute of Technology, Vellore 632014, Tamil Nadu, India
2
Translational Health Science and Technology Institute, Faridabad 121001, Haryana, India
3
Department of Pharmaceutical Chemistry, Shri Vile Parle Kelavani Mandal’s Institute of Pharmacy, Dhule 424001, Maharashtra, India
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(5), 2315; https://doi.org/10.3390/molecules28052315
Submission received: 16 January 2023 / Revised: 6 February 2023 / Accepted: 12 February 2023 / Published: 2 March 2023
(This article belongs to the Special Issue Advances in Computational Chemistry for Drug Design, Discovery and Screening)
Browse Figures
Versions Notes

Abstract

:
Mutations in homodimeric isocitrate dehydrogenase (IDH) enzymes at specific arginine residues result in the abnormal activity to overproduce D-2 hydroxyglutarate (D-2HG), which is often projected as solid oncometabolite in cancers and other disorders. As a result, depicting the potential inhibitor for D-2HG formation in mutant IDH enzymes is a challenging task in cancer research. The mutation in the cytosolic IDH1 enzyme at R132H, especially, may be associated with higher frequency of all types of cancers. So, the present work specifically focuses on the design and screening of allosteric site binders to the cytosolic mutant IDH1 enzyme. The 62 reported drug molecules were screened along with biological activity to identify the small molecular inhibitors using computer-aided drug design strategies. The designed molecules proposed in this work show better binding affinity, biological activity, bioavailability, and potency toward the inhibition of D-2HG formation compare to the reported drugs in the in silico approach.

Graphical Abstract

1. Introduction

The wild-type cytosolic and mitochondrial homodimeric isocitrate dehydrogenase (IDH1/2) enzymes often catalyze the reversible oxidative decarboxylation of isocitrate into α-ketoglutarate (α-KG) using nicotinamide adenine dinucleotide phosphate (NADP+/NADPH) as a cofactor [1]. This biochemical reaction depicts the fundamental biotransformation reaction to maintain the post-translation modifications, DNA repair mechanism, cell signaling process, lipid synthesis, antioxidant formation, and control the redox potential [2,3,4]. Modification in the cellular metabolic process is a dynamic process for cancer progression [5]. Likewise, the alteration in enzyme nature is largely associated with many biological modifications including oncogene activation [6]. The mutations in IDH enzymes, especially, have an unexpected role in the genesis and progression of human malignancies [7]. Various clinical studies also state that the somatic point mutation in the mutant IDH (mIDH) enzymes causes a broad range of cancers [8]. Frequent experimental reports confirm that the mutations in IDH1/2 are the central grounds for gliomas [9], glioblastomas [10], medulloblastomas [11], acute myeloid leukemia [12], melanoma and sporadically in melanoma [13], intrahepatic cholangiocarcinoma [14], angioimmunoblastic T cell lymphoma [15], chondrosarcoma [16], prostate cancer [17], and sporadically in thyroid, breast, stomach, and pancreatic cancers and diseases including Ollier and Maffucci syndromes [7]. Mutations in IDH1 and 2 are the fundamental hallmarks of brain cancers and they are reported up to ≥80% in WHO grade II/III astrocytomas, oligodendrogliomas, glioblastomas, and oligoastrocytomas [18]. The mutations in IDH1/2 are somatic and heterozygous missense substituted, which are recognized to be playing a dominant role in gain of function over the wild-type enzyme. Intriguingly, the recurrent mutations in IDH1 and IDH2 occur exclusively at precise arginine residues to change a different amino acid in active sites of particular hotspots namely, R100Q, R132H, R132S, R132C, R132G, and R132L in IDH1, and R140Q and R172K in IDH2. This dictates the distinct catalytic properties which activate the oncogene pathway [19]. With respect to IDH1, the most accepted variation is R132H, which contributes ≥89% of mutation and plays a paramount role in the mutation to pledge the cancer progression extensively. All the remaining disparities are addressed at low frequencies and represented as passenger mutations [20]. The specific amino acid sequence R132 in IDH1 plays a prevailing role in the active site catalytic point of the isocitrate binding pocket for changing NADP+ into NADPH which yields α-KG [21]. The mutations of this enzyme change the amino acid sequence, resulting in the loss of the normal oxidative decarboxylation of isocitrate to α-KG process and leading to the gain of the neomorphic catalytic activity that reduces α-KG to an optically active endogenous oncometabolite 2HG in specific to D form as a product by consuming NADPH [22,23] (Scheme 1). Mutations in these enzymes harbor an exclusive gain of function to vary the normal cell route to produce specific enantiomer D-2HG [24], which is reported as an important biomarker in many cancers via the disruption of DNA and histone demethylation [25,26].
To date, huge numbers of mIDH1 inhibitors are reported worldwide, and most of the inhibitors are in the pre-clinical stage, except AG-120 and IDH-305. AG-120, also known as Ivosidenib, is a potential candidate for various IDH1 mutants and the only inhibitor which successfully entered into the phase III level of various clinical tests including Cholangiocarcinoma, Leukemia, Acute Myeloid Cancer, and solid tumors [27]. IDH-305 is also a potential candidate for the mIDH1 enzyme and has successfully entered the phase II level. Because of safety issues, its use was terminated for the clinical trials of mIDH1-associated Myeloid Leukemia, Grade II and Grade III gliomas, and neurological cancers [28,29]. Given that, it is important to continue designing small-molecular inhibitors or drugs that are active organic compounds modulating cellular pathways by targeting specific proteins, which have a low molecular weight and easy mobility in the biological environment. Small-molecular drug development is a challenging task that calls for a wide range of skills and several methodologies [30]. These small molecules can be obtained by target-based drug discovery via utilizing the structure and activity relationship, which typically entails target identification, target validation, assay development, hit identification, hit-to-lead, lead optimization, candidate selection, and later development. The exploration of the structure and activity relationship in detail yields potential, safe, and efficient inhibitors for drug discovery pathways [31]. In this direction, our underlying aim is to design a robust inhibitor to identify as a potential candidate to suppress the overproduction of D-2HG via the deactivation of the mIDH1 enzyme. Based on the literature review, 62 active reported inhibitors were identified [29,32,33,34], which specifically focus on the allosteric binding pocket of the mutant IDH1 enzyme.
Computer-aided drug discovery (CADD) tools are commonly used in drug discovery pathways which is a relatively lower-cost method to deliver potential drug candidates [35]. In CADD, the most commonly used methodology is the Quantitative Structure Activity Relationship (QSAR) which has been applied for decades in the development of relationships between physicochemical properties of chemical substances and their biological activities to obtain a well-grounded statistical model to design the potential drugs [36]. The three-dimensional QSAR (3D-QSAR) model is a statistical analysis method of the structure–activity relationship through the sequence of structural analogs to examine the relationship between the ligand and the protein [37]. Likewise, Molecular Docking is the foremost method in CADD used in a hit-to-lead optimization for drug discovery pathways to predict the binding energy conformations of the optimized ligand, which helps to realize the structural interaction between the ligand and its target protein [38]. The Molecular Dynamics (MD) simulation is the technique to describe the stability of the ligand and protein complex over time, under dynamic conditions [39]. Absorption, Distribution, Metabolism, and Excretion (ADME) are the fundamental properties to confirm the organic compound as the lead molecule in the biological condition via the pharmacokinetic properties such as drug-likeness, bioavailability, and synthetic accessibility [40]. So, the periodic combination of 3D-QSAR, Molecular Docking, MD simulation, and ADME studies is considered a standard method to investigate the binding nature between the ligand and its target receptor. In the present work, the research findings are systematically validated using 3D-QSAR, Molecular Docking, MD simulation, and ADME to figure out the potential inhibitors of the potent R132H in the mIDH enzyme at the allosteric site. To address our hypothesis, uniform analyses are performed to identify the hit molecules via the structure-activity relationship of Imidazole-Pyrimidine-Oxazolidin-based derivatives [29,32,33,34] as discussed earlier using 3D-QSAR modeling to build the effective predicted half maximal inhibitory concentration (pIC50) value and identify the required structural modification via the generated 3D contour maps. Based on the results, 3459 compounds are designed and shortlisted into 600 molecules to check their binding energy using molecular docking analysis and authenticate the 3D-QSAR results to discover the ligand–receptor interactions. Based on the docking results, the designed molecules are further subjected to MD simulation and ADME studies to validate the designed molecules that are potential candidates for the mutant IDH1 enzyme at the allosteric binding pocket.

2. Result and Discussion

2.1. 3D-QSAR

The fundamental principle underlying 3D-QSAR is that the difference in structural properties (3D-descriptors) is responsible for the variations in the biological activities (IC50) of the compounds; this leads to optimizing a series of compounds and identifying the target molecules. The result of the 3D-QSAR helps to enhance the specificity and potency of the small chemical moieties via altering the structural chemical changes of the compounds which enhances the affinity of the compounds for its target to preserve its original binding mode while forming additional and advantageous connections [36,41]. Based on these techniques and incorporated with the use of rational computational approaches, the neutralized reported compounds are split into a test set (17 compounds = 27.41%) and a training set (45 compounds = 72.59%) using the ratio of 70:30 to build a better Field-Based-3D-QSAR (FB-3D-QSAR) model. The FB-3D-QSAR is an advanced technique and is an execution of the CoMFA/CoMSIA methods with a specific set of parameters. It begins with the known activity of aligned compounds and predicts the biologically active or inactive ligands via steric, electrostatic, hydrophobic, and hydrogen bond donor and acceptor fields [42]. The resulting test set and a training set of reported molecules along with their biological activity and the parameters are shown in Table 1.
From the FB-3D-QSAR results, the obtained statistical data are well correlated and validated based on the references, as projected in Table S1 as the supporting information. The partial-least-square (PLS) factor 7 shows better results compared to all other factors where it contains minimal standard deviation (S.D = 0.1984), higher regression coefficient (R2 < 0.6 = 0.90), lower variance ratio, and root-mean-square error (P = 4.40 × 10−17 and RMSE = 0.71), and acceptable cross-validated correlation coefficient (Q2 < 0.5 = 0.5123) [43].
The reference ligands are aligned via the Tanimoto similarities [44] and clustered [45] based on the related central nucleus with limited disparities, as shown in Figure 1, along with the physicochemical properties. The most biologically active compound is selected to modify the favorable functional node based on the R-group enumeration method [46] to design the potential candidate.
Among the five physiochemical properties, the steric and hydrophobic functions are playing a major role to determine and enhance the predicted biological activity of the designed compounds. Here, the steric function contributes roughly 45% and the hydrophobic factor contributes up to 29%, which is shown in Table S2 in the Supplementary Materials. Based on contour map analysis, the favorable modification in the selected molecules can result in improved activity via the functional group modifications. The most anticipated modifications of favorable and non-favorable functional groups for steric (S and S′) and hydrophobic (H and H′) properties are shown in Figure 2.
The most expected favorable and non-favorable modifications of steric and hydrophobic functional nodes are highly available in the oxazolidine, followed by imidazole and pyrimidine moieties. The favorable steric modification (S3) in the imidazole moiety can observe the steric and hydrophobic repulsion through the neighboring methyl group. To overcome this effect, the hydrophilic modification or non-favorable hydrophobic modification (H2′), is required. Likewise, the favorable (H3, H4, and S4) and non-favorable (H3′, S3′, H4′, and S4′) modifications of steric and hydrophobic functional groups enhance the activity of the designed compounds in a minimal range. Similarly, the favorable modification of steric and hydrophobic properties (S1, H1, and H2) in the oxazolidine moiety highly augments the activity of the designed compounds due to the free rotation away from the steric groups. The steric modification (S1), especially, will boost the activity of the designed compound; specifically, below the plane (dative bond), not above the plane, and it is incorporated in the ligand design. Likewise, the favorable hydrophobic (H2) modification will increase the activity only above the plane and not below. The non-favorable steric modification (S2′) is highly correlated with the favorable steric point (S2) of the pyrimidine moiety due to the pyrimidine–oxazolidine rotation (N-C), and the non-favorable hydrophobic modification (H1′) arises due to the availability of neighboring methyl group. The non-favorable steric modification S1′ avails above the plane to destabilize the oxazolidine moiety. From the FB-3D-QSAR, one can conclude that the small modification in the selected groups with specific functional nodes can dramatically increase the activity of the designed compound. Based on the contour map results, the high biological activity of molecule RC-01 was taken as the template molecule for further studies.

2.2. Molecular Docking

One of the most essential and widely used tools to deliver the potential drug candidate using low cost computer-assisted drug designing is molecular docking. The core of molecular docking is to place the small molecules in the active region of the target enzyme and use scoring functions to estimate a compound’s biological activity to predict the ligand accessibility in the biological environment. Molecular docking provides essential information about the three-dimensional structure possibilities of the specific ligand in the various regions of its target protein connectivity to list out via the scoring functions [47]. Based on this scoring value, one can determine the binding affinity of the ligand–protein interactions. From the FB-3D-QSAR results, the template molecule can allow modification with the selected nodes of the specific functional groups to check their binding affinity using the R-group enumeration method. Based on the available space in each functional node, it is possible to generate more than 2750 entries using the diverse R-group functionalities [46] such as steric effects, electronegativity, hydrophobicity, and hydrogen bond donor and acceptor nature. The limited entries were specifically selected based on the contour map results and designed 3459 ligands, especially using steric and hydrophobic functional nodes. A total of 936 designed molecules were selected and further allowed for the ligand preparation analysis [48] and finally, 600 compounds were shortlisted based on the highest predicted biological activity, which is further examined in the molecular docking study [49]. The molecular docking results are shown in Table S3 (reference compounds) and Table S4A,B (designed compounds) are in Supplementary Materials. Based on the molecular docking results, two ligands were chosen as reference compounds from the reported drugs, where one has the highest biological activity (RC-01) and another one contains the highest binding energy (RC-02), they are considered the control group to validate the newly designed compounds in this work. From the docking results, the top five designed ligands are selected for further studies (DC-01 to DC-05). The selected designed ligands and their parameters, along with the reference molecules, are shown in Table 2.
Based on the FB-3D-QSAR modeling, the specific steric active modification, especially in the fourth position of the oxazolidine moiety, results in the enhanced predicted biological activity which improves the binding affinity of the ligand–receptor interaction. The pyrrolidine, piperidine, and butyl functional groups enhance the steric nature (highlighted in pink), especially in DC-01 (S-3-fluoro-3-methylpyrrolidine), DC-02 (R-3-methylpyrrolidin-3-ol), DC-03 (3,3-dimethylpyrrolidine), DC-04 (piperidin-4-ol), and DC-05 (S-2,3-dimethylbutan-1-aminium). From the molecular docking results, it is very clear that the reference and designed ligands are perfectly fitted in the allosteric binding site and are shown in Figure 3. The 2D interactions of receptor–ligands obtained from the molecular docking results are shown in Figure S1A–G in Supplementary Materials.
From the molecular docking results, it is deduced that the functional node modification using R-group enumeration is the potential target to augment the hydrophobic interactions via pi–pi and alkyl–alkyl interactions due to the steric and hydrophobic functional group modifications. The reference ligands contain the least hydrophobic interaction in the oxazolidine moiety except for one alkyl–alkyl interaction with the VAL-281 residue. The main aim is to improve the activity of the designed compound without disturbing the existing interactions in the reference compounds via the FB-3D-QSAR approach. Thus, the modified ligands showcase the reference ligand interactions along with the enhanced interactions, especially in the oxazolidine site; this is one of the more robust predictions of the rational computational approach. The designed ligands are well enriched with hydrophobic and steric interactions in the preferable site to increase the activity and affinity. The molecular docking results of receptor–ligand interactions with specific residues of reference compounds along with designed compounds are shown in Table S5 in Supplementary Materials. The strong hydrophobic connectivity of selective residues ILE-113, ILE-117, TYR-285, MET-290, and CYS-379 enhance the alkyl–alkyl and pi–alkyl interactions to improve the binding affinity of the receptor–ligand complex. Similarly, the ILE-112, ILE-117, and ILE-118 residues involve a strong hydrogen bond, halogen-induced hydrogen bond (Fluorine-HB), and non-classical hydrogen bond interaction. Based on the molecular docking results, one can conclude that the high binding affinity molecules show higher interaction with their target protein, which is the potential key for robust drug delivery.

2.3. MD Simulation

Molecular dynamics (MD) simulation is the foremost method used and it creates a paradigm in computer-assisted drug discovery pathways to conclude if the designed compounds are biologically active or inactive. The fundamental principle of MD simulation is to analyze the physical movement of atoms in the molecule by using intermolecular interactions under dynamic conditions over a time period to provide abundant information such as enzymatic favorable reactions, chemical pathways, thermodynamic and kinetic stability, and so on. For the flair drug discovery identification, the designed compounds must fulfill all the conditions in the MD simulation study which allow interaction with its target protein under the biological condition over a period to achieve results confirming whether the specific compound is biologically stable or unstable [50]. The selected designed compounds DC (01 to 05) and reference compounds are further subjected to an MD-simulation study [51] to identify the stability of the receptor–ligand complex in the biological environment. The results depict the ligand stability through the root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), the radius of gyration (Rg), and hydrogen bond (HB) plots, which all state that the designed compounds have better stability for 50 ns compared to the reference compounds.

2.3.1. Root-Mean-Square Deviation (RMSD)

The RMSD of carbon alpha (C-alpha) of all the investigated systems was calculated to determine the convergence, i.e., stability of the trajectory. The average RMSD values of the designed compounds DC-03 (0.18 nm), and DC-04 (0.19 nm) are slightly lower than the reference compounds RC-01 (0.25 nm) and RC-02 (0.20 nm), which indicates that the designed complex is comparably more stable than the reference complex which is shown in Figure 4. However, the average RMSD value of DC-01 (0.20 nm) displayed similar stability compared to RC-02 and DC-02 (0.24 nm) and displayed higher stability than RC-01 and lower stability than RC-02. While the apo-protein average RMSD value is (0.50 nm) and the designed compound DC-05 (0.90 nm) shows an extremely high value, which indicates that it is unstable in the protein–ligand environment. As the DC-05 complex system is unstable, it is excluded from further studies.

2.3.2. Root-Mean-Square Fluctuation (RMSF)

The inherent fluctuation of the receptor–ligand is frequently linked to its function. The residue-wise protein fluctuation factor is determined, which is directly related to the biological function of the specific ligand. In Figure 5, the overall RMSF value of the designed compound DC-01 shows a slightly high fluctuation in two specific regions (140–180 and 280–290) compared to the reference compound RC-01. Whereas, DC-02 shows identical fluctuation compared to RC-01, except it is in a different region (70–90). The overall RMSF value of DC-03 was considerably high in the zones of 70–90, 140–150, and 160–180 as compared to the RC-01. Similarly, DC-04 also shows an acceptable high fluctuation in specific regions (70–90, 110–120, and 270–280) as compared to the RC-01 where the average RMSF value of the RC-02 is highest among all the complex systems. This indicates that all of the designed compounds (DC-01 to DC-04) have comparable stability to the reference compounds RC-01 and RC-02.

2.3.3. Radius of Gyration (Rg)

The Rg is a very important parameter for measuring the structural compactness, overall folding, and shape of the protein. Figure 6 exhibits the Rg plot for all the complexes (DC and RC). The average Rg values of DC-01 (2.30 nm), DC-03 (2.26 nm), and DC-04 (2.28 nm) are lower than RC-01 (2.32 nm), which indicates that they have higher compactness, and DC-02 (2.33 nm) compactness is nearly similar to RC-01. The average Rg value of RC-02 is (2.36 nm) shows the least compactness in comparison to all other complexes, which designates that the designed compounds, DC-01 to DC-04, are much more stable and active in the biological environment compared to the reference compounds RC-01 and RC-02.

2.3.4. Hydrogen Bond (HB) Interaction

The HB interaction is very essential to determine the receptor–ligand interaction in dynamic conditions over time. Figure 7 shows the HB plot for all the investigated systems including reference and designed compounds. The average number of HB interactions for DC-01(1.833), DC-02(2.428), DC-03(2.164), and DC-04(2.901) are lower compare to RC-01(4.059) and RC-02(3.191) up to 50 ns. Comparably DC-02, DC-03, and DC-04 have acceptable HB average values but are not higher compared to the reference molecules. To overcome this consequence, the ligand–receptor interaction is analyzed at the final frame of the overall simulation (50 ns).
Figure 8 shows the final frame (50 ns) of receptor–ligand interaction after the MD simulation using the Groningen Machine for Chemical Simulation (GROMACS v2022.4) package [51]. The 2D representation of receptor–ligand interaction after the MD simulation at the final frame (50 ns) is shown in Figure S2A–F in Supplementary Materials.
From the receptor–ligand interaction after the simulation at 50 ns as shown in Figure 8A–F, it is very clear that the designed ligands have better interaction compared to the reference compounds. The total number of residues involved in ligand–receptor interaction after the simulation is shown in Table 3.
Here, the reference compound RC-01 has 12 interactions including 9 hydrophobic (alkyl–alkyl and pi–alkyl) and 3 strong HBs with the help of eight active residues around the ligand. Similarly, RC-02 has a total of 13 interactions which include 8 hydrophobic (alkyl–alkyl, pi–alkyl,) 3 strong HBs, one pi–pi stacking, and pi–sulfur with the help of eight active residues around the ligand. However, the designed compounds contain a higher number of interactions over the reference compounds. DC-01 has a total of 20 interactions which include 3 classical and non-classical HBs, one halogen-induced HB (Fluorine-HB), 2 pi–pi stacking, and 11 hydrophobic interactions (alkyl–alkyl and pi–alkyl) with the help of 12 residues around the ligand pocket site. Similarly, DC-02 has 20 interactions which include 5 strong HBs, 3 pi–pi stacking, and 12 hydrophobic (alkyl–alkyl and pi-alkyl) with the help of 12 active residues. Furthermore, DC-03 contains 24 interactions which include 3 classical and non-classical HBs, one pi–sigma and pi–sulfur, and 16 hydrophobic (alkyl–alkyl and pi–alkyl) with the help of 14 active residues around the ligand pocket. DC-04 contains 18 interactions which include 4 strong classical HBs, one non-classical HB, one pi–sigma and pi–pi, and 11 hydrophobic (alkyl–alkyl and pi–alkyl) with the help of nine active residues. Among every designed compound, DC-03 had the highest number of interactions with 24, DC-01 and DC-02 had the same number of interactions with 20, and DC-04 had the lowest number of interactions with 18. Thus, all the designed compounds have a higher number of interactions compared to the reference compounds, RC-01(12) and RC-02(13). This dictates that our designed compounds are much more stable in the dynamic condition up to 50 ns.

2.4. ADME

A fundamental segment of drug discovery and development is addressing the pharmacokinetics and metabolism of the designed drug, which is often denoted as ADME properties [40]. The therapeutic effect of an effective drug is not only dependent on its biological activity, but also has to have proper ADME properties. Thus, it is very important to address the pharmacokinetic properties of the designed compound using in silico modeling to validate it as a therapeutic candidate. The ADME study in CADD tools is an essential and effective method to reduce the number of animal models in in vivo tests. The biologically active small molecules or drug candidates must follow the specific set of parameters which are widely denoted as pharmacokinetic properties. These biologically active designed molecules are called preclinical candidates and it is must have the affinity to be absorbed into the specific target surface, completely distributed to its target binding site, appropriately metabolized in the liver, and excreted after functioning. For this activity, a potent drug compound must follow Lipinski’s rule of five (Ro5) for excellent drug-likeness and absorption activity. The physicochemical properties on which Lipinski’s based drug-likeness comprise the number of H-bond donors (must be ≤5) and acceptors (must be ≤10), molecular weight (not more than 500 Da), and partition coefficient or lipophilicity as denoted as Log P (must be ≤5). Ligands, phytochemicals, and bioactive compounds derived from natural or synthetic sources or designed using CADD tools that would violate more than one of these rules are considered to have poor absorption [52]. The online free SwissADME server [53] is used to predict the drug bioavailability via the ADME property and confirm the designed drug molecules which obey Lipinski’s rule of five [54] in the biological environment. Overall, the ADME properties of our designed drugs along with reported drugs are shown in Table 4 with the required parameters.
From Table 4, one can conclude that our designed compounds have passed 90% in the drug-likeness test, and all have acceptable ranges except the MW, which is allowed in the drug discovery pathway [55]. Based on the results, the designed compound, DC-04 has highest solubility and PSA value, followed by DC-02, compared to all the other compounds. Similarly, DC-01 and DC-03 have higher PSA values and lower solubility compared to the reference compounds.

3. Conclusions

Abnormal activity in the mutated metabolic pathway leads to oncometabolite formation, which is the hallmark of cancers. Identifying potential drug candidates in cancer treatment is a challenging task for ongoing cancer research. The mIDH enzyme, especially, has gained huge attention in recent years to treat various cancers and other disease pathways. Specifically, the cytosolic mIDH1 enzyme gains the abnormal activity to be involved in the non-favorable reaction that produces the enantiomeric selective metabolite D-2HG. This metabolite is projected as a distinguished oncometabolite primarily in gliomas, glioblastomas, medulloblastomas, acute myeloid leukemia, melanoma and sporadically in melanoma, intrahepatic cholangiocarcinoma, angioimmunoblastic T cell lymphoma, and chondrosarcoma, and accountably in prostate, thyroid, breast, stomach, and pancreatic cancers and diseases including Ollier and Maffucci syndromes. The mutations in IDH enzymes are heterozygous somatic and recurrent modifications that take place in specific arginine residues to change the catalytic activity and endorse the abnormal activity involving the oncogene pathway. The R132H is the most accepted variation in the mIDH1 enzyme and it acts as a fuel to stimulate the oncogene pathway in more than 89% of human cancers and diseases. To date, the drug discovery pathway of the mIDH1 enzyme is quite complicated and only two inhibitors have been identified as potential candidates, namely, AG-120 and IDH-305. The AG-120 is the only inhibitor that successfully cleared phase III trials and is an accepted drug candidate for Cholangiocarcinoma, Leukemia, Acute Myeloid Cancer, and solid tumors; however, IDH-305 failed in phase III and the further developmental process is ongoing in IDH-305 modification. Examination of the combination of 3D-QSAR, molecular docking, MD simulation, and ADME studies yields the potential, safe, and efficient inhibitors for drug discovery pathways. This present work with the aid of CADD strategies aimed to identify the potential inhibitor for mIDH1, especially in the allosteric site to suppress the overproduction of D-2HG. Based on this direction, the Schrodinger suite 2022-3 was used to identify 62 allosteric site-reported inhibitors to build better FB-3D-QSAR modeling which obeys the required statistical parameters. From the FB-3D-QSAR method, the positive functional nodes are recognized and modified using the R-group enumeration method to design the robust derivatives of a highly active reference compound (RC-01) to predict biological activity. The designed compounds show high predicted biological activity compared to the reference compounds and selected designed compounds were further tested to understand the receptor–ligand interaction via molecular docking analysis to calculate the binding energy using the Glide package. The designed compounds DC-01 (−13.336 kcal/mol), DC-02 (−13.175 kcal/mol), DC-03 (−13.159 kcal/mol), DC-04 (−13.131 kcal/mol), and DC-05 (−12.952 kcal/mol) show high binding affinity compared to the reference compounds RC-01 (−11.800 kcal/mol) and RC-02 (−12.403 kcal/mol). These compounds were further tested to understand the stability via the MD simulation using the GROMACS v2022.4 package. From the MD simulation results, DC-01 to DC-04 show the highest stability compared to the reference compounds. The average RMSD value of DC-03 is 0.18 nm and DC-04 is 0.19 nm, which is a very low deviation compared to the reference compounds; RC-01 is 0.25 nm and RC-02 is 0.20 nm. Similarly, DC-01 has a value of 0.20 nm and DC-02 is 0.24 nm, which are an acceptable deviation as per the references. The designed compound DC-05 shows a high deviation in the RMSD of 0.90 nm compared to the apo-protein, which exhibits 0.50 nm. The RMSF plot of the designed compounds (except DC-05) shows acceptable fluctuation compared to the reference compound RC-01, where RC-02 shows higher fluctuation compared to all the other investigated systems. Similarly, the Rg plot shows very good compactness of designed compounds compared to the reference compounds. DC-01 has a value of 2.30 nm, DC-03 is 2.26 nm, and DC-04 is 2.28 nm and these are a higher compactness of ligands in the receptor region compared to RC-01 (2.32 nm) and RC-02 (2.36 nm). DC-02 has a value of 2.33 nm, which represents having similar and higher compactness compared to RC-01 and RC-02. The hydrogen bonding interaction analysis reveals that the average number of HBs is lower in all the designed compounds compared to the reference compounds. However, newly designed compounds are much augmented with hydrophobic and steric functional nodes. Thus, the types of interactions available in every compound after the simulation are explored, and surprisingly, our designed compounds have a higher number of interactions compared to the reference compounds, which drives the high stability over the reference compounds. The highest number of interactions observed for the investigated compounds after the simulation up to 50 ns is 24 in DC-03, followed by 20 in DC-01 and DC-02, 18 in DC-04, 12 in RC-01, and 13 in RC-02. The pharmacokinetic properties (ADME and Lipinski’s rule) are well correlated (except MW) to our designed compounds along with reference compounds. They exhibit pharmacokinetic, bioavailability, and drug-likeness properties. The results obtained from the rational computational approaches for the designed compounds show better stability, activity, bioavailability, and binding affinity compared to the reference compounds. This delineates that our designed compounds with in silico approaches using CADD tools could be potential candidates for an mIDH1 enzyme drug at the allosteric site to control the overproduction of D-2HG. The present work delineates the new strategies that would minimize the experimental cost involved in drug design and screening in laboratory research. This work provides new insights for researchers to synthesize the proposed class of molecules in the laboratory. Further, biochemical experimental investigations would be required to validate the designed compounds that could be therapeutic drug candidates for mIDH1 enzymes in various cancers.

4. Methodology

4.1. Data Set

The 62 specific allosteric binding inhibitors with similar molecular nuclei were taken from the existing literature and used in the present work. The 62 molecules showed a similar assay technique with a substantial disparity in their geometry and potency. The values of half maximal inhibitory concentration (IC50) varied from 49µM to 1nM, which were converted into pIC50 values using the following equation and converted into relatable values from 5.000 to 9.000.
pIC 50 = − log 10 IC 50
All the computational calculations were conducted using the Schrodinger suite 2022-3, which we acquired from www.schrodinger.com and used from 6 September 2022 to 27 September 2022 as a trial package, and the 3D geometries of the existing compounds were generated via the builder panel Maestro 13.3 [56] and consequently neutralized using Ligand preparation with Wizard segment [48]. We limited the reported ligand size to 500 atoms with default options such as desalt, determined the chirality from the 3D structure and disabled the tautomer generation, and executed the energy neutralization using the OPLS-2005 force field [57] for better results.

4.2. 3D-QSAR

The optimized compounds were treated in fingerprint analysis using Canvas similarities [45] and clustered to identify the mismatched molecular nucleus using Tanimoto analysis [44] to build a better 3D-QSAR model. We used the Field-Based (FB) 3D-QSAR techniques [42] to incorporate the 62 reported compounds into 70:30 ratios in the test set and training set for better results. We built the FB-3D-QSAR model using the Gaussian filed style and unchecked the use of input partial charges to generate electrostatic interactions. The pIC50 values of mIDH1 were used as the dependent variables and generated 3D molecular field descriptors were used as independent variables. Partial-least-square (PLS) procedure was used to correlate linearly in the FB-3D-QSAR method to limit the maximum PLS factor up to 7 to avoid over-fitting and Leave-one-out (LOO) cross-validation method is used to determine the cross-validation coefficient (R2), correlation coefficient (Q2), root-mean-square error (RMSE), the standard deviation of the regression (SD), and Pearson coefficient (Pearson-R). From these results, we explored the 5 physicochemical properties namely Gaussian Steric (GS), Gaussian Electrostatic (GE), Gaussian Hydrophobic (GH), Gaussian Hydrogen Bond Acceptor (GHA), and Gaussian Hydrogen Bond Donor (GHD) fields, which were visualized in the QSAR visualization module.

4.3. Molecular Docking

The Protein Preparation Wizard [58] in Maestro 13.3 was used to prepare the crystal structure of mIDH1 enzyme complex with IDH-305 ligand (PDB: 6B0Z). This was obtained from the RCSB protein data bank to remove water and heteroatoms, add hydrogen and partial charges, assign protonated state at pH 7.0 ± 2.0 using Epik mode [59], fill the missing loops and side chains using Prime [60], optimize with appropriate pH, and finally, minimize the energy of the complex using OPLS4 force field [61]. Then we created the grid box size 10Å × 10Å × 10Å (x, y, z coordinates 27.16, 25.74, 29.27) to delineate the binding site of mIDH1 using the ligand position using the Receptor Grid Generation module in Maestro13.3 with default settings. All the ligands were optimized with the aforementioned neutralization methodology to generate an ionization state at pH 7.0 ± 2.0 in Epik mode [59] with the OPLS-2005 force field [57]. Finally, we performed molecular docking using Glide version 2022-3 [49] standard precision method with default settings. All the ligand–receptor interactions obtained from molecular docking were visualized with the help of the Discovery Studio—2021 package [62].

4.4. MD Simulation

To examine the essential residues interaction in ligand binding, we performed the MD simulation for mIDH1 protein and ligands based on the molecular docking results using the GROMACS v2022.4 package [51]. The highest biological activity and binding energy of reported compounds were taken as the reference compounds to test our designed compounds. The high binding affinity and predicted biological activity of selected protein–ligand complexes from the results were exported. They were then used to create the topology parameters and coordinate files of receptor and hit molecules using the CHARMM36 force field [63] in GROMACS and SwissParam [64], respectively. The docked protein–ligand complex was solvated using a transferable intermolecular potential 3 point (TIP3P) water solvation model [65] enclosed within the dodecahedron box of 1 Å and neutralized via sufficient numbers of Na+ and Cl− ions added via the Monte Carlo ion-placing method [66]. Further, the energy minimization was performed by using the steepest descent approach by applying a maximum force of 10 kJ/mol to avoid steric hindrance up to 500 ps. Berendsen coupling [67] and the Parrinello–Rahman method [68,69] were employed to regulate temperature (NVT) and pressure (NPT) at 300 k and 1 bar, respectively, inside the box up to 100 ps. The overall simulation time was fixed as 50 ns for the reference and the designed complexes to determine the ligand–protein stability.

4.5. ADME

To address the pharmacokinetics and drug-likeness properties, we utilized the online free SwissADME [53] server to check the ADME properties along with the Lipinski rule of five to confirm the drug ability of our designed and reported compounds.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules28052315/s1, Figure S1A: The molecular docking results for reference compound (RC-01) interaction with the receptor; Figure S1B: The molecular docking results for reference compound (RC-02) interaction with the receptor; Figure S1C: The molecular docking results for designed compound (DC-01) interaction with the receptor; Figure S1D: The molecular docking results for designed compound (DC-02) interaction with the receptor; Figure S1E: The molecular docking results for designed compound (DC-03) interaction with the receptor; Figure S1F: The molecular docking results for designed compound (DC-04) interaction with the receptor; Figure S1G: The molecular docking results for designed compound (DC-05) interaction with the receptor; Figure S2A: The MD simulation results for reference compound (RC-01) interaction with the receptor at the final frame (50 ns); Figure S2B: The MD simulation results for reference compound (RC-02) interaction with the receptor at the final frame (50 ns); Figure S2C: The MD simulation results for designed compound (DC-01) interaction with the receptor at the final frame (50 ns); Figure S2D: The MD simulation results for designed compound (DC-02) interaction with the receptor at the final frame (50 ns); Figure S2E: The MD simulation results for designed compound (DC-03) interaction with the receptor at the final frame (50 ns); Figure S2F: The MD simulation results for designed compound (DC-04) interaction with the receptor at the final frame (50 ns); Table S1: The PLS statistical data of the homogenous 3D-QSAR model using the Field-Based method; Table S2: The physiochemical property contribution of reported compounds using the FB-3D-QSAR method; Table S3: The molecular docking results for 62 reported compounds using Glide package in Schrodinger suite; Table S4: The molecular docking results for 229 designed compounds using Glide package in Schrodinger suite; Table S4A: The molecular docking results for 371 designed compounds using Glide package in Schrodinger suite; Table S5: The specific residues involved in the receptor-ligand interaction of reference and designed compounds using molecular docking analysis.

Author Contributions

Conceptualization, M.T.; Data curation, M.T., A.K.A. and P.G.; Formal analysis, M.T., A.K.A. and P.G.; Funding acquisition, K.T.; Investigation, M.T.; Methodology, M.T.; Resources, K.T.; Software, M.T.; Supervision, K.T.; Validation, M.T.; Visualization, M.T.; Writing—original draft, M.T.; Writing—review and editing, K.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Department of Science and Technology (DST) India under National Supercomputing Mission (NSM) DST/NSM/R&D-HPC-Applications/2021/03.27.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary material.

Acknowledgments

We acknowledge Vinod Devaraji, Senior Scientist at Schrödinger, Bengaluru, Karnataka, India for his extended support. The support and the resources provided by PARAM Sanganak under the National Supercomputing Mission, Government of India at the Indian Institute of Technology, Kanpur are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Cairns, R.A.; Mak, T.W. Oncogenic isocitrate dehydrogenase mutations: Mechanisms, models, and clinical opportunities. Cancer Discov. 2013, 3, 730–741. [Google Scholar] [CrossRef] [Green Version]
  2. Judge, A.; Dodd, M.S. Metabolism. Essays Biochem. 2020, 64, 607–647. [Google Scholar] [CrossRef]
  3. Badur, M.G.; Muthusamy, T.; Parker, S.J.; Ma, S.; McBrayer, S.K.; Cordes, T.; Magana, J.H.; Guan, K.L.; Metallo, C.M. Oncogenic R132 IDH1 Mutations Limit NADPH for De Novo Lipogenesis through (D)2-Hydroxyglutarate Production in Fibrosarcoma Sells. Cell Rep. 2018, 25, 1018–1026.e4. [Google Scholar] [CrossRef] [Green Version]
  4. Xiong, Y.; Xu, W.; Yang, H.; Wang, P.; Ma, S.; Zhao, S.; Ye, D.; Guan, K. The mechanisms of oncometabolites in epigenetic control, DNA repair, neural development and gliomagenesis. Neurol. Oncol. 2017, 19 (Suppl. S3), iii68. [Google Scholar] [CrossRef] [Green Version]
  5. Thamim, M.; Thirumoorthy, K. Computational studies of selective N-methylation in nicotinamide: Epigenetic reprogramming in cancer. Comput. Theor. Chem. 2021, 1194, 113058. [Google Scholar] [CrossRef]
  6. Oermann, E.K.; Wu, J.; Guan, K.L.; Xiong, Y. Alterations of metabolic genes and metabolites in cancer. Semin. Cell Dev. Biol. 2012, 23, 370–380. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, F.; Travins, J.; De La Barre, B.; Penard-Lacronique, V.; Schalm, S.; Hansen, E.; Straley, K.; Kernytsky, A.; Liu, W.; Gliser, C.; et al. IDH mutations and tumorigenicity. Mol. Cancer Ther. 2013, 12 (Suppl. S11), PL02-04. [Google Scholar] [CrossRef]
  8. Prensner, J.R.; Chinnaiyan, A.M. Metabolism unhinged: IDH mutations in cancer. Nat. Med. 2011, 17, 291–293. [Google Scholar] [CrossRef]
  9. Al-Khallaf, H. Isocitrate dehydrogenases in physiology and cancer: Biochemical and molecular insight. Cell Biosci. 2017, 7, 37. [Google Scholar] [CrossRef]
  10. Mirchia, K.; Richardson, T.E. Beyond IDH-mutation: Emerging molecular diagnostic and prognostic features in adult diffuse gliomas. Cancers 2020, 12, 1817. [Google Scholar] [CrossRef]
  11. Faiyaz-Ul-Haque, M.; Al-Sayed, M.D.; Faqeih, E.; Jamil, M.; Saeed, A.; Amoudi, M.S.; Kaya, N.; Abalkhail, H.; Al-Abdullatif, A.; Rashed, M.; et al. Clinical, neuroimaging, and genetic features of L-2-hydroxyglutaric aciduria in Arab kindreds. Ann. Saudi Med. 2014, 34, 107–114. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, J.Y.; Lai, Y.S.; Tsai, H.J.; Kuo, C.C.; Yen, B.L.; Yeh, S.P.; Sun, H.S.; Hung, W.C. The oncometabolite R-2-hydroxyglutarate activates NF-κB-dependent tumor-promoting stromal niche for acute myeloid leukemia cells. Sci. Rep. 2016, 6, 32428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mullarky, E.; Mattaini, K.R.; Vander Heiden, M.G.; Cantley, L.C.; Locasale, J.W. PHGDH amplification and altered glucose metabolism in human melanoma. Pigment Cell Melanoma Res. 2011, 24, 1112–1115. [Google Scholar] [CrossRef] [PubMed]
  14. Molenaar, R.J.; Radivoyevitch, T.; Maciejewski, J.P.; van Noorden, C.J.F.; Bleeker, F.E. The driver and passenger effects of isocitrate dehydrogenase 1 and 2 mutations in oncogenesis and survival prolongation. Biochim. Biophys. Acta Rev. Cancer 2014, 1846, 326–341. [Google Scholar] [CrossRef]
  15. Qiu, Z.J.; Lin, A.P.; Jiang, S.; Elkashef, S.M.; Myers, J.; Srikantan, S.; Sasi, B.; Cao, J.Z.; Godley, L.A.; Rakheja, D.; et al. MYC Regulation of D2HGDH and L2HGDH Influences the Epigenome and Epitranscriptome. Cell Chem. Biol. 2020, 27, 538–550.e7. [Google Scholar] [CrossRef]
  16. Mohammad, N.; Wong, D.; Lum, A.; Lin, J.; Ho, J.; Lee, C.-H.; Yip, S. Characterisation of isocitrate dehydrogenase 1/isocitrate dehydrogenase 2 gene mutation and the d-2-hydroxyglutarate oncometabolite level in dedifferentiated chondrosarcoma. Histopathology 2020, 76, 722–730. [Google Scholar] [CrossRef]
  17. Ghiam, A.F.; Cairns, R.A.; Thoms, J.; Dal Pra, A.; Ahmed, O.; Meng, A.; Mak, T.W.; Bristow, R.G. IDH mutation status in prostate cancer. Oncogene 2012, 31, 3826. [Google Scholar] [CrossRef] [Green Version]
  18. Cohen, A.L.; Holmen, S.L.; Colman, H. IDH1 and IDH2 mutations in gliomas. Curr. Neurol. Neurosci. Rep. 2013, 13, 345. [Google Scholar] [CrossRef] [Green Version]
  19. Matteo, D.A.; Grunseth, A.J.; Gonzalez, E.R.; Anselmo, S.L.; Kennedy, M.A.; Moman, P.; Scott, D.A.; Hoang, A.; Sohl, C.D. Molecular mechanisms of isocitrate dehydrogenase 1 (IDH1) mutations identified in tumors: The role of size and hydrophobicity at residue 132 on catalytic efficiency. J. Biol. Chem. 2017, 292, 7971–7983. [Google Scholar] [CrossRef] [Green Version]
  20. Horbinski, C. What do we know about IDH1/2 mutations so far, and how do we use it? Acta Neuropathol. 2013, 125, 621–636. [Google Scholar] [CrossRef] [Green Version]
  21. Guo, C.; Pirozzi, C.J.; Lopez, G.Y.; Yan, H. Isocitrate dehydrogenase mutations in gliomas: Mechanisms, biomarkers and therapeutic target. Curr. Opin. Neurol. 2011, 24, 648–652. [Google Scholar] [CrossRef]
  22. Sasaki, M.; Knobbe, C.B.; Itsumi, M.; Elia, A.J.; Harris, I.S.; Chio, I.I.C.; Cairns, R.A.; Mccracken, S.; Wakeham, A.; Haight, J.; et al. D-2-hydroxyglutarate produced by mutant Idh1 perturbs collagen maturation and basement membrane function. Genes Dev. 2012, 26, 2038–2049. [Google Scholar] [CrossRef] [Green Version]
  23. Waitkus, M.S.; Diplas, B.H.; Yan, H. Isocitrate dehydrogenase mutations in gliomas. Neurol. Oncol. 2016, 18, 16–26. [Google Scholar] [CrossRef] [Green Version]
  24. Thamim, M.; Thirumoorthy, K. Chiral discrimination in a mutated IDH enzymatic reaction in cancer: A computational perspective. Eur. Biophys. J. 2020, 49, 549–559. [Google Scholar] [CrossRef]
  25. Krell, D.; Mulholland, P.; Frampton, A.E.; Krell, J.; Stebbing, J.; Bardella, C. IDH mutations in tumorigenesis and their potential role as novel therapeutic targets. Future Oncol. 2013, 9, 1923–1935. [Google Scholar] [CrossRef] [Green Version]
  26. Virani, S.; Virani, S.; Colacino, J.A.; Kim, J.H.; Rozek, L.S. Cancer epigenetics: A brief review. ILAR J. 2012, 53, 359–369. [Google Scholar] [CrossRef] [Green Version]
  27. Popovici-Muller, J.; Lemieux, R.M.; Artin, E.; Saunders, J.O.; Salituro, F.G.; Travins, J.; Cianchetta, G.; Cai, Z.; Zhou, D.; Cui, D.; et al. Discovery of AG-120 (Ivosidenib): A first-in-class mutant IDH1 inhibitor for the treatment of IDH1 mutant cancers. ACS Med. Chem. Lett. 2018, 9, 300–305. [Google Scholar] [CrossRef] [Green Version]
  28. DiNardo, C.D.; Hochhaus, A.; Frattini, M.G.; Yee, K.; Zander, T.; Krämer, A.; Chen, X.; Ji, Y.; Parikh, N.S.; Choi, J.; et al. A phase 1 study of IDH305 in patients with IDH1(R132)-mutant acute myeloid leukemia or myelodysplastic syndrome. J. Cancer Res. Clin. Oncol. 2022. [Google Scholar] [CrossRef]
  29. Yao, K.; Liu, P.; Liu, H.; Wei, Q.; Yang, J.; Cao, P.; Lai, Y. 3D-QSAR, molecular docking and molecular dynamics simulations study of 3-pyrimidin-4-yl-oxazolidin-2-one derivatives to explore the structure requirements of mutant IDH1 inhibitors. J. Mol. Struct. 2019, 1189, 187–202. [Google Scholar] [CrossRef]
  30. Ma, Y.-S.; Xin, R.; Yang, X.-L.; Shi, Y.; Zhang, D.-D.; Wang, H.-M.; Wang, P.-Y.; Liu, J.-B.; Chu, K.-J.; Fu, D. Paving the way for small-molecule drug discovery. Am. J. Transl. Res. 2021, 13, 853–870. [Google Scholar]
  31. Guha, R. On exploring structure-activity relationships. Methods Mol. Biol. 2013, 993, 81–94. [Google Scholar] [CrossRef] [Green Version]
  32. Zheng, Q.; Chen, Z.; Wan, H.; Tang, S.; Ye, Y.; Xu, Y.; Jiang, L.; Ding, J.; Geng, M.; Huang, M.; et al. Discovery and structure-activity-relationship study of novel imidazole cyclopropyl amine analogues for mutant isocitric dehydrogenase 1 (IDH1) inhibitors. Bioorganic Med. Chem. Lett. 2018, 28, 3808–3812. [Google Scholar] [CrossRef]
  33. Cao, H.; Zhu, G.; Sun, L.; Chen, G.; Ma, X.; Luo, X.; Zhu, J. Discovery of new small molecule inhibitors targeting isocitrate dehydrogenase 1 (IDH1) with blood-brain barrier penetration. Eur. J. Med. Chem. 2019, 183, 111694. [Google Scholar] [CrossRef]
  34. Zhao, Q.; Manning, J.R.; Sutton, J.; Costales, A.; Sendzik, M.; Shafer, C.M.; Levell, J.R.; Liu, G.; Caferro, T.; Cho, Y.S.; et al. Optimization of 3-Pyrimidin-4-yl-oxazolidin-2-ones as Orally Bioavailable and Brain Penetrant Mutant IDH1 Inhibitors. ACS Med. Chem. Lett. 2018, 9, 746–751. [Google Scholar] [CrossRef]
  35. Ou-Yang, S.; Lu, J.; Kong, X.; Liang, Z.; Luo, C.; Jiang, H. Computational drug discovery. Acta Pharmacol. Sin. 2012, 33, 1131–1140. [Google Scholar] [CrossRef] [Green Version]
  36. Verma, J.; Khedkar, V.M.; Coutinho, E.C. 3D-QSAR in drug design-a review. Curr. Top. Med. Chem. 2010, 10, 95–115. [Google Scholar] [CrossRef]
  37. Perkins, R.; Fang, H.; Tong, W.; Welsh, W.J. Quantitative structure-activity relationship methods: Perspectives on drug discovery and toxicology. Environ. Toxicol. Chem. An Int. J. 2003, 22, 1666–1679. [Google Scholar] [CrossRef]
  38. Ferreira, L.G.; Dos Santos, R.N.; Oliva, G.; Andricopulo, A.D. Molecular docking and structure-based drug design strategies. Molecules 2015, 20, 13384–13421. [Google Scholar] [CrossRef] [Green Version]
  39. De Vivo, M.; Masetti, M.; Bottegoni, G.; Cavalli, A. Role of molecular dynamics and related methods in drug discovery. J. Med. Chem. 2016, 59, 4035–4061. [Google Scholar] [CrossRef]
  40. Krämer, S.D.; Wunderli-Allenspach, H. Physicochemical properties in pharmacokinetic lead optimization. Farmaco 2001, 56, 145–148. [Google Scholar] [CrossRef]
  41. Gupta, P.; Sharma, A.; Garg, P.; Roy, N. QSAR study of curcumine derivatives as HIV-1 integrase inhibitors. Curr. Comput. Aided. Drug Des. 2013, 9, 141–150. [Google Scholar] [CrossRef]
  42. Schrödinger LCC. Schrödinger Release 2022-3: Field-based QSAR; Schrödinger LLC: New York, NY, USA, 2021; Available online: https://www.schrodinger.com/products/field-based-qsar (accessed on 27 September 2022).
  43. Alexander, D.L.J.; Tropsha, A.; Winkler, D.A. Beware of R 2: Simple, unambiguous assessment of the prediction accuracy of QSAR and QSPR models. J. Chem. Inf. Model. 2015, 55, 1316–1322. [Google Scholar] [CrossRef] [Green Version]
  44. Bajusz, D.; Rácz, A.; Héberger, K. Why is Tanimoto index an appropriate choice for fingerprint-based similarity calculations? J. Cheminform. 2015, 7, 20. [Google Scholar] [CrossRef] [Green Version]
  45. Duan, J.; Dixon, S.L.; Lowrie, J.F.; Sherman, W. Analysis and comparison of 2D fingerprints: Insights into database screening performance using eight fingerprint methods. J. Mol. Graph. Model. 2010, 29, 157–170. [Google Scholar] [CrossRef]
  46. Salam, N.K.; Nuti, R.; Sherman, W. Novel method for generating structure-based pharmacophores using energetic analysis. J. Chem. Inf. Model. 2009, 49, 2356–2368. [Google Scholar] [CrossRef]
  47. Pinzi, L.; Rastelli, G. Molecular docking: Shifting paradigms in drug discovery. Int. J. Mol. Sci. 2019, 20, 4331. [Google Scholar] [CrossRef] [Green Version]
  48. Schrödinger LLC. Schrödinger Release 2022-3: LigPrep; Schrödinger LLC: New York, NY, USA, 2021; Available online: https://www.schrodinger.com/products/ligprep (accessed on 27 September 2022).
  49. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef]
  50. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef] [Green Version]
  51. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef] [Green Version]
  52. Hadni, H.; Elhallaoui, M. 3D-QSAR, docking and ADMET properties of aurone analogues as antimalarial agents. Heliyon 2020, 6, e03580. [Google Scholar] [CrossRef]
  53. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2012, 64, 4–17. [Google Scholar] [CrossRef]
  55. Cassidy, J.; Newell, D.R.; Wedge, S.R.; Cummings, J. Pharmacokinetics of high molecular weight agents. Cancer Surv. 1993, 17, 315–341. [Google Scholar] [PubMed]
  56. Schrödinger LLC. Schrödinger Release 2022-3: Maestro; Schrödinger LLC: New York, NY, USA, 2021; Available online: https://www.schrodinger.com/products/maestro (accessed on 27 September 2022).
  57. Banks, J.L.; Beard, H.S.; Cao, Y.; Cho, A.E.; Damm, W.; Farid, R.; Felts, A.K.; Halgren, T.A.; Mainz, D.T.; Maple, J.R.; et al. Integrated modeling program, applied chemical theory (IMPACT). J. Comput. Chem. 2005, 26, 1752–1780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. Aided. Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  59. Shelley, J.C.; Cholleti, A.; Frye, L.L.; Greenwood, J.R.; Timlin, M.R.; Uchimaya, M. Epik: A software program for pK a prediction and protonation state generation for drug-like molecules. J. Comput. Aided. Mol. Des. 2007, 21, 681–691. [Google Scholar] [CrossRef]
  60. Jacobson, M.P.; Pincus, D.L.; Rapp, C.S.; Day, T.J.F.; Honig, B.; Shaw, D.E.; Friesner, R.A. A hierarchical approach to all-atom protein loop prediction. Proteins Struct. Funct. Bioinform. 2004, 55, 351–367. [Google Scholar] [CrossRef] [Green Version]
  61. Lu, C.; Wu, C.; Ghoreishi, D.; Chen, W.; Wang, L.; Damm, W.; Ross, G.A.; Dahlgren, M.K.; Russell, E.; Von Bargen, C.D.; et al. OPLS4: Improving force field accuracy on challenging regimes of chemical space. J. Chem. Theory Comput. 2021, 17, 4291–4300. [Google Scholar] [CrossRef]
  62. BIOVIA Discovery Studio. Biovia Discovery Studio 2021 Client; Dassault Systèmes: San Diego, CA, USA, 2021. [Google Scholar]
  63. Gutiérrez, I.S.; Lin, F.-Y.; Vanommeslaeghe, K.; Lemkul, J.A.; Armacost, K.A.; Brooks III, C.L.; MacKerell, A.D., Jr. Parametrization of halogen bonds in the CHARMM general force field: Improved treatment of ligand--protein interactions. Bioorganic Med. Chem. 2016, 24, 4812–4825. [Google Scholar] [CrossRef] [Green Version]
  64. Zoete, V.; Cuendet, M.A.; Grosdidier, A.; Michielin, O. SwissParam: A fast force field generation tool for small 150 organic molecules. J. Comput. Chem 2011, 32, 2359–2368. [Google Scholar] [CrossRef]
  65. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  66. Zhurkin, V.B.; Ulyanov, N.B.; Gorin, A.A.; Jernigan, R.L. Static and statistical bending of DNA evaluated by Monte Carlo simulations. Proc. Natl. Acad. Sci. USA 1991, 88, 7046–7050. [Google Scholar] [CrossRef] [Green Version]
  67. Berendsen, H.J.C.; van Postma, J.P.M.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  68. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  69. Nosé, S.; Klein, M.L. Constant pressure molecular dynamics for molecular systems. Mol. Phys. 1983, 50, 1055–1076. [Google Scholar] [CrossRef]
Molecules 28 02315 sch001 550
Scheme 1. The mechanistic pathway of the wild-type and mutant-type IDH1/2 reaction in oncometabolite formation (Blue—allowed reaction, Red—forbidden reaction).
Scheme 1. The mechanistic pathway of the wild-type and mutant-type IDH1/2 reaction in oncometabolite formation (Blue—allowed reaction, Red—forbidden reaction).
Molecules 28 02315 sch001
Molecules 28 02315 g001 550
Figure 1. The aligned structures of (A) 62 reported ligands using Tanimoto similarities and the FB-3D-QSAR modeled results of favorable and unfavorable contour maps for (B) steric (Green—Positive, Yellow—Negative), (C) electrostatic (Blue—Positive, Red—Negative), (D) hydrophobic (Green—Positive, Red—Negative), (E) hydrogen bond acceptor (Cyan—Positive, Pink—Negative) and (F) donor (Blue—Positive, Red—Negative) properties.
Figure 1. The aligned structures of (A) 62 reported ligands using Tanimoto similarities and the FB-3D-QSAR modeled results of favorable and unfavorable contour maps for (B) steric (Green—Positive, Yellow—Negative), (C) electrostatic (Blue—Positive, Red—Negative), (D) hydrophobic (Green—Positive, Red—Negative), (E) hydrogen bond acceptor (Cyan—Positive, Pink—Negative) and (F) donor (Blue—Positive, Red—Negative) properties.
Molecules 28 02315 g001
Molecules 28 02315 g002 550
Figure 2. The steric (S and S′) and hydrophobic (H and H′) modification of favorable and non-favorable functional points in reference compound (RC-01). The green represents the oxazolidine moiety, pink represents the pyrimidine moiety, and cyan represents the imidazole moiety.
Figure 2. The steric (S and S′) and hydrophobic (H and H′) modification of favorable and non-favorable functional points in reference compound (RC-01). The green represents the oxazolidine moiety, pink represents the pyrimidine moiety, and cyan represents the imidazole moiety.
Molecules 28 02315 g002
Molecules 28 02315 g003 550
Figure 3. The cumulative picture representation of reference compound (stick model) and designed compounds (ball and stick model).
Figure 3. The cumulative picture representation of reference compound (stick model) and designed compounds (ball and stick model).
Molecules 28 02315 g003
Molecules 28 02315 g004 550
Figure 4. The overall RMSD analysis of protein and ligands compounds for 50 ns.
Figure 4. The overall RMSD analysis of protein and ligands compounds for 50 ns.
Molecules 28 02315 g004
Molecules 28 02315 g005 550
Figure 5. The overall RMSF analysis of designed and reference compounds.
Figure 5. The overall RMSF analysis of designed and reference compounds.
Molecules 28 02315 g005
Molecules 28 02315 g006 550
Figure 6. The overall Rg plot for designed and reference compounds for 50 ns.
Figure 6. The overall Rg plot for designed and reference compounds for 50 ns.
Molecules 28 02315 g006
Molecules 28 02315 g007 550
Figure 7. The HB interaction plot of the designed and reference compounds for 50 ns.
Figure 7. The HB interaction plot of the designed and reference compounds for 50 ns.
Molecules 28 02315 g007
Molecules 28 02315 g008a 550Molecules 28 02315 g008b 550
Figure 8. The receptor–ligand 3D interaction of reference compound (A) (RC-01), (B) (RC-02), and the designed compounds (C) (DC-01), (D) (DC-02), (E) (DC-03) and (F) (DC-04) after the simulation at the final frame (50 ns). Green—Classical and non-classical hydrogen bond interactions, Cyan—halogen induced hydrogen bond interactions, Lavender—alkyl-alkyl and pi-alkyl interactions, Magenta—pi-pi interactions.
Figure 8. The receptor–ligand 3D interaction of reference compound (A) (RC-01), (B) (RC-02), and the designed compounds (C) (DC-01), (D) (DC-02), (E) (DC-03) and (F) (DC-04) after the simulation at the final frame (50 ns). Green—Classical and non-classical hydrogen bond interactions, Cyan—halogen induced hydrogen bond interactions, Lavender—alkyl-alkyl and pi-alkyl interactions, Magenta—pi-pi interactions.
Molecules 28 02315 g008aMolecules 28 02315 g008b
Table
Table 1. Test set and training set compounds with their parameters using FB-3D-QSAR method.
Table 1. Test set and training set compounds with their parameters using FB-3D-QSAR method.
S.NOReported CompoundsQSAR SetReported Biological
Activity		Predicted Biological Activity
RC-01Molecules 28 02315 i001Test9.0007.960
RC-02Molecules 28 02315 i002Training8.3988.068
RC-03Molecules 28 02315 i003Training8.6788.208
RC-04Molecules 28 02315 i004Test8.6748.379
RC-05Molecules 28 02315 i005Training8.6588.634
RC-06Molecules 28 02315 i006Training8.6388.113
RC-07Molecules 28 02315 i007Training8.5698.573
RC-08Molecules 28 02315 i008Test8.3987.945
RC-09Molecules 28 02315 i009Test8.3017.414
RC-10Molecules 28 02315 i010Training8.0978.054
RC-11Molecules 28 02315 i011Training7.8868.105
RC-12Molecules 28 02315 i012Training7.8247.850
RC-13Molecules 28 02315 i013Test7.8247.511
RC-14Molecules 28 02315 i014Test7.8108.263
RC-15Molecules 28 02315 i015Training7.7357.873
RC-16Molecules 28 02315 i016Training7.7177.687
RC-17Molecules 28 02315 i017Training7.6588.097
RC-18Molecules 28 02315 i018Training7.6237.902
RC-19Molecules 28 02315 i019Training7.6207.487
RC-20Molecules 28 02315 i020Test7.6187.527
RC-21Molecules 28 02315 i021Test7.6027.205
RC-22Molecules 28 02315 i022Training7.5857.636
RC-23Molecules 28 02315 i023Training7.5697.410
RC-24Molecules 28 02315 i024Training7.5237.335
RC-25Molecules 28 02315 i025Test7.4956.904
RC-26Molecules 28 02315 i026Training7.4957.505
RC-27Molecules 28 02315 i027Training7.4767.554
RC-28Molecules 28 02315 i028Training7.4717.599
RC-29Molecules 28 02315 i029Training7.4567.358
RC-30Molecules 28 02315 i030Training7.4447.325
RC-31Molecules 28 02315 i031Training7.4447.481
RC-32Molecules 28 02315 i032Training7.4097.420
RC-33Molecules 28 02315 i033Training7.3987.231
RC-34Molecules 28 02315 i034Training7.3877.371
RC-35Molecules 28 02315 i035Training7.3777.439
RC-36Molecules 28 02315 i036Training7.3777.248
RC-37Molecules 28 02315 i037Training7.3707.590
RC-38Molecules 28 02315 i038Training7.3547.414
RC-39Molecules 28 02315 i039Training7.3287.366
RC-40Molecules 28 02315 i040Training7.3217.482
RC-41Molecules 28 02315 i041Training7.3107.202
RC-42Molecules 28 02315 i042Training7.2947.447
RC-43Molecules 28 02315 i043Training7.2687.351
RC-44Molecules 28 02315 i044Test7.2527.044
RC-45Molecules 28 02315 i045Training7.2087.281
RC-46Molecules 28 02315 i046Training7.1877.457
RC-47Molecules 28 02315 i047Training7.0716.963
RC-48Molecules 28 02315 i048Training6.9519.967
RC-49Molecules 28 02315 i049Training6.9226.937
RC-50Molecules 28 02315 i050Training6.8916.956
RC-51Molecules 28 02315 i051Training6.8186.872
RC-52Molecules 28 02315 i052Test6.6586.398
RC-53Molecules 28 02315 i053Test6.6208.192
RC-54Molecules 28 02315 i054Test6.6026.843
RC-55Molecules 28 02315 i055Test6.5757.394
RC-56Molecules 28 02315 i056Training6.3206.070
RC-57Molecules 28 02315 i057Training6.2526.442
RC-58Molecules 28 02315 i058Test6.2156.964
RC-59Molecules 28 02315 i059Training6.1426.224
RC-60Molecules 28 02315 i060Training6.0865.996
RC-61Molecules 28 02315 i061Test6.0807.028
RC-62Molecules 28 02315 i062Test5.0005.909
Table
Table 2. The selected designed and reference compounds along with their biological activity and docking scores.
Table 2. The selected designed and reference compounds along with their biological activity and docking scores.
S.NOActive LigandsBiological ActivityDocking Score
kcal/mol		Glide Score
kcal/mol
RC-01Molecules 28 02315 i0639.000 a−11.800−12.988
RC-02Molecules 28 02315 i0648.398 a−12.403−12.497
DC-01Molecules 28 02315 i0659.272 b−13.336−13.424
DC-02Molecules 28 02315 i0669.180 b−13.175−13.270
DC-03Molecules 28 02315 i0679.675 b−13.159−13.256
DC-04Molecules 28 02315 i0689.065 b−13.131−13.218
DC-05Molecules 28 02315 i0699.154 b−12.952−13.045
a—Experimental value. b—Predicted value using FB-3D-QSAR.
Table
Table 3. The ligand–receptor interactions of designed and reference compounds after the MD simulation at the final frame.
Table 3. The ligand–receptor interactions of designed and reference compounds after the MD simulation at the final frame.
S.NOActive LigandsHydrophobic InteractionsHydrogen Bond InteractionsTotal Number of Interactions
alkyl–alkyl and pi–alkylpi–pipi–sigmapi–sulfurClassicalNon-Classical
(Carbon–Hydrogen)		Halogen
1RC-01ARG-109, ALA-111, ILE-113, ILE-130, VAL-276, MET-291NilNilNilILE-128, SER-278ILE-128Nil12
2RC-02ALA-111, ILE-128, ILE-130, VAL-255, TRP-267TRP-124NilMET-291LEU-120, ILE-128NilNil13
3DC-01ALA-111, ILE-117, TRP-124, PRO-127, ILE-128, VAL-255, ALA-258, MET-290, MET-291TRP-124, TYR-285NilNilALA-111, ILE-128ALA-111, ILE-112, VAL-255CYS-11420
4DC-02ALA-111, ILE-113, CYS-114, LEU-120, TRP-124, PRO-127, ILE-128, ILE-130, MET-259, MET-290TRP-124, TRP-267, TYR-285NilNilLEU-120, ILE-128, TYR-285NilNil20
5DC-03ARG-109, ALA-111, ILE-117, TRP-124, ILE-128, ILE-130, VAL-281, ALA-282, TYR-285, MET-290, MET-291NilALA-111MET-291ILE-128, SER-287,ALA-111, ILE-112, SER-278Nil24
6DC-04ARG-109, ALA-111, ILE-113, TRP-124, ILE-128, ILE-130, VAL-255, MET-259TRP-124ALA-111NilALA-111, ILE-112, ILE-128ALA-111Nil18
Table
Table 4. ADME properties of the designed drugs along with reference drugs.
Table 4. ADME properties of the designed drugs along with reference drugs.
S.NOActive CompoundsLipinski’s Rule of Five
PSA
Å		MW
g/mol		Log P
o/w		NRBHBAHBD
1.RC-0185.17486.493.73681
2.RC-0285.17468.493.73681
3.DC-0189.66546.543.93692
4.DC-02109.89544.553.57693
5.DC-0389.66542.583.77682
6.DC-04109.89544.553.42693
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Thamim, M.; Agrahari, A.K.; Gupta, P.; Thirumoorthy, K. Rational Computational Approaches in Drug Discovery: Potential Inhibitors for Allosteric Regulation of Mutant Isocitrate Dehydrogenase-1 Enzyme in Cancers. Molecules 2023, 28, 2315. https://doi.org/10.3390/molecules28052315

AMA Style

Thamim M, Agrahari AK, Gupta P, Thirumoorthy K. Rational Computational Approaches in Drug Discovery: Potential Inhibitors for Allosteric Regulation of Mutant Isocitrate Dehydrogenase-1 Enzyme in Cancers. Molecules. 2023; 28(5):2315. https://doi.org/10.3390/molecules28052315

Chicago/Turabian Style

Thamim, Masthan, Ashish Kumar Agrahari, Pawan Gupta, and Krishnan Thirumoorthy. 2023. "Rational Computational Approaches in Drug Discovery: Potential Inhibitors for Allosteric Regulation of Mutant Isocitrate Dehydrogenase-1 Enzyme in Cancers" Molecules 28, no. 5: 2315. https://doi.org/10.3390/molecules28052315

APA Style

Thamim, M., Agrahari, A. K., Gupta, P., & Thirumoorthy, K. (2023). Rational Computational Approaches in Drug Discovery: Potential Inhibitors for Allosteric Regulation of Mutant Isocitrate Dehydrogenase-1 Enzyme in Cancers. Molecules, 28(5), 2315. https://doi.org/10.3390/molecules28052315

Article Metrics

Back to TopTop