-
Bundling instability of lophotrichous bacteria
Authors:
Jeungeun Park,
Yongsam Kim,
Wanho Lee,
Veronika Pfeifer,
Valeriia Muraveva,
Carsten Beta,
Sookkyung Lim
Abstract:
We present a mathematical model of lophotrichous bacteria, motivated by Pseudomonas putida, which swim through fluid by rotating a cluster of multiple flagella extended from near one pole of the cell body. Although the flagella rotate individually, they are typically bundled together, enabling the bacterium to exhibit three primary modes of motility: push, pull, and wrapping. One key determinant o…
▽ More
We present a mathematical model of lophotrichous bacteria, motivated by Pseudomonas putida, which swim through fluid by rotating a cluster of multiple flagella extended from near one pole of the cell body. Although the flagella rotate individually, they are typically bundled together, enabling the bacterium to exhibit three primary modes of motility: push, pull, and wrapping. One key determinant of these modes is the coordination between motor torque and rotational direction of motors. The computational variations in this coordination reveal a wide spectrum of dynamical motion regimes, which are modulated by hydrodynamic interactions between flagellar filaments. These dynamic modes can be categorized into two groups based on the collective behavior of flagella, i.e., bundled and unbundled configurations. For some of these configurations, experimental examples from fluorescence microscopy recordings of swimming P. putida cells are also presented. Furthermore, we analyze the characteristics of stable bundles, such as push and pull, and investigate the dependence of swimming behaviors on the elastic properties of the flagella.
△ Less
Submitted 23 August, 2024;
originally announced August 2024.
-
NCIDiff: Non-covalent Interaction-generative Diffusion Model for Improving Reliability of 3D Molecule Generation Inside Protein Pocket
Authors:
Joongwon Lee,
Wonho Zhung,
Woo Youn Kim
Abstract:
Advancements in deep generative modeling have changed the paradigm of drug discovery. Among such approaches, target-aware methods that exploit 3D structures of protein pockets were spotlighted for generating ligand molecules with their plausible binding modes. While docking scores superficially assess the quality of generated ligands, closer inspection of the binding structures reveals the inconsi…
▽ More
Advancements in deep generative modeling have changed the paradigm of drug discovery. Among such approaches, target-aware methods that exploit 3D structures of protein pockets were spotlighted for generating ligand molecules with their plausible binding modes. While docking scores superficially assess the quality of generated ligands, closer inspection of the binding structures reveals the inconsistency in local interactions between a pocket and generated ligands. Here, we address the issue by explicitly generating non-covalent interactions (NCIs), which are universal patterns throughout protein-ligand complexes. Our proposed model, NCIDiff, simultaneously denoises NCI types of protein-ligand edges along with a 3D graph of a ligand molecule during the sampling. With the NCI-generating strategy, our model generates ligands with more reliable NCIs, especially outperforming the baseline diffusion-based models. We further adopted inpainting techniques on NCIs to further improve the quality of the generated molecules. Finally, we showcase the applicability of NCIDiff on drug design tasks for real-world settings with specialized objectives by guiding the generation process with desired NCI patterns.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Statistical Analysis by Semiparametric Additive Regression and LSTM-FCN Based Hierarchical Classification for Computer Vision Quantification of Parkinsonian Bradykinesia
Authors:
Youngseo Cho,
In Hee Kwak,
Dohyeon Kim,
Jinhee Na,
Hanjoo Sung,
Jeongjae Lee,
Young Eun Kim,
Hyeo-il Ma
Abstract:
Bradykinesia, characterized by involuntary slowing or decrement of movement, is a fundamental symptom of Parkinson's Disease (PD) and is vital for its clinical diagnosis. Despite various methodologies explored to quantify bradykinesia, computer vision-based approaches have shown promising results. However, these methods often fall short in adequately addressing key bradykinesia characteristics in…
▽ More
Bradykinesia, characterized by involuntary slowing or decrement of movement, is a fundamental symptom of Parkinson's Disease (PD) and is vital for its clinical diagnosis. Despite various methodologies explored to quantify bradykinesia, computer vision-based approaches have shown promising results. However, these methods often fall short in adequately addressing key bradykinesia characteristics in repetitive limb movements: "occasional arrest" and "decrement in amplitude."
This research advances vision-based quantification of bradykinesia by introducing nuanced numerical analysis to capture decrement in amplitudes and employing a simple deep learning technique, LSTM-FCN, for precise classification of occasional arrests. Our approach structures the classification process hierarchically, tailoring it to the unique dynamics of bradykinesia in PD.
Statistical analysis of the extracted features, including those representing arrest and fatigue, has demonstrated their statistical significance in most cases. This finding underscores the importance of considering "occasional arrest" and "decrement in amplitude" in bradykinesia quantification of limb movement. Our enhanced diagnostic tool has been rigorously tested on an extensive dataset comprising 1396 motion videos from 310 PD patients, achieving an accuracy of 80.3%. The results confirm the robustness and reliability of our method.
△ Less
Submitted 31 March, 2024;
originally announced April 2024.
-
DeepBioisostere: Discovering Bioisosteres with Deep Learning for a Fine Control of Multiple Molecular Properties
Authors:
Hyeongwoo Kim,
Seokhyun Moon,
Wonho Zhung,
Jaechang Lim,
Woo Youn Kim
Abstract:
Optimizing molecules to improve their properties is a fundamental challenge in drug design. For a fine-tuning of molecular properties without losing bio-activity validated in advance, the concept of bioisosterism has emerged. Many in silico methods have been proposed for discovering bioisosteres, but they require expert knowledge for their applications or are restricted to known databases. Here, w…
▽ More
Optimizing molecules to improve their properties is a fundamental challenge in drug design. For a fine-tuning of molecular properties without losing bio-activity validated in advance, the concept of bioisosterism has emerged. Many in silico methods have been proposed for discovering bioisosteres, but they require expert knowledge for their applications or are restricted to known databases. Here, we introduce DeepBioisostere, a deep generative model to design suitable bioisosteric replacements. Our model allows an end-to-end chemical replacement by intelligently selecting fragments for removal and insertion along with their attachment orientation. Through various scenarios of multiple property control, we showcase the model's capability to modulate specific properties, addressing the challenge in molecular optimization. Our model's innovation lies in its capacity to design a bioisosteric replacement reflecting the compatibility with the surroundings of the modification site, facilitating the control of sophisticated properties like drug-likeness. DeepBioisostere can also provide previously unseen bioisosteric replacements, highlighting its capability for exploring diverse chemical modifications rather than just mining them from known databases. Lastly, we employed DeepBioisostere to improve the sensitivity of a known SARS-CoV-2 main protease inhibitor to the E166V mutant that exhibits drug resistance to the inhibitor, demonstrating its potential application in lead optimization.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Association between Prefrontal fNIRS signals during Cognitive tasks and College scholastic ability test (CSAT) scores: Analysis using a quantum annealing approach
Authors:
Yeaju Kim,
Junggu Choi,
Bora Kim,
Yongwan Park,
Jihyun Cha,
Jongkwan Choi,
Sanghoon Han
Abstract:
Academic achievement is a critical measure of intellectual ability, prompting extensive research into cognitive tasks as potential predictors. Neuroimaging technologies, such as functional near-infrared spectroscopy (fNIRS), offer insights into brain hemodynamics, allowing understanding of the link between cognitive performance and academic achievement. Herein, we explored the association between…
▽ More
Academic achievement is a critical measure of intellectual ability, prompting extensive research into cognitive tasks as potential predictors. Neuroimaging technologies, such as functional near-infrared spectroscopy (fNIRS), offer insights into brain hemodynamics, allowing understanding of the link between cognitive performance and academic achievement. Herein, we explored the association between cognitive tasks and academic achievement by analyzing prefrontal fNIRS signals. A novel quantum annealer (QA) feature selection algorithm was applied to fNIRS data to identify cognitive tasks correlated with CSAT scores. Twelve features (signal mean, median, variance, peak, number of peaks, sum of peaks, slope, minimum, kurtosis, skewness, standard deviation, and root mean square) were extracted from fNIRS signals at two time windows (10- and 60-second) to compare results from various feature variable conditions. The feature selection results from the QA-based and XGBoost regressor algorithms were compared to validate the former's performance. In a three-step validation process using multiple linear regression models, correlation coefficients between the feature variables and the CSAT scores, model fitness (adjusted R2), and model prediction error (RMSE) values were calculated. The quantum annealer demonstrated comparable performance to classical machine learning models, and specific cognitive tasks, including verbal fluency, recognition, and the Corsi block tapping task, were correlated with academic achievement. Group analyses revealed stronger associations between Tower of London and N-back tasks with higher CSAT scores. Quantum annealing algorithms have significant potential in feature selection using fNIRS data, and represents a novel research approach. Future studies should explore predictors of academic achievement and cognitive ability.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
PharmacoNet: Accelerating Large-Scale Virtual Screening by Deep Pharmacophore Modeling
Authors:
Seonghwan Seo,
Woo Youn Kim
Abstract:
As the size of accessible compound libraries expands to over 10 billion, the need for more efficient structure-based virtual screening methods is emerging. Different pre-screening methods have been developed for rapid screening, but there is still a lack of structure-based methods applicable to various proteins that perform protein-ligand binding conformation prediction and scoring in an extremely…
▽ More
As the size of accessible compound libraries expands to over 10 billion, the need for more efficient structure-based virtual screening methods is emerging. Different pre-screening methods have been developed for rapid screening, but there is still a lack of structure-based methods applicable to various proteins that perform protein-ligand binding conformation prediction and scoring in an extremely short time. Here, we describe for the first time a deep-learning framework for structure-based pharmacophore modeling to address this challenge. We frame pharmacophore modeling as an instance segmentation problem to determine each protein hotspot and the location of corresponding pharmacophores, and protein-ligand binding pose prediction as a graph-matching problem. PharmacoNet is significantly faster than state-of-the-art structure-based approaches, yet reasonably accurate with a simple scoring function. Furthermore, we show the promising result that PharmacoNet effectively retains hit candidates even under the high pre-screening filtration rates. Overall, our study uncovers the hitherto untapped potential of a pharmacophore modeling approach in deep learning-based drug discovery.
△ Less
Submitted 18 December, 2023; v1 submitted 1 October, 2023;
originally announced October 2023.
-
PIGNet2: A Versatile Deep Learning-based Protein-Ligand Interaction Prediction Model for Binding Affinity Scoring and Virtual Screening
Authors:
Seokhyun Moon,
Sang-Yeon Hwang,
Jaechang Lim,
Woo Youn Kim
Abstract:
Prediction of protein-ligand interactions (PLI) plays a crucial role in drug discovery as it guides the identification and optimization of molecules that effectively bind to target proteins. Despite remarkable advances in deep learning-based PLI prediction, the development of a versatile model capable of accurately scoring binding affinity and conducting efficient virtual screening remains a chall…
▽ More
Prediction of protein-ligand interactions (PLI) plays a crucial role in drug discovery as it guides the identification and optimization of molecules that effectively bind to target proteins. Despite remarkable advances in deep learning-based PLI prediction, the development of a versatile model capable of accurately scoring binding affinity and conducting efficient virtual screening remains a challenge. The main obstacle in achieving this lies in the scarcity of experimental structure-affinity data, which limits the generalization ability of existing models. Here, we propose a viable solution to address this challenge by introducing a novel data augmentation strategy combined with a physics-informed graph neural network. The model showed significant improvements in both scoring and screening, outperforming task-specific deep learning models in various tests including derivative benchmarks, and notably achieving results comparable to the state-of-the-art performance based on distance likelihood learning. This demonstrates the potential of this approach to drug discovery.
△ Less
Submitted 17 July, 2023; v1 submitted 3 July, 2023;
originally announced July 2023.
-
Geomagnetic field influences probabilistic abstract decision-making in humans
Authors:
Kwon-Seok Chae,
In-Taek Oh,
Soo Hyun Jeong,
Yong-Hwan Kim,
Soo-Chan Kim,
Yongkuk Kim
Abstract:
To resolve disputes or determine the order of things, people commonly use binary choices such as tossing a coin, even though it is obscure whether the empirical probability equals to the theoretical probability. The geomagnetic field (GMF) is broadly applied as a sensory cue for various movements in many organisms including humans, although our understanding is limited. Here we reveal a GMF-modula…
▽ More
To resolve disputes or determine the order of things, people commonly use binary choices such as tossing a coin, even though it is obscure whether the empirical probability equals to the theoretical probability. The geomagnetic field (GMF) is broadly applied as a sensory cue for various movements in many organisms including humans, although our understanding is limited. Here we reveal a GMF-modulated probabilistic abstract decision-making in humans and the underlying mechanism, exploiting the zero-sum binary stone choice of Go game as a proof-of-principle. The large-scale data analyses of professional Go matches and in situ stone choice games showed that the empirical probabilities of the stone selections were remarkably different from the theoretical probability. In laboratory experiments, experimental probability in the decision-making was significantly influenced by GMF conditions and specific magnetic resonance frequency. Time series and stepwise systematic analyses pinpointed the intentionally uncontrollable decision-making as a primary modulating target. Notably, the continuum of GMF lines and anisotropic magnetic interplay between players were crucial to influence the magnetic field resonance-mediated abstract decision-making. Our findings provide unique insights into the impact of sensing GMF in decision-makings at tipping points and the quantum mechanical mechanism for manifesting the gap between theoretical and empirical probability in 3-dimensional living space.
△ Less
Submitted 22 July, 2023; v1 submitted 28 June, 2023;
originally announced June 2023.
-
CancerGPT: Few-shot Drug Pair Synergy Prediction using Large Pre-trained Language Models
Authors:
Tianhao Li,
Sandesh Shetty,
Advaith Kamath,
Ajay Jaiswal,
Xianqian Jiang,
Ying Ding,
Yejin Kim
Abstract:
Large pre-trained language models (LLMs) have been shown to have significant potential in few-shot learning across various fields, even with minimal training data. However, their ability to generalize to unseen tasks in more complex fields, such as biology, has yet to be fully evaluated. LLMs can offer a promising alternative approach for biological inference, particularly in cases where structure…
▽ More
Large pre-trained language models (LLMs) have been shown to have significant potential in few-shot learning across various fields, even with minimal training data. However, their ability to generalize to unseen tasks in more complex fields, such as biology, has yet to be fully evaluated. LLMs can offer a promising alternative approach for biological inference, particularly in cases where structured data and sample size are limited, by extracting prior knowledge from text corpora. Our proposed few-shot learning approach uses LLMs to predict the synergy of drug pairs in rare tissues that lack structured data and features. Our experiments, which involved seven rare tissues from different cancer types, demonstrated that the LLM-based prediction model achieved significant accuracy with very few or zero samples. Our proposed model, the CancerGPT (with $\sim$ 124M parameters), was even comparable to the larger fine-tuned GPT-3 model (with $\sim$ 175B parameters). Our research is the first to tackle drug pair synergy prediction in rare tissues with limited data. We are also the first to utilize an LLM-based prediction model for biological reaction prediction tasks.
△ Less
Submitted 17 April, 2023;
originally announced April 2023.
-
Development of deep biological ages aware of morbidity and mortality based on unsupervised and semi-supervised deep learning approaches
Authors:
Seong-Eun Moon,
Ji Won Yoon,
Shinyoung Joo,
Yoohyung Kim,
Jae Hyun Bae,
Seokho Yoon,
Haanju Yoo,
Young Min Cho
Abstract:
Background: While deep learning technology, which has the capability of obtaining latent representations based on large-scale data, can be a potential solution for the discovery of a novel aging biomarker, existing deep learning methods for biological age estimation usually depend on chronological ages and lack of consideration of mortality and morbidity that are the most significant outcomes of a…
▽ More
Background: While deep learning technology, which has the capability of obtaining latent representations based on large-scale data, can be a potential solution for the discovery of a novel aging biomarker, existing deep learning methods for biological age estimation usually depend on chronological ages and lack of consideration of mortality and morbidity that are the most significant outcomes of aging. Methods: This paper proposes a novel deep learning model to learn latent representations of biological aging in regard to subjects' morbidity and mortality. The model utilizes health check-up data in addition to morbidity and mortality information to learn the complex relationships between aging and measured clinical attributes. Findings: The proposed model is evaluated on a large dataset of general populations compared with KDM and other learning-based models. Results demonstrate that biological ages obtained by the proposed model have superior discriminability of subjects' morbidity and mortality.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
PINNet: a deep neural network with pathway prior knowledge for Alzheimer's disease
Authors:
Yeojin Kim,
Hyunju Lee
Abstract:
Identification of Alzheimer's Disease (AD)-related transcriptomic signatures from blood is important for early diagnosis of the disease. Deep learning techniques are potent classifiers for AD diagnosis, but most have been unable to identify biomarkers because of their lack of interpretability. To address these challenges, we propose a pathway information-based neural network (PINNet) to predict AD…
▽ More
Identification of Alzheimer's Disease (AD)-related transcriptomic signatures from blood is important for early diagnosis of the disease. Deep learning techniques are potent classifiers for AD diagnosis, but most have been unable to identify biomarkers because of their lack of interpretability. To address these challenges, we propose a pathway information-based neural network (PINNet) to predict AD patients and analyze blood and brain transcriptomic signatures using an interpretable deep learning model. PINNet is a deep neural network (DNN) model with pathway prior knowledge from either the Gene Ontology or Kyoto Encyclopedia of Genes and Genomes databases. Then, a backpropagation-based model interpretation method was applied to reveal essential pathways and genes for predicting AD. We compared the performance of PINNet with a DNN model without a pathway. Performances of PINNet outperformed or were similar to those of DNN without a pathway using blood and brain gene expressions, respectively. Moreover, PINNet considers more AD-related genes as essential features than DNN without a pathway in the learning process. Pathway analysis of protein-protein interaction modules of highly contributed genes showed that AD-related genes in blood were enriched with cell migration, PI3K-Akt, MAPK signaling, and apoptosis in blood. The pathways enriched in the brain module included cell migration, PI3K-Akt, MAPK signaling, apoptosis, protein ubiquitination, and t-cell activation. Collectively, with prior knowledge about pathways, PINNet reveals essential pathways related to AD.
△ Less
Submitted 27 November, 2022;
originally announced November 2022.
-
Human Cognition and Language Processing with Neural-Lexicon Hypothesis
Authors:
Zang-Hee Cho,
Sun-Ha Paek,
Young-Bo Kim,
Taigyoun Cho,
Hyejin Jeong,
Haigun Lee
Abstract:
Cognition and language seem closely related to the human cognitive process, although they have not been studied and investigated in detail. Our brain is too complex to fully comprehend the structures and connectivity, as well as its functions, with the currently available technology such as electro-encephalography, positron emission tomography, or functional magnetic resonance imaging, and neurobi…
▽ More
Cognition and language seem closely related to the human cognitive process, although they have not been studied and investigated in detail. Our brain is too complex to fully comprehend the structures and connectivity, as well as its functions, with the currently available technology such as electro-encephalography, positron emission tomography, or functional magnetic resonance imaging, and neurobiological data. Therefore, the exploration of neurobiological processes, such as cognition, requires substantially more related evidences, especially from in-vivo human experiments. Cognition and language are of inter-disciplinary nature and additional methodological support is needed from other disciplines, such as deep learning in the field of artificial intelligence, for example. In this paper, we have attempted to explain the neural mechanisms underlying "cognition and language processing" or "cognition or thinking" using a novel neural network model with several newly emerging developments such as neuronal resonance, in-vivo human fiber tractography or connectivity data, Engram and Hebbian hypothesis, human memory formation in the high brain areas, deep learning, and more recently developed neural memory concepts, the neural lexicon. The neural lexicon is developed via language by repeated exposure to the neural system, similar to multilayer signal processing in deep learning. We have derived a neural model to explain how human "cognition and language processing" or "cognition and thinking" works, with a focus on language, a universal medium of the human society. Although the proposed hypothesis is not fully based on experimental evidences, a substantial portion of the observations in this study is directly and indirectly supported by recent experimental findings and the theoretical bases of deep learning research.
△ Less
Submitted 25 October, 2022; v1 submitted 24 October, 2022;
originally announced October 2022.
-
Bayesian inference of scaled versus fractional Brownian motion
Authors:
Samudrajit Thapa,
Seongyu Park,
Yeongjin Kim,
Jae-Hyung Jeon,
Ralf Metzler,
Michael A. Lomholt
Abstract:
We present a Bayesian inference scheme for scaled Brownian motion, and investigate its performance on synthetic data for parameter estimation and model selection in a combined inference with fractional Brownian motion. We include the possibility of measurement noise in both models. We find that for trajectories of a few hundred time points the procedure is able to resolve well the true model and p…
▽ More
We present a Bayesian inference scheme for scaled Brownian motion, and investigate its performance on synthetic data for parameter estimation and model selection in a combined inference with fractional Brownian motion. We include the possibility of measurement noise in both models. We find that for trajectories of a few hundred time points the procedure is able to resolve well the true model and parameters. Using the prior of the synthetic data generation process also for the inference, the approach is optimal based on decision theory. We include a comparison with inference using a prior different from the data generating one.
△ Less
Submitted 12 May, 2022; v1 submitted 30 December, 2021;
originally announced January 2022.
-
Finding Cortical Subregions Regarding the Dorsal Language Pathway Based on the Structural Connectivity
Authors:
Young-Eun Hwang,
Young-Bo Kim,
Young-Don Son
Abstract:
Although the language-related fiber pathways in the human brain, such as superior longitudinal fasciculus (SLF) and arcuate fasciculus (AF), are already well-known, understanding more sophisticated cortical regions connected by the fiber tracts is essential to scrutinizing the structural connectivity of language circuits. With the regions of interest that were selected based on the Brainnetome atl…
▽ More
Although the language-related fiber pathways in the human brain, such as superior longitudinal fasciculus (SLF) and arcuate fasciculus (AF), are already well-known, understanding more sophisticated cortical regions connected by the fiber tracts is essential to scrutinizing the structural connectivity of language circuits. With the regions of interest that were selected based on the Brainnetome atlas, the fiber orientation distribution estimation method for tractography was used to produce further elaborate connectivity information. The results indicated that both fiber bundles had two distinct connections with the prefrontal corte (PFC). The SLF-II and dorsal AF are mainly connected to the rostrodorsal part of the inferior parietal cortex (IPC) and lateral part of the fusiform gyrus with the inferior frontal junction (IFJ), respectively. In contrast, the SLF-III and ventral AF were primary linked to the anterior part of the supramarginal gyrus and superior part of the temporal cortex with the inferior frontal cortex, including the Broca's area. Moreover, the IFJ in the PFC, which has rarely been emphasized as a language-related subretion, also had the strongest connectivity with the previously known language-related subregions among the PFC; consequently, we proposed that these specific regions are interconnected via the SLF and AF within the PFC, IPC, and temporal cortex as language-related circuitry.
△ Less
Submitted 1 December, 2021;
originally announced December 2021.
-
Relational graph convolutional networks for predicting blood-brain barrier penetration of drug molecules
Authors:
Yan Ding,
Xiaoqian Jiang,
Yejin Kim
Abstract:
Evaluating the blood-brain barrier (BBB) permeability of drug molecules is a critical step in brain drug development. Traditional methods for the evaluation require complicated in vitro or in vivo testing. Alternatively, in silico predictions based on machine learning have proved to be a cost-efficient way to complement the in vitro and in vivo methods. However, the performance of the established…
▽ More
Evaluating the blood-brain barrier (BBB) permeability of drug molecules is a critical step in brain drug development. Traditional methods for the evaluation require complicated in vitro or in vivo testing. Alternatively, in silico predictions based on machine learning have proved to be a cost-efficient way to complement the in vitro and in vivo methods. However, the performance of the established models has been limited by their incapability of dealing with the interactions between drugs and proteins, which play an important role in the mechanism behind the BBB penetrating behaviors. To address this limitation, we employed the relational graph convolutional network (RGCN) to handle the drug-protein interactions as well as the properties of each individual drug. The RGCN model achieved an overall accuracy of 0.872, an AUROC of 0.919 and an AUPRC of 0.838 for the testing dataset with the drug-protein interactions and the Mordred descriptors as the input. Introducing drug-drug similarity to connect structurally similar drugs in the data graph further improved the testing results, giving an overall accuracy of 0.876, an AUROC of 0.926 and an AUPRC of 0.865. In particular, the RGCN model was found to greatly outperform the LightGBM base model when evaluated with the drugs whose BBB penetration was dependent on drug-protein interactions. Our model is expected to provide high-confidence predictions of BBB permeability for drug prioritization in the experimental screening of BBB-penetrating drugs.
△ Less
Submitted 6 April, 2022; v1 submitted 4 July, 2021;
originally announced July 2021.
-
Bayesian inference of Lévy walks via hidden Markov models
Authors:
Seongyu Park,
Samudrajit Thapa,
Yeongjin Kim,
Michael A. Lomholt,
Jae-Hyung Jeon
Abstract:
The Lévy walk is a non-Brownian random walk model that has been found to describe anomalous dynamic phenomena in diverse fields ranging from biology over quantum physics to ecology. Recurrently occurring problems are to examine whether observed data are successfully quantified by a model classified as Lévy walks or not and extract the best model parameters in accordance with the data. Motivated by…
▽ More
The Lévy walk is a non-Brownian random walk model that has been found to describe anomalous dynamic phenomena in diverse fields ranging from biology over quantum physics to ecology. Recurrently occurring problems are to examine whether observed data are successfully quantified by a model classified as Lévy walks or not and extract the best model parameters in accordance with the data. Motivated by such needs, we propose a hidden Markov model for Lévy walks and computationally realize and test the corresponding Bayesian inference method. We introduce a Markovian decomposition scheme to approximate a renewal process governed by a power-law waiting time distribution. Using this, we construct the likelihood function of Lévy walks based on a hidden Markov model and the forward algorithm. With the Lévy walk trajectories simulated at various conditions, we perform the Bayesian inference for parameter estimation and model classification. We show that the power-law exponent of the flight-time distribution can be successfully extracted even at the condition that the mean-squared displacement does not display the expected scaling exponent due to the noise or insufficient trajectory length. It is also demonstrated that the Bayesian method performs remarkably inferring the Lévy walk trajectories from given unclassified trajectory data set if the noise level is moderate.
△ Less
Submitted 8 July, 2021;
originally announced July 2021.
-
Objective comparison of methods to decode anomalous diffusion
Authors:
Gorka Muñoz-Gil,
Giovanni Volpe,
Miguel Angel Garcia-March,
Erez Aghion,
Aykut Argun,
Chang Beom Hong,
Tom Bland,
Stefano Bo,
J. Alberto Conejero,
Nicolás Firbas,
Ã’scar Garibo i Orts,
Alessia Gentili,
Zihan Huang,
Jae-Hyung Jeon,
Hélène Kabbech,
Yeongjin Kim,
Patrycja Kowalek,
Diego Krapf,
Hanna Loch-Olszewska,
Michael A. Lomholt,
Jean-Baptiste Masson,
Philipp G. Meyer,
Seongyu Park,
Borja Requena,
Ihor Smal
, et al. (9 additional authors not shown)
Abstract:
Deviations from Brownian motion leading to anomalous diffusion are ubiquitously found in transport dynamics, playing a crucial role in phenomena from quantum physics to life sciences. The detection and characterization of anomalous diffusion from the measurement of an individual trajectory are challenging tasks, which traditionally rely on calculating the mean squared displacement of the trajector…
▽ More
Deviations from Brownian motion leading to anomalous diffusion are ubiquitously found in transport dynamics, playing a crucial role in phenomena from quantum physics to life sciences. The detection and characterization of anomalous diffusion from the measurement of an individual trajectory are challenging tasks, which traditionally rely on calculating the mean squared displacement of the trajectory. However, this approach breaks down for cases of important practical interest, e.g., short or noisy trajectories, ensembles of heterogeneous trajectories, or non-ergodic processes. Recently, several new approaches have been proposed, mostly building on the ongoing machine-learning revolution. Aiming to perform an objective comparison of methods, we gathered the community and organized an open competition, the Anomalous Diffusion challenge (AnDi). Participating teams independently applied their own algorithms to a commonly-defined dataset including diverse conditions. Although no single method performed best across all scenarios, the results revealed clear differences between the various approaches, providing practical advice for users and a benchmark for developers.
△ Less
Submitted 14 May, 2021;
originally announced May 2021.
-
Drug repurposing for COVID-19 using graph neural network and harmonizing multiple evidence
Authors:
Kanglin Hsieh,
Yinyin Wang,
Luyao Chen,
Zhongming Zhao,
Sean Savitz,
Xiaoqian Jiang,
Jing Tang,
Yejin Kim
Abstract:
Amid the pandemic of 2019 novel coronavirus disease (COVID-19) infected by SARS-CoV-2, a vast amount of drug research for prevention and treatment has been quickly conducted, but these efforts have been unsuccessful thus far. Our objective is to prioritize repurposable drugs using a drug repurposing pipeline that systematically integrates multiple SARS-CoV-2 and drug interactions, deep graph neura…
▽ More
Amid the pandemic of 2019 novel coronavirus disease (COVID-19) infected by SARS-CoV-2, a vast amount of drug research for prevention and treatment has been quickly conducted, but these efforts have been unsuccessful thus far. Our objective is to prioritize repurposable drugs using a drug repurposing pipeline that systematically integrates multiple SARS-CoV-2 and drug interactions, deep graph neural networks, and in-vitro/population-based validations. We first collected all the available drugs (n= 3,635) involved in COVID-19 patient treatment through CTDbase. We built a SARS-CoV-2 knowledge graph based on the interactions among virus baits, host genes, pathways, drugs, and phenotypes. A deep graph neural network approach was used to derive the candidate representation based on the biological interactions. We prioritized the candidate drugs using clinical trial history, and then validated them with their genetic profiles, in vitro experimental efficacy, and electronic health records. We highlight the top 22 drugs including Azithromycin, Atorvastatin, Aspirin, Acetaminophen, and Albuterol. We further pinpointed drug combinations that may synergistically target COVID-19. In summary, we demonstrated that the integration of extensive interactions, deep neural networks, and rigorous validation can facilitate the rapid identification of candidate drugs for COVID-19 treatment. This is a post-peer-review, pre-copyedit version of an article published in Scientific Reports The final authenticated version is available online at: https://www.nature.com/articles/s41598-021-02353-5
△ Less
Submitted 1 February, 2022; v1 submitted 23 September, 2020;
originally announced September 2020.
-
PIGNet: A physics-informed deep learning model toward generalized drug-target interaction predictions
Authors:
Seokhyun Moon,
Wonho Zhung,
Soojung Yang,
Jaechang Lim,
Woo Youn Kim
Abstract:
Recently, deep neural network (DNN)-based drug-target interaction (DTI) models were highlighted for their high accuracy with affordable computational costs. Yet, the models' insufficient generalization remains a challenging problem in the practice of in-silico drug discovery. We propose two key strategies to enhance generalization in the DTI model. The first is to predict the atom-atom pairwise in…
▽ More
Recently, deep neural network (DNN)-based drug-target interaction (DTI) models were highlighted for their high accuracy with affordable computational costs. Yet, the models' insufficient generalization remains a challenging problem in the practice of in-silico drug discovery. We propose two key strategies to enhance generalization in the DTI model. The first is to predict the atom-atom pairwise interactions via physics-informed equations parameterized with neural networks and provides the total binding affinity of a protein-ligand complex as their sum. We further improved the model generalization by augmenting a broader range of binding poses and ligands to training data. We validated our model, PIGNet, in the comparative assessment of scoring functions (CASF) 2016, demonstrating the outperforming docking and screening powers than previous methods. Our physics-informing strategy also enables the interpretation of predicted affinities by visualizing the contribution of ligand substructures, providing insights for further ligand optimization.
△ Less
Submitted 13 December, 2021; v1 submitted 22 August, 2020;
originally announced August 2020.
-
A matter of time: Using dynamics and theory to uncover mechanisms of transcriptional bursting
Authors:
Nicholas C. Lammers,
Yang Joon Kim,
Jiaxi Zhao,
Hernan G. Garcia
Abstract:
Eukaryotic transcription generally occurs in bursts of activity lasting minutes to hours; however, state-of-the-art measurements have revealed that many of the molecular processes that underlie bursting, such as transcription factor binding to DNA, unfold on timescales of seconds. This temporal disconnect lies at the heart of a broader challenge in physical biology of predicting transcriptional ou…
▽ More
Eukaryotic transcription generally occurs in bursts of activity lasting minutes to hours; however, state-of-the-art measurements have revealed that many of the molecular processes that underlie bursting, such as transcription factor binding to DNA, unfold on timescales of seconds. This temporal disconnect lies at the heart of a broader challenge in physical biology of predicting transcriptional outcomes and cellular decision-making from the dynamics of underlying molecular processes. Here, we review how new dynamical information about the processes underlying transcriptional control can be combined with theoretical models that predict not only averaged transcriptional dynamics, but also their variability, to formulate testable hypotheses about the molecular mechanisms underlying transcriptional bursting and control.
△ Less
Submitted 17 December, 2020; v1 submitted 20 August, 2020;
originally announced August 2020.
-
Population stratification enables modeling effects of reopening policies on mortality and hospitalization rates
Authors:
Tongtong Huang,
Yan Chu,
Shayan Shams,
Yejin Kim,
Genevera Allen,
Ananth V Annapragada,
Devika Subramanian,
Ioannis Kakadiaris,
Assaf Gottlieb,
Xiaoqian Jiang
Abstract:
Objective: We study the influence of local reopening policies on the composition of the infectious population and their impact on future hospitalization and mortality rates. Materials and Methods: We collected datasets of daily reported hospitalization and cumulative morality of COVID 19 in Houston, Texas, from May 1, 2020 until June 29, 2020. These datasets are from multiple sources (USA FACTS, S…
▽ More
Objective: We study the influence of local reopening policies on the composition of the infectious population and their impact on future hospitalization and mortality rates. Materials and Methods: We collected datasets of daily reported hospitalization and cumulative morality of COVID 19 in Houston, Texas, from May 1, 2020 until June 29, 2020. These datasets are from multiple sources (USA FACTS, Southeast Texas Regional Advisory Council COVID 19 report, TMC daily news, and New York Times county level mortality reporting). Our model, risk stratified SIR HCD uses separate variables to model the dynamics of local contact (e.g., work from home) and high contact (e.g., work on site) subpopulations while sharing parameters to control their respective $R_0(t)$ over time. Results: We evaluated our models forecasting performance in Harris County, TX (the most populated county in the Greater Houston area) during the Phase I and Phase II reopening. Not only did our model outperform other competing models, it also supports counterfactual analysis to simulate the impact of future policies in a local setting, which is unique among existing approaches. Discussion: Local mortality and hospitalization are significantly impacted by quarantine and reopening policies. No existing model has directly accounted for the effect of these policies on local trends in infections, hospitalizations, and deaths in an explicit and explainable manner. Our work is an attempt to close this important technical gap to support decision making. Conclusion: Despite several limitations, we think it is a timely effort to rethink about how to best model the dynamics of pandemics under the influence of reopening policies.
△ Less
Submitted 10 August, 2020;
originally announced August 2020.
-
Efficient Reconstruction of Stochastic Pedigrees
Authors:
Younhun Kim,
Elchanan Mossel,
Govind Ramnarayan,
Paxton Turner
Abstract:
We introduce a new algorithm called {\sc Rec-Gen} for reconstructing the genealogy or \textit{pedigree} of an extant population purely from its genetic data. We justify our approach by giving a mathematical proof of the effectiveness of {\sc Rec-Gen} when applied to pedigrees from an idealized generative model that replicates some of the features of real-world pedigrees. Our algorithm is iterative…
▽ More
We introduce a new algorithm called {\sc Rec-Gen} for reconstructing the genealogy or \textit{pedigree} of an extant population purely from its genetic data. We justify our approach by giving a mathematical proof of the effectiveness of {\sc Rec-Gen} when applied to pedigrees from an idealized generative model that replicates some of the features of real-world pedigrees. Our algorithm is iterative and provides an accurate reconstruction of a large fraction of the pedigree while having relatively low \emph{sample complexity}, measured in terms of the length of the genetic sequences of the population. We propose our approach as a prototype for further investigation of the pedigree reconstruction problem toward the goal of applications to real-world examples. As such, our results have some conceptual bearing on the increasingly important issue of genomic privacy.
△ Less
Submitted 7 May, 2020;
originally announced May 2020.
-
Optimal transport and barycenters for dendritic measures
Authors:
Young-Heon Kim,
Brendan Pass,
David J. Schneider
Abstract:
We introduce and study a variant of the Wasserstein distance on the space of probability measures, specially designed to deal with measures whose support has a dendritic, or treelike structure with a particular direction of orientation. Our motivation is the comparison of and interpolation between plants' root systems. We characterize barycenters with respect to this metric, and establish that the…
▽ More
We introduce and study a variant of the Wasserstein distance on the space of probability measures, specially designed to deal with measures whose support has a dendritic, or treelike structure with a particular direction of orientation. Our motivation is the comparison of and interpolation between plants' root systems. We characterize barycenters with respect to this metric, and establish that the interpolations of root-like measures, using this new metric, are also root like, in a certain sense; this property fails for conventional Wasserstein barycenters. We also establish geodesic convexity with respect to this metric for a variety of functionals, some of which we expect to have biological importance.
△ Less
Submitted 17 October, 2019;
originally announced October 2019.
-
Scaffold-based molecular design using graph generative model
Authors:
Jaechang Lim,
Sang-Yeon Hwang,
Seungsu Kim,
Seokhyun Moon,
Woo Youn Kim
Abstract:
Searching new molecules in areas like drug discovery often starts from the core structures of candidate molecules to optimize the properties of interest. The way as such has called for a strategy of designing molecules retaining a particular scaffold as a substructure. On this account, our present work proposes a scaffold-based molecular generative model. The model generates molecular graphs by ex…
▽ More
Searching new molecules in areas like drug discovery often starts from the core structures of candidate molecules to optimize the properties of interest. The way as such has called for a strategy of designing molecules retaining a particular scaffold as a substructure. On this account, our present work proposes a scaffold-based molecular generative model. The model generates molecular graphs by extending the graph of a scaffold through sequential additions of vertices and edges. In contrast to previous related models, our model guarantees the generated molecules to retain the given scaffold with certainty. Our evaluation of the model using unseen scaffolds showed the validity, uniqueness, and novelty of generated molecules as high as the case using seen scaffolds. This confirms that the model can generalize the learned chemical rules of adding atoms and bonds rather than simply memorizing the mapping from scaffolds to molecules during learning. Furthermore, despite the restraint of fixing core structures, our model could simultaneously control multiple molecular properties when generating new molecules.
△ Less
Submitted 31 May, 2019;
originally announced May 2019.
-
From Brain Imaging to Graph Analysis: a study on ADNI's patient cohort
Authors:
Rui Zhang,
Luca Giancardo,
Danilo A. Pena,
Yejin Kim,
Hanghang Tong,
Xiaoqian Jiang
Abstract:
In this paper, we studied the association between the change of structural brain volumes to the potential development of Alzheimer's disease (AD). Using a simple abstraction technique, we converted regional cortical and subcortical volume differences over two time points for each study subject into a graph. We then obtained substructures of interest using a graph decomposition algorithm in order t…
▽ More
In this paper, we studied the association between the change of structural brain volumes to the potential development of Alzheimer's disease (AD). Using a simple abstraction technique, we converted regional cortical and subcortical volume differences over two time points for each study subject into a graph. We then obtained substructures of interest using a graph decomposition algorithm in order to extract pivotal nodes via multi-view feature selection. Intensive experiments using robust classification frameworks were conducted to evaluate the performance of using the brain substructures obtained under different thresholds. The results indicated that compact substructures acquired by examining the differences between patient groups were sufficient to discriminate between AD and healthy controls with an area under the receiver operating curve of 0.72.
△ Less
Submitted 14 May, 2019;
originally announced May 2019.
-
Discriminative Sleep Patterns of Alzheimer's Disease via Tensor Factorization
Authors:
Yejin Kim,
Xiaoqian Jiang,
Luyao Chen,
Xiaojin Li,
Licong Cui
Abstract:
Sleep change is commonly reported in Alzheimer's disease (AD) patients and their brain wave studies show decrease in dreaming and non-dreaming stages. Although sleep disturbance is generally considered as a consequence of AD, it might also be a risk factor of AD as new biological evidence shows. Leveraging National Sleep Research Resource (NSRR), we built a unique cohort of 83 cases and 331 contro…
▽ More
Sleep change is commonly reported in Alzheimer's disease (AD) patients and their brain wave studies show decrease in dreaming and non-dreaming stages. Although sleep disturbance is generally considered as a consequence of AD, it might also be a risk factor of AD as new biological evidence shows. Leveraging National Sleep Research Resource (NSRR), we built a unique cohort of 83 cases and 331 controls with clinical variables and EEG signals. Supervised tensor factorization method was applied for this temporal dataset to extract discriminative sleep patterns. Among the 30 patterns extracted, we identified 5 significant patterns (4 patterns for AD likely and 1 pattern for normal ones) and their visual patterns provide interesting linkage to sleep with repeated wakefulness, insomnia, epileptic seizure, and etc. This study is preliminary but findings are interesting, which is a first step to provide quantifiable evidences to measure sleep as a risk factor of AD.
△ Less
Submitted 14 May, 2019;
originally announced May 2019.
-
Repurposing of Ring Framework through TR screening
Authors:
Surendra Kumar,
Cheongyun Jang,
Lalita Subedi,
Sun Yeou Kim,
Mi-hyun Kim
Abstract:
In the drug repurposing approach, the chemically diverse and potentially safe molecules can be explored as therapeutic potential out of those originally targeted indications. However, the intellectual property rights, and competitive over-heating issue from sharing the same drug pipeline information are drawback of the approach. In current paradigm of drug discovery, despite approval of new chemic…
▽ More
In the drug repurposing approach, the chemically diverse and potentially safe molecules can be explored as therapeutic potential out of those originally targeted indications. However, the intellectual property rights, and competitive over-heating issue from sharing the same drug pipeline information are drawback of the approach. In current paradigm of drug discovery, despite approval of new chemical entity, the finding of a new ring system for drug discovery have stagnated. In addition, there are still rarely used ring systems from FDA approved drug that can be optimized and used for repurposing purpose. In the present work, these rarely used drug ring system were re-purposed with our bidirectional target and ring screening (TR- dual screening) approach within our defined objectives. The TR-dual screening suggested RARbeta and cyproheptadine as a pair of target and ring system. The selected ring system was virtually grown with predefined criterion, ranking of molecular docking and top score poses made us select test compounds for synthesis and biological proof. Among the tested compounds, 6c proved RARbeta agonist activity. The concept of TR dual optimization can contribute to picking either a compound having acceptable ADMET profile or a target for a selected drug scaffold.
△ Less
Submitted 11 April, 2019;
originally announced April 2019.
-
Waveforms of molecular oscillations reveal circadian timekeeping mechanisms
Authors:
Hang-Hyun Jo,
Yeon Jeong Kim,
Jae Kyoung Kim,
Mathias Foo,
David E. Somers,
Pan-Jun Kim
Abstract:
Circadian clocks play a pivotal role in orchestrating numerous physiological and developmental events. Waveform shapes of the oscillations of protein abundances can be informative about the underlying biochemical processes of circadian clocks. We derive a mathematical framework where waveforms do reveal hidden biochemical mechanisms of circadian timekeeping. We find that the cost of synthesizing p…
▽ More
Circadian clocks play a pivotal role in orchestrating numerous physiological and developmental events. Waveform shapes of the oscillations of protein abundances can be informative about the underlying biochemical processes of circadian clocks. We derive a mathematical framework where waveforms do reveal hidden biochemical mechanisms of circadian timekeeping. We find that the cost of synthesizing proteins with particular waveforms can be substantially reduced by rhythmic protein half-lives over time, as supported by previous plant and mammalian data, as well as our own seedling experiment. We also find that previously-enigmatic, cyclic expression of positive arm components within the mammalian and insect clocks allows both a broad range of peak time differences between protein waveforms and the symmetries of the waveforms about the peak times. Such various peak-time differences may facilitate tissue-specific or developmental stage-specific multicellular processes. Our waveform-guided approach can be extended to various biological oscillators, including cell-cycle and synthetic genetic oscillators.
△ Less
Submitted 26 November, 2018;
originally announced November 2018.
-
Harmonic holes as the submodules of brain network and network dissimilarity
Authors:
Hyekyoung Lee,
Moo K. Chung,
Hongyoon Choi,
Hyejin Kang,
Seunggyun Ha,
Yu Kyeong Kim,
Dong Soo Lee
Abstract:
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connec…
▽ More
Persistent homology has been applied to brain network analysis for finding the shape of brain networks across multiple thresholds. In the persistent homology, the shape of networks is often quantified by the sequence of $k$-dimensional holes and Betti numbers.The Betti numbers are more widely used than holes themselves in topological brain network analysis. However, the holes show the local connectivity of networks, and they can be very informative features in analysis. In this study, we propose a new method of measuring network differences based on the dissimilarity measure of harmonic holes (HHs). The HHs, which represent the substructure of brain networks, are extracted by the Hodge Laplacian of brain networks. We also find the most contributed HHs to the network difference based on the HH dissimilarity. We applied our proposed method to clustering the networks of 4 groups, normal control (NC), stable and progressive mild cognitive impairment (sMCI and pMCI), and Alzheimer's disease (AD). The results showed that the clustering performance of the proposed method was better than that of network distances based on only the global change of topology.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.
-
How Many Subpopulations is Too Many? Exponential Lower Bounds for Inferring Population Histories
Authors:
Younhun Kim,
Frederic Koehler,
Ankur Moitra,
Elchanan Mossel,
Govind Ramnarayan
Abstract:
Reconstruction of population histories is a central problem in population genetics. Existing coalescent-based methods, like the seminal work of Li and Durbin (Nature, 2011), attempt to solve this problem using sequence data but have no rigorous guarantees. Determining the amount of data needed to correctly reconstruct population histories is a major challenge. Using a variety of tools from informa…
▽ More
Reconstruction of population histories is a central problem in population genetics. Existing coalescent-based methods, like the seminal work of Li and Durbin (Nature, 2011), attempt to solve this problem using sequence data but have no rigorous guarantees. Determining the amount of data needed to correctly reconstruct population histories is a major challenge. Using a variety of tools from information theory, the theory of extremal polynomials, and approximation theory, we prove new sharp information-theoretic lower bounds on the problem of reconstructing population structure -- the history of multiple subpopulations that merge, split and change sizes over time. Our lower bounds are exponential in the number of subpopulations, even when reconstructing recent histories. We demonstrate the sharpness of our lower bounds by providing algorithms for distinguishing and learning population histories with matching dependence on the number of subpopulations. Along the way and of independent interest, we essentially determine the optimal number of samples needed to learn an exponential mixture distribution information-theoretically, proving the upper bound by analyzing natural (and efficient) algorithms for this problem.
△ Less
Submitted 8 May, 2019; v1 submitted 7 November, 2018;
originally announced November 2018.
-
Mapping the spectrum of 3D communities in human chromosome conformation capture data
Authors:
Sang Hoon Lee,
Yeonghoon Kim,
Sungmin Lee,
Xavier Durang,
Per Stenberg,
Jae-Hyung Jeon,
Ludvig Lizana
Abstract:
Several experiments show that the three dimensional (3D) organization of chromosomes affects genetic processes such as transcription and gene regulation. To better understand this connection, researchers developed the Hi-C method that is able to detect the pairwise physical contacts of all chromosomal loci. The Hi-C data show that chromosomes are composed of 3D compartments that range over a varie…
▽ More
Several experiments show that the three dimensional (3D) organization of chromosomes affects genetic processes such as transcription and gene regulation. To better understand this connection, researchers developed the Hi-C method that is able to detect the pairwise physical contacts of all chromosomal loci. The Hi-C data show that chromosomes are composed of 3D compartments that range over a variety of scales. However, it is challenging to systematically detect these cross-scale structures. Most studies have therefore designed methods for specific scales to study foremost topologically associated domains (TADs) and A/B compartments. To go beyond this limitation, we tailor a network community detection method that finds communities in compact fractal globule polymer systems. Our method allows us to continuously scan through all scales with a single resolution parameter. We found: (i) polymer segments belonging to the same 3D community do not have to be in consecutive order along the polymer chain. In other words, several TADs may belong to the same 3D community. (ii) CTCF proteins---a loop-stabilizing protein that is ascribed a big role in TAD formation---are well correlated with community borders only at one level of organization. (iii) TADs and A/B compartments are traditionally treated as two weakly related 3D structures and detected with different algorithms. With our method, we detect both by simply adjusting the resolution parameter. We therefore argue that they represent two specific levels of a continuous spectrum 3D communities, rather than seeing them as different structural entities.
△ Less
Submitted 2 May, 2019; v1 submitted 2 October, 2018;
originally announced October 2018.
-
Contribution of Data Categories to Readmission Prediction Accuracy
Authors:
Wendong Ge,
Hee Yeun Kim,
Sonali Desai,
Leonid Perlovsky,
Alexander Turchin
Abstract:
Identification of patients at high risk for readmission could help reduce morbidity and mortality as well as healthcare costs. Most of the existing studies on readmission prediction did not compare the contribution of data categories. In this study we analyzed relative contribution of 90,101 variables across 398,884 admission records corresponding to 163,468 patients, including patient demographic…
▽ More
Identification of patients at high risk for readmission could help reduce morbidity and mortality as well as healthcare costs. Most of the existing studies on readmission prediction did not compare the contribution of data categories. In this study we analyzed relative contribution of 90,101 variables across 398,884 admission records corresponding to 163,468 patients, including patient demographics, historical hospitalization information, discharge disposition, diagnoses, procedures, medications and laboratory test results. We established an interpretable readmission prediction model based on Logistic Regression in scikit-learn, and added the available variables to the model one by one in order to analyze the influences of individual data categories on readmission prediction accuracy. Diagnosis related groups (c-statistic increment of 0.0933) and discharge disposition (c-statistic increment of 0.0269) were the strongest contributors to model accuracy. Additionally, we also identified the top ten contributing variables in every data category.
△ Less
Submitted 22 March, 2018; v1 submitted 21 March, 2018;
originally announced March 2018.
-
BeWith: A Between-Within Method to Discover Relationships between Cancer Modules via Integrated Analysis of Mutual Exclusivity, Co-occurrence and Functional Interactions
Authors:
Phuong Dao,
Yoo-Ah Kim,
Damian Wojtowicz,
Sanna Madan,
Roded Sharan,
Teresa M. Przytycka
Abstract:
The analysis of the mutational landscape of cancer, including mutual exclusivity and co-occurrence of mutations, has been instrumental in studying the disease. We hypothesized that exploring the interplay between co-occurrence, mutual exclusivity, and functional interactions between genes will further improve our understanding of the disease and help to uncover new relations between cancer driving…
▽ More
The analysis of the mutational landscape of cancer, including mutual exclusivity and co-occurrence of mutations, has been instrumental in studying the disease. We hypothesized that exploring the interplay between co-occurrence, mutual exclusivity, and functional interactions between genes will further improve our understanding of the disease and help to uncover new relations between cancer driving genes and pathways. To this end, we designed a general framework, BeWith, for identifying modules with different combinations of mutation and interaction patterns. We focused on three different settings of the BeWith schema: (i) BeME-WithFun in which the relations between modules are enriched with mutual exclusivity while genes within each module are functionally related; (ii) BeME-WithCo which combines mutual exclusivity between modules with co-occurrence within modules; and (iii) BeCo-WithMEFun which ensures co-occurrence between modules while the within module relations combine mutual exclusivity and functional interactions. We formulated the BeWith framework using Integer Linear Programming (ILP), enabling us to find optimally scoring sets of modules. Our results demonstrate the utility of BeWith in providing novel information about mutational patterns, driver genes, and pathways. In particular, BeME-WithFun helped identify functionally coherent modules that might be relevant for cancer progression. In addition to finding previously well-known drivers, the identified modules pointed to the importance of the interaction between NCOR and NCOA3 in breast cancer. Additionally, an application of the BeME-WithCo setting revealed that gene groups differ with respect to their vulnerability to different mutagenic processes, and helped us to uncover pairs of genes with potentially synergetic effects, including a potential synergy between mutations in TP53 and metastasis related DCC gene.
△ Less
Submitted 28 April, 2017;
originally announced April 2017.
-
WeSME: Uncovering Mutual Exclusivity of Cancer Drivers and Beyond
Authors:
Yoo-Ah Kim,
Sanna Madan,
Teresa M. Przytycka
Abstract:
Mutual exclusivity is a widely recognized property of many cancer drivers. Knowledge about these relationships can provide important insights into cancer drivers, cancer-driving pathways, and cancer subtypes. It can also be used to predict new functional interactions between cancer driving genes and uncover novel cancer drivers. Currently, most of mutual exclusivity analyses are preformed focusing…
▽ More
Mutual exclusivity is a widely recognized property of many cancer drivers. Knowledge about these relationships can provide important insights into cancer drivers, cancer-driving pathways, and cancer subtypes. It can also be used to predict new functional interactions between cancer driving genes and uncover novel cancer drivers. Currently, most of mutual exclusivity analyses are preformed focusing on a limited set of genes in part due to the computational cost required to rigorously compute p-values. To reduce the computing cost and perform less restricted mutual exclusivity analysis, we developed an efficient method to estimate p-values while controlling the mutation rates of individual patients and genes similar to the permutation test. A comprehensive mutual exclusivity analysis allowed us to uncover mutually exclusive pairs, some of which may have relatively low mutation rates. These pairs often included likely cancer drivers that have been missed in previous analyses. More importantly, our results demonstrated that mutual exclusivity can also provide information that goes beyond the interactions between cancer drivers and can, for example, elucidate different mutagenic processes in different cancer groups. In particular, including frequently mutated, long genes such as TTN in our analysis allowed us to observe interesting patterns of APOBEC activity in breast cancer and identify a set of related driver genes that are highly predictive of patient survival. In addition, we utilized our mutual exclusivity analysis in support of a previously proposed model where APOBEC activity is the underlying process that causes TP53 mutations in a subset of breast cancer cases.
△ Less
Submitted 7 April, 2016;
originally announced April 2016.
-
NASCUP: Nucleic Acid Sequence Classification by Universal Probability
Authors:
Sunyoung Kwon,
Gyuwan Kim,
Byunghan Lee,
Jongsik Chun,
Sungroh Yoon,
Young-Han Kim
Abstract:
Motivated by the need for fast and accurate classification of unlabeled nucleotide sequences on a large scale, we developed NASCUP, a new classification method that captures statistical structures of nucleotide sequences by compact context-tree models and universal probability from information theory. NASCUP achieved BLAST-like classification accuracy consistently for several large-scale databases…
▽ More
Motivated by the need for fast and accurate classification of unlabeled nucleotide sequences on a large scale, we developed NASCUP, a new classification method that captures statistical structures of nucleotide sequences by compact context-tree models and universal probability from information theory. NASCUP achieved BLAST-like classification accuracy consistently for several large-scale databases in orders-of-magnitude reduced runtime, and was applied to other bioinformatics tasks such as outlier detection and synthetic sequence generation.
△ Less
Submitted 29 November, 2018; v1 submitted 16 November, 2015;
originally announced November 2015.
-
Single-molecule study on the orientation of coiled-coil proteins attached to gold nanostructures
Authors:
Jae-Byum Chang,
Yong Ho Kim,
Evan Thompson,
Young Hyun No,
Nam Hyeong Kim,
Jose Arrieta,
Vitor R. Manfrinato,
Amy E. Keating,
Karl K. Berggren
Abstract:
Methods for patterning biomolecules on a substrate at the single molecule level have been studied as a route to sensors with single-molecular sensitivity or as a way to probe biological phenomena at the single-molecule level. However, the arrangement and orientation of single biomolecules on substrates has been less investigated. Here, we examined the arrangement and orientation of two rod-like co…
▽ More
Methods for patterning biomolecules on a substrate at the single molecule level have been studied as a route to sensors with single-molecular sensitivity or as a way to probe biological phenomena at the single-molecule level. However, the arrangement and orientation of single biomolecules on substrates has been less investigated. Here, we examined the arrangement and orientation of two rod-like coiled-coil proteins, cortexillin and tropomyosin, around patterned gold nanostructures. The high aspect ratio of the coiled coils made it possible to study their orientations and to pursue a strategy of protein orientation via two-point attachment. The proteins were anchored to the surfaces using thiol groups, and the number of cysteine residues in tropomyosin was varied to test how this variation affected the structure and arrangement of the surface-attached proteins. Molecular dynamics studies were used to interpret the observed positional distributions. Based on initial studies of protein attachment to gold post structures, two 31-nm-long tropomyosin molecules were aligned between the two sidewalls of a trench with a width of 68 nm. Because the approach presented in this study uses one of twenty natural amino acids, this method provides a convenient way to pattern biomolecules on substrates using standard chemistry.
△ Less
Submitted 8 August, 2015;
originally announced August 2015.
-
Bacterial chemotaxis without gradient-sensing
Authors:
Changwook Yoon,
Yong-Jung Kim
Abstract:
Models for chemotaxis are based on gradient sensing of individual organisms. The key contribution of Keller and Segel is showing that erratic movements of individuals may result in an accurate chemotaxis phenomenon as a group. In this paper we provide another option to understand chemotactic behavior when individuals do not sense the gradient of chemical concentration by any means. We show that, i…
▽ More
Models for chemotaxis are based on gradient sensing of individual organisms. The key contribution of Keller and Segel is showing that erratic movements of individuals may result in an accurate chemotaxis phenomenon as a group. In this paper we provide another option to understand chemotactic behavior when individuals do not sense the gradient of chemical concentration by any means. We show that, if individuals increase their motility to find food when they are hungry, an accurate chemotactic behavior may obtained without sensing the gradient. Such a random dispersal has been suggested by Cho and Kim and is called starvation driven diffusion. This model is surprisingly similar to the original derivation of Keller-Segel model. A comprehensive picture of traveling band and front solutions is provided with numerical simulations.
△ Less
Submitted 29 July, 2013; v1 submitted 23 July, 2013;
originally announced July 2013.
-
Development of Hydrogen Bonding Magnetic Reaction-based Gene Regulation through Cyclic Electromagnetic DNA Simulation in Double-Stranded DNA
Authors:
Yeon Sook Kim,
Dae Gwan Lee,
Suk Keun Lee
Abstract:
The proton-magnetic reaction is commonly used in MRI machines with a strong magnetic field of over 1 T, while this study hypothesized that the electron magnetic reaction of hydrogen could affect the hydrogen bonds of double-stranded DNA (dsDNA) at a low magnetic field below 0.01 T. The goal is to develop a hydrogen bonding magnetic reaction-based gene regulation (HBMR-GR) system. The polarities of…
▽ More
The proton-magnetic reaction is commonly used in MRI machines with a strong magnetic field of over 1 T, while this study hypothesized that the electron magnetic reaction of hydrogen could affect the hydrogen bonds of double-stranded DNA (dsDNA) at a low magnetic field below 0.01 T. The goal is to develop a hydrogen bonding magnetic reaction-based gene regulation (HBMR-GR) system. The polarities of DNA base pairs are derived from the relative electrostatic charge between purines and pyrimidines, which become positively and negatively charged, respectively. The Pyu dsDNAs with pyrimidine(s)-purine(s) sequences, ds3T3A, ds3C3G, and ds3C3A, showed stronger DNA hybridization potential, increased infrared absorption at 3400-3200 cm-1, and a unique DNA conformation in HPLC analysis compared to the corresponding Puy dsDNAs. To target the three-dimensional structure of dsDNA based on the DNA base pair polarities, one can use cyclic electromagnetic DNA simulation (CEDS) with approximately 25% efficiency for randomly oriented dsDNAs. CEDS was found to induce sequence-specific hybridization of target oligo-dsDNAs in 0.005M NaCl solution and sequence-specific conformation of oligo-dsDNAs in 0.1M NaCl solution. It was found that the Pyu oligo-dsDNAs were more responsible for the hybridization and conformational changes by CEDS than the Puy oligo-dsDNAs. CEDS decreased ethidium bromide (EtBr) DNA intercalation and spermidine DNA condensation depending on CEDS time in the binding assay. The results also included that the Pyu oligo-dsDNAs were more responsible for CEDS by forming stable and unique conformation of oligo-dsDNA than the Therefore, it is postulated that the low-level HBMR-based CEDS can enhance the hybridization potential of oligo-dsDNAs and subsequently lead to the unique DNA conformation required for the initiation of various DNA functions.
△ Less
Submitted 29 June, 2024; v1 submitted 26 October, 2012;
originally announced October 2012.
-
Crowding induced entropy-enthalpy compensation in protein association equilibria
Authors:
Young C. Kim,
Jeetain Mittal
Abstract:
A statistical mechanical theory is presented to predict the effects of macromolecular crowding on protein association equilibria, accounting for both excluded volume and attractive interactions between proteins and crowding molecules. Predicted binding free energies are in excellent agreement with simulation data over a wide range of crowder sizes and packing fraction. It is shown that attractive…
▽ More
A statistical mechanical theory is presented to predict the effects of macromolecular crowding on protein association equilibria, accounting for both excluded volume and attractive interactions between proteins and crowding molecules. Predicted binding free energies are in excellent agreement with simulation data over a wide range of crowder sizes and packing fraction. It is shown that attractive interactions between proteins and crowding agents counteract the stabilizing effects of excluded volume interactions. A critical attraction strength, for which there is no net effect of crowding, is almost independent of the crowder packing fraction.
△ Less
Submitted 11 April, 2013; v1 submitted 27 September, 2012;
originally announced September 2012.
-
Innovative in silico approaches to address avian flu using grid technology
Authors:
V. Vincent Breton,
A. L. Da Costa,
P. De Vlieger,
L. Maigne,
D. Sarramia,
Y. -M. Kim,
D. Kim,
H. Q. Nguyen,
T. Solomonides,
Y. -T. Wu,
T. N. Hai
Abstract:
The recent years have seen the emergence of diseases which have spread very quickly all around the world either through human travels like SARS or animal migration like avian flu. Among the biggest challenges raised by infectious emerging diseases, one is related to the constant mutation of the viruses which turns them into continuously moving targets for drug and vaccine discovery. Another chal…
▽ More
The recent years have seen the emergence of diseases which have spread very quickly all around the world either through human travels like SARS or animal migration like avian flu. Among the biggest challenges raised by infectious emerging diseases, one is related to the constant mutation of the viruses which turns them into continuously moving targets for drug and vaccine discovery. Another challenge is related to the early detection and surveillance of the diseases as new cases can appear just anywhere due to the globalization of exchanges and the circulation of people and animals around the earth, as recently demonstrated by the avian flu epidemics. For 3 years now, a collaboration of teams in Europe and Asia has been exploring some innovative in silico approaches to better tackle avian flu taking advantage of the very large computing resources available on international grid infrastructures. Grids were used to study the impact of mutations on the effectiveness of existing drugs against H5N1 and to find potentially new leads active on mutated strains. Grids allow also the integration of distributed data in a completely secured way. The paper presents how we are currently exploring how to integrate the existing data sources towards a global surveillance network for molecular epidemiology.
△ Less
Submitted 23 December, 2008;
originally announced December 2008.
-
Metabolite essentiality elucidates robustness of Escherichia coli metabolism
Authors:
Pan-Jun Kim,
Dong-Yup Lee,
Tae Yong Kim,
Kwang Ho Lee,
Hawoong Jeong,
Sang Yup Lee,
Sunwon Park
Abstract:
Complex biological systems are very robust to genetic and environmental changes at all levels of organization. Many biological functions of Escherichia coli metabolism can be sustained against single-gene or even multiple-gene mutations by using redundant or alternative pathways. Thus, only a limited number of genes have been identified to be lethal to the cell. In this regard, the reaction-cent…
▽ More
Complex biological systems are very robust to genetic and environmental changes at all levels of organization. Many biological functions of Escherichia coli metabolism can be sustained against single-gene or even multiple-gene mutations by using redundant or alternative pathways. Thus, only a limited number of genes have been identified to be lethal to the cell. In this regard, the reaction-centric gene deletion study has a limitation in understanding the metabolic robustness. Here, we report the use of flux-sum, which is the summation of all incoming or outgoing fluxes around a particular metabolite under pseudo-steady state conditions, as a good conserved property for elucidating such robustness of E. coli from the metabolite point of view. The functional behavior, as well as the structural and evolutionary properties of metabolites essential to the cell survival, was investigated by means of a constraints-based flux analysis under perturbed conditions. The essential metabolites are capable of maintaining a steady flux-sum even against severe perturbation by actively redistributing the relevant fluxes. Disrupting the flux-sum maintenance was found to suppress cell growth. This approach of analyzing metabolite essentiality provides insight into cellular robustness and concomitant fragility, which can be used for several applications, including the development of new drugs for treating pathogens.
△ Less
Submitted 14 August, 2007;
originally announced August 2007.
-
Analyzing the Dynamic Relationship between Tumor Growth and Angiogenesis in a Two Dimensional Finite Element Model
Authors:
Eun Bo Shim,
Yoo Seok Kim,
Thomas S. Deisboeck
Abstract:
Employing a novel two-dimensional computational model we have simulated the feedback between angiogenesis and tumor growth dynamics. Analyzing vessel formation and elongation towards the concentration gradient of the tumor-derived angiogenetic basic fibroblast growth factor, bFGF, we assumed that prior to the blood vessels reaching the tumor surface, the resulting pattern of tumor growth is symm…
▽ More
Employing a novel two-dimensional computational model we have simulated the feedback between angiogenesis and tumor growth dynamics. Analyzing vessel formation and elongation towards the concentration gradient of the tumor-derived angiogenetic basic fibroblast growth factor, bFGF, we assumed that prior to the blood vessels reaching the tumor surface, the resulting pattern of tumor growth is symmetric, circular with a common center point. However, after the vessels reach the tumor surface, we assumed that the growth rate of that particular cancer region is accelerated compared to the tumor surface section that lacks neo-vascularization. Therefore, the resulting asymmetric tumor growth pattern is biased towards the site of the nourishing vessels. The simulation results show over time an increase in vessel density, a decrease in vessel branching length, and an increase in fracticality of the vascular branching architecture. Interestingly, over time the fractal dimension displayed a sigmoidal pattern with a reduced rate increase at earlier and later tumor growth stages due to distinct characteristics in vessel length and density. The finding that, at later stages, higher vascular fracticality resulted in a marked increase of tumor slice volume provides further in silico evidence for a functional impact of vascular patterns on cancer growth.
△ Less
Submitted 6 March, 2007;
originally announced March 2007.
-
Phylogenetic tree constructing algorithms fit for grid computing with SVD
Authors:
Young Rock Kim,
Oh-In Kwon,
Seong-Hun Paeng,
Chun-Jae Park
Abstract:
Erikkson showed that singular value decomposition(SVD) of flattenings determined a partition of a phylogenetic tree to be a split. In this paper, based on his work, we develop new statistically consistent algorithms fit for grid computing to construct a phylogenetic tree by computing SVD of flattenings with the small fixed number of rows.
Erikkson showed that singular value decomposition(SVD) of flattenings determined a partition of a phylogenetic tree to be a split. In this paper, based on his work, we develop new statistically consistent algorithms fit for grid computing to construct a phylogenetic tree by computing SVD of flattenings with the small fixed number of rows.
△ Less
Submitted 4 November, 2006;
originally announced November 2006.
-
Quartet consistency count method for reconstructing phylogenetic trees
Authors:
Jin-Hwan Cho,
Dosang Joe,
Young Rock Kim
Abstract:
Among the distance based algorithms in phylogenetic tree reconstruction, the neighbor-joining algorithm has been a widely used and effective method. We propose a new algorithm which counts the number of consistent quartets for cherry picking with tie breaking. We show that the success rate of the new algorithm is almost equal to that of neighbor-joining. This gives an explanation of the qualitat…
▽ More
Among the distance based algorithms in phylogenetic tree reconstruction, the neighbor-joining algorithm has been a widely used and effective method. We propose a new algorithm which counts the number of consistent quartets for cherry picking with tie breaking. We show that the success rate of the new algorithm is almost equal to that of neighbor-joining. This gives an explanation of the qualitative nature of neighbor-joining and that of dissimilarity maps from DNA sequence data. Moreover, the new algorithm always reconstructs correct trees from quartet consistent dissimilarity maps.
△ Less
Submitted 11 October, 2006;
originally announced October 2006.