In silico target fishing for rationalized ligand discovery exemplified on constituents of Ruta graveolens

The identification of targets whose interaction is likely to result in the successful treatment of a disease is of growing interest for natural product scientists. In the current study we performed an exemplary application of a virtual parallel screening approach to identify potential targets for 16 secondary metabolites isolated and identified from the aerial parts of the medicinal plant RUTA GRAVEOLENS L. Low energy conformers of the isolated constituents were simultaneously screened against a set of 2208 pharmacophore models generated in-house for the IN SILICO prediction of putative biological targets, i. e., target fishing. Based on the predicted ligand-target interactions, we focused on three biological targets, namely acetylcholinesterase (AChE), the human rhinovirus (HRV) coat protein and the cannabinoid receptor type-2 (CB (2)). For a critical evaluation of the applied parallel screening approach, virtual hits and non-hits were assayed on the respective targets. For AChE the highest scoring virtual hit, arborinine, showed the best inhibitory IN VITRO activity on AChE (IC (50) 34.7 muM). Determination of the anti-HRV-2 effect revealed 6,7,8-trimethoxycoumarin and arborinine to be the most active antiviral constituents with IC (50) values of 11.98 muM and 3.19 muM, respectively. Of these, arborinine was predicted virtually. Of all the molecules subjected to parallel screening, one virtual CB (2) ligand was obtained, i. e., rutamarin. Interestingly, in experimental studies only this compound showed a selective activity to the CB (2) receptor ( Ki of 7.4 muM) by using a radioligand displacement assay. The applied parallel screening paradigm with constituents of R. GRAVEOLENS on three different proteins has shown promise as an IN SILICO tool for rational target fishing and pharmacological profiling of extracts and single chemical entities in natural product research.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Cross-reactions vs co-sensitization evaluated by in silico motifs and in vitro IgE microarray testing

Using an in silico allergen clustering method, we have recently shown that allergen extracts are highly cross-reactive. Here we used serological data from a multi-array IgE test based on recombinant or highly purified natural allergens to evaluate whether co-reactions are true cross-reactions or co-sensitizations by allergens with the same motifs.

Besnoitia besnoiti protein disulfide isomerase (BbPDI): molecular characterization, expression and in silico modelling

Besnoitia besnoiti is an apicomplexan parasite responsible for bovine besnoitiosis, a disease with a high prevalence in tropical and subtropical regions and re-emerging in Europe. Despite the great economical losses associated with besnoitiosis, this disease has been underestimated and poorly studied, and neither an effective therapy nor an efficacious vaccine is available. Protein disulfide isomerase (PDI) is an essential enzyme for the acquisition of the correct three-dimensional structure of proteins. Current evidence suggests that in Neosporacaninum and Toxoplasmagondii, which are closely related to B. besnoiti, PDI play an important role in host cell invasion, is a relevant target for the host immune response, and represents a promising drug target and/or vaccine candidate. In this work, we present the nucleotide sequence of the B. besnoiti PDI gene. BbPDI belongs to the thioredoxin-like superfamily (cluster 00388) and is included in the PDI_a family (cluster defined cd02961) and the PDI_a_PDI_a'_c subfamily (cd02995). A 3D theoretical model was built by comparative homology using Swiss-Model server, using as a template the crystallographic deduced model of Tapasin-ERp57 (PDB code 3F8U chain C). Analysis of the phylogenetic tree for PDI within the phylum apicomplexa reinforces the close relationship among B. besnoiti, N. caninum and T. gondii. When subjected to a PDI-assay based on the polymerisation of reduced insulin, recombinant BbPDI expressed in E. coli exhibited enzymatic activity, which was inhibited by bacitracin. Antiserum directed against recombinant BbPDI reacted with PDI in Western blots and by immunofluorescence with B. besnoiti tachyzoites and bradyzoites.

Engineering Tocopherol Selectivity in alpha-TTP: A Combined In Vitro/In Silico Study

