Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Design of potent inhibitors of human RAD51 recombinase based on BRC motifs of BRCA2 protein: modeling and experimental validation of a chimera peptide.

We have previously shown that a 28-amino acid peptide derived from the BRC4 motif of BRCA2 tumor suppressor inhibits selectively human RAD51 recombinase (HsRad51). With the aim of designing better inhibitors for cancer treatment, we combined an in silico docking approach with in vitro biochemical testing to construct a highly efficient chimera peptide from eight existing human BRC motifs. We built a molecular model of all BRC motifs complexed with HsRad51 based on the crystal structure of the BRC4 motif-HsRad51 complex, computed the interaction energy of each residue in each BRC motif, and selected the best amino acid residue at each binding position. This analysis enabled us to propose four amino acid substitutions in the BRC4 motif. Three of these increased the inhibitory effect in vitro, and this effect was found to be additive. We thus obtained a peptide that is about 10 times more efficient in inhibiting HsRad51-ssDNA complex formation than the original peptide.

Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

Interaction between calcium intake and menarcheal age on bone mass gain: an eight-year follow-up study from prepuberty to postmenarche.

Both late menarcheal age and low calcium intake (Ca intake) during growth are risk factors for osteoporosis, probably by impairing peak bone mass. We investigated whether lasting gain in areal bone mineral density (aBMD) in response to increased Ca intake varies according to menarcheal age and, conversely, whether Ca intake could influence menarcheal age. In an initial study, 144 prepubertal girls were randomized in a double-blind controlled trial to receive either a Ca supplement (Ca-suppl.) of 850 mg/d or placebo from age 7.9-8.9 yr. Mean aBMD gain determined by dual energy x-ray absorptiometry at six sites (radius metaphysis, radius diaphysis, femoral neck, trochanter, femoral diaphysis, and L2-L4) was significantly (P = 0.004) greater in the Ca-suppl. than in the placebo group (27 vs. 21 mg/cm(2)). In 122 girls followed up, menarcheal age was recorded, and aBMD was determined at 16.4 yr of age. Menarcheal age was lower in the Ca-suppl. than in the placebo group (P = 0.048). Menarcheal age and Ca intake were negatively correlated (r = -0.35; P < 0.001), as were aBMD gains from age 7.9-16.4 yr and menarcheal age at all skeletal sites (range: r = -0.41 to r = -0.22; P < 0.001 to P = 0.016). The positive effect of Ca-suppl. on the mean aBMD gain from baseline remained significantly greater in girls below, but not in those above, the median of menarcheal age (13.0 yr). Early menarcheal age (12.1 +/- 0.5 yr): placebo, 286 +/- 36 mg/cm(2); Ca-suppl., 317 +/- 46 (P = 0.009); late menarcheal age (13.9 +/- 0.5 yr): placebo, 284 +/- 58; Ca-suppl., 276 +/- 50 (P > 0.05). The level of Ca intake during prepuberty may influence the timing of menarche, which, in turn, could influence long-term bone mass gain in response to Ca supplementation. Thus, both determinants of early menarcheal age and high Ca intake may positively interact on bone mineral mass accrual.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Human topoisomerase II-DNA interaction study by using atomic force microscopy.

Type II topoisomerases (Topo II) are unique enzymes that change the DNA topology by catalyzing the passage of two double-strands across each other by using the energy from ATP hydrolysis. In vitro, human Topo II relaxes positive supercoiled DNA around 10-fold faster than negative supercoiled DNA. By using atomic force microscopy (AFM) we found that human Topo II binds preferentially to DNA cross-overs. Around 50% of the DNA crossings, where Topo II was bound to, presented an angle in the range of 80-90°, suggesting a favored binding geometry in the chiral discrimination by Topo II. Our studies with AFM also helped us visualize the dynamics of the unknotting action of Topo II in knotted molecules.

Molecular basis of interaction between Na,K-ATPase and FXYD proteins.

