Next Article in Journal
Customizing a Cognitive Stimulation Program for Individuals with Dementia through a Participatory Design Approach
Next Article in Special Issue
Modification of 6,7-Dichloro-5,8-Quinolinedione at C2 Position: Synthesis, Quantum Chemical Properties, and Activity against DT-Diaphorase Enzyme
Previous Article in Journal
Goal-Conditioned Reinforcement Learning within a Human-Robot Disassembly Environment
Previous Article in Special Issue
A Look at the Importance of Chirality in Drug Activity: Some Significative Examples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probing New Antileukemia Agents That Target FLT3 and BCL-2 from Traditional Concoctions through a Combination of Mass Spectrometry Analysis and Consensus Docking Methods

by
Adebayo A. Adeniyi
1,2,3,*,
Joy Nkechinyere Adeniyi
4,
Manimbulu Nlooto
5 and
Parvesh Singh
2,*
1
Department of Industrial Chemistry, Federal University Oye Ekiti, Oye 371104, Nigeria
2
School of Chemistry and Physics, University of KwaZulu-Natal, P/Bag X54001, Westville, Durban 4000, South Africa
3
Department of Chemistry, University of the Free State, P.O. Box 339, Bloemfontein 9300, South Africa
4
Traditional Medicine Laboratory, School of Nursing and Public Health, College of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
5
Department of Pharmacy, School of Health Care Sciences, University of Limpopo, Sovenga 0727, South Africa
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(22), 11611; https://doi.org/10.3390/app122211611
Submission received: 27 September 2022 / Revised: 7 November 2022 / Accepted: 7 November 2022 / Published: 15 November 2022

Abstract

:
The search for new chemotherapeutics against leukemia is of great interest to researchers, owing to the limitation of the current drugs. In this research, new drug candidates against leukemia were probed through liquid chromatography-mass spectrometer (LC-MS) analysis of three traditional herbal concoctions, that provide the phytochemical profile of the samples. The identified compounds from the LC-MS were modeled for the analysis of their antileukemia activities, by using five different consensus methods, to combine the seven docking scores. The consensus methods are used to combine the docking scores to avoid losing promising drug candidates, due to a poor reproducibility of the docking scores across the different packages, due to differences in the scoring functions and training sets across the docking packages. The libraries of the potential drug candidates from the concoctions were constructed by searching the NIST database for molecules with a similar MS fragmentation. Venetoclax and gilteritinib, that target FLT3 and BCL-2 were ranked among the top hits, indicating the efficiency of this protocol without missing any potential drug. The results ranked rescinnamine and bisacodyl as new potential antileukemia agents that targets FLAT3, and BCL-2, including the mutated BCL-2 G101V receptor, that is known to be resistant to treatment with venetoclax.

1. Introduction

The traditional experimental method for screening a large number of molecules towards biological targets is a high-throughput screening technique, but its demand for a large number of resources and time lead to the introduction of computer-aided virtual screening methods to decrease the time and economic costs [1]. The docking programs are used to filter potentially active compounds towards a protein target from a large library of compounds. These docking programs are computationally efficient, fast and of great impact in drug discovery, to search for the favorable position, orientation, and assumed conformation of molecules upon the binding and assigning docking scores to estimate the binding affinity [1,2,3,4,5]. The major challenge of relying on docking methods to screen a library of molecules against the protein target, is the possibility of providing false binding and ending up eliminating even the potent drug candidates from the list of top hits that will be considered. The reasons for the failure of the docking scores in some instances, are related to their training-set dependencies and scoring-function parameterization [1]. The approach that has been found successful in limiting the possibility of throwing potential drugs into the trash bin of nonactive compounds and achieving a higher success rate, is to combine the results of multiple docking scores and then apply consensus ranking methods [1,6,7]. A higher success rate in the virtual screening processes has been obtained by combining results from different docking programs (i.e., consensus scoring) to obtain a final rank [1]. In addition to many traditional consensus methods that have been introduced, which make use of the average scores or the ranking of each molecule from different programs, a new exponential consensus ranking (ECR) was introduced recently that promises a good performance where the traditional consensus fails, due to the poor performance of one or more of the docking programs [1]. This study also employs the consensus docking approach in probing potential antileukemia agents, in combination with LC-MS. The three leukemia targets that are selected in this study are the Fms-like tyrosine kinase 3 (FLT3), B-cell lymphoma-2 (BCL-2), and the mutated BCL-2 G101V.
The abnormalities in FLT3, such as the over-expression of the gene or FLT3 mutations from internal tandem duplications (ITDs) within the FLT3 gene are linked to the poor prognosis factors in acute myeloblastic leukemia (AML) and acute lymphoblastic leukemia (ALL) [8,9,10]. FLT3 is considered a promising drug target in AML patients with FLT3 mutations [11]. The antileukemia drug that targets FLT3 is gilteritinib, which is known for the treatment of patients with relapsed/refractory FLT3-mutated AML [11,12]. It is classified as an ATP-competitive type I inhibitor and a selective FLT3 inhibitor that inhibits both FLT3-ITD and FLT3-TKD mutations [11,13].
BCL-2 is a founding member of the pro-survival class of proteins and functions by interacting with homology 3 (BH3) motifs of pro-apoptotic molecules to restrain its pro-apoptotic function [14]. The overexpression of BCL-2 proteins helps the cancer cells to evade treatment, therefore an effective step toward the treatment of cancers is the down-regulation of the anti-apoptotic BCL-2 proteins [15]. Venetoclax is the first BCL-2 antagonist approved, in 2016, for cancer treatment that is targeting the binding groove of BCL-2, that is usually used to interact with BH3 [14,16]. The binding of venetoclax to BCL-2 will prevent the BCL-2 BH3 protein-protein interaction and then promote the apoptosis in malignant cells, dependent on BCL-2 [14]. The first novel BCL-2 mutation G101V has been found in patients that initially responded to venetoclax treatment, but developed chronic lymphocytic leukemia (CLL) [14,17]. The BCL-2 G101V mutation provides resistance to the therapy by leading to some 180-fold reduction in the binding affinity of venetoclax and yet maintains an affinity for the BH3 motif of pro-apoptotic proteins [14].
This study considers the inhibitory activities of the proposed ligands from three samples of traditional concoctions that are named Emelia M, Mshikazi and Delosma H decoctions [18], against three leukemia targets BCL-2 (PDB 6o0k) and mutated BCL-2 G101V (PDB: 6o0l). To efficiently probe the inhibitory activities of the potential drug candidates, seven scoring functions and five consensus methods were used to combine the results. The library of potential drug candidates was generated from three traditional herbal samples that are locally recommended by traditional healers to possess anticancer agents. The potential bioactive compounds in the three traditional concoctions were predicted, based on their experimental MS fragmentations that match the MS fragmentation of molecules found in the NIST database. The combination of LC-MS and docking has been previously used by Mohammed et al. (2022), for the phytochemical profiling of Suaeda vermiculata extracts [19]. The authors made use of the molecular docking of the compounds identified in the LC-MS chemical profiling against three ATP-binding cassette proteins to identify the active anti-hepatocarcinoma compounds [19] In this work, we instead targeted the key fragments present in the three herbal samples using LC-MS and multiple docking instead of single molecular docking. This presents a fast and efficient method of analyzing traditional concoctions without necessarily tampering with the non-disclosure agreements made with the traditional healers as the source of their products is often kept secret to avoid knowledge theft. The study also included venetoclax and gilteritinib, which are known antileukemia drugs that target BCL-2 and FLAT3 (PDB: 6il3), respectively. The two drug candidates serve as a benchmark for the docking scores and consensus methods on their ability to rank active compounds among the top hits. The interest is to use the co-crystalized ligands (venetoclax and gilteritinib) as the benchmark for the scoring functions and consensus methods to provide the correct binding site orientation and to rank the active compounds among the top 10 lead hits.

