Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA.

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

BioMOBS: A multi-omics visual analytics workflow for biomolecular insight generation.

One of the challenges in multi-omics data analysis for precision medicine is the efficient exploration of undiscovered molecular interactions in disease processes. We present BioMOBS, a workflow consisting of two data visualization tools integrated with an open-source molecular information database to perform clinically relevant analyses (https://github.com/driesheylen123/BioMOBS). We performed exploratory pathway analysis with BioMOBS and demonstrate its ability to generate relevant molecular hypotheses, by reproducing recent findings in type 2 diabetes UK biobank data. The central visualisation tool, where data-driven and literature-based findings can be integrated, is available within the github link as well. BioMOBS is a workflow that leverages information from multiple data-driven interactive analyses and visually integrates it with established pathway knowledge. The demonstrated use cases place trust in the usage of BioMOBS as a procedure to offer clinically relevant insights in disease pathway analyses on various types of omics data.

Intravesicular Genomic DNA Enriched by Size Exclusion Chromatography Can Enhance Lung Cancer Oncogene Mutation Detection Sensitivity.

Extracellular vesicles (EVs) are cell-derived structures surrounded by a lipid bilayer that carry RNA and DNA as potential templates for molecular diagnostics, e.g., in cancer genotyping. While it has been established that DNA templates appear on the outside of EVs, no consensus exists on which nucleic acid species inside small EVs (<200 nm, sEVs) are sufficiently abundant and accessible for developing genotyping protocols. We investigated this by extracting total intravesicular nucleic acid content from sEVs isolated from the conditioned cell medium of the human NCI-H1975 cell line containing the epidermal growth factor (EGFR) gene mutation T790M as a model system for non-small cell lung cancer. We observed that mainly short genomic DNA (<35−100 bp) present in the sEVs served as a template. Using qEV size exclusion chromatography (SEC), significantly lower yield and higher purity of isolated sEV fractions were obtained as compared to exoEasy membrane affinity purification and ultracentrifugation. Nevertheless, we detected the EGFR T790M mutation in the sEVs' lumen with similar sensitivity using digital PCR. When applying SEC-based sEV separation prior to cell-free DNA extraction on spiked human plasma samples, we found significantly higher mutant allele frequencies as compared to standard cell-free DNA extraction, which in part was due to co-purification of circulating tumor DNA. We conclude that intravesicular genomic DNA can be exploited next to ctDNA to enhance EGFR T790M mutation detection sensitivity by adding a fast and easy-to-use sEV separation method, such as SEC, upstream of standard clinical cell-free DNA workflows.

Towards Building a Quantitative Proteomics Toolbox in Precision Medicine: A Mini-Review.

Precision medicine as a framework for disease diagnosis, treatment, and prevention at the molecular level has entered clinical practice. From the start, genetics has been an indispensable tool to understand and stratify the biology of chronic and complex diseases in precision medicine. However, with the advances in biomedical and omics technologies, quantitative proteomics is emerging as a powerful technology complementing genetics. Quantitative proteomics provide insight about the dynamic behaviour of proteins as they represent intermediate phenotypes. They provide direct biological insights into physiological patterns, while genetics accounting for baseline characteristics. Additionally, it opens a wide range of applications in clinical diagnostics, treatment stratification, and drug discovery. In this mini-review, we discuss the current status of quantitative proteomics in precision medicine including the available technologies and common methods to analyze quantitative proteomics data. Furthermore, we highlight the current challenges to put quantitative proteomics into clinical settings and provide a perspective to integrate proteomics data with genomics data for future applications in precision medicine.

Constrained Standardization of Count Data from Massive Parallel Sequencing.

In high-throughput omics disciplines like transcriptomics, researchers face a need to assess the quality of an experiment prior to an in-depth statistical analysis. To efficiently analyze such voluminous collections of data, researchers need triage methods that are both quick and easy to use. Such a normalization method for relative quantitation, CONSTANd, was recently introduced for isobarically-labeled mass spectra in proteomics. It transforms the data matrix of abundances through an iterative, convergent process enforcing three constraints: (I) identical column sums; (II) each row sum is fixed (across matrices) and (III) identical to all other row sums. In this study, we investigate whether CONSTANd is suitable for count data from massively parallel sequencing, by qualitatively comparing its results to those of DESeq2. Further, we propose an adjustment of the method so that it may be applied to identically balanced but differently sized experiments for joint analysis. We find that CONSTANd can process large data sets at well over 1 million count records per second whilst mitigating unwanted systematic bias and thus quickly uncovering the underlying biological structure when combined with a PCA plot or hierarchical clustering. Moreover, it allows joint analysis of data sets obtained from different batches, with different protocols and from different labs but without exploiting information from the experimental setup other than the delineation of samples into identically processed sets (IPSs). CONSTANd's simplicity and applicability to proteomics as well as transcriptomics data make it an interesting candidate for integration in multi-omics workflows.

CONSTANd: An Efficient Normalization Method for Relative Quantification in Small- and Large-Scale Omics Experiments in R BioConductor and Python.

For differential expression studies in all omics disciplines, data normalization is a crucial step that is often subject to a balance between speed and effectiveness. To keep up with the data produced by high-throughput instruments, researchers require fast and easy-to-use yet effective methods that fit into automated analysis pipelines. The CONSTANd normalization method meets these criteria, so we have made its source code available for R/BioConductor and Python. We briefly review the method and demonstrate how it can be used in different omics contexts for experiments of any scale. Widespread adoption across omics disciplines would ease data integration in multiomics experiments.

PRiSM: A prototype for exhaustive, restriction-free database searching for mass spectrometry-based identification.

RATIONALE: The current methods for identifying peptides in mass spectral product ion data still struggle to do so for the majority of spectra. Based on the experimental setup and other assumptions, such methods restrict the search space to speed up computations, but at the cost of creating blind spots. The proteomics community would greatly benefit from a method that is capable of covering the entire search space without using any restrictions, thus establishing a baseline for identification. METHODS: We conceived the "mass pattern paradigm" (MPP) that enables the creation of such an identification method, and we implemented it into a prototype database search engine "PRiSM" (PRotein-Spectrum Matching). We then assessed its operational characteristics by applying it to publicly available high-precision mass spectra of low and high identification difficulty. We used those characteristics to gain theoretical insights into trade-offs between sensitivity and speed when trying to establish a baseline for identification. RESULTS: Of 100 low difficulty spectra, PRiSM and SEQUEST agree on 84 identifications (of which 75 are statistically significant). Of 15 of 100 spectra not identified in a previous study (using SEQUEST), 13 are considered reliable after visual inspection and represent 3 proteins (out of 9 in total) not detected previously. CONCLUSIONS: Despite leaving noise intact, the simple PRiSM prototype can make statistically reliable identifications, while controlling the false discovery rate by fitting a null distribution. It also identifies some spectra previously unidentifiable in an "extremely open" SEQUEST search, paving the way to establishing a baseline for identification in proteomics.

Author Correction: Leveraging European infrastructures to access 1 million human genomes by 2022.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Leveraging European infrastructures to access 1 million human genomes by 2022.

Human genomics is undergoing a step change from being a predominantly research-driven activity to one driven through health care as many countries in Europe now have nascent precision medicine programmes. To maximize the value of the genomic data generated, these data will need to be shared between institutions and across countries. In recognition of this challenge, 21 European countries recently signed a declaration to transnationally share data on at least 1 million human genomes by 2022. In this Roadmap, we identify the challenges of data sharing across borders and demonstrate that European research infrastructures are well-positioned to support the rapid implementation of widespread genomic data access.

Enhancing the Performance of DNA Surface-Hybridization Biosensors through Target Depletion.

DNA surface-hybridization biosensors utilize the selective hybridization of target sequences in solution to surface-immobilized probes. In this process, the target is usually assumed to be in excess, so that its concentration does not significantly vary while hybridizing to the surface-bound probes. If the target is initially at low concentrations and/or if the number of probes is very large, and they have high affinity for the target, the DNA in solution may become depleted. In this paper we analyze the equilibrium and kinetics of hybridization of DNA biosensors in the case of strong target depletion, by extending the Langmuir adsorption model. We focus, in particular, on the detection of a small amount of a single-nucleotide "mutant" sequence (concentration c2) in a solution, which differs by one or more nucleotides from an abundant "wild-type" sequence (concentration c1 ≫ c2). We show that depletion can give rise to a strongly enhanced sensitivity of the biosensors. Using representative values of rate constants and hybridization free energies, we find that in the depletion regime one could detect relative concentrations c2/c1 that are up to 3 orders of magnitude smaller than in the conventional approach. The kinetics is surprisingly rich and exhibits a nonmonotonic adsorption with no counterpart in the no-depletion case. Finally, we show that, alongside enhanced detection sensitivity, this approach offers the possibility of sample enrichment, by substantially increasing the relative amount of the mutant over the wild-type sequence.

MIND: A Double-Linear Model To Accurately Determine Monoisotopic Precursor Mass in High-Resolution Top-Down Proteomics.

Top-down proteomics approaches are becoming ever more popular, due to the advantages offered by knowledge of the intact protein mass in correctly identifying the various proteoforms that potentially arise due to point mutation, alternative splicing, post-translational modifications, etc. Usually, the average mass is used in this context; however, it is known that this can fluctuate significantly due to both natural and technical causes. Ideally, one would prefer to use the monoisotopic precursor mass, but this falls below the detection limit for all but the smallest proteins. Methods that predict the monoisotopic mass based on the average mass are potentially affected by imprecisions associated with the average mass. To address this issue, we have developed a framework based on simple, linear models that allows prediction of the monoisotopic mass based on the exact mass of the most-abundant (aggregated) isotope peak, which is a robust measure of mass, insensitive to the aforementioned natural and technical causes. This linear model was tested experimentally, as well as in silico, and typically predicts monoisotopic masses with an accuracy of only a few parts per million. A confidence measure is associated with the predicted monoisotopic mass to handle the off-by-one-Da prediction error. Furthermore, we introduce a correction function to extract the "true" (i.e., theoretically) most-abundant isotope peak from a spectrum, even if the observed isotope distribution is distorted by noise or poor ion statistics. The method is available online as an R shiny app: https://valkenborg-lab.shinyapps.io/mind/.

Role of nanoparticle size and sialic acids in the distinct time-evolution profiles of nanoparticle uptake in hematopoietic progenitor cells and monocytes.

BACKGROUND: Human hematopoietic progenitor cells (HPCs) are important for cell therapy in cancer and tissue regeneration. In vitro studies have shown a transient association of 40 nm polystyrene nanoparticles (PS NPs) with these cells, which is of interest for intelligent design and application of NPs in HPC-based regenerative protocols. In this study, we aimed to investigate the involvement of nanoparticles' size and membrane-attached glycan molecules in the interaction of HPCs with PS NPs, and compared it with monocytes. Human cord blood-derived HPCs and THP-1 cells were exposed to fluorescently labelled, carboxylated PS NPs of 40, 100 and 200 nm. Time-dependent nanoparticle membrane association and/or uptake was observed by measuring fluorescence intensity of exposed cells at short time intervals using flow cytometry. By pretreating the cells with neuraminidase, we studied the possible effect of membrane-associated sialic acids in the interaction with NPs. Confocal microscopy was used to visualize the cell-specific character of the NP association. RESULTS: Confocal images revealed that the majority of PS NPs was initially observed to be retained at the outer membrane of HPCs, while the same NPs showed immediate internalization by THP-1 monocytic cells. After prolonged exposure up to 4 h, PS NPs were also observed to enter the HPCs' intracellular compartment. Cell-specific time courses of NP association with HPCs and THP-1 cells remained persistent after cells were enzymatically treated with neuraminidase, but significantly increased levels of NP association could be observed, suggesting a role for membrane-associated sialic acids in this process. CONCLUSIONS: We conclude that the terminal membrane-associated sialic acids contribute to the NP retention at the outer cell membrane of HPCs. This retention behavior is a unique characteristic of the HPCs and is independent of NP size.

QCQuan: A Web Tool for the Automated Assessment of Protein Expression and Data Quality of Labeled Mass Spectrometry Experiments.

In the context of omics disciplines and especially proteomics and biomarker discovery, the analysis of a clinical sample using label-based tandem mass spectrometry (MS) can be affected by sample preparation effects or by the measurement process itself, resulting in an incorrect outcome. Detection and correction of these mistakes using state-of-the-art methods based on mixed models can use large amounts of (computing) time. MS-based proteomics laboratories are high-throughput and need to avoid a bottleneck in their quantitative pipeline by quickly discriminating between high- and low-quality data. To this end we developed an easy-to-use web-tool called QCQuan (available at qcquan.net ) which is built around the CONSTANd normalization algorithm. It automatically provides the user with exploratory and quality control information as well as a differential expression analysis based on conservative, simple statistics. In this document we describe in detail the scientifically relevant steps that constitute the workflow and assess its qualitative and quantitative performance on three reference data sets. We find that QCQuan provides clear and accurate indications about the scientific value of both a high- and a low-quality data set. Moreover, it performed quantitatively better on a third data set than a comparable workflow assembled using established, reliable software.

Molecular drug susceptibility testing and strain typing of tuberculosis by DNA hybridization.

In Mycobacterium tuberculosis (Mtb) the detection of single nucleotide polymorphisms (SNPs) is of high importance both for diagnostics, since drug resistance is primarily caused by the acquisition of SNPs in multiple drug targets, and for epidemiological studies in which strain typing is performed by SNP identification. To provide the necessary coverage of clinically relevant resistance profiles and strain types, nucleic acid-based measurement techniques must be able to detect a large number of potential SNPs. Since the Mtb problem is pressing in many resource-poor countries, requiring low-cost point-of-care biosensors, this is a non-trivial technological challenge. This paper presents a proof-of-concept in which we chose simple DNA-DNA hybridization as a sensing principle since this can be transferred to existing low-cost hardware platforms, and we pushed the multiplex boundaries of it. With a custom designed probe set and a physicochemical-driven data analysis it was possible to simultaneously detect the presence of SNPs associated with first- and second-line drug resistance and Mtb strain typing. We have demonstrated its use for the identification of drug resistance and strain type from a panel of phylogenetically diverse clinical strains. Furthermore, reliable detection of the presence of a minority population (<5%) of drug-resistant Mtb was possible.

DNA capture into the ClyA nanopore: diffusion-limited versus reaction-limited processes.

The capture and translocation of biomolecules through nanometer-scale pores are processes with a potentially large number of applications, and hence they have been intensively studied in recent years. The aim of this paper is to review existing models of the capture process by a nanopore, together with some recent experimental data of short single- and double-stranded DNA captured by the Cytolysin A (ClyA) nanopore. ClyA is a transmembrane protein of bacterial origin which has been recently engineered through site-specific mutations, to allow the translocation of double- and single-stranded DNA. A comparison between theoretical estimations and experiments suggests that for both cases the capture is a reaction-limited process. This is corroborated by the observed salt dependence of the capture rate, which we find to be in quantitative agreement with the theoretical predictions.

Thermodynamic framework to assess low abundance DNA mutation detection by hybridization.

The knowledge of genomic DNA variations in patient samples has a high and increasing value for human diagnostics in its broadest sense. Although many methods and sensors to detect or quantify these variations are available or under development, the number of underlying physico-chemical detection principles is limited. One of these principles is the hybridization of sample target DNA versus nucleic acid probes. We introduce a novel thermodynamics approach and develop a framework to exploit the specific detection capabilities of nucleic acid hybridization, using generic principles applicable to any platform. As a case study, we detect point mutations in the KRAS oncogene on a microarray platform. For the given platform and hybridization conditions, we demonstrate the multiplex detection capability of hybridization and assess the detection limit using thermodynamic considerations; DNA containing point mutations in a background of wild type sequences can be identified down to at least 1% relative concentration. In order to show the clinical relevance, the detection capabilities are confirmed on challenging formalin-fixed paraffin-embedded clinical tumor samples. This enzyme-free detection framework contains the accuracy and efficiency to screen for hundreds of mutations in a single run with many potential applications in molecular diagnostics and the field of personalised medicine.

Direct detection of nano-scale extracellular vesicles derived from inflammation-triggered endothelial cells using surface plasmon resonance.

A major conceptual breakthrough in cell signaling has been the finding of EV as new biomarker shuttles in body fluids. Now, one of the major challenges in using these nanometer-sized biological entities as diagnostic marker is the development of translational methodologies to profile them. SPR offers a promising label-free and real time platform with a high potential for biomarker detection. Therefore, we aimed to develop a uniform SPR methodology to detect specific surface markers on EV derived from patient with CHD. EVs having an approximate size range between 30 and 100 nm (~48.5%) and 100-300 nm (~51.5%) were successfully isolated. The biomarker profile of EV was verified using immunogold labeling, ELISA and SPR. Using SPR, we demonstrated an increased binding of EV derived from patients with CHD to anti-ICAM-1 antibodies as compared to EV from healthy donors. Our current findings open up novel opportunities for in-depth and label-free investigation of EV.

Transient loading of CD34+ hematopoietic progenitor cells with polystyrene nanoparticles.

CD34+ hematopoietic progenitor cells (HPCs) offer great opportunities to develop new treatments for numerous malignant and non-malignant diseases. Nanoparticle (NP)-based strategies can further enhance this potential, and therefore a thorough understanding of the loading behavior of HPCs towards NPs is essential for a successful application. The present study focusses on the interaction kinetics of 40 nm sized carboxylated polystyrene (PS) NPs with HPCs. Interestingly, a transient association of the NPs with HPCs is observed, reaching a maximum within 1 hour and declining afterwards. This behavior is not seen in dendritic cells (CD34-DCs) differentiated from HPCs, which display a monotonic increase in NP load. We demonstrate that this transient interaction requires an energy-dependent cellular process, suggesting active loading and release of NPs by HPCs. This novel observation offers a unique approach to transiently equip HPCs. A simple theoretical approach modeling the kinetics of NP loading and release is presented, contributing to a framework of describing this phenomenon.

Nonequilibrium Hybridization Enables Discrimination of a Point Mutation within 5-40 °C.

Detection of point mutations and single nucleotide polymorphisms in DNA and RNA has a growing importance in biology, biotechnology, and medicine. For the application at hand, hybridization assays are often used. Traditionally, they differentiate point mutations only at elevated temperatures (>40 °C) and in narrow intervals (ΔT = 1-10 °C). The current study demonstrates that a specially designed multistranded DNA probe can differentiate point mutations in the range of 5-40 °C. This unprecedentedly broad ambient-temperature range is enabled by a controlled combination of (i) nonequilibrium hybridization conditions and (ii) a mismatch-induced increase of equilibration time in respect to that of a fully matched complex, which we dub "kinetic inversion".

CONSTANd : A Normalization Method for Isobaric Labeled Spectra by Constrained Optimization.

In quantitative proteomics applications, the use of isobaric labels is a very popular concept as they allow for multiplexing, such that peptides from multiple biological samples are quantified simultaneously in one mass spectrometry experiment. Although this multiplexing allows that peptide intensities are affected by the same amount of instrument variability, systematic effects during sample preparation can also introduce a bias in the quantitation measurements. Therefore, normalization methods are required to remove this systematic error. At present, a few dedicated normalization methods for isobaric labeled data are at hand. Most of these normalization methods include a framework for statistical data analysis and rely on ANOVA or linear mixed models. However, for swift quality control of the samples or data visualization a simple normalization technique is sufficient. To this aim, we present a new and easy-to-use data-driven normalization method, named CONSTANd. The CONSTANd method employs constrained optimization and prior information about the labeling strategy to normalize the peptide intensities. Further, it allows maintaining the connection to any biological effect while reducing the systematic and technical errors. As a result, peptides can not only be compared directly within a multiplexed experiment, but are also comparable between other isobaric labeled datasets from multiple experimental designs that are normalized by the CONSTANd method, without the need to include a reference sample in every experimental setup. The latter property is especially useful when more than six, eight or ten (TMT/iTRAQ) biological samples are required to detect differential peptides with sufficient statistical power and to optimally make use of the multiplexing capacity of isobaric labels.