Next Article in Journal
2-Methyloxolane as a Bio-Based Solvent for Green Extraction of Aromas from Hops (Humulus lupulus L.)
Previous Article in Journal
Preoperative and Intraoperative Methods of Parathyroid Gland Localization and the Diagnosis of Parathyroid Adenomas
Previous Article in Special Issue
Multiclass Classifier for P-Glycoprotein Substrates, Inhibitors, and Non-Active Compounds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Studies Aimed at Finding FLT3 Inhibitors and a Promising Compound and Molecular Pattern with Dual Aurora B/FLT3 Activity

by
Ítalo Antônio Fernandes
1,
Déborah Braga Resende
2,
Teodorico Castro Ramalho
1,3,
Kamil Kuca
3,* and
Elaine Fontes Ferreira da Cunha
1,*
1
Department of Chemistry, Federal University of Lavras, P.O. Box 3037, Lavras-MG 37200-000, Brazil
2
Department of Veterinary Medicine, Federal University of Lavras, P.O. Box 3037, Lavras-MG 37200-000, Brazil
3
Faculty of Science, Department of Chemistry, University of Hradec Kralove, Rokitanskeho 62, 500 03 Hradec Kralove, Czech Republic
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(7), 1726; https://doi.org/10.3390/molecules25071726
Submission received: 2 February 2020 / Revised: 5 April 2020 / Accepted: 6 April 2020 / Published: 9 April 2020
(This article belongs to the Special Issue Receptor-Dependent QSAR Methods)

Abstract

:
FLT3 and dual Aurora B/FLT3 inhibitors have shown relevance in the search for promising new anticancer compounds, mainly for acute myeloid leukemia (AML). This study was designed to investigate the interactions between human FLT3 in the kinase domain with several indolin−2-one derivatives, structurally similar to Sunitinib. Molegro Virtual Docker (MVD) software was utilized in docking analyses. The predicted model of the training group, considering nineteen amino acid residues, performed in Chemoface, achieved an R2 of 0.82, suggesting that the binding conformations of the ligands with FLT3 are reasonable, and the data can be used to predict the interaction energy of other FLT3 inhibitors with similar molecular patterns. The MolDock Score for energy for compound 1 showed more stable interaction energy (–233.25 kcal mol−1) than the other inhibitors studied, while Sunitinib presented as one of the least stable (–160.94 kcal mol−1). Compounds IAF70, IAF72, IAF75, IAF80, IAF84, and IAF88 can be highlighted as promising derivatives for synthesis and biological evaluation against FLT3. Furthermore, IAF79 can be considered to be a promising dual Aurora B/FLT3 inhibitor, and its molecular pattern can be exploited synthetically to search for new indolin−2-one derivatives that may become drugs used in the treatment of cancers, including AML.

1. Introduction

Feline McDonough Sarcoma (FMS)-like tyrosine kinase 3 (FLT3, EC: 2.7.10.1) is a class III receptor protein tyrosine kinase, wherein five immunoglobulin-like domains are present in the extracellular region [1]. Under normal conditions, this enzyme is expressed in the membranes of precursor hematopoietic cells and is important for cellular differentiation, proliferation, and multiplication. In acute myeloid leukemia (AML), blasts express FLT3, and deregulation of its signaling pathway leads to exacerbated cellular proliferation, either by hyperstimulation or mutations, both of which contribute to AML onset [2,3].
The FLT3 gene is highly affected in this disease, with mutations being observed in about 40% of AML cases [2]. Justamembranar domain mutations lead to a FLT3/ITD (internal tandem duplication) mutant allele, which has been associated with a decrease in leukemia remission of patients undergoing chemotherapy with a low response to standard cytotoxic agents [4]. Evidence shows that mutations in a justamembranar region interfere with the self-inhibitory mechanism, leading to a conformational modification and constitutive activation of this enzyme [1]. Bavetsias and Moore [5,6] showed that dual Aurora/FLT3 inhibitors might be more effective than selective FLT3 inhibitors in cases of AML mutations.
In the search for more efficient chemotherapic agents, FLT3/ITD inhibitors [7], such as Sunitinib (IC50 = 34.0 nM), Quizartinib (IC50 = 1.1 nM) [8], Midostaurin (IC50 = 9.3 nM), Lestaurtinib (IC50 = 8.6 nM), and Sorafenib (IC50 = 18.5 nM) [9] (Figure 1) have been discovered, and although presenting tolerability and toxicity problems, this has also inspired new studies in the search for new promising anticancer agents.
We reported on a study between indolin−2-one derivatives and human Aurora B kinase through molecular docking [10]. Now, in this work, we apply the same methodology looking forward to new potent FLT3 and Dual Aurora B/FLT3 inhibitors. Docking studies were performed in order to understand the binding mode of indolin−2-one derivatives, previously reported by Chern et al. [11], inside kinase domain FLT3. This was done to aid in the search for new Sunitinib analogues that may be more efficient and promising for the treatment of different cancer types, including AML.

2. Results

Crystal coordinates of the FLT3 enzyme Quizartinib and the crystallographic water molecules were taken from the Protein Data Bank (PDB) [12]. At the interaction site, Quizartinib was redocked, and the conformation (Figure 2A) with the most appropriate spatial arrangement and most favorable energy parameters was selected (RMSD = 0.99Å) and applied as a reference compound for our docking studies. The analyzed site in docking studies was defined as a subset region around the center of the Quizartinib chemical structure in the human FLT3 enzyme, and the amino acid residues are shown in Figure 2B.
Compounds 141 and Sunitinib (Table 1) were docked into the FLT3 binding site using the MVD. The ligands and amino acid residues close to 5Å from the interaction site (Figure 2B) were considered flexible during the docking simulations. Thus, a candidate solution was encoded by an array of real-valued numbers representing the ligand position, conformation, and orientation as cartesian coordinates for ligand translation, whereas there are four variables specifying the ligand orientation and one angle for each flexible torsion angle in the ligand. MolDock Score energy values were obtained from each of the selected poses for all evaluated compounds, including Quizartinib and Sunitinib, and the values are expressed in kcal mol−1 (Table S1); energy values for pose–protein interactions and hydrogen bonds were also determined. The interaction modes of the ligand with the interaction site were determined to represent the lowest energy scored protein–ligand complex used during docking, and the conformers of each compound were mostly associated with conformations of redocked Quizartinib, which was the co-crystalized compound in then interaction site. Thus, it was used as the reference compound in our docking studies. The choice of the conformation of the ligand into active site was determined as the conformation with lowest energy scored of the protein–ligand complex used during docking and the position of each compound mostly associated with conformations of redocked Quizartinib (co-crystalized compound into interaction site). The total energy between each selected ligand conformation and the amino acid residues present at a distance of 5Å (those that exhibited some type of interaction: hydrogen bond, electrostatic or steric) for the training group was calculated and the results are shown in Table 2.
The model predicted was developed based on the training set (28 compound) and predictive potential was checked for the test set (12 compounds). The Kennard–Stone algorithm was employed to split the test dataset. Multivariate analyses such as partial least square were carried out for correlating descriptors (energy values) with observed activity (pIC50). The model was validated by validation techniques and the results are showed in Table 3.
The pIC50 values for the training and test set were computed and residual values were calculated and are reported in Table S2. The standard deviations (SD) of the residual values were 0.22 (training) and 0.31 (test). To establish outlier compounds, the residuals which were more than twice the SD of the residual of fit were observed. Our data showed that the training set has two compounds outlier and test set does not have outlier. In addition, the applicability domain of the model was assessed, using William’s plot, to identify possible atypical behavior (Figure 3). In this plot, the horizontal and vertical straight lines indicate the limits of normal values: the first for the outliers and second for influential compounds. A compound will be considered outside of the applicability domain when the leverage value is higher than the critical value of 3p/n (0.44), where p is the number of model variables plus 1 and n is the number of compounds used to develop the model. A value of 2.5 for standardized residual is commonly used as a cut-off value for accepting predictions. There are no compounds which lie outside of the applicability domain of the model.

