Molecular mechanism of myosin Va recruitment to dense core secretory granules.

The brain-spliced isoform of Myosin Va (BR-MyoVa) plays an important role in the transport of dense core secretory granules (SGs) to the plasma membrane in hormone and neuropeptide-producing cells. The molecular composition of the protein complex that recruits BR-MyoVa to SGs and regulates its function has not been identified to date. We have identified interaction between SG-associated proteins granuphilin-a/b (Gran-a/b), BR-MyoVa and Rab27a, a member of the Rab family of GTPases. Gran-a/b-BR-MyoVa interaction is direct, involves regions downstream of the Rab27-binding domain, and the C-terminal part of Gran-a determines exon specificity. MyoVa and Gran-a/b are partially colocalised on SGs and disruption of Gran-a/b-BR-MyoVa binding results in a perinuclear accumulation of SGs which augments nutrient-stimulated hormone secretion in pancreatic beta-cells. These results indicate the existence of at least another binding partner of BR-MyoVa that was identified as rabphilin-3A (Rph-3A). BR-MyoVa-Rph-3A interaction is also direct and enhanced when secretion is activated. The BR-MyoVa-Rph-3A and BR-MyoVa-Gran-a/b complexes are linked to a different subset of SGs, and simultaneous inhibition of these complexes nearly completely blocks stimulated hormone release. This study demonstrates that multiple binding partners of BR-MyoVa regulate SG transport, and this molecular mechanism is universally used by neuronal, endocrine and neuroendocrine cells.

Molecular insights for optimizing T cell receptor specificity against cancer.

Cytotoxic CD8 T cells mediate immunity to pathogens and they are able to eliminate malignant cells. Immunity to viruses and bacteria primarily involves CD8 T cells bearing high affinity T cell receptors (TCRs), which are specific to pathogen-derived (non-self) antigens. Given the thorough elimination of high affinity self/tumor-antigen reactive T cells by central and peripheral tolerance mechanisms, anti-cancer immunity mostly depends on TCRs with intermediate-to-low affinity for self-antigens. Because of this, a promising novel therapeutic approach to increase the efficacy of tumor-reactive T cells is to engineer their TCRs, with the aim to enhance their binding kinetics to pMHC complexes, or to directly manipulate the TCR-signaling cascades. Such manipulations require a detailed knowledge on how pMHC-TCR and co-receptors binding kinetics impact the T cell response. In this review, we present the current knowledge in this field. We discuss future challenges in identifying and targeting the molecular mechanisms to enhance the function of natural or TCR-affinity optimized T cells, and we provide perspectives for the development of protective anti-tumor T cell responses.

Identification, using a siRNA screen approach, and molecular characterization of a new Fas pathway modulator: STK11/

