Coral reefs, among the most diverse ecosystems on Earth, currently face major threats from pollution, unsustainable fishing practices , and perturbations in environmental parameters brought on by climate change. Corals also sustain regular wounding from other sea life and human activity. Recent reef restoration practices have even involved intentional wounding by systematically breaking coral fragments and relocating them to revitalize damaged reefs, a practice known as microfragmentation. Despite its importance, very little research has explored the inner mechanisms of wound healing in corals. Some reef-building corals have been observed to initiate an immunological response to wounding similar to that observed in mammalian species. Utilizing prior models of wound healing in mammalian species as the mathematical basis, we formulated a mechanistic model of wound healing, including observations of the immune response and tissue repair in scleractinian corals for the species Pocillopora damicornis. The model consists of four differential equations which track changes in remaining wound debris, number of cells involved in inflammation, number of cells involved in proliferation, and amount of wound closure through re-epithelialization. The model is fit to experimental wound size data from linear and circular shaped wounds on a live coral fragment. Mathematical methods, including numerical simulations and local sensitivity analysis, were used to analyze the resulting model. The parameter space was also explored to investigate drivers of other possible wound outcomes. This model serves as a first step in generating mathematical models for wound healing in corals that will not only aid in the understanding of wound healing as a whole, but also help optimize reef restoration practices and predict recovery behavior after major wounding events.

As full-scale detailed hemodynamic simulations of the entire vasculature are not feasible, numerical analysis should be focused on specific regions of the cardiovascular system, which requires the identification of lumped parameters to represent the patient behavior outside the simulated computational domain. We present a novel technique for estimating cardiovascular model parameters using gappy Proper Orthogonal Decomposition (g-POD). A POD basis is constructed with FSI simulations for different values of the lumped model parameters, and a linear operator is applied to retain information that can be compared to the available patient measurements. Then, the POD coefficients of the reconstructed solution are computed either by projecting patient measurements or by solving a minimization problem with constraints. The POD reconstruction is then used to estimate the model parameters. In the first test case, the parameter values of a 3-element Windkessel model are approximated using artificial patient measurements, obtaining a relative error of less than 4.2%. In the second case, 4 sets of 3-element Windkessel are approximated in a patient\'s aorta geometry, resulting in an error of less than 8% for the flow and less than 5% for the pressure. The method shows accurate results even with noisy patient data. It automatically calculates the delay between measurements and simulations and has flexibility in the types of patient measurements that can handle (at specific points, spatial or time averaged). The method is easy to implement and can be used in simulations performed in general-purpose FSI software.

Sulfide oxidizing bacteria are used in industrial biodesulfurization processes to convert sulfide to sulfur. These bacteria can spatially separate sulfide removal from terminal electron transfer, thereby acting as sulfide shuttles. The mechanisms underlying sulfide shuttling are not yet clear. In this work, newly obtained sulfide removal data were used to develop a new model for anaerobic sulfide removal and this model was shown to be an improvement over two previously published models. The new model describes a fast chemical step and a consecutive slow enzymatic step. The improved model includes the effect of pH, with higher total sulfide removal at increasing pH, as well as partial sulfide removal at higher sulfide concentrations. The two-stage model is supported by recent developments in anaerobic sulfide removal research and contributes to a better understanding of the underlying mechanisms. The model is a step toward accurately modelling anaerobic sulfide removal in industrial systems.

