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.

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.

Immersive Projection for biochemical pattern matching

In all biological processes, protein molecules and other small molecules interact to function and form transient macromolecular complexes. This interaction of two or more molecules can be described by a docking event. Docking is an important phase for structure-based drug design strategies, as it can be used as a method to simulate protein-ligand interactions. Various docking programs exist that allow automated docking, but most of them have limited visualization and user interaction. It would be advantageous if scientists could visualize the molecules participating in the docking process, manipulate their structures and manually dock them before submitting the new conformations to an automated docking process in an immersive environment, which can help stimulate the design/docking process. This also could greatly reduce docking time and resources. To achieve this, we propose a new virtual modelling/docking program, whereby the advantages of virtual modelling programs and the efficiency of the algorithms in existing docking programs will be merged.

Effect of sulphation on the oestrogen agonist activity of the phytoestrogens genistein and daidzein in MCF-7 human breast cancer cells

The phytoestrogens genistein, daidzein and the daidzein metabolite equol have been shown previously to possess oestrogen agonist activity. However, following consumption of soya diets, they are found in the body not only as aglycones but also as metabolites conjugated at their 4'- and 7-hydroxyl groups with sulphate. This paper describes the effects of monosulphation on the oestrogen agonist properties of these three phytoestrogens in MCF-7 human breast cancer cells in terms of their relative ability to compete with [H-3]oestradiol for binding to oestrogen receptor (ER), to induce a stably transfected oestrogen-responsive reporter gene (ERE-CAT) and to stimulate cell growth. In no case did sulphation abolish activity. The 4'-sulphadon of genistein reduced oestrogen agonist activity to a small extent in whole-cell assays but increased the relative binding affinity to ER. The 7-sulphation of genistein, and also of equol, reduced oestrogen agonist activity substantially in all assays. By contrast, the position of monosulphation of daidzein acted in an opposing manner on oestrogen agonist activity. Sulphation at the 4'-position of daidzein resulted in a modest reduction in oestrogen agonist activity but sulphation of daidzein at the 7-position resulted in an increase in oestrogen agonist activity. Molecular modelling and docking studies suggested that the inverse effects of sulphation could be explained by the binding of daidzein into the ligand-binding domain of the ER in the opposite orientation compared with genistein and equol. This is the first report of sulphation enhancing activity of an isoflavone and inverse effects of sulphation between individual phytoestrogens.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Titanocene anticancer complexes and their binding mode of action to human serum albumin: a computational study

Due to the pivotal role played by human serum albumin (HSA) in the transport and cytotoxicity of titanocene complexes, a docking study has been performed on a selected set of titanocene complexes to aid in the current understanding of the potential mode of action of these titanocenes upon binding HSA. Analysis of the docking results has revealed potential binding at the known drug binding sites in HSA and has provided some explanation for the specificity and subsequent cytotoxicity of these titanocenes. Additionally, a new alternative binding site for these titanocenes has been postulated.

A novel, non-canonical mechanism of regulation of MST3 (mammalian Sterile20-related kinase 3)

The canonical pathway of regulation of the germinal centre kinase (GCK) III subgroup member, mammalian Sterile20-related kinase 3 (MST3), involves a caspase-mediated cleavage between N-terminal catalytic and C-terminal regulatory domains with possible concurrent autophosphorylation of the activation loop MST3(Thr178-), induction of Ser-/Thr-protein kinase activity and nuclear localisation. We identified an alternative ‘non-canonical’ pathway of MST3 activation (regulated primarily through dephosphorylation) which may also be applicable to other GCKIII (and GCKVI) subgroup members. In the basal state, inactive MST3 co-immunoprecipitated with the Golgi protein, GOLGA2/gm130. Activation of MST3 by calyculin A (a protein Ser-/Thr- phosphatase 1/2A inhibitor) stimulated (auto)phosphorylation of MST3(Thr178-) in the catalytic domain with essentially simultaneous cis-autophosphorylation of MST3(Thr328-) in the regulatory domain, an event also requiring the MST3(341-376) sequence which acts as a putative docking domain. MST3(Thr178-) phosphorylation increased MST3 kinase activity but this activity was independent of MST3(Thr328-) phosphorylation. Interestingly, MST3(Thr328-) lies immediately C-terminal to a STRAD pseudokinase-like site recently identified as being involved in binding of GCKIII/GCKVI members to MO25 scaffolding proteins. MST3(Thr178- /Thr328-) phosphorylation was concurrent with dissociation of MST3 from GOLGA2/gm130 and association of MST3 with MO25, and MST3(Thr328-) phosphorylation was necessary for formation of the activated MST3-MO25 holocomplex.

