The multiple phenotypes of allosteric receptor mutants.

Channel-linked neurotransmitter receptors are membrane-bound heterooligomers made up of distinct, although homologous, subunits. They mediate chemo-electrical signal transduction and its regulation via interconversion between multiple conformations that exhibit distinct pharmacological properties and biological activities. The large diversity of functional properties and the widely pleiotropic phenotypes, which arise from point mutations in their subunits (or from subunit substitutions), are interpreted in terms of an allosteric model that incorporates multiple discrete conformational states. The model predicts that three main categories of phenotypes may result from point mutations, altering selectively one (or more) of the following features: (i) the properties of individual binding sites (K phenotype), (ii) the biological activity of the ion channel (gamma phenotype) of individual conformations, or (iii) the isomerization constants between receptor conformations (L phenotype). Several nicotinic acetylcholine and glycine receptor mutants with complex phenotypes are quantitatively analyzed in terms of the model, and the analogies among phenotypes are discussed.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Développement d'un biosenseur BRET permettant le criblage de drogues qui causent l'activation de canaux Kir3 via les récepteurs couplés aux protéines G

Les récepteurs couplés aux protéines G forment des complexes multimériques comprenant protéines G et effecteurs. Nous cherchons à caractériser de tels complexes comprenant les récepteurs opioïdes delta (DOR) et les canaux Kir3, qui nous sont d’intérêt vu leur implication dans l’analgésie des opioïdes. Des expériences d’immunopurification, de BRET et de liaison GTPgS ont été réalisées à l’intérieur de cellules HEK293 transfectées. Les canaux Kir3 ont été co-immunopurifiés avec les DOR, suggérant une interaction spontanée entre récepteur et effecteur. Des essais BRET ont corroboré que l’interaction était présente dans des cellules vivantes et nous ont permis d’identifier une interaction spontanée et spécifique entre DOR/Gg et Gg/Kir3, indiquant leur coexistence en un même complexe. Puisque l’activation du récepteur implique la présence de changements conformationnels à l’intérieur de celui-ci, nous étions intéressés à vérifier si l’information conformationnelle circule à partir du récepteur lié au ligand jusqu’à l’effecteur en aval. Ainsi, nous avons déterminé l’effet de différents ligands sur le signal BRET généré par les paires suivantes : DOR/Gbg, DOR/Kir3 et Kir3/Gbg. Nous avons constaté une modulation de l’interaction DOR/Gbg et Gbg/Kir3 suivant l’ordre d’efficacité des ligands à stimuler la protéine G, ce que nous n’avons pas observé entre DOR et Kir3. Donc, nous concluons que l’information conformationnelle circule du récepteur au canal Kir3 via la protéine Gbg. Ces résultats nous ont permis de développer un biosenseur BRET (EYFP-Gg2/Kir3.1-Rluc) qui pourrait être utilisé dans le criblage à haut débit afin de détecter de nouvelles molécules ayant une grande efficacité à activer les canaux Kir3.

Crystal structure of Schistosoma purine nucleoside phosphorylase complexed with a novel monocyclic inhibitor

A novel inhibitor of Schistosoma PNP was identified using an ""in silico"" approach allied to enzyme inhibition assays. The compound has a monocyclic structure which has not been previously described for PNP inhibitors The crystallographic structure of the complex was determined and used to elucidate the binding mode within the active site Furthermore, the predicted pose was very similar to that determined crystallographically, validating the methodology The compound Sm_VS1, despite its low molecular weight, possesses an IC(50) of 1 3 mu M, surprisingly low when compared with purine analogues This is presumably due to the formation of eight hydrogen bonds with key residues in the active site E203, N245 and T244. The results of this study highlight the importance of the use of multiple conformations for the target during virtual screening. Indeed the Sm_VS1 compound was only identified after flipping the N245 side chain It is expected that the structure will be of use in the development of new highly active non-purine based compounds against the Sclustosoma enzyme. (c) 2010 Elsevier B V. All rights reserved

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

