QM/MM and molecular dynamics investigation of the mechanism of covalent inhibition of TAK1 kinase.

TAK1 is a serine/threonine kinase which is involved in the moderation of cell survival and death via the TNFα signalling pathway. It is also implicated in a range of cancer and anti-inflammatory diseases. Drug discovery efforts on this target have focused on both traditional reversible ATP-binding site inhibitors and increasingly popular irreversible covalent binding inhibitors. Irreversible inhibitors can offer benefits in terms of potency, selectivity and PK/PD meaning they are increasingly pursued where the strategy exists. TAK1 kinase differs from the better-known kinase EGFR in that the reactive cysteine nucleophile targeted by electrophilic inhibitors is located towards the back of the ATP binding site, not at its mouth. While a wealth of structural and computational effort has been spent exploring EGFR, only limited studies on TAK1 have been reported. In this work we report the first QM/MM study on TAK1 aiming to better understand aspects of covalent adduct formation. Our goal is to identify the general base in the catalytic reaction, whether the process proceeds via a stepwise or concerted pathway, and how the highly flexible G-loop and A-loop affect the catalytic cysteine located nearby.

Preparation, biological & cheminformatics-based assessment of N2,N4-diphenylpyrimidine-2,4-diamine as potential Kinase-targeted antimalarials.

Twenty eight new N2,N4-diphenylpyrimidine-2,4-diamines have been prepared in order to expand our understanding of the anti-malarial SAR of the scaffold. The aim of the study was to make structural modifications to improve the overall potency, selectivity and solubility of the series by varying the anilino groups attached to the 2- and 4-position. We evaluated the activity of the compounds against Plasmodium falciparum (Pf) 3D7, cytotoxicity against HepG2, % inhibition at a panel of 10 human kinases, solubility, permeability and lipophilicity, and human and rat in vitro clearance. 11 was identified as a potent anti-malarial with an IC50 of 0.66 µM at the 3D7 strain and a selectivity (SI) of ~ 40 in terms of cytotoxicity against the HepG2 cell line. It also displayed low experimental logD7.4 (2.27), reasonable solubility (124 µg/ml), good metabolic stability, but low permeability. A proteo-chemometric workflow was employed to identify putative Pf targets of the most promising compounds. Ligand-based similarity searching of the ChEMBL database led to the identification of most probable human targets. These were then used as input for sequence-based searching of the Pf proteome. Homology modelling and molecular docking were used to evaluate whether compounds could indeed bind to these targets with valid binding modes. In vitro biological testing against close human analogs of these targets was subsequently undertaken. This allowed us to identify potential Pf targets and human anti-targets that could be exploited in future development.

Theoretical Investigation of the Enantioselective [4 + 2] Cycloaddition Reaction of o-Hydroxystyrene and Azlactone.

Theoretical studies have been undertaken to rationalize the origin of the enantioselective Diels-Alder reaction (DA) of o-hydroxystyrene and azlactone catalyzed by (a) chiral BINOL-phosphoric acid (CPA) and (b) CPA and chiral guanidine (TBO). The sequence of events leading to increased enantioselectivity under the latter conditions have been studied using density functional theory (DFT) methods. The computational results indicate that both the mono- and co-catalytic processes proceed via stepwise [4 + 2] cycloaddition reactions involving three steps, which are (1) C-C bond formation, (2) C-O bond formation, and (3) the opening of the azlactone ring. This results in the formation of an oxygenous cycle with one chiral center. The origin of greater enantioselectivity under the latter catalytic conditions are discussed in terms of the structural characteristics and energetics of the intermediates and transition states formed on the potential energy surface of the competing reactions.

Design, synthesis and evaluation of N2,N4-diaminoquinazoline based inhibitors of phosphodiesterase type 5.