Prediction of novel and analogous folds using fragment assembly and fold recognition

A number of new and newly improved methods for predicting protein structure developed by the Jones–University College London group were used to make predictions for the CASP6 experiment. Structures were predicted with a combination of fold recognition methods (mGenTHREADER, nFOLD, and THREADER) and a substantially enhanced version of FRAGFOLD, our fragment assembly method. Attempts at automatic domain parsing were made using DomPred and DomSSEA, which are based on a secondary structure parsing algorithm and additionally for DomPred, a simple local sequence alignment scoring function. Disorder prediction was carried out using a new SVM-based version of DISOPRED. Attempts were also made at domain docking and “microdomain” folding in order to build complete chain models for some targets.

Characterization of Grb2-binding proteins in human platelets activated by Fc gamma RIIA cross-linking.

Glutathione-S-transferase (GST)-Grb2 fusion proteins have been used to identify the potential role of Grb2-binding proteins in platelet activation by the platelet low-affinity IgG receptor, Fc gamma RIIA. Two tyrosine phosphoproteins of 38 and 63 kD bind to the SH2 domain of Grb2 following Fc gamma RIIA stimulation of platelets. Both are located in the particulate fraction following platelet activation and are also able to bind to a GST-construct containing the SH2 and SH3 domains of phospholipase C gamma 1. p38 also forms a complex with the tyrosine kinase csk in stimulated cells and is a substrate for the kinase. The SH3 domains of Grb2 form a stable complex with SOS1 and two proteins of 75 kD and 120 kD, which undergo tyrosine phosphorylation in Fc gamma RIIA stimulated cells. The 75-kD protein is recognized by antibodies to SLP-76, which has recently been isolated from T cells and sequenced. Tyrosine phosphorylation of p38 and p63 is also observed in platelets stimulated by the tyrosine kinase-linked receptor agonist collagen and by the G protein-coupled receptor agonist thrombin, although phosphorylation of SLP-76 is only observed in collagen-stimulated platelets. p38 and p63 may provide a docking site for Grb2, thereby linking Grb2 SH3-binding proteins SOS1, SLP-76, and p120 to downstream signalling events.

A novel transport mechanism for MOMP in Chlamydophila pneumoniae and its putative role in immune-therapy

Major outer membrane proteins (MOMPs) of Gram negative bacteria are one of the most intensively studied membrane proteins. MOMPs are essential for maintaining the structural integrity of bacterial outer membranes and in adaptation of parasites to their hosts. There is evidence to suggest a role for purified MOMP from Chlamydophila pneumoniae and corresponding MOMP-derived peptides in immune-modulation, leading to a reduced atherosclerotic phenotype in apoE−/− mice via a characteristic dampening of MHC class II activity. The work reported herein tests this hypothesis by employing a combination of homology modelling and docking to examine the detailed molecular interactions that may be responsible. A three-dimensional homology model of the C. pneumoniae MOMP was constructed based on the 14 transmembrane β-barrel crystal structure of the fatty acid transporter from Escherichia coli, which provides a plausible transport mechanism for MOMP. Ligand docking experiments were used to provide details of the possible molecular interactions driving the binding of MOMP-derived peptides to MHC class II alleles known to be strongly associated with inflammation. The docking experiments were corroborated by predictions from conventional immuno-informatic algorithms. This work supports further the use of MOMP in C. pneumoniae as a possible vaccine target and the role of MOMP-derived peptides as vaccine candidates for immune-therapy in chronic inflammation that can result in cardiovascular events.

CLEC-2 is not required for platelet aggregation at arteriolar shear