2. Methodology

The method employed in this work is best summarized as a LS-M analysis of the samples, followed by the identification of the ligands from the spectra, the docking of the identified ligands to the target, and the application of the consensus method to the docking score.

2.1. Generating Potential Drug Candidates from Traditional Herbs

Three samples of traditional herb medicine concoctions (labeled as A, B, and C) that are suggested to have antileukemia activities were collected from traditional healers. In this study, the samples were analyzed using LC-MS (ABSCIEX 4000 QTRAP hybrid triple quadrupole mass spectrometer with Shimadzu’s front end). Twenty microliters of each sample were injected onto a Discovery C18 reverse-phase column (150 × 2.1 mm, Supelco) and separated using 0% to 95% linear water (Solvent A) and methanol (Solvent B) gradient over 30 min at 0.3 mL/min. The eluting analytes were first analyzed in the positive and then the negative ionization mode, each time, using an information-dependent acquisition (IDA) method where the ions between 200 and 1000 Da with intensities above 100,000 counts per second (cps) originating from an enhanced MS (EMS) survey scan, were selected and fragmented in the collision cell, and the fragments were recorded following an enhanced production (EPI) scan. The obtained LC-MS fragmentation of the samples was further analyzed to identify drug candidates that have a similar fragmentation from the NIST database (NIST 17 MS) that best matches the fragmentation, using the automatic mass spectral deconvolution and identification system (AMDIS 2.73) package that is designed for automating the complex process of extracting data from GC/MS data. It is worth pointing out that the identified molecules are not the exact molecules from the concoctions but the best representative target molecules from the NIST database. In all three samples, a total of 271 unique molecules (excluding all repetition) from the MS retention time, ranging from 0.0806 to 44.8792 min, were identified in the NIST database. The structures of the target molecules were downloaded from the NIST database, using their CAS numbering and were optimized, using the semi-empirical method PM7 for their virtual screening as potential inhibitors of FLAT3, BCL-2, and mutated BCL-2 G101V.

2.2. Docking and Scoring-Function

The downloaded 271 target structures from the NIST database were first optimized with the G16 package [20], using the semi-empirical method PM7 before their virtual screening, using different docking programs. Following the exclusion of the molecules with inorganic metals, such as silicon, from the list, a total of 262 molecules were screened against the three targets. A total of seven docking scores from five docking programs were used to analyze the inhibitory activity of the 262 identified molecules. The CAS number and the names of all of the 262 final molecules that were docked against the three target receptors are added to the Supplementary Table S1. In the four docking software packages that were used (Autodock, Molegro, Haddock, and Schrodinger), the binding site residues of the proteins were first identified and the ligands binding was restricted to the grid area of the selected binding site residues. The first two docking scores (Moledock and Rerank scores) are from the Molegro package [21], the third score (Vina score) is from Autodock vina [22], the fourth score (Glide score) is from Schrodinger package [23,24], and the last three are from the Haddock docking package [25,26], named according to the docking criteria (rigid Haddock(rigid)), flexible refinement (Haddock(flex)) and flexible refinement in water (Haddock(flex, water)). The default settings were used in all docking simulations with little modifications. The search area distance of 25 Angstroms from the selected point on the binding site was used in Autodock-vina. The parameterization of the ligands for the Haddock docking was prepared using Acpype packages [27,28] and we chose Gasteiger charges for all of the atoms. The haddock was set to generate 50 rigid structures, 10 flexible refined poses, and finally 10 final flexible refined poses in water media. The Glide docking was carried out, using Schrodinger Python version 68013, using the default settings for the grid generation and docking.

2.3. Consensus Ranking of the Combined Docking Scores

A total of five different consensus methods were used to combine the results of the seven docking scores and to generate a better ranking of the top lead hits. The five consensus methods that were used are:
Rank by number (RbN) that combine the docking score of each molecule (sj) using the average values over the number of the docking programs (n) [1,29].
R b N i = 1 n j s i j
Rank by rank (RbR) makes use of the rank of each molecule (rj) and averages it over all of the docking programs (n) for each of the molecules [1,29].
R b R i = 1 n j r i j
Average of auto-scaled scores (AASS), that first normalizes each score between 0 and 1, using the minimum (sjmin) and maximum (sjmax) scores per each program and then averages over the number of the docking programs [1,30].
A A S S i = 1 n j s i j s m i n j s m a x j s m i n j
Z-score (ZS) that scales the score of each molecule (sj), using the average (μj) and standard deviation (σj) of the scores per each program and averages over the number of the docking programs used [1,31].
Z S i = 1 n j s i j μ j σ j
Exponential Consensus ranking (ECR) makes use of the exponential score of the rank of each molecule (rj) and the expected value of the exponential distribution (σ) and them sums over all of the docking programs (n) for each of the molecules [1].
E C R i = 1 σ j exp r i j σ

3. Results and Discussion

3.1. Identification of the Drug Candidates from the Concoctions Using LC-MS