Résumé au large public Notre corps est constitué de différents types de cellules. La condition minimale ou primordiale pour la survie des cellules est d'avoir de l'énergie. Cette tâche est assumée en partie par une protéine qui se situe dans la membrane de chaque cellule. Nommé Na, K¬ATPase ou pompe à sodium, c'est une protéine pressente dans toutes les cellules chez les mammifères est composée de deux sous-unités, α et β. En transportant 3 ions de sodium hors de la cellule et 2 ions de potassium à l'intérieur de la cellule, elle transforme l'énergie chimique sous forme de l'ATP en énergie motrice, qui permet aux cellules par la suite d'échanger des matériaux entre l'espace intracellulaire et extracellulaire ainsi que d'ingérer des nutriments provenant de son environnement. Le manque de cette protéine chez la souris entraîne la mort de l'embryon. Des défauts fonctionnels de cette protéine sont responsables de plusieurs maladies humaines comme par exemple, un type de migraine. En dehors de sa fonction vitale, cette protéine est également engagée dans diverses activités physiologiques comme la contractilité musculaire, l'activité nerveuse et la régulation du volume sanguin. Vue l'importance de cette protéine, sa découverte par Jens C. Skou en 1957 a été honorée d'un Prix Noble de chimie quarante ans plus tard. Depuis lors, nous connaissons de mieux en mieux les mécanismes de fonctionnement de la Na, K-ATPase. Entre autre, sa régulation par une famille de protéines appelées protéines FXYD. Cette famille contient 7 membres (FXYD 1-7). L'un d'entre eux nommé FXYD 2 est lié à une maladie héréditaire connue sous le nom de hypomagnesemia. Nous disposons actuellement d'informations concernant les conséquences de la régulation par les protéines FXYD sur activité de la Na, K-ATPase, mais nous savons très peu sur le mode d'interaction entre les protéines FXYD et la Na, K-ATPase. Dans ce travail de thèse, nous avons réussi à localiser des zones d'interaction dans la sous- unité a de la Na, K-ATPase et dans FXYD 7. En même temps, nous avons déterminé un 3ème site de liaison spécifique au sodium de la Na, K-ATPase. Une partie de ce site se situe à l'intérieur d'un domaine protéique qui interagit avec les protéines FXYD. De plus, ce site a été démontré comme responsable d'un mécanisme de transport de la Na, K-ATPase caractérisé par un influx ionique. En conclusion, les résultats de ce travail de thèse fournissent de nouvelles preuves sur les régions d'interaction entre la Na, K-ATPase et les protéines FXYD. La détermination d'un 3ème site spécifique au sodium et sa relation avec un influx ionique offrent la possibilité 1) d'explorer les mécanismes avec lesquels les protéines FXYD régulent l'activité de la Na, ATPase et 2) de localiser un site à sodium qui est essentielle pour mieux comprendre l'organisation et le fonctionnement de la Na, K-ATPase. Résumé Les gradients de concentration de Na+ et de K+ à travers la membrane plasmatique des cellules animales sont cruciaux pour la survie et l'homéostasie de cellules. De plus, des fonctions cellulaires spécifiques telles que la reabsorption de Na dans le rein et le côlon, la contraction musculaire et l'excitabilité nerveuse dépendent de ces gradients. La Na, K¬ATPase ou pompe à sodium est une protéine membranaire ubiquitaire. Elle crée et maintient ces gradients en utilisant l'énergie obtenu par l'hydrolyse de l'adénosine triphosphate. L'unité fonctionnelle minimale de cette protéine se compose d'une sous-unité catalytique α et d'une sous-unité régulatrice β. Récemment, il a été montré que des membres de la famille FXYD, sont des régulateurs tissu-spécifiques de la Na, K-ATPase qui influencent ses propriétés de transport. Cependant, on connaît peu de chose au sujet de la nature moléculaire de l'interaction entre les protéines FXYD et la Na, K-ATPase. Dans cette étude, nous fournissons, pour la première fois, l'évidence directe que des résidus du domaine transmembranaire (TM) 9 de la sous-unité α de la Na, K-ATPase sont impliqués dans l'interaction fonctionnelle et structurale avec les protéines FXYD. De plus nous avons identifié des régions dans le domaine transmembranaire de FXYD 7 qui sont importantes pour l'association stable avec la Na, K-ATPase et une série de résidus responsables des régulations fonctionnelles. Nous avons aussi montré les contributions fonctionnelles du TM 9 de la Na, K-ATPase à la translocation de Na + en déterminant un 3ème site spécifique au Na+. Ce site se situe probablement dans un espace entre TM 9, TM 6 et TM 5 de la sous-unité α de la pompe à sodium. De plus, nous avons constaté que le 3ème site de Na + est fonctionnellement lié à un courant entrant de la pompe sensible à l'ouabaïne et activé par le pH acide. En conclusion, ce travail donne de nouvelles perspectives de l'interaction structurale et fonctionnelle entre les protéines FXYD et la Na, K-ATPase. En outre, les contributions fonctionnelles de TM 9 offrent de nouvelles possibilités pour explorer le mécanisme par lequel les protéines FXYD régulent les propriétés fonctionnelles de la Na, K-ATPase. La détermination du 3ème site au Na + fournit une compréhension avancée du site spécifique au Na + de la Na, K-ATPase et du mécanisme de transport de la Na, K-ATPase. Summary The Na+ and K+ gradients across the plasma membrane of animal cells are crucial for cell survival and homeostasis. Moreover, specific tissue functions such as Na+ reabsorption in kidney and colon, muscle contraction and nerve excitability depend on the maintenance of these gradients. Na, K-ATPase or sodium pump, an ubiquitous membrane protein, creates and maintains these gradients by using the energy from the hydrolysis of ATP. The minimal functional unit of this protein is composed of a catalytic α subunit and a regulatory β subunit. Recently, members of the FXYD family, have been reported to be tissue-specific regulators of Na, K-ATPase by influencing its transport properties. However, little is known about the molecular nature of the interaction between FXYD proteins and Na, K-ATPase. In this study, we provide, for the first time, direct evidence that residues from the transmembrane (TM) domain 9 of the α subunit of Na, K-ATPase are implicated in the functional and structural interaction with FXYD proteins. Moreover, we have identified regions in the TM domain of FXYD 7 important for the stable association with Na, K-ATPase and a stretch of residues responsible for the functional regulations. We have further revealed the functional contributions of TM 9 of the Na, K-ATPase α subunit to the Na+ translocation by determining a 3rd Na+-specific cation binding site. This site is likely in a space between TM 9, TM 6 and TM 5 of the a subunit of the sodium pump. Moreover, we have found that the 3rd Na+ binding site is functionally linked to an acidic pH- activated ouabain-sensitive inward pump current. In conclusion, this work gives new insights into the structural and functional interaction between FXYD proteins and Na, K-ATPase. Functional contributions of TM 9 offer new possibilities to explore the mechanism by which FXYD proteins regulate functional properties of Na, K-ATPase. The determination of the 3rd Na+ binding site provides an advanced understanding concerning the Na+ -specific binding site of Na, K-ATPase and the 3rd Na+ site related transport mechanism.

