Phosphatidylinositol 4,5-bisphosphate clusters act as molecular beacons for vesicle recruitment.

Synaptic-vesicle exocytosis is mediated by the vesicular Ca(2+) sensor synaptotagmin-1. Synaptotagmin-1 interacts with the SNARE protein syntaxin-1A and acidic phospholipids such as phosphatidylinositol 4,5-bisphosphate (PIP2). However, it is unclear how these interactions contribute to triggering membrane fusion. Using PC12 cells from Rattus norvegicus and artificial supported bilayers, we show that synaptotagmin-1 interacts with the polybasic linker region of syntaxin-1A independent of Ca(2+) through PIP2. This interaction allows both Ca(2+)-binding sites of synaptotagmin-1 to bind to phosphatidylserine in the vesicle membrane upon Ca(2+) triggering. We determined the crystal structure of the C2B domain of synaptotagmin-1 bound to phosphoserine, allowing development of a high-resolution model of synaptotagmin bridging two different membranes. Our results suggest that PIP2 clusters organized by syntaxin-1 act as molecular beacons for vesicle docking, with the subsequent Ca(2+) influx bringing the vesicle membrane close enough for membrane fusion.

Phytohormone collaboration: zooming in on auxin-brassinosteroid interactions.

Similar to animal hormones, classic plant hormones are small organic molecules that regulate physiological and developmental processes. In development, this often involves the regulation of growth through the control of cell size or division. The plant hormones auxin and brassinosteroid modulate both cell expansion and proliferation and are known for their overlapping activities in physiological assays. Recent molecular genetic analyses in the model plant Arabidopsis suggest that this reflects interdependent and often synergistic action of the two hormone pathways. Such pathway interactions probably occur through the combinatorial regulation of common target genes by auxin- and brassinosteroid-controlled transcription factors. Moreover, auxin and brassinosteroid signaling and biosynthesis and auxin transport might be linked by an emerging upstream connection involving calcium-calmodulin and phosphoinositide signaling.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Structural analysis of TCR-ligand interactions studied on H-2Kd-restricted cloned CTL specific for a photoreactive peptide derivative.

To study the interaction of the TCR with its ligand, the complex of a MHC molecule and an antigenic peptide, we modified a TCR contact residue of a H-2Kd-restricted antigenic peptide with photoreactive 4-azidobenzoic acid. The photoreactive group was a critical component of the epitope recognized by CTL clones derived from mice immunized with such a peptide derivative. The majority of these clones expressed V beta 1-encoded beta chains that were paired with J alpha TA28-encoded alpha chains. For one of these TCR, the photoaffinity labeled sites were mapped on the alpha chain as a J alpha TA28-encoded tryptophan and on the beta chain as a residue of the C' strand of V beta 1. Molecular modeling of this TCR suggested the presence of a hydrophobic pocket that harbors this tryptophan as well as a tyrosine on the C' strand of V beta 1 between which the photoreactive side chain inserts. It is concluded that this avid binding principle may account for the preferential selection of V beta 1 and J alpha TA28-encoded TCR.

Functional interactions of peroxisome proliferator-activated receptor, retinoid-X receptor, and Sp1 in the transcriptional regulation of the acyl-coenzyme-A oxidase promoter.

Peroxisome proliferator-activated receptor (PPARs) are members of the nuclear receptor superfamily. For transcriptional activation of their target genes, PPARs heterodimerize with the retinoid-X receptor (RXR). The convergence of the PPAR and RXR signaling pathways has been shown to have an important function in lipid metabolism. The promoter of the gene encoding the acyl-coenzyme-A oxidase (ACO), the rate-limiting enzyme in peroxisomal beta-oxidation of fatty acids, is a target site of PPAR action. In this study, we examined the role and the contribution of both cis-and trans-acting factors in the transcriptional regulation of this gene using transient transfections in insect cells. We identified several functional cis-acting elements present in the promoter of the ACO gene and established that PPAR-dependent as well as PPAR-independent mechanisms can activate the ACO promoter in these cells. We show that the PPAR/RXR heterodimer exerts its effect through two response elements within the ACO promoter, in synergy with the transcription factor Sp1 via five Sp1-binding sites. Furthermore, this functional interaction also occurs when Sp1 is co-expressed with PPAR or RXR alone, indicating that activation can occur independently of PPAR/RXR heterodimers.

Probing the interactions of NK cell receptors with ligand expressed in trans and cis.

Certain receptors on natural killer (NK) cells, which are specific for MHC class I (MHC-I) molecules, do not only interact with ligand expressed on opposing cell membranes (in trans) but also interact with those on the same cell membrane (in cis). Cis interactions have been demonstrated for only a small number of cell surface receptors. However, this has not been tested systematically, raising the possibility that additional receptors may be able to bind ligand expressed in cis. Here we describe a number of approaches to evaluate trans and cis binding of the Ly49A NK cell receptor to its H-2D(d) ligand. These procedures should facilitate the investigation of cis/trans interactions of other receptor-ligand pairs and simplify the analysis of NK cell receptor variants.

