Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biological system 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale Modeling Meets Machine Learning: What Can We Learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learn- ing to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

machine learning, multiscale modeling, physics-bas↗
DOE OSTI2023

Standards, dissemination, and best practices in systems biology

In this study, the reproducibility of scientific research is crucial to the success of the scientific method. Here, we review the current best practices when publishing mechanistic models in systems biology. We recommend, where possible, to use software engineering strategies such as testing, verification, validation, documentation, versioning, iterative development, and continuous integration. In addition, adhering to the Findable, Accessible, Interoperable, and Reusable modeling principles allows other scientists to collaborate and build off of each other’s work. Existing standards such as Systems Biology Markup Language, CellML, or Simulation Experiment Description Markup Language can greatly improve the likelihood that a published model is reproducible, especially if such models are deposited in well-established model repositories. Where models are published in executable programming languages, the source code and their data should be published as open-source in public code repositories together with any documentation and testing code. For complex models, we recommend container-based solutions where any software dependencies and the run-time context can be easily replicated.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Investigating the Unique Ability of Trichodesmium To Fix Carbon and Nitrogen Simultaneously Using MiMoSA

The open ocean is an extremely competitive environment, partially due to the dearth of nutrients. Trichodesmium erythraeum, a marine diazotrophic cyanobacterium, is a keystone species in the ocean due to its ability to fix nitrogen and leak 30 to 50% into the surrounding environment, providing a valuable source of a necessary macronutrient to other species. While there are other diazotrophic cyanobacteria that play an important role in the marine nitrogen cycle, Trichodesmium is unique in its ability to fix both carbon and nitrogen simultaneously during the day without the use of specialized cells called heterocysts to protect nitrogenase from oxygen. Here, we use the advanced modeling framework called multiscale multiobjective systems analysis (MiMoSA) to investigate how Trichodesmium erythraeum can reduce dimolecular nitrogen to ammonium in the presence of oxygen. Our simulations indicate that nitrogenase inhibition is best modeled as Michealis-Menten competitive inhibition and that cells along the filament maintain microaerobia using high flux through Mehler reactions in order to protect nitrogenase from oxygen. We also examined the effect of location on metabolic flux and found that cells at the end of filaments operate in distinctly different metabolic modes than internal cells despite both operating in a photoautotrophic mode. These results give us important insight into how this species is able to operate photosynthesis and nitrogen fixation simultaneously, giving it a distinct advantage over other diazotrophic cyanobacteria because they can harvest light directly to fuel the energy demand of nitrogen fixation.

09 BIOMASS FUELS↗
DOE OSTI2020

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo

Stochastic reaction network models are often used to explain and predict the dynamics of gene regulation in single cells.These models usually involve several parameters, such as the kinetic rates of chemical reactions, that are not directly measurable and must be inferred from experimental data. Bayesian inference provides a rigorous probabilistic frame-work for identifying these parameters by finding a posterior parameter distribution that captures their uncertainty.Traditional computational methods for solving inference problems such as Markov Chain Monte Carlo methods based on classical Metropolis-Hastings algorithm involve numerous serial evaluations of the likelihood function, which in turn requires expensive forward solutions of the chemical master equation (CME). We propose an alternate approach based on a multifidelity extension of the Sequential Tempered Markov Chain Monte Carlo (ST-MCMC) sampler. This algorithm is built upon Sequential Monte Carlo and solves the Bayesian inference problem by decomposing it into a sequence of efficiently solved subproblems that gradually increase both model fidelity and the influence of the observed data. We reformulate the finite state projection (FSP) algorithm, a well-known method for solving the CME, to produce a hierarchy of surrogate master equations to be used in this multifidelity scheme. To determine the appropriate fidelity, we introduce a novel information-theoretic criteria that seeks to extract the most information about the ultimate Bayesian posterior from each model in the hierarchy without inducing significant bias. This novel sampling scheme is tested with high performance computing resources using biologically relevant problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Advocating Feedback Control for Human-Earth System Applications