We describe the design, synthesis and evaluation of a series of N2,N4-diaminoquinazoline analogs as PDE5 inhibitors. Twenty compounds were prepared and these were assessed in terms of their PDE5 and PDE6 activity, ex-vivo vasodilation response, mammalian cytotoxicity and aqueous solubility. Molecular docking was used to determine the binding mode of the series and this was demonstrated to be consistent with the observed SAR. Compound 15 was the most active PDE5 inhibitor (IC50 = 0.072 ±â€¯0.008 µM) and exhibited 4.6-fold selectivity over PDE6. Ex-vivo assessment of 15 and 22 in a rat pulmonary artery vasodilation model demonstrated EC50s of 1.63 ±â€¯0.72 µM and 2.28 ±â€¯0.74 µM respectively.

Insights into the EGFR SAR of N-phenylquinazolin-4-amine-derivatives using quantum mechanical pairwise-interaction energies.

Protein kinases are an important class of enzymes that play an essential role in virtually all major disease areas. In addition, they account for approximately 50% of the current targets pursued in drug discovery research. In this work, we explore the generation of structure-based quantum mechanical (QM) quantitative structure-activity relationship models (QSAR) as a means to facilitate structure-guided optimization of protein kinase inhibitors. We explore whether more accurate, interpretable QSAR models can be generated for a series of 76 N-phenylquinazolin-4-amine inhibitors of epidermal growth factor receptor (EGFR) kinase by comparing and contrasting them to other standard QSAR methodologies. The QM-based method involved molecular docking of inhibitors followed by their QM optimization within a ~ 300 atom cluster model of the EGFR active site at the M062X/6-31G(d,p) level. Pairwise computations of the interaction energies with each active site residue were performed. QSAR models were generated by splitting the datasets 75:25 into a training and test set followed by modelling using partial least squares (PLS). Additional QSAR models were generated using alignment dependent CoMFA and CoMSIA methods as well as alignment independent physicochemical, e-state indices and fingerprint descriptors. The structure-based QM-QSAR model displayed good performance on the training and test sets (r2 ~ 0.7) and was demonstrably more predictive than the QSAR models built using other methods. The descriptor coefficients from the QM-QSAR models allowed for a detailed rationalization of the active site SAR, which has implications for subsequent design iterations.

Elucidation of Vasodilation Response and Structure Activity Relationships of N²,N4-Disubstituted Quinazoline 2,4-Diamines in a Rat Pulmonary Artery Model.

Pulmonary arterial hypertension (PAH) is a rare and progressive disease arising from various etiologies and pathogenesis. PAH decreases life expectancy due to pulmonary vascular remodeling, elevation of mean pulmonary arterial pressure, and ultimately progresses to heart failure. While clinical treatments are available to reduce the associated symptoms, a complete cure has yet to be found. Phosphodiesterase-5 (PDE-5) inhibition has been identified as a possible intervention point in PAH treatment. The functional vasodilation response to N²,N4-diamino quinazoline analogues with differing PDE-5 inhibitory activities and varying physicochemical properties were assessed in both endothelium-intact and denuded rat pulmonary arteries to gain greater insight into their mode of action. All analogues produced vasorelaxant effects with EC50s ranging from 0.58 ± 0.22 µM to ˃30 µM. It was observed that vasodilation response in intact vessels was highly correlated with that of denuded vessels. The ~10% drop in activity is consistent with a loss of the nitric oxide mediated cyclic guanosine monophosphate (NO/cGMP) pathway in the latter case. A moderate correlation between the vasodilation response and PDE-5 inhibitory activity in the intact vessels was observed. Experimental protocol using the alpha-adrenergic (α1) receptor agonist, phenylephrine (PE), was undertaken to assess whether quinazoline derivatives showed competitive behavior similar to the α1 receptor blocker, prazosin, itself a quinazoline derivative, or to the PDE-5 inhibitor, sildenafil. Competitive experiments with the α1-adrenergic receptor agonist point to quinazoline derivatives under investigation here act via PDE-5 inhibition and not the former. The pre-incubation of pulmonary arterial rings with quinazoline test compounds (10 µM) reduced the contractile response to PE around 40⁻60%. The most promising compound (9) possessed ~32 folds higher selectivity in terms of vasodilation to its mammalian A549 cell cytotoxicity. This study provides experi0 0mental basis for PDE-5 inhibition as the mode of action for vasodilation by N²,N4-diamino quinazoline analogues along with their safety studies that may be beneficial in the treatment of various cardiovascular pathologies.

