T cell receptor-ligand interactions: a conformational preequilibrium or an induced fit.

Kinetic parameters of T cell receptor (TCR) interactions with its ligand have been proposed to control T cell activation. Analysis of kinetic data obtained has so far produced conflicting insights; here, we offer a consideration of this problem. As a model system, association and dissociation of a soluble TCR (sT1) and its specific ligand, an azidobenzoic acid derivative of the peptide SYIPSAEK-(ABA)I (residues 252-260 from Plasmodium berghei circumsporozoite protein), bound to class I MHC H-2K(d)-encoded molecule (MHCp) were studied by surface plasmon resonance. The association time courses exhibited biphasic patterns. The fast and dominant phase was assigned to ligand association with the major fraction of TCR molecules, whereas the slow component was attributed to the presence of traces of TCR dimers. The association rate constant derived for the fast phase, assuming a reversible, single-step reaction mechanism, was relatively slow and markedly temperature-dependent, decreasing from 7.0 x 10(3) at 25 degrees C to 1.8 x 10(2) M(-1).s(-1) at 4 degrees C. Hence, it is suggested that these observed slow rate constants are the result of unresolved elementary steps of the process. Indeed, our analysis of the kinetic data shows that the time courses of TCR-MHCp interaction fit well to two different, yet closely related mechanisms, where an induced fit or a preequilibrium of two unbound TCR conformers are operational. These mechanisms may provide a rationale for the reported conformational flexibility of the TCR and its unusual ligand recognition properties, which combine high specificity with considerable crossreactivity.

Stable masking by H-2Dd cis ligand limits Ly49A relocalization to the site of NK cell/target cell contact.

Ly49A is an inhibitory receptor, which counteracts natural killer (NK) cell activation on the engagement with H-2D(d) (D(d)) MHC class I molecules (MHC-I) on target cells. In addition to binding D(d) on apposed membranes, Ly49A interacts with D(d) ligand expressed in the plane of the NK cells' membrane. Indeed, multivalent, soluble MHC-I ligand binds inefficiently to Ly49A unless the NK cells' D(d) complexes are destroyed. However, it is not known whether masked Ly49A remains constitutively associated with cis D(d) also during target cell interaction. Alternatively, it is possible that Ly49A has to be unmasked to significantly interact with its ligand on target cells. These two scenarios suggest distinct roles of Ly49A/D(d) cis interaction for NK cell function. Here, we show that Ly49A contributes to target cell adhesion and efficiently accumulates at synapses with D(d)-expressing target cells when NK cells themselves lack D(d). When NK cells express D(d), Ly49A no longer contributes to adhesion, and ligand-driven recruitment to the cellular contact site is strongly reduced. The destruction of D(d) complexes on NK cells, which unmasks Ly49A, is necessary and sufficient to restore Ly49A adhesive function and recruitment to the synapse. Thus, cis D(d) continuously sequesters a considerable fraction of Ly49A receptors, preventing efficient Ly49A recruitment to the synapse with D(d)+ target cells. The reduced number of Ly49A receptors that can functionally interact with D(d) on target cells explains the modest inhibitory capacity of Ly49A in D(d) NK cells. This property renders Ly49A NK cells more sensitive to react to diseased host cells.

The Oncogene Metadherin Modulates the Apoptotic Pathway Based on the Tumor Necrosis Factor Superfamily Member TRAIL (Tumor Necrosis Factor-related Apoptosis-inducing Ligand) in Breast Cancer.

