Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical biology”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Creative-Dynamics Approach To Neural Intelligence

Paper discusses approach to mathematical modeling of artificial neural networks exhibiting complicated behaviors reminiscent of creativity and intelligence of biological neural networks. Neural network treated as non-Lipschitzian dynamical system - as described in "Non-Lipschitzian Dynamics For Modeling Neural Networks" (NPO-17814). System serves as tool for modeling of temporal-pattern memories and recognition of complicated spatial patterns.

Zak, Michail A.↗
NASA NTRS2013

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗
DOE OSTI2026

Towards modeling phage therapy

Patients infected with life-threatening multi-drug resistant (MDR) bacteria have been treated with cocktails of bacteriophages. This is a complicated form of personalized medicine as the phages given to a patient have to be selected beforehand on the basis of their lytic capacity of the infecting bacteria. Because bacteria rapidly become resistant, the evolution of resistance to a diverse cocktail of phages is a complicated dynamical process, during which competing bacterial strains replace one another by accumulating several resistance mechanisms, each of which may involve a fitness cost. As a consequence, it is typically not known why a particular phage therapy succeeded or failed, and how one can optimize the composition of the cocktails to maximize the rate of success. To improve upon this, we extend an existing in vivo -calibrated mouse model into a novel mathematical model for the human situation, and include multiple phages infecting multiple bacterial strains, differing in their resistance to each of the phages. We adjust several parameter estimates of the bacterial model to the human situation, and use the model to describe a successful case of phage therapy involving several cocktails, each containing several phages. In the model, treatment success crucially depended on pretreatment resistance levels, and on the diversity and the timing of the cocktails. Once an appropriate cocktail is found, it is less important to further optimize the infection rates of the phages. Resistant bacterial strains expand rapidly when sensitive strains decline, and the higher the infectivity of the phages, the faster resistant strains expand. Because resistance evolves rapidly, it is best to provide a diverse set of phages right from the start of therapy, i.e., to hit hard and early, and create a high genetic barrier to bacterial resistance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dynamics and growth rate implications of ribosomes and mRNAs interaction in E. coli

Understanding how cells grow and adapt under various nutrient conditions is pivotal in the study of biological stoichiometry. Recent studies provide empirical evidence that cells use multiple strategies to maintain an optimal protein production rate under different nutrient conditions. Mathematical models can provide a solid theoretical foundation that can explain experimental observations and generate testable hypotheses to further our understanding of the growth process. In this study, we generalize a modeling framework that centers on the translation process and study its asymptotic behaviors to validate algebraic manipulations involving the steady states. Using experimental results on the growth of E. coli under C-, N-, and P-limited environments, we simulate the expected quantitative measurements to show the feasibility of using the model to explain empirical evidence. Our results support the findings that cells employ multiple strategies to maintain a similar protein production rate across different nutrient limitations. Moreover, we find that the previous study underestimates the significance of certain biological rates, such as the binding rate of ribosomes to mRNA and the transition rate between different ribosomal stages. Furthermore, our simulation shows that the strategies used by cells under C- and P-limitations result in a faster overall growth dynamics than under N-limitation. In conclusion, the general modeling framework provides a valuable platform to study cell growth under different nutrient supply conditions, which also allows straightforward extensions to the coupling of transcription, translation, and energetics to deepen our understanding of the growth process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗
DOE OSTI2021

Cric searchable image database as a public platform for conventional pap smear cytology data

Amidst the current health crisis and social distancing, telemedicine has become an important part of mainstream of healthcare, and building and deploying computational tools to support screening more efficiently is an increasing medical priority. The early identification of cervical cancer precursor lesions by Pap smear test can identify candidates for subsequent treatment. However, one of the main challenges is the accuracy of the conventional method, often subject to high rates of false negative. While machine learning has been highlighted to reduce the limitations of the test, the absence of high-quality curated datasets has prevented strategies development to improve cervical cancer screening. The Center for Recognition and Inspection of Cells (CRIC) platform enables the creation of CRIC Cervix collection, currently with 400 images (1,376 × 1,020 pixels) curated from conventional Pap smears, with manual classification of 11,534 cells. This collection has the potential to advance current efforts in training and testing machine learning algorithms for the automation of tasks as part of the cytopathological analysis in the routine work of laboratories.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Molecular basis of crosstalk in nuclear receptors: heterodimerization between PXR and CAR and the implication in gene regulation

