Updates to Binding MOAD (Mother of All Databases): Polypharmacology Tools and Their Utility in Drug Repurposing

The goal of Binding MOAD is to provide users with a data set focused on high-quality x-ray crystal structures that have been solved with biologically relevant ligands bound. Where available, experimental binding affinities (K_a, K_d, K_i, IC₅₀) are provided from the primary literature of the crystal structure. The database has been updated regularly since 2005, and this most recent update has added nearly 7000 new structures (growth of 21%). MOAD currently contains 32,747 structures, composed of 9117 protein families and 16,044 unique ligands. The data are freely available on www.BindingMOAD.org. This paper outlines updates to the data in Binding MOAD as well as improvements made to both the website and its contents. The NGL viewer has been added to improve visualization of the ligands and protein structures. MarvinJS has been implemented, over the outdated MarvinView, to work with JChem for small molecule searching in the database. To add tools for predicting polypharmacology, we have added information about sequence, binding-site, and ligand similarity between entries in the database. A main premise behind polypharmacology is that similar binding sites will bind similar ligands. The large amount of protein–ligand information available in Binding MOAD allows us to compute pairwise ligand and binding-site similarities. Lists of similar ligands and similar binding sites have been added to allow users to identify potential polypharmacology pairs. To show the utility of the polypharmacology data, we detail a few examples from Binding MOAD of drug repurposing targets with their respective similarities.     

DRoP: Automated detection of conserved solvent-binding sites on proteins

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.     

Landscape of protein–small ligand binding modes

Elucidating the mechanisms of specific small‐molecule (ligand) recognition by proteins is a long‐standing conundrum. While the structures of these molecules, proteins and ligands, have been extensively studied, protein–ligand interactions, or binding modes, have not been comprehensively analyzed. Although methods for assessing similarities of binding site structures have been extensively developed, the methods for the computational treatment of binding modes have not been well established. Here, we developed a computational method for encoding the information about binding modes as graphs, and assessing their similarities. An all‐against‐all comparison of 20,040 protein–ligand complexes provided the landscape of the protein–ligand binding modes and its relationships with protein‐ and chemical spaces. While similar proteins in the same SCOP Family tend to bind relatively similar ligands with similar binding modes, the correlation between ligand and binding similarities was not very high (R² = 0.443). We found many pairs with novel relationships, in which two evolutionally distant proteins recognize dissimilar ligands by similar binding modes (757,474 pairs out of 200,790,780 pairs were categorized into this relationship, in our dataset). In addition, there were an abundance of pairs of homologous proteins binding to similar ligands with different binding modes (68,217 pairs). Our results showed that many interesting relationships between protein–ligand complexes are still hidden in the structure database, and our new method for assessing binding mode similarities is effective to find them.     

Ligand binding to proteins: the binding landscape model.

Models of ligand binding are often based on four assumptions: (1) steric fit: that binding is determined mainly by shape complementarity; (2) native binding: that ligands mainly bind to native states; (3) locality: that ligands perturb protein structures mainly at the binding site; and (4) continuity: that small changes in ligand or protein structure lead to small changes in binding affinity. Using a generalization of the 2D HP lattice model, we study ligand binding and explore these assumptions. We first validate the model by showing that it reproduces typical binding behaviors. We observe ligand-induced denaturation, ANS and heme-like binding, and "lock-and-key" and "induced-fit" specific binding behaviors characterized by Michaelis-Menten or more cooperative types of binding isotherms. We then explore cases where the model predicts violations of the standard assumptions. For example, very different binding modes can result from two ligands of identical shape. Ligands can sometimes bind highly denatured states more tightly than native states and yet have Michaelis-Menten isotherms. Even low-population binding to denatured states can cause changes in global stability, hydrogen-exchange rates, and thermal B-factors, contrary to expectations, but in agreement with experiments. We conclude that ligand binding, similar to protein folding, may be better described in terms of energy landscapes than in terms of simpler mass-action models.     

Determination of Amino Acid Residues Responsible for Specific Interaction of Protein Kinases with Small Molecule Inhibitors