3. Discussion

By analyzing the overlapping of compounds 141 and Sunitinib (Figure 4), the distribution pattern of the selected poses for these compounds can be verified, all following the alignment for the co-crystallized Quizartinib compound. Compound 1 showed a more stable interaction energy (−233.25 kcal mol−1) than the other inhibitors studied, while Sunitinib presented as one of the least stable (−160.94 kcal mol−1). The overlapping conformations for compound 1 and Sunitinib, compared to co-crystallized Quizartinib, can be visualized in Figure 5. Both exhibited good alignment with the reference compound, highlighting the fact that compound 1 fills, similarly, most of the positions occupied by Quizartinib.
Considering the results of Table 2, the great majority interact via hydrogen bonds, with Glu692 and Cys694 (–NH and carbonyl indolinone moieties, respectively). The following situations were observed: thirty-eight compounds with Cys694; thirty compounds plus Sunitinib with Glu692; and twenty-six compounds with both. Besides those, thirty-four compounds plus Sunitinib interacted with Cys694 or Leu616 residues through an NH–pyrrolic ring. Thus, Glu692 and Cys694 residues have an important contribution to the potency of these compounds, wherein indolin−2-one and pyrrolic ring groups are important subunits. In addition, hydrogen bond interactions were observed between the following:
i) Compounds 3, 10, 13, 15, 16, 17, 20, 21, 24, 25, 28, 30, 31, 34, 36, 38, 40 and, Asp829 through pyridinone carbonyl;
ii) Other compounds also interacted with Asp829 residue through different groups: N-imidazolidinone (37, 41); amidic carbonyl (9, 20); NH between pyridinone and indolinone groups (26, 29); and amidic carbonyl located next to the terminal benzene (12, 18, 27);
iii) Compounds 4, 10, 13, 21, 26, 30, 38 and Glu661 through NH-pyridinone;
iv) Compounds 2, 6, 8, 15 and 18 with Leu616 and 4, 11, 27 and 35 with Tyr693 and Cys695 through the –COOH terminal group;
v) Other compounds also interacted with the Tyr693 residue through different groups: carbonyl group (3 and Sunitinib) attached to the pyrrol ring; and N-pyrrolidine (compound 23); N-morpholine (compound 32);
vi) Compounds 9 and Cys828 interacted through amidic carbonyl between pyridinone and indolinone groups;
vii) Compounds 1, 5, 6, 14, 19, 22, 35 and 40 interacted with Asp829 and 12, 18, 23, 27 and 33 interacted with Glu661 through the –NH group located between two acyclic carbonyls.
As observed by Zorn et al. [12], regarding redocked Quizartinib, hydrogen bond interactions were identified through NH-diaryl urea (Figure 2B)—two interactions with Glu661 (−1.65 kcal mol−1; −0.67 kcal mol−1) and another with Asp829 (−0.93 kcal mol−1) residues.
Analyzing the amino acid residues that presented the greatest energy ranges (Table 2) for compounds 141, Quizartinib, and Sunitinib, Met665 and Glu661 can be highlighted. Sunitinib showed no interaction for both residues. For Met665, the most unfavorable energy values were observed for compounds 2 (70.68 kcal mol−1) and 7 (77.63 kcal mol−1), while for Glu661, unfavorable values were observed for compounds 18 (44.19 kcal mol−1) and 23 (32.16 kcal mol−1). These can be explained by the high structural proximity of one of the terminal moieties of these compounds, which generates high steric impairment levels as well as unfavorable electrostatic interactions. Observing the level of steric impairment, it can be verified that these are the two amino acid residues that lead the greatest steric impairment, of which Met665 is the largest, followed by Glu661. In contrast, by analyzing the most favorable values, compounds 30 (−4.64 kcal mol−1) and 41 (−4.84 kcal mol−1) for Met665 and compounds 12 (−25.33 kcal mol−1) and 33 (−21.55 kcal mol−1) for Glu661 were registered, orienting more appropriately and reflecting directly in interaction energy values, being considerably better for Glu661 due to having lower steric impairment and being a charged amino acid and thus achieving a greater number of favorable electrostatic interactions compared with Met665, an uncharged amino acid.
Another highlighted residue is Asp829, which showed a range of energy values between −8.60 kcal mol−1 (compound 2) and −28.56 kcal mol−1 (compound 5), but if Sunitinib was also considered, it would show the least favorable value (−2.84 kcal mol−1) among all compounds present in Table 1. This lower value for Sunitinib, but still negative and favorable, can be explained by the repulsion of charges between this negatively charged amino acid and the fluorine atom, despite the low steric impairment. For compounds 2 and 5, it can be noticed that compound 2 showed a considerably more unfavorable orientation than compound 5, which achieved high steric impairment and hydrogen bond absence.
For Phe691, the range of values varied from −20.81 kcal mol−1 (compound 20) to +5.56 kcal mol−1 (compound 31), which can be explained by the better overlap of the pyridinone subunit, present in both compounds, in compound 20, with the benzenic subunit of this aromatic amino acid, thereby generating a more favorable energy value for compound 22 in contrast with compound 31, which has been shown to be sterically hindered by this amino acid due to torsion in the pyridinone subunit.
The model used for the prediction of activities of FLT3 inhibitors depends on statistical significance and predictive ability. Table 3 showed the results generated, an R2 value of 0.81 (RMSEc = 0.29) was achieved from calibration. The R2 values greater than 0.8 indicates that the model is correlated and may be considered to represent the training set in the same manner. A q2 of 0.60 (RMSEcv = 0.40) reveals that the model can be a useful tool for predicting affinities of new compounds based on these structures. In respect to y-randomization, the model has R2 higher than R2rand (0.22) and R2p (0.62) higher than 0.5, thereby assuring the inexistence of chance correlation. Therefore, it can be concluded that the developed model has good performance, with suitable goodness of fit, robustness and predictivity.
The equation coefficients from the regression analysis may provide useful information: the terms with positive coefficient signs decrease the predicted potency of a compound. The same concept can be applied to terms with negative coefficient signs that increase the predicted potency. In both cases, the predicted potency values are directly proportional to the magnitude of each coefficient term of this equation. The predicted pIC50 values for training were computed and are shown in Table S2. The standard deviation (SD) of the residual values was 0.26. To establish outlier compounds, we observed which residuals were more than twice the SD of the residual of the fit. Analyses of the data showed one outlier, compound 1, but due to its better activity compared with the other compounds and its use as a molecular pattern for searching for new promising analogues, it was considered for the predictive model.
Lys644 contributes to a decrease in the interaction potency of compounds and presents the highest positive regression coefficient (0.132). This positively charged amino acid is located near the terminal subunit, wherein the majority of compounds 141 have a substituted benzene or pyridinone. In the case of Quizartinib, it is next to the isoxazole group, and Sunitinib has no interaction. It is verified that repulsive electrostatic interactions occur with positively charged atoms in the molecules, as does steric impairment. Its distance is higher when compared with other residues in the interaction site and can be observed in Figure 2A, thus generating a lower influence on interactions for the evaluated compounds. Residues that also showed positive regression coefficient values were Phe691 (0.030), Glu661 (0.025), and Met665 (0.017), which all showed interaction energy values in a wide range (Table 2).
Val624 was shown to contribute to the increase in the interaction potency of the compounds and presented the most negative regression coefficient (–0.204). This amino acid is close to the indolin−2-one and pyrrol groups (compounds 145 and Sunitinib), except in compound 22 which has a furane, and the central phenyl-benzoimidazothiazole subunit for the Quizartinib compound (Figure 2A). All of these compounds exhibited favorable hydrophobic interactions with this uncharged, non-polar amino acid residue and there was an absence of highly unfavorable steric or electrostatic interactions. Two other residues that can be pointed to with negative regression coefficient values were Glu692 (–0.176) and Leu818 (–0.146).
By structurally analyzing pyridinone, imidazolidinone, and acyclic diamide in compounds 145, it was verified that pyridinone generates a greater steric impairment between the three groups, therefore it is a group that can be structurally optimized. However, an acyclic subunit was maintained, as in the case of the diaryl urea group in Quizartinib. Furthermore, imidazolidinone was not as sterically unfavorable as pyridinone, but the acyclic diamide or diaryl urea still seem to be the most favorable, possibly due to their greater conformational freedom.
When compound 1 and Quizartinib are compared using the interaction energy values (kcal mol−1) of each amino acid residue, the widest energy variations were shown in Asp829, Cys694, Cys695, Leu616, Phe691, Phe830, Tyr693, Tyr696, Val624, and Val675. All nineteen important amino acid residues appeared in both compound 1 and Quizartinib. When comparing compound 1 and Sunitinib, since compounds 141 can all be considered structural optimizations of this one, the widest energy variations were shown in Asp829, Cys695, Cys828, Glu661, Leu616, Lys644, Met665, Phe691, Tyr696, and Val675 residues, of which amino acid residues Glu661, Lys644, and Met665 showed no interactions.
Through observations from docking studies, the same one hundred new indolinonic derivatives proposed by Fernandes et al. [10] were evaluated by applying our current predictive model, as well as carrying out observations of related parameters (MolDock scores, pose–protein interactions, and hydrogen bond energy values). Among the one hundred indolinonic derivatives, satisfactory poses were obtained in eighty-one derivatives, and all of them showed interactions in the nineteen defined amino acid residues. Six compounds (IAF70, IAF72, IAF75, IAF80, IAF84, and IAF88) can be highlighted and considered to be more promising for future synthesis and biological evaluations against the human FLT3 enzyme. These six compounds have a benzyle moiety attached to an sp2 carbon atom located between indolin−2-one and a pyrrolic ring, suggesting the filling of an important verified vacancy, leading to additional interactions that would considerably contribute to the inhibitory activity of these compounds. It is also noted that the presence of a benzyle moiety would be more important than the presence of a phenyl moiety for these compounds, as verified for the Aurora B target [10]. Furthermore, the IAF75 compound presented the most satisfactory set of parameters among the six most promising compounds, as follows: predictive activity (pIC50 = 10.24), MolDock Score (−257.27 kcal mol−1), Lig-Prot. (−254.38 kcal mol−1), and hydrogen bond (−13.95 kcal mol−1). This is shown in Table 4. Seven hydrogen bond interactions were observed for IAF75: Asn701, Asp698, Asp829, Cys694, Glu692, Gly697, and Tyr693 residues. IAF84 can also be highlighted, as it presented the lowest MolDock Score value (−261.42 kcal mol−1) and pose–protein interaction energy (−258.18 kcal mol−1), in addition to a considerable hydrogen bond energy (−11.16 kcal mol−1) for seven amino acid residues—Asn701, Asp698, Asp829, Cys694, Glu692, Leu616, and Tyr693—although it did not show the highest predicted activity.
Considering studies conducted by our research group in the search for new promising inhibitors of Aurora B kinase [10] and the current results, the IAF79 compound (Table 5) can be considered to be a promising dual Aurora B/FLT3 inhibitor aimed at the treatment of different cancer types, including myeloid leukemia, as it presents considerable predicted activity and energy values for both targets studied. As shown by its chemical structure, this compound, an analogous structural variation of compound 1 (Table 1), has an insertion of 4-ethylphenol or the p-cresol moiety attached in the sp2 carbon atom located between the indolin−2-one and pyrrolic groups, leading to additional and favorable interactions in both Aurora B and FLT3 therapeutic targets; thus, it fills an important vacancy in both enzymes. Another structural difference observed was the absence of two methyl groups attached to the pyrrole ring that generate steric hindrance to the 4-ethylphenol group inserted, which could limit its conformational freedom during interactions with the targets in question. Thus, its molecular pattern (shown in Figure 6) can be exploited synthetically at different positions, especially at X, Y, W, and Z, and may aid in the search for new indolin−2-one derivatives with dual Aurora B/FLT3 activity that may become drugs used in the treatment of different types of cancer in the future.