The role of Sse1 in the de novo formation and variant determination of the [PSI+] prion.

Yeast prions are a group of non-Mendelian genetic elements transmitted as altered and self-propagating conformations. Extensive studies in the last decade have provided valuable information on the mechanisms responsible for yeast prion propagation. How yeast prions are formed de novo and what cellular factors are required for determining prion "strains" or variants--a single polypeptide capable of existing in multiple conformations to result in distinct heritable phenotypes--continue to defy our understanding. We report here that Sse1, the yeast ortholog of the mammalian heat-shock protein 110 (Hsp110) and a nucleotide exchange factor for Hsp70 proteins, plays an important role in regulating [PSI+] de novo formation and variant determination. Overproduction of the Sse1 chaperone dramatically enhanced [PSI+] formation whereas deletion of SSE1 severely inhibited it. Only an unstable weak [PSI+] variant was formed in SSE1 disrupted cells whereas [PSI+] variants ranging from very strong to very weak were formed in isogenic wild-type cells under identical conditions. Thus, Sse1 is essential for the generation of multiple [PSI+] variants. Mutational analysis further demonstrated that the physical association of Sse1 with Hsp70 but not the ATP hydrolysis activity of Sse1 is required for the formation of multiple [PSI+] variants. Our findings establish a novel role for Sse1 in [PSI+] de novo formation and variant determination, implying that the mammalian Hsp110 may likewise be involved in the etiology of protein-folding diseases.

Structural flexibility in transcription complex formation revealed by protein–DNA photocrosslinking

The Oct-1 POU domain binds diverse DNA-sequence elements and forms a higher-order regulatory complex with the herpes simplex virus coregulator VP16. The POU domain contains two separate DNA-binding domains joined by a flexible linker. By protein–DNA photocrosslinking we show that the relative positioning of the two POU DNA-binding domains on DNA varies depending on the nature of the DNA target. On a single VP16-responsive element, the POU domain adopts multiple conformations. To determine the structure of the Oct-1 POU domain in a multiprotein complex with VP16, we allowed VP16 to interact with previously crosslinked POU-domain–DNA complexes and found that VP16 can associate with multiple POU-domain conformations. These results reveal the dynamic potential of a DNA-binding domain in directing transcriptional regulatory complex formation.

Multiple four-stranded conformations of human telomere sequence d(CCCTAA) in solution

By detailed NMR analysis of a human telomere repeating unit, d(CCCTAA), we have found that three distinct tetramers, each of which consists of four symmetric single-strands, slowly exchange in a slightly acidic solution. Our new finding is a novel i-motif topology (T-form) where T4 is intercalated between C1 and C2 of the other duplex. The other two tetramers have a topology where C1 is intercalated between C2 and C3 of the other parallel duplex, resulting in the non-stacking T4 residues (R-form), and a topology where C1 is stacked between C3 and T4 of the other duplex (S-form). From the NMR denaturation profile, the R-form is the most stable of the three structures in the temperature range of 15–50°C, the S-form the second and the T-form the least stable. The thermodynamic parameters indicate that the T-form is the most enthalpically driven and entropically opposed, and its population is increased with decreasing temperature. The T-form structure determined by restrained molecular dynamics calculation suggests that inter-strand van der Waals contacts in the narrow grooves should contribute to the enthalpic stabilization of the T-form.

Distinct conformations of Ly49 natural killer cell receptors mediate MHC class I recognition in trans and cis.