Horizontal, but not vertical, biotic interactions affect fine-scale plant distribution patterns in a low-energy system

Studies of species range determinants have traditionally focused on abiotic variables (typically climatic conditions), and therefore the recent explicit consideration of biotic interactions represents an important advance in the field. While these studies clearly support the role of biotic interactions in shaping species distributions, most examine only the influence of a single species and/or a single interaction, failing to account for species being subject to multiple concurrent interactions. By fitting species distribution models (SDMs), we examine the influence of multiple vertical (i.e., grazing, trampling, and manuring by mammalian herbivores) and horizontal (i.e., competition and facilitation; estimated from the cover of dominant plant species) interspecific interactions on the occurrence and cover of 41 alpine tundra plant species. Adding plant-plant interactions to baseline SDMs (using five field-quantified abiotic variables) significantly improved models' predictive power for independent data, while herbivore-related variables had only a weak influence. Overall, abiotic variables had the strongest individual contributions to the distribution of alpine tundra plants, with the importance of horizontal interaction variables exceeding that of vertical interaction variables. These results were consistent across three modeling techniques, for both species occurrence and cover, demonstrating the pattern to be robust. Thus, the explicit consideration of multiple biotic interactions reveals that plant-plant interactions exert control over the fine-scale distribution of vascular species that is comparable to abiotic drivers and considerably stronger than herbivores in this low-energy system.

