scispace - formally typeset
Search or ask a question

Showing papers on "Docking (molecular) published in 2012"


Journal ArticleDOI
TL;DR: Of the various structure-based approaches that have been applied to soluble protein targets, such as proteases and kinases, in silico docking is among the most ready applicable to GPCRs.

266 citations


Journal ArticleDOI
TL;DR: Based on this Capture-SELEX procedure, the successful DNA aptamer selection for the aminoglycoside antibiotic kanamycin A as a small molecule target is described.
Abstract: Small organic molecules are challenging targets for an aptamer selection using the SELEX technology (SELEX—Systematic Evolution of Ligans by EXponential enrichment). Often they are not suitable for immobilization on solid surfaces, which is a common procedure in known aptamer selection methods. The Capture-SELEX procedure allows the selection of DNA aptamers for solute targets. A special SELEX library was constructed with the aim to immobilize this library on magnetic beads or other surfaces. For this purpose a docking sequence was incorporated into the random region of the library enabling hybridization to a complementary oligo fixed on magnetic beads. Oligonucleotides of the library which exhibit high affinity to the target and a secondary structure fitting to the target are released from the beads for binding to the target during the aptamer selection process. The oligonucleotides of these binding complexes were amplified, purified, and immobilized via the docking sequence to the magnetic beads as the starting point of the following selection round. Based on this Capture-SELEX procedure, the successful DNA aptamer selection for the aminoglycoside antibiotic kanamycin A as a small molecule target is described.

201 citations


Journal ArticleDOI
TL;DR: Molecular dynamics simulation was used to examine the binding of nine substrates, two inhibitors, and two nonsubstrates to the distal binding pocket of AcrB, identified earlier by X-ray crystallography and found that all substrates tested bound to the pocket, whereas the binding to this site was not preferred for the nonssubstrates.
Abstract: Multidrug resistance in Gram-negative bacteria, to which multidrug efflux pumps such as the AcrB transporter makes a major contribution, is becoming a major public health problem Unfortunately only a few compounds have been cocrystallized with AcrB, and thus computational approaches are essential in elucidating the interaction between diverse ligands and the pump protein We used molecular dynamics simulation to examine the binding of nine substrates, two inhibitors, and two nonsubstrates to the distal binding pocket of AcrB, identified earlier by X-ray crystallography This approach gave us more realistic views of the binding than the previously used docking approach, as the explicit water molecules contributed to the process and the flexible binding site was often seen to undergo large structural changes We analyzed the interaction in detail in terms of the binding energy, hydrophobic surface-matching, and the residues involved in the process We found that all substrates tested bound to the pocket, whereas the binding to this site was not preferred for the nonsubstrates Interestingly, both inhibitors [Phe-Arg-β-naphthylamide and 1-(1-naphtylmethyl)-piperazine] tended to move out of the pocket at least partially, getting into contact with a glycine-rich loop that separates the distal pocket from the more proximal region of the protein and is thought to control the access of substrates to the distal pocket

194 citations


Journal ArticleDOI
TL;DR: It is shown that, surprisingly, in spite of the limited number of protein–protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the knownprotein–protein interactions, provided the components themselves have a known structure or can be homology-built.
Abstract: Traditional approaches to protein–protein docking sample the binding modes with no regard to similar experimentally determined structures (templates) of protein–protein complexes. Emerging template-based docking approaches utilize such similar complexes to determine the docking predictions. The docking problem assumes the knowledge of the participating proteins’ structures. Thus, it provides the possibility of aligning the structures of the proteins and the template complexes. The progress in the development of template-based docking and the vast experience in template-based modeling of individual proteins show that, generally, such approaches are more reliable than the free modeling. The key aspect of this modeling paradigm is the availability of the templates. The current common perception is that due to the difficulties in experimental structure determination of protein–protein complexes, the pool of docking templates is insignificant, and thus a broad application of template-based docking is possible only at some future time. The results of our large scale, systematic study show that, surprisingly, in spite of the limited number of protein–protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the known protein–protein interactions, provided the components themselves have a known structure or can be homology-built. About one-third of the templates are of good quality when they are compared to experimental structures in test sets extracted from the Protein Data Bank and would be useful starting points in modeling the complexes. This finding dramatically expands our ability to model protein interactions, and has far-reaching implications for the protein docking field in general.

