The roles of contact residue disorder and domain composition in characterizing protein-ligand binding specificity and promiscuity

Most protein chains interact with only one ligand but a small number of protein chains can bind several ligands, and many examples are available in the protein-ligand complex database of PDB. Among these proteins, some show preferences for the ligands or types of ligands they bind; however, so far we have only poor understanding of what determines protein-ligand binding and its specificity. Here we investigate the structural and functional properties of proteins in protein-ligand complexes. Analysis of the protein-ligand complex dataset from the PDB structure database reveals that proteins with more interactions have more disordered contact residues. Those proteins containing few disordered contact residues that bind multiple ligands have a tendency to consist of several domains. Analysis of physicochemical properties of hub contact residues binding multiple ligands indicates that they are enriched for hydrophilic, charged, polar and His-Asp catalytic triad residues. Finally, in order to differentiate proteins binding different classes of ligands, we mapped the three most prominent classes of ligands onto different superfamily domains. Our results demonstrate that contact residue disorder and ordered multiple domains are complementary factors that play a crucial role in determining ligand binding specificity and promiscuity.     

A new protein binding pocket similarity measure based on comparison of clouds of atoms in 3D: application to ligand prediction

Background  

Predicting which molecules can bind to a given binding site of a protein with known 3D structure is important to decipher the protein function, and useful in drug design. A classical assumption in structural biology is that proteins with similar 3D structures have related molecular functions, and therefore may bind similar ligands. However, proteins that do not display any overall sequence or structure similarity may also bind similar ligands if they contain similar binding sites. Quantitatively assessing the similarity between binding sites may therefore be useful to propose new ligands for a given pocket, based on those known for similar pockets.     

Bile acid interactions with rabbit ileal lipid binding protein and an engineered helixless variant reveal novel ligand binding properties of a versatile beta-clam shell protein scaffold

The intracellular ileal lipid binding proteins (ILBPs) are involved in the transport and enterohepatic circulation of bile acids. ILBPs from different species show high sequence and structural homology and have been shown to bind multiple bile acid ligands with differing degrees of selectivity and positive co-operativity. Human ILBP binds bile acid derivatives in a well-characterised 2:1 ligand:protein complex, however, we show that the highly homologous rabbit ILBP (82% sequence identity) with seven conservative substitutions preferentially binds multiple conjugated deoxycholate ligands in a novel 3:1 binding mode essentially within the same beta-clam shell structure. We have extended these studies to investigate the role of the alpha-helical capping motif (residues 9-35) in controlling the dimensions of the binding cavity and ligand uptake. Substituting the alpha-helical motif (residues 9-35) with a short Gly-Gly-Ser-Gly linker dramatically affects the protein stability such that under physiological conditions the mutant (Deltaalpha-ILBP) is highly disordered. However, we show that the inability of the mutant to adopt a stable three-dimensional structure under these conditions is no barrier to binding ligands with near-native affinity. These structural modifications not only demonstrate the possibility of strong coupling between ligand binding and protein folding, but result in changes in bile acid selectivity and binding stoichiometry, which we characterise in detail using isothermal calorimetry and mass spectrometry.     

Determination of the binding frame within a physiological ligand for the chaperone SecB.

The hallmark of the class of proteins called chaperones is the amazing ability to bind tightly to a wide array of polypeptide ligands that have no consensus in sequence; chaperones recognize non-native structure. As a step in the elucidation of the molecular mechanism of such remarkable binding, we have characterized complexes between the bacterial chaperone SecB and a series of ligands related to maltose-binding protein. SecB interacts at multiple sites on its polypeptide ligand. The entire binding region covers approximately half of the primary sequence of maltose-binding protein and comprises contiguous sites positioned around the center of the sequence.     

Use of the multiple copy simultaneous search (MCSS) method to design a new class of picornavirus capsid binding drugs

A combinatorial ligand design approach based on the multiple copy simultaneous search (MCSS) method and a simple scheme for joining MCSS functional group sites was applied to the binding pocket of P3/Sabin poliovirus and rhinovirus 14. The MCSS method determines where specific functional (chemical) groups have local potential energy minima in the binding site. Before the virus application, test calculations were run to determine the optimal set of input parameters to be used in evaluating the MCSS results. The MCSS minima are analyzed and selected minima are connected with (CH₂) _n linkers to form candidate ligands, whose structures are optimized in the binding site. Estimates of the binding strength were made for the ligands and compared with those for known drugs. The results indicate that the proposed ligands should bind to P3/Sabin poliovirus at least as well as the best of the existing drugs, and that they should also bind to P1/Mahoney poliovirus and rhinovirus 14. A detailed comparison of the poliovirus and rhinovirus binding pockets and an analysis of drug binding specificity is presented. Proteins 29:32–58, 1997. © 1997 Wiley-Liss, Inc.     

Binding specificities and affinities of egf domains for ErbB receptors

ErbB receptor activation is a complex process and is dependent upon the type and number of receptors expressed on a given cell. Previous studies with defined combinations of ErbB receptors expressed in mammalian cells have helped elucidate specific biological responses for many of the recognized gene products that serve as ligands for these receptors. However, no study has examined the binding of these ligands in a defined experimental system. To address this issue, the relative binding affinities of the egf domains of eleven ErbB ligands were measured on six ErbB receptor combinations using a soluble receptor-ligand binding format. The ErbB2/4 heterodimer was shown to bind all ligands tested with moderate to very high affinity. In contrast, ErbB3 showed much more restrictive ligand binding specificity and measurable binding was observed only with heregulin, neuregulin2beta, epiregulin and the synthetic heregulin/egf chimera, biregulin. These studies also revealed that ErbB2 preferentially enhances ligand binding to ErbB3 or ErbB4 and to a lesser degree to ErbB1.     

Peptide interactions with G-protein coupled receptors.

Peptide recognition by G-protein coupled receptors (GPCRs) is reviewed with an emphasis on the indirect approach used to determine the receptor-bound conformation of peptide ligands. This approach was developed in response to the lack of detailed structural information available for these receptors. Recent advances in the structural determination of rhodopsin (the GPCR of the visual system) by crystallography have provided a scaffold for homology modeling of the inactive state of a wide variety of GPCRs that interact with peptide messages. Additionally, the ability to mutate GPCRs and assay compounds of similar chemical structure to test a common binding site on the receptor provides a firm experimental basis for structure-activity studies. Recognition motifs, common in other well-studied systems such as proteolytic enzymes and major histocompatibility class receptors (MHC) are reviewed briefly to provide a basis of comparison. Finally, the development of true peptidomimetics is contrasted with nonpeptide ligands, discovered through combinatorial chemistry. In many systems, the evidence suggests that the peptide ligands bind at the interface between the transmembrane segments and the extracellular loops, while nonpeptide antagonists bind within the transmembrane segments. Plausible models of GPCRs and the mechanism by which they activate G-proteins on binding peptides are beginning to emerge.     

Amyloidogenic Propensity of a Natural Variant of Human Apolipoprotein A-I: Stability and Interaction with Ligands

A number of naturally occurring mutations of human apolipoprotein A-I (apoA-I) have been associated with hereditary amyloidoses. The molecular mechanisms involved in amyloid-associated pathology remain largely unknown. Here we examined the effects of the Arg173Pro point mutation in apoA-I on the structure, stability, and aggregation propensity, as well as on the ability to bind to putative ligands. Our results indicate that the mutation induces a drastic loss of stability, and a lower efficiency to bind to phospholipid vesicles at physiological pH, which could determine the observed higher tendency to aggregate as pro-amyloidogenic complexes. Incubation under acidic conditions does not seem to induce significant desestabilization or aggregation tendency, neither does it contribute to the binding of the mutant to sodium dodecyl sulfate. While the binding to this detergent is higher for the mutant as compared to wt apoA-I, the interaction of the Arg173Pro variant with heparin depends on pH, being lower at pH 5.0 and higher than wt under physiological pH conditions. We suggest that binding to ligands as heparin or other glycosaminoglycans could be key events tuning the fine details of the interaction of apoA-I variants with the micro-environment, and probably eliciting the toxicity of these variants in hereditary amyloidoses.     

The Structure,Stability and Pheromone Binding of the Male Mouse Protein Sex Pheromone Darcin

Mouse urine contains highly polymorphic major urinary proteins that have multiple functions in scent communication through their abilities to bind, transport and release hydrophobic volatile pheromones. The mouse genome encodes for about 20 of these proteins and are classified, based on amino acid sequence similarity and tissue expression patterns, as either central or peripheral major urinary proteins. Darcin is a male specific peripheral major urinary protein and is distinctive in its role in inherent female attraction. A comparison of the structure and biophysical properties of darcin with MUP11, which belongs to the central class, highlights similarity in the overall structure between the two proteins. The thermodynamic stability, however, differs between the two proteins, with darcin being much more stable. Furthermore, the affinity of a small pheromone mimetic is higher for darcin, although darcin is more discriminatory, being unable to bind bulkier ligands. These attributes are due to the hydrophobic ligand binding cavity of darcin being smaller, caused by the presence of larger amino acid side chains. Thus, the physical and chemical characteristics of the binding cavity, together with its extreme stability, are consistent with darcin being able to exert its function after release into the environment.     

Organization of ligand binding sites at the acetylcholine receptor: a study with monoclonal antibodies

We have studied 20 monoclonal antibodies directed against both the solubilized and the membrane-bound receptor from Torpedo marmorata. We find the following: (i) Six of the antibodies compete with cholinergic ligands for receptor binding and, hence, are directed against the ligand binding regions. (ii) Of these six antibodies, two cross-react with receptor from Electrophorus electricus, rat myotubes, and chicken sympathetic ganglia. These two antibodies therefore define a preserved structure within the ligand binding regions. The other four antibodies bind to structures not common between the receptor preparations tested. (iii) From competition binding studies using internally 3H-labeled antibodies, nine nonoverlapping antigenic regions were defined at the surface of the receptor. Three of these regions overlap with the ligand binding regions. Since two of these three regions do not overlap with each other, two structurally distinct ligand binding regions must exist at the receptor. (iv) From competition binding studies with representative cholinergic ligands, the antibodies directed against the ligand binding regions can be subdivided into three groups: one group competes with all ligands tested; the second group competes with all ligands except the bismethonium compounds; the third group competes with all ligands except the bismethonium compounds and tubocurarine. The results are summarized in a model of the organization of ligand binding sites at the receptor: There are two ligand binding regions differing in their antigenic properties. Furthermore, either there exists separate sites for distinct groups of ligands within each of these binding regions or some ligands produce conformational changes of the receptor that reversibly abolish some antigenic sites. In any case, the cholinergic ligands must interact with the receptor by more and/or other structural determinants than are provided by the structure of acetylcholine.     

Structural basis for the binding of compatible solutes by ProX from the hyperthermophilic archaeon Archaeoglobus fulgidus

Compatible solutes such as glycine betaine and proline betaine serve as protein stabilizers because of their preferential exclusion from protein surfaces. To use extracellular sources of this class of compounds as osmo-, cryo-, or thermoprotectants, Bacteria and Archaea have developed high affinity uptake systems of the ATP-binding cassette type. These transport systems require periplasmic- or extracellular-binding proteins that are able to bind the transported substance with high affinity. Therefore, binding proteins that bind compatible solutes have to avoid the exclusion of their ligands within the binding pocket. In the present study we addressed the question to how compatible solutes can be effectively bound by a protein at temperatures around 83 degrees C as this is done by the ligand-binding protein ProX from the hyperthermophilic archaeon Archaeoglobus fulgidus. We solved the structures of ProX without ligand and in complex with both of its natural ligands glycine betaine and proline betaine, as well as in complex with the artificial ligand trimethylammonium. Cation-pi interactions and non-classical hydrogen bonds between four tyrosine residues, a main chain carbonyl oxygen, and the ligand have been identified to be the key determinants in binding the quaternary amines of the three investigated ligands. The comparison of the ligand binding sites of ProX from A. fulgidus and the recently solved structure of ProX from Escherichia coli revealed a very similar solution for the problem of compatible solute binding, although both proteins share only a low degree of sequence identity. The residues involved in ligand binding are functionally equivalent but not conserved in the primary sequence.     

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca2? dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 ?, 5.5 ?, 6 ?, 6.5 ?, 7.0 ? and 7.5 ? were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.     

The dual-function disabled-1 PTB domain exhibits site independence in binding phosphoinositide and peptide ligands