Abstract : Apoptosis is an evolutionarily conserved cellular suicide mechanism that can be triggered by activation of various pathways, such as the Fas-Pathway. Upon stimulation by its specific ligand (FasL), present at the surface of Cytotoxic Τ lymphocytes, the death receptor Fas initiates a signaling cascade culminating in the activation of cellular caspases, leading thus to cell death of the target cell (e.g. transformed cell). Dysregulation of apoptosis in general, and of Fas pathway in particular, was shown to contribute to pathogenesis of cancers and many human diseases. Even though, during the last decades the molecular mechanisms of apoptosis have been widely studied, it is important to better understand the mechanisms leading to apoptosis, to improve our understanding of pathological processes, and generate more subtle apoptosis-modulating therapies to fight cancer and other diseases. In order to identify new components of the Fas signaling pathway, a screen based on the mechanism of RNA interference was undertaken. After a first and a second manual whole-kinome screen, we identified several strong positive hits that showed a protection against Fas ligand-induced apoptosis with distinct siRNAs, notably STK11, an interesting tumor suppressor mutated in several sporadic and inherited cancers. The STK11 functional characterization reveals that this kinase represents an apically acting general pro-apoptotic modulator of the extrinsic pathway (FasL, TRAIL, TNF-induced apoptosis), but not of the intrinsic apoptotic pathway. The STK11 action on the Fas pathway was shown to be dependent on its kinase activity, but independent of AMPK, a well-characterized STK11 downstream substrate. Furthermore, STK11 was shown to interact with caspase-8, a major mediator of the extrinsic pathway, and modulate its activity through an unclear mechanism that may involve an STK11-dependant caspase-8 phosphorylation. This modification may allow a proper caspase-8 polyubiquitination and activation in p62 sequestosmes aggregates, but may also increase the activation of caspase-8 at the DISC level. In addition, we observed that STK11 modulate not only the apoptotic pathway induced by Fas engagement, but also FasL-induced JNK and NF- KB, sustaining an upstream role of this kinase in the pathway. In conclusion, our report reveals that STK11 is an important pro-apoptotic modulator of the Fas pathway in particular, and extrinsic pathway in general. Our finding could explain, at least partially, why inactivating mutations of the kinase leads to cancer, by allowing resistance to apoptosis and accordingly evasion of immune surveillance. Résumé : L'apoptose est un mécanisme de suicide cellulaire, conservé dans diverses espèces, et qui au niveau moléculaire est déclenché par différentes voies de signalisation, comme par exemple lors de l'activation du récepteur Fas. La liaison du ligand FasL au récepteur de la mort Fas, induit une cascade de signalisation qui conduit à l'activation des caspases. Les lymphocytes Τ cytotoxiques peuvent utiliser la voie Fas pour induire la mort et se débarrasser de cellules dangereuses pour le reste de l'organisme, tel que les cellules transformées. La dysrégulation de l'apoptose en général, et de la voie Fas en particulier, peut contribuer à diverses maladies telles que le cancer. Même si ces dernières décennies, les mécanismes moléculaires conduisant à l'apoptose ont été extensivement étudiés, il reste néanmoins important de mieux comprendre le phénomène d'apoptose, pour améliorer notre compréhension des processus pathologiques, mais surtout dans le but de développer de nouvelles thérapies ciblant l'apoptose contre le cancer et d'autres pathologies. Pour identifier de nouveau constituants de la voie Fas, un criblage génétique basé sur l'interférence à l'ARN a été entrepris. Après un premier et un deuxième criblage d'une librairie du kinome, nous avons identifié différentes protéines qui pourraient jouer un rôle positif dans la voie Fas, et en particulier la protéine suppresseur de tumeur STK11, qui est fréquemment mutée dans divers cancers sporadiques et héréditaires. La caractérisation fonctionnelle de STK11 a révélé que cette kinase était un modulateur apical de la voie extrinsèque de l'apoptose en général (Fas, TNF, TRAIL), mais pas de la voie intrinsèque. L'action de STK11 sur la voie Fas est dépendante de sa fonction kinase, mais indépendante de l'AMPK, un substrat bien caractérisé de STK11. De plus, STK11 interagît avec la caspase-8, un constituant majeur de la voie Fas, et module son activité, par un mécanisme encore peu clair qui pourrait impliquer une phosphorylation de la caspase-8 par STK11. Cette modification pourrait permettre une activation optimale de la caspase-8 en jouant un rôle dans le processus de polyubiquitination de la caspase-8, phénomène qui semble être important pour l'activation de la caspase-8 dans des agrégats protéiques avec p62, mais qui pourrait aussi augmenter son activation au niveau du DISC. Finalement, nous avons observé que STK11 modulait non seulement la voie apoptotique déclenchée par l'activation de Fas, mais aussi les voies non-apoptotiques de Fas, comme JNK et NF-KB. En conclusion notre étude, révèle que STK11 est un important modulateur pro- apoptotique de la voie Fas, et de la voie extrinsèque en général. Cette découverte pourrait expliquer, du moins partiellement, pourquoi les mutations inactivatrices de STK11 conduisent au cancer, par une augmentation de la résistance à l'apoptose et donc par l'évasion de la surveillance immunitaire.

Cytogenetic and molecular relationships between Zarudny's rock shrew (Crocidura zarudnyi; Mammalia : Soricomorpha) and Eurasian taxa

