Three dimensional pharmacophore modelling for c-Kit receptor tyrosine kinase inhibitors

Practical Synthesis of 7-Bromo-4-chloro-1H-indazol-3-amine: An Important Intermediate to Lenacapavir

7-Bromo-4-chloro-1H-indazol-3-amine is a heterocyclic fragment used in the synthesis of Lenacapavir, a potent capsid inhibitor for the treatment of HIV-1 infections. In this manuscript, we describe a new approach to synthesizing 7-bromo-4-chloro-1H-indazol-3-amine from inexpensive 2,6-dichlorobenzonitrile. This synthetic method utilizes a two-step sequence including regioselective bromination and heterocycle formation with hydrazine to give the desired product in an overall isolated yield of 38-45%. The new protocol has been successfully demonstrated on hundred-gram scales without the need for column chromatography purification. This new synthesis provides a potential economical route to the large-scale production of this heterocyclic fragment of Lenacapavir.

Carica papaya Reduces Muscle Insulin Resistance via IR/GLUT4 Mediated Signaling Mechanisms in High Fat Diet and Streptozotocin-Induced Type-2 Diabetic Rats

In the management of type 2 diabetes, oral antidiabetic drugs have several side effects, which in turn have led the pharmaceutical industry to search for good therapeutic, non-toxic and reliable drugs. Carica papaya (C. papaya) is one of several plants in nature that have been found to possess anti-diabetic properties. Despite studies being focused on the antidiabetic activity of C. papaya, the molecular mechanism against high fat diet induced insulin resistance is yet to be identified. The role of C. papaya was evaluated on insulin signaling molecules, such as the insulin receptor (IR) and glucose transporter-4 (GLUT4) in high fat, diet-streptozotocin induced type 2 diabetic rats, and analyzed the bioactive compounds of C. papaya against IR and GLUT4 via molecular docking and dynamics. The ethanolic extract of C. papaya leaves (600 mg/kg of body weight) was given daily to male wistar rats for 45 days and we observed the various biochemical parameters, gene expression analysis and histopathology of skeletal muscle. Molecular docking and dynamics were undertaken to understand the bioactive compounds with the greatest hit rate. C. papaya treatment was able to control blood glucose levels, the lipid profile and serum insulin, but it facilitated tissue antioxidant enzymes and IR and GLUT4 levels. The in-silico study showed that kaempferol, quercitin and transferulic acid were the top three ligands with the greatest hit rate against the protein targets. Our preliminary findings, for the first time, showed that C. papaya reinstates the glycemic effect in the diabetic skeletal muscle by accelerating the expression of IR and GLUT4.

Discovery of novel potential KIT inhibitors for the treatment of gastrointestinal stromal tumor

Numerous inhibitors of tyrosine-protein kinase KIT, a receptor tyrosine kinase, have been explored as a viable therapy for the treatment of gastrointestinal stromal tumor (GIST). However, drug resistance due to acquired mutations in KIT makes these drugs almost useless. The present study was designed to screen the novel inhibitors against the activity of the KIT mutants through pharmacophore modeling and molecular docking. The best two pharmacophore models were established using the KIT mutants’ crystal complexes and were used to screen the new compounds with possible KIT inhibitory activity against both activation loop and ATP-binding mutants. As a result, two compounds were identified as potential candidates from the virtual screening, which satisfied the potential binding capabilities, molecular modeling characteristics, and predicted absorption, distribution, metabolism, excretion, toxicity (ADMET) properties. Further molecular docking simulations showed that two compounds made strong hydrogen bond interaction with different KIT mutant proteins. Our results indicated that pharmacophore models based on the receptor–ligand complex had excellent ability to screen KIT inhibitors, and two compounds may have the potential to develop further as the future KIT inhibitors for GIST treatment.

Beta-sitosterol attenuates insulin resistance in adipose tissue via IRS-1/Akt mediated insulin signaling in high fat diet and sucrose induced type-2 diabetic rats