Working with traditional concoctions is very challenging for the scientist because most of the traditional herbal practitioners withhold the names of the plants and other materials that were used in the preparation and the method of preparation is concealed for fear that someone else can steal their ideas. In this study, we demonstrate the fast means of identifying the biologically active compounds from the traditional concoctions without tampering with the non-disclosure agreements made with the traditional healers.
In this study, we analyzed three traditional concoctions collected from traditional healers in the Kwazulu-Natal province of South Africa, using LC-MS fragmentation. The identification of the closely similar drug candidates that best represent the LC-MS fragmentation was carried out using the automatic mass spectral deconvolution and identification system (AMDIS) package to analyze the fragmentation and the closely related molecules were identified from the NIST database. The three samples were analyzed to identify the potential antileukemia compounds, because of the unsubstantiated claims by the traditional healers on the anticancer activity of their products. A total of 271 molecules were generated from the MS fragmentation of the three traditional herbs samples that were analyzed with the AMDIS program and found to have related target structures in the NIST database that match best with the LC-MS fragmentation. A typical example of the LC-MS fragmentation that was used to identify 865043, which was found among the active compounds in this study, is shown in Figure 1.

3.2. Docking Study

A total number of 262 molecules from the 270 identified molecules were used for the docking study after excluding those that contained inorganic metals, such as silicon. The docking study was conducted using the docking programs Molgro, Autodock Vina, Schrodinger Glide, and the Haddock programs. A total of seven docking scores were generated, two from Molegro (Moledock and Rerank scores), Vina, Glide, and the last three from Haddock, according to the scores obtained for the rigid, flexible, and flexible with water refinement, as described in the methodology.
The 262 potential drug candidates from the traditional herbs were screened for their antileukemia properties against two major leukemia targets BCL-2 and FLT3. The study also included the G101V mutant of BCL-2 (BCL-2-G101V) that has been recently reported to display a resistance to therapy [14]. The docking study included the co-crystallized ligands venetoclax [14,16] and gilteritinib [11,12] which are known antileukemia drugs against the target BCL-2 and FLT3, respectively, to serve as a benchmark for our study. The benchmark is used to determine the ability of the docking programs to effectively rank the experimentally proven inhibitor, among the top 10 lead hits and also the possibility of reproducing the experimental orientation. The results of the seven docking scores are analyzed using five consensus ranking methods RbN, RbR, AAS, ZS, and ECR, as described in the methodology.

3.3. Optimizing Ranking to Select the Top Hit Molecules

The list of the top 10 hit molecules from each of the seven docking scores used, is shown in Table 1, while those of the five consensus methods that were used to combine the seven scores are shown in Table 2. One of the most significant applications of the docking programs is to use it to screen the large library of molecules for the top hits that can bind strongly to the target. This has been of great use in the drug discovery research and has led to a significant decrease in the time and cost of discovering new drugs, compared to the traditional high-throughput screening. However, the challenge with the docking scores is that they are significantly different, based on their training set and parametrization method, which consequently leads to a different ranking of the screened molecules. This implies that no single docking score is completely reliable as it is possible to obtain a false hit and lose the promising molecules in the process.

3.4. Performance of the Docking Scores with Known Drugs

The ranking of known antileukemia drugs venetoclax and gilteritinib that target BCL-2 and FLAT3 are examined. The geometries of BCL-2 (PDB: 6o0k) and the mutated BCL-2 G101V (PDB: 6o0l) have a co-crystallized structure of venetoclax binding to their binding site. The results of the top 10 hits of seven docking scores are shown in Table 1. It is obvious from the results that the ranking of the molecules is significantly different from each other across the docking scores. It is obvious that gilteritinib, that is a known inhibitor of FLT3 (6il3), was ranked among the top 10 hit molecules only by the Moledock and Haddock(flex), in the absence of these two docking scores, this potent drug candidate would have been thrown into the trash bin of the non-active, by the other five docking scores. The performances of the docking scores were better for a known drug against BCL-2 (6o0k), which is venetoclax. Six of the seven docking scores, Moledock, Glide, and the three Haddock, ranked venetoclax among the best 10 hits, only the Rerank score got it wrong. The interaction of venetoclax with the mutated BCL-2 G101V is expected to be lower than the activity toward BCL-2, based on the experimental information, among the six that found venetoclax as a promising BCL-2 inhibitor, Vina, Haddock(rigid), Haddock(flex) and Haddock(flex, water) ranked its inhibitory activity towards the mutated BCL-2 G101V lower than that of BCL-2, while Moledock and Glide gave it a better inhibitory activity towards the mutated BCL-2 G101V than BCL-2.
Using the consensus method to combine the results, the ECR ranked gilteritinib among the top 10 hits while it is excluded from the list of top 10 hits obtained from other consensus methods. Four of the consensus methods RbR, AASS, ZS, and ECR ranked venetoclax among the top 10 hits and even predicted it has the best inhibitor of the mutated BCL-2 G101V, by these four consensus methods. Furthermore, the known drug gilteritinib against FLAT3, is ranked among the top best inhibitors of the mutated BCL-2 G101V, by the docking scores (Moledock, Rerank, and Haddock(flex)), and also by all of the five consensus methods.

3.5. The Distribution of the Docking Scores

To understand the distribution of the scores obtained from the seven docking scores used for the ranking of the 262 potential drug candidates, from the three samples of the traditional concoctions and the two known drugs with the level of overlap of the docking scores in the ranking of the molecules, we normalized all of the docking scores using the normalizing expression explained for the consensus AASS. The plots of the distribution of the normalized docking scores of all of the molecules across the seven docking scores showing their level of overlap, are shown in Figure 2. The distributions for each of the seven docking scores are added to the Supplementary Figure S1.
The results of the distribution show that the scores of the ligands have a relatively normal distribution in the seven docking scores. The normal distribution is an indication that the molecules are grouped into three major groups, comprised of a small number of most active molecules on the extreme left, a large number of molecules as less active in the middle, and a small number of non-active molecules on the extreme right.
The pick of the distribution from the Rerank score is found towards the left of the others and is higher than others, indicating that it ranked many of the molecules as active, more than the other docking scores. The mean value of the docking scores from Moledock, Vina, Glide, and the three Haddock scores, are closer to each other than to the Rerank scores, indicating a much more similar distribution among the six excluding, the Rerank scores (Figure 2). Furthermore, the plots of the multi-variance relationship that show the level of correlation among the seven docking scores in the ranking of the molecules are shown in Figure 3. The plots of the multivariance analysis in Figure 3 clearly show that in the three targets, the ranking of the ligands by Moledock and Haddock(flex) have a significant level of correlation which also included Haddock(flex, water), in the target BCL-2 and its mutated BCL-2 G101V. The scores from Haddock(flex) and Haddock(flex, water) are also significantly correlated and there is a relatively fair correlation between Glide and Vina.
The distribution of the five consensus methods that were used to combine the docking scores, is shown in Figure S2. The methods RbN, RbR, AASS, and ZP also follow the normal distribution, indicating a few numbers of the molecules as active compounds on the extreme left side of the distribution. The distribution using the consensus ECR does not show a normal distribution; it has a higher number of molecules on the extreme left. Since the ranking of the ECR is different from the other consensus methods, in the sense that the active molecule will have the highest score values, it also clearly shows that the active compounds are a small proportion of the molecules on the extreme right of the distribution, in relation with the small number of active molecules found on the extreme left of the other consensus methods.

