An integrated numerical and experimental framework for modeling of CTB and GD1b ganglioside binding kinetics

Systems Engineering Approach to Modeling and Analysis of Chronic Obstructive Pulmonary Disease Part II: Extension for Variable Metabolic Rates

Recently, we developed a systems engineering model of the human cardiorespiratory system [Kurian et al. ACS Omega2023, 8 (23), 20524-20535. DOI: 10.1021/acsomega.3c00854] based on existing models of physiological processes and adapted it for chronic obstructive pulmonary disease (COPD)-an inflammatory lung disease with multiple manifestations and one of the leading causes of death in the world. This control engineering-based model is extended here to allow for variable metabolic rates established at different levels of physical activity. This required several changes to the original model: the model of the controller was enhanced to include the feedforward loop that is responsible for cardiorespiratory control under varying metabolic rates (activity level, characterized as metabolic equivalent of the task-Rm-and normalized to one at rest). In addition, a few refinements were made to the cardiorespiratory mechanics, primarily to introduce physiological processes that were not modeled earlier but became important at high metabolic rates. The extended model is verified by analyzing the impact of exercise (Rm > 1) on the cardiorespiratory system of healthy individuals. We further formally justify our previously proposed adaptation of the model for COPD patients through sensitivity analysis and refine the parameter tuning through the use of a parallel tempering stochastic global optimization method. The extended model successfully replicates experimentally observed abnormalities in COPD-the drop in arterial oxygen tension and dynamic hyperinflation under high metabolic rates-without being explicitly trained on any related data. It also supports the prospects of remote patient monitoring in COPD.

Nanocube-Based Fluidic Glycan Array

The nature of cell membrane fluidity permits glycans, which are attached to membrane proteins and lipids, to freely diffuse on cell surfaces. Through such two-dimensional motion, some weakly binding glycans can participate in lectin binding processes, eventually changing lectin binding behaviors. This chapter discusses a plasmonic nanocube sensor that allows users to detect lectin binding kinetics in a cell membrane mimicking environment. This assay only requires standard laboratory spectrometers, including microplate readers. We describe the basics of the technology in detail, including sensor fabrication, sensor calibration, data processing, a general protocol for detecting lectin-glycan interactions, and a troubleshooting guide.

A Framework for Economically Optimal Operation of Explosive Waste Incineration Process to Reduce NOx Emission Concentration

Explosives, especially those used for military weapons, have a short lifespan and their performance noticeably deteriorates over time. These old explosives need to be disposed of safely. Fluidized bed incinerators (FBIs) are safe for disposal of explosive waste (such as TNT) and produce fewer gas emissions compared to conventional methods, such as the rotary kiln. However, previous studies on this FBI process have only focused on minimizing the amount of NOx emissions without considering the operating and unitality costs (i.e., total cost) associated with the process. It is important to note that, in general, a number of different operating conditions are available to achieve a target NOx emission concentration and, thus, it requires a significant computational requirement to compare the total costs among those candidate operating conditions using a computational fluid dynamics simulation. To this end, a novel framework is proposed to quickly determine the most economically viable FBI process operating condition for a target NOx concentration. First, a surrogate model was developed to replace the high-fidelity model of an FBI process, and utilized to determine a set of possible operating conditions that may lead to a target NOx emission concentration. Second, the candidate operating conditions were fed to the Aspen Plus™ process simulation program to determine the most economically competitive option with respect to its total cost. The developed framework can provide operational guidelines for a clean and economical incineration process of explosive waste.

Development of a hybrid model for a partially known intracellular signaling pathway through correction term estimation and neural network modeling

Developing an accurate first-principle model is an important step in employing systems biology approaches to analyze an intracellular signaling pathway. However, an accurate first-principle model is difficult to be developed since it requires in-depth mechanistic understandings of the signaling pathway. Since underlying mechanisms such as the reaction network structure are not fully understood, significant discrepancy exists between predicted and actual signaling dynamics. Motivated by these considerations, this work proposes a hybrid modeling approach that combines a first-principle model and an artificial neural network (ANN) model so that predictions of the hybrid model surpass those of the original model. First, the proposed approach determines an optimal subset of model states whose dynamics should be corrected by the ANN by examining the correlation between each state and outputs through relative order. Second, an L2-regularized least-squares problem is solved to infer values of the correction terms that are necessary to minimize the discrepancy between the model predictions and available measurements. Third, an ANN is developed to generalize relationships between the values of the correction terms and the system dynamics. Lastly, the original first-principle model is coupled with the developed ANN to finalize the hybrid model development so that the model will possess generalized prediction capabilities while retaining the model interpretability. We have successfully validated the proposed methodology with two case studies, simplified apoptosis and lipopolysaccharide-induced NFκB signaling pathways, to develop hybrid models with in silico and in vitro measurements, respectively.

