Traumatic brain injury (TBI) is a leading cause of disability and death. Thus, timely and effective secondary brain injury intervention is crucial, with potential to improve the prognosis of TBI. Oxidative stress contributes to post-traumatic secondary cognitive impairment, and the reduction of post-traumatic oxidative stress effectively enhances cognitive function.
View Article and Find Full Text PDFWe introduce an exploratory active learning (AL) algorithm using Gaussian process regression and marginalized graph kernel (GPR-MGK) to sample chemical compound space (CCS) at minimal cost. Targeting 251,728 enumerated alkane molecules with 4-19 carbon atoms, we applied the AL algorithm to select a diverse and representative set of molecules and then conducted high-throughput molecular simulations on these selected molecules. To demonstrate the power of the AL algorithm, we built directed message-passing neural networks (D-MPNN) using simulation data as the training set to predict liquid densities, heat capacities, and vaporization enthalpies of the CCS.
View Article and Find Full Text PDFJ Chem Inf Model
August 2023
Marginalized graph kernels have shown competitive performance in molecular machine learning tasks but currently lack measures of interpretability, which are important to improve trust in the models, detect biases, and inform molecular optimization campaigns. We here conceive and implement two interpretability measures for Gaussian process regression using a marginalized graph kernel (GPR-MGK) to quantify (1) the contribution of specific training data to the prediction and (2) the contribution of specific nodes of the graph to the prediction. We demonstrate the applicability of these interpretability measures for molecular property prediction.
View Article and Find Full Text PDFThe solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, p and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using the Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model for electrostatic interaction, the nonpolar interaction contribution of the solvent and the conformational entropy of the solute in the solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of the conformational entropy of the solute are tested.
View Article and Find Full Text PDFJ Chem Inf Model
November 2021
This work proposes a state-of-the-art hybrid kernel to calculate molecular similarity. Combined with Gaussian process models, the performance of the hybrid kernel in predicting molecular properties is comparable to that of the directed message-passing neural network (D-MPNN). The hybrid kernel consists of a marginalized graph kernel (MGK) and a radial basis function (RBF) kernel that operate on molecular graphs and global molecular features, respectively.
View Article and Find Full Text PDFThis work presents a Gaussian process regression (GPR) model on top of a novel graph representation of chemical molecules that predicts thermodynamic properties of pure substances in single, double, and triple phases. A transferable molecular graph representation is proposed as the input for a marginalized graph kernel, which is the major component of the covariance function in our GPR models. Radial basis function kernels of temperature and pressure are also incorporated into the covariance function when necessary.
View Article and Find Full Text PDFData-driven prediction of molecular properties presents unique challenges to the design of machine learning methods concerning data structure/dimensionality, symmetry adaption, and confidence management. In this paper, we present a kernel-based pipeline that can learn and predict the atomization energy of molecules with high accuracy. The framework employs Gaussian process regression to perform predictions based on the similarity between molecules, which is computed using the marginalized graph kernel.
View Article and Find Full Text PDFMolecular fingerprints, i.e., feature vectors describing atomistic neighborhood configurations, is an important abstraction and a key ingredient for data-driven modeling of potential energy surface and interatomic force.
View Article and Find Full Text PDFMesoscopic numerical simulations provide a unique approach for the quantification of the chemical influences on red blood cell functionalities. The transport Dissipative Particles Dynamics (tDPD) method can lead to such effective multiscale simulations due to its ability to simultaneously capture mesoscopic advection, diffusion, and reaction. In this paper, we present a GPU-accelerated red blood cell simulation package based on a tDPD adaptation of our red blood cell model, which can correctly recover the cell membrane viscosity, elasticity, bending stiffness, and cross-membrane chemical transport.
View Article and Find Full Text PDFWe present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment-simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles-using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation.
View Article and Find Full Text PDFWe analyze hydrodynamic fluctuations of a hybrid simulation under shear flow. The hybrid simulation is based on the Navier-Stokes (NS) equations on one domain and dissipative particle dynamics (DPD) on the other. The two domains overlap, and there is an artificial boundary for each one within the overlapping region.
View Article and Find Full Text PDFSickle-cell anaemia (SCA) is an inherited blood disorder exhibiting heterogeneous cell morphology and abnormal rheology, especially under hypoxic conditions. By using a multiscale red blood cell (RBC) model with parameters derived from patient-specific data, we present a mesoscopic computational study of the haemodynamic and rheological characteristics of blood from SCA patients with hydroxyurea (HU) treatment (on-HU) and those without HU treatment (off-HU). We determine the shear viscosity of blood in health as well as in different states of disease.
View Article and Find Full Text PDFBackground: The identification of inversions of DNA segments shorter than read length (e.g., 100 bp), defined as micro-inversions (MIs), remains challenging for next-generation sequencing reads.
View Article and Find Full Text PDFChem Commun (Camb)
July 2015
We present a non-isothermal mesoscopic model for investigation of the phase transition dynamics of thermoresponsive polymers. Since this model conserves energy in the simulations, it is able to correctly capture not only the transient behavior of polymer precipitation from solvent, but also the energy variation associated with the phase transition process. Simulations provide dynamic details of the thermally induced phase transition and confirm two different mechanisms dominating the phase transition dynamics.
View Article and Find Full Text PDFWe present large-scale simulation results on the self-assembly of amphiphilic systems in bulk solution and under soft confinement. Self-assembled unilamellar and multilamellar vesicles are formed from amphiphilic molecules in bulk solution. The system is simulated by placing amphiphilic molecules inside large unilamellar vesicles (LUVs) and the dynamic soft confinement-induced self-assembled vesicles are investigated.
View Article and Find Full Text PDF