3.6. Potential Inhibitors BCL-2

In addition to the known drug venetoclax, the best-ranked inhibitors of BCL-2, by at least three of the seven docking scores among the 10 top hits, are 486475, 603509, 85734, and 865043. They are all captured among the top 10 hits of the five consensus methods except 486475, which is excluded in the RbR. Selecting the top 8 hits from the ECR and the top one from the other consensus methods, as as carried out for the target FLAT3, the promising potential inhibitors of BCL-2 from the three samples of the traditional concoctions are 865043, 603509, 85734, 486475, 58208, 113155, and 81232 of which, the known drug venetoclax is ranked best of all. The most promising, among the potential BCL-2 inhibitors from the traditional concoction, are 865043 and 603509, which are ranked best by many of the consensus methods and also by the docking scores.
Those ranked as potential inhibitors of the mutated BCL-2 G101V by at least three of the docking scores among the top 10 hits, are 24815245, 50146, 604320, and 865043 from the traditional concoctions and gilteritinib, which is a known inhibitor of FLAT3. The most promising potential inhibitors from the traditional concoctions that are included in the top 8 hits obtained from the ECR and the top 1 from the other consensus methods, are 865043, 604320, 24815245, 50146, 28782425, and 84640. The known inhibitor of BCL-2 (venetoclax,) is ranked best of all and also the known inhibitor of FLAT3 (gilteritinib) is ranked third. The best ranked across the consensus methods as the best inhibitors of the mutated BCL-2 G101V, are 865043 and 24815245, which were also ranked best for the native BCL-2 and FLAT3 targets, respectively.

3.7. The Difference in the Binding Site Orientation

The binding orientation obtained for the docking of known antileukemia drugs venetoclax and gilteritinib, that target FLAT3 and BCL-2 are shown in Figure 4 for the six docking scores (the seventh docking score Rank has the same binding site orientation as Moledock). Just as there is a significant difference in the ranking of the molecules, there is also a significant difference in their binding site orientation. In addition to the influence of the different score parametrization and training set, that leads to a difference in the docking scores [1], our result further shows that the difference in the docking poses plays a significant role in the observed disparity in the ranking of the molecules. Concerning the crystal structure binding orientation, none of the docking scores predicts the crystal structure binding orientation of gilteritinib to FLAT3. The binding orientation of venetoclax to BCL-2 was predicted by only Vina as the best pose, while the other scores predicted different orientations (Figure 4b).

3.8. The Top-Ranked Molecules from the Traditional Concoctions

The results of the analysis of the seven docking scores and five consensus methods show that the two potential inhibitors of FLAT3, from the traditional concoction are 24815245 (rescinnamine) and 101757 (benzenamine, N-phenyl-4-(phenylazo)). The two most prominent inhibitors of BCL-2 are 865043 (methoserpidine) and 603509 (bisacodyl), while that of the mutated BCL-2 G101V are 865043, and 24815245, which are already found to be the best inhibitors of BCL-2 and FLAT3, respectively. The two molecules methoserpidine (865043) and rescinnamine (24815245) that are found to be the most promising potential inhibitor of the mutated BCL-2 G101V, are very similar in structure with just a minor difference (Figure 5), but they are significantly different in function and toxicity, as explained below. The structures of these four molecules are shown in Figure 5. These molecules have been considered for biological applications, especially rescinnamine (that targets FLAT3 and BCL-2 G101V in this study) and bisacodyl (that targets BCL-2 in this study) are known drugs in clinical applications. Rescinnamine is a vinca alkaloid that is derived from Rauvolfia serpentina and is used as an antihypertensive drug [32,33]. Bisacodyl is a stimulant laxative drug that is being used to treat constipation and found to be both safe and effective for long-term use [34,35]. Benzenamine, N-phenyl-4-(phenylazo), which can also be called benzeneazodiphenylamine, is part of the diphenylamine derivatives that are well known to act as an antioxidant, but not considered seriously for use in edible products for the fear of toxicity [36]. Methoserpidine is a synthetic isomer of reserpine that is being used as an antihypertensive drug but found to have a side effect that results in mental depression [37,38].
The binding site residues interactions of rescinnamine with FLAT3 and the mutated BCL-2 G101V, and the interaction of bisacodyl with BCL-2 are shown in Figure 6 in an effort to compare their residue interactions with the residue interactions that defined the activity of the gilteritnib interaction with FLAT3 and the venetoclax interaction with BCL-2 and the mutated BCL-2 G101V, using their crystal interaction structures. The crystal interaction of venetoclax with BCL-2 and the mutated BCL-2 G101V are very similar in terms of protein residue interaction with only a minor difference that involves the Arg146 interaction in BCL-2, but Phe 153 interaction in the mutated BCL-2 G101V. The interaction of bisacodyl occupies some of the binding pockets that defined the activity of venetoclax, in its interaction with BCL-2, as they both share many residue interactions in common, such as Phe 104, Arg146, Tyr 108, Leu 137, Glue 136, Met 115, and Asp 111. The interaction of rescinnamine with the mutated BCL-2 G101V is also very similar to what was observed for the interaction of bisacodyl with BCL-2. The rescinnamine interaction with the mutated BCL-2 G101V was found to occupy a greater portion of the binding pockets that define the activity of venetoclax, as they shared many common residue interactions, such as Tyr 202, Val 148, Ala 100, Tyr 108, Phe 104, Asp 111, and Leu 137. The rescinnamine interaction with FLAT3 also shared some similar residue interactions that define the activity of the gilteritinib interaction with FLAT3, such as residues Cys 695, Gly 697, Lys 614, and Asp 698.