Mechanistic models are powerful tools for chromatographic process development and optimization. However, hydrophobic interaction chromatography (HIC) mechanistic models lack an effective and logical parameter estimation method, especially for multi-component system. In this study, a parameter-by-parameter method for multi-component system (called as mPbP-HIC) was derived based on the retention mechanism to estimate the six parameters of the Mollerup isotherm for HIC. The linear parameters (ks,i and keq,i) and nonlinear parameters (ni and qmax,i) of the isotherm can be estimated by the linear regression (LR) and the linear approximation (LA) steps, respectively. The remaining two parameters (kp,i and kkin,i) are obtained by the inverse method (IM). The proposed method was verified with a two-component model system. The results showed that the model could accurately predict the protein elution at a loading of 10 g/L. However, the elution curve fitting was unsatisfactory for high loadings (12 g/L and 14 g/L), which is mainly attributed to the demanding experimental conditions of the LA step and the potential large estimation error of the parameter qmax. Therefore, the inverse method was introduced to further calibrate the parameter qmax, thereby reducing the estimation error and improving the curve fitting. Moreover, the simplified linear approximation (SLA) was proposed by reasonable assumption, which provides the initial guess of qmax without solving any complex matrix and avoids the problem of matrix unsolvable. In the improved mPbP-HIC method, qmax would be initialized by the SLA and finally determined by the inverse method, and this strategy was named as SLA+IM. The experimental validation showed that the improved mPbP-HIC method has a better curve fitting, and the use of SLA+IM reduces the error accumulation effect. In process optimization, the parameters estimated by the improved mPbP-HIC method provided the model with excellent predictive ability and reasonable extrapolation. In conclusion, the SLA+IM strategy makes the improved mPbP-HIC method more rational and can be easily applied to the practical separation of protein mixture, which would accelerate the process development for HIC in downstream of biopharmaceuticals.

In-field human motion capture (HMC) is drawing increasing attention due to the multitude of application areas. Plenty of research is currently invested in camera-based (markerless) HMC, with the advantage of no infrastructure being required on the body, and additional context information being available from the surroundings. However, the inherent drawbacks of camera-based approaches are the limited field of view and occlusions. In contrast, inertial HMC (IHMC) does not suffer from occlusions, thus being a promising approach for capturing human motion outside the laboratory. However, one major challenge of such methods is the necessity of spatial registration. Typically, during a predefined calibration sequence, the orientation and location of each inertial sensor are registered with respect to the underlying skeleton model. This work contributes to calibration-free IHMC, as it proposes a recursive estimator for the simultaneous online estimation of all sensor poses and joint positions of a kinematic chain model like the human skeleton. The full derivation from an optimization objective is provided. The approach can directly be applied to a synchronized data stream from a body-mounted inertial sensor network. Successful evaluations are demonstrated on noisy simulated data from a three-link chain, real lower-body walking data from 25 young, healthy persons, and walking data captured from a humanoid robot. The estimated and derived quantities, global and relative sensor orientations, joint positions, and segment lengths can be exploited for human motion analysis and anthropometric measurements, as well as in the context of hybrid markerless visual-inertial HMC.

Chromatographic separation processes are most often modeled in the form of partial differential equations (PDEs) to describe the complex adsorption equilibria and kinetics. However, identifying parameters in such a model requires substantial computational effort. In this work, a novel parameter estimation approach using a Physics-informed Neural Network (PINN) model is developed and tested for a binary component system. Numerical accuracy of our PINN model is confirmed by validating its simulations against those of the finite element method (FEM). Furthermore, model parameters in the kinetic model are estimated by the PINN model with sufficient accuracy from the observed data at the column outlet, where parameter fitting error can be reduced by up to 35.0 % from the conventional method. In a comparison with the conventional numerical method, our approach can reduce the computational time by up to 95 %. The robustness of the PINN model has also been demonstrated by estimating model parameters from noisy artificial experimental data.

UNASSIGNED: The Philadelphia chromosome-negative myeloproliferative neoplasms are a group of slowly progressing haematological malignancies primarily characterised by an overproduction of myeloid blood cells. Patients are treated with various drugs, including the JAK1/2 inhibitor ruxolitinib. Mathematical modelling can help propose and test hypotheses of how the treatment works.
UNASSIGNED: We present an extension of the Cancitis model, which describes the development of myeloproliferative neoplasms and their interactions with inflammation, that explicitly models progenitor cells and can account for treatment with ruxolitinib through effects on the malignant stem cell response to cytokine signalling and the death rate of malignant progenitor cells. The model has been fitted to individual patients\' data for the JAK2 V617F variant allele frequency from the COMFORT-II and RESPONSE studies for patients who had substantial reductions (20 percentage points or 90% of the baseline value) in their JAK2 V617F variant allele frequency (n = 24 in total).
UNASSIGNED: The model fits very well to the patient data with an average root mean square error of 0.0249 (2.49%) when allowing ruxolitinib treatment to affect both malignant stem and progenitor cells. This average root mean square error is much lower than if allowing ruxolitinib treatment to affect only malignant stem or only malignant progenitor cells (average root mean square errors of 0.138 (13.8%) and 0.0874 (8.74%), respectively).
UNASSIGNED: Systematic simulation studies and fitting of the model to the patient data suggest that an initial reduction of the malignant cell burden followed by a monotonic increase can be recapitulated by the model assuming that ruxolitinib affects only the death rate of malignant progenitor cells. For patients exhibiting a long-term reduction of the malignant cells, the model predicts that ruxolitinib also affects stem cell parameters, such as the malignant stem cells\' response to cytokine signalling.