Metadherin (MTDH), the newly discovered gene, is overexpressed in more than 40% of breast cancers. Recent studies have revealed that MTDH favors an oncogenic course and chemoresistance. With a number of breast cancer cell lines and breast tumor samples, we found that the relative expression of MTDH correlated with tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) sensitivity in breast cancer. In this study, we found that knockdown of endogenous MTDH cells sensitized the MDA-MB-231 cells to TRAIL-induced apoptosis both in vitro and in vivo. Conversely, stable overexpression of MTDH in MCF-7 cells enhanced cell survival with TRAIL treatment. Mechanically, MTDH down-regulated caspase-8, decreased caspase-8 recruitment into the TRAIL death-inducing signaling complex, decreased caspase-3 and poly(ADP-ribose) polymerase-2 processing, increased Bcl-2 expression, and stimulated TRAIL-induced Akt phosphorylation, without altering death receptor status. In MDA-MB-231 breast cancer cells, sensitization to TRAIL upon MTDH down-regulation was inhibited by the caspase inhibitor Z-VAD-fmk (benzyloxycarbonyl-VAD-fluoromethyl ketone), suggesting that MTDH depletion stimulates activation of caspases. In MCF-7 breast cancer cells, resistance to TRAIL upon MTDH overexpression was abrogated by depletion of Bcl-2, suggesting that MTDH-induced Bcl-2 expression contributes to TRAIL resistance. We further confirmed that MTDH may control Bcl-2 expression partly by suppressing miR-16. Collectively, our results point to a protective function of MTDH against TRAIL-induced death, whereby it inhibits the intrinsic apoptosis pathway through miR-16-mediated Bcl-2 up-regulation and the extrinsic apoptosis pathway through caspase-8 down-regulation.

The role of Receptor Activator of NFKB Ligand (RANKL) in mammary gland development