While typical intracellular protein modules have only one ligand-binding site, there are rare examples of single modules that bind two different ligands at distinct binding sites. Here we present a detailed mutational and energetic analysis of one such domain, the phosphotyrosine binding (PTB) domain of Disabled-1 (Dab1), which binds to both peptide and phosphoinositide (PI) ligands simultaneously at structurally distinct binding sites. Through the techniques of isothermal titration calorimetry (ITC), analysis of Dab1 PTB domain mutants, and nuclear magnetic resonance (NMR), we have evaluated the characteristics of binding of the Dab1 PTB domain to various peptide and PI ligands. These studies reveal that the presence of saturating concentrations of one ligand has little effect on the binding constant for a second ligand at the other site. In addition, proteins with single-point mutations in the peptide-binding site retain native affinity for PI ligands, while proteins with mutations that prevent PI binding retain native affinity for peptide. NMR titrations show that the final structure of the ternary complex is the same independent of the order of addition of the two ligands. Together, these studies show that binding of peptide and PI ligands is energetically independent and noncooperative.     

HLA class I allele promiscuity revisited

The peptide repertoire presented on human leukocyte antigen (HLA) class I molecules is largely determined by the structure of the peptide binding groove. It is expected that the molecules having similar grooves (i.e., belonging to the same supertype) might present similar/overlapping peptides. However, the extent of promiscuity among HLA class I ligands remains controversial: while in many studies T cell responses are detected against epitopes presented by alternative molecules across HLA class I supertypes and loci, peptide elution studies report minute overlaps between the peptide repertoires of even related HLA molecules. To get more insight into the promiscuous peptide binding by HLA molecules, we analyzed the HLA peptide binding data from the large epitope repository, Immune Epitope Database (IEDB), and further performed in silico analysis to estimate the promiscuity at the population level. Both analyses suggest that an unexpectedly large fraction of HLA ligands (>50%) bind two or more HLA molecules, often across supertype or even loci. These results suggest that different HLA class I molecules can nevertheless present largely overlapping peptide sets, and that “functional” HLA polymorphism on individual and population level is probably much lower than previously anticipated.     

Recognition of RNA duplex by a neomycin–Hoechst 33258 conjugate

DNA minor groove binding drugs such as Hoechst 33258 have been shown to bind to a number of RNA structures. Similarly, RNA binding ligands such as neomycin have been shown by us to bind to a number of A-form DNA structures. A neomycin–Hoechst 33258 conjugate was recently shown to bind B-DNA, where Hoechst exhibits high affinity for the minor groove of A/T tract DNA and neomycin docks into the major groove. Further studies now indicate that the Hoechst moiety of the conjugate can be driven to bind RNA duplex as a consequence of neomycin binding in the RNA major groove. This is the first example of Hoechst 33258 binding to RNA duplex not containing bulges or loop motifs.     

Fluorophore labeled kinase detects ligands that bind within the MAPK insert of p38α kinase

The vast majority of small molecules known to modulate kinase activity, target the highly conserved ATP-pocket. Consequently, such ligands are often less specific and in case of inhibitors, this leads to the inhibition of multiple kinases. Thus, selective modulation of kinase function remains a major hurdle. One of the next great challenges in kinase research is the identification of ligands which bind to less conserved sites and target the non-catalytic functions of protein kinases. However, approaches that allow for the unambiguous identification of molecules that bind to these less conserved sites are few in number. We have previously reported the use of fluorescent labels in kinases (FLiK) to develop direct kinase binding assays that exclusively detect ligands which stabilize inactive (DFG-out) kinase conformations. Here, we present the successful application of the FLiK approach to develop a high-throughput binding assay capable of directly monitoring ligand binding to a remote site within the MAPK insert of p38α mitogen-activated protein kinase (MAPK). Guided by the crystal structure of an initially identified hit molecule in complex with p38α, we developed a tight binding ligand which may serve as an ideal starting point for further investigations of the biological function of the MAPK insert in regulating the p38α signaling pathway.     

Joint interaction of ethidium bromide and methylene blue with DNA. The effect of ionic strength on binding thermodynamic parameters

Large amount of data of experimental and theoretical studies have shown that ethidium bromide (EtBr) and methylene blue (MB) may bind to nucleic acids via three modes: intercalation between two adjacent base pairs, insertion into the plane between neighboring bases in the same strand (semi-intercalation), and outside binding with negatively charged backbone phosphate groups. The aim of the given research is to examine the behavior of these two ligands at both separate and joint DNA binding. The obtained experimental data show that the effect of simultaneous binding of EtBr and MB on double-stranded DNA has a non-additive effect of separate binding. The analyses of the melting thermodynamic parameters of DNA complexes with two bound ligands suggest competitive mechanism of interaction.     

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.