181 citations


Journal ArticleDOI
TL;DR: A coherent structural model for MAPK docking specificity is suggested that reveals how short linear motifs binding to a common kinase docking groove can mediate diverse interaction patterns and contribute to correct MAPK partner selection in signaling networks.
Abstract: Mitogen-activated protein kinases (MAPKs) have a docking groove that interacts with linear "docking" motifs in binding partners. To determine the structural basis of binding specificity between MAPKs and docking motifs, we quantitatively analyzed the ability of 15 docking motifs from diverse MAPK partners to bind to c-Jun amino-terminal kinase 1 (JNK1), p38α, and extracellular signal-regulated kinase 2 (ERK2). Classical docking motifs mediated highly specific binding only to JNK1, and only those motifs with a sequence pattern distinct from the classical MAPK binding docking motif consensus differentiated between the topographically similar docking grooves of ERK and p38α. Crystal structures of four complexes of MAPKs with docking peptides, representing JNK-specific, ERK-specific, or ERK- and p38-selective binding modes, revealed that the regions located between consensus positions in the docking motifs showed conformational diversity. Although the consensus positions in the docking motifs served as anchor points that bound to common MAPK surface features and mostly contributed to docking in a nondiscriminatory fashion, the conformation of the intervening region between the anchor points mostly determined specificity. We designed peptides with tailored MAPK binding profiles by rationally changing the length and amino acid composition of intervening regions located between anchor points. These results suggest a coherent structural model for MAPK docking specificity that reveals how short linear motifs binding to a common kinase docking groove can mediate diverse interaction patterns and contribute to correct MAPK partner selection in signaling networks.

146 citations


Journal ArticleDOI
TL;DR: A web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with recently developed scoring functions based on robust regression analysis and quantum chemical charge models.
Abstract: Identification of possible protein targets of small chemical molecules is an important step for unravelling their underlying causes of actions at the molecular level. To this end, we construct a web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with our recently developed scoring functions based on robust regression analysis and quantum chemical charge models. Affinity profiles of the protein targets are used to provide the confidence levels of prediction. The divide-and-conquer docking approach uses adaptively constructed small overlapping grids to constrain the searching space, thereby achieving better docking efficiency. Unlike previous approaches that screen against a specific class of targets or a limited number of targets, idTarget screen against nearly all protein structures deposited in the Protein Data Bank (PDB). We show that idTarget is able to reproduce known off-targets of drugs or drug-like compounds, and the suggested new targets could be prioritized for further investigation. idTarget is freely available as a web-based server at http://idtarget.rcas.sinica.edu.tw.

143 citations


Journal ArticleDOI
TL;DR: It is reported that peroxisome proliferator-activated receptor γ (PPARγ) modulators can sample multiple binding modes manifesting in multiple receptor conformations in slow conformational exchange and that ligand and receptor dynamics affect the graded transcriptional output of PPARγ modulators.

136 citations


Journal ArticleDOI
TL;DR: The potency and efficiency of these molecules has few precedents among protein–protein interface inhibitors, and supports structure-based efforts to discover leads for chemokine GPCRs.
Abstract: G-protein–coupled receptors (GPCRs) are key signaling molecules and are intensely studied. Whereas GPCRs recognizing small-molecules have been successfully targeted for drug discovery, protein-recognizing GPCRs, such as the chemokine receptors, claim few drugs or even useful small molecule reagents. This reflects both the difficulties that attend protein–protein interface inhibitor discovery, and the lack of structures for these targets. Imminent structure determination of chemokine receptor CXCR4 motivated docking screens for new ligands against a homology model and subsequently the crystal structure. More than 3 million molecules were docked against the model and then against the crystal structure; 24 and 23 high-scoring compounds from the respective screens were tested experimentally. Docking against the model yielded only one antagonist, which resembled known ligands and lacked specificity, whereas the crystal structure docking yielded four that were dissimilar to previously known scaffolds and apparently specific. Intriguingly, several were potent and relatively small, with IC50 values as low as 306 nM, ligand efficiencies as high as 0.36, and with efficacy in cellular chemotaxis. The potency and efficiency of these molecules has few precedents among protein–protein interface inhibitors, and supports structure-based efforts to discover leads for chemokine GPCRs.

