Search NASA⌕ Search

SEARCH · Search NASA

Results for “cell biology modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Design patterns of biological cells

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Multi-scale models of whole cells: progress and challenges

Whole-cell modeling is “the ultimate goal” of computational systems biology and “a grand challenge for 21st century”. These complex, highly detailed models account for the activity of every molecule in a cell and serve as comprehensive knowledgebases for the modeled system. Their scope and utility far surpass those of other systems models. In fact, whole-cell models (WCMs) are an amalgam of several types of “system” models. The models are simulated using a hybrid modeling method where the appropriate mathematical methods for each biological process are used to simulate their behavior. Given the complexity of the models, the process of developing and curating these models is labor-intensive and to date only a handful of these models have been developed. While whole-cell models provide valuable and novel biological insights, and to date have identified some novel biological phenomena, their most important contribution has been to highlight the discrepancy between available data and observations that are used for the parametrization and validation of complex biological models. Another realization has been that current whole-cell modeling simulators are slow and to run models that mimic more complex (e.g., multi-cellular) biosystems, those need to be executed in an accelerated fashion on high-performance computing platforms. In this manuscript, we review the progress of whole-cell modeling to date and discuss some of the ways that they can be improved.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Altered microRNA expression in animal models of Huntington’s disease and potential therapeutic strategies

A review of recent animal models of Huntington’s disease showed many microRNAs had altered expression levels in the striatum and cerebral cortex, and which were mostly downregulated. Among the altered microRNAs were miR-9/9*, miR-29b, miR-124a, miR-132, miR-128, miR-139, miR-122, miR-138, miR-23b, miR-135b, miR-181 (all downregulated) and miR-448 (upregulated), and similar changes had been previously found in Huntington’s disease patients. In the animal cell studies, the altered microRNAs included miR-9, miR-9*, miR-135b, miR-222 (all downregulated) and miR-214 (upregulated). In the animal models, overexpression of miR-155 and miR-196a caused a decrease in mutant huntingtin mRNA and protein level, lowered the mutant huntingtin aggregates in striatum and cortex, and improved performance in behavioral tests. Improved performance in behavioral tests also occurred with overexpression of miR-132 and miR-124. In the animal cell models, overexpression of miR-22 increased the viability of rat primary cortical and striatal neurons infected with mutant huntingtin and decreased huntingtin -enriched foci of ≥ 2 µm. Also, overexpression of miR-22 enhanced the survival of rat primary striatal neurons treated with 3-nitropropionic acid. Exogenous expression of miR-214, miR-146a, miR-150, and miR-125b decreased endogenous expression of huntingtin mRNA and protein in Hdh Q111 /Hdh Q111 cells. Further studies with animal models of Huntington’s disease are warranted to validate these findings and identify specific microRNAs whose overexpression inhibits the production of mutant huntingtin protein and other harmful processes and may provide a more effective means of treating Huntington’s disease in patients and slowing its progression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Spatial modeling algorithms for reactions and transport in biological cells

Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Here we present Spatial Modeling Algorithms for Reactions and Transport (SMART), a software package that takes in high-level user specifications about cell signaling networks and then assembles and solves the associated mathematical systems. SMART uses state-of-the-art finite element analysis, via the FEniCS Project software, to efficiently and accurately resolve cell signaling events over discretized cellular and subcellular geometries. We demonstrate its application to several different biological systems, including yes-associated protein (YAP)/PDZ-binding motif (TAZ) mechanotransduction, calcium signaling in neurons and cardiomyocytes, and ATP generation in mitochondria. Throughout, we utilize experimentally derived realistic cellular geometries represented by well-conditioned tetrahedral meshes. These scenarios demonstrate the applicability, flexibility, accuracy and efficiency of SMART across a range of temporal and spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Online Machine Learning for Accelerating Molecular Dynamics Modeling of Cells