We karyotyped and sequenced 1,140 base pairs of the mitochondrial DNA cytochrome b of a specimen of Zarudny's rock shrew (Crocidura zarudnyi) from Baluchestan, southeastern Iran, to clarify its cytogenetic and molecular relationships with other Eurasian species of Crocidura. According to the karyotype (2N = 40, FN = 50), Zarudny's rock shrew belongs to the group of the lesser white-toothed shrew (C. suaveolens), which is different from other known crocidurine karyotypes, considering the combination of the diploid and fundamental number of chromosomes. Molecular results revealed that C. zarudnyi is included in a monophyletic clade with the C. suaveolens group, where it is a sister taxon to the others (mean Kimura 2-parameter distance = 9.7%).

The early Miocene rotation of the Corso-Sardinian block. New paleomagnetic constraints for the end of the motion

The paleomagnetic investigations carried out in the 70's on Oligo-Miocene volcanics of Sardinia have demonstrated that the island was turned by 35-30 degrees clockwise from 33 Ma up to 3-1-20.5 Ma and rotated counterclockwise in a few million years [De Jong et al., 1969, 1973; Bobier et Coulon, 1970; Coulon et al., 1974; Manzoni, 1974, 1975; Bellon rr nl.. 1977: Edel et Lortscher, 1977; Edel, 1979, 1980]. Since then, the end of the rotation fixed at 19 Ma by Montigny er al. [1981] was the subject of discussions and several studies associating paleomagnetism and radiometric dating were undertaken [Assorgia er al., 1994: Vigliotti et Langenheim, 1995: Deino et al., 1997; Gattacceca rt Deino, 1999]. This is a contribution to this debate that is hampered by thr important secular variation recorded in the volcanics. The only way to get our of this problem is to sample series of successive flows as completely as possible, and to reduce the effect of secular variation by the calculation of means. Sampling was performed north of Bonorva in 5 pyroclastic flows that belong to the upper ignimbritic series SI2 according to Coulon rr nl. [1974] or LBLS according to Assorgia et al. [1997] (fig. I). Ar-40/Ar-39 dating of biotites from the debris flow (MDF) has yielded an age or 18.35 +/- 0.03 Ma [Dubois, 2000]. Five of the investigated sites are located beneath the debris flow ITV, TVB, TVD, SPM85, SPM86), one site was cured in the matrix of the debris flow (MDF) and one in 4 metric blocks included in the flow (DFC). Another site was sampled in the upper ash flow (PDM) that marks the end of the pyroclastic activity, just before the marine transgression. According to micropaleontological and radiometric dating this transgression has occurred between 18.35 and 17.6 Ma [Dubois, 2000]. After removal of a soft viscous component, the thermal demagnetization generally shows a univectorial behaviour of the remanent magnetization (fig. 2a). The maximum unblocking temperatures of 580-620 degrees (tab. I) and a rapid saturation below 100 mT (fig. 3) indicate that the carrier of the characteristic magnetization is magnetite. The exception comes: from the upper site PDM in which were found two characteristic components, one with a normal polarity and low unblocking temperatures up to 350 degreesC and one with a reversed polarity and maximum unblocking temperatures at 580-600 degreesC of magnetite. After calculation of a mean direction for each flow, the mean << Al >> direction 4 degrees /57 degrees (alpha (95) = 13 degrees) computed with the mean directions for the 5 flows may be considered as weakly affected by secular variation. But the results require a more careful examination. The declinations are N to NNW beneath the debris flow. NNW in the debris flow. and NNE (or SSW) above the debris flow, The elongated distribution of the directions obtained at sites TVB and TVD. scattered from the mean direction of TV to the mean direction of MDF is interpreted as due to partial overprinting during the debris How volcanic episode, The low temperature component PDMa is likely related to the alteration seen on thin sections and is also viewed as an overprint. As NNE/SSW directions occur as well below (mean direction << B >> : 5 degrees /58 degrees) as above the debris flow (PDMb : 200 degrees/-58 degrees). the NNW directions (<< C >> : 337 degrees /64 degrees) associated with the debris flow volcanism may be interpreted as resulting from a magnetic field excursion. According to the polarity scale of Cande and Kent [1992, 1995] and the radiometric age of MDF, the directions with normal polarity (TV, TVB, TVD, SPM85. SPM86a. MDF. DFC) may represent the period 5En. while the directions with reversed polarity PDMb and SPM86b were likely acquired during the period 5Dr. Using the mean << Al >> direction, the mean << B >>, or the PDM direction (tab. I). the deviation in declination with the direction of stable Europe 6.4 degrees /58.7 degrees (alpha (95) = 8 degrees) for a selection of 4 middle Tertiary poles by Besse et Courtillot [1991] or 7 degrees /56 degrees (alpha (95) = 3 degrees) for 19 poles listed by Edel [1980] can be considered as negligible. Using the results from the uppermost ignimbritic layer of Anglona also emplaced around 18.3 Ma [Odin rt al.. 1994]. the mean direction << E >> (3 degrees /51.5 degrees) leads to the same conclusion. On the contrary, when taking into account all dated results available for the period 5En (mean direction << D >> 353 degrees /56 degrees for 45 sites) (tab. II). the deviation 13 degrees is much more significant. As the rotation of Sardinia started around 21-20.5 Ma. the assumption of a constant velocity of rotation and the deviations of the Sardinia directions with respect to the stable Europe direction locate the end of the motion between 18.3 and 17.2 or 16.7 Ma (fig. 4). During the interval 18.35-17.5 Ma, the marine transgression took place. At the same period a NE-SW shortening interpreted as resulting from the collision of Sardinia with Apulia affected different parts of the island [Letouzey et al., 1982]. Consequently, the new paleomagnetic results and the tectono-sedimentary evolution are in favour of an end of the rotation at 17.5-18 Ma.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.

