NNVA: Neural Network Assisted Visual Analysis of Yeast Cell Polarization Simulation

Complex computational models are often designed to simulate real-world physical phenomena in many scientific disciplines. However, these simulation models tend to be computationally very expensive and involve a large number of simulation input parameters which need to be analyzed and properly calibrated before the models can be applied for real scientific studies. We propose a visual analysis system to facilitate interactive exploratory analysis of high-dimensional input parameter space for a complex yeast cell polarization simulation. The proposed system can assist the computational biologists, who designed the simulation model, to visually calibrate the input parameters by modifying the parameter values and immediately visualizing the predicted simulation outcome without having the need to run the original expensive simulation for every instance. Our proposed visual analysis system is driven by a trained neural network-based surrogate model as the backend analysis framework. Surrogate models are widely used in the field of simulation sciences to efficiently analyze computationally expensive simulation models. In this work, we demonstrate the advantage of using neural networks as surrogate models for visual analysis by incorporating some of the recent advances in the field of uncertainty quantification, interpretability and explainability of neural network-based models. We utilize the trained network to perform interactive parameter sensitivity analysis of the original simulation at multiple levels-of-detail as well as recommend optimal parameter configurations using the activation maximization framework of neural networks. We also facilitate detail analysis of the trained network to extract useful insights about the simulation model, learned by the network, during the training process.Comment: Published at IEEE Transactions on Visualization and Computer Graphic

Noise filtering tradeoffs in spatial gradient sensing and cell polarization response

<p>Abstract</p> <p>Background</p> <p>Cells sense chemical spatial gradients and respond by polarizing internal components. This process can be disrupted by gradient noise caused by fluctuations in chemical concentration.</p> <p>Results</p> <p>We investigated how external gradient noise affects spatial sensing and response focusing on noise-filtering and the resultant tradeoffs. First, using a coarse-grained mathematical model of gradient-sensing and cell polarity, we characterized three negative consequences of noise: Inhibition of the extent of polarization, degradation of directional accuracy, and production of a noisy output polarization. Next, we explored filtering strategies and discovered that a combination of positive feedback, multiple signaling stages, and time-averaging produced good results. There was an important tradeoff, however, because filtering resulted in slower polarization. Simulations demonstrated that a two-stage filter-amplifier resulted in a balanced outcome. Then, we analyzed the effect of noise on a mechanistic model of yeast cell polarization in response to gradients of mating pheromone. This analysis showed that yeast cells likely also combine the above three filtering mechanisms into a filter-amplifier structure to achieve impressive spatial-noise tolerance, but with the consequence of a slow response time. Further investigation of the amplifier architecture revealed two positive feedback loops, a fast inner and a slow outer, both of which contributed to noise-tolerant polarization. This model also made specific predictions about how orientation performance depended upon the ratio between the gradient slope (signal) and the noise variance. To test these predictions, we performed microfluidics experiments measuring the ability of yeast cells to orient to shallow gradients of mating pheromone. The results of these experiments agreed well with the modeling predictions, demonstrating that yeast cells can sense gradients shallower than 0.1% μm^-1, approximately a single receptor-ligand molecule difference between front and back, on par with motile eukaryotic cells.</p> <p>Conclusions</p> <p>Spatial noise impedes the extent, accuracy, and smoothness of cell polarization. A combined filtering strategy implemented by a filter-amplifier architecture with slow dynamics was effective. Modeling and experimental data suggest that yeast cells employ these elaborate mechanisms to filter gradient noise resulting in a slow but relatively accurate polarization response.</p

Computations with parametric equations

Abstract. We present a complete method of implicitization for general rational parametric equations. We also present a method to decide whether the parameters of a set of parametric equations are independent, and if not, reparameterize the parametric equations so that the new parametric equations have independent parameters. We give a method to compute the inversion maps of parametric equations with independent parameters, and as a consequence, we can decide whether the parametric equations are proper. A new method to find a proper reparameterization for a set of improper parametric equations of algebraic curves is presented.

Robust Spatial Sensing of Mating Pheromone Gradients by Yeast Cells

