Putative Ligand Binding Sites Research Articles

Computational Resources for Protein Modelling and Drug Discovery Applications

The design of new medications is an intensive, time-consuming and costly process. Over the years, a rational approach that exploits the structural knowledge of a biological target has led to many successes. This procedure can be expedited using computer-aided modelling techniques. The structure-based approach to drug design relies on knowing the three-dimensional structure of the target macromolecule. If an experimental structure has not been determined yet, a good approximation of the protein target structure can be obtained through computational modelling, provided that some structures of its homologues are available to serve as templates. The vast majority of drugs currently on the market act by disrupting the interaction between a protein and its physiological ligand(s). Hence, once a molecular model is available, the next step is to identify and study its putative ligand-binding sites. Molecular "docking" may then be performed in silico to predict the modes of interaction between the ligand and the target. In this review, a list of computational resources for structure-based drug design has been compiled. It is hoped that readers who do not have much experiences will be equipped with the appropriate tools to make a first attempt at protein modelling and in silico ligand docking exercises.

Thermostabilization of the Neurotensin Receptor NTS1

Structural studies on G-protein-coupled receptors have been hampered for many years by their instability in detergent solution and by the number of potential conformations that receptors can adopt. Recently, the structures of the β1 and β2 adrenergic receptors and the adenosine A2a receptor were determined in the antagonist-bound state, a receptor conformation that is thought to be more stable than the agonist-bound state. In contrast to these receptors, the neurotensin (NT) receptor NTS1 is much less stable in detergent solution. We have therefore used a systematic mutational approach coupled with activity assays to identify receptor mutants suitable for crystallization, both alone and in complex with the peptide agonist NT. The best receptor mutant NTS1-7m contained four point mutations. It showed increased stability compared to the wild-type receptor, in the absence of ligand, after solubilization with a variety of detergents. In addition, NTS1-7m bound to NT was more stable than unliganded NTS1-7m. Of the four thermostabilizing mutations, only one residue (A86L) is predicted to be in the lipid environment. In contrast, I260A appears to be buried within the transmembrane helix bundle, F342A may form a distant part of the putative ligand-binding site, whereas F358A is likely to be in a region that is important for receptor activation. NTS1-7m binds NT with a similar affinity for the wild-type receptor. However, agonist dissociation was slower, and NTS1-7m activated G-proteins poorly. The affinity of NTS1-7m for the antagonist SR48692 was also lower than that of the wild-type receptor. Thus, we have successfully stabilized NTS1 in an agonist-binding conformation that does not efficiently couple to G-proteins.

Homology Model Adjustment and Ligand Screening with a Pseudoreceptor of the Human Histamine H4 Receptor

A computer-assisted method for the generation of pseudoreceptor models is presented together with two practical applications. From a three-dimensional alignment of known histamine H(4) receptor ligands, a pseudoreceptor model of the putative ligand binding site was constructed and used for virtual screening of a large collection of commercially available compounds. Two bioactive chemotypes were retrieved, demonstrating the general applicability of the approach. The pseudoreceptor model was also used to find the putative ligand binding pocket within the transmembrane receptor domain. For each frame of a molecular dynamics simulation of a homology-based H(4) receptor model, we automatically extracted potential ligand binding pockets and used their compatibility with the pseudoreceptor as a selection criterion. The best-matching pocket fits perfectly with existing mutation data and previously published hypotheses suggesting Glu182(5.46) as the preferred binding partner of a positively charged moiety of H(4) receptor ligands. This new pseudoreceptor approach has demonstrated its suitability for both structure-based prioritization of protein receptor models, and ligand-based virtual screening with the aim to perform scaffold hopping.

SITEHOUND-web: a server for ligand binding site identification in protein structures