4. Materials and Methods

4.1. Data Set

Enzyme crystal coordinates (organism: Homo sapiens) were downloaded from the Protein Data Bank (PDB code: 4XUF) [12]. Compounds 145 [11], Quizartinib, and Sunitinib FLT3 inhibitor (Table 1) were docked into the FLT3 kinase domain using the Molegro Virtual Docker (MVD) [13,14,15]. MVD is a program for predicting the most likely conformation of how a ligand will bind to a macromolecule. The MolDock Scoring Function (MolDock Score) employed by the MVD program is regulated by a new hybrid search algorithm, called guided differential evolution. This algorithm combines the differential evolution optimization technique with a cavity prediction algorithm during the searching procedure, which allows fast and accurate recognition of potential binding modes—the poses. It is derived from the Piecewise Linear Potential (PLP), a simplified potential whose parameters are fit with protein–ligand structures and binding data scoring functions [15] and further extended in the GEMDOCK program (Generic Evolutionary Method for molecular DOCK) with a new hydrogen bond term and new charge schemes. The docking scoring function, Escore, is defined by the following energy terms:
E score   =   E inter +   E intra
where the ligand–protein interaction energy, Einter, is given by
E inter = i ligand j protein [ E PLP ( r ij ) + 332.0   q i   q j 4 r ij 2 ]
The EPLP term is a “piecewise linear potential” that uses two different set of parameters: one set for approximating the steric term between atoms (van der Waals) and another stronger potential for hydrogen bonds. The second term describes the electrostatic interactions between charged atoms. It is a Coulomb potential with a distance-dependent dielectric constant, D(r) = 4r. The numerical value of 332.0 fixes the units of the electrostatic energy to kilocalories per mole (Molegro ApS).
The Eintra terms represent the internal energy value of the ligand:
E intra = i ligand j ligand   E PLP ( r ij ) + flexiblebonds   A   [ 1 cos ( m . θ     θ 0 ) ] +   E clash
The first term, double summation, is the energy between all atom pairs in the ligand, excluding atom pairs, which are connected by two bonds. The second term is a torsional energy term, where θ is the torsional angle of the bond. The average torsional energy bond contribution is used if several torsions can be determined. The last term, Eclash, assigns a penalty of 1000 if the distance between two heavy atoms is less than 2.0 Å, punishing infeasible ligand conformations (Molegro ApS). The docking search algorithm used in MVD is based on interaction optimization techniques conducted by Darwinian Evolution Theory (evolutionary algorithms, EA). A population of individuals is exposed to competitive selection that weeds out poor solutions. Recombination and mutation are used to create new solutions [15,16,17,18,19,20,21,22].