A model for physical activity, bahaviour and energy expenditure assessment in real-life condition

RésuméL'origine de l'obésité, qui atteint des proportions épidémiques, est complexe. Elle est liée au mode de vie et au comportement des individus par rapport à l'activité physique, expression des choix individuels et de l'interaction avec l'environnement. Les mesures du comportement au niveau de l'activité physique des individus face à leur environnement, la répartition des types d'activité physique, la durée, la fréquence, l'intensité, et la dépense énergétique sont d'une grande importance. Aujourd'hui, il y a un manque de méthodes permettant une évaluation précise et objective de l'activité physique et du comportement des individus. Afin de compléter les recherches relatives à l'activité physique, à l'obésité et à certaines maladies, le premier objectif du travail de thèse était de développer un modèle pour l'identification objective des types d'activité physique dans des conditions de vie réelles et l'estimation de la dépense énergétique basée sur une combinaison de 2 accéléromètres et 1 GPS. Le modèle prend en compte qu'une activité donnée peut être accomplie de différentes façons dans la vie réelle. Les activités quotidiennes ont pu être classées en 8 catégories, de sédentaires à actives, avec une précision de 1 min. La dépense énergétique a pu peut être prédite avec précision par le modèle. Après validation du modèle, le comportement des individus de l'activité physique a été évalué dans une seconde étude. Nous avons émis l'hypothèse que, dans un environnement caractérisé par les pentes, les personnes obèses sont tentées d'éviter les pentes raides et de diminuer la vitesse de marche au cours d'une activité physique spontanée, ainsi que pendant les exercices prescrits et structurés. Nous avons donc caractérisé, par moyen du modèle développé, le comportement des individus obèses dans un environnement vallonné urbain. La façon dont on aborde un environnement valloné dans les déplacements quotidiens devrait également être considérée lors de la prescription de marche supplémentaire afin d'augmenter l'activité physique.SummaryOrigin of obesity, that reached epidemic proportion, is complex and may be linked to different lifestyle and physical activity behaviour. Measurement of physical activity behaviour of individuals towards their environment, the distribution of physical activity in terms of physical activity type, volume, duration, frequency, intensity, and energy expenditure is of great importance. Nowadays, there is a lack of methods for accurate and objective assessment of physical activity and of individuals' physical activity behaviour. In order to complement the research relating physical activity to obesity and related diseases, the first aim of the thesis work was to develop a model for objective identification of physical activity types in real-life condition and energy expenditure based on a combination of 2 accelerometers and 1 GPS device. The model takes into account that a given activity can be achieved in many different ways in real life condition. Daily activities could be classified in 8 categories, as sedentary to active physical activity, within 1 min accuracy, and physical activity patterns determined. The energy expenditure could be predicted accurately with an accuracy below 10%. Furthermore, individuals' physical activity behaviour is expression of individual choices and their interaction with the neighbourhood environment. In a second study, we hypothesized that, in an environment characterized by inclines, obese individuals are tempted to avoid steep positive slopes and to decrease walking speed during spontaneous outdoor physical activity, as well as during prescribed structured bouts of exercise. Finally, we characterized, by mean of the developed model, the physical activity behaviour of obese individuals in a hilly urban environment. Quantifying how one tackles hilly environment or avoids slope in their everyday displacements should be also considered while prescribing extra walking in free-living conditions in order to increase physical activity.

Rise in plasma lactate concentrations with psychosocial stress: a possible sign of cerebral energy demand.