The C-type lectin receptor CLEC-2 is expressed primarily on the surface of platelets, where it is present as a dimer, and is found at low level on a subpopulation of other hematopoietic cells, including mouse neutrophils [1–4] Clustering of CLEC-2 by the snake venom toxin rhodocytin, specific antibodies or its endogenous ligand, podoplanin, elicits powerful activation of platelets through a pathway that is similar to that used by the collagen receptor glycoprotein VI (GPVI) [4–6]. The cytosolic tail of CLEC-2 contains a conserved YxxL sequence preceded by three upstream acidic amino acid residues, which together form a novel motif known as a hemITAM. Ligand engagement induces tyrosine phosphorylation of the hemITAM sequence providing docking sites for the tandem-SH2 domains of the tyrosine kinase Syk across a CLEC-2 receptor dimer [3]. Tyrosine phosphorylation of Syk by Src family kinases and through autophosphorylation leads to stimulation of a downstream signaling cascade that culminates in activation of phospholipase C γ2 (PLCγ2) [4,6]. Recently, CLEC-2 has been proposed to play a major role in supporting activation of platelets at arteriolar rates of flow [1]. Injection of a CLEC-2 antibody into mice causes a sustained depletion of the C-type lectin receptor from the platelet surface [1]. The CLEC-2-depleted platelets were unresponsive to rhodocytin but underwent normal aggregation and secretion responses after stimulation of other platelet receptors, including GPVI [1]. In contrast, there was a marked decrease in aggregate formation relative to controls when CLEC-2-depleted blood was flowed at arteriolar rates of shear over collagen (1000 s−1 and 1700 s−1) [1]. Furthermore, antibody treatment significantly increased tail bleeding times and mice were unable to occlude their vessels after ferric chloride injury [1]. These data provide evidence for a critical role for CLEC-2 in supporting platelet aggregation at arteriolar rates of flow. The underlying mechanism is unclear as platelets do not express podoplanin, the only known endogenous ligand of CLEC-2. In the present study, we have investigated the role of CLEC-2 in platelet aggregation and thrombus formation using platelets from a novel mutant mouse model that lacks functional CLEC-2.

A global proteomics approach identifies novel phosphorylated signaling proteins in GPVI-activated platelets: involvement of G6f, a novel platelet Grb2-binding membrane adapter.

Collagen-related peptide (CRP) stimulates powerful activation of platelets through the glycoprotein VI (GPVI)-FcR gamma-chain complex. We have combined proteomics and traditional biochemistry approaches to study the proteome of CRP-activated platelets, focusing in detail on tyrosine phosphorylation. In two separate approaches, phosphotyrosine immunoprecipitations followed by 1-D-PAGE, and 2-DE, were used for protein separation. Proteins were identified by MS. By following these approaches, 96 proteins were found to undergo PTM in response to CRP in human platelets, including 11 novel platelet proteins such as Dok-1, SPIN90, osteoclast stimulating factor 1, and beta-Pix. Interestingly, the type I transmembrane protein G6f was found to be specifically phosphorylated on Tyr-281 in response to platelet activation by CRP, providing a docking site for the adapter Grb2. G6f tyrosine phoshporylation was also found to take place in response to collagen, although not in response to the G protein-coupled receptor agonists, thrombin and ADP. Further, we also demonstrate for the first time that Grb2 and its homolog Gads are tyrosine-phosphorylated in CRP-stimulated platelets. This study provides new insights into the mechanism of platelet activation through the GPVI collagen receptor, helping to build the basis for the development of new drug targets for thrombotic disease.

Comparison of CLEC-2 and GPVI signaling in platelets: the role of adaptor proteins

GPVI activates platelets through an ITAM pathway by activation of Src and Syk kinases leading to activation of PLCy2. CLEC-2 has been shown to activate platelets using an ITAM-like sequence in its cytoplasmic tail that is also dependent on Src and Syk kinases, but shows a partial rather than an absolute dependence on adapter SLP-76 for activation of PLCy2. The aim of this thesis is to understand some of the key differences in these signalling pathways. GPVI is in complex with FcRwhich contains the ITAM sequence (Yxx(L/I)x6−12Yxx(L/I)). These two tyrosines provide a docking site for the tandem-SH2 domains of Syk. In this thesis I show that CLEC-2 signalling through Syk is mediated by phosphorylation of the CLEC-2 YxxL sequence, receptor dimerisation and cross-linking by the Syk SH2 domains. I also show that the differential requirement for SLP-76 is not mediated by Gads. Both signalling pathways also show partial dependency for LAT. I also show that a novel protein, G6f, is not able to substitute for LAT in this signalling pathway and also exclude the LAT-family proteins PAG, LIME, LAX and NTAL as potential LAT replacements in platelet activation by GPVI. These results extend our understanding of platelet activation by CLEC-2.