Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Isoelectric Point Determination for Glossoscolex paulistus Extracellular Hemoglobin: Oligomeric Stability in Acidic pH and Relevance to Protein-Surfactant Interactions

The extracellular hemoglobin from Glossoscolex paulistus (HbGp) has a molecular mass of 3.6 M Da, It has a high oligomeric stability at pH 7.0 and low autoxidation rates, as compared to vertebrate hemoglobins. In this work, fluorescence and light scattering experiments were performed with the three oxidation forms of HbGp exposed to acidic pH. Our focus is on the HbGp stability at acidic pH and also on the determination of the isoelectric point (pI) of the protein. Our results show that the protein in the cyanomet form is more stable than in the other two forms, in the whole range. Our zeta-potential data are consistent with light scattering results. Average values apt obtained by different techniques were 5.6 +/- 0.5, 5.4 +/- 0.2 and 5.2 +/- 0.5 for the oxy, met, and cyanomet forms. Dynamic light scattering (DLS) experiments have shown that, at pH 6.0, the aggregation (oligomeric) state of oxy-, met- and cyanomet-HbGp remains the same as that at 7.0. The interaction between the oxy-HbGp and ionic surfactants at pH 5.0 and 6.0 was also monitored in the present study. At pH 5,0, below the protein pI, the effects of sodium dodecyl sulfate (SDS) and cetyltrimethylammonium chloride (CTAC) are inverted when compared to pH 7.0. For CTAC, in acid pH 5.0, no precipitation is observed, while for SDS an intense light scattering appears due to a precipitation process. HbGp interacts strongly with the cationic surfactant at pH 7.0 and with the anionic one at pH 5.0. This effect is due to the predominance, in the protein surface, of residues presenting opposite charges to the surfactant headgroups. This information can be relevant for the development of extracellular hemoglobin-based artificial blood substitutes.

Ric-8A, a G alpha protein guanine nucleotide exchange factor potentiates taste receptor signaling

Taste receptors for sweet, bitter and umami tastants are G-protein-coupled receptors (GPCRs). While much effort has been devoted to understanding G-protein-receptor interactions and identifying the components of the signalling cascade downstream of these receptors, at the level of the G-protein the modulation of receptor signal transduction remains relatively unexplored. In this regard a taste-specific regulator of G-protein signaling (RGS), RGS21, has recently been identified. To study whether guanine nucleotide exchange factors (GEFs) are involved in the transduction of the signal downstream of the taste GPCRs we investigated the expression of Ric-8A and Ric-8B in mouse taste cells and their interaction with G-protein subunits found in taste buds. Mammalian Ric-8 proteins were initially identified as potent GEFs for a range of G alpha subunits and Ric-8B has recently been shown to amplify olfactory signal transduction. We find that both Ric-8A and Ric-8B are expressed in a large portion of taste bud cells and that most of these cells contain IP3R-3 a marker for sweet, umami and bitter taste receptor cells. Ric-8A interacts with G alpha-gustducin and G alpha i2 through which it amplifies the signal transduction of hTas2R16, a receptor for bitter compounds. Overall, these findings are consistent with a role for Ric-8 in mammalian taste signal transduction.

Cloning, expression, molecular modelling and docking analysis of glutathione transferase from Saccharum officinarum

Sugarcane yield and quality are affected by a number of biotic and abiotic stresses. In response to such stresses, plants may increase the activities of some enzymes such as glutathione transferase (GST), which are involved in the detoxification of xenobiotics. Thus, a sugarcane GST was modelled and molecular docked using the program LIGIN to investigate the contributions of the active site residues towards the binding of reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB). As a result, W13 and I119 were identified as key residues for the specificity of sugarcane GSTF1 (SoGSTF1) towards CDNB. To obtain a better understanding of the catalytic specificity of sugarcane GST (SoGSTF1), two mutants were designed, W13L and I119F. Tertiary structure models and the same docking procedure were performed to explain the interactions between sugarcane GSTs with GSH and CDNB. An electron-sharing network for GSH interaction was also proposed. The SoGSTF1 and the mutated gene constructions were cloned and expressed in Escherichia coli and the expressed protein purified. Kinetic analyses revealed different Km values not only for CDNB, but also for GSH. The Km values were 0.2, 1.3 and 0.3 mM for GSH, and 0.9, 1.2 and 0.5 mM for CDNB, for the wild type, W13L mutant and I119F mutant, respectively. The V(max) values were 297.6, 224.5 and 171.8 mu mol min(-1) mg(-1) protein for GSH, and 372.3, 170.6 and 160.4 mu mol min(-1) mg(-1) protein for CDNB.

The Xanthomonas citri effector protein PthA interacts with citrus proteins involved in nuclear transport, protein folding and ubiquitination associated with DNA repair

P>Xanthomonas axonopodis pv. citri utilizes the type III effector protein PthA to modulate host transcription to promote citrus canker. PthA proteins belong to the AvrBs3/PthA family and carry a domain comprising tandem repeats of 34 amino acids that mediates protein-protein and protein-DNA interactions. We show here that variants of PthAs from a single bacterial strain localize to the nucleus of plant cells and form homo- and heterodimers through the association of their repeat regions. We hypothesize that the PthA variants might also interact with distinct host targets. Here, in addition to the interaction with alpha-importin, known to mediate the nuclear import of AvrBs3, we describe new interactions of PthAs with citrus proteins involved in protein folding and K63-linked ubiquitination. PthAs 2 and 3 preferentially interact with a citrus cyclophilin (Cyp) and with TDX, a tetratricopeptide domain-containing thioredoxin. In addition, PthAs 2 and 3, but not 1 and 4, interact with the ubiquitin-conjugating enzyme complex formed by Ubc13 and ubiquitin-conjugating enzyme variant (Uev), required for K63-linked ubiquitination and DNA repair. We show that Cyp, TDX and Uev interact with each other, and that Cyp and Uev localize to the nucleus of plant cells. Furthermore, the citrus Ubc13 and Uev proteins complement the DNA repair phenotype of the yeast Delta ubc13 and Delta mms2/uev1a mutants, strongly indicating that they are also involved in K63-linked ubiquitination and DNA repair. Notably, PthA 2 affects the growth of yeast cells in the presence of a DNA damage agent, suggesting that it inhibits K63-linked ubiquitination required for DNA repair.

Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

Crystallographic Structure of Phosphofructokinase-2 from Escherichia coli in Complex with Two ATP Molecules. Implications for Substrate Inhibition

Phosphofructokinase-1 and -2 (Pfk-1 and Pfk-2, respectively) from Escherichia coli belong to different homologous superfamilies. However, in spite of the lack of a common ancestor, they share the ability to catalyze the same reaction and are inhibited by the substrate MgATP. Pfk-2, an ATP-dependent 6-phosphofructokinase member of the ribokinase-like superfamily, is a homodimer of 66 kDa subunits whose oligomerization state is necessary for catalysis and stability. The presence of MgATP favors the tetrameric form of the enzyme. In this work, we describe the structure of Pfk-2 in its inhibited tetrameric form, with each subunit bound to two ATP molecules and two Mg ions. The present structure indicates that substrate inhibition occurs due to the sequential binding of two MgATP molecules per subunit, the first at the usual site occupied by the nucleotide in homologous enzymes and the second at the allosteric site, making a number of direct and Mg-mediated interactions with the first. Two configurations are observed for the second MgATP, one of which involves interactions with Tyr23 from the adjacent subunit in the dimer and the other making an unusual non-Watson-Crick base pairing with the adenine in the substrate ATP. The oligomeric state observed in the crystal is tetrameric, and some of the structural elements involved in the binding of the Substrate and allosteric ATPs are also participating in the dimer-dimer interface. This structure also provides the grounds to compare analogous features of the nonhomologous phosphofructokinases from E. coli. (C) 2008 Elsevier Ltd. All rights reserved.

Crystallographic structure and substrate-binding interactions of the molybdate-binding protein of the phytopathogen Xanthomonas axonopodis pv. citri

In Xanthomonas axonopodis pv. citri (Xac or X citri), the modA gene codes for a periplasmic protein (ModA) that is capable of binding molybdate and tungstate as part of the ABC-type transporter required for the uptake of micronutrients. In this study, we report the crystallographic structure of the Xac ModA protein with bound molybdate. The Xac ModA structure is similar to orthologs with known three-dimensional structures and consists of two nearly symmetrical domains separated by a hinge region where the oxyanion-binding site lies. Phylogenetic analysis of different ModA orthologs based on sequence alignments revealed three groups of molybdate-binding proteins: bacterial phytopathogens, enterobacteria and soil bacteria. Even though the ModA orthologs are segregated into different groups, the ligand-binding hydrogen bonds are mostly conserved, except for Archaeglobus fulgidus ModA. A detailed discussion of hydrophobic interactions in the active site is presented and two new residues, Ala(38) and Ser(151), are shown to be part of the ligand-binding pocket. (c) 2007 Elsevier B.V All rights reserved.

The Function and Three-Dimensional Structure of a Thromboxane A(2)/Cysteinyl Leukotriene-Binding Protein from the Saliva of a Mosquito Vector of the Malaria Parasite

The highly expressed D7 protein family of mosquito saliva has previously been shown to act as an anti-inflammatory mediator by binding host biogenic amines and cysteinyl leukotrienes (CysLTs). In this study we demonstrate that AnSt-D7L1, a two-domain member of this group from Anopheles stephensi, retains the CysLT binding function seen in the homolog AeD7 from Aedes aegypti but has lost the ability to bind biogenic amines. Unlike any previously characterized members of the D7 family, AnSt-D7L1 has acquired the important function of binding thromboxane A(2) (TXA(2)) and its analogs with high affinity. When administered to tissue preparations, AnSt-D7L1 abrogated Leukotriene C(4) (LTC(4))-induced contraction of guinea pig ileum and contraction of rat aorta by the TXA(2) analog U46619. The protein also inhibited platelet aggregation induced by both collagen and U46619 when administered to stirred platelets. The crystal structure of AnSt-D7L1 contains two OBP-like domains and has a structure similar to AeD(7). In AnSt-D7L1, the binding pocket of the C-terminal domain has been rearranged relative to AeD7, making the protein unable to bind biogenic amines. Structures of the ligand complexes show that CysLTs and TXA(2) analogs both bind in the same hydrophobic pocket of the N-terminal domain. The TXA(2) analog U46619 is stabilized by hydrogen bonding interactions of the omega-5 hydroxyl group with the phenolic hydroxyl group of Tyr 52. LTC(4) and occupies a very similar position to LTE(4) in the previously determined structure of its complex with AeD7. As yet, it is not known what, if any, new function has been acquired by the rearranged C-terminal domain. This article presents, to our knowledge, the first structural characterization of a protein from mosquito saliva that inhibits collagen mediated platelet activation.

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations

Background: Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP(II)), which is the most variable segment of the protein. Methods: To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP(II) in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP(II), and T-and B-cell epitopes were localized on the 3-D structure. Results: The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP(II), and (ii) PvDBP(II) appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium. Conclusions: This study shows that some polymorphisms of PvDBP(II) are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

The Importance of Protein-Protein Interactions on the pH-Induced Conformational Changes of Bovine Serum Albumin: A Small-Angle X-Ray Scattering Study

The combined effects of concentration and pH on the conformational states of bovine serum albumin (BSA) are investigated by small-angle x-ray scattering. Serum albumins, at physiological conditions, are found at concentrations of similar to 35-45 mg/mL (42 mg/mL in the case of humans). In this work, BSA at three different concentrations (10, 25, and 50 mg/mL) and pH values (2.0-9.0) have been studied. Data were analyzed by means of the Global Fitting procedure, with the protein form factor calculated from human serum albumin (HSA) crystallographic structure and the interference function described, considering repulsive and attractive interaction potentials within a random phase approximation. Small-angle x-ray scattering data show that BSA maintains its native state from pH 4.0 up to 9.0 at all investigated concentrations. A pH-dependence of the absolute net protein charge is shown and the charge number per BSA is quantified to 10(2), 8(l), 13(2), 20(2), and 26(2) for pH values 4.0, 5.4, 7.0, 8.0, and 9.0, respectively. The attractive potential diminishes as BSA concentration increases. The coexistence of monomers and dimers is observed at 50 mg/mL and pH 5.4, near the BSA isoelectric point. Samples at pH 2.0 show a different behavior, because BSA overall shape changes as a function of concentration. At 10 mg/mL, BSA is partially unfolded and a strong repulsive protein-protein interaction occurs due to the high amount of exposed charge. At 25 and 50 mg/mL, BSA undergoes some refolding, which likely results in a molten-globule state. This work concludes by confirming that the protein concentration plays an important role on the pH-unfolded BSA state, due to a delicate compromise between interaction forces and crowding effects.

Protein lethality investigated in terms of long range dynamical interactions

The relationship between network structure/dynamics and biological function constitutes a fundamental issue in systems biology. However, despite many related investigations, the correspondence between structure and biological functions is not yet fully understood. A related subject that has deserved particular attention recently concerns how essentiality is related to the structure and dynamics of protein interactions. In the current work, protein essentiality is investigated in terms of long range influences in protein-protein interaction networks by considering simulated dynamical aspects. This analysis is performed with respect to outward activations, an approach which models the propagation of interactions between proteins by considering self-avoiding random walks. The obtained results are compared to protein local connectivity. Both the connectivity and the outward activations were found to be strongly related to protein essentiality.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

The binding of synthetic triiodo L-thyronine analogs to human transthyretin: Molecular basis of cooperative and non-cooperative ligand recognition

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. Several classes of small molecules can bind to TTR delaying its amyloid fibril formation, thus being promising drug candidates to treat TTR amyloidoses. In the present study, we characterized the interactions of the synthetic triiodo L-thyronine analogs and thyroid hormone nuclear receptor TR beta-selecfive agonists GC-1 and GC-24 with the wild type and V30M variant of human transthyretin (TTR). To achieve this aim, we conducted in vitro TTR acid-mediated aggregation and isothermal titration calorimetry experiments and determined the TTR:GC-1 and TTR:GC-24 crystal structures. Our data indicate that both GC-1 and GC-24 bind to TTR in a non-cooperative manner and are good inhibitors of TTR aggregation, with dissociation constants for both hormone binding sites (HBS) in the low micromolar range. Analysis of the crystal structures of TTRwt:GC-1(24) complexes and their comparison with the TTRwt X-ray structure bound to its natural ligand thyroxine (T4) suggests, at the molecular level, the basis for the cooperative process displayed by T4 and the non-cooperative process provoked by both GC-1 and GC-24 during binding to TTR. (C) 2010 Elsevier Inc. All rights reserved.