Monte Carlo algorithms based on the number of potential moves

We discuss Monte Carlo dynamics based on _E, the (microcanonical) average number of potential moves which increase the energy by Delta E in a single spin flip. The microcanonical average can be sampled using Monte Carlo dynamics of a single spin flip with a transition rate min(1, _E' / _E) from energy E to E'. A cumulative average (over Monte Carlo steps) can be used as a first approximation to the exact microcanonical average in the flip rate. The associated histogram is a constant independent of the energy. The canonical distribution of energy can be obtained from the transition matrix Monte Carlo dynamics. This second dynamics has fast relaxation time - at the critical temperature the relaxation time is proportional to specific heat. The dynamics are useful in connection with reweighting methods for computing thermodynamic quantities.Comment: 8 pages, 4 figures, invited talk on CCP99 conference, 20-27 May 99, Atlanta, G

Heat Conduction in two-dimensional harmonic crystal with disorder

We study the problem of heat conduction in a mass-disordered two-dimensional harmonic crystal. Using two different stochastic heat baths, we perform simulations to determine the system size (L) dependence of the heat current (J). For white noise heat baths we find that J ~ 1/L^a with $a \approx 0.59$ while correlated noise heat baths gives $a \approx 0.51$. A special case with correlated disorder is studied analytically and gives a=3/2 which agrees also with results from exact numerics.Comment: Revised version. 4 pages, 3 figure

A new approach to the study of the ground-state properties of 2D Ising spin glass

A new approach known as flat histogram method is used to study the +/-J Ising spin glass in two dimensions. Temperature dependence of the energy, the entropy, and other physical quantities can be easily calculated and we give the results for the zero-temperature limit. For the ground-state energy and entropy of an infinite system size, we estimate e0 = -1.4007 +/- 0.0085 and s0 = 0.0709 +/- 0.006, respectively. Both of them agree well with previous calculations. The time to find the ground-states as well as the tunneling times of the algorithm are also reported and compared with other methods.Comment: 11 pages, 4 figure

Flat histogram simulation of lattice polymer systems

We demonstrate the use of a new algorithm called the Flat Histogram sampling algorithm for the simulation of lattice polymer systems. Thermodynamics properties, such as average energy or entropy and other physical quantities such as end-to-end distance or radius of gyration can be easily calculated using this method. Ground-state energy can also be determined. We also explore the accuracy and limitations of this method. Key words: Monte Carlo algorithms, flat histogram sampling, HP model, lattice polymer systemsComment: 7 RevTeX two-column page

Identification of a Blood-Based Protein Biomarker Panel for Lung Cancer Detection

Lung cancer is the deadliest cancer worldwide, mainly due to its advanced stage at the time of diagnosis. A non-invasive method for its early detection remains mandatory to improve patients’ survival. Plasma levels of 351 proteins were quantified by Liquid Chromatography-Parallel Reaction Monitoring (LC-PRM)-based mass spectrometry in 128 lung cancer patients and 93 healthy donors. Bootstrap sampling and least absolute shrinkage and selection operator (LASSO) penalization were used to find the best protein combination for outcome prediction. The PanelomiX platform was used to select the optimal biomarker thresholds. The panel was validated in 48 patients and 49 healthy volunteers. A 6-protein panel clearly distinguished lung cancer from healthy individuals. The panel displayed excellent performance: area under the receiver operating characteristic curve (AUC) = 0.999, positive predictive value (PPV) = 0.992, negative predictive value (NPV) = 0.989, specificity = 0.989 and sensitivity = 0.992. The panel detected lung cancer independently of the disease stage. The 6-protein panel and other sub-combinations displayed excellent results in the validation dataset. In conclusion, we identified a blood-based 6-protein panel as a diagnostic tool in lung cancer. Used as a routine test for high- and average-risk individuals, it may complement currently adopted techniques in lung cancer screening.publishedVersio

Complexity of the microRNA repertoire revealed by next-generation sequencing

MicroRNAs (miRNAs) have been implicated to play key roles in normal physiological functions, and altered expression of specific miRNAs has been associated with a number of diseases. It is of great interest to understand their roles and a prerequisite for such study is the ability to comprehensively and accurately assess the levels of the entire repertoire of miRNAs in a given sample. It has been shown that some miRNAs frequently have sequence variations termed isomirs. To better understand the extent of miRNA sequence heterogeneity and its potential implications for miRNA function and measurement, we conducted a comprehensive survey of miRNA sequence variations from human and mouse samples using next generation sequencing platforms. Our results suggest that the process of generating this isomir spectrum might not be random and that heterogeneity at the ends of miRNA affects the consistency and accuracy of miRNA level measurement. In addition, we have constructed a database from our sequencing data that catalogs the entire repertoire of miRNA sequences (http://galas.systemsbiology.net/cgi-bin/isomir/find.pl). This enables users to determine the most abundant sequence and the degree of heterogeneity for each individual miRNA species. This information will be useful both to better understand the functions of isomirs and to improve probe or primer design for miRNA detection and measurement

MaRiMba: A Software Application for Spectral Library-Based MRM Transition List Assembly

Multiple reaction monitoring mass spectrometry (MRM-MS) is a targeted analysis method that has been increasingly viewed as an avenue to explore proteomes with unprecedented sensitivity and throughput. We have developed a software tool, called MaRiMba, to automate the creation of explicitly defined MRM transition lists required to program triple quadrupole mass spectrometers in such analyses. MaRiMba creates MRM transition lists from downloaded or custom-built spectral libraries, restricts output to specified proteins or peptides, and filters based on precursor peptide and product. ion properties. MaRiMba can also create MRM lists containing corresponding transitions for isotopically heavy peptides, for which the precursor and product ions are adjusted according to user specifications. This open-source application is operated through a graphical user interface incorporated into the Trans-Proteomic Pipeline, and it outputs the final MRM list to a text file for upload to MS instruments. To illustrate the use of MaRiMba, we used the tool to design and execute an MRM-MS experiment in which we targeted the proteins of a well-defined and previously published standard mixture

Monitoring of Circulating Tumor DNA Improves Early Relapse Detection After Axicabtagene Ciloleucel Infusion in Large B-Cell Lymphoma: Results of a Prospective Multi-Institutional Trial

Although the majority of patients with relapsed or refractory large B-cell lymphoma respond to axicabtagene ciloleucel (axi-cel), only a minority of patients have durable remissions. This prospective multicenter study explored the prognostic value of circulating tumor DNA (ctDNA) before and after standard-of-care axi-cel for predicting patient outcomes. Lymphoma-specific variable, diversity, and joining gene segments (VDJ) clonotype ctDNA sequences were frequently monitored via next-generation sequencing from the time of starting lymphodepleting chemotherapy until progression or 1 year after axi-cel infusion. We assessed the prognostic value of ctDNA to predict outcomes and axi-cel-related toxicity. A tumor clonotype was successfully detected in 69 of 72 (96%) enrolled patients. Higher pretreatment ctDNA concentrations were associated with progression after axi-cel infusion and developing cytokine release syndrome and/or immune effector cell-associated neurotoxicity syndrome. Twenty-three of 33 (70%) durably responding patients versus 4 of 31 (13%) progressing patients demonstrated nondetectable ctDNA 1 week after axi-cel infusion ( < .0001). At day 28, patients with detectable ctDNA compared with those with undetectable ctDNA had a median progression-free survival and OS of 3 months versus not reached ( < .0001) and 19 months versus not reached ( = .0080), respectively. In patients with a radiographic partial response or stable disease on day 28, 1 of 10 patients with concurrently undetectable ctDNA relapsed; by contrast, 15 of 17 patients with concurrently detectable ctDNA relapsed ( = .0001). ctDNA was detected at or before radiographic relapse in 29 of 30 (94%) patients. All durably responding patients had undetectable ctDNA at or before 3 months after axi-cel infusion. Noninvasive ctDNA assessments can risk stratify and predict outcomes of patients undergoing axi-cel for the treatment of large B-cell lymphoma. These results provide a rationale for designing ctDNA-based risk-adaptive chimeric antigen receptor T-cell clinical trials