In our previous study, we have shown that β-sitosterol (SIT) enhances glycemic control by increasing the activation of insulin receptor (IR) and glucose transporter 4 (GLUT4) proteins in adipose tissue. However, the possible role of SIT on the regulation of post-receptor insulin signal transduction is not known. Hence, the study was aimed to assess the effects of SIT on IRS-1/Akt mediated insulin signaling molecules in high-fat diet and sucrose induced type-2 diabetic rats. An oral effective dose of SIT (20 mg/kg b.wt) was given for 30 days to high fat-fed type-2 diabetic rats to find out whether SIT regulates IRS-1/Akt pathway of insulin signaling. The results showed that SIT attenuated the insulin receptor substrate-1 serine phosphorylation (p-IRS-1^Ser636) (P = 0.0003). However, it up-regulated the mRNA expression of IR (P = 0.0036) and post-receptor insulin signaling molecules such as IRS-1 (P < 0.0001), β-arrestin-2 (P < 0.0058), Akt (P = 0.0008), AS160 (P = 0.0030) and GLUT4 (P < 0.0001) with a concomitant increase in the levels of IRS-1(P < 0.0001), p-IRS1-1^Tyr632 (P = 0.0014), Akt (P < 0.0001), p-Akt^{Ser473/Thr308} (P = 0.0006; P < 0.0001), AS160 and p-AS160^Thr642 (P < 0.0001) compared with type-2 diabetic rats. In Silico analysis was also performed and it showed that SIT possesses the greater binding affinity with β-arrestin-2, c-Src, and IRS-1 as well as Akt proteins and proved to attenuate insulin resistance as this study coincides with in vivo findings. Our present study clearly shows that SIT attenuates high fat diet-induced detrimental changes in adipose tissue. Therefore, it is concluded from the present findings that, SIT could be used as potential therapeutic phytomedicine for the management of type-2 diabetes.

Identifying dual leucine zipper kinase (DLK) inhibitors using e-pharamacophore screening and molecular docking

Alzheimer's is a neural disorder causing gradual loss in structure and function of nerve cell. To treat such disorders, c-Jun N-terminal Kinase (JNK) Pathway inhibitors were developed by representing chemical compounds that were used to inhibit the JNK signaling pathways. DLK is the stress sensor and implicating as regulatory factor in JNK pathway. Therefore, in the present investigation, pharmacophore screening was tried to identify the chemical compounds that involving inhibition of DLK proteins. To explore the pharmacophore region and mode of binding with DLK protein, N- (I H-pyrazol-3-y l) pyridin-2-aminer inhibitors were docked with DLK. Results reveal the information on the interaction mechanism of protein and ligand with chemical characteristics required to inhibit DLK protein. Such predicted information (AAAARH) was used as query to find out potential novel lead compounds sourced from public database. As an outcome of 65 compounds were listed based on the fitness score (2≥), and were subjected to glide HTVS.SP and XP. Best performing 5 lead compounds were shortlisted for dynamic simulations. This exhibited a constant RMSD over 20 ns of timescale.

Discovery of Dual Inhibitors for Wild Type and D816V Mutant of c-KIT Kinase through Virtual and Biochemical Screening of Natural Products

Although stem cell factor receptor (c-KIT) kinase is responsible for various malignant human cancers, the presence of constitutively active gain-of-function mutants has made it difficult to discover new anticancer agents using c-KIT as the target protein. To identify the common inhibitors of wild-type c-KIT and the most abundant gain-of-function mutant (D816V), the virtual screening of natural products was performed for the two target proteins in parallel with the scoring function improved by implementing a sophisticated solvation free energy term. As a result, four common inhibitors of natural origin are found with biochemical potencies ranging from low micromolar to submicromolar levels. The results of extensive docking simulations show that although the natural-product inhibitors establish weaker hydrophobic interactions with the D816V mutant than with the wild type, they exhibit a little higher inhibitory activity for the former than the latter by strengthening the hydrogen-bond interactions to a sufficient extent. Of the four natural-product inhibitors, (Z)-6-hydroxy-2-(4-methoxybenzylidene)benzofuran-3(2H)-one (3) is anticipated to serve as a new molecular core for the structure-activity relationship studies to optimize the biochemical potencies because it exhibits good inhibitory activity against both the wild type and D816V mutant despite its low molecular weight (268.3 amu).

