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.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Shikimate Kinase (EC 2.7.1.71) from Mycobacterium tuberculosis: Kinetics and Structural Dynamics of a Potential Molecular Target for Drug Development

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

Whydration effects on DNA double helix stability modulates ligand binding to natural DNA in response to changes in water activity

In this work we present evidence that water molecules are actively involved on the control of binding affinity and binding site discrimination of a drug to natural DNA. In a previous study, the effect of water activity (a(w)) on the energetic parameters of actinomycin-D intercalation to natural DNA was determined using the osmotic stress method (39). This earlier study has shown evidence that water molecules act as an allosteric regulator of ligand binding to DNA via the effect of water activity on the long-range stability of the DNA secondary structure. In this work we have carried out DNA circularization experiments using the plasmid pUC18 in the absence of drugs and in the presence of different neutral solutes to evaluate the contribution of water activity to the energetics of DNA helix unwinding. The contribution of water to these independent reactions were made explicit by the description of how the changes in the free energy of ligand binding to DNA and in the free energy associated with DNA helix torsional deformation are linked to a(w) via changes in structural hydration. Taken together, the results of these studies reveal an extensive linkage between ligand binding affinity and site binding discrimination, and long range helix conformational changes and DNA hydration, This is strong evidence that water molecules work as a classical allosteric regulator of ligand binding to the DNA via its contribution to the stability of the double helix secondary structure, suggesting a possible mechanism by which the biochemical machinery of DNA processing takes advantage of the low activity of water into the cellular milieu.

Mapping eIF5A binding sites for Dys1 and Lia1: In vivo evidence for regulation of eIF5A hypusination

The evolutionarily conserved factor eIF5A is the only protein known to undergo hypusination, a unique posttranslational modification triggered by deoxyhypusine synthase (Dys1). Although eIF5A is essential for cell viability, the function of this putative translation initiation factor is still obscure. To identify eIF5A-binding proteins that could clarify its function, we screened a two-hybrid library and identified two eIF-5A partners in S. cerevisiae: Dys1 and the protein encoded by the gene YJR070C, named Lia1 (Ligand of eIF5A). The interactions were confirmed by GST pulldown. Mapping binding sites for these proteins revealed that both eIF5A domains can bind to Dys1, whereas the C-terminal domain is sufficient to bind Lia1. We demonstrate for the first time in vivo that the N-terminal α-helix of Dys1 can modulate enzyme activity by inhibiting eIF5A interaction. We suggest that this inhibition be abrogated in the cell when hypusinated and functional eIF5A is required. © 2003 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

Probing the binding of the coumarin drugs using peptide fragments of DNA gyrase B protein

Bacterial DNA gyrase, has been identified as the target of several antibacterial agents, including the coumarin drugs. The coumarins inhibit the gyrase action by competitive binding to the ATP-binding site of DNA gyrase B (GyrB) protein. The high in vitro inhibitory potency of coumarins against DNA gyrase reactions has raised interest in studies on coumarin-gyrase interactions. In this context, a series of low-molecular weight peptides, including the coumarin resistance-determining region of subunit B of Escherichia coli gyrase, has been designed and synthesized. The first peptide model was built using the natural fragment 131-146 of GyrB and was able to bind to novobiocin (K a = 1.8 ± 0.2 × 105/M) and ATP (Ka = 1.9 ± 0.4 × 103/M). To build the other sequences, changes in the Arg136 residue were introduced so that the binding to the drug was progressively reduced with the hydrophobicity of this residue (Ka = 1.3 ± 0.1 × 105/M and 1.0 ± 0.2 × 105/M for Ser and His, respectively). No binding was observed for the change Arg136 to Leu. In contrast, the binding to ATP was not altered, independently of the changes promoted. On the contrary, for peptide-coumarin and peptide-ATP complexes, Mg2+ appears to modulate the binding process. Our results demonstrate the crucial role of Arg 136 residue for the stability of coumarin-gyrase complex as well as suggest a different binding site for ATP and in both cases the interactions are mediated by magnesium ions. Copyright Blackwell Munksgaard, 2005.

Few substitutions affect the bioluminescence spectra of Phrixotrix (Coleoptera: Phengodidae) luciferases: a site-directed mutagenesis survey

Phrixotrix (railroad worm) luciferases produce bioluminescence in the green and red regions of the spectrum, depending on the location of the lanterns, and are the only luciferases naturally producing red bioluminescence. Comparison of the luciferase sequences showed a set of substitutions that could be involved in bioluminescence colour determination: (a) unique substitutions in the red luciferase replacing otherwise invariant residues; (b) conserved basic residues in the green-yellow emitting luciferases; and (c) an additional R353 residue in red-emitting luciferase (Viviani et al., 1999). To investigate whether these sites have a functional role in bioluminescence colour determination, we performed a site-directed mutagenesis. Natural substitutions in the region 220-344 and residues in the putative luciferin-binding site were also investigated. With the exception of the previously identified substitution of R215 and T226 (Viviani et al., 2002), which display dramatic red-shift effects on the spectrum of green-yellow-emitting luciferases, only a few substitutions had a moderate effect on the spectrum of the green-emitting luciferase. In contrast, no single substitution affected the spectrum of the red-emitting luciferase. The results suggest that the identity of the active site residues is not so critical for determining red bioluminescence in PxRE luciferase. Rather, the conformation assumed during the emitting step could be critical to set up proper interactions with excited oxyluciferin. Copyright ©2007 John Wiley & Sons, Ltd.

Interaction of Bacillus thuringiensis Cry1 and Vip3A proteins with Spodoptera frugiperda midgut binding sites

Vip3Aa, Vip3Af, Cry1Ab, and Cry1Fa were tested for their toxicities and binding interactions. Vip3A proteins were more toxic than Cry1 proteins. Binding assays showed independent specific binding sites for Cry1 and Vip3A proteins. Cry1Ab and Cry1Fa competed for the same binding sites, whereas Vip3Aa competed for those of Vip3Af. Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Novel aryl beta-aminocarbonyl derivatives as inhibitors of Trypanosoma cruzi trypanothione reductase: binding mode revised by docking and GRIND2-based 3D-QSAR procedures

Trypanothione reductase has long been investigated as a promising target for chemotherapeutic intervention in Chagas disease, since it is an enzyme of a unique metabolic pathway that is exclusively present in the pathogen but not in the human host, which has the analog Glutathione reductase. In spite of the present data-set includes a small number of compounds, a combined use of flexible docking, pharmacophore perception, ligand binding site prediction, and Grid-Independent Descriptors GRIND2-based 3D-Quantitative Structure-Activity Relationships (QSAR) procedures allowed us to rationalize the different biological activities of a series of 11 aryl beta-aminocarbonyl derivatives, which are inhibitors of Trypanosoma cruzi trypanothione reductase (TcTR). Three QSAR models were built and validated using different alignments, which are based on docking with the TcTR crystal structure, pharmacophore, and molecular interaction fields. The high statistical significance of the models thus obtained assures the robustness of this second generation of GRIND descriptors here used, which were able to detect the most important residues of such enzyme for binding the aryl beta-aminocarbonyl derivatives, besides to rationalize distances among them. Finally, a revised binding mode has been proposed for our inhibitors and independently supported by the different methodologies here used, allowing further optimization of the lead compounds with such combined structure- and ligand-based approaches in the fight against the Chagas disease.

DNA-Interactive Properties of Crotamine, a Cell-Penetrating Polypeptide and a Potential Drug Carrier

Crotamine, a 42-residue polypeptide derived from the venom of the South American rattlesnake Crotalus durissus terrificus, has been shown to be a cell-penetrating protein that targets chromosomes, carries plasmid DNA into cells, and shows specificity for actively proliferating cells. Given this potential role as a nucleic acid-delivery vector, we have studied in detail the binding of crotamine to single- and double-stranded DNAs of different lengths and base compositions over a range of ionic conditions. Agarose gel electrophoresis and ultraviolet spectrophotometry analysis indicate that complexes of crotamine with long-chain DNAs readily aggregate and precipitate at low ionic strength. This aggregation, which may be important for cellular uptake of DNA, becomes less likely with shorter chain length. 25-mer oligonucleotides do not show any evidence of such aggregation, permitting the determination of affinities and size via fluorescence quenching experiments. The polypeptide binds non-cooperatively to DNA, covering about 5 nucleotide residues when it binds to single (ss) or (ds) double stranded molecules. The affinities of the protein for ss-vs. ds-DNA are comparable, and inversely proportional to salt levels. Analysis of the dependence of affinity on [NaCl] indicates that there are a maximum of,3 ionic interactions between the protein and DNA, with some of the binding affinity attributable to non-ionic interactions. Inspection of the three-dimensional structure of the protein suggests that residues 31 to 35, Arg-Trp-Arg-Trp-Lys, could serve as a potential DNA-binding site. A hexapeptide containing this sequence displayed a lower DNA binding affinity and salt dependence as compared to the full-length protein, likely indicative of a more suitable 3D structure and the presence of accessory binding sites in the native crotamine. Taken together, the data presented here describing crotamine-DNA interactions may lend support to the design of more effective nucleic acid drug delivery vehicles which take advantage of crotamine as a carrier with specificity for actively proliferating cells. Citation: Chen P-C, Hayashi MAF, Oliveira EB, Karpel RL (2012) DNA-Interactive Properties of Crotamine, a Cell-Penetrating Polypeptide and a Potential Drug Carrier. PLoS ONE 7(11): e48913. doi:10.1371/journal.pone.0048913

Relative positioning of diazepam in the benzodiazepine-binding-pocket of GABA receptors

GABA(A) receptors are the major inhibitory neurotransmitter receptors in the brain. Some of them are targets of benzodiazepines that are widely used in clinical practice for their sedative/hypnotic, anxiolytic, muscle relaxant and anticonvulsant effects. In order to rationally separate these different drug actions, we need to understand the interaction of such compounds with the benzodiazepine-binding pocket. With this aim, we mutated residues located in the benzodiazepine-binding site individually to cysteine. These mutated receptors were combined with benzodiazepine site ligands carrying a cysteine reactive group in a defined position. Proximal apposition of reaction partners will lead to a covalent reaction. We describe here such proximity-accelerated chemical coupling reactions of alpha(1)S205C and alpha(1)T206C with a diazepam derivative modified at the C-3 position with a reactive isothiocyanate group (-NCS). We also provide new data that identify alpha(1)H101C and alpha(1)N102C as exclusive sites of the reaction of a diazepam derivative where the -Cl atom is replaced by a -NCS group. Based on these observations we propose a relative positioning of diazepam within the benzodiazepine-binding site of alpha(1)beta(2)gamma(2) receptors.

Molecular analysis of the site for 2-arachidonylglycerol (2-AG) on the β 2 subunit of GABA A receptors

2-arachidonyl glycerol (2-AG) allosterically potentiates GABAA receptors via a binding site located in transmembrane segment M4 of the β2 subunit. Two amino acid residues have been described that are essential for this effect. With the aim to further describe this potential drug target, we performed a cysteine scanning of the entire M4 and part of M3. All four residues in M4 affecting the potentiation here and the two already identified residues locate to the same side of the α-helix. This side is exposed to M3, where further residues were identified. From the fact that the important residues span > 18 Å, we conclude that the hydrophobic tail of the bound 2-AG molecule must be near linear and that the site mainly locates to the inner leaflet but stretches far into the membrane. The influence of the structure of the head group of the ligand molecule on the activity of the molecule was also investigated. We present a model of 2-AG docked to the GABAA receptor.

Structural studies of cytochrome P4501A1: Identification of substrate and cumene hydroperoxide binding regions in the active site of P4501A1 and characterization of the P4501A1 interaction with cytochrome P450 reductase

Cytochromes P450 are a superfamily of heme-thiolate proteins that function in a concert with another protein, cytochrome P450 reductase, as terminal oxidases of an enzymatic system catalyzing the metabolism of a variety of foreign compounds and endogenous substrates. In order to better understand P450s catalytic mechanism and substrate specificity, information about the structure of the active site is necessary. Given the lack of a crystal structure of mammalian P450, other methods have been used to elucidate the substrate recognition and binding site structure in the active center. In this project I utilized the photoaffinity labeling technique and site-directed mutagenesis approach to gain further structural insight into the active site of mammalian cytochrome P4501AI and examine the role of surface residues in the interaction of P4501A1 with the reductase. ^ Four crosslinked peptides were identified by photoaffinity labeling using diazido benzphetamine as a substrate analog. Alignment of the primary structure of cytochrome P4501A1 with that of bacterial cytochrome P450102 (the crystal structure of which is known) revealed that two of the isolated crosslinked peptides can be placed in the vicinity of heme (in the L helix region and β10-β11 sheet region of cytochrome P450102) and could be involved in substrate binding. The other two peptides were located on the surface of the protein with the label bound specifically to Lys residues that were proposed to be involved in reductase-P450 interaction. ^ Alternatively, it has been shown that some of the organic hydroperoxides can support P450 catalyzed reactions in the absence of NADPH, O2 and reductase. By means of photoaffinity labeling the cumene hydroperoxide binding region was identified. Using azidocumene as the photoaffinity label, the tripeptide T501-L502-K503 was shown to be the site where azidocumene covalently binds to P4501A1. The sequence alignment of cytochrome P4501A1 with cytochrome P450102 predicts that this region might correspond to β-sheet structure localized on the distal side of the heme ring near the I helix and the oxygen binding pocket. The role of Thr501 in the cumene hydroperoxide binding was confirmed by mutations of this residue and kinetic analysis of the effects of the mutations. ^ In addition, the role of two lysine residues, Lys271 and Lys279, in the interaction with reductase was examined by means of site-directed mutagenesis. The lysine residues were substituted with isoleucine and enzymatic activity of the wild type and the mutants were compared in reductase- and cumene hydroperoxide-supported systems. The lysine 279 residue has been shown to play a critical role in the P4501A1-reductase interaction. ^

Synthesis of photo-crosslinking probes and their application for the site-selective chemical modification of the 5-HT3 receptor

The 5-HT3 receptor (5-HT3R) is an important ion channel responsible for the transmission of nerve impulses in the CNS and PNS that is activated by the endogenous agonist serotonin (5-hydroxytryptamine, 5-HT). 5-HT3R is the only serotonin receptor belonging to the Cys-loop superfamily of neurotransmitter receptors. Different structural biology approaches can be applied, such as crystallization and x-ray analysis. Nonetheless, characterizing the exact ligand binding site(s) of these dynamic receptors is still challenging. The use of photo-crosslinking probes is an alternative validated approach allowing identification of regions in the protein that are important for the binding of small molecules. We designed our probes based on the core structure of the 5-HT3R antagonist granisetron, a FDA approved drug used for the treatment of chemotherapy-induced nausea and vomiting. We synthesized a small library of photo-crosslinking probes by conjugating diazirines and benzophenones via various linkers to granisetron. We were able to obtain several compounds with diverse linker lengths and different photo-crosslinking moieties that show nanomolar binding affinity for the orthosteric binding site. Furthermore we established a stable h5-HT3R expressing cell line and a purification protocol to yield the receptor in a high purity. Several experiments showed unambiguously that we are able to photo-crosslink our probes with the receptor site-specifically. The functionalised protein was analysed by Western blot and MS-analysis. This yielded the exact covalent modification site, corroborating current ligand binding models derived from mutagenesis and docking studies.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.