4. Conclusions

A total of 262 potential antileukemia candidates were identified from the three samples, using LC-MS and were screened against leukemia targets FLAT3, BCL-2, and the mutated BCL-2 G101V, after excluding those that have inorganic metals, such as silicon.
The results from the ECR offered a better ranking of the potential drug candidates without excluding the benchmark molecules venetoclax and gilteritinib, which are known antileukemia drugs that target BCL-2 and FLAT3, respectively.
Among the top hit inhibitors from both the docking and consensus results, are 24815245 (rescinnamine) and 101757 (benzenamine, N-phenyl-4-(phenylazo)) that target FLAT3, as well as 865043 (methoserpidine) and 603509 (bisacodyl), that target BCL-2.
Methoserpidine and rescinnamine, which were found to be potential inhibitors of BCL-2 and FLAT3, respectively, were also predicted to be the top two promising inhibitors of the mutated BCL-2 G101V.
The binding interactions of the top hit potential drug candidates from the traditional concoctions are demonstrated to have similar residue interactions as the known drugs that target FLAT3 and BCL-2. This is an indication that the proposed potential drugs explored have similar binding pockets that define the activity of the known inhibitors of these targets. A different study on the experimental anticancer activities of the concoctions using leukemia cancer cell lines is also under consideration.

Supplementary Materials

The following supporting information can be download at: https://www.mdpi.com/article/10.3390/app122211611/s1.

Author Contributions

Conceptualization, A.A.A.; Formal analysis, J.N.A.; Investigation, A.A.A. and M.N.; Project administration, M.N.; Supervision, P.S.; Writing—original draft, A.A.A.; Writing—review & editing, P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Foundation (Grant No. 121276).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All of the necessary data that can help in the reproducibility of the results are included in the manuscript and also in the Supplementary Materials. The experimental LC-MS raw files were analyzed using the automatic mass spectral deconvolution and identification system (AMDIS) package version 2.73. The docking study was accomplished using Moledock version 2013.6.0, Autodock vina version 1.1.2, Glide in Schrodinger package version 2020, and Haddock docking package version 2.2. In all software used, the default parameters were used and all of the modifications are clearly stated in the methodology. Most of the data analysis and plotting are accomplished using Gnuplot version 5.2.

Acknowledgments