Structure-based de novo design and identification of D816V mutant-selective c-KIT inhibitors

To identify potent and selective inhibitors of D816V, the most common gain-of-function c-KIT mutant, we carried out structure-based de novo design using 7-azaindole as the core and the scoring function improved by implementing an accurate solvation free energy term. This approach led to the identification of new c-KIT inhibitors specific for the D816V mutant. The 3-(3,4-dimethoxyphenyl)-7-azaindole scaffold was optimized and represents a lead structure for the design of the potent and specific inhibitors of the D816V mutant. The results of molecular dynamics simulations indicate that hydrogen bonding interactions between the 7-azadindole moiety and the backbone groups of Cys673 are the most significant determinant for the potency and selectivity of c-KIT inhibitors.

High throughput virtual screening and E-pharmacophore filtering in the discovery of new BACE-1 inhibitors

Alzheimer's disease is a progressive neurodegenerative disorder, which is characterized by amyloid β peptide deposition in the brain. Aβ peptide, the major component of amyloid plaques is generated by the sequential processing of a larger protein called amyloid Precursor Protein by β-amyloid cleaving enzyme (BACE-1). In this study, we appllied computer assisted methodology unifying molecular docking and pharmacophore filtering to identify potent inhibitors against BACE-1. In order to inspect the pharmacophore region and binding mode of BACE-1 135 reported co-crystallized ligands of BACE-1 were docked into the active site using Glide XP. The present molecular docking studies provided critical information on protein ligand interactions that revealed imminent information on chemical features essential to inhibiting BACE-1. Based on the docking results we proposed structure based pharmacophore features that hold well as potent BACE-1 inhibitors. A huge set of compounds was docked into the active site of BACE-1 and the hits from the docking were filtered to match the chemical features of the pharmacophore model. The compounds resulting from the pharmacophore filtering were again re-docked into the active site of BACE-1 and the three hits bound well into the active sites and matched the pharmacophore models which were identified as possible potential inhibitors of BACE-1. Molecular dynamics simulation reveals that lead 3 shows constant RMSD and the number of hydrogen bonding with the protein among the identified three lead molecules.

Pharmacophore-based virtual screening and density functional theory approach to identifying novel butyrylcholinesterase inhibitors

AIM

To identify the critical chemical features, with reliable geometric constraints, that contributes to the inhibition of butyrylcholinesterase (BChE) function.

METHODS

Ligand-based pharmacophore modeling was used to identify the critical chemical features of BChE inhibitors. The generated pharmacophore model was validated using various techniques, such as Fischer's randomization method, test set, and decoy set. The best pharmacophore model was used as a query in virtual screening to identify novel scaffolds that inhibit BChE. Compounds selected by the best hypothesis in the virtual screening were tested for drug-like properties, and molecular docking study was applied to determine the optimal orientation of the hit compounds in the BChE active site. To find the reactivity of the hit compounds, frontier orbital analysis was carried out using density functional theory.

RESULTS

Based on its correlation coefficient (0.96), root mean square (RMS) deviation (1.01), and total cost (105.72), the quantitative hypothesis Hypo1 consisting of 2 HBA, 1 Hy-Ali, and 1 Hy-Ar was selected as the best hypothesis. Thus, Hypo1 was used as a 3D query in virtual screening of the Maybridge and Chembridge databases. The hit compounds were filtered using ADMET, Lipinski's Rule of Five, and molecular docking to reduce the number of false positive results. Finally, 33 compounds were selected based on their critical interactions with the significant amino acids in BChE's active site. To confirm the inhibitors' potencies, the orbital energies, such as HOMO and LUMO, of the hit compounds and 7 training set compounds were calculated. Among the 33 hit compounds, 10 compounds with the highest HOMO values were selected, and this set was further culled to 5 compounds based on their energy gaps important for stability and energy transfer. From the overall results, 5 hit compounds were confirmed to be potential BChE inhibitors that satisfied all the pharmacophoric features in Hypo1.

CONCLUSION

This study pinpoints important chemical features with geometric constraints that contribute to the inhibition of BChE activity. Five compounds are selected as the best hit BchE-inhibitory compounds.