Identifying amino acid positions that determine the specific interaction of proteins with small molecule ligands, is required for search of pharmaceutical targets, drug design, and solution of other biotechnology problems. We studied applicability of an original method SPrOS (specificity projection on sequence) developed to recognize functionally significant positions in amino acid sequences. The method allows residues specific to functional subgroups to be determined within the protein family based on their local surroundings in amino acid sequences. The efficiency of the method has been estimated on the protein kinase family. The residues associated with the protein specificity to inhibitors have been predicted. The results have been verified using 3D structures of protein–ligand complexes. Three small molecule inhibitors have been tested. Residues predicted with SPrOS either in contacted the inhibitor or influenced the conformation of the ligand–binding area. Excluding close homologues from the studied set makes it possible to decrease the number of difficult to interpret positions. The expediency of this procedure was determined by the relationship between an inhibitory spectrum and phylogenic partition. Thus, the method efficiency has been confirmed by matching the prediction results with the protein 3D structures.     

Ligand-protein inverse docking and its potential use in the computer search of protein targets of a small molecule

Ligand-protein docking has been developed and used in facilitating new drug discoveries. In this approach, docking single or multiple small molecules to a receptor site is attempted to find putative ligands. A number of studies have shown that docking algorithms are capable of finding ligands and binding conformations at a receptor site close to experimentally determined structures. These algorithms are expected to be equally applicable to the identification of multiple proteins to which a small molecule can bind or weakly bind. We introduce a ligand-protein inverse-docking approach for finding potential protein targets of a small molecule by the computer-automated docking search of a protein cavity database. This database is developed from protein structures in the Protein Data Bank (PDB). Docking is conducted with a procedure involving multiple-conformer shape-matching alignment of a molecule to a cavity followed by molecular-mechanics torsion optimization and energy minimization on both the molecule and the protein residues at the binding region. Scoring is conducted by the evaluation of molecular-mechanics energy and, when applicable, by the further analysis of binding competitiveness against other ligands that bind to the same receptor site in at least one PDB entry. Testing results on two therapeutic agents, 4H-tamoxifen and vitamin E, showed that 50% of the computer-identified potential protein targets were implicated or confirmed by experiments. The application of this approach may facilitate the prediction of unknown and secondary therapeutic target proteins and those related to the side effects and toxicity of a drug or drug candidate. Proteins 2001;43:217-226.     

Cleavage of the iron-methionine bond in c-type cytochromes: Crystal structure of oxidized and reduced cytochrome c2 from Rhodopseudomonas palustris and its ammonia complex

The three-dimensional structures of the native cytochrome c(2) from Rhodopseudomonas palustris and of its ammonia complex have been obtained at pH 4.4 and pH 8.5, respectively. The structure of the native form has been refined in the oxidized state at 1.70 A and in the reduced state at 1.95 A resolution. These are the first high-resolution crystal structures in both oxidation states of a cytochrome c(2) with relatively high redox potential (+350 mV). The differences between the two oxidation states of the native form, including the position of internal water molecules, are small. The unusual six-residue insertion Gly82-Ala87, which precedes the heme binding Met93, forms an isolated 3(10)-helix secondary structural element not previously observed in other c-type cytochromes. Furthermore, this cytochrome shows an external methionine residue involved in a strained folding near the exposed edge of the heme. The structural comparison of the present cytochrome c(2) with other c-type cytochromes has revealed that the presence of such a residue, with torsion angles phi and psi of approximately -140 and -130 degrees, respectively, is a typical feature of this family of proteins. The refined crystal structure of the ammonia complex, obtained at 1.15 A resolution, shows that the sulphur atom of the Met93 axial ligand does not coordinate the heme iron atom, but is replaced by an exogenous ammonia molecule. This is the only example so far reported of an X-ray structure with the heme iron coordinated by an ammonia molecule. The detachment of Met93 is accompanied by a very localized change in backbone conformation, involving mainly the residues Lys92, Met93, and Thr94. Previous studies under typical denaturing conditions, including high-pH values and the presence of exogenous ligands, have shown that the detachment of the Met axial ligand is a basic step in the folding/unfolding process of c-type cytochromes. The ammonia adduct represents a structural model for this important step of the unfolding pathway. Factors proposed to be important for the methionine dissociation are the strength of the H-bond between the Met93 and Tyr66 residues that stabilizes the native form, and the presence in this bacterial cytochrome c(2) of the rare six-residue insertion in the helix 3(10) conformation that increases Met loop flexibility.     

In search of allosteric modulators of a7-nAChR by solvent density guided virtual screening

Nicotinic acetylcholine receptors (nAChR) are pentameric ligand gated ion channels whose activity can be modulated by endogenous neurotransmitters as well as by synthetic ligands that bind the same or distinct sites from the natural ligand. The subtype of α7 nAChR has been considered as a potenial therapeutic target for Alzheimer's disease, schizophrenia and other neurological and psychiatric disorders. Here we have developed a homology model of α7 nAChR based on two high resolution crystal structures with Brookhaven Protein Data Bank (PDB) codes 2QC1 and 2WN9 for threading on one monomer and then for building a pentamer, respectively. A number of small molecule binding sites are identified using Pocket Finder (J. An, M. Tortov, and R. Abagyan, Molecular & Cellular Proteomics, 4.6, 752-761 (2005)) of Internal Coordinate Mechanics (ICM). Remarkably, these computer-identified sites match perfectly with ordered solvent densities found in the high-resolution crystal structure of α1 nAChR, suggesting that the surface cavities in the α7 nAChR model are likely binding sites of small molecules. A high throughput virtual screening by flexible ligand docking of 5008 small molecule compounds was performed at three potential allosteric modulator (AM) binding sites of α7 nAChR using Molsoft ICM software (R. Abagyan, M. Tortov and D. Kuznetsov, J Comput Chem 15, 488-506, (1994)). Some experimentally verified allosteric modulators of α7 like CCMI comp-6, LY 7082101, 5-HI, TQS, PNU-120596, genistein, and NS-1738 ranked among top 100 compounds, while the rest of the compounds in the list could guide further search for new allosteric modulators.     

BI-D1870 is a specific inhibitor of the p90 RSK (ribosomal S6 kinase) isoforms in vitro and in vivo

Lipocalins are extracellular proteins (17-25 kDa) that bind and transport small lipophilic molecules. The three-dimensional structure of the first lipocalin from a metatherian has been determined at different values of pH both with and without bound ligands. Trichosurin, a protein from the milk whey of the common brushtail possum, Trichosurus vulpecula, has been recombinantly expressed in Escherichia coli, refolded from inclusion bodies, purified and crystallized at two different pH values. The three-dimensional structure of trichosurin was solved by X-ray crystallography in two different crystal forms to 1.9 A (1 A=0.1 nm) and 2.6 A resolution, from crystals grown at low and high pH values respectively. Trichosurin has the typical lipocalin fold, an eight-stranded anti-parallel beta-barrel but dimerizes in an orientation that has not been seen previously. The putative binding pocket in the centre of the beta-barrel is well-defined in both high and low pH structures and is occupied by water molecules along with isopropanol molecules from the crystallization medium. Trichosurin was also co-crystallized with a number of small molecule ligands and structures were determined with 2-naphthol and 4-ethylphenol bound in the centre of the beta-barrel. The binding of phenolic compounds by trichosurin provides clues to the function of this important marsupial milk protein, which is highly conserved across metatherians.     

Structural basis of non-specific lipid binding in maize lipid-transfer protein complexes revealed by high-resolution X-ray crystallography

Non-specific lipid-transfer proteins (nsLTPs) are involved in the movement of phospholipids, glycolipids, fatty acids, and steroids between membranes. Several structures of plant nsLTPs have been determined both by X-ray crystallography and nuclear magnetic resonance. However, the detailed structural basis of the non-specific binding of hydrophobic ligands by nsLTPs is still poorly understood. In order to gain a better understanding of the structural basis of the non-specific binding of hydrophobic ligands by nsLTPs and to investigate the plasticity of the fatty acid binding cavity in nsLTPs, seven high-resolution (between 1.3 A and 1.9 A) crystal structures have been determined. These depict the nsLTP from maize seedlings in complex with an array of fatty acids.A detailed comparison of the structures of maize nsLTP in complex with various ligands reveals a new binding mode in an nsLTP-oleate complex which has not been seen before. Furthermore, in the caprate complex, the ligand binds to the protein cavity in two orientations with equal occupancy. The volume of the hydrophobic cavity in the nsLTP from maize shows some variation depending on the size of the bound ligands.The structural plasticity of the ligand binding cavity and the predominant involvement of non-specific van der Waals interactions with the hydrophobic tail of the ligands provide a structural explanation for the non-specificity of maize nsLTP. The hydrophobic cavity accommodates various ligands from C10 to C18. The C18:1 ricinoleate with its hydroxyl group hydrogen bonding to Ala68 possibly mimics cutin monomer binding which is of biological importance. Some of the myristate binding sites in human serum albumin resemble the maize nsLTP, implying the importance of a helical bundle in accommodating the non-specific binding of fatty acids.     

Dodecins: a family of lumichrome binding proteins

Dodecin is a small dodecameric flavoprotein from Halobacterium salinarum that contains two flavins stacked between two tryptophan residues to form an aromatic tetrade. The functional properties of heterologously expressed dodecin were investigated by fluorescence spectroscopy, which allowed the determination of dissociation constants for a number of protein-ligand complexes. The values obtained were in the nanomolar to micromolar range and correlate positively with the ligand size. These data were supplemented by X-ray crystal structures of the apododecin and holocomplexes with lumichrome, lumiflavin, riboflavin and FMN at resolutions between 1.55 to 1.95 A to unravel a gating mechanism as the structural basis for the preferential binding of the small ligands lumichrome and lumiflavin. The detailed analysis of the dodecin manifold for preferential binding of lumichrome and lumiflavin provides insight on a subatom level into a protein's strategy to gain selectivity for low molecular mass compounds by steric restrictions rather than specific interactions. Investigations on the ligand composition of a wild-type dodecin crystal (1.32 A resolution) support conclusions of functional and structural investigations on heterologously expressed dodecin, and strongly suggest that lumichrome, a molecule associated with the flavin metabolism, is a ligand of dodecin in vivo. Studies on mutant protein and a Halorhodospira halophila homologue spread the idea of a lumichrome binding system as a possible "waste"-trapping device, widely distributed in prokaryotes.     

Pre-docking filter for protein and ligand 3D structures

Virtual drug screening using protein-ligand docking techniques is a time-consuming process, which requires high computational power for binding affinity calculation. There are millions of chemical compounds available for docking. Eliminating compounds that are unlikely to exhibit high binding affinity from the screening set should speed-up the virtual drug screening procedure. We performed docking of 6353 ligands against twenty-one protein X-ray crystal structures. The docked ligands were ranked according to their calculated binding affinities, from which the top five hundred and the bottom five hundred were selected. We found that the volume and number of rotatable bonds of the top five hundred docked ligands are similar to those found in the crystal structures and corresponded with the volume of the binding sites. In contrast, the bottom five hundred set contains ligands that are either too large to enter the binding site, or too small to bind with high specificity and affinity to the binding site. A pre-docking filter that takes into account shapes and volumes of the binding sites as well as ligand volumes and flexibilities can filter out low binding affinity ligands from the screening sets. Thus, the virtual drug screening procedure speed is increased.     

The nature of inhibition of DNA gyrase by the coumarins and the cyclothialidines revealed by X-ray crystallography.

This study describes the first crystal structures of a complex between a DNA topoisomerase and a drug. We present the structures of a 24 kDa N-terminal fragment of the Escherichia coli DNA gyrase B protein in complexes with two different inhibitors of the ATPase activity of DNA gyrase, namely the coumarin antibiotic, novobiocin, and GR122222X, a member of the cyclothialidine family. These structures are compared with the crystal structure of the complex with an ATP analogue, adenylyl-beta-gamma-imidodiphosphate (ADPNP). The likely mechanism, by which mutant gyrase B proteins become resistant to inhibition by novobiocin are discussed in light of these comparisons. The three ligands are quite dissimilar in chemical structure and bind to the protein in very different ways, but their binding is competitive because of a small degree of overlap of their binding sites. These crystal structures consequently describe a chemically well characterized ligand binding surface and provide useful information to assist in the design of novel ligands.     

Crystallographic and calorimetric analysis of peptide binding to OppA protein.

Isothermal titration calorimetry has been used to study the binding of 20 different peptides to the peptide binding protein OppA, and the crystal structures of the ligand complexes have been refined. This periplasmic binding protein, part of the oligopeptide permease system of Gram negative bacteria, has evolved to bind and enclose small peptides of widely varying sequences. The peptides used in this study have the sequence Lys-X-Lys, where X is any of the 20 commonly occurring amino acids. The various side-chains found at position 2 on the ligand fit into a hydrated pocket. The majority of side-chains are restrained to particular conformations within the pocket. Water molecules act as flexible adapters, matching the hydrogen-bonding requirements of the protein and ligand and shielding charges on the buried ligand. This use of water by OppA to broaden the repertoire of its binding site is not unique, but contrasts sharply with other proteins which use water to help bind ligands highly selectively. Predicting the thermodynamics of binding from the structure of the complexes is highly complicated by the influence of water on the system.     

Preference of small molecules for local minimum conformations when binding to proteins

It is well known that small molecules (ligands) do not necessarily adopt their lowest potential energy conformations when binding to proteins. Analyses of protein-bound ligand crystal structures have reportedly shown that many of them do not even adopt the conformations at local minima of their potential energy surfaces (local minimum conformations). The results of these analyses raise a concern regarding the validity of virtual screening methods that use ligands in local minimum conformations. Here we report a normal-mode-analysis (NMA) study of 100 crystal structures of protein-bound ligands. Our data show that the energy minimization of a ligand alone does not automatically stop at a local minimum conformation if the minimum of the potential energy surface is shallow, thus leading to the folding of the ligand. Furthermore, our data show that all 100 ligand conformations in their protein-bound ligand crystal structures are nearly identical to their local minimum conformations obtained from NMA-monitored energy minimization, suggesting that ligands prefer to adopt local minimum conformations when binding to proteins. These results both support virtual screening methods that use ligands in local minimum conformations and caution about possible adverse effect of excessive energy minimization when generating a database of ligand conformations for virtual screening.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

Crystallographic Study of Novel Transthyretin Ligands Exhibiting Negative-Cooperativity between Two Thyroxine Binding Sites

Background

Transthyretin (TTR) is a homotetrameric serum and cerebrospinal fluid protein that transports thyroxine (T4) and retinol by binding to retinol binding protein. Rate-limiting tetramer dissociation and rapid monomer misfolding and disassembly of TTR lead to amyloid fibril formation in different tissues causing various amyloid diseases. Based on the current understanding of the pathogenesis of TTR amyloidosis, it is considered that the inhibition of amyloid fibril formation by stabilization of TTR in native tetrameric form is a viable approach for the treatment of TTR amyloidosis.

Methodology and Principal Findings

We have examined interactions of the wtTTR with a series of compounds containing various substitutions at biphenyl ether skeleton and a novel compound, previously evaluated for binding and inhibiting tetramer dissociation, by x-ray crystallographic approach. High resolution crystal structures of five ligands in complex with wtTTR provided snapshots of negatively cooperative binding of ligands in two T4 binding sites besides characterizing their binding orientations, conformations, and interactions with binding site residues. In all complexes, the ligand has better fit and more potent interactions in first T4 site i.e. (AC site) than the second T4 site (BD site). Together, these results suggest that AC site is a preferred ligand binding site and retention of ordered water molecules between the dimer interfaces further stabilizes the tetramer by bridging a hydrogen bond interaction between Ser117 and its symmetric copy.

Conclusion

Novel biphenyl ether based compounds exhibit negative-cooperativity while binding to two T4 sites which suggests that binding of only single ligand molecule is sufficient to inhibit the TTR tetramer dissociation.     

Interactions of cyclic and non-cyclic naphthalene diimide derivatives with different nucleic acids

Recently, strategy based on stabilization of G-quadruplex telomeric DNA by small organic molecule has been realized by naphthalene diimide derivatives (NDIs). At the same time NDIs bind to DNA duplex as threading intercalators. Here we present cyclic derivative of naphthalene diimide (ligand 1) as DNA-binding ligand with ability to recognition of different structures of telomeric G-quadruplexes and ability to bis-intercalate to double-stranded helixes. The results have been compared to non-cyclic derivative (ligand 2) and revealed that preferential binding of ligands to nucleic acids strongly depends on their topology and structural features of ligands.