Abstract : The female reproductive hormones estrogen, progesterone and prolactin control postnatal breast development and are important to breast carcinogenesis. The mechanisms by which they elicit proliferation and morphogenesis remain poorly understood. Using the mouse as a model to study the molecular mechanisms through which hormones elicit morphogenetic changes in the mammary gland in vivo, we found the Receptor Activator of NFκB Ligand, a Tumor Necrosis Factor family member, to be strongly induced by progesterone. Recent publications suggested that hormone dependant RANKURANK signals are involved in the terminal differentiation of mammary gland alveolar buds into lobulo-alveolar structures competent for lactation. I show that in the absence of epithelial RANKL a distinct earlier stage of mammary gland development, side branch formation, is blocked. RANKL acts as a major mediator downstream of progesterone; it is required for progesterone-induced paracrine proliferation and completely rescues the mutant phenotype when ectopically expressed in progesterone receptor (PR) KO mammary epithelia. RANKL is not required for cell autonomous division of estrogen receptor alpha (ERa) /PR positive cells. Cyclin D1, previously implicated as a mediator of RANKL, is not affected by ablation of RANKL and is not required for RANKL-induced paracrine proliferation but for the cell autonomous proliferation. Gene expression arrays to find specific RANKL downstream targets have identified Id4, ElfS and one secreted metalloprotease (Adamtsl8) as potential candidates validated by Q-RT-PCR. Interestingly, Id4 and Adamtsl8 are expressed by the myoepithelial cells. Their expression additionally coincides with RANKL mRNA expression at mid pregnancy, possibly implying a functional contribution of both genes to RANKL mediated sidebranch formation. ElfS in contrast, is found to be strongly expressed by the end of pregnancy supporting recent findings of a prolactin mediated regulation. As for RANKL, this gene was in particular induced in luminal cells. Taken together, I report that progesterone is the major proliferative stimulus in the adult mammary gland eliciting proliferation of ERaJPR positive cells by a cell autonomous, cyclin D1-dependent and a paracrine RANKL-dependent mechanism. My work moreover suggests, that RANKL acts as a major orchestrator affecting different downstream mediators, through which progesterone exerts its effects concomitantly on different cellular compartments. Résumé : Les hormones sexuelles telles que l'oestrogène, la progestérone et la prolactine contrôlent le développement postnatal du sein et sont impliquées dans la cazcinogenèse. Les mécanismes par lesquels elles induisent la prolifération et la morphogénèse demeurent incompris. En utilisant la souris comme modèle, J'ai trouvé que le ligand activateur du récepteur de NFκB, une protéine de la famille du facteur de nécrose des tumeurs, peut être fortement induit par la progestérone. Les publications récentes ont suggéré que cette protéine est nécessaire à la fin de la grossesse, quand les cellules sécrétrices du lait apparaissent. Par des techniques de transplantation d'épithélium, je montre contrairement aux études précédentes, qu'en l'absence de RANKL dans l'épithélium une partie distincte du développement mammaire, la formation de branches latérales, est bloquée. La progestérone agit de manière pazacrine par l'intermédiaire de 12ANKL pour induire la prolifération tandis que la mort cellulaire n'est pas affectée. De plus, l'injection d'une protéine recombinante RANKL dans une souris mutante pour le récepteur à la progestérone induit la prolifération des cellules épithéliales en l'absence de grossesse ; la surexpression de RANKL dans ces mêmes mutants mène à une réversion complète du phénotype. Mes expériences démontrent que la progestérone induit deux types distincts de prolifération. Un premier type direct dans laquelle les cellules positives au récepteur à la progestérone prolifèrent. Cette division cellulaire est alors indépendante de RANKL mais dépendante de la cycline D1. Le second type de prolifération est induit par un mécanisme pazacrine et dépend de RANKL mais pas de la cycline D1. Ici, les cellules négatives au récepteur à la progestérone prolifèrent. Pour détecter des gènes cibles de la voie de signalisation du RANKL, un profil d'expression des gènes a été généré. Les facteurs de transcription Id4, EIf5 et une métalloprotéase sécrétée (Adamtsl8) ont été identifiés en tant que cibles potentielles. D'autres analyses de validation démontrent qu'Id4, Adamtsl8, RANKL mais pas E1f5 sont fortement exprimés au cours de la grossesse, coïncidant avec la formation de branchements latéraux induit par progestérone. EIf5 s'est avéré être exprimé vers la fin de la grossesse appuyant des résultats récents proposant une régulation par la prolactine. Le système canalaire mammaire se compose de couches cellulaires: une couche interne de cellules luminales et une externe de cellules myoépithéliale. Les expériences génétiques d'expression ont révélé que RANKL. et E1f5 sont exprimés dans la partie luminale tandis qu'Id4 et Adamtsl8 sont dans les cellules myoépithéliales. En conclusion, je prouve que la progestérone est le stimulus principal induisant la prolifération dans la glande mammaire d'adulte. Deux mécanismes de prolifération sont impliqués: l'un direct dépendant de la cycline Dl et l'autre paracrine dépendant de RANKI.. Mon travail suggère par ailleurs que RANKL agit en tant que médiateur important, par lequel la progestérone exerce ses effets sur différents compartiments cellulaires tels que la coordination de la prolifération des cellules épithéliales avec la réorganisation de la matrice extracellulaire et de la membrane basale exigées pour la morphogénèse du système canalaire latéral.

Characterization of Fas (Apo-1, CD95)-Fas ligand interaction.

The death-inducing receptor Fas is activated when cross-linked by the type II membrane protein Fas ligand (FasL). When human soluble FasL (sFasL, containing the extracellular portion) was expressed in human embryo kidney 293 cells, the three N-linked glycans of each FasL monomer were found to be essential for efficient secretion. Based on the structure of the closely related lymphotoxin alpha-tumor necrosis factor receptor I complex, a molecular model of the FasL homotrimer bound to three Fas molecules was generated using knowledge-based protein modeling methods. Point mutations of amino acid residues predicted to affect the receptor-ligand interaction were introduced at three sites. The F275L mutant, mimicking the loss of function murine gld mutation, exhibited a high propensity for aggregation and was unable to bind to Fas. Mutants P206R, P206D, and P206F displayed reduced cytotoxicity toward Fas-positive cells with a concomitant decrease in the binding affinity for the recombinant Fas-immunoglobulin Fc fusion proteins. Although the cytotoxic activity of mutant Y218D was unaltered, mutant Y218R was inactive, correlating with the prediction that Tyr-218 of FasL interacts with a cluster of three basic amino acid side chains of Fas. Interestingly, mutant Y218F could induce apoptosis in murine, but not human cells.