"We present a combined in vitro/in silico study to determine the molecular origin of the selectivity of a-tocopherol transfer" "protein (a-TTP) towards a-tocopherol. Molecular dynamics simulations combined to free energy perturbation calculations predict a binding free energy for a-tocopherol to a-TTP 8.26+2.13 kcal mol{1 lower than that of c-tocopherol. Our calculations show that c-tocopherol binds to a-TTP in a significantly distorted geometry as compared to that of the natural ligand. Variations in the hydration of the binding pocket and in the protein structure are found as well. We propose a mutation, A156L, which significantly modifies the selectivity properties of a-TTP towards the two tocopherols. In particular, our simulations predict that A156L binds preferentially to c-tocopherol, with striking structural similarities to the wild-type- a-tocopherol complex. The affinity properties are confirmed by differential scanning fluorimetry as well as in vitro competitive binding assays. Our data indicate that residue A156 is at a critical position for determination of the selectivity of a-TTP. The engineering of TTP mutants with modulating binding properties can have potential impact at industrial level for easier purification of single tocopherols from vitamin E mixtures coming from natural oils or synthetic processes. Moreover," "the identification of a c-tocopherol selective TTP offers the possibility to challenge the hypotheses for the evolutionary development of a mechanism for a-tocopherol selection in omnivorous animals."

Identification and in silico analyses of novel TGFBR1 and TGFBR2 mutations in Marfan syndrome-related disorders

Very recently, heterozygous mutations in the genes encoding transforming growth factor beta receptors I (TGFBR1) and II (TGFBR2) have been reported in Loeys-Dietz aortic aneurysm syndrome (LDS). In addition, dominant TGFBR2 mutations have been identified in Marfan syndrome type 2 (MFS2) and familial thoracic aortic aneurysms and dissections (TAAD). In the past, mutations of these genes were associated with atherosclerosis and several human cancers. Here, we report a total of nine novel and one known heterozygous sequence variants in the TGFBR1 and TGFBR2 genes in nine of 70 unrelated individuals with MFS-like phenotypes who previously tested negative for mutations in the gene encoding the extracellular matrix protein fibrillin-1 (FBN1). To assess the pathogenic impact of these sequence variants, in silico analyses were performed by the PolyPhen, SIFT, and Fold-X algorithms and by means of a 3D homology model of the TGFBR2 kinase domain. Our results showed that in all but one of the patients the pathogenic effect of at least one sequence variant is highly probable (c.722C > T, c.799A > C, and c.1460G > A in TGFBR1 and c.773T > G, c.1106G > T, c.1159G > A, c.1181G > A, and c.1561T > C in TGFBR2). These deleterious alleles occurred de novo or segregated with the disease in the families, indicating a causative association between the sequence variants and clinical phenotypes. Since TGFBR2 mutations found in patients with MFS-related disorders cannot be distinguished from heterozygous TGFBR2 mutations reported in tumor samples, we emphasize the importance of segregation analysis in affected families. In order to be able to find the mutation that is indeed responsible for a MFS-related phenotype, we also propose that genetic testing for sequence alterations in TGFBR1 and TGFBR2 should be complemented by mutation screening of the FBN1 gene.

Pseudodominant Inheritance of Goitrous Congenital Hypothyroidism Caused by TPO Mutations: Molecular and in Silico Studies

Context and Objective: Most cases of goitrous congenital hypothyroidism (CH) from thyroid dyshormonogenesis 1) follow a recessive mode of inheritance and 2) are due to mutations in the thyroid peroxidase gene (TPO). We report the genetic mechanism underlying the apparently dominant inheritance of goitrous CH in a nonconsanguineous family of French Canadian origin. Design, Setting, and Participants: Two brothers identified by newborn TSH screening had severe hypothyroidism and a goiter with increased (99m)Tc uptake. The mother was euthyroid, but the father and two paternal uncles had also been diagnosed with goitrous CH. After having excluded PAX8 gene mutations, we hypothesized that the underlying defect could be TPO mutations. Results: Both compound heterozygous siblings had inherited a mutant TPO allele carried by their mother (c.1496delC; p.Pro499Argfs2X), and from their father, one brother had inherited a missense mutation (c.1978C-->G; p.Gln660Glu) and the other an insertion (c.1955insT; p.Phe653Valfs15X). The thyroid gland of one uncle who is a compound heterozygote for TPO mutations (p.Phe653Valfs15X/p.Gln660Glu) was removed because of concurrent multiple endocrine neoplasia type 2A. Immunohistochemistry revealed normal TPO staining, implying that Gln660Glu TPO is expressed properly. Modeling of this mutant in silico suggests that its three-dimensional structure is conserved, whereas the electrostatic binding energy between the Gln660Glu TPO and its heme group becomes repulsive. Conclusion: We report a pedigree presenting with pseudodominant goitrous CH due to segregation of three different TPO mutations. Although goitrous CH generally follows a recessive mode of inheritance, the high frequency of TPO mutations carriers may lead to pseudodominant inheritance.