In this work, we introduce a compartmental model of ovarian follicle development all along lifespan, based on ordinary differential equations. The model predicts the changes in the follicle numbers in different maturation stages with aging. Ovarian follicles may either move forward to the next compartment (unidirectional migration) or degenerate and disappear (death). The migration from the first follicle compartment corresponds to the activation of quiescent follicles, which is responsible for the progressive exhaustion of the follicle reserve (ovarian aging) until cessation of reproductive activity. The model consists of a data-driven layer embedded into a more comprehensive, knowledge-driven layer encompassing the earliest events in follicle development. The data-driven layer is designed according to the most densely sampled experimental dataset available on follicle numbers in the mouse. Its salient feature is the nonlinear formulation of the activation rate, whose formulation includes a feedback term from growing follicles. The knowledge-based, coating layer accounts for cutting-edge studies on the initiation of follicle development around birth. Its salient feature is the co-existence of two follicle subpopulations of different embryonic origins. We then setup a complete estimation strategy, including the study of structural identifiability, the elaboration of a relevant optimization criterion combining different sources of data (the initial dataset on follicle numbers, together with data in conditions of perturbed activation, and data discriminating the subpopulations) with appropriate error models, and a model selection step. We finally illustrate the model potential for experimental design (suggestion of targeted new data acquisition) and in silico experiments.

Modeling studies of household transmission data have helped characterize the role of children in influenza and COVID-19 epidemics. However, estimates from these studies may be biased since they do not account for the heterogeneous nature of household contacts. Here, we quantified the impact of contact heterogeneity between household members on the estimation of child relative susceptibility and infectivity. We simulated epidemics of SARS-CoV-2-like and influenza-like infections in a synthetic population of 1,000 households assuming heterogeneous contact levels. Relative contact frequencies were derived from a household contact study according to which contacts are more frequent in the father-mother pair, followed by the child-mother, child-child, and finally child-father pairs. Child susceptibility and infectivity were then estimated while accounting for heterogeneous contacts or not. When ignoring contact heterogeneity, child relative susceptibility was underestimated by approximately 20% in the two disease scenarios. Child relative infectivity was underestimated by 20% when children and adults had different infectivity levels. These results are sensitive to our assumptions of European-style household contact patterns; but they highlight that household studies collecting both disease and contact data are needed to assess the role of complex household contact behavior on disease transmission and improve estimation of key biological parameters.

Contrary to most vessels, the ascending thoracic aorta (ATA) not only distends but also elongates in the axial direction. The purpose of this study is to investigate the biomechanical behavior of the ascending thoracic aorta (ATA) in response to dynamic axial stretching during the cardiac cycle. In addition, the implications of neglecting this dynamic axial stretching when estimating the constitutive model parameters of the ATA are investigated. The investigations were performed through in silico simulations by assuming a Gasser-Ogden-Holzapfel (GOH) constitutive model representative of ATA tissue material. The GOH model parameters were obtained from biaxial tests performed on four human ATA tissues in a previous study. Pressure-diameter curves were simulated as synthetic data to assess the effect of neglecting dynamic axial stretching on estimating constitutive model parameters. Our findings reveal a significant increase in axial stress (~ 16%) and stored strain energy (~ 18%) in the vessel when dynamic axial stretching is considered, as opposed to assuming a fixed axial stretch. All but one artery showed increased volume compliance while considering a dynamic axial stretching condition. Furthermore, we observe a notable difference in the estimated constitutive model parameters when dynamic axial stretching of the ATA is neglected, compared to the ground truth model parameters. These results underscore the critical importance of accounting for axial deformations when conducting in vivo biomechanical characterization of the ascending thoracic aorta.