4.2. Predicted Model

The prediction model was constructed considering the interaction energies from docking for the nineteen most relevant amino acid residues located at a distance of 5 Å from the evaluated interaction site and the biological activity against FLT3 according to Chern et al. [11]. Subsequently, Chemoface software version 1.61 was used toward the statistical analysis:
i)
Coefficient of determination (R2): is the proportion of the variance in the dependent variable that is predictable from the independent variable.
ii)
Leave-one-out cross-validation (LOOcv) correlation coefficient (q2): estimating the performance of a predictive model.
iii)
Y-Randomization (R2rand): consists of the random exchange of the independent variable values. Thus, the R2rand value must be less than the correlation coefficient of the non-randomized models.
iv)
R2p: penalizes the model for the difference between R2 of randomized models and the R2 of the non-randomized model:
R P 2 = R 2 × R 2 R rand 2
v)
Correlation coefficient of external validation set (R2pred): reflects the degree of correlation between the observed (YExp(test)) and predicted (YPred(test)) activity data of the test set:
R Pred   2 = 1 1 n     ( Y Exp ( test ) Y Pred ( test ) ) 2 1 n   ( Y Exp ( test )   Y ¯ training ) 2
where Y ¯ training is average value for the dependent variable for the training set.
vi)
Modified R2 (R2m(test)): equation determining the proximity between the observed and predicted values with the zero axis intersection:
R m 2 = r 2 ( 1 | R 2   R 0 2 | )

4.3. Searching for More Potent Indolin−2-one FLT3 Inhibitors

In our previous study, we performed a search for new promising inhibitors of human Aurora B kinase [10]. Using the same methodology and the same one hundred compounds, we looked for promising inhibitors for human FLT3, considering the orientations of compounds inside the interaction site compared with the co-crystallized compound Quizartinib. These derivatives were evaluated in silico, being submitted to docking studies in MVD, and, subsequently, a pre-established predictive model was applied. The selected compounds for the predictive model were all those that presented adequate poses, comparing to redocked Quizartinib with a satisfactory RMSD value. Thus, by getting the Quizartinib redocked conformation in the interaction site under study, the compounds were selected, considering satisfactory conformations and energy parameters. Interactions in most relevant amino acid residues occurring at a distance of 5 Å from the site of enzyme interaction were considered. By searching for potential promising inhibitors for the human FLT3 enzyme, besides the predictive model, the following energetic docking parameters were considered: MolDock Scores, protein–ligand interactions, and hydrogen bonds.

5. Conclusions

We carried out molecular docking studies in order to understand the interactions of a variety of indolin−2-one derivatives with the human FLT3 enzyme. The different substituent groups were exploited at positions 3 and 6 of the indolinone ring and are shown in Table 1. Docked structures were evaluated based on their binding energy and electrostatic and hydrogen bond interactions. Our molecular docking results, combined with experimental data for human FLT3 enzyme inhibition, suggest the presence of important empty space around the indolin−2-one moiety. Thus, more favorable substituents attached at the 3-position of the pyrrolic ring and in the sp2 carbon atom located between the indolin−2-one and pyrrolic ring groups may increase the affinity for the human FLT3 enzyme, which was verified to be of a similar form to the Aurora B kinase by Fernandes et al. [10] using the same benzyle moiety. For new IAF1IAF100 compounds, six of them could be detached (IAF70, IAF72, IAF75, IAF80, IAF84, and IAF88). These all have a benzyle group attached in the sp2 carbon atom located between the indolin−2-one and pyrrolic ring groups, and IAF75 presented the most satisfactory set of parameters, while IAF84 presented the lowest MolDock score and pose–protein interaction values. Moreover, considerable hydrogen bond energy was added to seven different amino acid residues.
The present study, together with studies by Fernandes et al. [10] on this class of indolinone compounds, seem to suggest that future obtainment of a candidate for an innovative drug for the treatment of different cancer types, including myeloid leukemia, is possible, since some compounds, particularly IAF79, have been shown to exhibit possible promising dual activity against Aurora B/FLT3 and have promising molecular patterns that can be further investigated synthetically in the search for new anticancer drug candidates to build on studies that have already been carried out.