QM/MM Investigation for Protonation States in a Bilin Reductase PcyA-Biliverdin IXα Complex.

Herein we report quantum mechanical/molecular mechanical (QM/MM) studies to investigate the most probable protonation states of active site amino acids and bound substrate based on a recently reported neutron diffraction structure of phycocyanobilin:ferredoxin oxidoreductase (PcyA) by Unno et al. This structure was considered to be bound in its initial state of biliverdin IXα (BV), which has the C-pyrrole ring in the deprotonated state. The protonation state of BV suggested by neutron and spectroscopic studies is a stable, two-electron reduced complex with a bound hydronium ion. Several ambiguities in the neutron structure were observed which prompted a further theoretical analysis of the structure. This structural investigation provides new understanding of the PcyA and BV protonation states not previously reported in the literature. Our calculations suggest that the hydronium ion (H3 O+ ) is energetically unfavorable, preferentially protonating the neighboring His88 residue and that the C-ring of BV is not protonated.

Elucidation of the catalytic mechanism of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase using QM/MM calculations.

The folate pathway is a recognized intervention point for treating parasitic and bacterial infections in humans. However, the efficacy of treatments targeting dihydropteroate synthase (DHPS) and dihydrofolate reductase (DHFR) has reduced due to disease-related mutations. This has prompted interest in other enzyme targets on this clinically validated pathway, including 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK). A challenge in the design of molecules to target this enzyme is that the precise mechanism of the reaction and the role of the active site residues are not fully understood. In this study, we report the first theoretical analysis of the catalytic pathway of the natural substrate using hybrid quantum mechanical/molecular mechanical (QM/MM) methods. The reaction profiles associated with three proposed general bases have been investigated, as well as the profile for two mutant enzymes, namely R92A and R82A. We identified R92 as the general base in the wildtype reaction. The predicted barriers are in good agreement with the observed experimental kcat values obtained for wildtype and mutant proteins.

Synthesis and evaluation of the NSCLC anti-cancer activity and physical properties of 4-aryl-N-phenylpyrimidin-2-amines.

Reported herein are efforts to profile 4-aryl-N-phenylpyrimidin-2-amines in terms of their anti-cancer activity towards non small-cell lung carcinoma (NSCLC) cells. We have synthesized new 4-aryl-N-phenylpyrimidin-2-amines and assessed them in terms of their cytotoxicity (A549, NCI-H187, MCF7, Vero & KB) and physicochemical properties (logD7.4 and solubility). 13f and 13c demonstrated potent anti-cancer activity in A549 cells (0.2µM), compared to 0.4µM for the NSCLC drug Doxorubicin. 13f also displayed low experimental logD7.4 (2.9) and the best solubility (∼40µM). Compounds 13b and 13d showed the best balance of A549 anti-cancer activity and selectivity. 13g showed good activity and selectivity comparable with the anti-cancer drug Doxorubicin.

Exploring the catalytic mechanism of dihydropteroate synthase: elucidating the differences between the substrate and inhibitor.

Dihydropteroate synthase (DHPS) catalyzes the condensation of 6-hydroxymethyl-7,8-dihydropterin pyrophosphate (DHPPP) with p-aminobenzoic acid (pABA) and is a well validated target for anti-malarial and anti-bacterial drugs. However, in recent years its utility as a therapeutic target has diminished considerably due to multiple mutations. As such, considerable structural biology and medicinal chemistry effort has been expended to understand and overcome this issue. To date no detailed computational analysis of the protein mechanism has been made despite the detailed crystal structures and multiple mechanistic proposals being made. In this study the mechanistic proposals for DHPS have been systematically investigated using a hybrid QM/MM method. We aimed to compare the energetics associated with SN1 and SN2 processes, whether the SN1 process involves a carbocation or neutral DHP intermediate, uncover the identity of the general base in the catalytic mechanism, and understand the differences in substrate vs. inhibitor reactivity. Our results suggest a reaction that follows an SN1 process with the rate determining step being C-O bond breaking to give a carbocation intermediate. Comparative studies on the inhibitor STZ confirm the experimental observations that it is also a DHPS substrate.