SITEHOUND-web (http://sitehound.sanchezlab.org) is a binding-site identification server powered by the SITEHOUND program. Given a protein structure in PDB format SITEHOUND-web will identify regions of the protein characterized by favorable interactions with a probe molecule. These regions correspond to putative ligand binding sites. Depending on the probe used in the calculation, sites with preference for different ligands will be identified. Currently, a carbon probe for identification of binding sites for drug-like molecules, and a phosphate probe for phosphorylated ligands (ATP, phoshopeptides, etc.) have been implemented. SITEHOUND-web will display the results in HTML pages including an interactive 3D representation of the protein structure and the putative sites using the Jmol java applet. Various downloadable data files are also provided for offline data analysis.

Gating Current Measurements Reveal Ligand-Selective Conformational Changes in the M2 Muscarinic Receptor

In addition to external stimuli, membrane voltage modulates the activity of several G-protein-coupled receptors (GPCRs). For example, membrane depolarization reduces the affinity of the M2 muscarinic receptor (M2R) for acetylcholine (ACh). The intrinsic capacity of muscarinic receptors to “sense” transmembrane potential was confirmed by the recording of “gating” currents that reflect reorientation of charges within the receptor in response to changes in voltage. We studied the effects of membrane voltage on agonist activation of M2R in atrial myocytes and on gating currents in oocytes heterologously expressing M2R. Depolarization decreased the apparent affinity of M2R for ACh, but increased the apparent affinity and efficacy of pilocarpine (PILO), as measured by acetylcholine-activated K+ current, IKACh. Voltage-induced conformational changes in M2R (as measured by gating currents) were modified in a ligand-selective manner. For example, ACh reduced gating charge displacement while PILO increased the amount of charge displaced. Thus, these ligands manifest opposite voltage-dependent IKACh modulation and exert opposite effects on M2R gating charge displacement. Finally, mutations in the putative ligand binding site perturbed the movement of the M2R voltage sensor and confirmed a link between voltage sensing and changes in agonist affinity. Our data suggest that voltage modulates the conformation of the ligand binding site, while ligand binding induces ligand-specific conformational changes in the receptor. Voltage-dependent GPCR modulation has important implications for cellular signaling in excitable tissues. Gating current measurement allows for the tracking of subtle conformational changes in the receptor that accompany ligand binding and changes in membrane voltage.

Structures of Two Major Allergens, Bla g 4 and Per a 4, from Cockroaches and Their IgE Binding Epitopes

Inhalant allergens from cockroaches are an important cause of asthma to millions of individuals worldwide. Here we report for the first time the structures of two major cockroach allergens, Bla g 4 and Per a 4, that adopt a typical lipocalin fold but with distinct structural features as compared with other known lipocalin allergens. Both Bla g 4 and Per a 4 contain two long-range disulfide bonds linking the N and C termini to a beta-barrel. The C-terminal helix of Bla g 4 is bent and greatly extended toward the N terminus. Bla g 4 is found to be a monomer, whereas Per a 4 exists as a dimer in solution with a novel dimeric interface involving residues from loops at the top and bottom of the beta-barrel. Putative ligand binding sites of both allergens are determined by docking of the juvenile hormone III inside the beta-barrel and found to interact with the ligand using non-conserved residues. Bla g 4 and Per a 4 are found to be cross-reactive in sera IgE binding, at least in the Singaporean Chinese population tested. A major IgE binding epitope unique to Per a 4 is found on the loops at the bottom of the beta-barrel that may aid the development of hypoallergens for immunotherapy.

Structural insights into phenylethanolamines high-affinity binding site in NR2B from binding and molecular modeling studies

BackgroundPhenylethanolamines selectively bind to NR2B subunit-containing N-methyl-D-aspartate-subtype of ionotropic glutamate receptors and negatively modulate receptor activity. To investigate the structural and functional properties of the ifenprodil binding domain on the NR2B protein, we have purified a soluble recombinant rat NR2B protein fragment comprising the first ~400 amino acid amino-terminal domain (ATD2B) expressed in E. coli. Spectral measurements on refolded ATD2B protein demonstrated specific binding to ifenprodil. We have used site-directed mutagenesis, circular dichroism spectroscopy and molecular modeling to obtain structural information on the interactions between critical amino acid residues and ifenprodil of our soluble refolded ATD2B proteins. Ligand-induced changes in protein structure were inferred from changes in the circular dichroism spectrum, and the concentration dependence of these changes was used to determine binding constants for ifenprodil and its analogues.ResultsLigand binding of ifenprodil, RO25,6981 and haloperidol on soluble recombinant ATD2B determined from circular dichroism spectroscopy yielded low-to-high micromolar equilibrium constants which concurred with functional IC50 measurement determined in heterologously expressed NR1/NR2B receptors in Xenopus oocytes. Amino acid residue substitutions of Asp101, Ile150 and Phe176 with alanine residue within the ATD2B protein altered the recombinant protein dissociation constants for ifenprodil, mirroring the pattern of their functional phenotypes. Molecular modeling of ATD2B as a clam-shell-like structure places these critical residues near a putative ligand binding site.ConclusionWe report for the first time biochemical measurements show that the functional measurements actually reflect binding to the ATD of NR2B subunit. Insights gained from this study help advance the theory that ifenprodil is a ligand for the ATD of NR2B subunit.

What induces pocket openings on protein surface patches involved in protein–protein interactions?

We previously showed for the proteins BCL-X(L), IL-2, and MDM2 that transient pockets at their protein-protein binding interfaces can be identified by applying the PASS algorithm to molecular dynamics (MD) snapshots. We now investigated which aspects of the natural conformational dynamics of proteins induce the formation of such pockets. The pocket detection protocol was applied to three different conformational ensembles for the same proteins that were extracted from MD simulations of the inhibitor bound crystal conformation in water and the free crystal/NMR structure in water and in methanol. Additional MD simulations studied the impact of backbone mobility. The more efficient CONCOORD or normal mode analysis (NMA) techniques gave significantly smaller pockets than MD simulations, whereas tCONCOORD generated pockets comparable to those observed in MD simulations for two of the three systems. Our findings emphasize the influence of solvent polarity and backbone rearrangements on the formation of pockets on protein surfaces and should be helpful in future generation of transient pockets as putative ligand binding sites at protein-protein interfaces.

The Structure and Binding Behavior of the Bacterial Cell Surface Layer Protein SbsC

Surface layers (S-layers) comprise the outermost cell envelope component of most archaea and many bacteria. Here we present the structure of the bacterial S-layer protein SbsC from Geobacillus stearothermophilus, showing a very elongated and flexible molecule, with strong and specific binding to the secondary cell wall polymer (SCWP). The crystal structure of rSbsC((31-844)) revealed a novel fold, consisting of six separate domains, which are connected by short flexible linkers. The N-terminal domain exhibits positively charged residues regularly spaced along the putative ligand binding site matching the distance of the negative charges on the extended SCWP. Upon SCWP binding, a considerable stabilization of the N-terminal domain occurs. These findings provide insight into the processes of S-layer attachment to the underlying cell wall and self-assembly, and also accommodate the observed mechanical strength, the polarity of the S-layer, and the pronounced requirement for surface flexibility inherent to cell growth and division.

Structure and ligand binding of the soluble domain of a Thermotoga maritima membrane protein of unknown function TM1634

As a part of the Joint Center for Structural Genomics (JCSG) biological targets, the structures of soluble domains of membrane proteins from Thermotoga maritima were pursued. Here, we report the crystal structure of the soluble domain of TM1634, a putative membrane protein of 128 residues (15.1 kDa) and unknown function. The soluble domain of TM1634 is an alpha-helical dimer that contains a single tetratrico peptide repeat (TPR) motif in each monomer where each motif is similar to that found in Tom20. The overall fold, however, is unique and a DALI search does not identify similar folds beyond the 38-residue TPR motif. Two different putative ligand binding sites, in which PEG200 and Co(2+) were located, were identified using crystallography and NMR, respectively.

Crystal Structure of the MACPF Domain of Human Complement Protein C8α in Complex with the C8γ Subunit

Human C8 is one of five complement components (C5b, C6, C7, C8, and C9) that assemble on bacterial membranes to form a porelike structure referred to as the “membrane attack complex” (MAC). C8 contains three genetically distinct subunits (C8α, C8β, C8γ) arranged as a disulfide-linked C8α–γ dimer that is noncovalently associated with C8β. C6, C7 C8α, C8β, and C9 are homologous. All contain N- and C-terminal modules and an intervening 40-kDa segment referred to as the membrane attack complex/perforin (MACPF) domain. The C8γ subunit is unrelated and belongs to the lipocalin family of proteins that display a β-barrel fold and generally bind small, hydrophobic ligands. Several hundred proteins with MACPF domains have been identified based on sequence similarity; however, the structure and function of most are unknown. Crystal structures of the secreted bacterial protein Plu-MACPF and the human C8α MACPF domain were recently reported and both display a fold similar to those of the bacterial pore-forming cholesterol-dependent cytolysins (CDCs). In the present study, we determined the crystal structure of the human C8α MACPF domain disulfide-linked to C8γ (αMACPF-γ) at 2.15 Å resolution. The αMACPF portion has the predicted CDC-like fold and shows two regions of interaction with C8γ. One is in a previously characterized 19-residue insertion (indel) in C8α and fills the entrance to the putative C8γ ligand-binding site. The second is a hydrophobic pocket that makes contact with residues on the side of the C8γ β-barrel. The latter interaction induces conformational changes in αMACPF that are likely important for C8 function. Also observed is structural conservation of the MACPF signature motif Y/W-G-T/S-H-F/Y-X 6-G-G in αMACPF and Plu-MACPF, and conservation of several key glycine residues known to be important for refolding and pore formation by CDCs.

Crystal structure of the PH1932 protein, a unique archaeal ArsR type winged‐HTH transcription factor from Pyrococcus horikoshii OT3

Crystal structure of the PH1932 protein, a unique archaeal ArsR type winged‐HTH transcription factor from <i>Pyrococcus horikoshii</i> OT3

Crystal structure of complement protein C8γ in complex with a peptide containing the C8γ binding site on C8α: Implications for C8γ ligand binding

Human C8 is one of five complement components (C5b, C6, C7, C8 and C9) that interact to form the cytolytic membrane attack complex. It contains three genetically distinct subunits; C8α and C8γ form a disulfide-linked C8α–γ heterodimer that is noncovalently associated with C8β. The C8α subunit is homologous to C8β, C6, C7 and C9 and together they form the MAC family of proteins. By contrast, C8γ is the only lipocalin in the complement system. Like other lipocalins, it has a core β-barrel structure forming a calyx with a distinct binding pocket for a small and as yet unidentified ligand. The binding site on C8α for C8γ was previously localized to a 19-residue segment which contains an insertion (indel) that is unique to C8α. Included in the indel is C8α Cys 164 which links to Cys 40 in C8γ. In the present study, C8γ containing a C40A substitution was co-crystallized with a synthetic indel peptide containing the equivalent of a C8α C164A substitution. The X-ray crystal structure shows that the indel peptide completely fills the upper portion of the putative C8γ ligand binding pocket and is in contact with all four loops at the calyx entrance. The lower part of the C8γ cavity is either unoccupied or contains disordered solvent. The validity of the structure is supported by the spatial arrangement of C8α Ala 164 in the peptide and C8γ Ala 40, which are within disulfide-bonding distance of each other. Corresponding studies in solution indicate the C8γ ligand binding site is also occupied by the indel segment of C8α in whole C8. These results suggest a role for C8α in regulating access to the putative C8γ ligand binding site.

Detection of pockets on protein surfaces using small and large probe spheres to find putative ligand binding sites

One of the simplest ways to predict ligand binding sites is to identify pocket-shaped regions on the protein surface. Many programs have already been proposed to identify these pocket regions. Examination of their algorithms revealed that a pocket intrinsically has two arbitrary properties, "size" and "depth". We proposed a new definition for pockets using two explicit adjustable parameters that correspond to these two arbitrary properties. A pocket region is defined as a space into which a small probe can enter, but a large probe cannot. The radii of small and large probe spheres are the two parameters that correspond to the "size" and "depth" of the pockets, respectively. These values can be adjusted individual putative ligand molecule. To determine the optimal value of the large probe spheres radius, we generated pockets for thousands of protein structures in the database, using several size of large probe spheres, examined the correspondence of these pockets with known binding site positions. A new measure of shallowness, a minimum inaccessible radius, R(inaccess), indicated that binding sites of coenzymes are very deep, while those for adenine/guanine mononucleotide have only medium shallowness and those for short peptides and oligosaccharides are shallow. The optimal radius of large probe spheres was 3-4 A for the coenzymes, 4 A for adenine/guanine mononucleotides, and 5 A or more for peptides/oligosaccharides. Comparison of our program with two other popular pocket-finding programs showed that our program had a higher performance of detecting binding pockets, although it required more computational time.

Action of Nucleosides and Nucleotides at 7 Transmembrane-Spanning Receptors

Ribose ring-constrained nucleosides and nucleotides to act at cell-surface purine recesptors have been designed and synthesized. At the P2Y1 nucleotide receptor and the A3 adenosine receptor (AR) the North envelope conformation of ribose is highly preferred. We have applied mutagenesis and rhodopsin-based homology modeling to the study of purine receptors and used the structural insights gained to assist in the design of novel ligands. Two subgroups of P2Y receptors have been defined, containing different sets of cationic residues for coordinating the phosphate groups. Modeling/mutagenesis of adenosine receptors has focused on determinants of intrinsic efficacy in adenosine derivatives and on a conserved Trp residue (6.48) which is involved in the activation process. The clinical use of adenosine agonists as cytoprotective agents has been limited by the widespread occurrence of ARs, thus, leading to undesirable side effects of exogenously administered adenosine derivatives. In order to overcome the inherent nonselectivity of activating the native receptors, we have introduced the concept of neoceptors. By this strategy, intended for eventual use in gene therapy, the putative ligand binding site of a G protein-coupled receptor is reengineered for activation by synthetic agonists (neoligands) built to have a structural complementarity. Using a rational design process we have identified neoceptor-neoligand pairs which are pharmacologically orthogonal with respect to the native species.

The transthyretin-related protein: Structural investigation of a novel protein family

The transthyretin-related protein (TRP) family comprises proteins predicted to be structurally related to the homotetrameric transport protein transthyretin (TTR). The function of TRPs is not yet fully established, but recent data suggest that they are involved in purine catabolism. We have determined the three-dimensional structure of the Escherichia coli TRP in two crystal forms; one at 1.65 Å resolution in the presence of zinc, and the other at 2.1 Å resolution in the presence of zinc and bromide. The structures revealed five zinc-ion-binding sites per monomer. Of these, the zinc ions bound at sites I and II are coordinated in tetrahedral geometries to the side chains of residues His9, His96, His98, Ser114, and three water molecules at the putative ligand-binding site. Of these four residues, His9, His98, and Ser114 are conserved. His9 and His98 bind the central zinc (site I) together with two water molecules. The side chain of His98 also binds to the zinc ion at site II. Bromide ions bind at site I only, replacing one of the water molecules coordinated to the zinc ion. The C-terminal four amino acid sequence motif Y-[RK]-G-[ST] constitutes the signature sequence of the TRP family. Two Tyr111 residues form direct hydrogen bonds to each other over the tetramer interface at the area, which in TTR constitutes the rear part of its thyroxine-binding channel. The putative substrate/ligand-binding channel of TRP is consequently shallower and broader than its counterpart in TTR.

Structures of the Carbohydrate Recognition Domain of Ca2+-independent Cargo Receptors Emp46p and Emp47p

Emp46p and Emp47p are type I membrane proteins, which cycle between the endoplasmic reticulum (ER) and the Golgi apparatus by vesicles coated with coat protein complexes I and II (COPI and COPII). They are considered to function as cargo receptors for exporting N-linked glycoproteins from the ER. We have determined crystal structures of the carbohydrate recognition domains (CRDs) of Emp46p and Emp47p of Saccharomyces cerevisiae, in the absence and presence of metal ions. Both proteins fold as a beta-sandwich, and resemble that of the mammalian ortholog, p58/ERGIC-53. However, the nature of metal binding is distinct from that of Ca(2+)-dependent p58/ERGIC-53. Interestingly, the CRD of Emp46p does not bind Ca(2+) ion but instead binds K(+) ion at the edge of a concave beta-sheet whose position is distinct from the corresponding site of the Ca(2+) ion in p58/ERGIC-53. Binding of K(+) ion to Emp46p appears essential for transport of a subset of glycoproteins because the Y131F mutant of Emp46p, which cannot bind K(+) ion fails to rescue the transport in disruptants of EMP46 and EMP47 genes. In contrast the CRD of Emp47p binds no metal ions at all. Furthermore, the CRD of Emp46p binds to glycoproteins carrying high mannosetype glycans and the is promoted by binding not the addition of Ca(2+) or K(+) ion in These results suggest that Emp46p can be regarded as a Ca(2+)-independent intracellular lectin at the ER exit sites.

Potencies of Human Immunodeficiency Virus Protease Inhibitors In Vitro againstPlasmodium falciparumand In Vivo against Murine Malaria

Parasite resistance to antimalarial drugs is a serious threat to human health, and novel agents that act on enzymes essential for parasite metabolism, such as proteases, are attractive targets for drug development. Recent studies have shown that clinically utilized human immunodeficiency virus (HIV) protease inhibitors can inhibit the in vitro growth of Plasmodium falciparum at or below concentrations found in human plasma after oral drug administration. The most potent in vitro antimalarial effects have been obtained for parasites treated with saquinavir, ritonavir, or lopinavir, findings confirmed in this study for a genetically distinct P. falciparum line (3D7). To investigate the potential in vivo activity of antiretroviral protease inhibitors (ARPIs) against malaria, we examined the effect of ARPI combinations in a murine model of malaria. In mice infected with Plasmodium chabaudi AS and treated orally with ritonavir-saquinavir or ritonavir-lopinavir, a delay in patency and a significant attenuation of parasitemia were observed. Using modeling and ligand docking studies we examined putative ligand binding sites of ARPIs in aspartyl proteases of P. falciparum (plasmepsins II and IV) and P. chabaudi (plasmepsin) and found that these in silico analyses support the antimalarial activity hypothesized to be mediated through inhibition of these enzymes. In addition, in vitro enzyme assays demonstrated that P. falciparum plasmepsins II and IV are both inhibited by the ARPIs saquinavir, ritonavir, and lopinavir. The combined results suggest that ARPIs have useful antimalarial activity that may be especially relevant in geographical regions where HIV and P. falciparum infections are both endemic.

Grand Canonical Monte Carlo Simulation of Ligand−Protein Binding

A new application of the grand canonical thermodynamics ensemble to compute ligand-protein binding is described. The described method is sufficiently rapid that it is practical to compute ligand-protein binding free energies for a large number of poses over the entire protein surface, thus identifying multiple putative ligand binding sites. In addition, the method computes binding free energies for a large number of poses. The method is demonstrated by the simulation of two protein-ligand systems, thermolysin and T4 lysozyme, for which there is extensive thermodynamic and crystallographic data for the binding of small, rigid ligands. These low-molecular-weight ligands correspond to the molecular fragments used in computational fragment-based drug design. The simulations correctly identified the experimental binding poses and rank ordered the affinities of ligands in each of these systems.

Identification of Determinants of Ligand Binding Affinity and Selectivity in the Prostaglandin D2 Receptor CRTH2

The chemoattractant receptor-homologous molecule expressed on Th2 cells (CRTH2) is a G protein-coupled receptor that mediates the pro-inflammatory effects of prostaglandin D(2) (PGD(2)) generated in allergic inflammation. The CRTH2 receptor shares greatest sequence similarity with chemoattractant receptors compared with prostanoid receptors. To investigate the structural determinants of CRTH2 ligand binding, we performed site-directed mutagenesis of putative mCRTH2 ligand-binding residues, and we evaluated mutant receptor ligand binding and functional properties. Substitution of alanine at each of three residues in the transmembrane (TM) helical domains (His-106, TM III; Lys-209, TM V; and Glu-268, TM VI) and one in extracellular loop II (Arg-178) decreased PGD(2) binding affinity, suggesting that these residues play a role in binding PGD(2). In contrast, the H106A and E268A mutants bound indomethacin, a nonsteroidal anti-inflammatory drug, with an affinity similar to the wild-type receptor. HEK293 cells expressing the H106A, K209A, and E268A mutants displayed reduced inhibition of intracellular cAMP and chemotaxis in response to PGD(2), whereas the H106A and E268A mutants had functional responses to indomethacin similar to the wild-type receptor. Binding of PGE(2) by the E268A mutant was enhanced compared with the wild-type receptor, suggesting that Glu-268 plays a role in determining prostanoid ligand selectivity. Replacement of Tyr-261 with phenylalanine did not affect PGD(2) binding but decreased the binding affinity for indomethacin. These results provided the first details of the ligand binding pocket of an eicosanoid-binding chemoattractant receptor.