The authors would like to acknowledge the University of Kwazulu-Natal financial support and the Centre for High-Performance Computing (CHPC), Cape Town (South Africa), for the supercomputing facilities. P.S. thanks the National Research Foundation (NRF) South Africa for a Competitive Grant for unrated Researchers (Grant No. 121276).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Palacio-Rodríguez, K.; Lans, I.; Cavasotto, C.N.; Cossio, P. Exponential consensus ranking improves the outcome in docking and receptor ensemble docking. Sci. Rep. 2019, 9, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jorgensen, W.L. The Many Roles of Computation in Drug Discovery. Science 2004, 303, 1813–1818. [Google Scholar] [CrossRef] [PubMed]
  3. Rognan, D. The impact of in silico screening in the discovery of novel and safer drug candidates. Pharmacol. Ther. 2017, 175, 47–66. [Google Scholar] [CrossRef] [PubMed]
  4. Jorgensen, W.L. Efficient drug lead discovery and optimization. Acc. Chem. Res. 2009, 42, 724–733. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, Y.C. Beware of docking! Trends Pharmacol. Sci. 2015, 36, 78–95. [Google Scholar] [CrossRef]
  6. Charifson, P.S.; Corkery, J.J.; Murcko, M.A.; Walters, W.P. Consensus scoring: A method for obtaining improved hit rates from docking databases of three-dimensional structures into proteins. J. Med. Chem. 1999, 42, 5100–5109. [Google Scholar] [CrossRef]
  7. Ren, X.; Shi, Y.S.; Zhang, Y.; Liu, B.; Zhang, L.H.; Peng, Y.B.; Zeng, R. Novel Consensus Docking Strategy to Improve Ligand Pose Prediction. J. Chem. Inf. Model. 2018, 58, 1662–1668. [Google Scholar] [CrossRef]
  8. Perelló-Reus, C.M.; Català, A.; Caviedes-Cárdenas, L.; García-Vega, N.; Camós, M.; Pérez-Torras, S.; Pastor-Anglada, M. FMS-like tyrosine kinase 3 (FLT3) modulates key enzymes of nucleotide metabolism implicated in cytarabine responsiveness in pediatric acute leukemia. Pharmacol. Res. 2020, 151, 104556. [Google Scholar] [CrossRef]
  9. Ozeki, K.; Kiyoi, H.; Hirose, Y.; Iwai, M.; Ninomiya, M.; Kodera, Y.; Miyawaki, S.; Kuriyama, K.; Shimazaki, C.; Akiyama, H.; et al. Biologic and clinical significance of the FLT3 transcript level in acute myeloid leukemia. Blood 2004, 103, 1901–1908. [Google Scholar] [CrossRef] [Green Version]
  10. Stam, R.W.; Den Boer, M.L.; Schneider, P.; Nollau, P.; Horstmann, M.; Beverloo, H.B.; Van Der Voort, E.; Valsecchi, M.G.; De Lorenzo, P.; Sallan, S.E.; et al. Targeting FLT3 in primary MLL-gene-rearranged infant acute lymphoblastic leukemia. Blood 2005, 106, 2484–2490. [Google Scholar] [CrossRef]
  11. Kawase, T.; Nakazawa, T.; Eguchi, T.; Tsuzuki, H.; Ueno, Y.; Amano, Y.; Suzuki, T.; Mori, M.; Yoshida, T. Effect of Fms-like tyrosine kinase 3 (FLT3) ligand (FL) on antitumor activity of gilteritinib, a FLT3 inhibitor, in mice xenografted with FL-overexpressing cells. Oncotarget 2019, 10, 6111–6123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gorcea, C.M.; Burthem, J.; Tholouli, E. ASP2215 in the treatment of relapsed/refractory acute myeloid leukemia with FLT3 mutation: Background and design of the ADMIRAL trial. Futur. Oncol. 2018, 14, 1995–2004. [Google Scholar] [CrossRef] [PubMed]
  13. Mori, M.; Kaneko, N.; Ueno, Y.; Yamada, M.; Tanaka, R.; Saito, R.; Shimada, I.; Mori, K.; Kuromitsu, S. Gilteritinib, a FLT3/AXL inhibitor, shows antileukemic activity in mouse models of FLT3 mutated acute myeloid leukemia. Investig. New Drugs 2017, 35, 556–565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Birkinshaw, R.W.; Gong, J.; Luo, C.S.; Lio, D.; White, C.A.; Anderson, M.A.; Blombery, P.; Lessene, G.; Majewski, I.J.; Thijssen, R.; et al. Structures of BCL-2 in complex with venetoclax reveal the molecular basis of resistance mutations. Nat. Commun. 2019, 10, 2385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Liu, T.; Wu, Z.; He, Y.; Xiao, Y.; Xia, C. Single and dual target inhibitors based on Bcl-2: Promising anti-tumor agents for cancer therapy. Eur. J. Med. Chem. 2020, 201, 112446. [Google Scholar] [CrossRef]
  16. Roberts, A.W.; Davids, M.S.; Pagel, J.M.; Kahl, B.S.; Puvvada, S.D.; Gerecitano, J.F.; Kipps, T.J.; Anderson, M.A.; Brown, J.R.; Gressick, L.; et al. Targeting BCL2 with venetoclax in relapsed chronic lymphocytic leukemia. N. Engl. J. Med. 2016, 374, 311–322. [Google Scholar] [CrossRef]
  17. Blombery, P.; Anderson, M.A.; Gong, J.; Thijssen, R.; Birkinshaw, R.W.; Thompson, E.; Teh, C.E.; Nguyen, T.; Xu, Z.; Flensburg, C.; et al. Acquisition of the Recurrent Gly101Val Mutation in BCL2 Confers Resistance to Venetoclax in Patients with Progressive Chronic Lymphocytic Leukemia. Blood 2018, 132, LBA-7. [Google Scholar] [CrossRef]
  18. Nkechinyere Adeniyi, J.; Nlooto, M.; Ngcobo, M.; Moodley, R.; Gomo, E. Phytochemical profile and in vitro antioxidant activity of Emelia M (EMB), Mshikazi and Delosma H herbal medicines as demonstrated in THP-1 and Jurkat leukaemia cell lines. Afr. Health Sci. 2021, 21, 1924–1937. [Google Scholar] [CrossRef]
  19. Mohammed, H.A.; Almahmoud, S.A.; Arfeen, M.; Srivastava, A.; El-Readi, M.Z.; Ragab, E.A.; Shehata, S.M.; Mohammed, S.A.A.; Mostafa, E.M.; El-khawaga, H.A.; et al. Phytochemical profiling, molecular docking, and in vitro anti-hepatocellular carcinoid bioactivity of Suaeda vermiculata extracts. Arab. J. Chem. 2022, 15, 103950. [Google Scholar] [CrossRef]
  20. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, X.; et al. Gaussian 09, Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  21. Thomsan, R.; Christensen, M.H. MolDock: A new technique for high-accuracy molecular docking. J. Med. Chem. 2006, 49, 3315. [Google Scholar] [CrossRef]
  22. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [Green Version]
  24. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A New Approach for Rapid, Accurate Docking and Scoring 2. Enrichment Factors in Database Screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef] [PubMed]
  25. Dominguez, C.; Boelens, R.; Bonvin, A.M.J.J. HADDOCK: A protein-protein docking approach based on biochemical or biophysical information. J. Am. Chem. Soc. 2003, 125, 1731–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. De Vries, S.J.; Van Dijk, A.D.J.; Krzeminski, M.; Van Dijk, M.; Thureau, A.; Hsu, V.; Wassenaar, T.; Bonvin, A.M.J.J. HADDOCK versus HADDOCK: New features and performance of HADDOCK2.0 on the CAPRI targets. Proteins Struct. Funct. Genet. 2007, 69, 726–733. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic atom type and bond type perception in molecular mechanical calculations. J. Mol. Graph. Model. 2006, 25, 247–260. [Google Scholar] [CrossRef]
  28. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general Amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  29. Wang, R.; Wang, S. How Does Consensus Scoring Work for Virtual Library Screening? An Idealized Computer Experiment. J. Chem. Inf. Comput. Sci. 2001, 41, 1422–1426. [Google Scholar] [CrossRef]
  30. Oda, A.; Tsuchida, K.; Takakura, T.; Yamaotsu, N.; Hirono, S. Comparison of consensus scoring strategies for evaluating computational models of protein-ligand complexes. J. Chem. Inf. Model. 2006, 46, 380–391. [Google Scholar] [CrossRef]
  31. Liu, S.; Fu, R.; Zhou, L.H.; Chen, S.P. Application of consensus scoring and principal component analysis for virtual screening against β-secretase (BACE-1). PLoS One 2012, 7, e38086. [Google Scholar] [CrossRef]
  32. Smirk, F.H.; Mcqueen, E.G. Comparison of Rescinnamine and Reserpine As Hypotensive Agents. Lancet 1955, 266, 115–116. [Google Scholar] [CrossRef]
  33. Fife, R.; Maclaurin, J.C.; Wright, J.H. Rescinnamine in treatment of hypertension in hospital clinic and in general practice. Br. Med. J. 1960, 2, 1848–1850. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Stiens, S.A.; Luttrel, W.; Binard, J.E. Polyethylene glycol versus vegetable oil based bisacodyl suppositories to initiate side-lying bowel care: A clinical trial in persons with spinal cord injury. Spinal Cord 1998, 36, 777–781. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wald, A. Constipation advances in diagnosis and treatment. JAMA J. Am. Med. Assoc. 2016, 315, 185–191. [Google Scholar] [CrossRef]
  36. Bickoff, E.M.; Livingston, A.L.; Guggolz, J.; Ray Thompson, C. Comparative evaluation of antioxidants for carotene. J. Am. Oil Chem. Soc. 1952, 29, 445–446. [Google Scholar] [CrossRef]
  37. Holt, M.C. Methoserpidine. Br. Med. J. 1961, 2, 415–418. [Google Scholar] [CrossRef]
  38. Forrester, T.M.; Shirriffs, G.G. Methyldopa and Methoserpidine in Hypertension. Lancet 1967, 289, 141–142. [Google Scholar] [CrossRef]