The 48 human nuclear receptors (NRs) form a superfamily of transcription factors that regulate major physiological and pathological processes. Emerging evidence suggests that NR crosstalk can fundamentally change our understanding of NR biology, but detailed molecular mechanisms of crosstalk are lacking. Here, we report the molecular basis of crosstalk between the pregnane X receptor (PXR) and constitutive androstane receptor (CAR), where they form a novel heterodimer, resulting in their mutual inhibition. PXR and CAR regulate drug metabolism and energy metabolism. Although they have been broadly perceived as functionally redundant, a growing number of reports suggests a mutual inhibitory relation, but their precise mode of coordinated action remains unknown. Using methods including RNA sequencing, small-angle X-ray scattering and crosslinking mass spectrometry we demonstrate that the mutual inhibition altered gene expression globally and is attributed to the novel PXR–CAR heterodimerization via the same interface used by each receptor to heterodimerize with its functional partner, retinoid X receptor (RXR). These findings establish an unexpected functional relation between PXR, CAR and RXR, change the perceived functional relation between PXR and CAR, open new perspectives on elucidating their role and designing approaches to regulate them, and highlight the importance to comprehensively investigate nuclear receptor crosstalk.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Arabidopsis cytochrome b 5 proteins support fatty acid ω-3 but not ω-6 desaturation

Fatty acids are primary components of lipids, which serve as major energy sources in cells and play essential roles in membrane structure, signaling, and metabolic regulation (Shanklin and Cahoon 1998). The degree of fatty acid unsaturation critically influences lipid physicochemical properties, thereby affecting membrane fluidity and biological function (Nguyen et al. 2019). In Arabidopsis thaliana, fatty acid desaturation occurs via 2 parallel pathways: the “prokaryotic pathway” in plastids, involving glycosylglycerides, such as monogalactosyldiacylglycerol (MGDG) and digalactosyldiacylglycerol (DGDG), and phospholipid phosphatidylglycerol (PG); and the “eukaryotic pathway” in the endoplasmic reticulum (ER), involving phosphatidylcholine (PC) (Lou et al. 2014) (Supplementary Figure S1). Seven fatty acid desaturases (FADs) in Arabidopsis differentially desaturate each glycerolipid class in the plastid and ER (Nguyen et al. 2019). FAD2, an ER-resident ω-6 fatty acid desaturase, catalyzes the conversion of oleic acid (18:1) to linoleic acid (18:2), which can be further desaturated to α-linolenic acid (18:3) by FAD3, an ER-resident ω-3 fatty acid desaturase. In plastids, FAD6 catalyzes the desaturation of 18:1/16:1 to produce 18:2/16:2, while FAD7 and FAD8 redundantly convert 18:2/16:2 to 18:3/16:3 (Li-Beisson et al. 2013; Nguyen et al. 2019). Additionally, fatty acids synthesized in the ER can also be reimported into plastids to their site of de novo synthesis (Xu et al. 2010). All FADs require reducing power, in the form of 2 electrons, for catalysis, but the sources of the electrons vary between their subcellular localizations. In the ER, FAD2 and FAD3 receive electrons from a cytochrome b 5 (CB5)-based electron transfer chain comprising cytochrome b 5 reductase (CBR) and CB5. In contrast, ferredoxin serves as the electron donor for plastid-localized FAD6, FAD7, and FAD8 (Ohlrogge and Browse 1995; Andreu et al. 2007). While the relative contributions of the 2 pathways to total cellular desaturation products vary across tissues and species, most polyunsaturated FA biosynthesis in seeds occurs via ER-resident FAD2 and FAD3 (Miquel and Browse 1992; Ohlrogge and Browse 1995).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluation of Potential Impacts of Free Chlorine during Washing of Fresh-Cut Leafy Greens on Escherichia coli O157:H7 Cross-Contamination and Risk of Illness