Supplementary Materials

The following are available online; Table S1: MolDock Score, Rerank Score, Pose–Protein Interaction, and Hydrogen Bond Energy values from each selected pose for compounds 140, Quizartinib, and Sunitinib; Table S2: Measured and predicted pIC50 and residual values for the training compounds.

Author Contributions

E.F.F.d.C., T.C.R. and K.K. designed this study; Í.A.F. and D.B.R. did computational experiments; Í.A.F. and E.F.F.d.C. did data interpretation, Í.A.F., T.C.R. and E.F.F.d.C. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors thank the Brazilian agencies “Conselho Nacional de Desenvolvimento Científico e Tecnológico” (CNPq) and “Coordenação de Aperfeiçoamento de Pessoal de Nível Superior” (CAPES) for their financial support. This study was supported by the Excelence Project UHK.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berenstein, R. Class III Receptor Tyrosine Kinases in Acute Leukemia—Biological Functions and Modern Laboratory Analysis. Biomark. Insights 2015, 10, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kayser, S.; Schlenk, R.F.; Londono, M.C.; Breitenbuecher, F.; Wittke, K.; Du, J.; Groner, S.; Späth, D.; Krauter, J.; Ganser, A.; et al. Insertion of FLT3 internal tandem duplication in the tyrosine kinase domain-1 is associated with resistance to chemotherapy and inferior outcome. Blood 2009, 114, 2386–2392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Port, M.; Böttcher, M.; Thol, F.; Ganser, A.; Schlenk, R.; Wasem, J.; Neumann, A.; Pouryamout, L. Prognostic significance of FLT3 internal tandem duplication, nucleophosmin 1, and CEBPA gene mutations for acute myeloid leukemia patients with normal karyotype and younger than 60 years: A systematic review and meta-analysis. Ann. Hematol. 2014, 93, 1279–1286. [Google Scholar] [CrossRef] [PubMed]
  4. Pastore, D. FLT3 inhibitors in acute myeloid leukemia. J. Hematol. Oncol. 2018, 11, 133–144. [Google Scholar] [CrossRef]
  5. Bavetsias, V.; Crumpler, S.; Sun, C.; Avery, S.; Atrash, B.; Faisal, A.; Moore, A.S.; Kosmopoulou, M.; Brown, N.; Sheldrake, P.W.; et al. Optimization of Imidazo [4–b]pyridine-Based Kinase Inhibitors: Identification of a Dual FLT3/Aurora Kinase Inhibitor as an Orally Bioavailable Preclinical Development Candidate for the Treatment of Acute Myeloid Leukemia. J. Med. Chem. 2012, 55, 8721–8734. [Google Scholar] [CrossRef]
  6. Moore, A.S.; Faisal, A.; de Castro, G.D.; Bavetsias, V.; Sun, C.; Atrash, B.; Valenti, M.; de Haven, B.A.; Avery, S.; Mair, D.; et al. Selective FLT3 inhibition of FLT3-ITD+ acute myeloid leukaemia resulting in secondary D835Y mutation: A model for emerging clinical resistance patterns. Leukemia 2012, 26, 1462–1470. [Google Scholar] [CrossRef]
  7. Ma, F.; Liu, P.; Lei, M.; Liu, J.; Wang, H.; Zhao, S.; Hu, L. Design, synthesis and biological evaluation of indolin-2-one-based derivatives as potent, selective and efficacious inhibitors of FMS-like tyrosine kinase3 (FLT3). Eur. J. Med. Chem. 2017, 127, 72–86. [Google Scholar] [CrossRef]
  8. Zarrinkar, P.P.; Gunawardane, R.N.; Cramer, M.D.; Gardner, M.F.; Brigham, D.; Belli, B.; Karaman, M.W.; Pratz, K.W.; Pallares, G.; Chao, Q.; et al. AC220 is a uniquely potent and selective inhibitor of FLT3 for the treatment of acute myeloid leukemia (AML). Blood 2009, 114, 2984–2992. [Google Scholar] [CrossRef]
  9. Nguyen, B.; Williams, A.B.; Young, D.J.; Ma, H.; Li, L.; Levis, M.; Brown, P.; Small, D. FLT3 activating mutations display differential sensitivity to multiple tyrosine kinase inhibitors. Oncotarget 2017, 8, 10931–10944. [Google Scholar] [CrossRef] [Green Version]
  10. Fernandes, I.A.; Assis, T.M.; Rosa, I.A.; da Cunha, E.F.F. Indolin-2-one Derivatives: Theoretical Studies Aimed at Finding More Potent Aurora B Kinase Inhibitors. Lett. Drug Des. Discov. 2019, 16, 138–152. [Google Scholar] [CrossRef]
  11. Chern, J.W.; Jagtap, A.D.; Wang, H.C.; Chen, G.S. Indolin-2-one Derivatives as Protein Kinase Inhibitors. WO 2013/158373 A1, 24 October 2013. [Google Scholar]
  12. Zorn, J.A.; Wang, Q.; Fujimura, E.; Barros, T.; Kuriyan, J. Crystal structure of the FLT3 kinase domain bound to the inhibitor Quizartinib (AC220). PLoS ONE 2015, 10, e0121177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hopfinger, A.J.; Wang, S.; Tokarski, J.S.; Jin, B.; Albuquerque, M.; Madhav, P.J.; Duraiswami, C. Construction of 3D-QSAR Models Using the 4D-QSAR Analysis Formalism. J. Am. Chem. Soc. 1997, 119, 10509–10524. [Google Scholar] [CrossRef]
  14. Gehlhaar, D.K.; Verkhivker, G.M.; Rejto, P.A.; Sherman, C.J.; Fogel, D.R.; Fogel, L.J.; Freer, S.T. Molecular Recognition of the Inhibitor AG-1343 by HIV-1 Protease: Conformationally Flexible Docking by Evolutionary Programming. Chem. Biol. 1995, 2, 317–324. [Google Scholar] [CrossRef] [Green Version]
  15. Thomsen, R.; Christensen, M.H. MolDock:  A New Technique for High-Accuracy Molecular Docking. J. Med. Chem. 2006, 49, 3315–3321. [Google Scholar] [CrossRef]
  16. Guimarães, A.P.; França, T.C.C.; Ramalho, T.C.; Rennó, M.N.; da Cunha, E.F.F.; Matos, K.S.; Mancini, D.T.; Kuca, K. Docking studies and effects of syn-anti isomery of oximes derived from pyridine imidazol bicycled systems as potential human acetylcholinesterase reactivators. J. Appl. Biomed. 2011, 9, 163–171. [Google Scholar] [CrossRef] [Green Version]
  17. Ramalho, T.C.; França, T.C.C.; Renno, M.N.; Guimaraes, A.P.; da Cunha, E.F.F.; Kuca, K. Development of new acetylcholinesterase reactivators: Molecular modeling versus in vitro data. Chem. Biol. Interact. 2010, 185, 73–77. [Google Scholar] [CrossRef]
  18. Da Cunha, E.F.F.; Ramalho, T.C.; Maia, E.R.; de Alencastro, R.B. The search for new DHFR inhibitors: A review of patents, January 2001 February 2005. Expert Opin. Ther. Pat. 2005, 8, 967–986. [Google Scholar] [CrossRef]
  19. Ramalho, T.C.; da Cunha, E.F.F.; de Alencastro, R.B. A density functional study on the complexation of ethambutol with divalent cations. J. Mol. Struct. 2004, 676, 149–153. [Google Scholar] [CrossRef]
  20. Tavares, T.S.; da rocha, E.P.; Nogueira, F.G.E.; Torres, J.A.; Silva, M.C.; Kuca, K.; Ramalho, T.C. Δ-FeOOH as Support for Immobilization Peroxidase: Optimization via a Chemometric Approach. Molecules 2020, 25, 259. [Google Scholar] [CrossRef] [Green Version]
  21. Lam, S.Y.S.; Leung, A.Y.H. Overcoming Resistance to FLT3 Inhibitors in the Treatment of FLT3-Mutated AML. Int. J. Mol. Sci. 2020, 21, 1537. [Google Scholar] [CrossRef] [Green Version]
  22. Bohmer, A.; Barz, S.; Schwab, K.; Kolve, U.; Gabel, A.; Kirkpatrick, J.; Ohlenshlager, O.; Gorlach, M.; Bohmer, F.D. Modulation of FLT3 signal transduction through cytoplasmic cysteine residues indicates the potential for redox regulation. Redox Biol. 2020, 28, 101325. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemical structures of Sunitinib, Quizartinib, Sorafenib, Midostaurin, and Lestaurtinib.