Specific targeting of pro-death NMDA receptor signals with differing reliance on the NR2B PDZ ligand.

NMDA receptors (NMDARs) mediate ischemic brain damage, for which interactions between the C termini of NR2 subunits and PDZ domain proteins within the NMDAR signaling complex (NSC) are emerging therapeutic targets. However, expression of NMDARs in a non-neuronal context, lacking many NSC components, can still induce cell death. Moreover, it is unclear whether targeting the NSC will impair NMDAR-dependent prosurvival and plasticity signaling. We show that the NMDAR can promote death signaling independently of the NR2 PDZ ligand, when expressed in non-neuronal cells lacking PSD-95 and neuronal nitric oxide synthase (nNOS), key PDZ proteins that mediate neuronal NMDAR excitotoxicity. However, in a non-neuronal context, the NMDAR promotes cell death solely via c-Jun N-terminal protein kinase (JNK), whereas NMDAR-dependent cortical neuronal death is promoted by both JNK and p38. NMDAR-dependent pro-death signaling via p38 relies on neuronal context, although death signaling by JNK, triggered by mitochondrial reactive oxygen species production, does not. NMDAR-dependent p38 activation in neurons is triggered by submembranous Ca(2+), and is disrupted by NOS inhibitors and also a peptide mimicking the NR2B PDZ ligand (TAT-NR2B9c). TAT-NR2B9c reduced excitotoxic neuronal death and p38-mediated ischemic damage, without impairing an NMDAR-dependent plasticity model or prosurvival signaling to CREB or Akt. TAT-NR2B9c did not inhibit JNK activation, and synergized with JNK inhibitors to ameliorate severe excitotoxic neuronal loss in vitro and ischemic cortical damage in vivo. Thus, NMDAR-activated signals comprise pro-death pathways with differing requirements for PDZ protein interactions. These signals are amenable to selective inhibition, while sparing synaptic plasticity and prosurvival signaling.

Competition between SOCS36E and Drk modulates Sevenless receptor tyrosine kinase activity

Modulation of signalling pathways can trigger different cellular responses, including differences in cell fate. This modulation can be achieved by controlling the pathway activity with great precision to ensure robustness and reproducibility of the specification of cell fate. The development of the photoreceptor R7 in the Drosophila melanogasterretina has become a model in which to investigate the control of cell signalling. During R7 specification, a burst of Ras small GTPase (Ras) and mitogen-activated protein kinase (MAPK) controlled by Sevenless receptor tyrosine kinase (Sev) is required. Several cells in each ommatidium express sev. However, the spatiotemporal expression of the boss ligand and the action of negative regulators of the Sev pathway will restrict the R7 fate to a single cell. The Drosophila suppressor of cytokine signalling 36E (SOCS36E) protein contains an SH2 domain and acts as a Sev signalling attenuator. By contrast, downstream of receptor kinase (Drk), the fly homolog of the mammalian Grb2 adaptor protein, which also contains an SH2 domain, acts as a positive activator of the pathway. Here, we apply the Förster resonance energy transfer (FRET) assay to transfected Drosophila S2 cells and demonstrate that Sev binds directly to either the suppressor protein SOCS36E or the adaptor protein Drk. We propose a mechanistic model in which the competition between these two proteins for binding to the same docking site results in either attenuation of the Sev transduction in cells that should not develop R7 photoreceptors or amplification of the Ras-MAPK signal only in the R7 precursor.

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

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

Identification of Rare High-Avidity, Tumor-Reactive CD8+ T Cells by Monomeric TCR-Ligand Off-Rates Measurements on Living Cells.

