Decoding the coupled decision-making of the epithelial-mesenchymal transition and metabolic reprogramming in cancer.

Cancer metastasis relies on an orchestration of traits driven by different interacting functional modules, including metabolism and epithelial-mesenchymal transition (EMT). During metastasis, cancer cells can acquire a hybrid metabolic phenotype (W/O) by increasing oxidative phosphorylation without compromising glycolysis and they can acquire a hybrid epithelial/mesenchymal (E/M) phenotype by engaging EMT. Both the W/O and E/M states are associated with high metastatic potentials, and many regulatory links coupling metabolism and EMT have been identified. Here, we investigate the coupled decision-making networks of metabolism and EMT. Their crosstalk can exhibit synergistic or antagonistic effects on the acquisition and stability of different coupled metabolism-EMT states. Strikingly, the aggressive E/M-W/O state can be enabled and stabilized by the crosstalk irrespective of these hybrid states' availability in individual metabolism or EMT modules. Our work emphasizes the mutual activation between metabolism and EMT, providing an important step toward understanding the multifaceted nature of cancer metastasis.

Stochastic fluctuations promote ordered pattern formation of cells in the Notch-Delta signaling pathway.

The Notch-Delta signaling pathway mediates cell differentiation implicated in many regulatory processes including spatiotemporal patterning in tissues by promoting alternate cell fates between neighboring cells. At the multicellular level, this "lateral inhibition" principle leads to checkerboard patterns with alternation of Sender and Receiver cells. While it is well known that stochasticity modulates cell fate specification, little is known about how stochastic fluctuations at the cellular level propagate during multicell pattern formation. Here, we model stochastic fluctuations in the Notch-Delta pathway in the presence of two different noise types-shot and white-for a multicell system. Our results show that intermediate fluctuations reduce disorder and guide the multicell lattice toward checkerboard-like patterns. By further analyzing cell fate transition events, we demonstrate that intermediate noise amplitudes provide enough perturbation to facilitate "proofreading" of disordered patterns and cause cells to switch to the correct ordered state (Sender surrounded by Receivers, and vice versa). Conversely, high noise can override environmental signals coming from neighboring cells and lead to switching between ordered and disordered patterns. Therefore, in analogy with spin glass systems, intermediate noise levels allow the multicell Notch system to escape frustrated patterns and relax towards the lower energy checkerboard pattern while at large noise levels the system is unable to find this ordered base of attraction.

Contributions of higher-order proximal distribution functions to solvent structure around proteins.

The proximal distribution function (pDF) quantifies the probability of finding a solvent molecule in the vicinity of solutes. The approach constitutes a hierarchically organized theory for constructing approximate solvation structures around solutes. Given the assumption of universality of atom cluster-specific solvation, reconstruction of the solvent distribution around arbitrary molecules provides a computationally convenient route to solvation thermodynamics. Previously, such solvent reconstructions usually considered the contribution of the nearest-neighbor distribution only. We extend the pDF reconstruction algorithm to terms including next-nearest-neighbor contribution. As a test, small molecules (alanine and butane) are examined. The analysis is then extended to include the protein myoglobin in the P6 crystal unit cell. Molecular dynamics simulations are performed, and solvent density distributions around the solute molecules are compared with the results from different pDF reconstruction models. It is shown that the next-nearest-neighbor modification significantly improves the reconstruction of the solvent number density distribution in concave regions and between solute molecules. The probability densities are then used to calculate the solute-solvent non-bonded interaction energies including van der Waals and electrostatic, which are found to be in good agreement with the simulated values.

Towards decoding the coupled decision-making of metabolism and epithelial-to-mesenchymal transition in cancer.

Cancer cells have the plasticity to adjust their metabolic phenotypes for survival and metastasis. A developmental programme known as epithelial-to-mesenchymal transition (EMT) plays a critical role during metastasis, promoting the loss of polarity and cell-cell adhesion and the acquisition of motile, stem-cell characteristics. Cells undergoing EMT or the reverse mesenchymal-to-epithelial transition (MET) are often associated with metabolic changes, as the change in phenotype often correlates with a different balance of proliferation versus energy-intensive migration. Extensive crosstalk occurs between metabolism and EMT, but how this crosstalk leads to coordinated physiological changes is still uncertain. The elusive connection between metabolism and EMT compromises the efficacy of metabolic therapies targeting metastasis. In this review, we aim to clarify the causation between metabolism and EMT on the basis of experimental studies, and propose integrated theoretical-experimental efforts to better understand the coupled decision-making of metabolism and EMT.

Decoding the mechanisms underlying cell-fate decision-making during stem cell differentiation by random circuit perturbation.

Stem cells can precisely and robustly undergo cellular differentiation and lineage commitment, referred to as stemness. However, how the gene network underlying stemness regulation reliably specifies cell fates is not well understood. To address this question, we applied a recently developed computational method, random circuit perturbation (RACIPE), to a nine-component gene regulatory network (GRN) governing stemness, from which we identified robust gene states. Among them, four out of the five most probable gene states exhibit gene expression patterns observed in single mouse embryonic cells at 32-cell and 64-cell stages. These gene states can be robustly predicted by the stemness GRN but not by randomized versions of the stemness GRN. Strikingly, we found a hierarchical structure of the GRN with the Oct4/Cdx2 motif functioning as the first decision-making module followed by Gata6/Nanog. We propose that stem cell populations, instead of being viewed as all having a specific cellular state, can be regarded as a heterogeneous mixture including cells in various states. Upon perturbations by external signals, stem cells lose the capacity to access certain cellular states, thereby becoming differentiated. The new gene states and key parameters regulating transitions among gene states proposed by RACIPE can be used to guide experimental strategies to better understand differentiation and design reprogramming. The findings demonstrate that the functions of the stemness GRN is mainly determined by its well-evolved network topology rather than by detailed kinetic parameters.

Identifying trends in hydration behavior for modifications to the hydrophobicity of poly(n-isopropylacrylamide).

Poly(n-isopropylacrylamide), PNIPAM, is a thermo-responsive polymer that has been thoroughly studied for its many applications, such as drug delivery and actuators. Below the lower critical solution temperature (LCST), PNIPAM is well hydrated in the elongated conformation. The transition occuring at the LCST results in a less hydrated collapsed state above the LCST. This volume phase transition is dependent upon the hydration of the polymer and its hydrophobicity. Some research has been done on potential modifications of PNIPAM for applications, but until now there has not been a study of the hydration properties as a function of hydrophobicity. The work presented in this paper applies a Voronoi analysis of the hydration of PNIPAM, as well as PNIPAM with other alkyl substituents. We show from classical MD simulations that increasing hydrophobicity can increase the volume phase change, but there is a lower limit to this trend. Additionally, replica exchange molecular dynamics were conducted on PNIPAM showing a fluctuation between elongated and collapsed states near the LCST.