Interactions between cycloguanil derivatives and wild type and resistance-associated mutant Plasmodium falciparum dihydrofolate reductases

Application of molecular docking and PSO-SVR intelligent approaches in antimalarial activity prediction of enantiomeric cycloguanil analogues

A series of antifolate compounds, i.e. 1-(4-chlorophenyl)-6,6-dimethyl-1,3,5-triazine-2,4-diamine, or cycloguanil analogues, have shown effective inhibiting properties against Plasmodium falciparum dihydrofolate reductase (PfDHFR). In this work, the stereoselectivity of PfDHFR to the R and S enantiomer of cycloguanil analogues was obtained from molecular docking calculations and integrated into QSAR study to obtain a more accurate prediction model. Results indicate that PfDHFR can bind to cycloguanil analogues in the R and S enantiomers. Cycloguanil analogues with alkyl chain substituent prefer the R enantiomer over S because they do not experience steric hindrance with the Phe58 side chain, while cycloguanil analogues with phenol chain substituent prefer the S enantiomer over R because they do not experience steric hindrance with Leu46 and Met55 side chains. Particle swarm optimization and support vector regression were used to select relevant descriptors and generate the effective prediction model, with a high statistical significance level (r²_training = 0.941; r²_test = 0.884).

A novel prediction approach for antimalarial activities of Trimethoprim, Pyrimethamine, and Cycloguanil analogues using extremely randomized trees

Malaria is still one of the most serious diseases in tropical regions. This is due in part to the high resistance against available drugs for the inhibition of parasites, Plasmodium, the cause of the disease. New potent compounds with high clinical utility are urgently needed. In this work, we created a novel model using a regression tree to study structure-activity relationships and predict the inhibition constant, K_i of three different antimalarial analogues (Trimethoprim, Pyrimethamine, and Cycloguanil) based on their molecular descriptors. To the best of our knowledge, this work is the first attempt to study the structure-activity relationships of all three analogues combined. The most relevant descriptors and appropriate parameters of the regression tree are harvested using extremely randomized trees. These descriptors are water accessible surface area, Log of the aqueous solubility, total hydrophobic van der Waals surface area, and molecular refractivity. Out of all possible combinations of these selected parameters and descriptors, the tree with the strongest coefficient of determination is selected to be our prediction model. Predicted K_i values from the proposed model show a strong coefficient of determination, R²=0.996, to experimental K_i values. From the structure of the regression tree, compounds with high accessible surface area of all hydrophobic atoms (ASA_H) and low aqueous solubility of inhibitors (Log S) generally possess low K_i values. Our prediction model can also be utilized as a screening test for new antimalarial drug compounds which may reduce the time and expenses for new drug development. New compounds with high predicted K_i should be excluded from further drug development. It is also our inference that a threshold of ASA_H greater than 575.80 and Log S less than or equal to -4.36 is a sufficient condition for a new compound to possess a low K_i.

Ligand binding studies, preliminary structure-activity relationship and detailed mechanistic characterization of 1-phenyl-6,6-dimethyl-1,3,5-triazine-2,4-diamine derivatives as inhibitors of Escherichia coli dihydrofolate reductase

Gram-negative bacteria are implicated in the causation of life-threatening hospital-acquired infections. They acquire rapid resistance to multiple drugs and available antibiotics. Hence, there is the need to discover new antibacterial agents with novel scaffolds. For the first time, this study explores the 1,3,5-triazine-2,4-diamine and 1,2,4-triazine-2,4-diamine group of compounds as potential inhibitors of Escherichia coli DHFR, a pivotal enzyme in the thymidine and purine synthesis pathway. Using differential scanning fluorimetry, DSF, fifteen compounds with various substitutions on either the 3rd or 4th positions on the benzene group of 6,6-dimethyl-1-(benzene)-1,3,5-triazine-2,4-diamine were shown to bind to the enzyme with varying affinities. Then, the dose dependence of inhibition by these compounds was determined. Preliminary quantitative structure-activity relationship analysis and docking studies implicate the alkyl linker group and the sulfonyl fluoride group in increasing the potency of inhibition. 4-[4-[3-(4,6-diamino-2,2-dimethyl-1,3,5-triazin-1-yl)phenyl]butyl]benzenesulfonyl fluoride (NSC120927), the best hit from the study and a molecule with no reported inhibition of E. coli DHFR, potently inhibits the enzyme with a Ki value of 42.50 ± 5.34 nM, followed by 4-[6-[4-(4,6-diamino-2,2-dimethyl-1,3,5-triazin-1-yl)phenyl]hexyl]benzenesulfonyl fluoride (NSC132279), with a Ki value of 100.9 ± 12.7 nM. Detailed kinetic characterization of the inhibition brought about by five small-molecule hits shows that these inhibitors bind to the dihydrofolate binding site with preferential binding to the NADPH-bound binary form of the enzyme. Furthermore, in search of novel diaminotriazine scaffolds, it is shown that lamotrigine, a 1,2,4-triazine-3,5-diamine and a sodium-ion channel blocker class of antiepileptic drug, also inhibits E. coli DHFR. This is the first comprehensive study on the binding and inhibition brought about by diaminotriazines of a gram-negative prokaryotic enzyme and provides valuable insights into the SAR as an aid to the discovery of novel antibiotics.