The avidity of the T-cell receptor (TCR) for antigenic peptides presented by the peptide-MHC (pMHC) on cells is a key parameter for cell-mediated immunity. Yet a fundamental feature of most tumor antigen-specific CD8(+) T cells is that this avidity is low. In this study, we addressed the need to identify and select tumor-specific CD8(+) T cells of highest avidity, which are of the greatest interest for adoptive cell therapy in patients with cancer. To identify these rare cells, we developed a peptide-MHC multimer technology, which uses reversible Ni(2+)-nitrilotriacetic acid histidine tags (NTAmers). NTAmers are highly stable but upon imidazole addition, they decay rapidly to pMHC monomers, allowing flow-cytometric-based measurements of monomeric TCR-pMHC dissociation rates of living CD8(+) T cells on a wide avidity spectrum. We documented strong correlations between NTAmer kinetic results and those obtained by surface plasmon resonance. Using NTAmers that were deficient for CD8 binding to pMHC, we found that CD8 itself stabilized the TCR-pMHC complex, prolonging the dissociation half-life several fold. Notably, our NTAmer technology accurately predicted the function of large panels of tumor-specific T cells that were isolated prospectively from patients with cancer. Overall, our results demonstrated that NTAmers are effective tools to isolate rare high-avidity cytotoxic T cells from patients for use in adoptive therapies for cancer treatment.

Ligand binding to heparan sulfate proteoglycans induces their aggregation and distribution along actin cytoskeleton

Cell surface heparan sulfate proteoglycans (HSPGs) participate in molecular events that regulate cell adhesion, migration, and proliferation. The present study demonstrates that soluble heparin-binding proteins or cross-linking antibodies induce the aggregation of cell surface HSPGs and their distribution along underlying actin filaments. Immunofluorescence and confocal microscopy and immunogold and electron microscopy indicate that, in the absence of ligands, HSPGs are irregularly distributed on the fibroblast cell surface, without any apparent codistribution with the actin cytoskeleton. In the presence of ligand (lipoprotein lipase) or antibodies against heparan sulfate, HSPGs aggregate and colocalize with the actin cytoskeleton. Triton X-100 extraction and immunoelectron microscopy have demonstrated that in this condition HSPGs were clustered and associated with the actin filaments. Crosslinking experiments that use biotinylated lipoprotein lipase have revealed three major proteoglycans as binding sites at the fibroblast cell surface. These cross-linked proteoglycans appeared in the Triton X-100 insoluble fraction. Platinum/carbon replicas of the fibroblast surface incubated either with lipoprotein lipase or antiheparan sulfate showed large aggregates of HSPGs regularly distributed along cytoplasmic fibers. Quantification of the spacing between HSPGs by confocal microscopy confirmed that the nonrandom distribution of HSPG aggregates along the actin cytoskeleton was induced by ligand binding. When cells were incubated either with lipoprotein lipase or antibodies against heparan sulfate, the distance between immunofluorescence spots was uniform. In contrast, the spacing between HSPGs on fixed cells not incubated with ligand was more variable. This highly organized spatial relationship between actin and proteoglycans suggests that cortical actin filaments could organize the molecular machinery involved in signal transduction and molecular movements on the cell surface that are triggered by heparin-binding proteins.

Programmed death ligand 1 expression and tumor-infiltrating lymphocytes in glioblastoma.

