Sharing data from molecular simulations

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations has become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 (https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way toward more open, interoperable, and reproducible outputs coming from research studies using MD simulations

Endoscopic mucosal resection (EMR) versus endoscopic submucosal dissection (ESD) for resection of large distal non-pedunculated colorectal adenomas (MATILDA-trial): Rationale and design of a multicenter randomized clinical trial

Background: Endoscopic mucosal resection (EMR) is currently the most used technique for resection of large distal colorectal polyps. However, in large lesions EMR can often only be performed in a piecemeal fashion resulting in relatively low radical (R0)-resection rates and high recurrence rates. Endoscopic submucosal dissection (ESD) is a newer procedure that is more difficult resulting in a longer procedural time, but is promising due to the high en-bloc resection rates and the very low recurrence rates. We aim to evaluate the (cost-)effectiveness of ESD against EMR on both short (i.e. 6 months) and long-term (i.e. 36 months). We hypothesize that in the short-run ESD is more time consuming resulting in higher healthcare costs, but is (cost-) effective on the long-term due to lower patients burden, a higher number of R0-resections and lower recurrence rates with less need for repeated procedures. Methods: This is a multicenter randomized clinical trial in patients with a non-pedunculated polyp larger than 20 mm in the rectum, sigmoid, or descending colon suspected to be an adenoma by means of endoscopic assessment. Primary endpoint is recurrence rate at follow-up colonoscopy at 6 months. Secondary endpoints are R0-resection rate, perceived burden and quality of life, healthcare resources utilization and costs, surgical referral rate, complication rate and recurrence rate at 36 months. Quality-adjusted-life-year (QALY) will be estimated taking an area under the curve approach and using EQ-5D-indexes. Healthcare costs will be calculated by multiplying used healthcare services with unit prices. The cost-effectiveness of ESD against EMR will be expressed as incremental cost-effectiveness ratios (ICER) showing additional costs per recurrence free patient and as ICER showing additional costs per QALY. Discussion: If this trial confirms ESD to be favorable on the long-term, the burden of extra colonoscopies and repeated procedures can be prevented for future patients. Trial registration:NCT02657044(Clinicaltrials.gov), registered January 8, 2016

Outcome of the First wwPDB Hybrid / Integrative Methods Task Force Workshop

Structures of biomolecular systems are increasingly computed by integrative modeling that relies on varied types of experimental data and theoretical information. We describe here the proceedings and conclusions from the first wwPDB Hybrid/Integrative Methods Task Force Workshop held at the European Bioinformatics Institute in Hinxton, UK, on October 6 and 7, 2014. At the workshop, experts in various experimental fields of structural biology, experts in integrative modeling and visualization, and experts in data archiving addressed a series of questions central to the future of structural biology. How should integrative models be represented? How should the data and integrative models be validated? What data should be archived? How should the data and models be archived? What information should accompany the publication of integrative models

Determination of the Form Factors for the Decay B0 --> D*-l+nu_l and of the CKM Matrix Element |Vcb|

We present a combined measurement of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{cb}|$ and of the parameters $\rho^2$, $R_1$, and $R_2$, which fully characterize the form factors of the $B^0 \to D^{*-}\ell^{+}\nu_\ell$ decay in the framework of HQET, based on a sample of about 52,800 $B^0 \to D^{*-}\ell^{+}\nu_\ell$ decays recorded by the BABAR detector. The kinematical information of the fully reconstructed decay is used to extract the following values for the parameters (where the first errors are statistical and the second systematic): $\rho^2 = 1.156 \pm 0.094 \pm 0.028$, $R_1 = 1.329 \pm 0.131 \pm 0.044$, $R_2 = 0.859 \pm 0.077 \pm 0.022$, $\mathcal{F}(1)|V_{cb}| = (35.03 \pm 0.39 \pm 1.15) \times 10^{-3}$. By combining these measurements with the previous BABAR measurements of the form factors which employs a different technique on a partial sample of the data, we improve the statistical accuracy of the measurement, obtaining: $\rho^2 = 1.179 \pm 0.048 \pm 0.028, R_1 = 1.417 \pm 0.061 \pm 0.044, R_2 = 0.836 \pm 0.037 \pm 0.022,$ and $\mathcal{F}(1)|V_{cb}| = (34.68 \pm 0.32 \pm 1.15) \times 10^{-3}.$ Using the lattice calculations for the axial form factor $\mathcal{F}(1)$, we extract $|V_{cb}| =(37.74 \pm 0.35 \pm 1.25 \pm ^{1.23}_{1.44}) \times 10^{-3}$, where the third error is due to the uncertainty in $\mathcal{F}(1)$

Study of the Exclusive Initial-State Radiation Production of the DDˉD \bar D System

A study of exclusive production of the $D \bar D$ system through initial-state r adiation is performed in a search for charmonium states, where $D=D^0$ or $D^+$. The $D^0$ mesons are reconstructed in the $D^0 \to K^- \pi^+$, $D^0 \to K^- \pi^+ \pi^0$, and $D^0 \to K^- \pi^+ \pi^+ \pi^-$ decay modes. The $D^+$ is reconstructed through the $D^+ \to K^- \pi^+ \pi^+$ decay mode. The analysis makes use of an integrated luminosity of 288.5 fb$^{-1}$ collected by the BaBar experiment. The $D \bar D$ mass spectrum shows a clear $\psi(3770)$ signal. Further structures appear in the 3.9 and 4.1 GeV/$c^2$ regions. No evidence is found for Y(4260) decays to $D \bar D$, implying an up per limit \frac{\BR(Y(4260)\to D \bar D)}{\BR(Y(4260)\to J/\psi \pi^+ \pi^-)} < 7.6 (95 % confidence level)

A search for the decay B+→K+ννˉB^+ \to K^+ \nu \bar{\nu}

We search for the rare flavor-changing neutral-current decay $B^+ \to K^+ \nu \bar{\nu}$ in a data sample of 82 fb$^{-1}$ collected with the {\sl BABAR} detector at the PEP-II B-factory. Signal events are selected by examining the properties of the system recoiling against either a reconstructed hadronic or semileptonic charged-B decay. Using these two independent samples we obtain a combined limit of ${\mathcal B}(B^+ \to K^+ \nu \bar{\nu})<5.2 \times 10^{-5}$ at the 90% confidence level. In addition, by selecting for pions rather than kaons, we obtain a limit of ${\mathcal B}(B^+ \to \pi^+ \nu \bar{\nu})<1.0 \times 10^{-4}$ using only the hadronic B reconstruction method.Comment: 7 pages, 8 postscript figures, submitted to Phys. Rev. Let

High-reflectivity broadband distributed Bragg reflector lattice matched to ZnTe

We report on the realization of a high quality distributed Bragg reflector with both high and low refractive index layers lattice matched to ZnTe. Our structure is grown by molecular beam epitaxy and is based on binary compounds only. The high refractive index layer is made of ZnTe, while the low index material is made of a short period triple superlattice containing MgSe, MgTe, and ZnTe. The high refractive index step of Delta_n=0.5 in the structure results in a broad stopband and the reflectivity coefficient exceeding 99% for only 15 Bragg pairs.Comment: 4 pages, 3 figure