Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Chemodiversity and molecular plasticity: recognition processes as explored by property spaces.

In the last few years, a need to account for molecular flexibility in drug-design methodologies has emerged, even if the dynamic behavior of molecular properties is seldom made explicit. For a flexible molecule, it is indeed possible to compute different values for a given conformation-dependent property and the ensemble of such values defines a property space that can be used to describe its molecular variability; a most representative case is the lipophilicity space. In this review, a number of applications of lipophilicity space and other property spaces are presented, showing that this concept can be fruitfully exploited: to investigate the constraints exerted by media of different levels of structural organization, to examine processes of molecular recognition and binding at an atomic level, to derive informative descriptors to be included in quantitative structure--activity relationships and to analyze protein simulations extracting the relevant information. Much molecular information is neglected in the descriptors used by medicinal chemists, while the concept of property space can fill this gap by accounting for the often-disregarded dynamic behavior of both small ligands and biomacromolecules. Property space also introduces some innovative concepts such as molecular sensitivity and plasticity, which appear best suited to explore the ability of a molecule to adapt itself to the environment variously modulating its property and conformational profiles. Globally, such concepts can enhance our understanding of biological phenomena providing fruitful descriptors in drug-design and pharmaceutical sciences.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

A web-based tool to aid the identification of chemicals potentially posing a health risk through percutaneous exposure

Occupational hygiene practitioners typically assess the risk posed by occupational exposure by comparing exposure measurements to regulatory occupational exposure limits (OELs). In most jurisdictions, OELs are only available for exposure by the inhalation pathway. Skin notations are used to indicate substances for which dermal exposure may lead to health effects. However, these notations are either present or absent and provide no indication of acceptable levels of exposure. Furthermore, the methodology and framework for assigning skin notation differ widely across jurisdictions resulting in inconsistencies in the substances that carry notations. The UPERCUT tool was developed in response to these limitations. It helps occupational health stakeholders to assess the hazard associated with dermal exposure to chemicals. UPERCUT integrates dermal quantitative structure-activity relationships (QSARs) and toxicological data to provide users with a skin hazard index called the dermal hazard ratio (DHR) for the substance and scenario of interest. The DHR is the ratio between the estimated 'received' dose and the 'acceptable' dose. The 'received' dose is estimated using physico-chemical data and information on the exposure scenario provided by the user (body parts exposure and exposure duration), and the 'acceptable' dose is estimated using inhalation OELs and toxicological data. The uncertainty surrounding the DHR is estimated with Monte Carlo simulation. Additional information on the selected substances includes intrinsic skin permeation potential of the substance and the existence of skin notations. UPERCUT is the only available tool that estimates the absorbed dose and compares this to an acceptable dose. In the absence of dermal OELs it provides a systematic and simple approach for screening dermal exposure scenarios for 1686 substances.

Structure function relationships of ENaC and its role in sodium handling.

The epithelial sodium channel (ENaC) in the apical membrane of polarized epithelial cells is the rate-limiting step for Na entry into the cell; in series with the basolateral Na pump, it allows the vectorial transepithelial transport of Na ions. ENaC is expressed in different epithelia like the distal nephron or colon, and the airways epithelium. In the lung ENaC controls the composition and the amount of pulmonary fluid, whereas in the distal nephron ENaC under the control of aldosterone and vasopressin, is essential to adapt the amount of Na+ reabsorbed with the daily sodium intake. Activating mutations of ENaC cause severe disturbances of Na+ homeostasis leading to hypertension in human and in mouse models. Functional expression of ENaC in different cell systems allowed the identification of structural domains of the protein that are essential for channel function and/or modulation of channel activity. Site-directed mutations in specific domains of the channel protein lead to channel hyperactivity or channel loss of function. Knowledge about ENaC structure-function relationships opens new opportunities for development of pharmacological tools for controlling ENaC activity, such as channel activators of potential benefit in the treatment of pulmonary edema, or highly potent ENaC blockers with natriuretic effects.

Structure-function relationships of the alpha1b-adrenergic receptor.

