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.

Screening, Identification and Preliminary Investigation of Target Transporters in Pregnancy Pathologies

Screening, Identification and Preliminary Investigation of Target Transporters in Pregnancy Pathologies. INTRODUCTION: Pre-eclampsia (PE), intrauterine growth restriction (IUGR) and gestational diabetes mellitus (GDM) are major sources of clinical morbidity and mortality in pregnant women worldwide. The mechanisms underlying these gestational diseases are complex and not yet fully understood, but one factor contributing to their development is impaired maternal-fetal nutrient transport. Therefore, we aimed to identify candidate membrane transporters involved in transplacental nutrient transfer associated with PE/IUGR or GDM. METHODS: Using in silico strategies, we analysed various gene expression data sets generated on different platforms focusing on solute carriers, ABC transporters and TRP channels in order to identify transporters that are differently expressed between patients and gestational age-matched controls. These bioinformatic analyses were combined with literature data to define a catalogue of target transporters that could be involved in the development of PE/IUGR or GDM. Transporters of interest were then analysed for gene expression using qRT-PCR in placental tissues of patients and controls. For validating the results on protein and functional level, we started to establish an in vitro assay using freshly isolated primary cytotrophoblast cells polarized on the Transwell® system. RESULTS: Using bioinformatics approaches, we initially identified 37 target membrane proteins which were mainly associated with the transport of amino acids, vitamins, and trace elements. At the current state of analysis, the amino acid transporters SLC7A7, SLC38A2, SLC38A5, and the thiamine transporter SLC19A3 showed significant differences in placental mRNA expression between controls and patients affected by PE and/or IUGR. Subsequent gene expression analysis in our in-house GDM placental tissue bank is still ongoing. CONCLUSIONS: Based on our in silico analyses, literature data and first follow-up in vitro validations, we were able to define potentially interesting candidate transporters implicated in PE/IUGR or GDM. To date, additional newly defined candidate targets are being analysed on mRNA level in PE/IUGR and GDM. Subsequent analyses on protein and functional level will reveal whether these targets could be of diagnostic or therapeutical interest in these pregnancy pathologies.

miRNA-26b Overexpression in Ulcerative Colitis-associated Carcinogenesis

BACKGROUND Longstanding ulcerative colitis (UC) bears a high risk for development of UC-associated colorectal carcinoma (UCC). The inflammatory microenvironment influences microRNA expression, which in turn deregulates target gene expression. microRNA-26b (miR-26b) was shown to be instrumental in normal tissue growth and differentiation. Thus, we aimed to investigate the impact of miR-26b in inflammation-associated colorectal carcinogenesis. METHODS Two different cohorts of patients were investigated. In the retrospective group, a tissue microarray with 38 samples from 17 UC/UCC patients was used for miR-26b in situ hybridization and quantitative reverse transcription polymerase chain reaction analyses. In the prospective group, we investigated miR-26b expression in 25 fresh-frozen colon biopsies and corresponding serum samples of 6 UC and 15 non-UC patients, respectively. In silico analysis, Ago2-RNA immunoprecipitation, luciferase reporter assay, quantitative reverse transcription polymerase chain reaction examination, and miR-26b mimic overexpression were employed for target validation. RESULTS miR-26b expression was shown to be upregulated with disease progression in tissues and serum of UC and UCC patients. Using miR-26b and Ki-67 expression levels, an UCC was predicted with high accuracy. We identified 4 novel miR-26b targets (DIP1, MDM2, CREBBP, BRCA1). Among them, the downregulation of the E3 ubiquitin ligase DIP1 was closely related to death-associated protein kinase stabilization along the normal mucosa-UC-UCC sequence. In silico functional pathway analysis revealed that the common cellular pathways affected by miR-26b are highly related to cancerogenesis and the development of gastrointestinal diseases. CONCLUSIONS We suggest that miR-26b could serve as a biomarker for inflammation-associated processes in the gastrointestinal system. Because miR-26b expression is downregulated in sporadic colon cancer, it could discriminate between UCC and the sporadic cancer type.

Genome profiling of sterol synthesis shows convergent evolution in parasites and guides chemotherapeutic attack

Sterols are an essential class of lipids in eukaryotes, where they serve as structural components of membranes and play important roles as signaling molecules. Sterols are also of high pharmacological significance: cholesterol-lowering drugs are blockbusters in human health, and inhibitors of ergosterol biosynthesis are widely used as antifungals. Inhibitors of ergosterol synthesis are also being developed for Chagas's disease, caused by Trypanosoma cruzi. Here we develop an in silico pipeline to globally evaluate sterol metabolism and perform comparative genomics. We generate a library of hidden Markov model-based profiles for 42 sterol biosynthetic enzymes, which allows expressing the genomic makeup of a given species as a numerical vector. Hierarchical clustering of these vectors functionally groups eukaryote proteomes and reveals convergent evolution, in particular metabolic reduction in obligate endoparasites. We experimentally explore sterol metabolism by testing a set of sterol biosynthesis inhibitors against trypanosomatids, Plasmodium falciparum, Giardia, and mammalian cells, and by quantifying the expression levels of sterol biosynthetic genes during the different life stages of T. cruzi and Trypanosoma brucei. The phenotypic data correlate with genomic makeup for simvastatin, which showed activity against trypanosomatids. Other findings, such as the activity of terbinafine against Giardia, are not in agreement with the genotypic profile.

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.