This paper proposes a feedback control perspective for Human-Earth Systems (HESs) which essentially are complex systems that capture the interactions between humans and nature. Recent attention in HES research has been directed towards devising strategies for climate change mitigation and adaptation, aimed at achieving environmental and societal objectives. However, existing approaches heavily rely on HES models, which inherently suffer from inaccuracies due to the complexity of the system. Moreover, overly detailed models often prove impractical for optimization tasks. We propose a framework inheriting from feedback control strategies the robustness against model errors, because inaccuracies are mitigated using measurements retrieved from the field. The framework comprises two nested control loops. The outer loop computes the optimal inputs to the HES, which are then implemented by actuators controlled in the inner loop. Potential fields of applications are also identified and a numerical example is provided.

biological system modeling↗
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 OSTI2020

Logic-based analysis of gene expression data predicts association between TNF, TGFB1 and EGF pathways in basal-like breast cancer

For breast cancer, clinically important subtypes are well characterized at the molecular level in terms of gene expression profiles. In addition, signaling pathways in breast cancer have been extensively studied as therapeutic targets due to their roles in tumor growth and metastasis. However, it is challenging to put signaling pathways and gene expression profiles together to characterize biological mechanisms of breast cancer subtypes since many signaling events result from post-translational modifications, rather than gene expression differences. We designed a logic-based computational framework to explain the differences in gene expression profiles among breast cancer subtypes using Pathway Logic and transcriptional network information. Pathway Logic is a rewriting-logic-based formal system for modeling biological pathways including post-translational modifications. Our method demonstrated its utility by constructing subtype-specific path from key receptors (TNFR, TGFBR1 and EGFR) to key transcription factor (TF) regulators (RELA, ATF2, SMAD3 and ELK1) and identifying potential association between pathways via TFs in basal-specific paths, which could provide a novel insight on aggressive breast cancer subtypes. Lastly, codes and results are available at http://epigenomics.snu.ac.kr/PL/.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Genomes to Structure and Function Workshop Report 2022