Certain cell-surface receptors engage ligands expressed on juxtaposed cells and ligands on the same cell. The structural basis for trans versus cis binding is not known. Here, we showed that Ly49 natural killer (NK) cell receptors bound two MHC class I (MHC-I) molecules in trans when the two ligand-binding domains were backfolded onto the long stalk region. In contrast, dissociation of the ligand-binding domains from the stalk and their reorientation relative to the NK cell membrane allowed monovalent binding of MHC-I in cis. The distinct conformations (backfolded and extended) define the structural basis for cis-trans binding by Ly49 receptors and explain the divergent functional consequences of cis versus trans interactions. Further analyses identified specific stalk segments that were not required for MHC-I binding in trans but were essential for inhibitory receptor function. These data identify multiple distinct roles of stalk regions for receptor function.

Assessing the Reproducibility of Clustering of Molecular Dynamics Conformations on Self-Organizing Maps

The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.

A conditional probability analysis of cystic fibrosis transmembrane conductance regulator gating indicates that ATP has multiple effects during the gating cycle

ATP-binding cassette (ABC) transporters bind and hydrolyze ATP. In the cystic fibrosis transmembrane conductance regulator Cl− channel, this interaction with ATP generates a gating cycle between a closed (C) and two open (O1 and O2) conformations. To understand better how ATP controls channel activity, we examined gating transitions from the C to the O1 and O2 states and from these open states to the C conformation. We made three main observations. First, we found that the channel can open into either the O1 or O2 state, that the frequency of transitions to both states was increased by ATP concentration, and that ATP increased the relative proportion of openings into O1 vs. O2. These results indicate that ATP can interact with the closed state to open the channel in at least two ways, which may involve binding to nucleotide-binding domains (NBDs) NBD1 and NBD2. Second, ATP prolonged the burst duration and altered the way in which the channel closed. These data suggest that ATP also interacts with the open channel. Third, the channel showed runs of specific types of open–closed transitions. This finding suggests a mechanism with more than one cycle of gating transitions. These data suggest models to explain how ATP influences conformational transitions in cystic fibrosis transmembrane conductance regulator and perhaps other ABC transporters.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

Conformations and folding of lysozyme ions in vacuo.

Proton transfer reactivity of isolated charge states of the protein hen egg-white lysozyme shows that multiple distinct conformations of this protein are stable in the gas phase. The reactivities of the 9+ and 10+ charge state ions, formed by electrospray ionization of "native" (disulfide-intact) and "denatured" (disulfide-reduced) solutions, are consistent with values calculated for ions in their crystal structure and fully denatured conformations, respectively. Charge states below 8+ of both forms, formed by proton stripping, have similar or indistinguishable reactivities, indicating that the disulfide-reduced ions fold in the gas phase to a more compact conformation.

Significance of bound water to local chain conformations in protein crystals.

We examine how the polypeptide chain in protein crystal structures exploits the multivalent hydrogen-bonding potential of bound water molecules. This shows that multiple interactions with a single water molecule tend to occur locally along the chain. A distinctive internal-coordinate representation of the local water-binding segments reveals several consensus conformations. The fractional water occupancy of each was found by comparison of the total number of conformations in the database regardless of the presence or absence of bound water. The water molecule appears particularly frequently in type II beta-turn geometries and an N-terminal helix feature. This work constitutes a first step into assessing not only the generality but also the significance of specific water binding in globular proteins.

Self-condensation of a thiazole-peptide bearing a 21-membered loop into a library of giant macrocycles with multiple orthogonal loops

Tetrapeptide analogue H-[Glu-Ser-Lys(Thz)]-OH, containing a turn-inducing thiazole constraint, was used as a template to produce a 21-membered structurally characterized loop by linking Glu and Lys side chains with a Val-Ile dipeptide. This template was oligomerized in one pot to a library (cyclo-[1](n), n = 2-10) of giant symmetrical macrocycles (up to 120-membered rings), fused to 2-10 appended loops that were carried intact through multiple oligomerization (chain extension) and cyclization (chain terminating) reactions of the template. A three-dimensional solution structure for cyclo-[1](3) shows all three appended loops projecting from the same face of the macrocycle. This is a promising approach to separating pepticle motifs over large distances.