Search

Purpose. As a scientific field, bioinformatics has drawn remarkable attention from various fields, such as information technology, mathematics, and modern biological sciences, in recent years. The topic models originating from the field of natural language processing have become the focus of attention with the rapid accumulation of biological datasets. Thus, this research is aimed at modeling the topic content of the bioinformatics literature presented by Iranian researchers in the Scopus Citation Database. Methodology. This research was a descriptive-exploratory study, and the studied population included 3899 papers indexed in the Scopus database, which had been indexed in this database until March 9, 2022. The topic modeling was then performed on the abstracts and titles of the papers. A combination of LDA and TF-IDF was utilized for topic modeling. Findings. The data analysis with topic modeling resulted in identifying seven main topics "Molecular Modeling," "Gene Expression," "Biomarker," "Coronavirus," "Immunoinformatics," "Cancer Bioinformatics," and "Systems Biology." Moreover, "Systems Biology" and "Coronavirus" had the largest and smallest clusters, respectively. Conclusion. The present investigation demonstrated an acceptable performance for the LDA algorithm in classifying the topics included in this field. The extracted topic clusters indicated excellent consistency and topic connection with each other. [ABSTRACT FROM AUTHOR]

Malignant melanoma is one of the most threatening cancers to human health. Only 14% of patients with malignant melanoma have a remaining life span of 5 years. At present, there have been some studies looking for potential prognostic indicators of esophageal cancer from the level of genes and infiltrating immune cells, but there are still some problems that need to be resolved urgently. This paper proposes IGHG3 as the immune infiltration of malignant melanoma, which takes into account the computational mathematics. It aims to deduce the characteristics of immune cell infiltration in malignant melanoma and study the relationship between different immune cell infiltration characteristics and prognosis. The method in this article is to establish a computational mathematical model for the immunotherapy of melanoma, then study the method of identification of the affinity of the IGHG3 reagent, and finally obtain the gene expression of immune infiltration. The functions of these methods are, respectively, to predict the dynamic behavior of T cells with two different specificities through mathematical models and to test the matching degree of different concentrations of IGHG3 reagent with the human body. Then use the ssGSEA algorithm to obtain immune infiltration related data and calculate the difference between the weighted empirical cumulative distribution function of all genes in the effect of IGHG3 on melanoma that was carried out. The experimental results showed the computational mathematical method genome and all the remaining genes. In this study, a computational mathematical method to detect the IGHG3 gene expression had a significant inhibitory effect on A375 cells in the experimental group, and the knockdown efficiency reached 85.6%. [ABSTRACT FROM AUTHOR]

In this paper, an algorithm is introduced to solve the global optimization problem for melanoma skin cancer segmentation. The algorithm is based on the smoothing of an auxiliary function that is constructed using a known local minimizer and smoothed by utilising Bezier curves. This function achieves all filled function properties. The proposed optimization method is applied to find the threshold values in melanoma skin cancer images. The proposed algorithm is implemented on PH2, ISBI2016 challenge, and ISBI 2017 challenge datasets for melanoma segmentation. The results show that the proposed algorithm exhibits high accuracy, sensitivity, and specificity compared with other methods. [ABSTRACT FROM AUTHOR]

This paper proposes a method using multidomain features and support vector machine (SVM) for classifying normal and abnormal heart sound recordings. The database was provided by the PhysioNet/CinC Challenge 2016. A total of 515 features are extracted from nine feature domains, i.e., time interval, frequency spectrum of states, state amplitude, energy, frequency spectrum of records, cepstrum, cyclostationarity, high-order statistics, and entropy. Correlation analysis is conducted to quantify the feature discrimination abilities, and the results show that “frequency spectrum of state”, “energy”, and “entropy” are top domains to contribute effective features. A SVM with radial basis kernel function was trained for signal quality estimation and classification. The SVM classifier is independently trained and tested by many groups of top features. It shows the average of sensitivity, specificity, and overall score are high up to 0.88, 0.87, and 0.88, respectively, when top 400 features are used. This score is competitive to the best previous scores. The classifier has very good performance with even small number of top features for training and it has stable output regardless of randomly selected features for training. These simulations demonstrate that the proposed features and SVM classifier are jointly powerful for classifying heart sound recordings. [ABSTRACT FROM AUTHOR]