The goal of the U.S. Department of Energy (DOE) Biological and Environmental Research (BER) Program is to achieve a predictive understanding of complex biological, earth, and environmental systems with the aim of advancing the nation’s energy and infrastructure security. (https://www.energy.gov/science/ ber/biological-and-environmental-research). To pursue this goal, collaborations among experts in diverse research areas that lead to multidisciplinary projects are indispensable. The roles of DOE’s User Facilities, which offer unique and powerful resources for such research projects, are evolving, and expectations for the facilities are increasing. To respond to Users’ needs, the Joint Genome Institute (JGI) and Environmental Molecular Sciences Laboratory (EMSL) initiated the Facilities Integrating Collaborations for User Science (FICUS) program in 2014. This collaboration has grown into a popular and successful program, advancing more than 100 multidisciplinary projects to date. Similarly, the new interFacility collaborations among the JGI, EMSL, and User resources for BER structural biology and imaging at the Basic Energy Science (BES) Program’s synchrotron and neutron facilities are becoming essential for cutting-edge transdisciplinary science. To further explore the need for the BER research community to combine genomic, functional, and structural approaches to advance their research, an organizing committee was formed to develop and jointly host a 3-part workshop. The committee’s members represented seven DOE National Laboratory User Facilities (Appendix 1 lists the members). The “Genomes to Structure and Function” virtual workshop (see Appendices 2–5) was composed of three sessions. The first session, titled “Molecular Structures” (October 27– 28, 2021), highlighted diverse integrative experimental and computational approaches correlating structural data with sequencing and functional information, as well as predicting protein structures to model complex biological systems. The second session, “Intracellular Organization, and Material Synthesis and Decomposition” (December 15–16, 2021), covered imaging methods for observing, quantifying, and manipulating biosystems. The third session, “Imaging the Rhizosphere and Cellular Organization” (January 26–27, 2022) emphasized advanced and non-invasive imaging techniques applied to plant root-microbe-soil interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

In Situ Investigation of Upper Airway Occlusion in Sleep Disordered Breathing Using Ultrasonic Transducer Arrays

This work presents a novel application of ultrasound for the real-time, non-invasive investigation of occlusion of the upper airway during events of obstructive sleep apnea/hypopnea syndrome. It is hypothesized that ultrasonic pulses applied to the neck during apneic events produce spectral and temporal features that can detect apnea occurrence. Theoretical models of ultrasound propagation and an in vitro test were conducted to test this hypothesis in both transmission and reflection modes. Complete specifications and technical details of the system design and fabrication, which is mounted on each subject’s neck, are presented, including the methodology. Nine patients (seven male and two female, mean age of 42 years, with a range of 25 to 56 years, and body mass index 37.6 ± 6.6 kg/m2) were recruited for a full night study, which included simultaneous nocturnal polysomnography for the validation of the results. Nine temporal features and four spectral features were extracted from the envelope of the received pulse waveform. These were used to compute 26 metrics to quantify the changes in the ultrasonic waveforms between normal breathing and apneic events. The statistical analysis of the collected ultrasonic data showed that at least two or more of the proposed features could detect apneic events in all subjects. The findings establish the feasibility of the proposed method as a cost-effective and non-invasive OSAHS screening tool.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Transport and Assembly of Magnetic Surface Rotors**

This minireview focuses on recent advances with surface magnetic rotors, namely field-responsive spherical or anisotropic microparticles that translate close to, or are embedded in a confining surface. The application of external magnetic modulations allows these microscopic wheels to be remotely spun and steered while also tuning their interactions and inducing assembly from a collection of disordered, moving units. With optical microscopy one can observe and characterize the complex collective phenomena that emerge in dissipative colloidal systems driven far from equilibrium by external fields. From a technological point of view, magnetic surface rotors envisage implementation into microfluidic devices, and can be used for drug delivery, mixing as well as a model system for biological active matter.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Central transcriptional regulator controls photosynthetic growth and carbon storage in response to high light

Carbon capture and biochemical storage are some of the primary drivers of photosynthetic yield and productivity. To elucidate the mechanisms governing carbon allocation, we designed a photosynthetic light response test system for genetic and metabolic carbon assimilation tracking, using microalgae as simplified plant models. The systems biology mapping of high light-responsive photophysiology and carbon utilization dynamics between two variants of the same Picochlorum celeri species, TG1 and TG2 elucidated metabolic bottlenecks and transport rates of intermediates using instationary 13 C-fluxomics. Simultaneous global gene expression dynamics showed 73% of the annotated genes responding within one hour, elucidating a singular, diel-responsive transcription factor, closely related to the CCA1/LHY clock genes in plants, with significantly altered expression in TG2. Transgenic P. celeri TG1 cells expressing the TG2 CCA1/LHY gene, showed 15% increase in growth rates and 25% increase in storage carbohydrate content, supporting a coordinating regulatory function for a single transcription factor.

09 BIOMASS FUELS↗
DOE OSTI2022

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Budget-Constrained Sizing of Renewable and Energy Storage Systems for Farm-Scale Ammonia Production Within the Food-Energy-Water Nexus

In the transition toward sustainable agriculture, farms have emerged as eco-friendly pioneers, harnessing clean hybrid wind and solar systems to improve farm performance. A concern in this paradigm is the effective sizing of renewable energy systems to ensure optimal energy use within budget considerations. This research focuses on optimizing renewable energy sizing in small-scale ammonia production to meet specific farm demands and enhance local resilience, emphasizing the interplay between environmental and economic factors. These findings promise increased energy efficiency and sustainability in this innovative agricultural sector. Additionally, our approach considers small-scale ammonia plant needs and the dynamic relationships between ammonia, water, and farm demands. Simulations demonstrate substantial cost savings in farm electricity consumption. Specifically, scenarios with renewable energy integration in the farm can reduce at least 13% electricity cost compared to a grid-dependent system in the 15-year simulation.

ammonia↗
DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Biodiversity and Stoichiometric Plasticity Increase Pico‐Phytoplankton Contributions to Marine Net Primary Productivity and the Biological Pump

Earth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low-latitude oligotrophic systems harboring pico-sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon-climate system. We have added three pico-phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico-eukaryote populations that contribute significantly to productivity and export within the subtropical to mid-latitude transition zone, with the western subtropical regions of each basin supporting the most P-poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C-rich/P-poor stoichiometry, preferentially inhabited by pico-eukaryotes and diatoms. Collectively, pico-phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico-phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.

Environmental Sciences & Ecology↗