NO rebinding to myoglobin: a reactive molecular dynamics study.

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.

Integrating receptor interactions and dynamics and Interference with endocrine disruptors in the mechanisms of action of PPAR nuclear receptors

Abstract Peroxisome Proliferator-Activated Receptors (PPARs) form a family of three nuclear receptors regulating important cellular and metabolic functions. PPARs control gene expression by directly binding to target promoters as heterodimers with the Retinoid X Receptor (RXR), and their transcriptional activity is enhanced upon activation by natural or pharmacological ligands. The binding of PPAR/RXR heterodimers on target promoters allows the anchoring of a series of coactivators and corepressors involved in promoter remodeling and the recruitment of the transcription machinery. The transcriptional output finally depends on a complex interplay between (i) the respective expression levels of PPARs, RXRs and of other nuclear receptors competing for DNA binding and RXR recruitment, (ii) the availability and the nature of PPAR and RXR ligands, (iii) the expression levels and the nature of the different coactivators and corepressors and (iv) the sequence and the epigenetic status of the promoter. Understanding how all these factors and signals integrate and fine-tune transcription remains a challenge but is necessary to understand the specificity of the physiological functions regulated by PPARs. The work presented herein focuses on the molecular mechanisms of PPAR action and aims at understanding how the interactions and mobility of the receptor modulate transcription in the physiological context of a living cell: Such observations in vivo rely on the use of engineered fluorescent protein chimeras and require the development and the application of complementary imaging techniques such as Fluorescence Recovery After Photobleaching (FRAP), Fluorescence Resonance Energy Transfer (FRET) and Fluorescence Correlation Spectroscopy (FCS). Using such techniques, PPARs are shown to reside solely in the nucleus where they are constitutively associated with RXR but transcriptional activation by ligand binding -does not promote the formation of sub-nuclear structures as observed with other nuclear receptors. In addition, the engagement of unliganded PPARs in large complexes of cofactors in living cells provides a molecular basis for their ligand-independent activity. Ligand binding reduces receptor diffusion by promoting the recruitment of coactivators which further enlarge the size of PPAR complexes to acquire full transcriptional competence. Using these molecular approaches, we deciphered the molecular mechanisms through which phthalates, a class of pollutants from the plastic industry, interfere with PPARγ signaling. Mono-ethyl-hexyl-phthalate (MEHP) binding induces the recruitment of a specific subset of cofactors and translates into the expression of a specific subset of target genes, the transcriptional output being strongly conditioned by the differentiation status of the cell. This selective PPARγ modulation induces limited adipogenic effects in cellular models while exposure to phthalates in animal models leads to protective effects on glucose tolerance and diet-induced obesity. These results demonstrate that phthalates influence lipid and carbohydrate metabolism through complex mechanisms which most likely involve PPARγ but also probably PPARα and PPARß, Altogether, the molecular and physiological demonstration of the interference of pollutants with PPAR action outlines an important role of chemical exposure in metabolic regulations. Résumé Les PPARs (Peroxisome Proliferator-Activated Receptors) forment une famille de récepteurs nucléaires qui régulent des fonctions cellulaires et métaboliques importantes. Les PPARs contrôlent l'expression des gènes en se liant directement à leurs promoteurs sous forme d'hétérodimères avec les récepteurs RXR (Retinoid X Receptor), et leur activité transcriptionnelle est stimulée par la liaison de ligands naturels ou pharmacologiques. L'association des hétérodimères PPAR/RXR avec les promoteurs des gènes cibles permet le recrutement de coactivateurs et de corépresseurs qui vont permettre le remodelage de la chromatine et le recrutement de la machinerie transcriptionnelle. Les actions transcriptionnelles du récepteur dépendent toutefois d'interactions complexes qui sont régulées par (i) le niveau d'expression des PPARs, des RXRs et d'autres récepteurs nucléaires entrant en compétition pour la liaison à l'ADN et l'association avec RXR, (ii) la disponibilité et la nature de ligands de PPAR et de RXR, (iii) les niveaux d'expression et la nature des différents coactivateurs et corépresseurs et (iv) la séquence et le marquage épigénétique des promoteurs. La compréhension des mécanismes qui permettent d'intégrer ces aspects pour assurer une régulation fine de l'activité transcriptionnelle est un défi qu'il est nécessaire de relever pour comprendre la spécificité des fonctions physiologiques régulées par les PPARs. Ce travail concerne l'étude des mécanismes d'action moléculaire des PPARs et vise à mieux comprendre comment les interactions du récepteur avec d'autres protéines ainsi que la mobilité de ce dernier régulent son activité transcriptionnelle dans le contexte physiologique des cellules vivantes. De telles observations reposent sur l'emploi de protéines fusionnées à des protéines fluorescentes ainsi que sur le développement et l'utilisation de techniques d'imagerie complémentaires telles que le FRAP (Fluorescence Recovery After Photobleaching), le FRET (Fluorescence Resonance Energy Transfer) ou la FCS (Fluorescence Corrélation Spectroscopy). En appliquant ces méthodes, nous avons pu montrer que les PPARs résident toujours dans le noyau où ils sont associés de manière constitutive à RXR, mais que l'ajout de ligand n'induit pas la formation de structures sub-nucléaires comme cela a pu être décrit pour d'autres récepteurs nucléaires. De plus, les PPARs sont engagés dans de larges complexes protéiques de cofacteurs en absence de ligand, ce qui procure une explication moléculaire à leur activité ligand-indépendante. La liaison du ligand réduit la vitesse de diffusion du récepteur en induisant le recrutement de coactivateurs qui augmente encore plus la taille des complexes afin d'acquérir un potentiel d'activation maximal. En utilisant ces approches moléculaires, nous avons pu caractériser les mécanismes permettant aux phtalates, une classe de polluants provenant de l'industrie plastique, d'interférer avec PPARγ. La liaison du mono-ethyl-hexyl-phtalate (NERF) à PPARγ induit un recrutement sélectif de cofacteurs, se traduisant par l'induction spécifique d'un sous-ensemble de gènes qui varie en fonction du niveau de différentiation cellulaire. La modulation sélective de PPARγ par le MEHP provoque une adipogenèse modérée dans des modèles cellulaires alors que l'exposition de modèles animaux aux phtalates induit des effets bénéfiques sur la tolérance au glucose et sur le développement de l'obésité. Toutefois, les phtalates ont une action complexe sur le métabolisme glucido-lipidique en faisant intervenir PPARγ mais aussi probablement PPARα et PPARß. Cette démonstration moléculaire et physiologique de l'interférence des polluants avec les récepteurs nucléaires PPAR souligne un rôle important de l'exposition à de tels composés dans les régulations métaboliques.