Addition of chlorine-based antimicrobial substances to fresh-cut leafy green wash water is done to minimize microbial cross-contamination during processing. We developed the FDA Leafy Green Risk Assessment Model (FDA-LGRAM) to quantify the impact of free chlorine concentration in wash water during fresh-cut lettuce processing on the extent of water-mediated cross-contamination between shredded lettuce and the associated risk of illness due to exposure to Escherichia coli O157:H7. At different contamination prevalence and levels of E. coli O157:H7 on incoming lettuce heads, the model compared the predicted prevalence of contaminated fresh-cut lettuce packages and the risk of illness per serving between: (1) a scenario where fresh-cut lettuce was packaged without washing; and (2) scenarios involving washing fresh-cut lettuce with different levels of free chlorine (0 ppm, 5 ppm, 10 ppm, 15 ppm, and 20 ppm) prior to packaging. Our results indicate that the free chlorine level in wash water has a substantial impact on the predicted prevalence of contaminated fresh-cut lettuce packages and the risk of illness associated with E. coli O157:H7 in fresh-cut lettuce. Results showed that the required level of free chlorine that can minimize water-mediated cross-contamination and reduce the corresponding risk of illness depended on contamination prevalence and levels of E. coli O157:H7 on incoming lettuce heads. Our model also indicated that the pathogen inactivation rate in wash water via free chlorine was a key model parameter that had a significant impact on the extent of cross-contamination during washing and the predicted associated risk of illness.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Updated distribution maps of predominant Culex mosquitoes across the Americas

Abstract Background Estimates of the geographical distribution of Culex mosquitoes in the Americas have been limited to state and provincial levels in the United States and Canada and based on data from the 1980s. Since these estimates were made, there have been many more documented observations of mosquitoes and new methods have been developed for species distribution modeling. Moreover, mosquito distributions are affected by environmental conditions, which have changed since the 1980s. This calls for updated estimates of these distributions to understand the risk of emerging and re-emerging mosquito-borne diseases. Methods We used contemporary mosquito data, environmental drivers, and a machine learning ecological niche model to create updated estimates of the geographical range of seven predominant Culex species across North America and South America: Culex erraticus , Culex nigripalpus , Culex pipiens , Culex quinquefasciatus , Culex restuans , Culex salinarius , and Culex tarsalis . Results We found that Culex mosquito species differ in their geographical range. Each Culex species is sensitive to both natural and human-influenced environmental factors, especially climate and land cover type. Some prefer urban environments instead of rural ones, and some are limited to tropical or humid areas. Many are found throughout the Central Plains of the USA. Conclusions Our updated contemporary Culex distribution maps may be used to assess mosquito-borne disease risk. It is critical to understand the current geographical distributions of these important disease vectors and the key environmental predictors structuring their distributions not only to assess current risk, but also to understand how they will respond to climate change. Since the environmental predictors structuring the geographical distribution of mosquito species varied, we hypothesize that each species may have a different response to climate change. Graphical abstract

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Transmission risk of Oropouche fever across the Americas

Abstract Background Vector-borne diseases (VBDs) are important contributors to the global burden of infectious diseases due to their epidemic potential, which can result in significant population and economic impacts. Oropouche fever, caused by Oropouche virus (OROV), is an understudied zoonotic VBD febrile illness reported in Central and South America. The epidemic potential and areas of likely OROV spread remain unexplored, limiting capacities to improve epidemiological surveillance. Methods To better understand the capacity for spread of OROV, we developed spatial epidemiology models using human outbreaks as OROV transmission-locality data, coupled with high-resolution satellite-derived vegetation phenology. Data were integrated using hypervolume modeling to infer likely areas of OROV transmission and emergence across the Americas. Results Models based on one-support vector machine hypervolumes consistently predicted risk areas for OROV transmission across the tropics of Latin America despite the inclusion of different parameters such as different study areas and environmental predictors. Models estimate that up to 5 million people are at risk of exposure to OROV. Nevertheless, the limited epidemiological data available generates uncertainty in projections. For example, some outbreaks have occurred under climatic conditions outside those where most transmission events occur. The distribution models also revealed that landscape variation, expressed as vegetation loss, is linked to OROV outbreaks. Conclusions Hotspots of OROV transmission risk were detected along the tropics of South America. Vegetation loss might be a driver of Oropouche fever emergence. Modeling based on hypervolumes in spatial epidemiology might be considered an exploratory tool for analyzing data-limited emerging infectious diseases for which little understanding exists on their sylvatic cycles. OROV transmission risk maps can be used to improve surveillance, investigate OROV ecology and epidemiology, and inform early detection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Structural models and functional annotations for the Sphagnum divinum proteome

This dataset contains the structural models for the primary transcripts of the Sphagnum divinum proteome. Additionally, for a subset of these proteins, sequence and structural alignment results are provided. This dataset represents the most thorough structural study of a Sphagnum species, also known as peat mosses, by providing three-dimensional atomic resolution structures of the majority of the encoded proteins as well as structural alignment results used in the application of annotating the proteome. References (DOI) AlphaFold v2 Monomer: https://doi.org/10.1038/s41586-021-03819-2. References (DOI) US-align2: https://doi.org/10.1038/s41592-022-01585-1

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural Models and Sequence Alignment Results of the Rhodospirillum rubrum Proteome

This dataset contains the structural models for the primary transcripts of the Rhodospirillum rubrum proteome as well as sequence alignment results for a subset of the encoded proteins. For each protein, the five models inferred from AlphaFold 2 are provided. The largest pTM-scoring model for each protein was energy minimized; this minimized structure as well as its AlphaFold pickle output file are also provided. This set of structures represent an alternate source of models for the R. rubrum proteome to those available in the AlphaFold Protein Structure Database. For proteins that have been annotated as hypothetical, sequence alignment results from the HHblits and SAdLSA alignment methods are provided. These methods are often more capable to resolve sequence homology than other methods. Therefore, the results from both HHblits and SAdLSA are provided to identify possible homologs for these challenging proteins. Numerous sequence databases are utilized for these alignments. References AlphaFold v2 Multimer: https://doi.org/10.1101/2021.10.04.463034. References HHBlits: https://doi.org/10.1186/s12859-019-3019-7. References SAdLSA: https://doi.org/10.3389/fbinf.2021.689960.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural Models and Sequence Alignment Results of the Desulfovibrio vulgaris Proteome

This dataset contains the structural models for the primary transcripts of the Desulfovibrio vulgaris proteome as well as sequence alignment results for a subset of the encoded proteins. For each protein, the five models inferred from AlphaFold 2 are provided. The largest pTM-scoring model for each protein was energy minimized; this minimized structure as well as its AlphaFold pickle output file are also provided. This set of structures represent an alternate source of models for the D. vulgaris proteome to those available in the AlphaFold Protein Structure Database (AFDB). This is a bit more complicated since the proteins reporting in the AFDB originate from an outdated form of the D. vulgaris sequence. The different versions of the D. vulgaris gene annotation are collected in the Chronology subdirectory; further consideration of these changes on the structural space of the proteome are currently underway. For proteins that have been annotated as hypothetical, sequence alignment results from the HHblits and SAdLSA alignment methods are provided. These methods are often more capable to resolve sequence homology than other methods. Therefore, the results from both HHblits and SAdLSA are provided to identify possible homologs for these challenging proteins. Numerous sequence databases are utilized for these alignments. References AlphaFold v2 Multimer: https://doi.org/10.1101/2021.10.04.463034. References HHblits: hhtps://doi.org/10.1186/s12859-019-3019-7. References SAdLSA: hhtps://doi.org/10.3389/fbinf.2021.689960.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural Models of the Rhodopseudomonas palustris Proteome

This dataset contains the structural models for the primary transcripts of the Rhodopseudomonas palustris proteome. For each protein, the five models inferred from AlphaFold 2 are provided. The largest pTM-scoring model for each protein was energy minimized; this minimized structure as well as its AlphaFold pickle output file are also provided. This set of structures represent an alternate source of models for the R. palustris proteome to those available in the AlphaFold Protein Structure Database.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1 , a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

2019 LDRD Annual Report

This summary report provides an overview of all LDRD projects at Argonne that concluded in Fiscal Year 2019. Many projects are funded for multiple years, the initial fiscal year for each project is indicated by the first four digits of the LDRD project number. Focus areas were: Advanced Computing; Biological and Environmental Science Capability Development; Hard X-ray Sciences; Manufacturing Science and Engineering; Materials and Chemistry; National and Global Security; Photon Sciences; Securing Energy and Critical Resources; and, The Universe as Our Laboratory (ULab).

36 MATERIALS SCIENCE↗
DOE OSTI2025

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