136 citations


Journal ArticleDOI
20 Apr 2012-PLOS ONE
TL;DR: This work successfully determined the amino acid residues responsible for binding to sweeteners in the cleft of hT1R2 ATD and suggests that individual ligands have sets of specific residues for binding in correspondence with the chemical structures and other residuesresponsible for interacting with multiple ligands.
Abstract: One of the most distinctive features of human sweet taste perception is its broad tuning to chemically diverse compounds ranging from low-molecular-weight sweeteners to sweet-tasting proteins. Many reports suggest that the human sweet taste receptor (hT1R2–hT1R3), a heteromeric complex composed of T1R2 and T1R3 subunits belonging to the class C G protein–coupled receptor family, has multiple binding sites for these sweeteners. However, it remains unclear how the same receptor recognizes such diverse structures. Here we aim to characterize the modes of binding between hT1R2–hT1R3 and low-molecular-weight sweet compounds by functional analysis of a series of site-directed mutants and by molecular modeling–based docking simulation at the binding pocket formed on the large extracellular amino-terminal domain (ATD) of hT1R2. We successfully determined the amino acid residues responsible for binding to sweeteners in the cleft of hT1R2 ATD. Our results suggest that individual ligands have sets of specific residues for binding in correspondence with the chemical structures and other residues responsible for interacting with multiple ligands.

126 citations


Journal ArticleDOI
18 May 2012-PLOS ONE
TL;DR: The high selectivity and low concentration dependent inhibitory effect of Auphen (in the nanomolar range) together with its high water solubility makes the compound a suitable drug lead for future in vivo studies.
Abstract: Aquaporins (AQPs) are membrane channels that conduct water and small solutes such as glycerol and are involved in many physiological functions. Aquaporin-based modulator drugs are predicted to be of broad potential utility in the treatment of several diseases. Until today few AQP inhibitors have been described as suitable candidates for clinical development. Here we report on the potent inhibition of AQP3 channels by gold(III) complexes screened on human red blood cells (hRBC) and AQP3-transfected PC12 cells by a stopped-flow method. Among the various metal compounds tested, Auphen is the most active on AQP3 (IC50 = 0.8±0.08 µM in hRBC). Interestingly, the compound poorly affects the water permeability of AQP1. The mechanism of gold inhibition is related to the ability of Au(III) to interact with sulphydryls groups of proteins such as the thiolates of cysteine residues. Additional DFT and modeling studies on possible gold compound/AQP adducts provide a tentative description of the system at a molecular level. The mapping of the periplasmic surface of an homology model of human AQP3 evidenced the thiol group of Cys40 as a likely candidate for binding to gold(III) complexes. Moreover, the investigation of non-covalent binding of Au complexes by docking approaches revealed their preferential binding to AQP3 with respect to AQP1. The high selectivity and low concentration dependent inhibitory effect of Auphen (in the nanomolar range) together with its high water solubility makes the compound a suitable drug lead for future in vivo studies. These results may present novel metal-based scaffolds for AQP drug development.

123 citations