Skin sensitization prediction using quantum chemical calculations: a theoretical model for the S(N)Ar domain.

It is widely accepted that skin sensitization begins with the sensitizer in question forming a covalent adduct with a protein electrophile or nucleophile. We investigate the use of quantum chemical methods in an attempt to rationalize the sensitization potential of chemicals of the S(N)Ar reaction domain. We calculate the full reaction profile for 23 chemicals with experimental sensitization data. For all quantitative measurements, we find that there is a good correlation between the reported pEC3 and the calculated barrier to formation of the low energy product or intermediate (r(2) = 0.64, N = 12) and a stronger one when broken down by specific subtype (r(2) > 0.9). Using a barrier cutoff of ∼10 kcal/mol allows us to categorize 100% (N = 12) of the sensitizers from the nonsensitizers (N = 11), with just 1 nonsensitizer being mispredicted as a weak sensitizer (9%). This model has an accuracy of ∼96%, with a sensitivity of 100% and a specificity of ∼91%. We find that the kinetic and thermodynamic information provided by the complete profile can help in the rationalization process, giving additional insight into a chemical's potential for skin sensitization.

Comparison of feline and human immunodeficiency virus reverse transcriptase enzymes through chemical screening and computational analysis.

Feline immunodeficiency virus (FIV) is a common infection found in domesticated and wild cats worldwide. Despite the wealth of therapeutic understanding of the disease in humans, considerably less information exists regarding the treatment of the disease in felines. Current treatment relies on drugs developed for the related human immunodeficiency virus (HIV) and includes compounds of the popular non-nucleotide reverse transcriptase (NNRTI) class. This is despite FIV-RT being only 67% similar to HIV-1 RT at the enzyme level, increasing to 88% for the allosteric pocket targeted by NNRTIs. The goal of this project was to try to quantify how well the more extensive pharmacological knowledge available for human disease translates to felines. To this end we screened known NNRTIs and 10 diverse pyrimidine analogs identified virtually. We use this chemo-centric probe approach to (a) assess the similarity between the two related RT targets based on the observed experimental inhibition values, (b) try to identify more potent inhibitors at FIV, and (c) gain a better appreciation of the structure-activity relationships (SAR). We found the correlation between IC50s at the two targets to be strong (r2 = 0.87) and identified compound 1 as the most potent inhibitor of FIV with IC50 of 0.030 µM ± 0.009. This compared to FIV IC50 values of 0.22 ± 0.17 µM, 0.040 ± 0.010 µM and >160 µM for known anti HIV-1 RT drugs Efavirenz, Rilpivirine, and Nevirapine, respectively. This knowledge, along with an understanding of the structural origin that give rise to any differences could improve the way HIV drugs are repurposed for FIV.

Synthesis of substituted 2-arylindanes from E-(2-stilbenyl)methanols via Lewis acid-mediated cyclization and nucleophililc transfer from trialkylsilyl reagents.

A preparative method for the synthesis of functionalized 2-arylindanes has been developed via the Lewis acid-mediated ring closure of stilbenyl methanols followed by nucleophilic transfer from trialkylsilyl reagents. The reactions gave the corresponding products in moderate to high yields and diastereoselectivity. The solvent as well as the nucleophile played an important role in determining the type(s) of product arising either from nucleophilic addition (indanes) or loss of a proton ß to the indanyl-type carbocations (indenes). Electron-donating groups on the fused aromatic ring (Y and Z = OMe) or the presence of electron-withdrawing groups (NO2) on the nonfused Ar ring facilitate the cyclization. In contrast, the presence of electron-donating groups (OMe) on the nonfused Ar ring impedes the process. In the case of Cl on the nonfused Ar ring, temperature modulates the resonance versus inductive field effects on the overall reaction pathways involving cyclization to form the indanyl-type cation. Quantum chemical calculations supported the intermediacy of the carbocation species and the transfer of hydride from triethylsilane (Nu = H) to the indanyl-type cations to form the trans-1,2-disubstituted indane as the single diastereomer product.