The alpha1b-adrenergic receptor (AR) is a member of the large superfamily of seven transmembrane domain (TMD) G protein-coupled receptors (GPCR). Combining site-directed mutagenesis of the alpha1b-AR with computational simulations of receptor dynamics, we have explored the conformational changes underlying the process of receptor activation, i.e. the transition between the inactive and active states. Our findings suggest that the structural constraint stabilizing the alpha1b-AR in the inactive form is a network of H-bonding interactions amongst conserved residues forming a polar pocket and R143 of the DRY sequence at the end of TMDIII. We have recently reported that point mutations of D142, of the DRY sequence and of A293 in the distal portion of the third intracellular loop resulted in ligand-independent (constitutive) activation of the alpha1b-AR. These constitutively activating mutations could induce perturbations resulting in the shift of R143 out of the polar pocket. The main role of R143 may be to mediate receptor activation by triggering the exposure of several basic amino acids of the intracellular loops towards the G protein. Our investigation has been extended also to the biochemical events involved in the desensitization process of alpha1b-AR. Our results indicate that immediately following agonist-induced activation, the alpha1b-AR can undergo rapid agonist-induced phosphorylation and desensitization. Different members of the G protein coupled receptor kinase family can play a role in agonist-induced regulation of the alpha1b-AR. In addition, constitutively active alpha1b-AR mutants display different phosphorylation and internalization features. The future goal is to further elucidate the molecular mechanism underlying the complex equilibrium between activation and inactivation of the alpha1b-AR and its regulation by pharmacological substances. These findings can help to elucidate the mechanism of action of various agents displaying properties of agonists or inverse agonists at the adrenergic system.

The activation process of the alpha1B-adrenergic receptor: potential role of protonation and hydrophobicity of a highly conserved aspartate.

In this study, a quantitative approach was used to investigate the role of D142, which belongs to the highly conserved E/DRY sequence, in the activation process of the alpha1B-adrenergic receptor (alpha1B-AR). Experimental and computer-simulated mutagenesis were performed by substituting all possible natural amino acids at the D142 site. The resulting congeneric set of proteins together with the finding that all the receptor mutants show various levels of constitutive (agonist-independent) activity enabled us to quantitatively analyze the relationships between structural/dynamic features and the extent of constitutive activity. Our results suggest that the hydrophobic/hydrophilic character of D142, which could be regulated by protonation/deprotonation of this residue, is an important modulator of the transition between the inactive (R) and active (R*) state of the alpha1B-AR. Our study represents an example of quantitative structure-activity relationship analysis of the activation process of a G protein-coupled receptor.

PPAR alpha structure-function relationships derived from species-specific differences in responsiveness to hypolipidemic agents.

The nuclear receptor PPAR alpha is a key regulatory transcription factor in lipid homeostasis, some liver detoxification processes and the control of inflammation. Recent findings suggest that many hypolipidemic drugs and anti-inflammatory agents can potentially act by binding to PPAR alpha and inducing its activity. Here, we identify some structure-function relationships in PPAR alpha, by using the species-specific responsiveness to the two hypolipidemic agents, Wy 14,643 and 5,8,11,14-eicosatetraynoic acid (ETYA). We first show that the species-specific differences are mediated primarily via the ligand binding domain of the receptor and that these two drugs are indeed ligands of PPAR alpha. By mutagenesis analyses we identify amino acid residues in the ligand binding domains of Xenopus, mouse and human PPAR alpha, that confer preferential responsiveness to ETYA and Wy 14,643. These findings will aid in the development of new synthetic PPAR alpha ligands as effective therapeutics for lipid-related diseases and inflammatory disorders.

Evolutionary aspects of population structure for molecular and quantitative traits in the freshwater snail Radix balthica.

