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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Forecasting Multi-Wave Epidemics Through Bayesian Inference

We present a simple, near-real-time Bayesian method to infer and forecast a multiwave outbreak, and demonstrate it on the COVID-19 pandemic. The approach uses timely epidemiological data that has been widely available for COVID-19. It provides short-term forecasts of the outbreak’s evolution, which can then be used for medical resource planning. The method postulates one- and multiwave infection models, which are convolved with the incubation-period distribution to yield competing disease models. The disease models’ parameters are estimated via Markov chain Monte Carlo sampling and information-theoretic criteria are used to select between them for use in forecasting. The method is demonstrated on two- and three-wave COVID-19 outbreaks in California, New Mexico and Florida, as observed during Summer-Winter 2020. We find that the method is robust to noise, provides useful forecasts (along with uncertainty bounds) and that it reliably detected when the initial single-wave COVID-19 outbreaks transformed into successive surges as containment efforts in these states failed by the end of Spring 2020.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Development and performance of a HemeLB GPU code for human-scale blood flow simulation

In recent years, it has become increasingly common for high performance computers (HPC) to possess some level of heterogeneous architecture - typically in the form of GPU accelerators. In some machines these are isolated within a dedicated partition, whilst in others they are integral to all compute nodes - often with multiple GPUs per node - and provide the majority of a machine's compute performance. In light of this trend, it is becoming essential that codes deployed on HPC are updated to execute on accelerator hardware. Here, in this paper, we introduce a GPU implementation of the 3D blood flow simulation code HemeLB that has been developed using CUDA C++. We demonstrate how taking advantage of NVIDIA GPU hardware can achieve significant performance improvements compared to the equivalent CPU only code on which it has been built whilst retaining the excellent strong scaling characteristics that have been repeatedly demonstrated by the CPU version. With HPC positioned on the brink of the exascale era, we use HemeLB as a motivation to provide a discussion on some of the challenges that many users will face when deploying their own applications on upcoming exascale machines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A time-varying vulnerability index for COVID-19 in New Mexico, USA using generalized propensity scores

The coronavirus disease (COVID-19) pandemic has highlighted systemic inequities in the United States and resulted in a larger burden of negative social outcomes for marginalized communities. New Mexico, a state in the southwestern US, has a unique population with a large racial minority population and a high rate of poverty that may make communities more vulnerable to negative social outcomes from COVID-19. To identify which communities may be at the highest relative risk, we created a county-level vulnerability index. After the first COVID-19 case was reported in New Mexico on March 11, 2020, we fit a generalized propensity score model that incorporates sociodemographic factors to predict county-level viral exposure and thus, the generic risk to negative social outcomes such as unemployment or mental health impacts. We used four static sociodemographic covariates important for the state of New Mexico—population, poverty, household size, and minority population—and weekly cumulative case counts to iteratively run our model each week and normalize the exposure score to create a time-varying vulnerability index. We found the relative vulnerability between counties varied in the first eight weeks from the initial COVID-19 case before stabilizing. This framework for creating a location-specific vulnerability index in response to an ongoing disaster may be used as a quick, deployable metric to inform health policy decisions such as allocating state resources to the county level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Continuous-time probabilistic models for longitudinal electronic health records

Analysis of longitudinal Electronic Health Record (EHR) data is an important goal for precision medicine. Difficulty in applying Machine Learning (ML) methods, either predictive or unsupervised, stems in part from the heterogeneity and irregular sampling of EHR data. Here, we present an unsupervised probabilistic model that captures nonlinear relationships between variables over continuous-time. This method works with arbitrary sampling patterns and captures the joint probability distribution between variable measurements and the time intervals between them. Inference algorithms are derived that can be used to evaluate the likelihood of future using under a trained model. As an example, we consider data from the United States Veterans Health Administration (VHA) in the areas of diabetes and depression. Likelihood ratio maps are produced showing the likelihood of risk for moderate-severe vs minimal depression as measured by the Patient Health Questionnaire-9 (PHQ-9).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Unsupervised probabilistic models for sequential Electronic Health Records

We develop an unsupervised probabilistic model for heterogeneous Electronic Health Record (EHR) data. Utilizing a mixture model formulation, our approach directly models sequences of arbitrary length, such as medications and laboratory results. This allows for subgrouping and incorporation of the dynamics underlying heterogeneous data types. The model consists of a layered set of latent variables that encode underlying structure in the data. These variables represent subject subgroups at the top layer, and unobserved states for sequences in the second layer. We train this model on episodic data from subjects receiving medical care in the Kaiser Permanente Northern California integrated healthcare delivery system. The resulting properties of the trained model generate novel insight from these complex and multifaceted data. In addition, we show how the model can be used to analyze sequences that contribute to assessment of mortality likelihood.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Deep learning uncertainty quantification for clinical text classification