Figure 1. Chemical structures of Sunitinib, Quizartinib, Sorafenib, Midostaurin, and Lestaurtinib.
Molecules 25 01726 g001
Figure 2. Quizartinib into the Feline McDonough Sarcoma (FMS)-like tyrosine kinase 3 (FLT3) interaction site. (A) The structure with yellow carbon atoms is the redocked Quizartinib. (B) Hydrogen bonds provide evidence for redocked Quizartinib and amino acid residues are located at a distance of 5Å.
Figure 2. Quizartinib into the Feline McDonough Sarcoma (FMS)-like tyrosine kinase 3 (FLT3) interaction site. (A) The structure with yellow carbon atoms is the redocked Quizartinib. (B) Hydrogen bonds provide evidence for redocked Quizartinib and amino acid residues are located at a distance of 5Å.
Molecules 25 01726 g002
Figure 3. William’s plot.
Figure 3. William’s plot.
Molecules 25 01726 g003
Figure 4. Overlapped conformations of FLT3 inhibitors. Quizartinib is shown in yellow, and Sunitinib is shown in purple.
Figure 4. Overlapped conformations of FLT3 inhibitors. Quizartinib is shown in yellow, and Sunitinib is shown in purple.
Molecules 25 01726 g004
Figure 5. Overlapped conformations of compound 1, Quizartinib (yellow), and Sunitinib (purple).
Figure 5. Overlapped conformations of compound 1, Quizartinib (yellow), and Sunitinib (purple).
Molecules 25 01726 g005
Figure 6. Promising indolin−2-one molecular pattern for dual Aurora B/FLT3 activity.
Figure 6. Promising indolin−2-one molecular pattern for dual Aurora B/FLT3 activity.
Molecules 25 01726 g006
Table 1. Chemical structures, IC50 and pIC50 values of 141 [11], Quizartinib and Sunitinib FLT3 inhibitor compounds. Test set compounds are marked with an asterisk.
Table 1. Chemical structures, IC50 and pIC50 values of 141 [11], Quizartinib and Sunitinib FLT3 inhibitor compounds. Test set compounds are marked with an asterisk.
CompoundChemical StructureIC50 (nM) FLT3pIC50 FLT3
1 Molecules 25 01726 i0010.59.301
2 Molecules 25 01726 i0021.38.886
3 Molecules 25 01726 i0031.48.854
4* Molecules 25 01726 i0041.68.796
5* Molecules 25 01726 i0052.48.620
6* Molecules 25 01726 i0062.78.569
7 Molecules 25 01726 i0072.78.569
8 Molecules 25 01726 i0082.78.569
9 Molecules 25 01726 i0092.98.538
10 Molecules 25 01726 i0102.98.538
11 Molecules 25 01726 i0113.58.456
12 Molecules 25 01726 i0124.48.357
13* Molecules 25 01726 i0134.88.319
14* Molecules 25 01726 i0146.18.215
15 Molecules 25 01726 i0157.38.137
16* Molecules 25 01726 i0168.18.092
17 Molecules 25 01726 i0178.68.066
18 Molecules 25 01726 i01810.08.000
19* Molecules 25 01726 i01910.08.000
20 Molecules 25 01726 i02010.37.987
21 Molecules 25 01726 i02111.67.936
22* Molecules 25 01726 i02215.47.812
23 Molecules 25 01726 i02316.27.790
24 Molecules 25 01726 i02424.67.609
25 Molecules 25 01726 i02524.97.604
26 Molecules 25 01726 i02627.77.558
27 Molecules 25 01726 i02731.77.499
28 Molecules 25 01726 i02834.97.457
29 Molecules 25 01726 i02937.57.426
30 Molecules 25 01726 i03037.97.421
31* Molecules 25 01726 i03138.17.419
32 Molecules 25 01726 i03239.47.404
33 Molecules 25 01726 i03345.77.340
34* Molecules 25 01726 i03448.47.315
35 Molecules 25 01726 i03554.27.266
36 Molecules 25 01726 i03664.27.192
37 Molecules 25 01726 i03774.17.130
38* Molecules 25 01726 i03886.67.062
39* Molecules 25 01726 i039115.36.938
40 Molecules 25 01726 i040148.46.828
Quizartinib Molecules 25 01726 i0414.2 a1
1.1 a2
8.377
8.959
Sunitinib Molecules 25 01726 i0429.9 a1
34.0 a2
8.004
7.468
a1: FLT3/Wild [8]; a2: FLT3/ITD [8]. * Test set compounds.
Table 2. Interaction energy values (kcal mol−1) between amino acid residues and ligand for training and test set.
Table 2. Interaction energy values (kcal mol−1) between amino acid residues and ligand for training and test set.
Cpd/aaAla
642
Asp
698
Asp
829
Cys
694
Cys
695
Cys
828
Glu
661
Glu
692
Gly
697
Leu
616
Leu
818
Lys
644
Met
665
Phe
691
Phe
830
Tyr
693
Tyr
696
Val
624
Val
675
1−1.90−1.91−28.29−12.11−3.75−18.82−10.10−4.10−10.57−23.07−11.32−1.914.22−14.37−8.25−14.65−4.19−3.53−10.21
2−3.47−2.52−8.60−11.61−3.78−8.4611.60−4.30−9.89−23.55−12.07−9.6270.68−12.02−8.44−12.14−3.83−5.74−7.10
3−3.18−2.03−14.63−10.25−8.59−8.33−21.14−1.75−10.32−21.16−13.13−1.0527.99−12.19−9.04−13.02−10.30−4.48−4.60
4−2.06−1.81−12.19−11.28−8.05−9.41−15.47−1.70−9.61−16.22−14.06−3.231.95−4.00−12.28−13.12−6.78−5.76−8.14
5−3.93−4.79−28.56−10.96−3.05−24.34−7.84−4.92−11.72−14.22−12.61−2.692.79−14.80−9.70−10.11−5.02−5.04−8.56
6−3.72−3.25−19.77−12.07−2.87−21.27−4.86−4.94−9.70−24.01−11.81−2.577.50−12.93−8.96−10.96−2.66−5.69−10.26
7−3.58−5.82−9.43−12.55−3.54−10.520.52−2.20−10.90−17.27−11.75−7.1877.63−8.56−7.08−10.75−6.27−5.58−4.75
8−4.03−5.13−15.80−12.87−2.02−23.44−2.34−5.47−8.53−20.14−12.68−1.619.50−14.83−8.84−9.66−2.13−5.05−7.23
9−5.14−2.17−11.28−8.87−3.82−6.943.47−3.14−10.13−18.14−9.58−6.5365.61−12.54−7.25−11.71−4.78−8.17−5.69
10−2.30−1.62−12.27−11.99−3.68−12.31−9.30−5.20−9.17−21.41−11.99−4.3538.66−3.53−7.94−12.27−3.17−3.58−8.62
11−2.39−1.70−13.05−12.58−9.24−10.14−16.87−2.00−9.45−11.96−14.02−2.973.85−2.70−11.56−15.15−7.67−5.44−8.52
12−3.22−9.77−17.90−10.84−4.04−11.92−25.33−3.83−15.35−15.56−12.73−4.4115.24−14.71−9.85−10.81−7.82−5.72−6.56
13−1.79−1.59−13.55−10.83−6.14−11.91−8.55−2.41−9.81−15.30−13.45−5.53−4.10−2.98−10.12−15.00−6.59−5.87−7.99
14−2.15−1.80−25.45−12.12−3.52−19.77−8.72−5.03−9.32−21.43−11.89−1.713.03−16.56−8.81−12.82−3.25−3.42−7.48
15−2.66−2.27−10.52−12.48−4.21−10.40−1.77−3.69−10.22−24.67−12.05−6.2054.98−1.02−7.62−11.22−4.02−3.64−8.03
16−2.23−1.92−13.13−13.07−3.91−10.60−1.06−4.51−10.82−17.85−11.47−5.1239.50−7.34−7.64−11.77−4.50−3.50−9.17
17−3.33-0.91−15.54−14.26−2.48−12.128.11−4.44−6.15−12.24−11.87−7.2234.52−3.82−6.98−12.19−1.25−5.27−7.94
18−3.48−5.04−20.27−11.13−2.16−12.3244.19−1.95−9.89−22.54−10.40−2.13−2.50−18.51−11.80−7.48−2.21−4.43−7.12
19−2.15−1.32−25.70−14.60−2.07−18.54−9.08−6.09−6.94−17.34−11.95−2.472.22−16.83−7.05−12.42−1.64−4.73−9.06
20−2.94−1.97−12.01−12.66−6.93−11.42−17.81−4.69−14.44−16.58−10.88-0.3133.24−20.81−8.59−11.82−14.06−2.895.41
21−1.80−1.48−10.81−11.82−7.05−11.32−12.80−1.94−6.00−19.99−13.32−4.11-0.97−3.64−11.52−6.08−3.90−6.15−7.82
22−1.70−1.31−12.20−14.35−1.81−23.16−2.52−6.30−4.04−10.48−11.75−2.028.27−13.68−6.37−10.84−1.64−5.55−6.99
23−2.59−1.66−18.36−11.64−3.54−12.6332.16−3.09−9.04−13.48−12.99−7.69−2.75−19.03−8.37−14.20−3.84−3.67−8.31
24−4.49−2.12−13.80−9.24−3.64−11.077.10−3.21−10.41−18.62−11.74−9.5627.49−10.24−8.61−11.81−4.27−6.65−7.45
25−4.44−1.16−23.81−11.81−3.31−11.652.30−2.16−4.03−10.48−12.26−8.10−1.2311.81−8.98−11.76−2.85−6.58−8.21
26−3.46−1.01−17.06−12.26−2.70−15.765.74−1.79−5.98−10.60−12.48−3.4050.76−13.89−7.20−10.54−2.08−6.03−1.60
27−3.04−1.25−21.42−10.60−3.94−11.96−21.37−4.68−8.62−15.67−11.82−4.61−1.05−14.40−9.63−16.59−3.29−6.16−7.27
28−1.88−1.12−16.49−8.57−6.61−10.78−7.95-0.93−7.19−21.41−10.72−6.684.370.54−9.77−12.12−4.89−8.58−7.27
29−2.58-0.85−19.75−12.76−3.21−10.3212.68−1.89−6.22−13.12−11.82−4.0829.51−12.09−6.24−12.57−2.05−6.41−4.56
30−1.89-0.79−14.09−13.14−2.10−11.77−6.60−5.49−5.51−5.58−12.82−4.48−4.64−8.71−11.04−11.76−1.12−4.29−7.49
31−3.00-0.67−20.12−13.76−2.70−11.78−4.00−4.92−5.82−10.84−11.69−6.47−2.565.56−7.77−12.88−1.14−5.24−8.90
32−3.59−2.33−17.04−10.53−4.09−15.48−15.18−1.60−10.97−18.08−11.57−1.906.44−16.41−8.20−9.10−5.18−7.444.76
33−2.96−2.01−14.00−11.59−3.68−12.27−21.55−3.87−10.00−19.56−12.71−3.8121.37−16.47−8.23−12.49−4.03−4.28−5.11
34−2.04−1.24−14.85−14.31−2.46−12.116.84−5.96−5.13−8.97−12.32−5.8127.07−12.89−6.44−11.51−1.79−4.35−8.87
35−2.42−1.55−25.07−11.15−9.38−18.03−8.01−2.73-0.36−13.27−12.19−1.68−2.38−15.94−7.62−16.61−6.66−4.76−4.90
36−2.07−1.14−14.32−14.53−1.87−11.2810.54−5.39−5.15−8.90−12.44−7.022.20−10.78−7.14−10.95−1.30−4.78−8.37
37−2.23−1.22−17.23−14.97−1.41−20.78−5.56−6.80−5.08−10.34−12.09−2.184.18−17.93−6.69−10.49−1.07−4.24-0.69
38−2.24−1.51−14.39−10.47−3.64−11.66−8.32−5.08−8.20−11.91−12.40−4.53−4.61−5.64−9.53−11.34−2.65−3.94−7.51
39−3.53−2.64−19.41−8.51−2.35−17.95−11.74−1.44−9.49−17.55−10.16−1.618.11−16.42−13.80−6.79−3.21−8.192.66
40−2.77-0.98−22.08−12.72−3.36−7.884.48−1.86−6.34−12.98−12.08−4.22−1.32−10.22−5.81−12.47−2.26−6.27−7.50
Quiz−4.48−1.55−18.63−8.59−7.63−19.56−11.13−1.64−10.93−15.17−9.37−4.922.19−9.99−12.34−10.03−13.10−7.47−6.01
Sunit−2.11−1.65−2.84−10.80−8.63−0.98-−3.42−9.61−19.49−12.80--−4.01−10.82−16.27−19.13−5.61−3.63
Table 3. Statistical parameters evaluated in the analysis.
Table 3. Statistical parameters evaluated in the analysis.
ParametersAccept ValuesObtained Values
R2>0.80.80
RMSEc-0.29
q2>0.50.60
RMSEcv-0.46
LV-5
R2rand< r20.22
RMSE y-rand-0.60
R2pred>0.80.80
RMSEpred-0.31
r2m(test)>0.50.68
R2p>0.50.61
R2: Coefficient of determination; RMSEc: root mean square error in calibration; q2: Leave-one-out cross-validation correlation coefficient; RMSEcv: root mean square error in validation; LV: latent variable; R2pred: Correlation coefficient of external validation; r2m(test): equation 3; R2rand: Y-randomization; RMSE y-rand: root mean square error in Y-randomization; R2p: equation 1.
Table 4. Chemical structures, predicted pIC50 values, MolDock Scores (kcal mol−1), pose-protein interaction (kcal mol−1) and hydrogen bond interactions (kcal mol−1) for the most promising derivatives against FLT3.
Table 4. Chemical structures, predicted pIC50 values, MolDock Scores (kcal mol−1), pose-protein interaction (kcal mol−1) and hydrogen bond interactions (kcal mol−1) for the most promising derivatives against FLT3.
CompdChemical StructurepIC50 predMolDock ScoreInter EnergyH Bond Energy
IAF70 Molecules 25 01726 i04310.06−236.96−249.26−9.37
IAF72 Molecules 25 01726 i04410.25−238.17−251.73−6.74
IAF75 Molecules 25 01726 i04510.24−257.27−254.38−13.95
IAF80 Molecules 25 01726 i04610.22−255.22−252.98−11.75
IAF84 Molecules 25 01726 i04710.13−261.42−258.18−11.16
IAF88 Molecules 25 01726 i04810.12−254.29−250.82−8.30
Table 5. Chemical structure, predicted pIC50, MolDock Scores (kcal mol−1) of ligand-protein interactions (kcal mol−1) and hydrogen bond interactions (kcal mol−1) for a promising dual Aurora B/FLT3 inhibitor.
Table 5. Chemical structure, predicted pIC50, MolDock Scores (kcal mol−1) of ligand-protein interactions (kcal mol−1) and hydrogen bond interactions (kcal mol−1) for a promising dual Aurora B/FLT3 inhibitor.
CompdChemical StructureEnzymepIC50 predMolDock ScoreInter EnergyH Bond Energy
IAF79 Molecules 25 01726 i049Aurora B11.39 b−236.71 b−240.59 b−15.61 b
FLT39.83−248.06−246.77−8.35
b [10].