An intracellular signaling pathway is often represented by a set of nonlinear ordinary differential equations, which translate our current knowledge about the signaling pathway into a testable mathematical model. However, predictions from such models are often subject to high uncertainty since many signaling pathways are only partially known beforehand. In this study, we propose a systematic approach to develop a hybrid model to improve model accuracy by combining machine learning and the first-principle modeling. Specifically, model correction terms are learned from discrepancy between model predictions and measurements, and these terms are added to the first-principle model to enhance the prediction accuracy. Once these correction terms are learned from the data, an artificial neural network (ANN) model is developed to find an empirical relation between the model and the correction terms so that the developed ANN can be used to posses improved predictive capabilities even in new operating conditions (i.e., generalizability). The final hybrid model is then constructed by coupling the first-principle model with the developed ANN.

Identification of a time-varying intracellular signalling model through data clustering and parameter selection: application to NF-[inline-formula removed]B signalling pathway induced by LPS in the presence of BFA

Developing a model for a signalling pathway requires several iterations of experimentation and model refinement to obtain an accurate model. However, the implementation of such an approach to model a signalling pathway induced by a poorly-known stimulus can become labour intensive because only limited information on the pathway is available beforehand to formulate an initial model. Therefore, a large number of iterations are required since the initial model is likely to be erroneous. In this work, a numerical scheme is proposed to construct a time-varying model for a signalling pathway induced by a poorly-known stimulus when its nominal model is available in the literature. Here, the nominal model refers to one that describes the signalling dynamics under a well-characterised stimulus. First, global sensitivity analysis is implemented on the nominal model to identify the most important parameters, which are assumed to be piecewise constants. Second, measurement data are clustered to determine temporal subdomains where the parameters take different values. Finally, a least-squares problem is solved to estimate the parameter values in each temporal subdomain. The effectiveness of this approach is illustrated by developing a time-varying model for NF-[inline-formula removed]B signalling dynamics induced by lipopolysaccharide in the presence of brefeldin A.

Evaluation of hetero-multivalent lectin binding using a turbidity-based emulsion agglutination assay

Lectin hetero-multivalency, binding to two or more different types of ligands, has been demonstrated to play a role in case of both LecA (a Pseudomonas aeruginosa adhesin) and Cholera Toxin subunit B (a Vibrio cholerae toxin). In order to screen the ligand candidates that involve in hetero-multivalent binding from large molecular libraries, we present a turbidity-based emulsion agglutination (TEA) assay that can be conducted in a high-throughput format using the standard laboratory instruments and reagents. The benefit of this assay is that it relies on the use of emulsions that can be formed using ultrasonication, minimizing the bottleneck of substrate surface functionalization. By measuring the change in turbidity, we could quantify the lectin-induced aggregation rate of oil droplets to determine the relative binding strength between different ligand combinations. The TEA results are consistent with our prior binding results using a nanocube sensor. The developed TEA assay can serve as a high-throughput and customizable tool to screen the potential ligands involved in hetero-multivalent binding.

Successive complementary model-based experimental designs for parameter estimation of fed-batch bioreactors

When a dynamic model is used for the description of (fed-)batch bioreactors, it is typical that the model parameters are highly correlated to each other. In this case, it is important to keep the parameter correlation as small as possible to obtain a reliable set of parameter estimates. In this study, we propose an anticorrelation parameter estimation scheme that can be best utilized when a number of different batch experiments are sequentially processed. The scheme iteratively performs parameter estimation and model-based design of experiment (MBDOE) at the beginning and between the batches. The important difference from the existing approaches is that the MBDOE objective is defined according to the system analysis performed a priori, so that each new batch supplements what is lacking from the previous batches combined, in terms of information. The use of the scheme is illustrated on a fed-batch bioreactor model.