J wave is the bulge generated in the descending slope of the terminal portion of the QRS complex in the electrocardiogram. The presence of J wave may lead to sudden death. However, the diagnosis of J wave variation only depends on doctor’s clinical experiences at present and missed diagnosis is easy to occur. In this paper, a new method is proposed to realize the automatic detection of J wave. First, the synchrosqueezed wavelet transform is used to obtain the precise time-frequency information of the ECG. Then, the inverse transformation of SST is computed to get the intrinsic mode function of the ECG. At last, the time-frequency features and SST-based and the entropy features based on modes are fed to Random forest to realize the automatic detection of J wave. As the experimental results shown, the proposed method has achieved the highest accuracy, sensitivity, and specificity compared with existing techniques. [ABSTRACT FROM AUTHOR]

This paper is concerned with the mathematical modeling of a severe and common congenital defect called hypoplastic left heart syndrome (HLHS). Surgical approaches are utilized for palliating this heart condition; however, a brain white matter injury called periventricular leukomalacia (PVL) occurs with high prevalence at or around the time of surgery, the exact cause of which is not known presently. Our main goal in this paper is to study the hemodynamic conditions under which HLHS physiology may lead to the occurrence of PVL. A lumped parameter model of the HLHS circulation has been developed integrating diffusion modeling of oxygen and carbon dioxide concentrations in order to study hemodynamic variables such as pressure, flow, and blood gas concentration. Results presented include calculations of blood pressures and flow rates in different parts of the circulation. Simulations also show changes in the ratio of pulmonary to systemic blood flow rates when the sizes of the patent ductus arteriosus and atrial septal defect are varied. These changes lead to unbalanced blood circulations and, when combined with low oxygen and carbon dioxide concentrations in arteries, result in poor oxygen delivery to the brain. We stipulate that PVL occurs as a consequence. [ABSTRACT FROM AUTHOR]

Fluorescence planar imaging (FPI) is failure to capture high resolution images of deep fluorochromes due to photon diffusion. This paper presents an image restoration method to deal with this kind of blurring. The scheme of this method is conceived based on a reconstruction method in fluorescence molecular tomography (FMT) with diffusion model. A new unknown parameter is defined through introducing the first mean value theorem for definite integrals. System matrix converting this unknown parameter to the blurry image is constructed with the elements of depth conversion matrices related to a chosen plane named focal plane. Results of phantom and mouse experiments show that the proposed method is capable of reducing the blurring of FPI image caused by photon diffusion when the depth of focal plane is chosen within a proper interval around the true depth of fluorochrome. This method will be helpful to the estimation of the size of deep fluorochrome. [ABSTRACT FROM AUTHOR]

GRAPPA (Generalized Autocalibrating Partially Parallel Acquisition) is a widely used parallel MRI reconstruction technique. The processing of data from multichannel receiver coils may increase the storage and computational requirements of GRAPPA reconstruction. Random projection on GRAPPA (RP-GRAPPA) uses random projection (RP) method to overcome the computational overheads of solving large linear equations in the calibration phase of GRAPPA, saving reconstruction time. However, RP-GRAPPA compromises the reconstruction accuracy in case of large reductions in the dimensions of calibration equations. In this paper, we present the implementation of GRAPPA reconstruction method using potential iterative solvers to estimate the reconstruction coefficients from the randomly projected calibration equations. Experimental results show that the proposed methods withstand the reconstruction accuracy (visually and quantitatively) against large reductions in the dimension of linear equations, when compared with RP-GRAPPA reconstruction. Particularly, the proposed method using conjugate gradient for least squares (CGLS) demonstrates more savings in the computational time of GRAPPA, without significant loss in the reconstruction accuracy, when compared with RP-GRAPPA. It is also demonstrated that the proposed method using CGLS complements the channel compression method for reducing the computational complexities associated with higher channel count, thereby resulting in additional memory savings and speedup. [ABSTRACT FROM AUTHOR]

This paper proposed to take advantages of resonant ultrasound spectroscopy (RUS) to measure the mechanical properties of human dentin specimen. The resonant spectroscopy of the dentin specimen was obtained between the frequency bands 155 and 575 kHz, and resonant frequencies were extracted by linear predictive filter and then by Levenberg-Marquardt method. By inverse problem approach, 13 experimental resonant frequencies progressively matched to the first 30 orders of theoretical resonant frequencies calculated by Lagrangian variational method. The full second-order elastic tensor of dentin specimen was adjusted. The whole set of human dentin engineering moduli, including Young’s moduli (E11=22.641 GPa, E33=13.637 GPa), shear moduli (G12=10.608 GPa, G23=7.742 Gpa), and Poisson’s ratios (ν12=0.067, ν31=0.378), were finally calculated. This study demonstrates that RUS can be successfully adapted to measure the mechanical properties of low quality factor biomaterials. [ABSTRACT FROM AUTHOR]

This paper opens up a new issue dealing with Luffa cylindrica (LC) lignocellulosic biomass recovery in order to produce 2G bioethanol. LC fibers are composed of three principal fractions, namely, α-cellulose (45.80%  ± 1.3), hemicelluloses (20.76%  ± 0.3), and lignins (13.15%  ± 0.6). The optimization of LC fibers hydrothermal and diluted acid pretreatments duration and temperature were achieved through the cubic central composite experimental design CCD. The pretreatments optimization was monitored via the determination of reducing sugars. Then, the 2G bioethanol process feasibility was tested by means of three successive steps, namely, LC fibers hydrothermal pretreatment performed at 96°C during 54 minutes, enzymatic saccharification carried out by means of a commercial enzyme AP2, and the alcoholic fermentation fulfilled with Saccharomyces cerevisiae. LC fibers hydrothermal pretreatment liberated 33.55 g/kg of reducing sugars. Enzymatic hydrolysis allowed achieving 59.4 g/kg of reducing sugars. The conversion yield of reducing sugar to ethanol was 88.66%. After the distillation step, concentration of ethanol was 1.58% with a volumetric yield about 70%. [ABSTRACT FROM AUTHOR]

Generalized autocalibrating partially parallel acquisition (GRAPPA) has been a widely used parallel MRI technique. However, noise deteriorates the reconstructed image when reduction factor increases or even at low reduction factor for some noisy datasets. Noise, initially generated from scanner, propagates noise-related errors during fitting and interpolation procedures of GRAPPA to distort the final reconstructed image quality. The basic idea we proposed to improve GRAPPA is to remove noise from a system identification perspective. In this paper, we first analyze the GRAPPA noise problem from a noisy input-output system perspective; then, a new framework based on errors-in-variables (EIV) model is developed for analyzing noise generation mechanism in GRAPPA and designing a concrete method—instrument variables (IV) GRAPPA to remove noise. The proposed EIV framework provides possibilities that noiseless GRAPPA reconstruction could be achieved by existing methods that solve EIV problem other than IV method. Experimental results show that the proposed reconstruction algorithm can better remove the noise compared to the conventional GRAPPA, as validated with both of phantom and in vivo brain data. [ABSTRACT FROM AUTHOR]

Tumors are a serious threat to human health. The oncolytic virus is a kind of tumor killer virus which can infect and lyse cancer cells and spread through the tumor, while leaving normal cells largely unharmed. Mathematical models can help us to understand the tumor-virus dynamics and find better treatment strategies. This paper gives a new mathematical model of tumor therapy with oncolytic virus and MEK inhibitor. Stable analysis was given. Because mitogen-activated protein kinase (MEK) can not only lead to greater oncolytic virus infection into cancer cells, but also limit the replication of the virus, in order to provide the best dosage of MEK inhibitors and balance the positive and negative effect of the inhibitors, we put forward an optimal control problem of the inhibitor. The optimal strategies are given by theory and simulation. [ABSTRACT FROM AUTHOR]

Although heart researches and acquirement of clinical and experimental data are progressively open to public use, cardiac biophysical functions are still not well understood. Due to the complex and fine structures of the heart, cardiac electrophysiological features of interest may be occluded when there is a necessity to demonstrate cardiac electrophysiological behaviors. To investigate cardiac abnormal electrophysiological features under the pathological condition, in this paper, we implement a human cardiac ischemic model and acquire the electrophysiological data of excitation propagation. A visualization framework is then proposed which integrates a novel depth weighted optic attenuation model into the pathological electrophysiological model. The hidden feature of interest in pathological tissue can be revealed from sophisticated overlapping biophysical information. Experiment results verify the effectiveness of the proposed method for intuitively exploring and inspecting cardiac electrophysiological activities, which is fundamental in analyzing and explaining biophysical mechanisms of cardiac functions for doctors and medical staff. [ABSTRACT FROM AUTHOR]

Constructing rooted phylogenetic networks from rooted phylogenetic trees has become an important problem in molecular evolution. So far, many methods have been presented in this area, in which most efficient methods are based on the incompatible graph, such as the CASS, the LNETWORK, and the BIMLR. This paper will research the commonness of the methods based on the incompatible graph, the relationship between incompatible graph and the phylogenetic network, and the topologies of incompatible graphs. We can find out all the simplest datasets for a topology G and construct a network for every dataset. For any one dataset C, we can compute a network from the network representing the simplest dataset which is isomorphic to C. This process will save more time for the algorithms when constructing networks. [ABSTRACT FROM AUTHOR]

This paper presents an ultrasound simulation model for pulsatile blood flow, modulated by the motion of a stenosed vessel wall. It aims at generating more realistic ultrasonic signals to provide an environment for evaluating ultrasound signal processing and imaging and a framework for investigating the behaviors of blood flow field modulated by wall motion. This model takes into account fluid-structure interaction, blood pulsatility, stenosis of the vessel, and arterial wall movement caused by surrounding tissue’s motion. The axial and radial velocity distributions of blood and the displacement of vessel wall are calculated by solving coupled Navier-Stokes and wall equations. With these obtained values, we made several different phantoms by treating blood and the vessel wall as a group of point scatterers. Then, ultrasound echoed signals from oscillating wall and blood in the axisymmetric stenotic-carotid arteries were computed by ultrasound simulation software, Field II. The results show better consistency with corresponding theoretical values and clinical data and reflect the influence of wall movement on the flow field. It can serve as an effective tool not only for investigating the behavior of blood flow field modulated by wall motion but also for quantitative or qualitative evaluation of new ultrasound imaging technology and estimation method of blood velocity. [ABSTRACT FROM AUTHOR]

Phylogenetic networks are a generalization of phylogenetic trees that allow for the representation of evolutionary events acting at the population level, such as recombination between genes, hybridization between lineages, and horizontal gene transfer. The researchers have designed several measures for computing the dissimilarity between two phylogenetic networks, and each measure has been proven to be a metric on a special kind of phylogenetic networks. However, none of the existing measures is a metric on the space of partly reduced phylogenetic networks. In this paper, we provide a metric, de-distance, on the space of partly reduced phylogenetic networks, which is polynomial-time computable. [ABSTRACT FROM AUTHOR]

Neurons are specialized, electrically excitable cells which use electrical to chemical signals to transmit and elaborate information. Understanding how the cooperation of a great many of neurons in a grid may modify and perhaps improve the information quality, in contrast to few neurons in isolation, is critical for the rational design of cell-materials interfaces for applications in regenerative medicine, tissue engineering, and personalized lab-on-a-chips. In the present paper, we couple an integrate-and-fire model with information theory variables to analyse the extent of information in a network of nerve cells. We provide an estimate of the information in the network in bits as a function of cell density and short-term depression time. In the model, neurons are connected through a Delaunay triangulation of not-intersecting edges; in doing so, the number of connecting synapses per neuron is approximately constant to reproduce the early time of network development in planar neural cell cultures. In simulations where the number of nodes is varied, we observe an optimal value of cell density for which information in the grid is maximized. In simulations in which the posttransmission latency time is varied, we observe that information increases as the latency time decreases and, for specific configurations of the grid, it is largely enhanced in a resonance effect. [ABSTRACT FROM AUTHOR]

Cancer is a class of diseases characterized by out-of-control cells’ growth which affect cells and make them damaged. Many treatment options for cancer exist. Chemotherapy as an important treatment option is the use of drugs to treat cancer. The anticancer drug travels to the tumor and then diffuses in it through capillaries. The diffusion of drugs in the solid tumor is limited by penetration depth which is different in case of different drugs and cancers. The computation of this depth is important as it helps physicians to investigate about treatment of infected tissue. Although many efforts have been made on studying and measuring drug penetration depth, less works have been done on computing this length from a mathematical point of view. In this paper, first we propose phase lagging model for diffusion of drug in the tumor. Then, using this model on one side and considering the classic diffusion on the other side, we compute the drug penetration depth in the solid tumor. This computed value of drug penetration depth is corroborated by comparison with the values measured by experiments. [ABSTRACT FROM AUTHOR]

Protein remote homology detection is one of the central problems in bioinformatics. Although some computational methods have been proposed, the problem is still far from being solved. In this paper, an ensemble classifier for protein remote homology detection, called SVM-Ensemble, was proposed with a weighted voting strategy. SVM-Ensemble combined three basic classifiers based on different feature spaces, including Kmer, ACC, and SC-PseAAC. These features consider the characteristics of proteins from various perspectives, incorporating both the sequence composition and the sequence-order information along the protein sequences. Experimental results on a widely used benchmark dataset showed that the proposed SVM-Ensemble can obviously improve the predictive performance for the protein remote homology detection. Moreover, it achieved the best performance and outperformed other state-of-the-art methods. [ABSTRACT FROM AUTHOR]

The monoexponential model is widely used in quantitative biomedical imaging. Notable applications include apparent diffusion coefficient (ADC) imaging and pharmacokinetics. The application of ADC imaging to the detection of malignant tissue has in turn prompted several studies concerning optimal experiment design for monoexponential model fitting. In this paper, we propose a new experiment design method that is based on minimizing the determinant of the covariance matrix of the estimated parameters (D-optimal design). In contrast to previous methods, D-optimal design is independent of the imaged quantities. Applying this method to ADC imaging, we demonstrate its steady performance for the whole range of input variables (imaged parameters, number of measurements, and range of b-values). Using Monte Carlo simulations we show that the D-optimal design outperforms existing experiment design methods in terms of accuracy and precision of the estimated parameters. [ABSTRACT FROM AUTHOR]

Objective. To show how random variables concern fatigue behaviour by a probabilistic finite element method. Methods. Uncertainties on material properties due to the existence of defects that cause material elastic constant are not the same in the whole dental implant the dimensions of the structural element and load history have a decisive influence on the fatigue process and therefore on the life of a dental implant. In order to measure these uncertainties, we used a method based on Markoff chains, Bogdanoff and Kozin cumulative damage model, and probabilistic finite elements method. Results. The results have been obtained by conventional and probabilistic methods. Mathematical models obtained the same result regarding fatigue life; however, the probabilistic model obtained a greater mean life but with more information because of the cumulative probability function. Conclusions. The present paper introduces an improved procedure to study fatigue behaviour in order to know statistics of the fatigue life (mean and variance) and its probability of failure (fatigue life versus probability of failure). [ABSTRACT FROM AUTHOR]