Share and Cite

MDPI and ACS Style

Fernandes, Í.A.; Braga Resende, D.; Ramalho, T.C.; Kuca, K.; da Cunha, E.F.F. Theoretical Studies Aimed at Finding FLT3 Inhibitors and a Promising Compound and Molecular Pattern with Dual Aurora B/FLT3 Activity. Molecules 2020, 25, 1726. https://doi.org/10.3390/molecules25071726

AMA Style

Fernandes ÍA, Braga Resende D, Ramalho TC, Kuca K, da Cunha EFF. Theoretical Studies Aimed at Finding FLT3 Inhibitors and a Promising Compound and Molecular Pattern with Dual Aurora B/FLT3 Activity. Molecules. 2020; 25(7):1726. https://doi.org/10.3390/molecules25071726

Chicago/Turabian Style

Fernandes, Ítalo Antônio, Déborah Braga Resende, Teodorico Castro Ramalho, Kamil Kuca, and Elaine Fontes Ferreira da Cunha. 2020. "Theoretical Studies Aimed at Finding FLT3 Inhibitors and a Promising Compound and Molecular Pattern with Dual Aurora B/FLT3 Activity" Molecules 25, no. 7: 1726. https://doi.org/10.3390/molecules25071726

APA Style

Fernandes, Í. A., Braga Resende, D., Ramalho, T. C., Kuca, K., & da Cunha, E. F. F. (2020). Theoretical Studies Aimed at Finding FLT3 Inhibitors and a Promising Compound and Molecular Pattern with Dual Aurora B/FLT3 Activity. Molecules, 25(7), 1726. https://doi.org/10.3390/molecules25071726

Article Metrics

Back to TopTop