Journal ArticleDOI
TL;DR: This work combines fluorescence spectroscopy, NMR experiments, docking and molecular dynamics simulations to study the configurational and recognition properties of IL-8 towards a series of HA and CS derivatives and DS, and the effects of GAG length and sulfation patterns in binding strength and specificity.
Abstract: The interactions between glycosaminoglycans (GAGs), important components of the extracellular matrix, and proteins such as growth factors and chemokines play critical roles in cellular regulation processes. Therefore, the design of GAG derivatives for the development of innovative materials with bio-like properties in terms of their interaction with regulatory proteins is of great interest for tissue engineering and regenerative medicine. Previous work on the chemokine interleukin-8 (IL-8) has focused on its interaction with heparin and heparan sulfate, which regulate chemokine function. However, the extracellular matrix contains other GAGs, such as hyaluronic acid (HA), dermatan sulfate (DS) and chondroitin sulfate (CS), which have so far not been characterized in terms of their distinct molecular recognition properties towards IL-8 in relation to their length and sulfation patterns. NMR and molecular modeling have been in great part the methods of choice to study the structural and recognition properties of GAGs and their protein complexes. However, separately these methods have challenges to cope with the high degree of similarity and flexibility that GAGs exhibit. In this work, we combine fluorescence spectroscopy, NMR experiments, docking and molecular dynamics simulations to study the configurational and recognition properties of IL-8 towards a series of HA and CS derivatives and DS. We analyze the effects of GAG length and sulfation patterns in binding strength and specificity, and the influence of GAG binding on IL-8 dimer formation. Our results highlight the importance of combining experimental and theoretical approaches to obtain a better understanding of the molecular recognition properties of GAG-protein systems.

Journal ArticleDOI
TL;DR: Modelling protein–protein association using the associative-memory, water-mediated, structure, and energy model (AWSEM), a coarse-grained protein folding model that has been optimized using energy-landscape theory, shows the ability of the model to predict interfaces and addresses the relative importance of monomer geometry, flexibility, and nonnative intermonomeric contacts in the association process for the homodimers.
Abstract: We investigate protein–protein association using the associative-memory, water-mediated, structure, and energy model (AWSEM), a coarse-grained protein folding model that has been optimized using energy-landscape theory. The potential was originally parameterized by enforcing a funneled nature for a database of dimeric interfaces but was later further optimized to create funneled folding landscapes for individual monomeric proteins. The ability of the model to predict interfaces was not tested previously. The present results show that simulated annealing of the model indeed is able to predict successfully the native interfaces of eight homodimers and four heterodimers, thus amounting to a flexible docking algorithm. We go on to address the relative importance of monomer geometry, flexibility, and nonnative intermonomeric contacts in the association process for the homodimers. Monomer surface geometry is found to be important in determining the binding interface, but it is insufficient. Using a uniform binding potential rather than the water-mediated potential results in sampling of misbound structures that are geometrically preferred but are nonetheless energetically disfavored by AWSEM, as well as in nature. Depending on the stability of the unbound monomers, nonnative contacts play different roles in the association process. For unstable monomers, thermodynamic states stabilized by nonnative interactions correspond to productive, on-pathway intermediates and can, therefore, catalyze binding through a fly-casting mechanism. For stable monomers, in contrast, states stabilized by nonnative interactions generally correspond to traps that impede binding.

Journal ArticleDOI
TL;DR: The mechanism of binding of a new series of potent inhibitors of Adenosine Deaminase is developed and revealed and the lowest energy binding modes of the most potent compound of the series, 4-decyl-pyrazolo[1,5-a]pyrimidin-7-one are identified.
Abstract: An exhaustive description of the molecular recognition mechanism between a ligand and its biological target is of great value because it provides the opportunity for an exogenous control of the related process. Very often this aim can be pursued using high resolution structures of the complex in combination with inexpensive computational protocols such as docking algorithms. Unfortunately, in many other cases a number of factors, like protein flexibility or solvent effects, increase the degree of complexity of ligand/protein interaction and these standard techniques are no longer sufficient to describe the binding event. We have experienced and tested these limits in the present study in which we have developed and revealed the mechanism of binding of a new series of potent inhibitors of Adenosine Deaminase. We have first performed a large number of docking calculations, which unfortunately failed to yield reliable results due to the dynamical character of the enzyme and the complex role of the solvent. Thus, we have stepped up the computational strategy using a protocol based on metadynamics. Our approach has allowed dealing with protein motion and solvation during ligand binding and finally identifying the lowest energy binding modes of the most potent compound of the series, 4-decyl-pyrazolo[1,5-a]pyrimidin-7-one.