Figure 1. The (a) chromatogram, (b) scan, (c) extracted and (d) suggested target methoserpidine (CAS number: 865043) mass spectra derived from one of the three samples.
Figure 1. The (a) chromatogram, (b) scan, (c) extracted and (d) suggested target methoserpidine (CAS number: 865043) mass spectra derived from one of the three samples.
Applsci 12 11611 g001
Figure 2. The distribution of the normalized docking scores of all of the molecules across the seven docking scores, showing the level of overlap and mean value (dash line).
Figure 2. The distribution of the normalized docking scores of all of the molecules across the seven docking scores, showing the level of overlap and mean value (dash line).
Applsci 12 11611 g002
Figure 3. The distribution of the individual seven docking scores along with the two principal components with the quality of the contribution (cos2) of each to the PCA computed from the multivariance analysis to show their level of correlation. When the docking scores appear within close range to each other and are of a similar contribution to the PCA (similar color), then it shows a good level of correlation.
Figure 3. The distribution of the individual seven docking scores along with the two principal components with the quality of the contribution (cos2) of each to the PCA computed from the multivariance analysis to show their level of correlation. When the docking scores appear within close range to each other and are of a similar contribution to the PCA (similar color), then it shows a good level of correlation.
Applsci 12 11611 g003
Figure 4. The binding site orientation of (a) gilteritinib to FLAT3 and (b) venetoclax to BCL-2 obtained from Molegro (green), Vina (cyan), Glide (magenta), Haddock(rigid) (yellow), Haddock(flex) (blue), and Haddock(flex, water) (red), in relation to their crystal structure crystal (orange).
Figure 4. The binding site orientation of (a) gilteritinib to FLAT3 and (b) venetoclax to BCL-2 obtained from Molegro (green), Vina (cyan), Glide (magenta), Haddock(rigid) (yellow), Haddock(flex) (blue), and Haddock(flex, water) (red), in relation to their crystal structure crystal (orange).
Applsci 12 11611 g004
Figure 5. The known antileukemia drugs venetoclax and gliteritinib that target BCL-2 and FLAT3, respectively, and the four most promising antileukemia candidates predicted from the docking study.
Figure 5. The known antileukemia drugs venetoclax and gliteritinib that target BCL-2 and FLAT3, respectively, and the four most promising antileukemia candidates predicted from the docking study.
Applsci 12 11611 g005
Figure 6. The docking pose binding site interaction of the rescinnamine (24815245) with FLAT 3 (6il3) and the mutated BCL-2 G101V (6o0l) and the interaction of bisacodyl (603509) with BCL-2 (6o0k), comparing the residue interaction with that of the crystal structure of the gilteritinib interaction with FLAT3 and the venetoclax interaction with BCL-2 and the mutated BCL-2 G101V.
Figure 6. The docking pose binding site interaction of the rescinnamine (24815245) with FLAT 3 (6il3) and the mutated BCL-2 G101V (6o0l) and the interaction of bisacodyl (603509) with BCL-2 (6o0k), comparing the residue interaction with that of the crystal structure of the gilteritinib interaction with FLAT3 and the venetoclax interaction with BCL-2 and the mutated BCL-2 G101V.
Applsci 12 11611 g006
Table 1. The top 10 hit molecules from the seven docking scores (Moledock, Rerank, Vina, Glide, Haddock(rigid), Haddock(flex) and Haddock(flex, water)) ranked according to their binding. The molecules from the traditional herbs are named, according to the CAS number that was used to extract them from the NIST database. The names of all of the molecules are added to Table S1 of the Supplementary Materials.
Table 1. The top 10 hit molecules from the seven docking scores (Moledock, Rerank, Vina, Glide, Haddock(rigid), Haddock(flex) and Haddock(flex, water)) ranked according to their binding. The molecules from the traditional herbs are named, according to the CAS number that was used to extract them from the NIST database. The names of all of the molecules are added to Table S1 of the Supplementary Materials.
MoleculeMoledockMoleculeRerankMoleculeVinaMoleculeGlideMoleculeHaddock (Rigid)MoleculeHaddock (Flex)MoleculeHaddock (Flex, Water)
6o0k
1989522−137.6281232−102.43venetoclax−12.10venetoclax−7.63venetoclax−21.34venetoclax−57.22venetoclax−75.17
113155−137.0257874−102.08215587−10.203342981−6.832642800−19.92gilteritinib−52.12865043−54.66
24815245−130.773648213−98.1671636−10.2067199660−6.7850293−19.60865043−51.33486475−53.25
7235407−130.72865043−96.37113155−10.106267023−6.713848519−18.8118713514−47.9878308−53.24
venetoclax−129.0158208−93.331662017−9.60215587−6.707173844−18.34603509−46.7085734−53.07
67970−128.311241947−88.363342981−9.1028322023−6.68614266−18.02629970−46.4418713514−51.34
486475−127.40101757−87.6977098−9.1072480−6.6422232719−17.83486475−45.08604320−51.02
58208−126.2985734−86.17213467−9.001771182−6.61603509−17.7178308−44.3632634687−50.77
71636−126.27603509−85.0756495−9.0086884−6.572227170−17.6985734−43.66603509−50.63
3648213−125.93563042−83.6081232−9.00213467−6.56604320−17.66630024−43.6352645531−50.25
6o0l
venetoclax−152.76venetoclax−110.00venetoclax−11.70venetoclax−8.45544763−19.3828782425−44.33venetoclax−55.58
24815245−146.21gilteritinib−102.5471636−10.20215587−6.69604320−18.82gilteritinib−43.56865043−48.77
1989522−142.783648213−98.85113155−10.00239010−6.6550146−17.86865043−43.091164165−47.11
604320−137.8957874−96.8081232−9.7086884−6.58531919−17.7767970−42.302432908−46.82
113155−136.4459052−93.50215587−9.5028322023−6.5439208009−17.4250146−41.8750146−46.36
7235407−134.60865043−90.17213467−9.401207698−6.47101757−16.83604320−40.96126147−45.92
gilteritinib−133.7624815245−89.951662017−9.2063077009−6.43117817−16.75venetoclax−40.0824815245−45.08
865043−132.21563042−89.8156495−9.101771182−6.3985698−16.687173844−39.51629970−44.27
85734−130.391330865−89.343342981−9.0072480−6.391120281−16.5424815245−39.2728782425−43.87
71636−129.9958208−88.92632520−8.80897789−6.392259963−16.32117817−39.0432634687−43.46
6il3
24815245 −143.83 59052 −110.66 56495 −10.50 897789 −9.48 133493 −15.78 24815245 −49.40 1617909 −49.58
1989522 −141.97 111013 −105.01 213467 −10.40 520285 −8.91 111013 −15.07 gilteritinib −45.96 59052 −47.70
venetoclax −141.78 865043 −101.04 63077009 −10.40 239010 −8.00 608935 −15.03 85723 −44.74 1166525 −46.69
111013 −139.62 55334515 −100.24 113155 −9.70 3457485 −7.97 1207698 −14.88 642659 −44.60 101757 −45.61
113155 −139.41 2432908 −98.73 1662017 −9.70 2440224 −7.82 116165 −14.64 117840 −44.09 1031932 −45.36
604320 −136.22 24815245 −94.75 129000 −9.60 101757 −7.56 56495 −14.58 215587 −43.39 563042 −44.39
865043 −132.98 83885 −93.28 3342981 −9.60 531919 −7.44 58208 −13.93 213467 −43.16 126147 −44.23
59052 −131.80 57114 −92.62 239010 −9.50 15597769 −7.37 111068 −13.78 81812 −42.76 28782425 −44.12
gilteritinib −131.66 84786 −92.60 venetoclax −9.50 40817081 −7.10 215587 −13.76 1617909 −42.63 93765 −43.45
2432908 −128.74 1843056 −91.88 40817081 −9.30 85734 −7.10 614266 −13.72 1843056 −42.61 85734 −43.03
Table 2. The top 10 hit molecules from the five consensus ranking methods that were used to combine the seven docking scores. The molecules from the traditional herbs are named according to the CAS number that was used to extract them from the NIST database. The names of all of the molecules are added to Table S1 of the Supplementary Materials.
Table 2. The top 10 hit molecules from the five consensus ranking methods that were used to combine the seven docking scores. The molecules from the traditional herbs are named according to the CAS number that was used to extract them from the NIST database. The names of all of the molecules are added to Table S1 of the Supplementary Materials.
MoleculeRbNMoleculeRbRMoleculeAASSMoleculeZ-ScoreMoleculeECR
6o0k
865043−50.83586504318.429venetoclax0.156venetoclax−2.015venetoclax0.513
85734−48.06260350920.8578650430.229865043−1.6768650430.246
58208−48.048384851934.8576035090.272603509−1.3986035090.200
633318−47.648166201737.286gilteritinib0.28085734−1.238857340.174
57874−46.8307709838.857857340.291gilteritinib−1.2254864750.174
3648213−46.6198573439.000582080.29258208−1.205582080.159
81232−46.535venetoclax39.42938485190.2981662017−1.1931131550.149
357573−46.2875820840.57116620170.3013848519−1.158812320.146
486475−45.935gilteritinib49.7145630420.316486475−1.0432155870.143
563042−45.6855264553151.5714864750.319563042−0.99516620170.138
865043−50.83586504318.429venetoclax0.156venetoclax−2.015venetoclax0.513
6o0l
24815245−56.941venetoclax6.857venetoclax0.027venetoclax−2.550venetoclax0.505
venetoclax−56.263gilteritinib18.857gilteritinib0.190gilteritinib−1.5628650430.256
gilteritinib−49.9462878242526.714248152450.20824815245−1.495gilteritinib0.241
865043−48.6132481524530.1678650430.217865043−1.4146043200.224
357573−48.304225996338.714287824250.22828782425−1.318248152450.217
3648213−46.1605264553140.000501460.24250146−1.142501460.212
28782425−45.07186504340.429679700.25367970−1.122287824250.172
1989522−44.95510175745.8575319190.27352645531−1.012846400.160
39208009−43.894116416547.857526455310.275531919−0.98236482130.147
28322023−43.3386797053.1436043200.2841164165−0.9482155870.144
24815245−56.941venetoclax6.857venetoclax0.027venetoclax−2.550venetoclax0.505
6il3
24815245 −57.671 101757 25.143 24815245 0.224 24815245 −1.458 24815245 0.230
604320 −51.421 52645531 30.857 101757 0.237 101757 −1.374 111013 0.218
59052 −49.836 56495 36.857 56495 0.248 56495 −1.336 59052 0.210
111013 −49.366 531919 37.429 63077009 0.271 63077009 −1.194 56495 0.168
67970 −47.673 642659 39.143 213467 0.274 213467 −1.162 101757 0.164
2432908 −47.184 63077009 39.429 52645531 0.278 52645531 −1.157 113155 0.141
113155 −46.806 24815245 39.500 531919 0.279 531919 −1.153 213467 0.138
865043 −46.134 1843056 40.857 642659 0.279 642659 −1.133 gilteritinib 0.133
633318 −45.547 67199660 42.429 897789 0.280 897789 −1.125 215587 0.129
83885 −45.322 51234287 45.714 59052 0.281 215587 −1.114 1617909 0.127
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adeniyi, A.A.; Adeniyi, J.N.; Nlooto, M.; Singh, P. Probing New Antileukemia Agents That Target FLT3 and BCL-2 from Traditional Concoctions through a Combination of Mass Spectrometry Analysis and Consensus Docking Methods. Appl. Sci. 2022, 12, 11611. https://doi.org/10.3390/app122211611

AMA Style

Adeniyi AA, Adeniyi JN, Nlooto M, Singh P. Probing New Antileukemia Agents That Target FLT3 and BCL-2 from Traditional Concoctions through a Combination of Mass Spectrometry Analysis and Consensus Docking Methods. Applied Sciences. 2022; 12(22):11611. https://doi.org/10.3390/app122211611

Chicago/Turabian Style

Adeniyi, Adebayo A., Joy Nkechinyere Adeniyi, Manimbulu Nlooto, and Parvesh Singh. 2022. "Probing New Antileukemia Agents That Target FLT3 and BCL-2 from Traditional Concoctions through a Combination of Mass Spectrometry Analysis and Consensus Docking Methods" Applied Sciences 12, no. 22: 11611. https://doi.org/10.3390/app122211611

APA Style

Adeniyi, A. A., Adeniyi, J. N., Nlooto, M., & Singh, P. (2022). Probing New Antileukemia Agents That Target FLT3 and BCL-2 from Traditional Concoctions through a Combination of Mass Spectrometry Analysis and Consensus Docking Methods. Applied Sciences, 12(22), 11611. https://doi.org/10.3390/app122211611

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop