Rapid analysis of sugars in fruit juices by FT-NIR spectroscopy.

A simple analytical procedure using FT-NIR and multivariate techniques for the rapid determination of individual sugars in fruit juices was evaluated. Different NIR detection devices and sample preparation methods were tested by using model solutions to determine their analytical performance. Aqueous solutions of sugar mixtures (glucose, fructose, and sucrose; 0-8% w/v) were used to develop a calibration model. Direct measurements were made by transflection using a reflectance accessory, by transmittance using a 0.5-mm cell, and by reflectance using a fiberglass paper filter. FT-NIR spectral data were transformed to the second derivative. Partial least-squares regression (PLSR) was used to create calibration models that were cross-validated (leave-one-out approach). The prediction ability of the models was evaluated on fruit juices and compared with HPLC and standard enzymatic techniques. The PLSR loading spectra showed characteristic absorption bands for the different sugars. Models generated from transmittance spectra gave the best performance with standard error of prediction (SEP) <0.10% and R(2) of 99.9% that accurately and precisely predicted the sugar levels in juices, whereas lower precision was obtained with models generated from reflectance spectra. FT-NIR spectroscopy allowed for the rapid ( approximately 3 min analysis time), accurate and non-destructive analysis of sugars in juices and could be applied in quality control of beverages or to monitor for adulteration or contamination.     

Multi-model recursive identification for nonlinear systems with non-uniformly sampling

A recursive least squares based on Multi-model is proposed for non-uniformly sampled-data nonlinear (NUSDN) systems. The corresponding state space model of an NUSDN system is derived using lifting technique. Taking advantage of the Fuzzy c-Mean Clustering algorithm, NUSDN is divided into several local models. The basic idea is that the NUSDN system is viewed as a model switching system under a given rule. Once the local models are identified, the global model is determined. A pH neutralization process validate the performance of the proposed algorithm.     

不同地形条件下植被盖度信息提取技术研究

为系统地研究特定区域的植被盖度信息提取技术, 在不同的地形条件下, 比较了目前流行的多种高光谱遥感植被盖度提取方法。结果表明: 最优高光谱归一化植被指数(NDVI₁)的建模和验证精度均高于其他两种归一化植被指数(NDVI), 直接采用NDVI建立的回归模型对研究区植被盖度的估测能力低于像元二分模型; 阴坡的最佳模型为基于一阶微分的偏最小二乘回归模型(PLSR模型), 其建模决定系数(R²)为0.810, 均方根误差(RMSE)为6.29, 验证R²为0.773, RMSE为8.85; 阳坡的最佳模型为基于二阶微分的PLSR模型, 其建模R²为0.823, RMSE为6.04, 验证R²为0.801, RMSE为7.35; 平原的最佳模型为全受限的线性光谱混合分解模型(FCLS), 其验证R²为0.852, RMSE为5.86。     

A mechanostatistical approach to cortical bone remodelling: an equine model

In this study, the development of a mechanostatistical model of three-dimensional cortical bone remodelling informed with in vivo equine data is presented. The equine model was chosen as it is highly translational to the human condition due to similar Haversian systems, availability of in vivo bone strain and biomarker data, and furthermore, equine models are recommended by the US Federal Drugs Administration for comparative joint research. The model was derived from micro-computed tomography imaged specimens taken from the equine third metacarpal bone, and the Frost-based ‘mechanostat’ was informed from both in vivo strain gauges and biomarkers to estimate bone growth rates. The model also described the well-known ‘cutting cone’ phenomena where Haversian canals tunnel and replace bone. In order to make this model useful in practice, a partial least squares regression (PLSR) surrogate model was derived based on training data from finite element simulations with different loads. The PLSR model was able to predict microstructure and homogenised Young’s modulus with errors less than 2.2 % and \(0.6\,\% \), respectively.     

Compressed images for affinity prediction-2 (CIFAP-2): an improved machine learning methodology on protein–ligand interactions based on a study on caspase 3 inhibitors

Abstract

The aim of this study is to propose an improved computational methodology, which is called Compressed Images for Affinity Prediction-2 (CIFAP-2) to predict binding affinities of structurally related protein–ligand complexes. CIFAP-2 method is established based on a protein–ligand model from which computational affinity information is obtained by utilizing 2D electrostatic potential images determined for the binding site of protein–ligand complexes. The quality of the prediction of the CIFAP-2 algorithm was tested using partial least squares regression (PLSR) as well as support vector regression (SVR) and adaptive neuro-fuzzy ?nference system (ANFIS), which are highly promising prediction methods in drug design. CIFAP-2 was applied on a protein–ligand complex system involving Caspase 3 (CASP3) and its 35 inhibitors possessing a common isatin sulfonamide pharmacophore. As a result, PLSR affinity prediction for the CASP3–ligand complexes gave rise to the most consistent information with reported empirical binding affinities (pIC₅₀) of the CASP3 inhibitors.     

Developing Subdomain Allocation Algorithms Based on Spatial and Communicational Constraints to Accelerate Dust Storm Simulation

Dust storm has serious disastrous impacts on environment, human health, and assets. The developments and applications of dust storm models have contributed significantly to better understand and predict the distribution, intensity and structure of dust storms. However, dust storm simulation is a data and computing intensive process. To improve the computing performance, high performance computing has been widely adopted by dividing the entire study area into multiple subdomains and allocating each subdomain on different computing nodes in a parallel fashion. Inappropriate allocation may introduce imbalanced task loads and unnecessary communications among computing nodes. Therefore, allocation is a key factor that may impact the efficiency of parallel process. An allocation algorithm is expected to consider the computing cost and communication cost for each computing node to minimize total execution time and reduce overall communication cost for the entire simulation. This research introduces three algorithms to optimize the allocation by considering the spatial and communicational constraints: 1) an Integer Linear Programming (ILP) based algorithm from combinational optimization perspective; 2) a K-Means and Kernighan-Lin combined heuristic algorithm (K&K) integrating geometric and coordinate-free methods by merging local and global partitioning; 3) an automatic seeded region growing based geometric and local partitioning algorithm (ASRG). The performance and effectiveness of the three algorithms are compared based on different factors. Further, we adopt the K&K algorithm as the demonstrated algorithm for the experiment of dust model simulation with the non-hydrostatic mesoscale model (NMM-dust) and compared the performance with the MPI default sequential allocation. The results demonstrate that K&K method significantly improves the simulation performance with better subdomain allocation. This method can also be adopted for other relevant atmospheric and numerical modeling.     

Advanced monitoring and control of pharmaceutical production processes with Pichia pastoris by using Raman spectroscopy and multivariate calibration methods

This contribution includes an investigation of the applicability of Raman spectroscopy as a PAT analyzer in cyclic production processes of a potential Malaria vaccine with Pichia pastoris. In a feasibility study, Partial Least Squares Regression (PLSR) models were created off‐line for cell density and concentrations of glycerol, methanol, ammonia and total secreted protein. Relative cross validation errors RMSE_cvrel range from 2.87% (glycerol) to 11.0% (ammonia). In the following, on‐line bioprocess monitoring was tested for cell density and glycerol concentration. By using the nonlinear Support Vector Regression (SVR) method instead of PLSR, the error RMSE_Prel for cell density was reduced from 5.01 to 2.94%. The high potential of Raman spectroscopy in combination with multivariate calibration methods was demonstrated by the implementation of a closed loop control for glycerol concentration using PLSR. The strong nonlinear behavior of exponentially increasing control disturbances was met with a feed‐forward control and adaptive correction of control parameters. In general the control procedure works very well for low cell densities. Unfortunately, PLSR models for glycerol concentration are strongly influenced by a correlation with the cell density. This leads to a failure in substrate prediction, which in turn prevents substrate control at cell densities above 16 g/L.     

Monitoring cellular C:N ratio in phytoplankton by means of FTIR‐spectroscopy

Statistical growth rate modelling can be applied in a variety of ecological and biotechnological applications. Such models are frequently based on Monod or Droop equations and, especially for the latter, require reliable determination of model input parameters such as C:N quotas. Besides growth rate modelling, a C:N quota quantification can be useful for monitoring and interpretation of physiological acclimation to abiotic and biotic disturbances (e.g., nutrient limitations). However, as high throughput C:N quota determination is difficult to perform, alternatives need to be established. Fourier‐transformed infrared (FTIR) spectroscopy is used to analyze a variety of biochemical, chemical, and physiological parameters in phytoplankton. Hence, a quantification of the C:N quota should also be feasible. Therefore, using FTIR spectroscopy, six phytoplankton species from among different phylogenetic groups have been analyzed to determine the effect of nutrient limitation on C:N quota patterns. The typical species‐specific response to increasing nitrogen limitation was an increase in the C:N quota. Irrespective of this species specificity, we were able to develop a reliable multi‐species C:N quota prediction model based on FTIR spectroscopy using the partial least square regression (PLSR) algorithm. Our data demonstrate that the PLSR approach is more robust in C:N quota quantification (R² = 0.93) than linear correlation of C:N quota versus growth rate (R² ranges from 0.74 to 0.86) or biochemical information based on FTIR spectra (R² ranges from 0.82 to 0.89). This accurate prediction of C:N values may support high throughput measurements in a broad range of future approaches.     

Evaluation of Multitype Mathematical Models for CFSE-Labeling Experiment Data

Carboxy-fluorescein diacetate succinimidyl ester (CFSE) labeling is an important experimental tool for measuring cell responses to extracellular signals in biomedical research. However, changes of the cell cycle (e.g., time to division) corresponding to different stimulations cannot be directly characterized from data collected in CFSE-labeling experiments. A number of independent studies have developed mathematical models as well as parameter estimation methods to better understand cell cycle kinetics based on CFSE data. However, when applying different models to the same data set, notable discrepancies in parameter estimates based on different models has become an issue of great concern. It is therefore important to compare existing models and make recommendations for practical use. For this purpose, we derived the analytic form of an age-dependent multitype branching process model. We then compared the performance of different models, namely branching process, cyton, Smith–Martin, and a linear birth–death ordinary differential equation (ODE) model via simulation studies. For fairness of model comparison, simulated data sets were generated using an agent-based simulation tool which is independent of the four models that are compared. The simulation study results suggest that the branching process model significantly outperforms the other three models over a wide range of parameter values. This model was then employed to understand the proliferation pattern of CD4+ and CD8+ T cells under polyclonal stimulation.     

不同土地利用类型下土壤有机碳含量的高光谱反演

不同土地利用类型下土壤光谱信息存在差异,了解不同土地利用类型下合适的建模方法可以高效准确地进行土壤有机碳含量反演。本研究以江西省奉新县中北部林地、耕地和园地3种土地利用类型共248个土壤样本为对象,首先对土壤原始光谱反射率曲线使用Savitzky-Golay(SG)滤波去噪并进行10 nm重采样减少数据冗余,之后采用偏最小二乘回归(PLSR)、基于网格搜索法的支持向量机回归(GRID-SVR)和基于粒子群算法的支持向量机回归(PSO-SVR)3种方法分别构建土壤有机碳含量的反演模型。结果表明: 构建单一土地利用类型反演模型时,PLSR方法在林地、耕地和园地的相对分析误差(RPD)分别为1.536、1.315和1.493,采用GRID-SVR方法时,其RPD分别提升0.150、0.183和0.502。采用PSO-SVR方法时精度最高,相较GRID-SVR方法,其林地、耕地和园地的RPD分别提高20.8%、10.0%和2.7%,林地和园地的RPD分别为2.036和2.049,可以极好地预测土壤有机碳含量,耕地的RPD为1.647,可以对土壤有机碳含量进行粗略估测。PSO-SVR方法对不同土地利用类型土壤有机碳反演效果最优,林地和园地土壤有机碳含量的反演精度相近且高于耕地。研究区不同土地利用类型对土壤有机碳含量的反演结果存在一定的影响,今后可以考虑在反演土壤有机碳时分不同土地利用类型进行建模。     

Fingerprint detection and process prediction by multivariate analysis of fed‐batch monoclonal antibody cell culture data

This work presents a sequential data analysis path, which was successfully applied to identify important patterns (fingerprints) in mammalian cell culture process data regarding process variables, time evolution and process response. The data set incorporates 116 fed‐batch cultivation experiments for the production of a Fc‐Fusion protein. Having precharacterized the evolutions of the investigated variables and manipulated parameters with univariate analysis, principal component analysis (PCA) and partial least squares regression (PLSR) are used for further investigation. The first major objective is to capture and understand the interaction structure and dynamic behavior of the process variables and the titer (process response) using different models. The second major objective is to evaluate those models regarding their capability to characterize and predict the titer production. Moreover, the effects of data unfolding, imputation of missing data, phase separation, and variable transformation on the performance of the models are evaluated. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1633–1644, 2015     

A genetic algorithm based molecular modeling technique for RNA stem-loop structures.

A new modeling technique for arriving at the three dimensional (3-D) structure of an RNA stem-loop has been developed based on a conformational search by a genetic algorithm and the following refinement by energy minimization. The genetic algorithm simultaneously optimizes a population of conformations in the predefined conformational space and generates 3-D models of RNA. The fitness function to be optimized by the algorithm has been defined to reflect the satisfaction of known conformational constraints. In addition to a term for distance constraints, the fitness function contains a term to constrain each local conformation near to a prepared template conformation. The technique has been applied to the two loops of tRNA, the anticodon loop and the T-loop, and has found good models with small root mean square deviations from the crystal structure. Slightly different models have also been found for the anticodon loop. The analysis of a collection of alternative models obtained has revealed statistical features of local variations at each base position.     

Prostate segmentation on T2 MRI using Optimal Surface Detection

This paper deals with a T2 MRI prostate segmentation method. We assume to have an initial surface mesh obtained interactively or after a first rough segmentation. The surface of the prostate is then searched within the initial mesh neighborhood using the Optimal Surface Detection algorithm (OSD). This algorithm is based on the construction of a directed graph from the information obtained around the initial mesh. The optimal surface is then obtained by a graph cut. Three different cost functions for the graph have been explored, one based on the local gradient, another on a statistical model of shape and a third on a model of gradient profile. The parameters of this method have been tuned on 33 different T2 MRI volumes.     

Modeling of the structure of protein–DNA complexes using the data from FRET and footprinting experiments

We discuss the question of constructing three-dimensional models of DNA in complex with proteins using computer modeling and indirect methods of studying the conformation of macromolecules. We consider the methods of interpreting the experimental data obtained by indirect methods of studying the three-dimensional structure of biomolecules. We discuss some aspects of integrating such data into the process of constructing the molecular models of DNA–protein complexes based on the geometric characteristics of DNA. We propose an algorithm for estimating conformations of such complexes based on the information about the local flexibility of DNA and on the experimental data obtained by Forster resonance energy transfer (FRET) and hydroxyl footprinting. Finally, we use this algorithm to predict the hypothetical configuration of DNA in a nucleosome bound with histone H1.     

Statistical estimation of statistical mechanical models: helix-coil theory and peptide helicity prediction.

Analysis of biopolymer sequences and structures generally adopts one of two approaches: use of detailed biophysical theoretical models of the system with experimentally-determined parameters, or largely empirical statistical models obtained by extracting parameters from large datasets. In this work, we demonstrate a merger of these two approaches using Bayesian statistics. We adopt a common biophysical model for local protein folding and peptide configuration, the helix-coil model. The parameters of this model are estimated by statistical fitting to a large dataset, using prior distributions based on experimental data. L(1)-norm shrinkage priors are applied to induce sparsity among the estimated parameters, resulting in a significantly simplified model. Formal statistical procedures for evaluating support in the data for previously proposed model extensions are presented. We demonstrate the advantages of this approach including improved prediction accuracy and quantification of prediction uncertainty, and discuss opportunities for statistical design of experiments. Our approach yields a 39% improvement in mean-squared predictive error over the current best algorithm for this problem. In the process we also provide an efficient recursive algorithm for exact calculation of ensemble helicity including sidechain interactions, and derive an explicit relation between homo- and heteropolymer helix-coil theories and Markov chains and (non-standard) hidden Markov models respectively, which has not appeared in the literature previously.     

Nonlinear Dispersive Model of Electroporation for Irregular Nucleated Cells

In this work, the electroporation phenomenon induced by pulsed electric field on different nucleated biological cells is studied. A nonlinear, non‐local, dispersive, and space–time multiphysics model based on Maxwell’s and asymptotic Smoluchowski’s equations has been developed to calculate the transmembrane voltage and pore density on both plasma and nuclear membrane perimeters. The irregular cell shape has been modeled by incorporating in the numerical algorithm the analytical functions pertaining to Gielis curves. The dielectric dispersion of the cell media has been modeled considering the multi‐relaxation Debye‐based relationship. Two different irregular nucleated cells have been investigated and their response has been studied applying both the dispersive and non‐dispersive models. By a comparison of the obtained results, differences can be highlighted confirming the need to make use of the dispersive model to effectively investigate the cell response in terms of transmembrane voltages, pore densities, and electroporation opening angle, especially when irregular cell shapes and short electric pulses are considered. Bioelectromagnetics. 2019;40:331–342. © 2019 Wiley Periodicals, Inc.     

Effects of Subsetting by Parent Materials on Prediction of Soil Organic Matter Content in a Hilly Area Using Vis–NIR Spectroscopy

Assessment and monitoring of soil organic matter (SOM) quality are important for understanding SOM dynamics and developing management practices that will enhance and maintain the productivity of agricultural soils. Visible and near-infrared (Vis–NIR) diffuse reflectance spectroscopy (350–2500 nm) has received increasing attention over the recent decades as a promising technique for SOM analysis. While heterogeneity of sample sets is one critical factor that complicates the prediction of soil properties from Vis–NIR spectra, a spectral library representing the local soil diversity needs to be constructed. The study area, covering a surface of 927 km² and located in Yujiang County of Jiangsu Province, is characterized by a hilly area with different soil parent materials (e.g., red sandstone, shale, Quaternary red clay, and river alluvium). In total, 232 topsoil (0–20 cm) samples were collected for SOM analysis and scanned with a Vis–NIR spectrometer in the laboratory. Reflectance data were related to surface SOM content by means of a partial least square regression (PLSR) method and several data pre-processing techniques, such as first and second derivatives with a smoothing filter. The performance of the PLSR model was tested under different combinations of calibration/validation sets (global and local calibrations stratified according to parent materials). The results showed that the models based on the global calibrations can only make approximate predictions for SOM content (RMSE (root mean squared error) = 4.23–4.69 g kg⁻¹; R² (coefficient of determination) = 0.80–0.84; RPD (ratio of standard deviation to RMSE) = 2.19–2.44; RPIQ (ratio of performance to inter-quartile distance) = 2.88–3.08). Under the local calibrations, the individual PLSR models for each parent material improved SOM predictions (RMSE = 2.55–3.49 g kg⁻¹; R² = 0.87–0.93; RPD = 2.67–3.12; RPIQ = 3.15–4.02). Among the four different parent materials, the largest R² and the smallest RMSE were observed for the shale soils, which had the lowest coefficient of variation (CV) values for clay (18.95%), free iron oxides (15.93%), and pH (1.04%). This demonstrates the importance of a practical subsetting strategy for the continued improvement of SOM prediction with Vis–NIR spectroscopy.     

Prevalence-adjusted optimisation of fuzzy habitat suitability models for aquatic invertebrate and fish species in New Zealand

For many years, habitat suitability models for aquatic species have been derived from ecological datasets by model optimisation. Previous research showed that optimisation of the predictive model performance did not necessarily lead to ecologically relevant models due to the impact of the dataset prevalence. Therefore, the adjusted average deviation was presented as a performance criterion that allowed incorporation of ecological relevance in the model optimisation process. This paper aims to analyse the relation between the adjusted average deviation (aAD) and the training set prevalence for three species in different New Zealand river systems: caddis flies Aoteapsyche spp., large brown trout Salmo trutta and rainbow trout Oncorhynchus mykiss. The aAD was implemented in a hill-climbing algorithm to optimise a fuzzy species distribution model for each species. Specifically, the hypotheses were tested that (1) similar relations between the aAD and the training set prevalence would be obtained, (2) training based on the aAD would lead to more accurate model predictions than training based on more frequently applied performance criteria such as CCI, and that (3) the final fuzzy model would produce a realistic model of habitat suitability. The approach in this paper may improve the transparency of the model training process and thus the insight into habitat suitability models. Consequently, this paper could lead to ecologically more relevant models and contribute to the implementation of these models in ecosystem management.     

Monitoring of complex industrial bioprocesses for metabolite concentrations using modern spectroscopies and machine learning: application to gibberellic acid production

Two rapid vibrational spectroscopic approaches (diffuse reflectance-absorbance Fourier transform infrared [FT-IR] and dispersive Raman spectroscopy), and one mass spectrometric method based on in vacuo Curie-point pyrolysis (PyMS), were investigated in this study. A diverse range of unprocessed, industrial fed-batch fermentation broths containing the fungus Gibberella fujikuroi producing the natural product gibberellic acid, were analyzed directly without a priori chromatographic separation. Partial least squares regression (PLSR) and artificial neural networks (ANNs) were applied to all of the information-rich spectra obtained by each of the methods to obtain quantitative information on the gibberellic acid titer. These estimates were of good precision, and the typical root-mean-square error for predictions of concentrations in an independent test set was <10% over a very wide titer range from 0 to 4925 ppm. However, although PLSR and ANNs are very powerful techniques they are often described as "black box" methods because the information they use to construct the calibration model is largely inaccessible. Therefore, a variety of novel evolutionary computation-based methods, including genetic algorithms and genetic programming, were used to produce models that allowed the determination of those input variables that contributed most to the models formed, and to observe that these models were predominantly based on the concentration of gibberellic acid itself. This is the first time that these three modern analytical spectroscopies, in combination with advanced chemometric data analysis, have been compared for their ability to analyze a real commercial bioprocess. The results demonstrate unequivocally that all methods provide very rapid and accurate estimates of the progress of industrial fermentations, and indicate that, of the three methods studied, Raman spectroscopy is the ideal bioprocess monitoring method because it can be adapted for on-line analysis.     

Evolving cell models for systems and synthetic biology

This paper proposes a new methodology for the automated design of cell models for systems and synthetic biology. Our modelling framework is based on P systems, a discrete, stochastic and modular formal modelling language. The automated design of biological models comprising the optimization of the model structure and its stochastic kinetic constants is performed using an evolutionary algorithm. The evolutionary algorithm evolves model structures by combining different modules taken from a predefined module library and then it fine-tunes the associated stochastic kinetic constants. We investigate four alternative objective functions for the fitness calculation within the evolutionary algorithm: (1) equally weighted sum method, (2) normalization method, (3) randomly weighted sum method, and (4) equally weighted product method. The effectiveness of the methodology is tested on four case studies of increasing complexity including negative and positive autoregulation as well as two gene networks implementing a pulse generator and a bandwidth detector. We provide a systematic analysis of the evolutionary algorithm’s results as well as of the resulting evolved cell models.