Whole-exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow-up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We utilized whole-exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mutant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow-up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

Biomechanical model of tumor growth: application to the European ContraCancrum project

This abstract presents the biomechanical model that is used in the European ContraCancrum project, aiming at simulating tumor evolution in the brain and lung. The construction of the finite element model as well as a simulation of tumor growth are shown. The construction of the mesh is fully automatic and is therefore compatible with a clinical application. This biomechanical model will be later combined to a cellular level simulator also developed in the project.

Atomistic simulations of 2D bicomponent self-assembly: from molecular recognition to self-healing

Supramolecular two-dimensional engineering epitomizes the design of complex molecular architectures through recognition events in multicomponent self-assembly. Despite being the subject of in-depth experimental studies, such articulated phenomena have not been yet elucidated in time and space with atomic precision. Here we use atomistic molecular dynamics to simulate the recognition of complementary hydrogen-bonding modules forming 2D porous networks on graphite. We describe the transition path from the melt to the crystalline hexagonal phase and show that self-assembly proceeds through a series of intermediate states featuring a plethora of polygonal types. Finally, we design a novel bicomponent system possessing kinetically improved self-healing ability in silico, thus demonstrating that a priori engineering of 2D self-assembly is possible.

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohn's disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10 ? ? ). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohn's disease.

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.

Block of the hERG channel by bupivacaine: Electrophysiological and modeling insights towards stereochemical optimization

The hERG voltage-gated potassium channel mediates the cardiac I(Kr) current, which is crucial for the duration of the cardiac action potential. Undesired block of the channel by certain drugs may prolong the QT interval and increase the risk of malignant ventricular arrhythmias. Although the molecular determinants of hERG block have been intensively studied, not much is known about its stereoselectivity. Levo-(S)-bupivacaine was the first drug reported to have a higher affinity to block hERG than its enantiomer. This study strives to understand the principles underlying the stereoselectivity of bupivacaine block with the help of mutagenesis analyses and molecular modeling simulations. Electrophysiological measurements of mutated hERG channels allowed for the identification of residues involved in bupivacaine binding and stereoselectivity. Docking and molecular mechanics simulations for both enantiomers of bupivacaine and terfenadine (a non-stereoselective blocker) were performed inside an open-state model of the hERG channel. The predicted binding modes enabled a clear depiction of ligand-protein interactions. Estimated binding affinities for both enantiomers were consistent with electrophysiological measurements. A similar computational procedure was applied to bupivacaine enantiomers towards two mutated hERG channels (Tyr652Ala and Phe656Ala). This study confirmed, at the molecular level, that bupivacaine stereoselectively binds the hERG channel. These results help to lay the foundation for structural guidelines to optimize the cardiotoxic profile of drug candidates in silico.

Age-modulated association between prefrontal NAA and the BDNF gene

Brain-derived neurotrophic factor (BDNF) has been implicated in the pathophysiology of psychiatric and neurological disorders and in the mechanisms of antidepressant pharmacotherapy. Psychiatric and neurological conditions have also been associated with reduced brain levels of N-acetyl-aspartate (NAA), which has been used as a putative marker of neural integrity. However, few studies have explored the relationship between BDNF polymorphisms and NAA levels directly. Here, we present data from a single-voxel proton magnetic resonance spectroscopy study of 64 individuals and explore the relationship between BDNF polymorphisms and prefrontal NAA level. Our results indicate an association between a single nucleotide polymorphism (SNP) within BDNF, known as rs1519480, and reduced NAA level (p = 0.023). NAA levels were further predicted by age and Asian ancestry. There was a significant rs1519480 × age interaction on NAA level (p = 0.031). Specifically, the effect of rs1519480 on NAA level became significant at age ⩾34.17 yr. NAA level decreased with advancing age for genotype TT (p = 0.001) but not for genotype CT (p = 0.82) or CC (p = 0.34). Additional in silico analysis of 142 post-mortem brain samples revealed an association between the same SNP and reduced BDNF mRNA expression in the prefrontal cortex. The rs1519480 SNP influences BDNF mRNA expression and has an impact on prefrontal NAA level over time. This genetic mechanism may contribute to inter-individual variation in cognitive performance seen during normal ageing, as well as contributing to the risk for developing psychiatric and neurological conditions.

A human-horse comparative map based on equine BAC end sequences

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.