NR2E3 mutations in enhanced S-cone sensitivity syndrome (ESCS), Goldmann-Favre syndrome (GFS), clumped pigmentary retinal degeneration (CPRD), and retinitis pigmentosa (RP).

NR2E3, also called photoreceptor-specific nuclear receptor (PNR), is a transcription factor of the nuclear hormone receptor superfamily whose expression is uniquely restricted to photoreceptors. There, its physiological activity is essential for proper rod and cone photoreceptor development and maintenance. Thirty-two different mutations in NR2E3 have been identified in either homozygous or compound heterozygous state in the recessively inherited enhanced S-cone sensitivity syndrome (ESCS), Goldmann-Favre syndrome (GFS), and clumped pigmentary retinal degeneration (CPRD). The clinical phenotype common to all these patients is night blindness, rudimental or absent rod function, and hyperfunction of the "blue" S-cones. A single p.G56R mutation is inherited in a dominant manner and causes retinitis pigmentosa (RP). We have established a new locus-specific database for NR2E3 (www.LOVD.nl/eye), containing all reported mutations, polymorphisms, and unclassified sequence variants, including novel ones. A high proportion of mutations are located in the evolutionarily-conserved DNA-binding domains (DBDs) and ligand-binding domains (LBDs) of NR2E3. Based on homology modeling of these NR2E3 domains, we propose a structural localization of mutated residues. The high variability of clinical phenotypes observed in patients affected by NR2E3-linked retinal degenerations may be caused by different disease mechanisms, including absence of DNA-binding, altered interactions with transcriptional coregulators, and differential activity of modifier genes.

Theoretical study on receptor-G protein recognition: new insights into the mechanisms of the alpha1b-adrenergic recptor activation

This work compares the structural/dynamics features of the wild-type alb-adrenergic receptor (AR) with those of the D142A active mutant and the agonist-bound state. The two active receptor forms were compared in their isolated states as well as in their ability to form homodimers and to recognize the G alpha q beta 1 gamma 2 heterotrimer. The analysis of the isolated structures revealed that, although the mutation- and agonist-induced active states of the alpha 1b-AR are different, they, however, share several structural peculiarities including (a) the release of some constraining interactions found in the wild-type receptor and (b) the opening of a cytosolic crevice formed by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. Accordingly, also their tendency to form homodimers shows commonalties and differences. In fact, in both the active receptor forms, helix 6 plays a crucial role in mediating homodimerization. However, the homodimeric models result from different interhelical assemblies. On the same line of evidence, in both of the active receptor forms, the cytosolic opened crevice recognizes similar domains on the G protein. However, the docking solutions are differently populated and the receptor-G protein preorientation models suggest that the final complexes should be characterized by different interaction patterns.

The fourth transmembrane segment of the Na,K-ATPase alpha subunit: a systematic mutagenesis study.

The Na,K-ATPase is a major ion-motive ATPase of the P-type family responsible for many aspects of cellular homeostasis. To determine the structure of the pathway for cations across the transmembrane portion of the Na,K-ATPase, we mutated 24 residues of the fourth transmembrane segment into cysteine and studied their function and accessibility by exposure to the sulfhydryl reagent 2-aminoethyl-methanethiosulfonate. Accessibility was also examined after treatment with palytoxin, which transforms the Na,K-pump into a cation channel. Of the 24 tested cysteine mutants, seven had no or a much reduced transport function. In particular cysteine mutants of the highly conserved "PEG" motif had a strongly reduced activity. However, most of the non-functional mutants could still be transformed by palytoxin as well as all of the functional mutants. Accessibility, determined as a 2-aminoethyl-methanethiosulfonate-induced reduction of the transport activity or as inhibition of the membrane conductance after palytoxin treatment, was observed for the following positions: Phe(323), Ile(322), Gly(326), Ala(330), Pro(333), Glu(334), and Gly(335). In accordance with a structural model of the Na,K-ATPase obtained by homology modeling with the two published structures of sarcoplasmic and endoplasmic reticulum calcium ATPase (Protein Data Bank codes 1EUL and 1IWO), the results suggest the presence of a cation pathway along the side of the fourth transmembrane segment that faces the space between transmembrane segments 5 and 6. The phenylalanine residue in position 323 has a critical position at the outer mouth of the cation pathway. The residues thought to form the cation binding site II ((333)PEGL) are also part of the accessible wall of the cation pathway opened by palytoxin through the Na,K-pump.

