Fatty-acid-binding protein from the flight muscle of Locusta migratoria: evolutionary variations in fatty acid binding

Intracellular lipid-binding proteins have evolved from a common ancestral gene with the appearance of mitochondrial oxidation, to guarantee, for example, transport of fatty acids through the aqueous cytosol to their site of utilization. The mammalian forms of these lipid carriers are structurally well-characterized and have been categorized, on the basis of sequence similarities and several typical ligand-binding features, into four subfamilies. Only a single complex structure of an invertebrate fatty-acid-binding protein (FABP) has been reported to date, which reveals a unique ligand-binding arrangement yet unknown in vertebrate FABPs. In the present study, the structure of a second invertebrate FABP (locust muscle) complexed with a fatty acid has been determined on the basis of intermolecular NOE connectivities between the protein and the uniformly (13)C-enriched oleate ligand. The resulting ligand conformation, although resembling the closely related mammalian heart- and adipocyte-type FABPs, is characterized by certain binding features that differ significantly from the typical hairpin-turn ligand shapes of the latter forms. This is primarily due to an alanine-to-leucine substitution in locust FABPs that produces a steric hindrance for ligand binding. A comparison with an FABP from tobacco hornworm larvae furthermore demonstrates that certain amino acid substitutions that appear to be specific for invertebrates decidedly influence the binding arrangement inside the protein cavity. Hence, as a result of these evolutionary variations, invertebrate FABPs may display a much greater diversity in intracellular lipid binding than observed for the mammalian transport proteins, thus possibly providing new insights for the design of modified lipid carriers.     

Computational detection of the binding-site hot spot at the remodeled human growth hormone-receptor interface

A hierarchical computational approach is used to identify the engineered binding-site cavity at the remodeled intermolecular interface between the mutants of human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp). Multiple docking simulations are conducted with the remodeled hGH-hGHbp complex for a panel of potent benzimidazole-containing inhibitors that can restore the binding affinity of the wild-type complex, and for a set of known nonactive small molecules that contain different heterocyclic motifs. Structural clustering of ligand-bound conformations and binding free-energy calculations, using the AMBER force field and a continuum solvation model, can rapidly locate and screen numerous ligand-binding modes on the protein surface and detect the binding-site hot spot at the intermolecular interface. Structural orientation of the benzimidazole motif in the binding-site cavity closely mimics the position of the hot spot residue W104 in the crystal structure of the wild-type complex, which is recognized as an important structural requirement for restoring binding affinity. Despite numerous pockets on the protein surface of the mutant hGH-hGHbp complex, the binding-site cavity presents the energetically favorable hot spot for the benzimidazole-containing inhibitors, whereas for a set of nonactive molecules, the lowest energy ligand conformations do not necessarily bind in the engineered cavity. The results reveal a dominant role of the intermolecular van der Waals interactions in providing favorable ligand-protein energetics in the redesigned interface, in agreement with the experimental and computational alanine scanning of the hGH-hGHbp complex.     

Ornithine and Homocitrulline Impair Mitochondrial Function,Decrease Antioxidant Defenses and Induce Cell Death in Menadione-Stressed Rat Cortical Astrocytes: Potential Mechanisms of Neurological Dysfunction in HHH Syndrome

Hyperornithinemia–hyperammonemia–homocitrullinuria (HHH) syndrome is caused by deficiency of ornithine translocase leading to predominant tissue accumulation and high urinary excretion of ornithine (Orn), homocitrulline (Hcit) and ammonia. Although affected patients commonly present neurological dysfunction manifested by cognitive deficit, spastic paraplegia, pyramidal and extrapyramidal signs, stroke-like episodes, hypotonia and ataxia, its pathogenesis is still poorly known. Although astrocytes are necessary for neuronal protection. Therefore, in the present study we investigated the effects of Orn and Hcit on cell viability (propidium iodide incorporation), mitochondrial function (thiazolyl blue tetrazolium bromide—MTT—reduction and mitochondrial membrane potential—ΔΨ_m), antioxidant defenses (GSH) and pro-inflammatory response (NFkB, IL-1β, IL-6 and TNF-α) in unstimulated and menadione-stressed cortical astrocytes that were previously shown to be susceptible to damage by neurotoxins. We first observed that Orn decreased MTT reduction, whereas both amino acids decreased GSH levels, without altering cell viability and the pro-inflammatory factors in unstimulated astrocytes. Furthermore, Orn and Hcit decreased cell viability and ΔΨ_m in menadione-treated astrocytes. The present data indicate that the major compounds accumulating in HHH syndrome impair mitochondrial function and reduce cell viability and the antioxidant defenses in cultured astrocytes especially when stressed by menadione. It is presumed that these mechanisms may be involved in the neuropathology of this disease.     

AMPA receptors and bacterial periplasmic amino acid-binding proteins share the ionic mechanism of ligand recognition.

In order to identify key structural determinants for ligand recognition, we subjected the ligand-binding domain of the alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA)-selective glutamate receptor GluR-D subunit to site-directed mutagenesis. Based on the analysis of the [3H]AMPA-binding properties of the mutated binding sites, we constructed a revised three-dimensional model of the ligand-binding site, different in many respects from previously published models. In particular, our results indicate that the residues Arg507 and Glu727 represent the structural and functional correlates of Arg77 and Asp161 in the homologous bacterial lysine/ornithine/arginine-binding protein and histidine-binding protein, and directly interact with the alpha-carboxyl and alpha-amino group of the bound ligand, respectively. In contrast, Glu424, implicated previously in ionic interactions with the alpha-amino group of the agonist, is unlikely to have such a role in ligand binding. Our results indicate that glutamate receptors share with the bacterial polar amino acid-binding proteins the fundamental mechanism of amino acid recognition.     

Crystal structures of the choline/acetylcholine substrate-binding protein ChoX from Sinorhizobium meliloti in the liganded and unliganded-closed states

The ATP-binding cassette transporter ChoVWX is one of several choline import systems operating in Sinorhizobium meliloti. Here fluorescence-based ligand binding assays were used to quantitate substrate binding by the periplasmic ligand-binding protein ChoX. These data confirmed that ChoX recognizes choline and acetylcholine with high and medium affinity, respectively. We also report the crystal structures of ChoX in complex with either choline or acetylcholine. These structural investigations revealed an architecture of the ChoX binding pocket and mode of substrate binding similar to that reported previously for several compatible solute-binding proteins. Additionally the ChoX-acetylcholine complex permitted a detailed structural comparison with the carbamylcholine-binding site of the acetylcholine-binding protein from the mollusc Lymnaea stagnalis. In addition to the two liganded structures of ChoX, we were also able to solve the crystal structure of ChoX in a closed, substrate-free conformation that revealed an architecture of the ligand-binding site that is superimposable to the closed, ligand-bound form of ChoX. This structure is only the second of its kind and raises the important question of how ATP-binding cassette transporters are capable of distinguishing liganded and unliganded-closed states of the binding protein.     

Functional studies on the ligand-binding domain of Ultraspiracle from Drosophila melanogaster

The functional insect ecdysteroid receptor is comprised of the ecdysone receptor (EcR) and Ultraspiracle (USP). The ligand-binding domain (LBD) of USP was fused to the GAL4 DNA-binding domain (GAL4-DBD) and characterized by analyzing the effect of site-directed mutations in the LBD. Normal and mutant proteins were tested for ligand and DNA binding, dimerization, and their ability to induce gene expression. The presence of helix 12 proved to be essential for DNA binding and was necessary to confer efficient ecdysteroid binding to the heterodimer with the EcR (LBD), but did not influence dimerization. The antagonistic position of helix 12 is indispensible for interaction between the fusion protein and DNA, whereas hormone binding to the EcR (LBD) was only partially reduced if fixation of helix 12 was disturbed. The mutation of amino acids, which presumably bind to a fatty acid evoked a profound negative influence on transactivation ability, although enhanced transactivation potency and ligand binding to the ecdysteroid receptor was impaired to varying degrees by mutation of these residues. Mutations of one fatty acid-binding residue within the ligand-binding pocket, 1323, however, evoked enhanced transactivation. The results confirmed that the LBD of Ultraspiracle modifies ecdysteroid receptor function through intermolecular interactions and demonstrated that the ligand-binding pocket of USP modifies the DNA-binding and transactivation abilities of the fusion protein.     

Ligand displaces heat shock protein 90 from overlapping binding sites within the aryl hydrocarbon receptor ligand-binding domain

Hsp90 (heat shock protein of 90 kDa) is often found associated with functional domains of client proteins, including those for ligand binding, dimerization, DNA binding, and enzymatic activity. Although Hsp90 can maintain the conformation of functionally important domains prior to activation of the client protein, its specific binding site and the mechanism(s) of Hsp90 dissociation during activation are unknown. Here, we have identified and characterized residues involved in Hsp90 binding within the aryl hydrocarbon receptor (AhR) ligand-binding domain and demonstrate that they overlap with those involved in ligand binding. In agreement with this spatial model, ligand binding results in Hsp90 dissociation from the AhR Per-ARNT-Sim B fragment. Interestingly, whereas Hsp90-binding residues within the ligand-binding domain were not involved in Hsp90-dependent AhR protein stability, several of these residues are important for ligand-dependent AhR activation, and their mutation resulted in conversion of two AhR antagonists/partial agonists into full AhR agonists. These studies reveal co-localization of a tentative Hsp90-binding site with that for AhR ligand binding and provide the first molecular mechanism for Hsp90 dissociation in the activation of a client protein.     

A family-based approach reveals the function of residues in the nuclear receptor ligand-binding domain

Literature studies, 3D structure data, and a series of sequence analysis techniques were combined to reveal important residues in the structure and function of the ligand-binding domain of nuclear hormone receptors. A structure-based multiple sequence alignment allowed for the seamless combination of data from many different studies on different receptors into one single functional model. It was recently shown that a combined analysis of sequence entropy and variability can divide residues in five classes; (1) the main function or active site, (2) support for the main function, (3) signal transduction, (4) modulator or ligand binding and (5) the rest. Mutation data extracted from the literature and intermolecular contacts observed in nuclear receptor structures were analyzed in view of this classification and showed that the main function or active site residues of the nuclear receptor ligand-binding domain are involved in cofactor recruitment. Furthermore, the sequence entropy-variability analysis identified the presence of signal transduction residues that are located between the ligand, cofactor and dimer sites, suggesting communication between these regulatory binding sites. Experimental and computational results agreed well for most residues for which mutation data and intermolecular contact data were available. This allows us to predict the role of the residues for which no functional data is available yet. This study illustrates the power of family-based approaches towards the analysis of protein function, and it points out the problems and possibilities presented by the massive amounts of data that are becoming available in the "omics era". The results shed light on the nuclear receptor family that is involved in processes ranging from cancer to infertility, and that is one of the more important targets in the pharmaceutical industry.     

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.     

Molecular modeling of protein function regions

Experimental protein structures often provide extensive insight into the mode and specificity of small molecule binding, and this information is useful for understanding protein function and for the design of drugs. We have performed an analysis of the reliability with which ligand-binding information can be deduced from computer model structures, as opposed to experimentally derived ones. Models produced as part of the CASP experiments are used. The accuracy of contacts between protein model atoms and experimentally determined ligand atom positions is the main criterion. Only comparative models are included (i.e., models based on a sequence relationship between the protein of interest and a known structure). We find that, as expected, contact errors increase with decreasing sequence identity used as a basis for modeling. Analysis of the causes of errors shows that sequence alignment errors between model and experimental template have the most deleterious effect. In general, good, but not perfect, insight into ligand binding can be obtained from models based on a sequence relationship, providing there are no alignment errors in the model. The results support a structural genomics strategy based on experimental sampling of structure space so that all protein domains can be modeled on the basis of 30% or higher sequence identity.     

Trivalent recognition unit of innate immunity system: crystal structure of trimeric human M-ficolin fibrinogen-like domain

Ficolins are a kind of pathogen-recognition molecule in the innate immune systems. To investigate the discrimination mechanism between self and non-self by ficolins, we determined the crystal structure of the human M-ficolin fibrinogen-like domain (FD1), which is the ligand-binding domain, at 1.9A resolution. Although the FD1 monomer shares a common fold with the fibrinogen gamma fragment and tachylectin-5A, the Asp-282-Cys-283 peptide bond, which is the predicted ligand-binding site on the C-terminal P domain, is a normal trans bond, unlike the cases of the other two proteins. The trimeric formation of FD1 results in the separation of the three P domains, and the spatial arrangement of the three predicted ligand-binding sites on the trimer is very similar to that of the trimeric collectin, indicating that such an arrangement is generally required for pathogen-recognition. The ligand binding study of FD1 in solution indicated that the recombinant protein binds to N-acetyl-d-glucosamine and the peptide Gly-Pro-Arg-Pro and suggested that the ligand-binding region exhibits a conformational equilibrium involving cis-trans isomerization of the Asp-282-Cys-283 peptide bond. The crystal structure and the ligand binding study of FD1 provide an insight of the self- and non-self discrimination mechanism by ficolins.     

Molecular basis for the subtype discrimination of the estrogen receptor-beta-selective ligand,diarylpropionitrile

Although the two subtypes of the human estrogen receptor (ER), ERalpha and ERbeta, share only 56% amino acid sequence identity in their ligand binding domain (LBD), the residues that surround the ligand are nearly identical; nevertheless, subtype-selective ligands are known. To understand the molecular basis by which diarylpropionitrile (DPN), an ERbeta-selective ligand, is able to discriminate between the two ERs, we examined its activity on ER mutants and chimeric constructs generated by DNA shuffling. The N-terminal region of the ERbeta LBD (through helix 6) appears to be fully responsible for the ERbeta selectivity of DPN. In fact, a single ERalpha point mutation (L384M) was largely sufficient to switch the DPN response of this ER to that of the ERbeta type, but residues in helix 3 are also important in achieving the full ERbeta selectivity of DPN. Using molecular modeling, we found an energetically favorable fit for the S-DPN enantiomer in ERbeta, in which the proximal phenol mimics the A ring of estradiol, and the nitrile engages in stabilizing interactions with residues in the ligand-binding pocket of ERbeta. Our findings highlight that a limited number of critical interactions of DPN with the ERbeta ligand-binding pocket underlie its ER subtype-selective character.     

Improving accuracy and efficiency of blind protein-ligand docking by focusing on predicted binding sites

The use of predicted binding sites (binding sites calculated from the protein structure alone) is evaluated here as a tool to focus the docking of small molecule ligands into protein structures, simulating cases where the real binding sites are unknown. The resulting approach consists of a few independent docking runs carried out on small boxes, centered on the predicted binding sites, as opposed to one larger blind docking run that covers the complete protein structure. The focused and blind approaches were compared using a set of 77 known protein-ligand complexes and 19 ligand-free structures. The focused approach is shown to: (1) identify the correct binding site more frequently than blind docking; (2) produce more accurate docking poses for the ligand; (3) require less computational time. Additionally, the results show that very few real binding sites are missed in spite of focusing on only three predicted binding sites per target protein. Overall the results indicate that, by improving the sampling in regions that are likely to correspond to binding sites, the focused docking approach increases accuracy and efficiency of protein ligand docking for those cases where the ligand-binding site is unknown. This is especially relevant in applications such as reverse virtual screening and structure-based functional annotation of proteins.     

Crystal structure of the human TbetaR2 ectodomain--TGF-beta3 complex

Transforming growth factor-beta (TGF-beta) is the prototype of a large family of structurally related cytokines that play key roles in maintaining cellular homeostasis by signaling through two classes of functionally distinct Ser/Thr kinase receptors, designated as type I and type II. TGF-beta initiates receptor assembly by binding with high affinity to the type II receptor. Here, we present the 2.15 A crystal structure of the extracellular ligand-binding domain of the human TGF-beta type II receptor (ecTbetaR2) in complex with human TGF-beta3. ecTbetaR2 interacts with homodimeric TGF-beta3 by binding identical finger segments at opposite ends of the growth factor. Relative to the canonical 'closed' conformation previously observed in ligand structures across the superfamily, ecTbetaR2-bound TGF-beta3 shows an altered arrangement of its monomeric subunits, designated the 'open' conformation. The mode of TGF-beta3 binding shown by ecTbetaR2 is compatible with both ligand conformations. This, in addition to the predicted mode for TGF-beta binding to the type I receptor ectodomain (ecTbetaR1), suggests an assembly mechanism in which ecTbetaR1 and ecTbetaR2 bind at adjacent positions on the ligand surface and directly contact each other via protein--protein interactions.     

Evidence for DNA-binding domain--ligand-binding domain communications in the androgen receptor