Analysis of independent microarray datasets of renal biopsies identifies a robust transcript signature of acute allograft rejection

Transcriptomics could contribute significantly to the early and specific diagnosis of rejection episodes by defining 'molecular Banff' signatures. Recently, the description of pathogenesis-based transcript sets offered a new opportunity for objective and quantitative diagnosis. Generating high-quality transcript panels is thus critical to define high-performance diagnostic classifier. In this study, a comparative analysis was performed across four different microarray datasets of heterogeneous sample collections from two published clinical datasets and two own datasets including biopsies for clinical indication, and samples from nonhuman primates. We characterized a common transcriptional profile of 70 genes, defined as acute rejection transcript set (ARTS). ARTS expression is significantly up-regulated in all AR samples as compared with stable allografts or healthy kidneys, and strongly correlates with the severity of Banff AR types. Similarly, ARTS were tested as a classifier in a large collection of 143 independent biopsies recently published by the University of Alberta. Results demonstrate that the 'in silico' approach applied in this study is able to identify a robust and reliable molecular signature for AR, supporting a specific and sensitive molecular diagnostic approach for renal transplant monitoring.

An Actor-Critic based controller for glucose regulation in type 1 diabetes

A novel adaptive approach for glucose control in individuals with type 1 diabetes under sensor-augmented pump therapy is proposed. The controller, is based on Actor-Critic (AC) learning and is inspired by the principles of reinforcement learning and optimal control theory. The main characteristics of the proposed controller are (i) simultaneous adjustment of both the insulin basal rate and the bolus dose, (ii) initialization based on clinical procedures, and (iii) real-time personalization. The effectiveness of the proposed algorithm in terms of glycemic control has been investigated in silico in adults, adolescents and children under open-loop and closed-loop approaches, using announced meals with uncertainties in the order of ±25% in the estimation of carbohydrates. The results show that glucose regulation is efficient in all three groups of patients, even with uncertainties in the level of carbohydrates in the meal. The percentages in the A+B zones of the Control Variability Grid Analysis (CVGA) were 100% for adults, and 93% for both adolescents and children. The AC based controller seems to be a promising approach for the automatic adjustment of insulin infusion in order to improve glycemic control. After optimization of the algorithm, the controller will be tested in a clinical trial.

Organotypic slice cultures from rat brain tissue: a new approach for Naegleria fowleri CNS infection in vitro

The free-living amoeba Naegleria fowleri is the aetiological agent of primary amoebic meningoencephalitis (PAM), a disease leading to death in the vast majority of cases. In patients suffering from PAM, and in corresponding animal models, the brain undergoes a massive inflammatory response, followed by haemorrhage and severe tissue necrosis. Both, in vivo and in vitro models are currently being used to study PAM infection. However, animal models may pose ethical issues, are dependent upon availability of specific infrastructural facilities, and are time-consuming and costly. Conversely, cell cultures lack the complex organ-specific morphology found in vivo, and thus, findings obtained in vitro do not necessarily reflect the situation in vivo. The present study reports infection of organotypic slice cultures from rat brain with N. fowleri and compares the findings in this culture system with in vivo infection in a rat model of PAM, that proved complementary to that of mice. We found that brain morphology, as present in vivo, is well retained in organotypic slice cultures, and that infection time-course including tissue damage parallels the observations in vivo in the rat. Therefore, organotypic slice cultures from rat brain offer a new in vitro approach to study N. fowleri infection in the context of PAM.