A virtual screen discovers novel, fragment-sized inhibitors of Mycobacterium tuberculosis InhA

Isoniazid (INH) is usually administered to treat latent Mycobacterium tuberculosis (Mtb) infections and is used in combination therapy to treat active tuberculosis (TB). Unfortunately, resistance to this drug is hampering its clinical effectiveness. INH is a prodrug that must be activated by Mtb catalase-peroxidase (KatG) before it can inhibit InhA (Mtb enoyl-acyl-carrier-protein reductase). Isoniazid-resistant cases of TB found in clinical settings usually involve mutations in or deletion of katG, which abrogate INH activation. Compounds that inhibit InhA without requiring prior activation by KatG would not be affected by this resistance mechanism and hence would display continued potency against these drug-resistant isolates of Mtb. Virtual screening experiments versus InhA in the GO Fight Against Malaria (GO FAM) project were designed to discover new scaffolds that display base-stacking interactions with the NAD cofactor. GO FAM experiments included targets from other pathogens, including Mtb, when they had structural similarity to a malaria target. Eight of the 16 soluble compounds identified by docking against InhA plus visual inspection were modest inhibitors and did not require prior activation by KatG. The best two inhibitors discovered are both fragment-sized compounds and displayed Ki values of 54 and 59 μM, respectively. Importantly, the novel inhibitors discovered have low structural similarity to known InhA inhibitors and thus help expand the number of chemotypes on which future medicinal chemistry efforts can be focused. These new fragment hits could eventually help advance the fight against INH-resistant Mtb strains, which pose a significant global health threat.

Design of novel FLT-3 inhibitors based on dual-layer 3D-QSAR model and fragment-based compounds in silico

FMS-like tyrosine kinase 3 (FLT-3) is strongly correlated with acute myeloid leukemia, but no FLT-3-inhibitor cocomplex structure is available to assist the design of therapeutic inhibitors. Hence, we propose a dual-layer 3D-QSAR model for FLT-3 that integrates the pharmacophore, CoMFA, and CoMSIA. We then coupled the model with the fragment-based design strategy to identify novel FLT-3 inhibitors. In the first layer, the previously established model, Hypo02, was evaluated in terms of its correlation coefficient (r), RMS, cost difference, and configuration cost, with values of 0.930, 1.24, 106.45, and 16.44, respectively. Moreover, Fischer's cross-validation test of data generated by Hypo02 yielded a 98% confidence level, and the validation of the testing set yielded a best r value of 0.87. The features of Hypo02 were separated into two parts and then used to screen the MiniMaybridge fragment compound database. Nine novel FLT-3 inhibitors were generated in this layer. In the second layer, Hypo02 was subjected to an alignment rule to generate CoMFA- and CoMSIA-based models, for which the partial least-squares validation method was utilized. The values of q(2), r(2), and predictive r(2) were 0.58, 0.98, and 0.76, respectively, derived from the CoMFA model with steric and electrostatic fields. The CoMSIA model with five different fields yielded values of 0.54, 0.97, and 0.76 for q(2), r(2), and predictive r(2), respectively. The CoMFA and CoMSIA models were used to constrain 3D structures of the nine novel FLT-3 inhibitors. This dual-layer 3D-QSAR model constitutes a valuable tool to easily and quickly screen and optimize novel potential FLT-3 inhibitors for the treatment of acute myeloid leukemia.

Binding modes of 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs to P. falciparum dihydrofolate reductase enzyme: Molecular docking studies

A molecular docking study was carried out on 28 compounds belonging to 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs using Glide, FlexX and GOLD programs and the X-ray crystallographic structures of the quadruple mutant (1J3K:pdb) and wild type (1J3I:pdb) Plasmodium falciparum dihydrofolate reductase enzyme. The experimental conformation the bound ligand WR99210 was precisely reproduced by the docking procedures as demonstrated by low (<2.00 Å) root-mean-square deviations. The results indicated that most of the compounds dock into the active sites of both the wild type and quadruple mutant P. falciparum dihydrofolate reductase enzymes. Visual inspection of the binding modes also demonstrated that most of the compounds could form H-bond interactions with the key amino acid residues (Asp54, Ile14 and Leu/Ile164) and with better docking scores than the bound compound (5). Their long side chains orient in the hydrophobic portion of the active site which is occupied by trichloro aryloxy side chain of WR99210 (5). Thus, avoid potential steric clashes with Asn108 (mutated from Ser108). Such a clash is known to be responsible for the resistance of the P. falciparum to pyrimethamine and cycloguanil.