Structure-based drug design studies of the interactions ofent-kaurane diterpenes derived from Wedelia paludosa with the Plasmodium falciparumsarco/endoplasmic reticulum Ca2+-ATPase PfATP6

Malaria is responsible for more deaths around the world than any other parasitic disease. Due to the emergence of strains that are resistant to the current chemotherapeutic antimalarial arsenal, the search for new antimalarial drugs remains urgent though hampered by a lack of knowledge regarding the molecular mechanisms of artemisinin resistance. Semisynthetic compounds derived from diterpenes from the medicinal plant Wedelia paludosawere tested in silico against the Plasmodium falciparumCa2+-ATPase, PfATP6. This protein was constructed by comparative modelling using the three-dimensional structure of a homologous protein, 1IWO, as a scaffold. Compound 21 showed the best docking scores, indicating a better interaction with PfATP6 than that of thapsigargin, the natural inhibitor. Inhibition of PfATP6 by diterpene compounds could promote a change in calcium homeostasis, leading to parasite death. These data suggest PfATP6 as a potential target for the antimalarial ent-kaurane diterpenes.

Molecular characterization of Unc5CL/ZUD

Résumé : Au cours de l'évolution, les organismes multicellulaires ont développé le système immunitaire afin de pouvoir se défendre contre les pathogènes tel que les bactéries, les virus, et les parasites. La réponse immunitaire doit être finement régulée par différentes voies de signalisation moléculaire, afin d'assurer une efficacité optimale, et d'éviter des dommages tissulaires indésirables. Les résultats expérimentaux décrits dans ce manuscrit, mettent en évidence que la protéine Unc5CL, qui contient un death domain (DD), est impliquée dans la régulation de la réponse immunitaire des muqueuses. Il a été démontré que cette protéine contient aussi un domaine transmembranaire de type III dans sa partie N-terminale, permettant ainsi de l'ancrer et d'exposer sa partie C-terminale dans le cytosol, un prérequis pour la signalisation dans ce compartiment cellulaire. De plus, cette protéine a la capacité d'activer le facteur de transcription NFxB, qui joue un rôle important dans le système immunitaire, ainsi que dans d'autres processus cellulaires essentiels. Le profil transcriptionnel révèle que l'activation de NF-κB induite par Unc5CL conduit principalement à une réponse inflammatoire, qui se caractérise par la production de diverses chimiokines (e.g. CXCL-1, IL-8 et CCL20). Il a également été démontré que Unc5CL requiert les mêmes molécules qui sont utilisées dans la voie de signalisation des récepteurs de la famille toll et de l'interleukine-1. De manière similaire à leur protéine adaptatrice MyD88, Unc5CL a la capacité de recruter, via une interaction homotypique DD-DD, les kinases IRAK1 et IRAK4 qui contiennent elles aussi un DD, permettant ainsi au signal d'être transmis. La production d'un anticorps polyclonal contre le DD de Unc5CL a permis d'identifier des lignées cellulaires et des tissus exprimant cette protéine, ainsi que de déterminer sa localisation sub-cellulaire. Unc5CL a été détecté dans les cellules de la muqueuse utérine et intestinale, ainsi que dans une lignée cellulaire issue d'un adénocarcinome colorectal humain, les CaCo-2. Dans chacun de ces cas, Unc5CL a été principalement détectée au niveau apical des cellules épithéliales polarisées. De manière similaire à PIDD, une protéine impliquée dans la réponse aux dommages à l'ADN, et au constituant des pores nucléaires Nup98, Unc5CL est constitutivement clivé de manière autoprotéolytique, au niveau d'un site HFS. Il est intéressant d'observer que les deux fragments ainsi générés restent fortement associés l'un à l'autre après clivage. Finalement, un criblage protéomique pour identifier un partenaire d'interaction, a mis en évidence l'ubiquitin ligase E3 ITCH, qui régule de manière négative Unc5CL en augmentant sa dégradation. Summary : Multicellular organisms have evolved the immune system in order to defend themselves against pathogens such as bacteria, viruses and eukaryotic parasites. Immune responses have to be tightly orchestrated by signaling mechanisms to achieve optimal effectiveness and minimal tissue damage. The experimental results in this thesis manuscript provide evidence that the death domain (DD)-containing protein Unc5CL might be involved in the regulation of mucosal immune responses. It could be shown that the protein contains an N-terminal type-III transmembrane domain that anchors the protein with its C-terminus exposed to the cytosol, a prerequisite for signaling events in this compartment. Furthermore, the protein has the capacity to activate the transcription factor NF-κB, which plays an important role in the immune system as well as in other essential cellular processes. Transcriptional profiling revealed that Unc5CL-mediated activation of NF-κB mainly leads to an inflammatory response, characterized by the production of chemokines (e.g. CXCL-l, IL-8 and CCL20). Furthermore, it could be shown that Unc5CL requires the same downstream signaling molecules as the evolutionarily ancient tolUinterleukin-1 receptor family. Similar to their adapter protein MyD88, Unc5CL has the capacity to recruit the DD-containing kinases IRAKI and IRAK4 for signaling and can interact with these proteins via homotypic DD-DD interactions. Generation of polyclonal antibodies raised against the DD of Unc5CL allowed the identification of cell lines and tissues that express the endogenous protein as well as to confine its subcellular localization. Unc5CL was detected in primary mucosal uterine and intestinal epithelial cells as well as in the human colorectal adenocarcinoma cell line CaCo-2. In all cases, the protein was mainly localized to the apical face of these polarized epithelial cells. Similar to PIDD, a protein critically involved in responses to DNA damage, and the nuclear pore component Nup98, Unc5CL is constitutively autoproteolytically processed at an HFS site. Interestingly, the two generated cleavage fragments remain tightly associated after processing. Finally, a proteomics screen for interaction partners identified the E3 ubiquitin ligase ITCH as a negative regulator of Unc5CL by targeting the protein for degradation.