Journal ArticleDOI
TL;DR: Molecular docking result, along with the biological assay data, suggested that compound 5c was a potential anti-inflammatory agent.

Journal ArticleDOI
14 May 2012-PLOS ONE
TL;DR: A new inhibitory mechanism of ACE by a peptide is provided, which broads the knowledge for drug designing against enzyme targets, and is supported by the structural model.
Abstract: Angiotensin I-converting enzyme (ACE) has an important function in blood pressure regulation. ACE-inhibitory peptides can lower blood pressure by inhibiting ACE activity. Based on the sequence of an ACE-inhibitory hexapeptide (TPTQQS) purified from yeast, enzyme kinetics experiments, isothermal titration calorimetry (ITC), and a docking simulation were performed. The hexapeptide was found to inhibit ACE in a non-competitive manner, as supported by the structural model. The hexapeptide bound to ACE via interactions of the N-terminal Thr1, Thr3, and Gln4 residues with the residues on the lid structure of ACE, and the C-terminal Ser6 attracted the zinc ion, which is vital for ACE catalysis. The displacement of the zinc ion from the active site resulted in the inhibition of ACE activity. The structural model based on the docking simulation was supported by experiments in which the peptide was modified. This study provides a new inhibitory mechanism of ACE by a peptide which broads our knowledge for drug designing against enzyme targets.

Journal ArticleDOI
TL;DR: Docking studies of selected analogues into a homology model of P-glycoprotein suggest that benzophenones show an interaction pattern similar to that previously identified for propafenone-type inhibitors.
Abstract: The drug efflux pump P-glycoprotein (P-gp) has been shown to promote multidrug resistance (MDR) in tumors as well as to influence ADME properties of drug candidates. Here we synthesized and tested a series of benzophenone derivatives structurally analogous to propafenone-type inhibitors of P-gp. Some of the compounds showed ligand efficiency and lipophilic efficiency (LipE) values in the range of compounds which entered clinical trials as MDR modulators. Interestingly, although lipophilicity plays a dominant role for P-gp inhibitors, all compounds investigated showed LipE values below the threshold for promising drug candidates. Docking studies of selected analogues into a homology model of P-glycoprotein suggest that benzophenones show an interaction pattern similar to that previously identified for propafenone-type inhibitors.

Journal ArticleDOI
TL;DR: A homology model of nsP2 protein based on the crystal structure of the nsP1 protein of Venezuelan equine encephalitis virus is presented, and the pharmacophore features that must be present in an inhibitor of ns P2 protease are proposed.
Abstract: To date, no suitable vaccine or specific antiviral drug is available to treat Chikungunya viral (CHIKV) fever. Hence, it is essential to identify drug candidates that could potentially impede CHIKV infection. Here, we present the development of a homology model of nsP2 protein based on the crystal structure of the nsP2 protein of Venezuelan equine encephalitis virus (VEEV). The protein modeled was optimized using molecular dynamics simulation; the junction peptides of a nonstructural protein complex were then docked in order to investigate the possible protein-protein interactions between nsP2 and the proteins cleaved by nsP2. The modeling studies conducted shed light on the binding modes, and the critical interactions with the peptides provide insight into the chemical features needed to inhibit the CHIK virus infection. Energy-optimized pharmacophore mapping was performed using the junction peptides. Based on the results, we propose the pharmacophore features that must be present in an inhibitor of nsP2 protease. The resulting pharmacophore model contained an aromatic ring, a hydrophobic and three hydrogen-bond donor sites. Using these pharmacophore features, we screened a large public library of compounds (Asinex, Maybridge, TOSLab, Binding Database) to find a potential ligand that could inhibit the nsP2 protein. The compounds that yielded a fitness score of more than 1.0 were further subjected to Glide HTVS and Glide XP. Here, we report the best four compounds based on their docking scores; these compounds have IDs of 27943, 21362, ASN 01107557 and ASN 01541696. We propose that these compounds could bind to the active site of nsP2 protease and inhibit this enzyme. Furthermore, the backbone structural scaffolds of these four lead compounds could serve as building blocks when designing drug-like molecules for the treatment of Chikungunya viral fever.