BACKGROUND: Immune checkpoint inhibitors targeting programmed cell death 1 (PD1) or its ligand (PD-L1) showed activity in several cancer types. METHODS: We performed immunohistochemistry for CD3, CD8, CD20, HLA-DR, phosphatase and tensin homolog (PTEN), PD-1, and PD-L1 and pyrosequencing for assessment of the O6-methylguanine-methyltransferase (MGMT) promoter methylation status in 135 glioblastoma specimens (117 initial resection, 18 first local recurrence). PD-L1 gene expression was analyzed in 446 cases from The Cancer Genome Atlas. RESULTS: Diffuse/fibrillary PD-L1 expression of variable extent, with or without interspersed epithelioid tumor cells with membranous PD-L1 expression, was observed in 103 of 117 (88.0%) newly diagnosed and 13 of 18 (72.2%) recurrent glioblastoma specimens. Sparse-to-moderate density of tumor-infiltrating lymphocytes (TILs) was found in 85 of 117 (72.6%) specimens (CD3+ 78/117, 66.7%; CD8+ 52/117, 44.4%; CD20+ 27/117, 23.1%; PD1+ 34/117, 29.1%). PD1+ TIL density correlated positively with CD3+ (P < .001), CD8+ (P < .001), CD20+ TIL density (P < .001), and PTEN expression (P = .035). Enrichment of specimens with low PD-L1 gene expression levels was observed in the proneural and G-CIMP glioblastoma subtypes and in specimens with high PD-L1 gene expression in the mesenchymal subtype (P = 5.966e-10). No significant differences in PD-L1 expression or TIL density between initial and recurrent glioblastoma specimens or correlation of PD-L1 expression or TIL density with patient age or outcome were evident. CONCLUSION: TILs and PD-L1 expression are detectable in the majority of glioblastoma samples but are not related to outcome. Because the target is present, a clinical study with specific immune checkpoint inhibitors seems to be warranted in glioblastoma.

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.

Protein pocket and ligand shape comparison and its application in virtual screening.

Understanding molecular recognition is one major requirement for drug discovery and design. Physicochemical and shape complementarity between two binding partners is the driving force during complex formation. In this study, the impact of shape within this process is analyzed. Protein binding pockets and co-crystallized ligands are represented by normalized principal moments of inertia ratios (NPRs). The corresponding descriptor space is triangular, with its corners occupied by spherical, discoid, and elongated shapes. An analysis of a selected set of sc-PDB complexes suggests that pockets and bound ligands avoid spherical shapes, which are, however, prevalent in small unoccupied pockets. Furthermore, a direct shape comparison confirms previous studies that on average only one third of a pocket is filled by its bound ligand, supplemented by a 50 % subpocket coverage. In this study, we found that shape complementary is expressed by low pairwise shape distances in NPR space, short distances between the centers-of-mass, and small deviations in the angle between the first principal ellipsoid axes. Furthermore, it is assessed how different binding pocket parameters are related to bioactivity and binding efficiency of the co-crystallized ligand. In addition, the performance of different shape and size parameters of pockets and ligands is evaluated in a virtual screening scenario performed on four representative targets.

Comparing assays for the measurement of AhR ligand activity

In this thesis (TFG) the results of the comparison of three assays for the measurement of AhR ligand activity are exposed. This study was part of a collaborative project aiming at the characterization of the AhR signaling activities of known naturally occurring compounds to explore the potential of using non-toxic compounds to treat inflammatory diseases via oral administration. The first goal of this project was to find an assay able to measure AhR-activity, so the comparison of different assays has been done in order to find the most convenient one according to the efficiency, sensitivity and precision. Moreover, other elements with operational nature such as price, toxicity of components or ease of use has been considered. From the use of compounds known from the literature to be AhR ligands, three assays have been tested: (1) P450-GloTM CYP1A2 Induction/Inhibition assay, (2) quantitative Polymerase Chain Reaction (qPCR) and (3) DR. CALUX® Bioassay. Moreover, a different experiment using the last assay was performed for the study in vivo of the transport of the compounds tested. The results of the TFG suggested the DR. CALUX® Bioassay as the most promising assay to be used for the screening of samples as AhR-ligands because it is quicker, easier to handle and less expensive than qPCR and more reproducible than the CYP1A2 Induction/Inhibition assay. Moreover, the use of this assay allowed having a first idea of which compounds are uptaken by the epithelial barrier and in with direction the transport happens.

rDock: A Fast, Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrodinger's Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/