The new quinazoline derivative (N2-methyl-N4-[(thiophen-2-yl)methyl]quinazoline-2,4-diamine) vasodilates isolated mesenteric arteries through endothelium-independent mechanisms and has acute hypotensive effects in Wistar rats.

During the screening of new N2,N4-disubstituted quinazoline 2,4-diamines as phosphodiesterase-5 inhibitors and pulmonary artery vasodilators, one N2-methyl-N4-[(thiophen-2-yl)methyl]quinazoline-2,4-diamine (compound 8) presented a greater selectivity for systemic than pulmonary vasculature. The present study aimed to characterize its vasorelaxant and hypotensive effects in Wistar rats. Vasorelaxant effects of compound 8 and underlying mechanisms were evaluated on isolated mesenteric arteries. Acute hypotensive effect was evaluated in anesthetized rats. Additionally, cell viability and cytochrome P450 (CYP) activities were studied in rat isolated hepatocytes. Nifedipine was used as a comparator. Compound 8 induced a strong vasorelaxant effect, similar to nifedipine. This was unaffected by endothelium removal but was decreased by inhibitors of guanylate cyclase (ODQ) and KCa channel (iberiotoxin). Compound 8 enhanced sodium nitroprusside-induced relaxation, but inhibited vasoconstriction evoked by α1-adrenergic receptor activation and extracellular Ca2+ influx via receptor-operated Ca2+ channels. Acute intravenous infusion of compound 8 (0.05 and 0.1 mg/kg) produced hypotension. It showed similar potency to nifedipine for lowering diastolic and mean arterial blood pressure, but less so for the effect on systolic blood pressure. Compound 8 had no effect on hepatocyte viability and CYP activities except at high concentration (10 µM) at which a weak inhibitory effect on CYP1A and 3A was observed. In conclusion, this study identified a N2-methyl-N4-[(thiophen-2-yl)methyl]quinazoline-2,4-diamine with a potent vasodilator effect on resistance vessels, leading to an acute hypotensive effect and a low risk of liver toxicity or drug-drug interactions. These vascular effects were mediated mainly through sGC/cGMP pathway, opening of KCa channels, and inhibition of calcium entry.

Evaluating the enthalpic contribution to ligand binding using QM calculations: effect of methodology on geometries and interaction energies.

As a result of research on ligand efficiency in the pharmaceutical industry, there is greater focus on optimizing the strength of polar interactions within receptors, so that the contribution of overall size and lipophilicity to binding can be decreased. A number of quantum mechanical (QM) methods involving simple probes are available to assess the H-bonding potential of different heterocycles or functional groups. However, in most receptors, multiple features are present, and these have distinct directionality, meaning very minimalist models may not be so ideal to describe the interactions. We describe how the use of gas phase QM models of kinase protein-ligand complex, which can more closely mimic the polar features of the active site region, can prove useful in assessing alterations to a core template, or different substituents. We investigate some practical issues surrounding the use of QM cluster models in structure based design (SBD). These include the choice of the method; semi-empirical, density functional theory or ab-initio, the choice of the basis set, whether to include implicit or explicit solvation, whether BSSE should be included, etc. We find a combination of the M06-2X method and the 6-31G* basis set is sufficiently rapid, and accurate, for the computation of structural and energetic parameters for this system.

Preparation, biological evaluation and QSAR analysis of urea substituted 2,4-diamino-pyrimidine anti-malarials.

The synthesis and evaluation of twenty six new phenylurea substituted 2,4-diamino-pyrimidines against Plasmodium falciparum (Pf) 3D7 are reported. Compounds were prepared to improve both anti-malarial activity and selectivity of the series previously reported by our group. Additional properties have been determined to assess their potential as anti-malarial leads including; HepG2 cytotoxicity, solubility, permeability, and lipophilicity, as well as in vitro stability in human and rat microsomes. We also assess their inhibition profile against a diverse set of 10 human kinases. Molecular docking, cheminformatics and bioinformatics analyses were also undertaken. Compounds 40 demonstrated the best anti-malarial activity at Pf 3D7 (0.09 µM), good selectivity with respect to mammalian cytotoxicity (SI = 54) and low microsomal clearance. Quantitative structure activity relationship (QSAR) analyses point to lipophilicity being a key driver of improved anti-malarial activity. The most active compounds in the series suffered from high lipophilicity, poor aqueous solubility and low permeability. The results provide useful information to guide further chemistry iterations.