Journal ArticleDOI
TL;DR: Target molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp and nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses.
Abstract: Multidrug resistance proteins that belong to the ATP-binding cassette family like the human P-glycoprotein (ABCB1 or Pgp) are responsible for many failed cancer and antiviral chemotherapies because these membrane transporters remove the chemotherapeutics from the targeted cells. Understanding the details of the catalytic mechanism of Pgp is therefore critical to the development of inhibitors that might overcome these resistances. In this work, targeted molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp. Crystal structures of homologues in four different conformations were used as intermediate targets in the dynamics simulations. Transitions from conformations that were wide open to the cytoplasm to transition state conformations that were wide open to the extracellular space were studied. Twenty-six nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses. Coupled movement of nucleotide binding domains (NBDs) and transmembrane domains (TMDs) that form the drug binding cavities were observed. Pronounced twisting of the NBDs as they approached each other as well as the quantification of a dramatic opening of the TMDs to the extracellular space as the ATP hydrolysis transition state was reached were observed. Docking interactions of 21 known transport ligands or inhibitors were analyzed with each of the 26 transitional structures. Many of the docking results obtained here were validated by previously published biochemical determinations. As the ATP hydrolysis transition state was approached, drug docking in the extracellular half of the transmembrane domains seemed to be destabilized as transport ligand exit gates opened to the extracellular space.

Journal ArticleDOI
TL;DR: The binding free energies obtained by the molecular mechanic-Poisson-Boltzmann surface area method indicate that curcumin displays higher binding affinity than nonsteroidal anti-inflammatory drugs.
Abstract: Binding of curcumin, naproxen, and ibuprofen to Aβ1–40 peptide and its fibrils is studied by docking method and all-atom molecular dynamics simulations. The Gromos96 43a1 force field and simple point charge model of water have been used for molecular dynamics simulations. It is shown that if the receptor is a monomer then naproxen and ibuprofen are bound to the same place that is different from the binding position of curcumin. However all of three ligands have the same binding pocket in fibrillar structures. The binding mechanism is studied in detail showing that the van der Waals interaction between ligand and receptor dominates over the electrostatic interaction. The binding free energies obtained by the molecular mechanic-Poisson–Boltzmann surface area method indicate that curcumin displays higher binding affinity than nonsteroidal anti-inflammatory drugs. Our results are in good agreement with the experiments.

Journal ArticleDOI
TL;DR: The identification of the first small molecules specifically designed to inhibit the HIV-1 replication by targeting the RNA binding site of the human DEAD-Box RNA helicase DDX3 is described.

Journal ArticleDOI
TL;DR: Clear correlations are observed between the docking accuracy and the similarity of the sequence of interest to the template, suggesting regimes in which docking can correctly identify ligand binding modes.
Abstract: The growing availability of novel structures for several G protein-coupled receptors (GPCRs) has provided new opportunities for structure-based drug design of ligands against this important class of targets. Here, we report a systematic analysis of the accuracy of docking small molecules into GPCR structures and homology models using both rigid receptor (Glide SP and Glide XP) and flexible receptor (Induced Fit Docking; IFD) methods. The ability to dock ligands into different structures of the same target (cross-docking) is evaluated for both agonist and inverse agonist structures of the A2A receptor and the β1- and β2-adrenergic receptors. In addition, we have produced homology models for the β1-adrenergic, β2-adrenergic, D3 dopamine, H1 histamine, M2 muscarine, M3 muscarine, A2A adenosine, S1P1, κ-opioid, and C-X-C chemokine 4 receptors using multiple templates and investigated the ability of docking to predict the binding mode of ligands in these models. Clear correlations are observed between the dock...