Detecting the action of selection in natural populations can be achieved using the QST-FST comparison that relies on the estimation of FST with neutral markers, and QST using quantitative traits potentially under selection. QST higher than FST suggests the action of directional selection and thus potential local adaptation. In this article, we apply the QST-FST comparison to four populations of the hermaphroditic freshwater snail Radix balthica located in a floodplain habitat. In contrast to most studies published so far, we did not detect evidence of directional selection for local optima for any of the traits we measured: QST calculated using three different methods was never higher than FST. A strong inbreeding depression was also detected, indicating that outcrossing is probably predominant over selfing in the studied populations. Our results suggest that in this floodplain habitat, local adaptation of R. balthica populations may be hindered by genetic drift, and possibly altered by uneven gene flow linked to flood frequency.

Peroxisome proliferator-activated receptor structures: ligand specificity, molecular switch and interactions with regulators.

Peroxisome proliferator-activated receptors (PPARs) compose a family of nuclear receptors that mediate the effects of lipidic ligands at the transcriptional level. In this review, we highlight advances in the understanding of the PPAR ligand binding domain (LBD) structure at the atomic level. The overall structure of PPARs LBD is described, and important protein ligand interactions are presented. Structure-activity relationships between isotypes structures and ligand specificity are addressed. It is shown that the numerous experimental three-dimensional structures available, together with in silico simulations, help understanding the role played by the activating function-2 (AF-2) in PPARs activation and its underlying molecular mechanism. The relation between the PPARs constitutive activity and the intrinsic stability of the active conformation is discussed. Finally, the interactions of PPARs LBD with co-activators or co-repressors, as well as with the retinoid X receptor (RXR) are described and considered in relation to PPARs activation.

Mutations of the Serine Protease CAP1/Prss8 Lead to Reduced Embryonic Viability, Skin Defects, and Decreased ENaC Activity.

CAP1/Prss8 is a membrane-bound serine protease involved in the regulation of several different effectors, such as the epithelial sodium channel ENaC, the protease-activated receptor PAR2, the tight junction proteins, and the profilaggrin polypeptide. Recently, the V170D and the G54-P57 deletion mutations within the CAP1/Prss8 gene, identified in mouse frizzy (fr) and rat hairless (fr(CR)) animals, respectively, have been proposed to be responsible for their skin phenotypes. In the present study, we analyzed those mutations, revealing a change in the protein structure, a modification of the glycosylation state, and an overall reduction in the activation of ENaC of the two mutant proteins. In vivo analyses demonstrated that both fr and fr(CR) mutant animals present analogous reduction of embryonic viability, similar histologic aberrations at the level of the skin, and a significant decrease in the activity of ENaC in the distal colon compared with their control littermates. Hairless rats additionally had dehydration defects in skin and intestine and significant reduction in the body weight. In conclusion, we provided molecular and functional evidence that CAP1/Prss8 mutations are accountable for the defects in fr and fr(CR) animals, and we furthermore demonstrate a decreased function of the CAP1/Prss8 mutant proteins. Therefore, fr and fr(CR) animals are suitable models to investigate the consequences of CAP1/Prss8 action on its target proteins in the whole organism.

Interaction of alpha-melanotropin with central dopamine systems: role of hormonal state and molecular structure

The tubero-infundibular and nigrostriatal DA neurone systems of rats respond to systemic (i.p.) injection of alpha-MSH (2-100 microgram/kg). The response of the tubero-infundibular (arcuate) DA neurones, an increase in cellular fluorescence intensity which can be interpreted as a sign of increased neuronal activity, is essentially the same in males, estrogen-progesterone-pretreated ovariectomized females and hypophysectomized males, whereas the type of response elicited by alpha-MSH in the nigral DA neurones depends upon the hormonal state of the animal. Differences between the two DA neurone groups exist also with regard to the effects of peptide fragments containing the two active sites of the alpha-MSH molecule. Results of lesion experiments in the lower brainstem (area postrema) and of blockade of muscarinic mechanisms by atropine further point to differences in the mechanisms underlying the peptide effects on the two neurone systems. The reaction of the tubero-infundibular DA system (which controls the pars intermedia of the pituitary) can be considered to reflect the activation of a feedback mechanism on MSH secretion, while the functional counterpart of the changes observed in the nigral system remains unknown at the present time.

