Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Novel synthesised flavone derivatives provide significant insight into the structural features required for enhanced anti-proliferative activity

With many cancers showing resistance to current chemotherapies, the search for novel anti-cancer agents is attracting considerable attention. Natural flavonoids have been identified as useful leads in such programmes. However, since an in-depth understanding of the structural requirements for optimum activity is generally lacking, further research is required before the full potential of flavonoids as anti-proliferative agents can be realised. Herein a broad library of 76 methoxy and hydroxy flavones, and their 4-thio analogues, was constructed and their structure-activity relationships for anti-proliferative activity against the breast cancer cell lines MCF-7 (ER+ve), MCF-7/DX (ER+ve, anthracycline resistant) and MDA-MB-231 (ER-ve) were probed. Within this library, 42 compounds were novel, and all compounds were afforded in good yields and > 95% purity. The most promising lead compounds, specifically the novel hydroxy 4-thioflavones 15f and 16f, were further evaluated for their anti-proliferative activities against a broader range of cancer cell lines by the National Cancer Institute (NCI), USA and displayed significant growth inhibition profiles (e.g Compound-15f: MCF-7 (GI50 = 0.18 μM), T-47D (GI50 = 0.03 μM) and MDA-MB-468 (GI50 = 0.47 μM) and compound-16f: MCF-7 (GI50 = 1.46 μM), T-47D (GI50 = 1.27 μM) and MDA-MB-231 (GI50 = 1.81 μM). Overall, 15f and 16f exhibited 7-46 fold greater anti-proliferative potency than the natural flavone chrysin (2d). A systematic structure-activity relationship study against the breast cancer cell lines highlighted that free hydroxyl groups and the B-ring phenyl groups were essential for enhanced anti-proliferative activities. Substitution of the 4-C=O functionality with a 4-C=S functionality, and incorporation of electron withdrawing groups at C4’ of the B-ring phenyl, also enhanced activity. Molecular docking and mechanistic studies suggest that the anti-proliferative effects of flavones 15f and 16f are mediated via ER-independent cleavage of PARP and downregulation of GSK-3β for MCF-7 and MCF-7/DX cell lines. For the MDA-MB-231 cell line, restoration of the wild-type p53 DNA binding activity of mutant p53 tumour suppressor gene was indicated.

GRID and docking analyses reveal a molecular basis for flavonoid inhibition of src-family kinase activity

Flavonoids reduce cardiovascular disease risk through anti-inflammatory, anti-coagulant and anti-platelet actions. One key flavonoid inhibitory mechanism is blocking kinase activity that drives these processes. Flavonoids attenuate activities of kinases including phosphoinositide-3-kinase (PI3K), Fyn, Lyn, Src, Syk, PKC, PIM1/2, ERK, JNK, and PKA. X-ray crystallographic analyses of kinase-flavonoid complexes show that flavonoid ring systems and their hydroxyl substitutions are important structural features for their binding to kinases. A clearer understanding of structural interactions of flavonoids with kinases is necessary to allow construction of more potent and selective counterparts. We examined flavonoid (quercetin, apigenin and catechin) interactions with Src-family kinases (Lyn, Fyn and Hck) applying the Sybyl docking algorithm and GRID. A homology model (Lyn) was used in our analyses to demonstrate that high quality predicted kinase structures are suitable for flavonoid computational studies. Our docking results revealed potential hydrogen bond contacts between flavonoid hydroxyls and kinase catalytic site residues. Identification of plausible contacts indicated that quercetin formed the most energetically stable interactions, apigenin lacked hydroxyl groups necessary for important contacts, and the non-planar structure of catechin could not support predicted hydrogen bonding patterns. GRID analysis using a hydroxyl functional group supported docking results. Based on these findings, we predicted that quercetin would inhibit activities of Src-family kinases with greater potency than apigenin and catechin. We validated this prediction using in vitro kinase assays. We conclude that our study can be used as a basis to construct virtual flavonoid interaction libraries to guide drug discovery using these compounds as molecular templates.

Molecular insights of p47phox phosphorylation dynamics in the regulation of NADPH oxidase activation and superoxide production

Phagocyte superoxide production by a multicomponent NADPH oxidase is important in host defense against microbial invasion. However inappropriate NADPH oxidase activation causes inflammation. Endothelial cells express NADPH oxidase and endothelial oxidative stress due to prolonged NADPH oxidase activation predisposes many diseases. Discovering the mechanism of NADPH oxidase activation is essential for developing novel treatment of these diseases. The p47phox is a key regulatory subunit of NADPH oxidase; however, due to the lack of full protein structural information, the mechanistic insight of p47phox phosphorylation in NADPH oxidase activation remains incomplete. Based on crystal structures of three functional domains, we generated a computational structural model of the full p47phox protein. Using a combination of in silico phosphorylation, molecular dynamics simulation and protein/protein docking, we discovered that the C-terminal tail of p47phox is critical for stabilizing its autoinhibited structure. Ser-379 phosphorylation disrupts H-bonds that link the C-terminal tail to the autoinhibitory region (AIR) and the tandem Src homology 3 (SH3) domains, allowing the AIR to undergo phosphorylation to expose the SH3 pocket for p22phox binding. These findings were confirmed by site-directed mutagenesis and gene transfection of p47phox_/_ coronary microvascular cells. Compared with wild-type p47phoxcDNAtransfected cells, the single mutation of S379A completely blocked p47phox membrane translocation, binding to p22phox and endothelial O2 . production in response to acute stimulation of PKC. p47phox C-terminal tail plays a key role in stabilizing intramolecular interactions at rest. Ser-379 phosphorylation is a molecular switch which initiates p47phox conformational changes and NADPH oxidase-dependent superoxide production by cells.