The relevance of the Laplacian of intracule and extracule density distributions for analyzing electron-electron interactions in molecules

A topological analysis of intracule and extracule densities and their Laplacians computed within the Hartree-Fock approximation is presented. The analysis of the density distributions reveals that among all possible electron-electron interactions in atoms and between atoms in molecules only very few are located rigorously as local maxima. In contrast, they are clearly identified as local minima in the topology of Laplacian maps. The conceptually different interpretation of intracule and extracule maps is also discussed in detail. An application example to the C2H2, C2H4, and C2H6 series of molecules is presented

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Development and validation of AMANDA, a new algorithm for selecting highly relevant regions in Molecular Interaction Fields

Descriptors based on Molecular Interaction Fields (MIF) are highly suitable for drug discovery, but their size (thousands of variables) often limits their application in practice. Here we describe a simple and fast computational method that extracts from a MIF a handful of highly informative points (hot spots) which summarize the most relevant information. The method was specifically developed for drug discovery, is fast, and does not require human supervision, being suitable for its application on very large series of compounds. The quality of the results has been tested by running the method on the ligand structure of a large number of ligand-receptor complexes and then comparing the position of the selected hot spots with actual atoms of the receptor. As an additional test, the hot spots obtained with the novel method were used to obtain GRIND-like molecular descriptors which were compared with the original GRIND. In both cases the results show that the novel method is highly suitable for describing ligand-receptor interactions and compares favorably with other state-of-the-art methods.

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.