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.

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.

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.

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.

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.

Polymorphism and chromosomal location of the MC4R (melanocortin-4 receptor) gene in the dog and red fox

The melanocortin-4 receptor (MC4R) is expressed in the hypothalamus and regulates energy intake and body weight. In silico screening of the canine chromosome 1 sequence and a comparison with the porcine MC4R sequence by BLAST were performed. The nucleotide sequence of the whole coding region and 3'- and 5'-flanking regions of the dog (1214 bp) and red fox (1177 bp) MC4R gene was established and high conservation of the nucleotide sequences was revealed (99%). Five sets of PCR primers were designed and a search for polymorphism was performed by the SSCP technique in a group of 31 dogs representing nineteen breeds and 35 farm red foxes. Sequencing of DNA fragments, representing the identified SSCP patterns, revealed three single nucleotide polymorphisms (including a missense one) in dogs and four silent SNPs in red foxes. An average SNP frequency was approx. 1/400 bp in the dog and 1/300 bp in the red fox. We mapped the MC4R gene by FISH to the canine chromosome 1 (CFA1q1.1) and to the red fox chromosome 5 (VVU5p1.2).

PCR and real-time PCR primers developed for detection and identification of Bifidobacterium thermophilum in faeces

BACKGROUND: Culture-independent methods based on the 16S ribosomal RNA molecule are nowadays widely used for assessment of the composition of the intestinal microbiota, in relation to host health or probiotic efficacy. Because Bifidobacterium thermophilum was only recently isolated from human faeces until now, no specific real-time PCR (qPCR) assay has been developed for detection of this species as component of the bifidobacterial community of the human intestinal flora. RESULTS: Design of specific primers and probe was achieved based on comparison of 108 published bifidobacterial 16S rDNA sequences with the recently published sequence of the human faecal isolate B. thermophilum RBL67. Specificity of the primer was tested in silico by similarity search against the sequence database and confirmed experimentally by PCR amplification on 17 Bifidobacterium strains, representing 12 different species, and two Lactobacillus strains. The qPCR assay developed was linear for B. thermophilum RBL67 DNA quantities ranging from 0.02 ng/microl to 200 ng/microl and showed a detection limit of 10(5) cells per gram faeces. The application of this new qPCR assay allowed to detect the presence of B. thermophilum in one sample from a 6-month old breast-fed baby among 17 human faecal samples tested. Additionally, the specific qPCR primers in combination with selective plating experiments led to the isolation of F9K9, a faecal isolate from a 4-month old breast-fed baby. The 16S rDNA sequence of this isolate is 99.93% similar to that of B. thermophilum RBL67 and confirmed the applicability of the new qPCR assay in faecal samples. CONCLUSION: A new B. thermophilum-specific qPCR assay was developed based on species-specific target nucleotides in the 16S rDNA. It can be used to further characterize the composition of the bifidobacterial community in the human gastrointestinal tract. Until recently, B. thermophilum was considered as a species of animal origin, but here we confirm with the application of this new PCR assay the presence of B. thermophilum strains in the human gut.

Novel autoantigens immunogenic in COPD patients

Background Chronic obstructive pulmonary disease (COPD) is a respiratory inflammatory condition with autoimmune features including IgG autoantibodies. In this study we analyze the complexity of the autoantibody response and reveal the nature of the antigens that are recognized by autoantibodies in COPD patients. Methods An array of 1827 gridded immunogenic peptide clones was established and screened with 17 sera of COPD patients and 60 healthy controls. Protein arrays were evaluated both by visual inspection and a recently developed computer aided image analysis technique. By this computer aided image analysis technique we computed the intensity values for each peptide clone and each serum and calculated the area under the receiver operator characteristics curve (AUC) for each clone and the separation COPD sera versus control sera. Results By visual evaluation we detected 381 peptide clones that reacted with autoantibodies of COPD patients including 17 clones that reacted with more than 60% of the COPD sera and seven clones that reacted with more than 90% of the COPD sera. The comparison of COPD sera and controls by the automated image analysis system identified 212 peptide clones with informative AUC values. By in silico sequence analysis we found an enrichment of sequence motives previously associated with immunogenicity. Conclusion The identification of a rather complex humoral immune response in COPD patients supports the idea of COPD as a disease with strong autoimmune features. The identification of novel immunogenic antigens is a first step towards a better understanding of the autoimmune component of COPD.