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.

Docking and small angle X-ray scattering studies of purine nucleoside phosphorylase

Docking simulations have been used to assess protein complexes with some success. Small angle X-ray scattering (SAXS) is a well-established technique to investigate protein spatial configuration. This work describes the integration of geometric docking with SAXS to investigate the quaternary structure of recombinant human purine nucleoside phosphorylase (PNP). This enzyme catalyzes the reversible phosphorolysis of N-ribosidic bonds of purine nucleosides and deoxynucleosides. A genetic deficiency due to mutations in the gene encoding for PNP causes gradual decrease in T-cell immunity. Inappropriate activation of T-cells has been implicated in several clinically relevant human conditions such as transplant rejection, rheumatoid arthritis, lupus, and T-cell lymphomas. PNP is therefore a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. The present analysis confirms the trimeric structure observed in the crystal. The potential application of the present procedure to other systems is discussed. (C) 2003 Elsevier B.V. All rights reserved.

Automatic design of decision-tree induction algorithms tailored to flexible-receptor docking data

Background: This paper addresses the prediction of the free energy of binding of a drug candidate with enzyme InhA associated with Mycobacterium tuberculosis. This problem is found within rational drug design, where interactions between drug candidates and target proteins are verified through molecular docking simulations. In this application, it is important not only to correctly predict the free energy of binding, but also to provide a comprehensible model that could be validated by a domain specialist. Decision-tree induction algorithms have been successfully used in drug-design related applications, specially considering that decision trees are simple to understand, interpret, and validate. There are several decision-tree induction algorithms available for general-use, but each one has a bias that makes it more suitable for a particular data distribution. In this article, we propose and investigate the automatic design of decision-tree induction algorithms tailored to particular drug-enzyme binding data sets. We investigate the performance of our new method for evaluating binding conformations of different drug candidates to InhA, and we analyze our findings with respect to decision tree accuracy, comprehensibility, and biological relevance. Results: The empirical analysis indicates that our method is capable of automatically generating decision-tree induction algorithms that significantly outperform the traditional C4.5 algorithm with respect to both accuracy and comprehensibility. In addition, we provide the biological interpretation of the rules generated by our approach, reinforcing the importance of comprehensible predictive models in this particular bioinformatics application. Conclusions: We conclude that automatically designing a decision-tree algorithm tailored to molecular docking data is a promising alternative for the prediction of the free energy from the binding of a drug candidate with a flexible-receptor.

Prediction of heparin binding sites in bone morphogenetic proteins (BMPs)

Heparin is a glycosaminoglycan known to bind bone morphogenetic proteins (BMPs) and the growth and differentiation factors (GDFs) and has strong and variable effects on BMP osteogenic activity. In this paper we report our predictions of the likely heparin binding sites for BMP-2 and 14. The N-terminal sequences upstream of TGF-β-type cysteine-knot domains in BMP-2, 7 and 14 contain the basic residues arginine and lysine, which are key components of the heparin/HS-binding sites, with these residues being highly non-conserved. Importantly, evolutionary conserved surfaces on the beta sheets are required for interactions with receptors and antagonists. Furthermore, BMP-2 has electropositive surfaces on two sides compared to BMP-7 and BMP-14. Molecular docking simulations suggest the presence of high and low affinity binding sites in dimeric BMP-2. Histidines were found to play a role in the interactions of BMP-2 with heparin; however, a pKa analysis suggests that histidines are likely not protonated. This is indicative that interactions of BMP-2 with heparin do not require acidic pH. Taken together, non-conserved amino acid residues in the N-terminus and residues protruding from the beta sheet (not overlapping with the receptor binding sites and the dimeric interface) and not C-terminal are found to be important for heparin–BMP interactions.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

A variant peptide of buffalo colostrum beta-lactoglobulin inhibits angiotensin I-converting enzyme activity

beta-lactoglobulin is a rich source of bioactive peptides. The LC-MS separated tryptic peptides of buffalo colostrum beta-lactoglobulin (BLG-col) were computed based on MS-MS fragmentation for de novo sequencing. Among the selected peptides (P1-P8), a variant was detected with methionine at position 74 instead of glutamate. The sequences of two peptides were identical to hypocholesterolemic peptides whereas the remaining peptides were in accordance with buffalo milk beta-lactoglobulin. Comparative sequence analysis of BLG-col to milk beta-lactoglobulin was carried out using CLUSTALW2 and a molecular model for BLG-col was constructed (PMDB ID-PM0076812). The synthesized variant pentapeptide (IIAMK, m/z-576 Da) was found to inhibit angiotensin I-converting enzyme (ACE) with an IC50 of 498 +/- 2 mu M, which was rationalized through docking simulations using Molgrow virtual docker. (C) 2012 Elsevier Masson SAS. All rights reserved.

Identification and characterization of novel indole based small molecules as anticancer agents through SIRT1 inhibition

In our pursuit to develop new potential anticancer leads, we designed a combination of structural units of indole and substituted triazole; and a library of 1-{1-methyl-2-4-phenyl-5-(propan-2-ylsulfanyl)-4H-1,2,4-triazol-3-yl ]-1H-indol-3-yl}methanamine derivatives was synthesized and characterized. Cytotoxic evaluations of these molecules over a panel of three human cancer cell lines were carried out. Few molecules exhibited potent growth inhibitory action against the treated cancer cell lines at lower micro molar concentration. An in vitro assay investigation of these active compounds using recombinant human SIRT1 enzyme showed that one of the compounds (IT-14) inhibited the deacetylation activity of the enzyme. The in vivo study of IT-14 exemplified its promising action by reducing the prostate weight to the body weight ratio in prostate hyperplasia animal models. A remarkable decrease in the disruption of histoarchitecture of the prostate tissues isolated from IT-14 treated animal compared to that of the positive control was observed. The molecular interactions with SIRT1 enzyme were also supported by molecular docking simulations. Hence this compound can act as a lead molecule to treat prostatic hyperplasia. (C) 2013 Elsevier Masson SAS. All rights reserved.

The electron transfer complex between nitrous oxide reductase and its electron donors

J Biol Inorg Chem (2011) 16:1241–1254 DOI 10.1007/s00775-011-0812-9

Molecular modeling databases: A new way in the search of protein targets for drug development

DBMODELING is a relational database of annotated comparative protein structure models and their metabolic, pathway characterization. It is focused on enzymes identified in the genomes of Mycobacterium tuberculosis and Xylella fastidiosa. The main goal of the present database is to provide structural models to be used in docking simulations and drug design. However, since the accuracy of structural models is highly dependent on sequence identity between template and target, it is necessary to make clear to the user that only models which show high structural quality should be used in such efforts. Molecular modeling of these genomes generated a database, in which all structural models were built using alignments presenting more than 30% of sequence identity, generating models with medium and high accuracy. All models in the database are publicly accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools. DBMODELING user interface provides users friendly menus, so that all information can be printed in one stop from any web browser. Furthermore, DBMODELING also provides a docking interface, which allows the user to carry out geometric docking simulation, against the molecular models available in the database. There are three other important homology model databases: MODBASE, SWISSMODEL, and GTOP. The main applications of these databases are described in the present article. © 2007 Bentham Science Publishers Ltd.

SKPDB: A structural database of shikimate pathway enzymes

Background: The functional and structural characterisation of enzymes that belong to microbial metabolic pathways is very important for structure-based drug design. The main interest in studying shikimate pathway enzymes involves the fact that they are essential for bacteria but do not occur in humans, making them selective targets for design of drugs that do not directly impact humans.Description: The ShiKimate Pathway DataBase (SKPDB) is a relational database applied to the study of shikimate pathway enzymes in microorganisms and plants. The current database is updated regularly with the addition of new data; there are currently 8902 enzymes of the shikimate pathway from different sources. The database contains extensive information on each enzyme, including detailed descriptions about sequence, references, and structural and functional studies. All files (primary sequence, atomic coordinates and quality scores) are available for downloading. The modeled structures can be viewed using the Jmol program.Conclusions: The SKPDB provides a large number of structural models to be used in docking simulations, virtual screening initiatives and drug design. It is freely accessible at http://lsbzix.rc.unesp.br/skpdb/. © 2010 Arcuri et al; licensee BioMed Central Ltd.

Triacontyl p-coumarate: An inhibitor of snake venom metalloproteinases

Snake venom metalloproteinases (SVMPs) participate in a number of important biological, physiological and pathophysiological processes and are primarily responsible for the local tissue damage characteristic of viperid snake envenomations. The use of medicinal plant extracts as antidotes against animal venoms is an old practice, especially against snake envenomations. Such plants are sources of many pharmacologically active compounds and have been shown to antagonize the effects of some venoms and toxins. The present study explores the activity of triacontyl p-coumarate (PCT), an active compound isolated from root bark of Bombacopsis glabra vegetal extract (Bg), against harmful effects of Bothropoides pauloensis snake venom and isolated toxins (SVMPs or phospholipase A2). Before inhibition assays, Bg or PCT was incubated with venom or toxins at ratios of 1:1 and 1:5 (w/w; venom or isolated toxins/PCT) for 30 min at 37 °C. Treatment conditions were also assayed to simulate snakebite with PCT inoculated at either the same venom or toxin site. PCT neutralized fibrinogenolytic activity and plasmatic fibrinogen depletion induced by B. pauloensis venom or isolated toxin. PCT also efficiently inhibited the hemorrhagic (3MDH-minimum hemorrhagic dose injected i.d into mice) and myotoxic activities induced by Jararhagin, a metalloproteinase from B. jararaca at 1:5 ratio (toxin: inhibitor, w/w) when it was previously incubated with PCT and injected into mice or when PCT was administered after toxin injection. Docking simulations using data on a metalloproteinase (Neuwiedase) structure suggest that the binding between the protein and the inhibitor occurs mainly in the active site region causing blockade of the enzymatic reaction by displacement of catalytic water. Steric hindrance may also play a role in the mechanism since the PCT hydrophobic tail was found to interact with the loop associated with substrate anchorage. Thus, PCT may provide a alternative to complement ophidian envenomation treatments. © 2012 Elsevier Ltd. All rights reserved.

Modelagem molecular da Chiquimato quinase de Mycobacterium leprae

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Bioinformática estrutural aplicada ao estudo de proteínas alvo do genoma do Mycobacterium tuberculosis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Characterization of Heparin-induced Glyceraldehyde-3-phosphate Dehydrogenase Early Amyloid-like Oligomers and Their Implication in alpha-Synuclein Aggregation

Lewy bodies and Lewy neurites, neuropathological hallmarks of several neurological diseases, are mainly made of filamentous assemblies of alpha-synuclein. However, other macromolecules including Tau, ubiquitin, glyceraldehyde-3-phosphate dehydrogenase, and glycosaminoglycans are routinely found associated with these amyloid deposits. Glyceraldehyde-3-phosphate dehydrogenase is a glycolytic enzyme that can form fibrillar aggregates in the presence of acidic membranes, but its role in Parkinson disease is still unknown. In this work, the ability of heparin to trigger the amyloid aggregation of this protein at physiological conditions of pH and temperature is demonstrated by infrared and fluorescence spectroscopy, dynamic light scattering, small angle x-ray scattering, circular dichroism, and fluorescence microscopy. Aggregation proceeds through the formation of short rod-like oligomers, which elongates in one dimension. Heparan sulfate was also capable of inducing glyceraldehyde-3-phosphate dehydrogenase aggregation, but chondroitin sulfates A, B, and C together with dextran sulfate had a negligible effect. Aided with molecular docking simulations, a putative binding site on the protein is proposed providing a rational explanation for the structural specificity of heparin and heparan sulfate. Finally, it is demonstrated that in vitro the early oligomers present in the glyceraldehyde-3-phosphate dehydrogenase fibrillation pathway promote alpha-synuclein aggregation. Taking into account the toxicity of alpha-synuclein prefibrillar species, the heparin-induced glyceraldehyde-3-phosphate dehydrogenase early oligomers might come in useful as a novel therapeutic strategy in Parkinson disease and other synucleinopathies.

Looking over Toxin-K+ Channel Interactions. Clues from the Structural and Functional Characterization of alpha-KTx Toxin Tc32, a Kv1.3 Channel Blocker

alpha-KTx toxin Tc32, from the Amazonian scorpion Tityus cambridgei, lacks the dyad motif; including Lys27, characteristic of the family and generally associated with channel blockage. The toxin has been cloned and expressed for the first time. Electrophysiological experiments, by showing that the recombinant form blocks Kv1.3 channels of olfactory bulb periglomerular cells like the natural Tc32 toxin, when tested on the Kv1.3 channel of human T lymphocytes, confirmed it is in an active fold. The nuclear magnetic resonance-derived structure revealed it exhibits an alpha/beta scaffold typical of the members of the alpha-KTx family. TdK2 and TdK3, all belonging to the same alpha-KTx 18 subfamily, share significant sequence identity with Tc32 but diverse selectivity and affinity for Kv1.3 and Kv1.1 channels. To gain insight into the structural features that may justify those differences, we used the recombinant Tc32 nuclear magnetic resonance-derived structure to model the other two toxins, for which no experimental structure is available. Their interaction with Kv1.3 and Kv1.1 has been investigated by means of docking simulations. The results suggest that differences in the electrostatic features of the toxins and channels, in their contact surfaces, and in their total dipole moment orientations govern the affinity and selectivity of toxins. In addition, we found that, regardless of whether the dyad motif is present, it is always a Lys side chain that physically blocks the channels, irrespective of its position in the toxin sequence.