We developed a biomechanics-informed online learning framework to learn the dynamics with ground truth generated with multiscale modeling simulation. It was built on Summit-like supercomputers, which were also used to benchmark and validate our framework on one physiologically significant modeling of deformable biological cells. We generalized the century-old equation of Jeffery orbits to a new equation of motion with additional parameters to account for the flow conditions and the cell deformability. Using simulation data at particle-based resolutions for flowing cells and the learned parameters from our framework, we validated the new equation by the motions, mostly rotations, of a human platelet in shear blood flow at various shear stresses and platelet deformability. Our online framework, which surrogates redundant computations in the conventional multiscale modeling by solutions of our learned equation, accelerates the conventional modeling by three orders of magnitude without visible loss of accuracy.

multiscale modeling↗
DOE OSTI2024

Formation and retrieval of cell assemblies in a biologically realistic spiking neural network model of area CA3 in the mouse hippocampus

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

hetero-association↗
DOE OSTI2020

Art and Science of the Cellular Mesoscale

Experimental information from microscopy, structural biology, and bioinformatics may be integrated to build structural models of entire cells with molecular detail. This integrative modeling is challenging in several ways: the intrinsic complexity of biology results in models with many closely packed and heterogeneous components; the wealth of available experimental data is scattered among multiple resources and must be gathered, reconciled, and curated; and computational infrastructure is only now gaining the capability of modeling and visualizing systems of this complexity. We present recent efforts to address these challenges, both with artistic approaches to depicting the cellular mesoscale, and development and application of methods to build quantitative models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Modular cell-free expression plasmids to accelerate biological design in cells

Industrial biotechnology aims to produce high-value products from renewable resources. This can be challenging because model microorganisms—organisms that are easy to use like Escherichia coli—often lack the machinery required to utilize desired feedstocks like lignocellulosic biomass or syngas. Non-model organisms, such as Clostridium, are industrially proven and have desirable metabolic features but have several hurdles to mainstream use. Namely, these species grow more slowly than conventional laboratory microbes, and genetic tools for engineering them are far less prevalent. To address these hurdles for accelerating cellular design, cell-free synthetic biology has matured as an approach for characterizing non-model organisms and rapidly testing metabolic pathways in vitro. Unfortunately, cell-free systems can require specialized DNA architectures with minimal regulation that are not compatible with cellular expression. In this work, we develop a modular vector system that allows for T7 expression of desired enzymes for cell-free expression and direct Golden Gate assembly into Clostridium expression vectors. Utilizing the Joint Genome Institute’s DNA Synthesis Community Science Program, we designed and synthesized these plasmids and genes required for our projects allowing us to shuttle DNA easily between our in vitro and in vivo experiments. We next validated that these vectors were sufficient for cell-free expression of functional enzymes, performing on par with the previous state-of-the-art. Lastly, we demonstrated automated six-part DNA assemblies for Clostridium autoethanogenum expression with efficiencies ranging from 68% to 90%. We anticipate this system of plasmids will enable a framework for facile testing of biosynthetic pathways in vitro and in vivo by shortening development cycles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Study of impacts of two types of cellular aging on the yeast bud morphogenesis

Understanding the mechanisms of the cellular aging processes is crucial for attempting to extend organismal lifespan and for studying age-related degenerative diseases. Yeast cells divide through budding, providing a classical biological model for studying cellular aging. With their powerful genetics, relatively short cell cycle, and well-established signaling pathways also found in animals, yeast cells offer valuable insights into the aging process. Recent experiments suggested the existence of two aging modes in yeast characterized by nucleolar and mitochondrial declines, respectively. By analyzing experimental data, this study shows that cells evolving into those two aging modes behave differently when they are young. While buds grow linearly in both modes, cells that consistently generate spherical buds throughout their lifespan demonstrate greater efficacy in controlling bud size and growth rate at young ages. A three-dimensional multiscale chemical-mechanical model was developed and used to suggest and test hypothesized impacts of aging on bud morphogenesis. Experimentally calibrated model simulations showed that during the early stage of budding, tubular bud shape in one aging mode could be generated by locally inserting new materials at the bud tip, a process guided by the polarized Cdc42 signal. Furthermore, the aspect ratio of the tubular bud could be stabilized during the late stage as observed in experiments in this work. The model simulation results suggest that the localization of new cell surface material insertion, regulated by chemical signal polarization, could be weakened due to cellular aging in yeast and other cell types, leading to the change and stabilization of the bud aspect ratio.

Tsai, Kevin↗
DOE OSTI2022

