ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Patient-Oriented SCORAD (PO-SCORAD): a new self-assessment scale in atopic dermatitis validated in Europe

BACKGROund: Patient-oriented medicine is an emerging concept, encouraged by the World Health Organization, to greater involvement of the patient in the management of chronic diseases. The Patient-Oriented SCORing Atopic Dermatitis (PO-SCORAD) index is a self-assessment score allowing the patient to comprehensively evaluate the actual course of atopic dermatitis (AD), using subjective and objective criteria derived mainly from the SCORAD, a validated AD severity clinical assessment tool.

Impact of complications in total ankle replacement and ankle arthrodesis analyzed with a validated outcome measurement

Major modifications in the design and techniques of total ankle replacement have challenged the perception that ankle arthrodesis is the treatment of choice for end-stage ankle arthritis. High complication and revision rates have been reported after both procedures.

Validated capillary electrophoresis assay for the simultaneous enantioselective determination of propafenone and its major metabolites in biological samples

A robust, inexpensive, and fully validated CE method for the simultaneous determination of the enantiomers of propafenone (PPF), 5-hydroxy-propafenone (5OH-PPF) and N-despropyl-propafenone (NOR-PPF) in serum and in in vitro media is described. It is based upon liquid-liquid extraction at alkaline pH followed by analysis of the reconstituted extract by CE in presence of a pH 2.0 running buffer composed of 100 mM sodium phosphate, 19% methanol, and 0.6% highly sulfated beta-CD. For each compound, the S-enantiomers are shown to migrate ahead of their antipodes, and the overall run time is about 30 min. Enantiomer levels between 25 and 1000 ng/mL provide linear calibration graphs, and the LOD for all enantiomers is between 10 and 12 ng/mL. The assay is shown to be suitable for the determination of the enantiomers of PPF and its metabolites in in vitro incubations comprising human liver microsomes or single CYP450 enzymes (SUPERSOMES). Incubations with CYP2D6 SUPERSOMES revealed, for the first time, the simultaneous formation of the enantiomers of 5OH-PPF and NOR-PPF with that enzyme. CE data can be used for the evaluation of the enzymatic N-dealkylation and hydroxylation rates.

Radiographic analysis of femoroacetabular impingement with Hip2Norm-reliable and validated

The purpose of this study was to validate the accuracy, consistency, and reproducibility/reliability of a new method for correction of pelvic tilt and rotation of radiographic hip parameters for pincer type of femoroacetabular impingement on an anteroposterior pelvic radiograph. Thirty cadaver hips and 100 randomized, blinded AP pelvic radiographs were used for investigation. To detect the software accuracy, the calculated femoral head coverage and classic hip parameters determined with our software were compared to reference measurements based on CT scans or conventional radiographs in a neutral orientation as gold standard. To investigate software consistency, differences among the different parameters for each cadaver pelvis were calculated when reckoned back from a random to the neutral orientation. Intra- and interobserver comparisons were used to analyze the reliability and reproducibility of all parameters. All but two parameters showed a good-to-very good accuracy with the reference measurements. No relevant systematic errors were detected in the Bland-Altman analysis. Software consistency was good-to-very good for all parameters. A good-to-very good reliability and reproducibility was found for a substantial number of the evaluated radiographic acetabular parameters. The software appears to be an accurate, consistent, reliable, and reproducible method for analysis of acetabular pathomorphologies.

A validated assay by liquid chromatography-tandem mass spectrometry for the simultaneous quantification of elvitegravir and rilpivirine in HIV positive patients

Because of the large variability in the pharmacokinetics of anti-HIV drugs, therapeutic drug monitoring in patients may contribute to optimize the overall efficacy and safety of antiretroviral therapy. An LC–MS/MS method for the simultaneous assay in plasma of the novel antiretroviral agents rilpivirine (RPV) and elvitegravir (EVG) has been developed to that endeavor. Plasma samples (100 μL) extraction is performed by protein precipitation with acetonitrile, and the supernatant is subsequently diluted 1:1 with 20-mM ammonium acetate/MeOH 50:50. After reverse-phase chromatography, quantification of RPV and EVG, using matrix-matched calibration samples, is performed by electrospray ionization–triple quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. The stable isotopic-labeled compounds RPV-13C6 and EVG-D6 were used as internal standards. The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effects variability (<6.4%), as well as EVG and RPV short and long-term stability in plasma. Calibration curves were validated over the clinically relevant concentrations ranging from 5 to 2500 ng/ml for RPV and from 50 to 5000 ng/ml for EVG. The method is precise (inter-day CV%: 3–6.3%) and accurate (3.8–7.2%). Plasma samples were found to be stable (<15%) in all considered conditions (RT/48 h, +4°C/48 h, −20°C/3 months and 60°C/1 h). Selected metabolite profiles analysis in patients' samples revealed the presence of EVG glucuronide, that was well separated from parent EVG, allowing to exclude potential interferences through the in-source dissociation of glucuronide to parent drug. This new, rapid and robust LCMS/MS assay for the simultaneous quantification of plasma concentrations of these two major new anti-HIV drugs EVG and RPV offers an efficient analytical tool for clinical pharmacokinetics studies and routine therapeutic drug monitoring service.

Rib fractures at postmortem computed tomography (PMCT) validated against the autopsy