Quinazoline-based human phosphodiesterase 5 inhibitors exhibited a selective vasorelaxant effect on rat isolated pulmonary arteries involving NO-sGC-cGMP pathway and calcium inhibitory effects.

Phosphodiesterase 5 (PDE5) inhibitors are an attractive option among the currently available therapies in the management of pulmonary arterial hypertension (PAH). Good selectivity for PDE5 is associated with reduced side effects and greater vasorelaxant effect on pulmonary arteries (PA). This study investigated the vasorelaxant effects of a series of quinazoline-based PDE5 inhibitors and their precise mechanisms action using rat isolated PA and aorta, as compared to sildenafil. Their effects on rat hepatocytes (viability and CYP activities) were also evaluated. Compounds 5 and 11 displayed lower human PDE5 IC50 of the analogs studied here and induced a greater relaxant effect on PA (EC50 0.94 ± 0.30 and 1.03 ± 0.23 µM, respectively). As compared to sildenafil (EC50 = 0.05 ± 0.02 µM on PA), the relaxant effect of 5 and 11 on PA was lower but their selectivity for PA compared to aorta was higher. The effects of 5 and 11 were reduced by NG-nitro-L-arginine methyl ester, 1H-[1,2,4]oxadiazolo[4,3-a]quinoxaline-1-one, but not by indomethacin or potassium channels blockers. They also enhanced the relaxant effect of sodium nitroprusside, and inhibited extracellular Ca2+ influx and intracellular Ca2+release. Compounds 5 and 11 did not reduce hepatocyte viability except at concentration > 10 µM, inhibited CYP3A at 10 µM, like sildenafil, but did not induce CYP1A. In conclusion, this study identified 2 quinazoline analogues with good PDE5 inhibitory activity and good selectivity for the pulmonary vasculature. Their relaxant effect involves both the potentiation of nitric oxide-sGC-cGMP pathway and calcium inhibition. These compounds are potential leads for developing new drugs for PAH.

The Effect of the EGFR - Targeting Compound 3-[(4-Phenylpyrimidin-2-yl) Amino] Benzene-1-Sulfonamide (13f) against Cholangiocarcinoma Cell Lines.

OBJECTIVE: Cholangiocarcinoma (CCA) is a noxious malignancy of epithelium of the bile duct with a low response rate to chemotherapy. The epidermal growth factor receptor (EGFR) signaling pathway is implicated in the development of cancerous cells, especially CCA. In this study, we report detailed biological profiling of 13f identified from our earlier hit expansion studies. The aim of this work was to expand our understanding of 13f via more detailed investigations of its mechanism of action against KKU-100, KKU-452 and KKU-M156 CCA cells, as well as in comparison to the EGFR inhibitor Gefitinib and non-specific chemotherapeutic agents such as Cisplatin. METHODS: Inhibiting EGFR-Kinase, cytotoxicity, clonogenic assay, wound healing and apoptosis were performed. Levels of total expression of EGFR and EGFR phosphorylation proteins were detected. RESULTS: 13f was confirmed as an inhibitor of EGFR with an IC50 value against the tyrosine kinase of EGFR of 22 nM and IC50 values for 48 h incubation period were 1.3 ± 1.9, 1.5 ± 0.4 and 1.7 ± 1.1 µM of KKU-100, KKU-452 and KKU-M156, respectively through dose- and time-dependent induction of early apoptosis of CCA cells. The compound also suppressed the clonogenic ability of KKU-100 and KKU-M156 cells stronger than Gefitinib, while potently inhibiting EGF-stimulated CCA cell migratory activity in KKU-452 cells. It was observed that under normal conditions EGFR was activated in CCA cells. EGF-stimulated basal expression of EGFR in KKU-452 cells was suppressed following 13f treatment, which was significantly greater than that of the marketed EGFR inhibitor Gefitinib. CONCLUSION: In summary, our study showed that 13f has potent anti-cancer activities including antiproliferation, clonogenic ability and migration through the modulation of EGFR signaling pathway in CCA for the first time. The compound represents an interesting starting point as a potential chemotherapeutic agent in ongoing efforts to improve response rate in CCA patients.
.