MicroRNA expression in animal models of amyotrophic lateral sclerosis and potential therapeutic approaches

A review of recent animal models of amyotrophic lateral sclerosis showed a large number of miRNAs had altered levels of expression in the brain and spinal cord, motor neurons of spinal cord and brainstem, and hypoglossal, facial, and red motor nuclei and were mostly upregulated. Among the miRNAs found to be upregulated in two of the studies were miR-21, miR-155, miR-125b, miR-146a, miR-124, miR-9, and miR-19b, while those downregulated in two of the studies included miR-146a, miR-29, miR-9, and miR-125b. A change of direction in miRNA expression occurred in some tissues when compared (e.g., miR-29b-3p in cerebellum and spinal cord of wobbler mice at 40 days), or at different disease stages (e.g., miR-200a in spinal cord of SOD1(G93A) mice at 95 days vs. 108 and 112 days). In the animal models, suppression of miR-129-5p resulted in increased lifespan, improved muscle strength, reduced neuromuscular junction degeneration, and tended to improve motor neuron survival in the SOD1(G93A) mouse model. Suppression of miR-155 was also associated with increased lifespan, while lowering of miR-29a tended to improve lifespan in males and increase muscle strength in SOD1(G93A) mice. Overexpression of members of miR-17~92 cluster improved motor neuron survival in SOD1(G93A) mice. Treatment with an artificial miRNA designed to target hSOD1 increased lifespan and improved muscle strength in SOD1(G93A) animals. Further studies with animal models of amyotrophic lateral sclerosis are warranted to validate these findings and identify specific miRNAs whose suppression or directed against hSOD1 results in increased lifespan, improved muscle strength, reduced neuromuscular junction degeneration, and improved motor neuron survival in SOD1(G93A) animals.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Solving the puzzle of Fe homeostasis by integrating molecular, mathematical, and societal models

To ensure optimal utilization and bioavailability, iron uptake, transport, subcellular localization, and assimilation are tightly regulated in plants. In this work, we examine recent advances in our understanding of cellular responses to Fe deficiency. We then use intracellular mechanisms of Fe homeostasis to discuss how formalizing cell biology knowledge via a mathematical model can advance discovery even when quantitative data is limited. Using simulation-based inference to identify plausible systems mechanisms that conform to known emergent phenotypes can yield novel, testable hypotheses to guide targeted experiments. However, this approach relies on the accurate encoding of domain-expert knowledge in exploratory mathematical models. We argue that this would be facilitated by fostering more “systems thinking” life scientists and that diversifying your research team may be a practical path to achieve that goal.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Investigating key cell types and molecules dynamics in PyMT mice model of breast cancer through a mathematical model

The most common kind of cancer among women is breast cancer. Understanding the tumor microenvironment and the interactions between individual cells and cytokines assists us in arriving at more effective treatments. Here, we develop a data-driven mathematical model to investigate the dynamics of key cell types and cytokines involved in breast cancer development. We use time-course gene expression profiles of a mouse model to estimate the relative abundance of cells and cytokines. We then employ a least-squares optimization method to evaluate the model’s parameters based on the mice data. The resulting dynamics of the cells and cytokines obtained from the optimal set of parameters exhibit a decent agreement between the data and predictions. We perform a sensitivity analysis to identify the crucial parameters of the model and then perform a local bifurcation on them. The results reveal a strong connection between adipocytes, IL6, and the cancer population, suggesting them as potential targets for therapies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A dynamic kinetic model captures cell-free metabolism for improved butanol production

Cell-free systems are useful tools for prototyping metabolic pathways and optimizing the production of various bioproducts. Mechanistically-based kinetic models are uniquely suited to analyze dynamic experimental data collected from cell-free systems and provide vital qualitative insight. However, to date, dynamic kinetic models have not been applied with rigorous biological constraints or trained on adequate experimental data to the degree that they would give high confidence in predictions and broadly demonstrate the potential for widespread use of such kinetic models. In this work, we construct a large-scale dynamic model of cell-free metabolism with the goal of understanding and optimizing butanol production in a cell-free system. Using a combination of parameterization methods, the resultant model captures experimental metabolite measurements across two experimental conditions for nine metabolites at timepoints between 0 and 24 hours. Here, we present analysis of the model predictions, provide recommendations for butanol optimization, and identify the aldehyde/alcohol dehydrogenase as the primary bottleneck in butanol production. Sensitivity analysis further reveals the extent to which various parameters are constrained, and our approach for probing valid parameter ranges can be applied to other modeling efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mechanistic insights explain the transforming potential of the T507K substitution in the protein-tyrosine phosphatase SHP2

The protein-tyrosine phosphatase SHP2 is an allosteric enzyme critical for cellular events downstream of growth factor receptors. Mutations in the SHP2 gene have been linked to many different types of human diseases, including developmental disorders, leukemia, and solid tumors. Unlike most SHP2-activating mutations, the T507K substitution in SHP2 is unique in that it exhibits oncogenic Ras-like transforming activity. However, the biochemical basis of how the SHP2/T507K variant elicits transformation remains unclear. By combining kinetic and biophysical methods, X-ray crystallography, and molecular modeling, as well as using cell biology approaches, in this work we uncovered that the T507K substitution alters both SHP2 substrate specificity and its allosteric regulatory mechanism. We found that although SHP2/T507K exists in the closed, autoinhibited conformation similar to the WT enzyme, the interactions between its N-SH2 and protein-tyrosine phosphatase domains are weakened such that SHP2/T507K possesses a higher affinity for the scaffolding protein Grb2-associated binding protein 1 (Gab1). We also discovered that the T507K substitution alters the structure of the SHP2 active site, resulting in a change in SHP2 substrate preference for Sprouty1, a known negative regulator of Ras signaling and a potential tumor suppressor. Our results suggest that SHP2/T507K's shift in substrate specificity coupled with its preferential association of SHP2/T507K with Gab1 enable the mutant SHP2 to more efficiently dephosphorylate Sprouty1 at pTyr-53. This dephosphorylation hyperactivates Ras signaling, which is likely responsible for SHP2/T507K's Ras-like transforming activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Generative diffusion model surrogates for mechanistic agent-based biological models

Mechanistic, multicellular, agent-based models are commonly used to investigate tissue, organ, and organism-scale biology at single-cell resolution. The Cellular-Potts Model (CPM) is a powerful and popular framework for developing and interrogating these models. CPMs become computationally expensive at large space- and time- scales making application and investigation of developed models difficult. Surrogate models may allow for the accelerated evaluation of CPMs of complex biological systems. However, the stochastic nature of these models means each set of parameters may give rise to different model configurations, complicating surrogate model development. In this work, we leverage denoising diffusion probabilistic models (DDPMs) to train a generative AI surrogate of a CPM used to investigate in vitro vasculogenesis. We describe the use of an image classifier to learn the characteristics that define unique areas of a 2-dimensional parameter space. We then apply this classifier to aid in surrogate model selection and verification. Our CPM model surrogate generates model configurations 20,000 timesteps ahead of a reference configuration and demonstrates approximately a 22x reduction in computational time as compared to native code execution. Our work represents a step towards the implementation of DDPMs to develop digital twins of stochastic biological systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Kullback–Leibler Divergence of an Open-Queuing Network of a Cell-Signal-Transduction Cascade

Queuing networks (QNs) are essential models in operations research, with applications in cloud computing and healthcare systems. However, few studies have analyzed the cell’s biological signal transduction using QN theory. This study entailed the modeling of signal transduction as an open Jackson’s QN (JQN) to theoretically determine cell signal transduction, under the assumption that the signal mediator queues in the cytoplasm, and the mediator is exchanged from one signaling molecule to another through interactions between the signaling molecules. Each signaling molecule was regarded as a network node in the JQN. The JQN Kullback–Leibler divergence (KLD) was defined using the ratio of the queuing time (λ) to the exchange time (μ), λ/μ. The mitogen-activated protein kinase (MAPK) signal-cascade model was applied, and the KLD rate per signal-transduction-period was shown to be conserved when the KLD was maximized. Our experimental study on MAPK cascade supported this conclusion. This result is similar to the entropy-rate conservation of chemical kinetics and entropy coding reported in our previous studies. Thus, JQN can be used as a novel framework to analyze signal transduction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