Scaffold attachment factor B1 directly interacts with nuclear receptors in living cells and represses transcriptional activity.

Transcriptional activity relies on coregulators that modify the chromatin structure and serve as bridging factors between transcription factors and the basal transcription machinery. Using the DE domain of human peroxisome proliferator-activated receptor gamma (PPARgamma) as bait in a yeast two-hybrid screen of a human adipose tissue library, we isolated the scaffold attachment factor B1 (SAFB1/HET/HAP), which was previously shown to be a corepressor of estrogen receptor alpha. We show here that SAFB1 has a very broad tissue expression profile in human and is also expressed all along mouse embryogenesis. SAFB1 interacts in pull-down assays not only with PPARgamma but also with all nuclear receptors tested so far, albeit with different affinities. The association of SAFB1 and PPARgamma in vivo is further demonstrated by fluorescence resonance energy transfer (FRET) experiments in living cells. We finally show that SAFB1 is a rather general corepressor for nuclear receptors. Its change in expression during the early phases of adipocyte and enterocyte differentiation suggests that SAFB1 potentially influences cell proliferation and differentiation decisions.

Qualitative and quantitative analysis of thalamo-cortical axons in the somatosensory cortex of normal and barrelless mice

Résumé Les rongeurs utilisent leurs moustaches (vibrisses) pour explorer le milieu environnant. Chaque moustache est mue par un système des muscles. Les récepteurs situés à sa base transmettent les informations au système nerveux central. La transmission vers l'écorce se fait via trois neurones de relais qui se trouvent au niveau du ganglion trigéminé, du tronc cérébral et du thalamus. La représentation corticale d'une vibrisse est une concentration des axones thalamo-corticaux (ATC) autour desquelles s'organisent leurs cibles, les cellules de la couche IV. La structure peut être identifiée histologiquement en coupes tangentielles et porte le nom de « barrel » (« tonneau »). Cette correspondance vibrisse - barrel fait de ce système un model idéal pour étudier l'influence de l'activité périphérique sur l'établissement et le maintien des cartes somatotopiques. Notre laboratoire dispose d'une souche de souris qui a subi une mutation spontanée pour le gène codant l'adenylyl cyclase I (ACI). Cette enzyme membranaire catalyse la formation de l'AMPc et joue un rôle important dans le guidage axonal, la libération des neurotransmetteurs et l'intégration des signaux postsynaptiques. Nous avons démontré dans un premier temps que cette souris adulte ne développe pas de barrels. Cela est dû à un manque d'organisation des ATC et aussi des cellules de la couche IV. De plus, les résultats électrophysiologiques montrent que les informations venant des vibrisses adjacentes ne sont pas intégrées d'une manière normale. Dans ce travail de thèse, j'ai analysé la morphologie des ATC révélés individuellement avec de la biocytine. L'analyse quantitative des ATC a mis en évidence les points suivants: 1. Les axones de la souris normale (NOR) quittent le thalamus, traversent la capsule interne et la substance blanche sous-corticale et pénètrent dans le cortex somato-sensoriel primaire. A l'intérieur de l'écorce ils traversent au maximum 3 colonnes corticales adjacentes dont une contient le barrel cible. En passant à travers les couches VI et V, ces axones arborisent et convergent progressivement vers le barrel dans lequel ils forment une riche arborisation. Un petit nombre des branches « errantes », pleines de boutons synaptiques, pénètrent dans les barrels voisins. Deux axones NOR provenant de corps cellulaires très proches dans le thalamus peuvent avoir un cheminement très divergent lors de la traversée de la capsule interne et de la substance blanche sous-corticale mais, à leur entrée dans le cortex, ils sont distants d'au maximum 2 colonnes corticales de la colonne qui contient le barrel cible et ils convergent progressivement vers ce barrel. 2. Les axones de la souris mutante (BRL) ont le même trajet sous-cortical que les axones NOR, mais leur entrée dans le cortex somato-sensoriel primaire est aléatoire. A l'interface entre la substance blanche sous-corticale et le cortex, l'axone principal se divise rapidement en troncs axonaux qui traversent les couches VI et V d'une manière divergente pour arriver dans la couche IV. Cela contraste beaucoup avec la trajectoire des NOR qui convergent graduellement vers leur barrel cible. Le nombre de branches radiales que les axones BRL utilisent pour entrer dans le cortex et dans la couche IV est double par rapport aux axones NOR. Parmi ces branches, seules quelques-unes donnent des arborisations, les autres ne sont pas développées et leur morphologie est semblable à celle des branches formées par les axones de la souris normale lors du développement. Deux axones BRL issus de corps cellulaires proches dans le thalamus peuvent avoir une trajectoire très divergente jusqu'à leur entrée dans la couche IV, mais à ce niveau ils sont réorientés pour se retrouver et faire un nombre maximal de branches et boutons synaptiques dans la même région corticale. Dans un cas extrême, un des axones observés est entré dans le cortex à la limite entre l'aire somatosensorielle primaire et secondaire et a parcouru une distance de 2 mm pour retrouver son partenaire thalamique et donner avec celui-ci un nombre maximal de branches dans la même région de la couche IV. 3. Les mesures quantitatives ont montré que les arborisations corticales des axones NOR ont une longueur moyenne de 18mm et sont formées par 200 segments qui portent 1200 boutons synaptiques. Par rapport à la souris NOR, les axones BRL ont en moyenne la même longueur, le même nombre de segments et boutons synaptiques, mais donnent deux fois plus de branches radiales. La surface tangentielle occupée par les arborisations BRL dans la couche IV est 2 fois plus grande que celle des NOR. Cela signifie que les 1000 boutons synaptiques qui caractérisent les arborisations NOR et BRL dans la couche IV sont disséminés sur une surface tangentielle double chez les derniers, et donc que la densité des boutons par unité de surface corticale est en moyenne plus faible. En effet, l'augmentation de la surface corticale tangentielle des BRL est due aux surfaces de faible et moyenne densité synaptique (0 - 8 boutons / 400pn2) qui augmentent 2 fois tandis que les surfaces de haute densité synaptiques (8 - 64 boutons / 4001.tm2) sont les mêmes. Nous émettons l'hypothèse selon laquelle, durant le développement, les ATC de la souris BRL divergent et forment un nombre exubérant de branches. Grâce à cette divergence et aux branches supranuméraires, ils trouvent l'endroit de l'écorce où se trouvent leurs voisins thalamiques et arborisent abondamment dans cette région. Cependant, le déficit en AGI ne leurs permet pas par la suite, sous influence de l'activité périphérique, de retirer les branches qui se trouvent dans les endroits inappropriés de l'écorce, avec de possibles conséquences sur la discrimination tactile.

The permeability P-glycoprotein: a focus on enantioselectivity and brain distribution.

IMPORTANCE OF THE FIELD: The permeability glycoprotein (P-gp) is an important protein transporter involved in the disposition of many drugs with different chemical structures, but few studies have examined a possible stereoselectivity in its activity. P-gp can have a major impact on the distribution of drugs in selected organs, including the brain. Polymorphisms of the ABCB1 gene, which encodes for P-gp, can influence the kinetics of several drugs. AREAS COVERED IN THIS REVIEW: A search including publications from 1990 up to 2009 was performed on P-gp stereoselectivity and on the impact of ABCB1 polymorphisms on enantiomer brain distribution. WHAT THE READER WILL GAIN: Despite stereoselectivity not being expected because of the large variability of chemical structures of P-gp substrates, structure-activity relationships suggest different P-gp-binding sites for enantiomers. Enantioselectivity in the activity of P-gp has been demonstrated by in vitro studies and in animal models (preferential transport of one enantiomer or different inhibitory potencies towards P-gp activity between enantiomers). There is also in vivo evidence of an enantioselective drug transport at the human blood-brain barrier. TAKE HOME MESSAGE: The significant enantioselective activity of P-gp might be clinically relevant and must be taken into account in future studies.