Journal ArticleDOI
TL;DR: PI-LZerD, a pairwise docking algorithm, which uses imperfect protein-protein binding interface prediction to improve docking accuracy, consistently showed better prediction accuracy over alternative methods in the series of benchmark experiments.
Abstract: Background Many important cellular processes are carried out by protein complexes. To provide physical pictures of interacting proteins, many computational protein-protein prediction methods have been developed in the past. However, it is still difficult to identify the correct docking complex structure within top ranks among alternative conformations.

Journal ArticleDOI
01 Jan 2012-Proteins
TL;DR: BSP‐SLIM was tested on 71 ligand–protein complexes from the Astex diverse set and was used in virtual screening for six target proteins, which prioritized actives of 25% and 50% in the top 9% and 17% of the library on average, respectively, demonstrating the usefulness of the template‐based coarse‐grained algorithms in the low‐resolution ligand‐protein docking and drug‐screening.
Abstract: We developed BSP-SLIM, a new method for ligand-protein blind docking using low-resolution protein structures. For a given sequence, protein structures are first predicted by I-TASSER; putative ligand binding sites are transferred from holo-template structures which are analogous to the I-TASSER models; ligand-protein docking conformations are then constructed by shape and chemical match of ligand with the negative image of binding pockets. BSP-SLIM was tested on 71 ligand-protein complexes from the Astex diverse set where the protein structures were predicted by I-TASSER with an average RMSD 2.92 A on the binding residues. Using I-TASSER models, the median ligand RMSD of BSP-SLIM docking is 3.99 A which is 5.94 A lower than that by AutoDock; the median binding-site error by BSP-SLIM is 1.77 A which is 6.23 A lower than that by AutoDock and 3.43 A lower than that by LIGSITE(CSC) . Compared to the models using crystal protein structures, the median ligand RMSD by BSP-SLIM using I-TASSER models increases by 0.87 A, while that by AutoDock increases by 8.41 A; the median binding-site error by BSP-SLIM increase by 0.69A while that by AutoDock and LIGSITE(CSC) increases by 7.31 A and 1.41 A, respectively. As case studies, BSP-SLIM was used in virtual screening for six target proteins, which prioritized actives of 25% and 50% in the top 9.2% and 17% of the library on average, respectively. These results demonstrate the usefulness of the template-based coarse-grained algorithms in the low-resolution ligand-protein docking and drug-screening. An on-line BSP-SLIM server is freely available at http://zhanglab.ccmb.med.umich.edu/BSP-SLIM.

Journal ArticleDOI
TL;DR: The data provide a mechanistic explanation for the allosteric action of antidepressants at SERT and suggest that the role of the vestibule is evolutionarily conserved among neurotransmitter:sodium symporter proteins as a binding pocket for small molecule ligands.

Journal ArticleDOI
TL;DR: Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCR structures and to analyze docking of other ligands in already crystallized GPCRs.
Abstract: Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A2A adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Journal ArticleDOI
Xian-Hui Yang1, Wen Qing1, Ting-Ting Zhao1, Jian Sun1, Xi Li1, Man Xing1, Xiang Lu1, Hai-Liang Zhu1 
TL;DR: Based on the preliminary results, compound 10e with potent inhibitory activity in tumor growth may be a potential anticancer agent.