OBJECTIVE: It is known that exogenous lactate given as an i.v. energy infusion is able to counteract a neuroglycopenic state that developed during psychosocial stress. It is unknown, however, whether the brain under stressful conditions can induce a rise in plasma lactate to satisfy its increased needs during stress. Since lactate is i) an alternative cerebral energy substrate to glucose and ii) its plasmatic concentration is influenced by the sympathetic nervous system, the present study aimed at investigating whether plasma lactate concentrations increase with psychosocial stress in humans. METHODS: 30 healthy young men participated in two sessions (stress induced by the Trier Social Stress Test and a non-stress control session). Blood samples were frequently taken to assess plasma lactate concentrations and stress hormone profiles. RESULTS: Plasma lactate increased 47% during psychosocial stress (from 0.9 ± 0.05 to 1.4 ± 0.1 mmol/l; interaction time × stress intervention: F = 19.7, p < 0.001). This increase in lactate concentrations during stress was associated with an increase in epinephrine (R(2) = 0.221, p = 0.02) and ACTH concentrations (R(2) = 0.460, p < 0.001). CONCLUSION: Plasma lactate concentrations increase during acute psychosocial stress in humans. This finding suggests the existence of a demand mechanism that functions to allocate an additional source of energy from the body towards the brain, which we refer to as 'cerebral lactate demand'.

Visual-gustatory interaction: orbitofrontal and insular cortices mediate the effect of high-calorie visual food cues on taste pleasantness.

Vision provides a primary sensory input for food perception. It raises expectations on taste and nutritional value and drives acceptance or rejection. So far, the impact of visual food cues varying in energy content on subsequent taste integration remains unexplored. Using electrical neuroimaging, we assessed whether high- and low-calorie food cues differentially influence the brain processing and perception of a subsequent neutral electric taste. When viewing high-calorie food images, participants reported the subsequent taste to be more pleasant than when low-calorie food images preceded the identical taste. Moreover, the taste-evoked neural activity was stronger in the bilateral insula and the adjacent frontal operculum (FOP) within 100 ms after taste onset when preceded by high- versus low-calorie cues. A similar pattern evolved in the anterior cingulate (ACC) and medial orbitofrontal cortex (OFC) around 180 ms, as well as, in the right insula, around 360 ms. The activation differences in the OFC correlated positively with changes in taste pleasantness, a finding that is an accord with the role of the OFC in the hedonic evaluation of taste. Later activation differences in the right insula likely indicate revaluation of interoceptive taste awareness. Our findings reveal previously unknown mechanisms of cross-modal, visual-gustatory, sensory interactions underlying food evaluation.

Mutations in the DNA-binding domain of NR2E3 affect in vivo dimerization and interaction with CRX.

BACKGROUND: NR2E3 (PNR) is an orphan nuclear receptor essential for proper photoreceptor determination and differentiation. In humans, mutations in NR2E3 have been associated with the recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS) and, more recently, with autosomal dominant retinitis pigmentosa (adRP). NR2E3 acts as a suppressor of the cone generation program in late mitotic retinal progenitor cells. In adult rod photoreceptors, NR2E3 represses cone-specific gene expression and acts in concert with the transcription factors CRX and NRL to activate rod-specific genes. NR2E3 and CRX have been shown to physically interact in vitro through their respective DNA-binding domains (DBD). The DBD also contributes to homo- and heterodimerization of nuclear receptors. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed NR2E3 homodimerization and NR2E3/CRX complex formation in an in vivo situation by Bioluminescence Resonance Energy Transfer (BRET(2)). NR2E3 wild-type protein formed homodimers in transiently transfected HEK293T cells. NR2E3 homodimerization was impaired in presence of disease-causing mutations in the DBD, except for the p.R76Q and p.R104W mutant proteins. Strikingly, the adRP-linked p.G56R mutant protein interacted with CRX with a similar efficiency to that of NR2E3 wild-type and p.R311Q proteins. In contrast, all other NR2E3 DBD-mutant proteins did not interact with CRX. The p.G56R mutant protein was also more effective in abolishing the potentiation of rhodospin gene transactivation by the NR2E3 wild-type protein. In addition, the p.G56R mutant enhanced the transrepression of the M- and S-opsin promoter, while all other NR2E3 DBD-mutants did not. CONCLUSIONS/SIGNIFICANCE: These results suggest different disease mechanisms in adRP- and ESCS-patients carrying NR2E3 mutations. Titration of CRX by the p.G56R mutant protein acting as a repressor in trans may account for the severe clinical phenotype in adRP patients.

