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.

Proteomic analysis of Plasmodium in the mosquito: progress and pitfalls

Here we discuss proteomic analyses of whole cell preparations of the mosquito stages of malaria parasite development (i.e. gametocytes, microgamete, ookinete, oocyst and sporozoite) of Plasmodium berghei. We also include critiques of the proteomes of two cell fractions from the purified ookinete, namely the micronemes and cell surface. Whereas we summarise key biological interpretations of the data, we also try to identify key methodological constraints we have met, only some of which we were able to resolve. Recognising the need to translate the potential of current genome sequencing into functional understanding, we report our efforts to develop more powerful combinations of methods for the in silico prediction of protein function and location. We have applied this analysis to the proteome of the male gamete, a cell whose very simple structural organisation facilitated interpretation of data. Some of the in silico predictions made have now been supported by ongoing protein tagging and genetic knockout studies. We hope this discussion may assist future studies.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

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.

An insulin infusion advisory system based on autotuning nonlinear model-predictive control

This paper aims at the development and evaluation of a personalized insulin infusion advisory system (IIAS), able to provide real-time estimations of the appropriate insulin infusion rate for type 1 diabetes mellitus (T1DM) patients using continuous glucose monitors and insulin pumps. The system is based on a nonlinear model-predictive controller (NMPC) that uses a personalized glucose-insulin metabolism model, consisting of two compartmental models and a recurrent neural network. The model takes as input patient's information regarding meal intake, glucose measurements, and insulin infusion rates, and provides glucose predictions. The predictions are fed to the NMPC, in order for the latter to estimate the optimum insulin infusion rates. An algorithm based on fuzzy logic has been developed for the on-line adaptation of the NMPC control parameters. The IIAS has been in silico evaluated using an appropriate simulation environment (UVa T1DM simulator). The IIAS was able to handle various meal profiles, fasting conditions, interpatient variability, intraday variation in physiological parameters, and errors in meal amount estimations.

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.

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.