Preclinical testing of strategies for therapeutic targeting of human T-cell trafficking in vivo

Naive T cells are migratory cells that continuously recirculate between blood and lymphoid tissues. Antigen-specific stimulation of T cells within the lymph nodes reprograms the trafficking properties of T cells by inducing a specific set of adhesion molecules and chemokine receptors on their surface which allow these activated and effector T cells to effectively and specifically home to extralymphoid organs. The observations of organ-specific homing of T cells initiated the development of therapeutic strategies targeting adhesion receptors for organ-specific inhibition of chronic inflammation. As most adhesion receptors have additional immune functions besides mediating leukocyte trafficking, these drugs may have additional immunomodulatory effects. Therapeutic targeting of T-cell trafficking to the central nervous system is the underlying concept of a novel treatment of relapsing remitting multiple sclerosis with the humanized anti-alpha-4-integrin antibody natalizumab. In this chapter, we describe a possible preclinical in vivo approach to directly visualize the therapeutic efficacy of a given drug in inhibiting T-cell homing to a certain organ at the example of the potential of natalizumab to inhibit the trafficking of human T cells to the inflamed central nervous system in an animal model of multiple sclerosis.

Genome-wide association study in German patients with attention deficit/hyperactivity disorder

The heritability of attention deficit hyperactivity disorder (ADHD) is approximately 0.8. Despite several larger scale attempts, genome-wide association studies (GWAS) have not led to the identification of significant results. We performed a GWAS based on 495 German young patients with ADHD (according to DSM-IV criteria; Human660W-Quadv1; Illumina, San Diego, CA) and on 1,300 population-based adult controls (HumanHap550v3; Illumina). Some genes neighboring the single nucleotide polymorphisms (SNPs) with the lowest P-values (best P-value: 8.38 × 10(-7)) have potential relevance for ADHD (e.g., glutamate receptor, metabotropic 5 gene, GRM5). After quality control, the 30 independent SNPs with the lowest P-values (P-values ≤ 7.57 × 10(-5) ) were chosen for confirmation. Genotyping of these SNPs in up to 320 independent German families comprising at least one child with ADHD revealed directionally consistent effect-size point estimates for 19 (10 not consistent) of the SNPs. In silico analyses of the 30 SNPs in the largest meta-analysis so far (2,064 trios, 896 cases, and 2,455 controls) revealed directionally consistent effect-size point estimates for 16 SNPs (11 not consistent). None of the combined analyses revealed a genome-wide significant result. SNPs in previously described autosomal candidate genes did not show significantly lower P-values compared to SNPs within random sets of genes of the same size. We did not find genome-wide significant results in a GWAS of German children with ADHD compared to controls. The second best SNP is located in an intron of GRM5, a gene located within a recently described region with an infrequent copy number variation in patients with ADHD.

Expression of ADAMTS1 in endothelial cells is induced by shear stress and suppressed in sprouting capillaries

ADAMTS1 inhibits capillary sprouting, and since capillary sprouts do not experience the shear stress caused by blood flow, this study undertook to clarify the relationship between shear stress and ADAMTS1. It was found that endothelial cells exposed to shear stress displayed a strong upregulation of ADAMTS1, dependent upon both the magnitude and duration of their exposure. Investigation of the underlying pathways demonstrated involvement of phospholipase C, phosphoinositide 3-kinase, and nitric oxide. Forkhead box protein O1 was identified as a likely inhibitor of the system, as its knockdown was followed by a slight increase in ADAMTS1 expression. In silico prediction displayed a transcriptional binding site for Forkhead box protein O1 in the promotor region of the ADAMTS1 gene, as well as sites for nuclear factor 1, SP1, and AP-1. The anti-angiogenic effects of ADAMTS1 were attributed to its cleavage of thrombospondin 1 into a 70-kDa fragment, and a significant enhancement of this fragment was indeed demonstrated by immunoblotting shear stress-treated cells. Accordingly, scratch wound closure displayed a slowdown in conditioned medium from shear stress-treated endothelial cells, an effect that could be completely blocked by a knockdown of thrombospondin 1 and partially blocked by a knockdown of ADAMTS1. Non-perfused capillary sprouts in rat mesenteries stained negative for ADAMTS1, while vessels in the microcirculation that had already experienced blood flow yielded the opposite results. The shear stress-dependent expression of ADAMTS1 in vitro was therefore also demonstrated in vivo and thereby confirmed as a mechanism connecting blood flow with the regulation of angiogenesis.