Development, evaluation and application of 3D QSAR Pharmacophore model in the discovery of potential human renin inhibitors

Background

Renin has become an attractive target in controlling hypertension because of the high specificity towards its only substrate, angiotensinogen. The conversion of angiotensinogen to angiotensin I is the first and rate-limiting step of renin-angiotensin system and thus designing inhibitors to block this step is focused in this study.

Methods

Ligand-based quantitative pharmacophore modeling methodology was used in identifying the important molecular chemical features present in the set of already known active compounds and the missing features from the set of inactive compounds. A training set containing 18 compounds including active and inactive compounds with a substantial degree of diversity was used in developing the pharmacophore models. A test set containing 93 compounds, Fischer randomization, and leave-one-out methods were used in the validation of the pharmacophore model. Database screening was performed using the best pharmacophore model as a 3D structural query. Molecular docking and density functional theory calculations were used to select the hit compounds with strong molecular interactions and favorable electronic features.

Results

The best quantitative pharmacophore model selected was made of one hydrophobic, one hydrogen bond donor, and two hydrogen bond acceptor features with high a correlation value of 0.944. Upon validation using an external test set of 93 compounds, Fischer randomization, and leave-one-out methods, this model was used in database screening to identify chemical compounds containing the identified pharmacophoric features. Molecular docking and density functional theory studies have confirmed that the identified hits possess the essential binding characteristics and electronic properties of potent inhibitors.

Conclusion

A quantitative pharmacophore model of predictive ability was developed with essential molecular features of a potent renin inhibitor. Using this pharmacophore model, two potential inhibitory leads were identified to be used in designing novel and future renin inhibitors as antihypertensive drugs.

Potent BACE-1 inhibitor design using pharmacophore modeling, in silico screening and molecular docking studies

BACKGROUND

Beta-site amyloid precursor protein cleaving enzyme (BACE-1) is a single-membrane protein belongs to the aspartyl protease class of catabolic enzymes. This enzyme involved in the processing of the amyloid precursor protein (APP). The cleavage of APP by BACE-1 is the rate-limiting step in the amyloid cascade leading to the production of two peptide fragments Aβ40 and Aβ42. Among two peptide fragments Aβ42 is the primary species thought to be responsible for the neurotoxicity and amyloid plaque formation that lead to memory and cognitive defects in Alzheimer's disease (AD). AD is a ravaging neurodegenerative disorder for which no disease-modifying treatment is currently available. Inhibition of BACE-1 is expected to stop amyloid plaque formation and emerged as an interesting and attractive therapeutic target for AD.

METHODS

Ligand-based computational approach was used to identify the molecular chemical features required for the inhibition of BACE-1 enzyme. A training set of 20 compounds with known experimental activity was used to generate pharmacophore hypotheses using 3D QSAR Pharmacophore Generation module available in Discovery studio. The hypothesis was validated by four different methods and the best hypothesis was utilized in database screening of four chemical databases like Maybridge, Chembridge, NCI and Asinex. The retrieved hit compounds were subjected to molecular docking study using GOLD 4.1 program.

RESULTS

Among ten generated pharmacophore hypotheses, Hypo 1 was chosen as best pharmacophore hypothesis. Hypo 1 consists of one hydrogen bond donor, one positive ionizable, one ring aromatic and two hydrophobic features with high correlation coefficient of 0.977, highest cost difference of 121.98 bits and lowest RMSD value of 0.804. Hypo 1 was validated using Fischer randomization method, test set with a correlation coefficient of 0.917, leave-one-out method and decoy set with a goodness of hit score of 0.76. The validated Hypo 1 was used as a 3D query in database screening and retrieved 773 compounds with the estimated activity value <100 nM. These hits were docked into the active site of BACE-1 and further refined based on molecular interactions with the essential amino acids and good GOLD fitness score.

CONCLUSION

The best pharmacophore hypothesis, Hypo 1, with high predictive ability contains chemical features required for the effective inhibition of BACE-1. Using Hypo 1, we have identified two compounds with diverse chemical scaffolds as potential virtual leads which, as such or upon further optimization, can be used in the designing of new BACE-1 inhibitors.