Cellular and molecular mechanisms underlying the effects of sleep loss on hippocampal synaptic plasticity

SUMMARY : The function of sleep for the organism is one of the most persistent and perplexing questions in biology. Current findings lead to the conclusion that sleep is primarily for the brain. In particular, a role for sleep in cognitive aspects of brain function is supported by behavioral evidence both in humans and animals. However, in spite of remarkable advancement in the understanding of the mechanisms underlying sleep generation and regulation, it has been proven difficult to determine the neurobiological mechanisms underlying the beneficial effect of sleep, and the detrimental impact of sleep loss, on learning and memory processes. In my thesis, I present results that lead to several critical steps forward in the link between sleep and cognitive function. My major result is the molecular identification and physiological analysis of a protein, the NR2A subunit of NMDA receptor (NMDAR), that confers sensitivity to sleep loss to the hippocampus, a brain structure classically involved in mnemonic processes. Specifically, I used a novel behavioral approach to achieve sleep deprivation in adult C57BL6/J mice, yet minimizing the impact of secondary factors associated with the procedure,.such as stress. By using in vitro electrophysiological analysis, I show, for the first time, that sleep loss dramatically affects bidirectional plasticity at CA3 to CA1 synapses in the hippocampus, a well established cellular model of learning and memory. 4-6 hours of sleep loss elevate the modification threshold for bidirectional synaptic plasticity (MT), thereby promoting long-term depression of CA3 to CA 1 synaptic strength after stimulation in the theta frequency range (5 Hz), and rendering long-term potentiation induction.more difficult. Remarkably, 3 hours of recovery sleep, after the deprivation, reset the MT at control values, thus re-establishing the normal proneness of synapses to undergo long-term plastic changes. At the molecular level, these functional changes are paralleled by a change in the NMDAR subunit composition. In particular, the expression of the NR2A subunit protein of NMDAR at CA3 to CA1 synapses is selectively and rapidly increased by sleep deprivation, whereas recovery sleep reset NR2A synaptic content to control levels. By using an array of genetic, pharmacological and computational approaches, I demonstrate here an obligatory role for NR2A-containing NMDARs in conveying the effect of sleep loss on CA3 to CAl MT. Moreover, I show that a genetic deletion of the NR2A subunit fully preserves hippocampal plasticity from the impact of sleep loss, whereas it does not alter sleepwake behavior and homeostatic response to sleep deprivation. As to the mechanism underlying the effects of the NR2A subunit on hippocampal synaptic plasticity, I show that the increased NR2A expression after sleep loss distinctly affects the contribution of synaptic and more slowly recruited NMDAR pools activated during plasticity-induction protocols. This study represents a major step forward in understanding the mechanistic basis underlying sleep's role for the brain. By showing that sleep and sleep loss affect neuronal plasticity by regulating the expression and function of a synaptic neurotransmitter receptor, I propose that an important aspect of sleep function could consist in maintaining and regulating protein redistribution and ion channel trafficking at central synapses. These findings provide a novel starting point for investigations into the connections between sleep and learning, and they may open novel ways for pharmacological control over hippocampal .function during periods of sleep restriction. RÉSUMÉ DU PROJET La fonction du sommeil pour l'organisme est une des questions les plus persistantes et difficiles dans la biologie. Les découvertes actuelles mènent à la conclusion que le sommeil est essentiel pour le cerveau. En particulier, le rôle du sommeil dans les aspects cognitifs est soutenu par des études comportementales tant chez les humains que chez les animaux. Cependant, malgré l'avancement remarquable dans la compréhension des mécanismes sous-tendant la génération et la régulation du sommeil, les mécanismes neurobiologiques qui pourraient expliquer l'effet favorable du sommeil sur l'apprentissage et la mémoire ne sont pas encore clairs. Dans ma thèse, je présente des résultats qui aident à clarifier le lien entre le sommeil et la fonction cognitive. Mon résultat le plus significatif est l'identification moléculaire et l'analyse physiologique d'une protéine, la sous-unité NR2A du récepteur NMDA, qui rend l'hippocampe sensible à la perte de sommeil. Dans cette étude, nous avons utilisé une nouvelle approche expérimentale qui nous a permis d'induire une privation de sommeil chez les souris C57BL6/J adultes, en minimisant l'impact de facteurs confondants comme, par exemple, le stress. En utilisant les techniques de l'électrophysiologie in vitro, j'ai démontré, pour la première fois, que la perte de sommeil est responsable d'affecter radicalement la plasticité bidirectionnelle au niveau des synapses CA3-CA1 de l'hippocampe. Cela correspond à un mécanisme cellulaire de l'apprentissage et de la mémoire bien établi. En particulier, 4-6 heures de privation de sommeil élèvent le seuil de modification pour la plasticité synaptique bidirectionnelle (SM). Comme conséquence, la dépression à long terme de la transmission synaptique est induite par la stimulation des fibres afférentes dans la bande de fréquences thêta (5 Hz), alors que la potentialisation à long terme devient plus difficile. D'autre part, 3 heures de sommeil de récupération sont suffisant pour rétablir le SM aux valeurs contrôles. Au niveau moléculaire, les changements de la plasticité synaptiques sont associés à une altération de la composition du récepteur NMDA. En particulier, l'expression synaptique de la protéine NR2A du récepteur NMDA est rapidement augmentée de manière sélective par la privation de sommeil, alors que le sommeil de récupération rétablit l'expression de la protéine au niveau contrôle. En utilisant des approches génétiques, pharmacologiques et computationnelles, j'ai démontré que les récepteurs NMDA qui expriment la sous-unité NR2A sont responsables de l'effet de la privation de sommeil sur le SM. De plus, nous avons prouvé qu'une délétion génétique de la sous-unité NR2A préserve complètement la plasticité synaptique hippocampale de l'impact de la perte de sommeil, alors que cette manipulation ne change pas les mécanismes de régulation homéostatique du sommeil. En ce qui concerne les mécanismes, j'ai .découvert que l'augmentation de l'expression de la sous-unité NR2A au niveau synaptique modifie les propriétés de la réponse du récepteur NMDA aux protocoles de stimulations utilisés pour induire la plasticité. Cette étude représente un pas en avant important dans la compréhension de la base mécaniste sous-tendant le rôle du sommeil pour le cerveau. En montrant que le sommeil et la perte de sommeil affectent la plasticité neuronale en régulant l'expression et la fonction d'un récepteur de la neurotransmission, je propose qu'un aspect important de la fonction du sommeil puisse être finalisé au règlement de la redistribution des protéines et du tracking des récepteurs aux synapses centraux. Ces découvertes fournissent un point de départ pour mieux comprendre les liens entre le sommeil et l'apprentissage, et d'ailleurs, ils peuvent ouvrir des voies pour des traitements pharmacologiques dans le .but de préserver la fonction hippocampale pendant les périodes de restriction de sommeil.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil.

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.

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.