Projecting or moving up a chemical gradient is a universal behavior of living organisms. We tested the ability of S. cerevisiae a-cells to sense and respond to spatial gradients of the mating pheromone α-factor produced in a microfluidics chamber; the focus was on bar1Δ strains, which do not degrade the pheromone input. The yeast cells exhibited good accuracy with the mating projection typically pointing in the correct direction up the gradient (∼80% under certain conditions), excellent sensitivity to shallow gradients, and broad dynamic range so that gradient-sensing was relatively robust over a 1000-fold range of average α-factor concentrations. Optimal directional sensing occurred at lower concentrations (5 nM) close to the Kd of the receptor and with steeper gradient slopes. Pheromone supersensitive mutations (sst2Δ and ste2300Δ) that disrupt the down-regulation of heterotrimeric G-protein signaling caused defects in both sensing and response. Interestingly, yeast cells employed adaptive mechanisms to increase the robustness of the process including filamentous growth (i.e. directional distal budding) up the gradient at low pheromone concentrations, bending of the projection to be more aligned with the gradient, and forming a more accurate second projection when the first projection was in the wrong direction. Finally, the cells were able to amplify a shallow external gradient signal of α-factor to produce a dramatic polarization of signaling proteins at the front of the cell. Mathematical modeling revealed insights into the mechanism of this amplification and how the supersensitive mutants can disrupt accurate polarization. Together, these data help to specify and elucidate the abilities of yeast cells to sense and respond to spatial gradients of pheromone

Modeling yeast cell polarization induced by pheromone gradients

Yeast cells respond to spatial gradients of mating pheromones by polarizing and projecting up the gradient toward the source. It is thought that they employ a spatial sensing mechanism in which the cell compares the concentration of pheromone at different points on the cell surface and determines the maximum point, where the projection forms. Here we constructed the first spatial mathematical model of the yeast pheromone response that describes the dynamics of the heterotrimeric and Cdc42p G-protein cycles, which are linked in a cascade. Two key performance objectives of this system are (1) amplification-converting a shallow external gradient of ligand to a steep internal gradient of protein components and (2) tracking-following changes in gradient direction. We used simulations to investigate amplification mechanisms that allow tracking. We identified specific strategies for regulating the spatial dynamics of the protein components (i.e. their changing location in the cell) that would enable the cell to achieve both objectives

Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson-Crick base pair at the two drug-binding interfaces can result in a single-stranded 'chair-shaped' DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a 'curved-shaped' duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs' potential synergetic use

REG3A overexpression functions as a negative predictive and prognostic biomarker in rectal cancer patients receiving CCRT

Background. Concurrent chemoradiotherapy (CCRT) is suggested before resection surgery in the control of rectal cancer. Unfortunately, treatment outcomes are widely variable and highly patientspecific. Notably, rectal cancer patients with distant metastasis generally have a much lower survival rate. Accordingly, a better understanding of the genetic background of patient cohorts can aid in predicting CCRT efficacy and clinical outcomes for rectal cancer before distant metastasis. Methods. A published transcriptome dataset (GSE35452) (n=46) was utilized to distinguish prospective genes concerning the response to CCRT. We recruited 172 rectal cancer patients, and the samples were collected during surgical resection after CCRT. Immunohistochemical (IHC) staining was performed to evaluate the expression level of regenerating family member 3 alpha (REG3A). Pearson's chi-squared test appraised the relevance of REG3A protein expression to clinicopathological parameters. The Kaplan-Meier method was utilized to generate survival curves, and the log-rank test was performed to compare the survival distributions between two given groups. Results. Employing a transcriptome dataset (GSE35452) and focusing on the inflammatory response (GO: 0006954), we recognized that REG3A is the most significantly upregulated gene among CCRT nonresponders (log2 ratio=1.2472, p=0.0079). Following IHC validation, high immunoexpression of REG3A was considerably linked to advanced post-CCRT tumor status (p<0.001), post-CCRT lymph node metastasis (p=0.042), vascular invasion (p=0.028), and low-grade tumor regression (p=0.009). In the multivariate analysis, high immunoexpression of REG3A was independently correlated with poor disease-specific survival (DSS) (p=0.004) and metastasis-free survival (MeFS) (p=0.045). The results of the bioinformatic analysis also supported the idea that REG3A overexpression is implicated in rectal carcinogenesis. Conclusion. In the current study, we demonstrated that REG3A overexpression is correlated with poor CCRT effectiveness and inferior patient survival in rectal cancer. The predictive and prognostic utility of REG3A expression may direct patient stratification and decisionmaking more accurately for those patients