Transcriptional Activation and Receptor Dimerization is Affected by Mutations in the NR2E3 Ligand-Binding Domain

Purpose: Mutations in the ligand-binding domain (LBD) of NR2E3 cause recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS), Goldmann-Favre syndrome (GFS) and clumped pigmentary retinal degeneration (CPRD). In addition to ligand binding, the LBD contains also essential amino acid sequences for the oligomerization of nuclear receptors. The aim of our studies is to characterize the impact of mutations in the LBD on receptor oligomerization and transcriptional activity of NR2E3. Methods: The different NR2E3 mutants were generated by QuickChange mutagenesis and analyzed in 293T-based transactivation studies and BRET2 (bioluminescence resonance electron transfer) assays. In silico homology modeling of mutant proteins was also performed using available crystallographic data of related nuclear receptors. Results: The mutants p.W234S, p.A256V, p.A256E, p.L263P, p.R309G, p.R311Q, p.R334G, p.L336P, p.L353V, p.R385P and p.M407K, all located in the LBD, showed impaired receptor dimerization at various degrees. Impaired repressor dimerization as assessed by BRET2 assays did not always correlate with impaired repressor function of NR2E3 as assessed by cell-based reporter assays. There were minor differences of transcriptional activity of mutant proteins on mouse S-opsin (opn1sw), mouse cone arrestin (arr3) and human cone arrestin, suggesting that the effect of LBD mutations was independent of the promoter context. Conclusions: Mutational analysis and homology modeling allowed the characterization of potential oligomerization interfaces of the NR2E3 LBD. Additionally, mutations in NR2E3 LBD may cause recessive retinal degenerations by different molecular mechanisms.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Local alignment refinement using structural assessment.

Homology modeling is the most commonly used technique to build a three-dimensional model for a protein sequence. It heavily relies on the quality of the sequence alignment between the protein to model and related proteins with a known three dimensional structure. Alignment quality can be assessed according to the physico-chemical properties of the three dimensional models it produces.In this work, we introduce fifteen predictors designed to evaluate the properties of the models obtained for various alignments. They consist of an energy value obtained from different force fields (CHARMM, ProsaII or ANOLEA) computed on residue selected around misaligned regions. These predictors were evaluated on ten challenging test cases. For each target, all possible ungapped alignments are generated and their corresponding models are computed and evaluated.The best predictor, retrieving the structural alignment for 9 out of 10 test cases, is based on the ANOLEA atomistic mean force potential and takes into account residues around misaligned secondary structure elements. The performance of the other predictors is significantly lower. This work shows that substantial improvement in local alignments can be obtained by careful assessment of the local structure of the resulting models.

Differential Dimerization of Variants Linked to Enhanced S-Cone Sensitivity Syndrome (ESCS) Located in the NR2E3 Ligand-Binding Domain.

NR2E3 encodes the photoreceptor-specific nuclear hormone receptor that acts as a repressor of cone-specific gene expression in rod photoreceptors, and as an activator of several rod-specific genes. Recessive variants located in the ligand-binding domain (LBD) of NR2E3 cause enhanced short wavelength sensitive- (S-) cone syndrome (ESCS), a retinal degeneration characterized by an excess of S-cones and non-functional rods. We analyzed the dimerization properties of NR2E3 and the effect of disease-causing LBD missense variants by bioluminescence resonance energy transfer (BRET(2) ) protein interaction assays. Homodimerization was not affected in presence of p.A256V, p.R039G, p.R311Q, and p.R334G variants, but abolished in presence of p.L263P, p.L336P, p.L353V, p.R385P, and p.M407K variants. Homology modeling predicted structural changes induced by NR2E3 LBD variants. NR2E3 LBD variants did not affect interaction with CRX, but with NRL and rev-erbα/NR1D1. CRX and NRL heterodimerized more efficiently together, than did either with NR2E3. NR2E3 did not heterodimerize with TLX/NR2E1 and RXRα/NR2C1. The identification of a new compound heterozygous patient with detectable rod function, who expressed solely the p.A256V variant protein, suggests a correlation between LBD variants able to form functional NR2E3 dimers and atypical mild forms of ESCS with residual rod function.

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D&#8201;=&#8201;&#8764;1&#8201;-&#8201;5&#8201;&#956;M. Beyond the affinity threshold at K D&#8201;<&#8201;1&#8201;&#956;M we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

Modelagem de proteínas por homologia

Modeling methods to derive 3D-structure of proteins have been recently developed. Protein homology-modeling, also known as comparative protein modeling, is nowadays the most accurate protein modeling method. This technique can produce useful models for about an order of magnitude more protein sequences than there have been structures determined by experiment in the same amount of time. All current protein homology-modeling methods consist of four sequential steps: fold assignment and template selection, template-target alignment, model building, and model evaluation. In this paper we discuss in some detail the protein-homology paradigm, its predictive power and its limitations.

In silico analysis identifies a C3HC4-RING finger domain of a putative E3 ubiquitin-protein ligase located at the C-terminus of a polyglutamine-containing protein

Almost identical polyglutamine-containing proteins with unknown structures have been found in human, mouse and rat genomes (GenBank AJ277365, AF525300, AY879229). We infer that an identical new gene (RING) finger domain of real interest is located in each C-terminal segment. A three-dimensional (3-D) model was generated by remote homology modeling and the functional implications are discussed. The model consists of 65 residues from terminal position 707 to 772 of the human protein with a total length of 796 residues. The 3-D model predicts a ubiquitin-protein ligase (E3) as a binding site for ubiquitin-conjugating enzyme (E2). Both enzymes are part of the ubiquitin pathway to label unwanted proteins for subsequent enzymatic degradation. The molecular contact specificities are suggested for both the substrate recognition and the residues at the possible E2-binding surface. The predicted structure, of a ubiquitin-protein ligase (E3, enzyme class number 6.3.2.19, CATH code 3.30.40.10.4) may contribute to explain the process of ubiquitination. The 3-D model supports the idea of a C3HC4-RING finger with a partially new pattern. The putative E2-binding site is formed by a shallow hydrophobic groove on the surface adjacent to the helix and one zinc finger (L722, C739, P740, P741, R744). Solvent-exposed hydrophobic amino acids lie around both zinc fingers (I717, L722, F738, or P765, L766, V767, V733, P734). The 3-D structure was deposited in the protein databank theoretical model repository (2B9G, RCSB Protein Data Bank, NJ).

Rôle des facteurs physico-chimiques du micro-environnement intestinal et des boucles inter-hélicales du Domaine I dans l’activité de la toxine insecticide Cry9Ca du bacille de Thuringe

Une fois ingérées par un insecte sensible, les toxines insecticides du bacille de Thuringe doivent être activées par les protéases intestinales de cet insecte. Leur premier domaine, un ensemble de sept hélices-α amphipathiques, est responsable de leur insertion dans la membrane luminale de certaines cellules de l’intestin médian, ce qui crée des pores peu sélectifs. La toxicité et la capacité à former des pores d’une telle toxine, la Cry9Ca, de ses mutants simples R164A et R164K et d’un fragment de 55 kDa résultant d’un clivage protéolytique au niveau de son résidu 164 ont été étudiées à l’aide d’une combinaison de modélisation par homologie, de bioessais, d’expériences de gonflement osmotique avec des vésicules de membrane en bordure en brosse de larves de sphinx du tabac et de mesures électrophysiologiques sur des intestins isolés. Ni les mutations simples ni le clivage protéolytique n’ont altéré la toxicité de la Cry9Ca. Dans une solution à faible force ionique, toutefois, la formation des pores dépend fortement du pH : une augmentation de celui-ci de 6,5 à 10,5 a entraîné une baisse irrégulière et par étapes successives de la perméabilité membranaire. Les quatre préparations de toxine ont néanmoins dépolarisé la membrane apicale d’intestins médians fraîchement isolés baignant dans une solution contenant 122 mM de KCl à pH 10,5. L’activité de la Cry9Ca, et des mutants R164A et R164K, a été grandement stimulée lorsque les expériences ont été effectuées en présence de suc intestinal, de lipides extraits d’un volume équivalent de suc intestinal ou d’un cocktail d’inhibiteurs de protéases solubles dans l’eau. De plus, le rôle des boucles inter-hélicales du Domaine I lors de l’insertion dans la membrane a été étudié avec des mutants doubles de la Cry9Ca dont les mutations introduisaient, neutralisaient ou renversaient une charge électrique. À l’exception de trois d’entres eux, tous ces mutants ont conservé une toxicité et une capacité à former des pores comparables à celles de la toxine parentale. L’ensemble de ces résultats suggère que le micro-environnement de l’intestin médian contribue à minimiser l’influence des charges de surface portées par les résidus des boucles inter-hélicales du Domaine I sur la capacité des toxines du bacille de Thuringe à former des pores. Il indique aussi que, d’une part, selon le site de clivage et les conditions expérimentales utilisées, des protéolyses supplémentaires de la toxine Cry9Ca activée peuvent soit stimuler, soit nuire à son activité et que, d’autre part, le suc intestinal du sphinx du tabac contient probablement un inhibiteur de protéases qui pourrait jouer un rôle important dans l’activité des toxines du bacille de Thuringe.

Calmodulin/KCa3.1 channel interactions as determinant to the KCa3.1 Ca2+ dependent gating : theoretical and experimental analyses

Differentes études ont montré que la sensibilité au Ca2+ du canal KCa3.1, un canal potassique indépendant du voltage, était conférée par la protéine calmoduline (CaM) liée de façon constitutive au canal. Cette liaison impliquerait la région C-lobe de la CaM et un domaine de $\ikca$ directement relié au segment transmembranaire S6 du canal. La CaM pourrait égalment se lier au canal de façon Ca2+ dépendante via une interaction entre un domaine de KCa3.1 du C-terminal (CaMBD2) et la région N-lobe de la CaM. Une étude fut entreprise afin de déterminer la nature des résidus responsables de la liaison entre le domaine CaMBD2 de KCa3.1 et la région N-lobe de la CaM et leur rôle dans le processus d'ouverture du canal par le Ca2+. Une structure 3D du complexe KCa3.1/CaM a d'abord été générée par modélisation par homologie avec le logiciel MODELLER en utilisant comme référence la structure cristalline du complexe SK2.2/CaM (PDB: 1G4Y). Le modèle ainsi obtenu de KCa3.1 plus CaM prévoit que le segment L361-S372 dans KCa3.1 devrait être responsable de la liaison dépendante du Ca2+ du canal avec la région N-lobe de la CaM via les résidus L361 et Q364 de KCa3.1 et E45, E47 et D50 de la CaM. Pour tester ce modèle, les résidus dans le segment L361-S372 ont été mutés en Cys et l'action du MTSET+ (chargé positivement) et MTSACE (neutre) a été mesurée sur l'activité du canal. Des enregistrements en patch clamp en configuration ``inside-out`` ont montré que la liaison du réactif chargé MTSET+ au le mutant Q364C entraîne une forte augmentation du courant, un effet non observé avec le MTSACE. De plus les mutations E45A et E47A dans la CaM, ont empêché l'augmentation du courant initié par MTSET+ sur le mutant Q364C. Une analyse en canal unitaire a confirmé que la liaison MTSET+ à Q364C cause une augmentation de la probabilité d'ouverture de KCa3.1 par une déstabilisation de l'état fermé du canal. Nous concluons que nos résultats sont compatibles avec la formation de liaisons ioniques entre les complexes chargés positivement Cys-MTSET+ à la position 364 de KCa3.1 et les résidus chargés négativement E45 et E47 dans la CaM. Ces données confirment qu'une stabilisation électrostatique des interactions CaM/KCa3.1 peut conduire à une augmentation de la probabilité d'ouverture du canal en conditions de concentrations saturantes de Ca2+.

Modulation de l’adressage membranaire et de la fonction du canal CaV2.3 par les résidus leucine du domaine guanylate kinase impliqués dans la liaison à forte affinité de CaVβ

Les canaux Ca2+ activés par le voltage (CaV) sont des protéines membranaires qui génèrent des courants Ca2+ dans les cellules excitables suite à une dépolarisation membranaire. Ces complexes oligomériques sont classifiés selon les propriétés structurelles de la sous-unité principale qui forme le pore du canal, soit la sous-unité CaVα1. La sous-unité auxiliaire CaVβ module l’expression membranaire et la dépendance au voltage du « gating » de la sous-unité CaVα1 des canaux HVA (« high-voltage-activated ») CaV1 et CaV2. La sous-unité CaVβ est formée par un domaine SH3 (« Src homology-3 ») connecté à un domaine GK (« guanylate kinase-like ») par le biais d’un domaine variable HOOK. Dans le but d’identifier les résidus dans la CaVβ3 qui sont responsables de la densité membranaire du CaV2.3, nous avons produit des mutants de la sous-unité auxiliaire le long de ses domaines fonctionnels. Cela dit, la délétion complète du domaine SH3 ainsi que la délétion du domaine HOOK n’ont pas modifié la densité membranaire de CaV2.3 ni ses propriétés d’activation. Cependant, la délétion de cinq résidus dans le domaine GK interrompt l’expression membranaire et l’expression fonctionnelle de CaV2.3. La mutation de résidus identifiés précédemment comme soutenant une affinité de liaison de l’ordre du nanomolaire dans le domaine GK de CaVβ n’a pas modifié de manière significative l’adressage membranaire de CaV2.3. Toutefois, les mutations de quatre résidus leucine dans les régions α3, α6, β10 et α9 du domaine GK ont grandement réduit l’adressage membranaire du canal CaV2.3. Nos résultats confirment que le domaine GK contient les déterminants moléculaires responsables de la fonction chaperone de CaVβ. Cela dit, l’adressage membranaire induit par CaVβ semble être déterminé par des éléments structuraux qui ne sont pas strictement dépendants d’une liaison à haute affinité de CaVβ sur CaVα1.

Propriétés électrophysiologiques des canaux ioniques formés par la toxine nématicide Cry5Ba du bacille de Thuringe dans les bicouches lipidiques planes

Les toxines Cry sont des protéines synthétisées sous forme de cristaux par la bactérie bacille de Thuringe pendant la sporulation. Elles sont largement utilisées comme agents de lutte biologique, car elles sont toxiques envers plusieurs espèces d’invertébrées, y compris les nématodes. Les toxines Cry5B sont actives contre certaines espèces de nématodes parasites, y compris Ankylostoma ceylanicum un parasite qui infeste le système gastro-intestinal des humains. Jusqu’au présent, le mode d’action des toxines Cry nématicides reste grandement inconnu, sauf que leurs récepteurs spécifiques sont des glycolipides et qu’elles causent des dommages importants aux cellules intestinales. Dans cette étude, on démontre pour la première fois que la toxine nématicide Cry5Ba, membre de la famille des toxines à trois domaines et produite par la bactérie bacille de Thuringe, forme des pores dans les bicouches lipidiques planes en absence de récepteurs. Les pores formés par cette toxine sont de sélectivité cationique, à pH acide ou alcalin. Les conductances des pores formés sous conditions symétriques de 150 mM de KCl varient entre 17 et 330 pS, à pH 6.0 et 9.0. Les niveaux des conductances les plus fréquemment observés diffèrent les uns des autres par environ 17 à 18 pS, ce qui est compatible avec l’existence d’arrangement d’un nombre différent de pores élémentaires similaires, activés de façon synchronisée, ou avec la présence d’oligomères de tailles variables et de différents diamètres de pores.

Catalytic mechanism of inulinase from Arthrobacter sp S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k(cat)/K-m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k(cat), but not due to variations in K-m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK(a) of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK(a). Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme. (C) 2008 Elsevier Inc. All rights reserved.