Journal ArticleDOI
01 Jul 2012-Proteins
TL;DR: A novel computational multiple protein docking algorithm, Multi‐LZerD, that builds models of multimeric complexes by effectively reusing pairwise docking predictions of component proteins, and was able to predict near‐native structures for multimerics complexes of various topologies.
Abstract: The tertiary structures of protein complexes provide a crucial insight about the molecular mechanisms that regulate their functions and assembly. However, solving protein complex structures by experimental methods is often more difficult than single protein structures. Here, we have developed a novel computational multiple protein docking algorithm, Multi-LZerD, that builds models of multimeric complexes by effectively reusing pairwise docking predictions of component proteins. A genetic algorithm is applied to explore the conformational space followed by a structure refinement procedure. Benchmark on eleven hetero-multimeric complexes resulted in near-native conformations for all but one of them (a root mean square deviation smaller than 2.5A). We also show that our method copes with unbound docking cases well, outperforming the methodology that can be directly compared with our approach. Multi-LZerD was able to predict near-native structures for multimeric complexes of various topologies.

Journal ArticleDOI
TL;DR: This approach, termed HDX‐DOCK, gave more tightly clustered docking poses than stand‐alone docking for all antigen–antibody interactions examined and improved docking results significantly for the cytochrome c–E8 interaction.
Abstract: Understanding antigen-antibody interactions at the sub-molecular level is of particular interest for scientific, regulatory, and intellectual property reasons, especially with increasing demand for monoclonal antibody therapeutic agents. Although various techniques are available for the determination of an epitope, there is no widely applicable, high-resolution, and reliable method available. Here, a combination approach using amide hydrogen/deuterium exchange coupled with proteolysis and mass spectrometry (HDX-MS) and computational docking was applied to investigate antigen-antibody interactions. HDX-MS is a widely applicable, medium-resolution, medium-throughput technology that can be applied to epitope identification. First, the epitopes of cytochrome c-E8, IL-13-CNTO607, and IL-17A-CAT-2200 interactions identified using the HDX-MS method were compared with those identified by X-ray co-crystal structures. The identified epitopes are in good agreement with those identified using high-resolution X-ray crystallography. Second, the HDX-MS data were used as constraints for computational docking. More specifically, the non-epitope residues of an antigen identified using HDX-MS were designated as binding ineligible during computational docking. This approach, termed HDX-DOCK, gave more tightly clustered docking poses than stand-alone docking for all antigen-antibody interactions examined and improved docking results significantly for the cytochrome c-E8 interaction.

Journal ArticleDOI
TL;DR: The vicinity of receptor-bound C2 groups was probed by homology modeling based on recent X-ray structure of an agonist-bound A(2A)AR, with a predicted helical rearrangement requiring an agonists-specific outward displacement of TM2 resembling opsin.
Abstract: (N)-Methanocarba adenosine 5′-methyluronamides containing known A3 AR (adenosine receptor)-enhancing modifications, i.e., 2-(arylethynyl)adenine and N6-methyl or N6-(3-substituted-benzyl), were nanomolar full agonists of human (h) A3AR and highly selective (Ki ∼0.6 nM, N6-methyl 2-(halophenylethynyl) analogues 13 and 14). Combined 2-arylethynyl-N6-3-chlorobenzyl substitutions preserved A3AR affinity/selectivity in the (N)-methanocarba series (e.g., 3,4-difluoro full agonist MRS5698 31, Ki 3 nM, human and mouse A3) better than that for ribosides. Polyaromatic 2-ethynyl N6-3-chlorobenzyl analogues, such as potent linearly extended 2-p-biphenylethynyl MRS5679 34 (Ki hA3 3.1 nM; A1, A2A, inactive) and fluorescent 1-pyrene adduct MRS5704 35 (Ki hA3 68.3 nM), were conformationally rigid; receptor docking identified a large, mainly hydrophobic binding region. The vicinity of receptor-bound C2 groups was probed by homology modeling based on recent X-ray structure of an agonist-bound A2AAR, with a predicted helica...

Journal ArticleDOI
TL;DR: The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosIn inhibition studies and rutin has the potential to be a potent antipigment agent.
Abstract: Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values We simulated docking (binding energies for AutoDock Vina: -91 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin The docking results suggest that rutin mostly interacts with histidine residues located in the active site A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin The inhibition was a typical competitive type with K(i) = 110±025 mM Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced Taken together, rutin has the potential to be a potent anti-pigment agent The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors