Inward facing conformations of the MetNI methionine ABC transporter: Implications for the mechanism of transinhibition

Two new crystal structures of the Escherichia coli high affinity methionine uptake ATP Binding Cassette (ABC) transporter MetNI, purified in the detergents cyclohexyl‐pentyl‐β‐D ‐maltoside (CY5) and n‐decyl‐β‐D ‐maltopyranoside (DM), have been solved in inward facing conformations to resolutions of 2.9 and 4.0 Å, respectively. Compared to the previously reported 3.7 Å resolution structure of MetNI purified in n‐dodecyl‐β‐D ‐maltopyranoside (DDM), the higher resolution of the CY5 data enabled significant improvements to the structural model in several regions, including corrections to the sequence registry, and identification of ADP in the nucleotide binding site. CY5 crystals soaked with selenomethionine established details of the methionine binding site in the C2 regulatory domain of the ABC subunit, including the displacement of the side chain of MetN residue methionine 301 by the exogenous ligand. When compared to the CY5 or DDM structures, the DM structure exhibits a significant repositioning of the dimeric C2 domains, including an unexpected register shift in the intermolecular β‐sheet hydrogen bonding between monomers, and a narrowing of the nucleotide binding space. The immediate proximity of the exogenous methionine binding site to the conformationally variable dimeric interface provides an indication of how methionine binding to the regulatory domains might mediate the phenomenon of transinhibition.     

Quantum mechanical analysis of oxygenated and deoxygenated states of hemocyanin: theoretical clues for a plausible allosteric model of oxygen binding.

Abstract: In this work with ab initio computations, we describe relevant interactions between protein active sites and ligands, using as a test case arthropod hemocyanins. A computational analysis of models corresponding to the oxygenated and deoxygenated forms of the hemocyanin active site is performed using the Density Functional Theory approach. We characterize the electron density distribution of the binding site with and without bound oxygen in relation to the geometry, which stems out of the crystals of three hemocyanin proteins, namely the oxygenated form from the horseshoe crab Limulus polyphemus, and the deoxygenated forms, respectively, from the same source and from another arthropod, the spiny lobster Panulirus interruptus. Comparison of the three available crystals indicate structural differences at the oxygen binding site, which cannot be explained only by the presence and absence of the oxygen ligand, since the geometry of the ligand site of the deoxygenated Panulirus hemocyanin is rather similar to that of the oxygenated Limulus protein. This finding was interpreted in the frame of a mechanism of allosteric regulation for oxygen binding. However, the cooperative mechanism, which is experimentally well documented, is only partially supported by crystallographic data, since no oxygenated crystal of Panulirus hemocyanin is presently available. We address the following question: is the local ligand geometry responsible for the difference of the dicopper distance observed in the two deoxygenated forms of hemocyanin or is it necessary to advocate the allosteric regulation of the active site conformations in order to reconcile the different crystal forms? We find that the difference of the dicopper distance between the two deoxygenated hemocyanins is not due to the small differences of ligand geometry found in the crystals and conclude that it must be therefore stabilized by the whole protein tertiary structure.     

Structural basis for the negative allostery between Ca(2+)- and Mg(2+)-binding in the intracellular Ca(2+)-receptor calbindin D9k.

The three-dimensional structures of the magnesium- and manganese-bound forms of calbindin D9k were determined to 1.6 A and 1.9 A resolution, respectively, using X-ray crystallography. These two structures are nearly identical but deviate significantly from both the calcium bound form and the metal ion-free (apo) form. The largest structural differences are seen in the C-terminal EF-hand, and involve changes in both metal ion coordination and helix packing. The N-terminal calcium binding site is not occupied by any metal ion in the magnesium and manganese structures, and shows little structural deviation from the apo and calcium bound forms. 1H-NMR and UV spectroscopic studies at physiological ion concentrations show that the C-terminal site of the protein is significantly populated by magnesium at resting cell calcium levels, and that there is a negative allosteric interaction between magnesium and calcium binding. Calcium binding was found to occur with positive cooperativity at physiological magnesium concentration.     

An autoantibody to single-stranded DNA: comparison of the three-dimensional structures of the unliganded Fab and a deoxynucleotide-Fab complex.

Crystal structures of the Fabs from an autoantibody (BV04-01) with specificity for single-stranded DNA have been determined in the presence and absence of a trinucleotide of deoxythymidylic acid, d(pT)3. Formation of the ligand-protein complex was accompanied by small adjustments in the orientations of the variable (VL and VH) domains. In addition, there were local conformational changes in the first hypervariable loop of the light chain and the third hypervariable loop of the heavy chain, which together with the domain shifts led to an improvement in the complementarity of nucleotide and Fab. The sugar-phosphate chain adopted an extended and "open" conformation, with the base, sugar, and phosphate components available for interactions with the protein. Nucleotide 1 (5'-end) was associated exclusively with the heavy chain, nucleotide 2 was shared by both heavy and light chains, and nucleotide 3 was bound by the light chain. The orientation of phosphate 1 was stabilized by hydrogen bonds with serine H52a and asparagine H53. Phosphate 2 formed an ion pair with arginine H52, but no other charge-charge interactions were observed. Insertion of the side chain of histidine L27d between nucleotides 2 and 3 resulted in a bend in the sugar-phosphate chain. The most dominant contacts with the protein involved the central thymine base, which was immobilized by cooperative stacking and hydrogen bonding interactions. This base was intercalated between a tryptophan ring (no. H100a) from the heavy chain and a tyrosine ring (no. L32) from the light chain. The resulting orientation of thymine was favorable for the simultaneous formation of two hydrogen bonds with the backbone carbonyl oxygen and the side chain hydroxyl group of serine L91 (the thymine atoms were the hydrogen on nitrogen 3 and keto oxygen 4).     

Computation of electrostatic complements to proteins: a case of charge stabilized binding.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.     

Understanding the basis of a class of paradoxical mutations in AraC through simulations

Most mutations at position 15 in the N‐terminal arm of the regulatory protein AraC leave the protein incapable of responding to arabinose and inducing the proteins required for arabinose catabolism. Mutations at other positions of the arm do not have this behavior. Simple energetic analysis of the interactions between the arm and bound arabinose do not explain the uninducibility of AraC with mutations at position 15. Extensive molecular dynamics (MD) simulations, carried out largely on the Open Science Grid, were done of the wild‐type protein with and without bound arabinose and of all possible mutations at position 15, many of which were constructed and measured for this work. Good correlation was found for deviation of arm position during the simulations and inducibility as measured in vivo of the same mutant proteins. Analysis of the MD trajectories revealed that preservation of the shape of the arm is critical to inducibility. To maintain the correct shape of the arm, the strengths of three interactions observed to be strong in simulations of the wild‐type AraC protein need to be preserved. These interactions are between arabinose and residue 15, arabinose and residues 8–9, and residue 13 and residue 15. The latter interaction is notable because residues L9, Y13, F15, W95, and Y97 form a hydrophobic cluster which needs to be preserved for retention of the correct shape. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The upper cytokine-binding module and the Ig-like domain of the leukaemia inhibitory factor (LIF) receptor are sufficient for a functional LIF receptor complex.

To elucidate the function of the two cytokine-binding modules (CBM) of the leukemia inhibitory factor receptor (LIFR), receptor chimeras of LIFR and the interleukin-6 receptor (IL-6R) were constructed. Either the NH(2)-terminal (chimera RILLIFdeltaI) or the COOH-terminal LIFR CBM (chimera RILLIFdeltaII) were replaced by the structurally related CBM of the IL-6R which does not bind LIF. Chimera RILLIFdeltaI is functionally inactive, whereas RILLIFdeltaII binds LIF and mediates signalling as efficiently as the wild-type LIFR. Deletion mutants of the LIFR revealed that both the NH(2)-terminal CBM and the Ig-like domain of the LIFR are involved in LIF binding, presumably via the LIF site III epitope. The main function of the COOH-terminal CBM of the LIFR is to position the NH(2)-terminal CBM and the Ig-like domain, so that these can bind to LIF. In analogy to a recently published model of the IL-6R complex, a model of the active LIFR complex is suggested which positions the COOH-terminal CBM at LIF site I and the NH(2)-terminal CBM and the Ig-like domain at site III. An additional contact is postulated between the Ig-like domain of gp130 and the NH(2)-terminal CBM of the LIFR.     

A study of a hierarchical structure of proteins and ligand binding sites of receptors using the triangular spatial relationship-based structure comparison method and development of a size-filtering feature designed for comparing different sizes of protein structures

The presence of receptors and the specific binding of the ligands determine nearly all cellular responses. Binding of a ligand to its receptor causes conformational changes of the receptor that triggers the subsequent signaling cascade. Therefore, systematically studying structures of receptors will provide insight into their functions. We have developed the triangular spatial relationship (TSR)-based method where all possible triangles are constructed with C_α atoms of a protein as vertices. Every triangle is represented by an integer denoted as a “key” computed through the TSR algorithm. A structure is thereby represented by a vector of integers. In this study, we have first defined substructures using different types of keys. Second, using different types of keys represents a new way to interpret structure hierarchical relations and differences between structures and sequences. Third, we demonstrate the effects of sequence similarity as well as sample size on the structure-based classifications. Fourth, we show identification of structure motifs, and the motifs containing multiple triangles connected by either an edge or a vertex are mapped to the ligand binding sites of the receptors. The structure motifs are valuable resources for the researchers in the field of signal transduction. Next, we propose amino-acid scoring matrices that capture “evolutionary closeness” information based on BLOSUM62 matrix, and present the development of a new visualization method where keys are organized according to evolutionary closeness and shown in a 2D image. This new visualization opens a window for developing tools with the aim of identification of specific and common substructures by scanning pixels and neighboring pixels. Finally, we report a new algorithm called as size filtering that is designed to improve structure comparison of large proteins with small proteins. Collectively, we provide an in-depth interpretation of structure relations through the detailed analyses of different types of keys and their associated key occurrence frequencies, geometries, and labels. In summary, we consider this study as a new computational platform where keys are served as a bridge to connect sequence and structure as well as structure and function for a deep understanding of sequence, structure, and function relationships of the protein family.     

Molecular basis of film formation from a soybean protein: comparison between the conformation of glycinin in aqueous solution and in films

Fourier transform infrared spectroscopy has been used to investigate the conformational changes of glycinin, a major storage protein of soybean seeds, upon film-forming. The results show that the secondary structure of glycinin is mainly composed of a β-sheet (48%) and unordered (49%) structures. The amide I band of glycinin in film-forming conditions, i.e. in alkaline media and in the presence of plasticizing agent, reveals the conversion of 18% of the secondary structure of the protein from the β-sheet (6%) and random coil (12%) to the -helical conformation due to the helicogenic effect of the ethylene glycol used as the plasticizing agent. Conformational changes also occur upon the film-forming process leading to the formation of intermolecular hydrogen-bonded β-sheet structures. Results obtained from other plant families indicate that, whatever the origin and conformation of protein, formation of films leads to the appearance of intermolecular hydrogen-bonded β-sheet structures, suggesting that this type of structure might be essential for the network formation in films. Thus, it is hypothesized that, in the film state, intermolecular hydrogen bonding between segments of β-sheet may act as junction zones in the film network. This study reveals for the first time that there is a close relationship between the conformation of proteins and the mechanical properties of films.     

SuperStar: comparison of CSD and PDB-based interaction fields as a basis for the prediction of protein-ligand interactions.

SuperStar is an empirical method for identifying interaction sites in proteins, based entirely on the experimental information about non-bonded interactions, present in the IsoStar database. The interaction information in IsoStar is contained in scatterplots, which show the distribution of a chosen probe around structure fragments. SuperStar breaks a template molecule (e.g. a protein binding site) into structural fragments which correspond to those in the scatterplots. The scatterplots are then superimposed on the corresponding parts of the template and converted into a composite propensity map.The original version of SuperStar was based entirely on scatterplots from the CSD. Here, scatterplots based on protein-ligand interactions are implemented in SuperStar, and validated on a test set of 122 X-ray structures of protein-ligand complexes. In this validation, propensity maps are compared with the experimentally observed positions of ligand atoms of comparable types. Although non-bonded interaction geometries in small molecule structures are similar to those found in protein-ligand complexes, their relative frequencies of occurrence are different. Polar interactions are more common in the first class of structures, while interactions between hydrophobic groups are more common in protein crystals. In general, PDB and CSD-based SuperStar maps appear equally successful in the prediction of protein-ligand interactions. PDB-based maps are more suitable to identify hydrophobic pockets, and inherently take into account the experimental uncertainties of protein atomic positions. If the protonation state of a histidine, aspartate or glutamate protein side-chain is known, specific CSD-based maps for that protonation state are preferred over PDB-based maps which represent an ensemble of protonation states.     

FoldGPCR: Structure prediction protocol for the transmembrane domain of G protein‐coupled receptors from class A

Building reliable structural models of G protein‐coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template‐based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1–3 Å Cα RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter‐residue contacts between the template and target, and exploits an all‐atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor‐ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human β₂‐adrenergic receptor (β₂AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X‐ray crystal structure of the β₂AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix‐helix packing interactions that distinguish rhodopsin from β₂AR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Profiling two indole‐2‐carboxamides for allosteric modulation of the CB1 receptor

Allosteric modulation of G‐protein coupled receptors (GPCRs) represents a novel approach for fine‐tuning GPCR functions. The cannabinoid CB1 receptor, a GPCR associated with the CNS, has been implicated in the treatment of drug addiction, pain, and appetite disorders. We report here the synthesis and pharmacological characterization of two indole‐2‐carboxamides:5‐chloro‐3‐ethyl‐1‐methyl‐N‐(4‐(piperidin‐1‐yl)phenethyl)‐1H‐indole‐2‐carboxamide (ICAM‐a) and 5‐chloro‐3‐pentyl‐N‐(4‐(piperidin‐1‐yl)phenethyl)‐1H‐indole‐2‐carboxamide (ICAM‐b). Although both ICAM‐a and ICAM‐b enhanced CP55, 940 binding, ICAM‐b exhibited the strongest positive cooperativity thus far demonstrated for enhancing agonist binding to the CB1 receptor. Although it displayed negative modulatory effects on G‐protein coupling to CB1, ICAM‐b induced β‐arrestin‐mediated downstream activation of extracellular signal‐regulated kinase (ERK) signaling. These results indicate that this compound represents a novel class of CB1 ligands that produce biased signaling via CB1.     

CLIX: a search algorithm for finding novel ligands capable of binding proteins of known three-dimensional structure.

A computer algorithm, CLIX, capable of searching a crystallographic data-base of small molecules for candidates which have both steric and chemical likelihood of binding a protein of known three-dimensional structure is presented. The algorithm is a significant advance over previous strategies which consider solely steric or chemical requirements for binding. The algorithm is shown to be capable of predicting the correct binding geometry of sialic acid to a mutant influenza-virus hemagglutinin and of proposing a number of potential new ligands to this protein.     

Modeling of the structure of the Haemophilus influenzae heme-binding protein suggests a mode of heme interaction.

The structure and function of the periplasmic heme-binding protein HbpA of Haemophilus influenzae were investigated. This protein is involved in the import of heme into the bacteria through the inner membrane, and thus is a key element of the organism's ability to survive in blood. A high degree of sequence similarity between HbpA and the dipeptide-binding protein of Escherichia coli is suggested to be the result of a functional relationship. An HbpA model built using the dipeptide-binding protein suggests a mode of heme binding that is distinct from those known in proteins of the human host. These results provide a starting point for rational drug design.     

Molecular strategies for protein stabilization: the case of a trehalose/maltose-binding protein from Thermus thermophilus

The trehalose/maltose-binding protein (MalE1) is one component of trehalose and maltose uptake system in the thermophilic organism Thermus thermophilus. MalE1 is a monomeric 48 kDa protein predominantly organized in alpha-helix conformation with a minor content of beta-structure. In this work, we used Fourier-infrared spectroscopy and in silico methodologies for investigating the structural stability properties of MalE1. The protein was studied in the absence and in the presence of maltose as well as in the absence and in the presence of SDS at different p(2)H values (neutral p(2)H and at p(2)H 9.8). In the absence of SDS, the results pointed out a high thermostability of the MalE1 alpha-helices, maintained also at basic p(2)H values. However, the obtained data also showed that at high temperatures the MalE1 beta-sheets underwent to structural rearrangements that were totally reversible when the temperature was lowered. At room temperature, the addition of SDS to the protein solution slightly modified the MalE1 secondary structure content by decreasing the protein thermostability. The infrared data, corroborated by molecular dynamics simulation experiments performed on the structure of MalE1, indicated that the protein hydrophobic interactions have an important role in the MalE1 high thermostability. Finally, the results obtained on MalE1 are also discussed in comparison with the data on similar thermostable proteins already studied in our laboratories.