Machine learning algorithms are expected to work side-by-side with humans in decision-making pipelines. Thus, the ability of classifiers to make reliable decisions is of paramount importance. Deep neural networks (DNNs) represent the state-of-the-art models to address real-world classification. Although the strength of activation in DNNs is often correlated with the network’s confidence, in-depth analyses are needed to establish whether they are well calibrated. In this paper, we demonstrate the use of DNN-based classification tools to benefit cancer registries by automating information extraction of disease at diagnosis and at surgery from electronic text pathology reports from the US National Cancer Institute (NCI) Surveillance, Epidemiology, and End Results (SEER) population-based cancer registries. In particular, we introduce multiple methods for selective classification to achieve a target level of accuracy on multiple classification tasks while minimizing the rejection amount—that is, the number of electronic pathology reports for which the model’s predictions are unreliable. We evaluate the proposed methods by comparing our approach with the current in-house deep learning-based abstaining classifier. Overall, all the proposed selective classification methods effectively allow for achieving the targeted level of accuracy or higher in a trade-off analysis aimed to minimize the rejection rate. On in-distribution validation and holdout test data, with all the proposed methods, we achieve on all tasks the required target level of accuracy with a lower rejection rate than the deep abstaining classifier (DAC). Interpreting the results for the out-of-distribution test data is more complex; nevertheless, in this case as well, the rejection rate from the best among the proposed methods achieving 97% accuracy or higher is lower than the rejection rate based on the DAC. We show that although both approaches can flag those samples that should be manually reviewed and labeled by human annotators, the newly proposed methods retain a larger fraction and do so without retraining—thus offering a reduced computational cost compared with the in-house deep learning-based abstaining classifier.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗
DOE OSTI2022

Machine Learning Prediction of the Critical Cooling Rate for Metallic Glasses from Expanded Datasets and Elemental Features

In this study, we use a random forest (RF) model to predict the critical cooling rate (R C ) for glass formation of various alloys from features of their constituent elements. The RF model was trained on a database that integrates multiple sources of direct and indirect R C data for metallic glasses to expand the directly measured R C database of less than 100 values to a training set of over 2000 values. The model error on 5-fold cross-validation (CV) is 0.66 orders of magnitude in K/s. The error on leave-out-one-group CV on alloy system groups is 0.59 log units in K/s when the target alloy constituents appear more than 500 times in training data. Using this model, we make predictions for the set of compositions with melt-spun glasses in the database and for the full set of quaternary alloys that have constituents which appear more than 500 times in training data. These predictions identify a number of potential new bulk metallic glass systems for future study, but the model is most useful for the identification of alloy systems likely to contain good glass formers rather than detailed discovery of bulk glass composition regions within known glassy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pose Classification Using Three-Dimensional Atomic Structure-Based Neural Networks Applied to Ion Channel–Ligand Docking

The identification of promising lead compounds showing pharmacological activities toward a biological target is essential in early stage drug discovery. With the recent increase in available small-molecule databases, virtual high-throughput screening using physics-based molecular docking has emerged as an essential tool in assisting fast and cost-efficient lead discovery and optimization. However, the best scored docking poses are often suboptimal, resulting in incorrect screening and chemical property calculation. We address the pose classification problem by leveraging data-driven machine learning approaches to identify correct docking poses from AutoDock Vina and Glide screens. To enable effective classification of docking poses, we present two convolutional neural network approaches: a three-dimensional convolutional neural network (3D-CNN) and an attention-based point cloud network (PCN) trained on the PDBbind refined set. We demonstrate the effectiveness of our proposed classifiers on multiple evaluation data sets including the standard PDBbind CASF-2016 benchmark data set and various compound libraries with structurally different protein targets including an ion channel data set extracted from Protein Data Bank (PDB) and an in-house KCa3.1 inhibitor data set. Our experiments show that excluding false positive docking poses using the proposed classifiers improves virtual high-throughput screening to identify novel molecules against each target protein compared to the initial screen based on the docking scores.

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

Behavior of Water Near Multimodal Chromatography Ligands and Its Consequences for Modulating Protein–Ligand Interactions

Multimodal chromatography is a powerful approach for purifying proteins that uses ligands containing multiple modes of interaction. Recent studies have shown that selectivity in multimodal chromatographic separations is a function of the ligand structure and geometry. Here, we performed molecular dynamics simulations to explore how the ligand structure and geometry affect ligand–water interactions and how these differences in solution affect the nature of protein–ligand interactions. Our investigation focused on three chromatography ligands: Capto MMC, Nuvia cPrime, and Prototype 4, a structural variant of Nuvia cPrime. First, the solvation characteristics of each ligand were quantified via three metrics: average water density, fluctuations, and residence time. We then explored how solvation was perturbed when the ligand was bound to the protein surface and found that the probability of the phenyl ring dewetting followed the order: Capto MMC > Prototype 4 > Nuvia cPrime. To explore how these differences in dewetting affect protein–ligand interactions, we calculated the probability of each ligand binding to different types of residues on the protein surface and found that the probability of binding to a hydrophobic residue followed the same order as the dewetting behavior. This study illustrates the role that wetting and dewetting play in modulating protein–ligand interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Tunable Heteroassembly of a Plant Pseudoenzyme–Enzyme Complex

Pseudoenzymes have emerged as key regulatory elements in all kingdoms of life despite being catalytically nonactive. Yet many factors defining why one protein is active while its homologue is inactive remain uncertain. For pseudoenzyme–enzyme pairs, the similarity of both subunits can often hinder conventional characterization approaches. In plants, a pseudoenzyme, PDX1.2, positively regulates vitamin B 6 production by association with its active catalytic homologues such as PDX1.3 through an unknown assembly mechanism. Here we used an integrative experimental approach to learn that such pseudoenzyme–enzyme pair associations result in heterocomplexes of variable stoichiometry, which are unexpectedly tunable. We also present the atomic structure of the PDX1.2 pseudoenzyme as well as the population averaged PDX1.2–PDX1.3 pseudoenzyme–enzyme pair. Finally, we dissected hetero-dodecamers of each stoichiometry to understand the arrangement of monomers in the heterocomplexes and identified symmetry-imposed preferences in PDX1.2–PDX1.3 interactions. In this work, our results provide a new model of pseudoenzyme–enzyme interactions and their native heterogeneity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