Theoretical Evaluation of the Reaction Mechanism of Serine Hydroxymethyltransferase.

Serine hydroxymethyltransferase (SHMT) is a pyridoxal phosphate (PLP)-dependent enzyme that catalyzes the reversible conversion of serine and tetrahydrofolate (THF) to glycine and 5,10-methylene THF. SHMT is a folate pathway enzyme and is therefore of considerable medical interest due to its role as an important intervention point for antimalarial, anticancer, and antibacterial treatments. Despite considerable experimental effort, the precise reaction mechanism of SHMT remains unclear. In this study, we explore the mechanism of SHMT to determine the roles of active site residues and the nature and the sequence of chemical steps. Molecular dynamics (MD) methods were employed to generate a suitable starting structure which then underwent analysis using hybrid quantum mechanical/molecular mechanical (QM/MM) simulations. The QM region consisted of 12 key residues, two substrates, and explicit solvent. Our results show that the catalytic reaction proceeds according to a retro-aldol synthetic process with His129 acting as the general base in the reaction. The rate-determining step involves the cleavage of the PLP-serine aldimine Cα-Cß bond and the formation of formaldehyde in line with experimental evidence. The pyridyl ring of the PLP-serine aldimine substrate exists in deprotonated form, being stabilized directly by Asp208 via a strong H-bond, as well as through interactions with Arg371, Lys237, and His211, and with the surrounding protein which was electrostatically embedded. This knowledge has the potential to impact the design and development of new inhibitors.

Synthesis, Plasmodium falciparum Inhibitory Activity, Cytotoxicity and Solubility of N2 ,N4 -Disubstituted Quinazoline-2,4-diamines.

BACKGROUND: Despite the development of extensive control strategies and treatment options, approximately 200 million malaria cases, leading to approximately 450,000 deaths, were reported in 2015. Due to issue of disease resistance, additional drug development efforts are needed to produce new, more effective treatments. Quinazoline-2,4-diamines were identified as antiparasitic compounds over three decades ago and have remained of interest to date in industry and academia. OBJECTIVE: An anti-malarial SAR evaluation of previously unreported N2 ,N4 -disubstituted quinazoline- 2,4-diamines have been undertaken in this study. We have synthesized and evaluated new derivatives against P. falciparum in our attempt to better characterize their biological activity and overall physical properties. METHODS: The synthesis of N2 ,N4 -disubstituted quinazoline-2,4-diamines inhibitors is reported along with activities in a radioactive labeled hypoxanthine incorporation assay against the f Plasmodium falciparum (Pf.) K1 strain. In addition, cytotoxicity was determined in the A549 and Vero cell lines using an MTT based. The aqueous solubility of key compounds was assessed at pH 7.4 using a shake flask-based approach. RESULTS: We identified compounds 1 and 6p as sub µM inhibitors of P. falciparum, having equivalent anti-malarial activity to Chloroquine. Compounds 1 and 6m are low µM inhibitors of P. falciparum with improved cytotoxicity profiles. Compound 6m displayed the best balance between P. falciparum Inhibitory activity (2 µM) and cytotoxicity, displaying >49 fold selectivity over A549 and Vero cell lines. CONCLUSION: Twenty one N2 ,N4 -Disubstituted Quinazoline-2,4-diamines have been prepared in our group and characterized in terms of their antimalarial activity, cytotoxicity and physical properties. Compounds with good activity and reasonable selectivity over mammalian cell lines have been identified. SAR analyses suggest further exploration is are necessary to improve the balance of P. falciparum Inhibitory activity, cytotoxicity and solubility.