DNA binding as well as ligand binding by nuclear receptors has been studied extensively. Both binding functions are attributed to isolated domains of which the structure is known. The crystal structure of a complete receptor in complex with its ligand and DNA-response element, however, has been solved only for the peroxisome proliferator-activated receptor γ (PPARγ)-retinoid X receptor α (RXRα) heterodimer. This structure provided the first indication of direct interactions between the DNA-binding domain (DBD) and ligand-binding domain (LBD). In this study, we investigated whether there is a similar interface between the DNA- and ligand-binding domains for the androgen receptor (AR). Despite the structural differences between the AR- and PPARγ-LBD, a combination of in silico modeling and docking pointed out a putative interface between AR-DBD and AR-LBD. The surfaces were subjected to a point mutation analysis, which was inspired by known AR mutations described in androgen insensitivity syndromes and prostate cancer. Surprisingly, AR-LBD mutations D695N, R710A, F754S, and P766A induced a decrease in DNA binding but left ligand binding unaffected, while the DBD-residing mutations K590A, K592A, and E621A lowered the ligand-binding but not the DNA-binding affinity. We therefore propose that these residues are involved in allosteric communications between the AR-DBD and AR-LBD.     

Crystal structure of Bordetella pertussis BugD solute receptor unveils the basis of ligand binding in a new family of periplasmic binding proteins

Periplasmic binding proteins of a new family particularly well represented in Bordetella pertussis have been called Bug receptors. One B.pertussis Bug protein is part of a tripartite tricarboxylate transporter while the functions of the other 77 are unknown. We report the first structure of a Bug receptor, BugD. It adopts the characteristic Venus flytrap motif observed in other periplasmic binding proteins, with two globular domains bisected by a deep cleft. BugD displays a closed conformation resulting from the fortuitous capture of a ligand, identified from the electron density as an aspartate. The structure reveals a distinctive alpha carboxylate-binding motif, involving two water molecules that bridge the carboxylate oxygen atoms to the protein. Both water molecules are hydrogen bonded to a common carbonyl group from Ala14, and each forms a hydrogen bond with one carboxylate oxygen atom of the ligand. Additional hydrogen bonds are found between the ligand alpha carboxylate oxygen atoms and protein backbone amide groups and with a threonine hydroxyl group. This specific ligand-binding motif is highly conserved in Bug proteins, indicating that they may all be receptors of amino acids or other carboxylated solutes, with a similar binding mode. The present structure thus unveils the bases of ligand binding in this large family of periplasmic binding proteins, several hundred members of which have been identified in various bacterial species.     

Three-dimensional structure and active site of three hydrophobic molecule-binding proteins with significant amino acid sequence similarity.

We review our work on bovine and human retinol-binding protein (RBP), bovine beta lactoglobulin (BLG), and bovine odorant-binding protein (OBP). These three proteins share a sequence similarity high enough to justify the proposal that their three-dimensional structure ought to be quite similar, and they also share the function of similar or even identical hydrophobic ligand binding, although with a very different degree of specificity. Thus they constitute an ideal system to exhaustively explore the question of three-dimensional structure prediction from sequence similarity and the related question of binding site prediction for similar ligands. We have used x-ray diffraction techniques on single crystals of human and bovine RBP, bovine milk BLG, and bovine nasal mucosa OBP to investigate this problem. The results of these crystallographic studies indicate that to the level of resolution so far attained, the three-dimensional structure of these three proteins is reasonably predicted from the sequence similarity. The fold is the same and structural differences are rather subtle. Finally, we present experimental evidence that the binding sites of RBP, BLG, and OBP are in different regions of the molecules. Thus, it appears that although sequence alignment has correctly predicted the protein fold, it has incorrectly predicted the hydrophobic ligand-binding sites.     

Structural evidence for ligand specificity in the binding domain of the human androgen receptor. Implications for pathogenic gene mutations

The crystal structures of the human androgen receptor (hAR) and human progesterone receptor ligand-binding domains in complex with the same ligand metribolone (R1881) have been determined. Both three-dimensional structures show the typical nuclear receptor fold. The change of two residues in the ligand-binding pocket between the human progesterone receptor and hAR is most likely the source for the specificity of R1881 to the hAR. The structural implications of the 14 known mutations in the ligand-binding pocket of the hAR ligand-binding domains associated with either prostate cancer or the partial or complete androgen receptor insensitivity syndrome were analyzed. The effects of most of these mutants could be explained on the basis of the crystal structure.