MM-GBSA binding free energy decomposition and T cell receptor engineering.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by class I major histocompatibility complexes (MHC) is the key event in the immune response against virus infected cells or tumor cells. The major determinant of T cell activation is the affinity of the TCR for the peptide-MHC complex, though kinetic parameters are also important. A study of the 2C TCR/SIYR/H-2Kb system using a binding free energy decomposition (BFED) based on the MM-GBSA approach had been performed to assess the performance of the approach on this system. The results showed that the TCR-p-MHC BFED including entropic terms provides a detailed and reliable description of the energetics of the interaction (Zoete and Michielin, 2007). Based on these results, we have developed a new approach to design sequence modifications for a TCR recognizing the human leukocyte antigen (HLA)-A2 restricted tumor epitope NY-ESO-1. NY-ESO-1 is a cancer testis antigen expressed not only in melanoma, but also on several other types of cancers. It has been observed at high frequencies in melanoma patients with unusually positive clinical outcome and, therefore, represents an interesting target for adoptive transfer with modified TCR. Sequence modifications of TCR potentially increasing the affinity for this epitope have been proposed and tested in vitro. T cells expressing some of the proposed TCR mutants showed better T cell functionality, with improved killing of peptide-loaded T2 cells and better proliferative capacity compared to the wild type TCR expressing cells. These results open the door of rational TCR design for adoptive transfer cancer therapy.

The DNA-binding domain (DBD) is essential for NR2E3 dimerization and interaction with Crx

Purpose:NR2E3 (PNR) is an orphan nuclear receptor essential for proper photoreceptor determination and differentiation. In humans, mutations in NR2E3 have been associated with the recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS) and, more recently, with autosomal dominant retinitis pigmentosa (adRP). NR2E3 acts in concert with the transcription factors Crx and Nrl to repress cone-specific genes and activate rod-specific genes. NR2E3 and Crx have been shown to physically interact by their DNA-binding domain (DBD), which may also be implicated in the dimerization process of the nuclear receptor. However, neither NR2E3 homodimerization nor NR2E3/Crx complex formation has been investigated in detail. Methods:In this present work, we analyzed the dimerization of the NR2E3 protein and its interaction with Crx by bioluminescence resonance energy transfer (BRET2) which utilizes Renilla luciferase (hRluc) protein and its substrate DeepBlueC as an energy donor and a mutant green fluorescent protein (GFP2) as the acceptor. We investigated, on whole intact cells, the role of NR2E3 DBD-mutations in dimerization and association with Crx. Results:We clearly showed that NR2E3 formed homodimers in HEK-293T cells. Moreover, all causative NR2E3 mutations present in the DBD of the protein showed an alteration in dimerization, except for the R76Q and the R104W mutants. Interestingly, the adRP-linked G56R mutant was the only DBD-NR2E3 mutant that showed a correct interaction with Crx. Finally, we observed a decrease in rhodospin gene transactivation for all DBD-NR2E3 mutants tested and no potentiation for the adRP-linked G56R mutant. In addition, the p.G56R mutant enhanced the transrepression of M-opsin promoter, while all other DBD-NR2E3 mutants did not repress M-opsin transactivation. Conclusions:A defect, either in the dimer formation or in the interaction of NR2E3 with Crx, leads to abnormal transcriptional activity on rhodopsin and M-opsin promoter and to an atypical retinal development; while the titration of Crx by p.G56R-NR2E3 leads to low levels of rhodopsin and M-opsin expression and may be responsible for the strong adRP phenotype.