ENU-induced phenovariance in mice: inferences from 587 mutations

Background We present a compendium of N-ethyl-N-nitrosourea (ENU)-induced mouse mutations, identified in our laboratory over a period of 10 years either on the basis of phenotype or whole genome and/or whole exome sequencing, and archived in the Mutagenetix database. Our purpose is threefold: 1) to formally describe many point mutations, including those that were not previously disclosed in peer-reviewed publications; 2) to assess the characteristics of these mutations; and 3) to estimate the likelihood that a missense mutation induced by ENU will create a detectable phenotype. Findings In the context of an ENU mutagenesis program for C57BL/6J mice, a total of 185 phenotypes were tracked to mutations in 129 genes. In addition, 402 incidental mutations were identified and predicted to affect 390 genes. As previously reported, ENU shows strand asymmetry in its induction of mutations, particularly favoring T to A rather than A to T in the sense strand of coding regions and splice junctions. Some amino acid substitutions are far more likely to be damaging than others, and some are far more likely to be observed. Indeed, from among a total of 494 non-synonymous coding mutations, ENU was observed to create only 114 of the 182 possible amino acid substitutions that single base changes can achieve. Based on differences in overt null allele frequencies observed in phenotypic vs. non-phenotypic mutation sets, we infer that ENU-induced missense mutations create detectable phenotype only about 1 in 4.7 times. While the remaining mutations may not be functionally neutral, they are, on average, beneath the limits of detection of the phenotypic assays we applied. Conclusions Collectively, these mutations add to our understanding of the chemical specificity of ENU, the types of amino acid substitutions it creates, and its efficiency in causing phenovariance. Our data support the validity of computational algorithms for the prediction of damage caused by amino acid substitutions, and may lead to refined predictions as to whether specific amino acid changes are responsible for observed phenotypes. These data form the basis for closer in silico estimations of the number of genes mutated to a state of phenovariance by ENU within a population of G3 mice.

Common obesity risk alleles in childhood attention-deficit/hyperactivity disorder

Children with attention-deficit/hyperactivity disorder (ADHD) have a higher rate of obesity than children without ADHD. Obesity risk alleles may overlap with those relevant for ADHD. We examined whether risk alleles for an increased body mass index (BMI) are associated with ADHD and related quantitative traits (inattention and hyperactivity/impulsivity). We screened 32 obesity risk alleles of single nucleotide polymorphisms (SNPs) in a genome-wide association study (GWAS) for ADHD based on 495 patients and 1,300 population-based controls and performed in silico analyses of the SNPs in an ADHD meta-analysis comprising 2,064 trios, 896 independent cases, and 2,455 controls. In the German sample rs206936 in the NUDT3 gene (nudix; nucleoside diphosphate linked moiety X-type motif 3) was associated with ADHD risk (OR: 1.39; P = 3.4 × 10(-4) ; Pcorr  = 0.01). In the meta-analysis data we found rs6497416 in the intronic region of the GPRC5B gene (G protein-coupled receptor, family C, group 5, member B; P = 7.2 × 10(-4) ; Pcorr  = 0.02) as a risk allele for ADHD. GPRC5B belongs to the metabotropic glutamate receptor family, which has been implicated in the etiology of ADHD. In the German sample rs206936 (NUDT3) and rs10938397 in the glucosamine-6-phosphate deaminase 2 gene (GNPDA2) were associated with inattention, whereas markers in the mitogen-activated protein kinase 5 gene (MAP2K5) and in the cell adhesion molecule 2 gene (CADM2) were associated with hyperactivity. In the meta-analysis data, MAP2K5 was associated with inattention, GPRC5B with hyperactivity/impulsivity and inattention and CADM2 with hyperactivity/impulsivity. Our results justify further research on the elucidation of the common genetic background of ADHD and obesity.