To evaluate the sensitivity of postmortem computed tomography (PMCT) in rib fracture detection validated against autopsy. Fifty-one forensic cases underwent a postmortem CT prior to forensic autopsy. Two image readers (radiologist and forensic pathologist) assessed high resolution CT data sets for rib fractures. Correct recognition rates (CRR), sensitivity and specificity values were calculated over all observations as well as individually for every rib and region. Additionally, for partial rib fractures the sensitivity of autopsy was calculated vice versa. 3876 entries in each study protocol (autopsy, PMCT radiologist and PMCT forensic pathologist) were investigated. A total of 690 fractures (autopsy), 491 (PMCT and radiologist) and 559 (PMCT and forensic pathologist) were detected. The CRR was 0.85. Sensitivity and specificity of PMCT for rib fracture detection were 0.63 (0.58 radiologist, 0.68 forensic pathologist) and 0.97 (both readers 0.97), respectively. Low CRR and sensitivity values were obtained for antero-lateral fractures. Partial rib fractures were better detected by PMCT. PMCT has a rather low sensitivity for rib fracture detection when validated against autopsy and indicates that clinical CT may also demonstrate a reasonable number of false negatives. Partial rib fractures often remain undetected at autopsy.

An experimentally validated finite element method for augmented vertebral bodies

Background Finite element models of augmented vertebral bodies require a realistic modelling of the cement infiltrated region. Most methods published so far used idealized cement shapes or oversimplified material models for the augmented region. In this study, an improved, anatomy-specific, homogenized finite element method was developed and validated to predict the apparent as well as the local mechanical behavior of augmented vertebral bodies. Methods Forty-nine human vertebral body sections were prepared by removing the cortical endplates and scanned with high-resolution peripheral quantitative CT before and after injection of a standard and a low-modulus bone cement. Forty-one specimens were tested in compression to measure stiffness, strength and contact pressure distributions between specimens and loading-plates. From the remaining eight, fourteen cylindrical specimens were extracted from the augmented region and tested in compression to obtain material properties. Anatomy-specific finite element models were generated from the CT data. The models featured element-specific, density-fabric-based material properties, damage accumulation, real cement distributions and experimentally determined material properties for the augmented region. Apparent stiffness and strength as well as contact pressure distributions at the loading plates were compared between simulations and experiments. Findings The finite element models were able to predict apparent stiffness (R2 > 0.86) and apparent strength (R2 > 0.92) very well. Also, the numerically obtained pressure distributions were in reasonable quantitative (R2 > 0.48) and qualitative agreement with the experiments. Interpretation The proposed finite element models have proven to be an accurate tool for studying the apparent as well as the local mechanical behavior of augmented vertebral bodies.

Coreleased orexin and glutamate evoke nonredundant spike outputs and computations in histamine neurons.

Stable wakefulness requires orexin/hypocretin neurons (OHNs) and OHR2 receptors. OHNs sense diverse environmental cues and control arousal accordingly. For unknown reasons, OHNs contain multiple excitatory transmitters, including OH peptides and glutamate. To analyze their cotransmission within computational frameworks for control, we optogenetically stimulated OHNs and examined resulting outputs (spike patterns) in a downstream arousal regulator, the histamine neurons (HANs). OHR2s were essential for sustained HAN outputs. OHR2-dependent HAN output increased linearly during constant OHN input, suggesting that the OHN→HAN(OHR2) module may function as an integral controller. OHN stimulation evoked OHR2-dependent slow postsynaptic currents, similar to midnanomolar OH concentrations. Conversely, glutamate-dependent output transiently communicated OHN input onset, peaking rapidly then decaying alongside OHN→HAN glutamate currents. Blocking glutamate-driven spiking did not affect OH-driven spiking and vice versa, suggesting isolation (low cross-modulation) of outputs. Therefore, in arousal regulators, cotransmitters may translate distinct features of OHN activity into parallel, nonredundant control signals for downstream effectors.

Towards the continuum limit in transport coefficient computations

The analytic continuation needed for the extraction of transport coefficients necessitates in principle a continuous function of the Euclidean time variable. We report on progress towards achieving the continuum limit for 2-point correlator measurements in thermal SU(3) gauge theory, with specific attention paid to scale setting. In particular, we improve upon the determination of the critical lattice coupling and the critical temperature of pure SU(3) gauge theory, estimating r0Tc ≃ 0.7470(7) after a continuum extrapolation. As an application the determination of the heavy quark momentum diffusion coefficient from a correlator of colour-electric fields attached to a Polyakov loop is discussed.

Accurate computations of the structures and binding energies of the imidazole ... benzene and pyrrole ... benzene complexes

Using explicitly-correlated coupled-cluster theory with single and double excitations, the intermolecular distances and interaction energies of the T-shaped imidazole⋯⋯benzene and pyrrole⋯⋯benzene complexes have been computed in a large augmented correlation-consistent quadruple-zeta basis set, adding also corrections for connected triple excitations and remaining basis-set-superposition errors. The results of these computations are used to assess other methods such as Møller–Plesset perturbation theory (MP2), spin-component-scaled MP2 theory, dispersion-weighted MP2 theory, interference-corrected explicitly-correlated MP2 theory, dispersion-corrected double-hybrid density-functional theory (DFT), DFT-based symmetry-adapted perturbation theory, the random-phase approximation, explicitly-correlated ring-coupled-cluster-doubles theory, and double-hybrid DFT with a correlation energy computed in the random-phase approximation.