Search

Profiling of circulating tumor DNA (ctDNA) may offer a non-invasive approach to monitor disease progression. Here, we develop a quantitative method, exploiting local tissue-specific cell-free DNA (cfDNA) degradation patterns, that accurately estimates ctDNA burden independent of genomic aberrations. Nucleosome-dependent cfDNA degradation at promoters and first exon-intron junctions is strongly associated with differential transcriptional activity in tumors and blood. A quantitative model, based on just 6 regulatory regions, could accurately predict ctDNA levels in colorectal cancer patients. Strikingly, a model restricted to blood-specific regulatory regions could predict ctDNA levels across both colorectal and breast cancer patients. Using compact targeted sequencing (, Circulating tumour DNA (ctDNA) represents a non-invasive option to monitor cancer progression. Here, the authors perform deep sequencing of plasma cell-free DNA, and find that nucleosome-dependent cfDNA degradation at 6 specific regulatory regions is predictive of ctDNA burden.

Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure‐based and physics‐based atomistic force field with an efficient sampling strategy is adopted to simulate a model di‐domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low‐energy structures and the minimum‐size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small‐angle X‐ray scattering data. It is illustrated that the regularizations of energy and ensemble‐size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high‐energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure‐ensemble optimizations with a topology‐based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates. MOE (Min. of Education, S’pore) Accepted version

An ensemble-modeling scheme incorporating coarse-grained simulations with experimental small-angle X-ray scattering (SAXS) data is applied to dengue virus 2 (DENV2) nonstructural protein 5 (NS5). NS5 serves a key role in viral replication through its two domains that are connected by a 10-residue polypeptide segment. A set of representative structures is generated from a simulated structure pool using SAXS data fitting by the non-negativity least squares (NNLS) or standard ensemble optimization method (EOM) based on a genetic algorithm (GA). It is found that a proper low-energy threshold of the structure pool is necessary to produce a conformational ensemble of two representative structures by both NNLS and GA that agrees well with the experimental SAXS profile. The stability of the constructed ensemble is validated also by molecular dynamics simulations with an all-atom force field. The constructed ensemble successfully revealed the domain–domain orientation and domain-contacting interface of DENV2 NS5. Using experimental data fitting and additional investigations with synthesized data, it is found that energy restraint on the conformational pool is necessary to avoid overinterpretation of experimental data by spurious conformational representations. MOE (Min. of Education, S’pore) Accepted version