-
Large language models surpass human experts in predicting neuroscience results
Authors:
Xiaoliang Luo,
Akilles Rechardt,
Guangzhi Sun,
Kevin K. Nejad,
Felipe Yáñez,
Bati Yilmaz,
Kangjoo Lee,
Alexandra O. Cohen,
Valentina Borghesani,
Anton Pashkov,
Daniele Marinazzo,
Jonathan Nicholas,
Alessandro Salatiello,
Ilia Sucholutsky,
Pasquale Minervini,
Sepehr Razavi,
Roberta Rocca,
Elkhan Yusifov,
Tereza Okalova,
Nianlong Gu,
Martin Ferianc,
Mikail Khona,
Kaustubh R. Patil,
Pui-Shee Lee,
Rui Mata
, et al. (14 additional authors not shown)
Abstract:
Scientific discoveries often hinge on synthesizing decades of research, a task that potentially outstrips human information processing capacities. Large language models (LLMs) offer a solution. LLMs trained on the vast scientific literature could potentially integrate noisy yet interrelated findings to forecast novel results better than human experts. To evaluate this possibility, we created Brain…
▽ More
Scientific discoveries often hinge on synthesizing decades of research, a task that potentially outstrips human information processing capacities. Large language models (LLMs) offer a solution. LLMs trained on the vast scientific literature could potentially integrate noisy yet interrelated findings to forecast novel results better than human experts. To evaluate this possibility, we created BrainBench, a forward-looking benchmark for predicting neuroscience results. We find that LLMs surpass experts in predicting experimental outcomes. BrainGPT, an LLM we tuned on the neuroscience literature, performed better yet. Like human experts, when LLMs were confident in their predictions, they were more likely to be correct, which presages a future where humans and LLMs team together to make discoveries. Our approach is not neuroscience-specific and is transferable to other knowledge-intensive endeavors.
△ Less
Submitted 21 June, 2024; v1 submitted 4 March, 2024;
originally announced March 2024.
-
Functional Linear Non-Gaussian Acyclic Model for Causal Discovery
Authors:
Tian-Le Yang,
Kuang-Yao Lee,
Kun Zhang,
Joe Suzuki
Abstract:
In causal discovery, non-Gaussianity has been used to characterize the complete configuration of a Linear Non-Gaussian Acyclic Model (LiNGAM), encompassing both the causal ordering of variables and their respective connection strengths. However, LiNGAM can only deal with the finite-dimensional case. To expand this concept, we extend the notion of variables to encompass vectors and even functions,…
▽ More
In causal discovery, non-Gaussianity has been used to characterize the complete configuration of a Linear Non-Gaussian Acyclic Model (LiNGAM), encompassing both the causal ordering of variables and their respective connection strengths. However, LiNGAM can only deal with the finite-dimensional case. To expand this concept, we extend the notion of variables to encompass vectors and even functions, leading to the Functional Linear Non-Gaussian Acyclic Model (Func-LiNGAM). Our motivation stems from the desire to identify causal relationships in brain-effective connectivity tasks involving, for example, fMRI and EEG datasets. We demonstrate why the original LiNGAM fails to handle these inherently infinite-dimensional datasets and explain the availability of functional data analysis from both empirical and theoretical perspectives. {We establish theoretical guarantees of the identifiability of the causal relationship among non-Gaussian random vectors and even random functions in infinite-dimensional Hilbert spaces.} To address the issue of sparsity in discrete time points within intrinsic infinite-dimensional functional data, we propose optimizing the coordinates of the vectors using functional principal component analysis. Experimental results on synthetic data verify the ability of the proposed framework to identify causal relationships among multivariate functions using the observed samples. For real data, we focus on analyzing the brain connectivity patterns derived from fMRI data.
△ Less
Submitted 17 January, 2024;
originally announced January 2024.
-
Russel and Rao Coefficient is a Suitable Substitute for Dice Coefficient in Studying Restriction Mapped Genetic Distances of Escherichia coli
Authors:
Zhu En Chay,
Chin How Lee,
Kun Cheng Lee,
Jack SH Oon,
Maurice HT Ling
Abstract:
Escherichia coli is one of many bacterial inhabitants found in human intestines and any adaptation as a result of mutations may affect its host. A commonly used technique employed to study these mutations is Restriction Fragment Length Polymorphism (RFLP) and is proceeded with a suitable distance coefficient to quantify genetic differences between 2 samples. Dice is considered a suitable distance…
▽ More
Escherichia coli is one of many bacterial inhabitants found in human intestines and any adaptation as a result of mutations may affect its host. A commonly used technique employed to study these mutations is Restriction Fragment Length Polymorphism (RFLP) and is proceeded with a suitable distance coefficient to quantify genetic differences between 2 samples. Dice is considered a suitable distance coefficient in RFLP analyses, while others were left unstudied in its suitability for use. Hence, this study aims to identify substitutes for Dice. Experimental data was obtained by subculturing E. coli for 72 passages in 8 different adaptation media and RFLP profiles analyzed using 20 distance coefficients. Our results suggest that Dennis, Fossum, Matching and Russel and Rao to work as well or better than Dice. Dennis, Matching and Fossum coefficients had highest discriminatory abilities but are limited by the lack of upper or lower boundaries. Russel and Rao coefficient is highly correlated with Dice coefficient (r2 = 0.998), with both higher and lower boundaries, suggesting that Russel and Rao coefficient can be used to substitute Dice coefficient in studying genetic distances in E. coli.
△ Less
Submitted 18 February, 2023;
originally announced February 2023.
-
Decoding Neural Correlation of Language-Specific Imagined Speech using EEG Signals
Authors:
Keon-Woo Lee,
Dae-Hyeok Lee,
Sung-Jin Kim,
Seong-Whan Lee
Abstract:
Speech impairments due to cerebral lesions and degenerative disorders can be devastating. For humans with severe speech deficits, imagined speech in the brain-computer interface has been a promising hope for reconstructing the neural signals of speech production. However, studies in the EEG-based imagined speech domain still have some limitations due to high variability in spatial and temporal inf…
▽ More
Speech impairments due to cerebral lesions and degenerative disorders can be devastating. For humans with severe speech deficits, imagined speech in the brain-computer interface has been a promising hope for reconstructing the neural signals of speech production. However, studies in the EEG-based imagined speech domain still have some limitations due to high variability in spatial and temporal information and low signal-to-noise ratio. In this paper, we investigated the neural signals for two groups of native speakers with two tasks with different languages, English and Chinese. Our assumption was that English, a non-tonal and phonogram-based language, would have spectral differences in neural computation compared to Chinese, a tonal and ideogram-based language. The results showed the significant difference in the relative power spectral density between English and Chinese in specific frequency band groups. Also, the spatial evaluation of Chinese native speakers in the theta band was distinctive during the imagination task. Hence, this paper would suggest the key spectral and spatial information of word imagination with specialized language while decoding the neural signals of speech.
△ Less
Submitted 15 April, 2022;
originally announced April 2022.
-
Bivalent Kinetics: Insights from Many Body Physics
Authors:
Richard E. Spinney,
Lawrence K. Lee,
Richard G. Morris
Abstract:
Bivalency confers several concentration-dependent phenomena, including avidity, competitive exchange and multi-site competitive exchange. Since these concepts are crucial for a wide variety of topics in cell and molecular biology, their extension, modification and/or re-purposing is also increasingly important for the design and construction of de-novo synthetic systems at the nanoscale. In this c…
▽ More
Bivalency confers several concentration-dependent phenomena, including avidity, competitive exchange and multi-site competitive exchange. Since these concepts are crucial for a wide variety of topics in cell and molecular biology, their extension, modification and/or re-purposing is also increasingly important for the design and construction of de-novo synthetic systems at the nanoscale. In this context, we draw upon classical techniques of statistical physics to revisit bivalency, highlighting that receptor site geometry offers a design modality independent of the chemistry of the individual binding interfaces themselves. Recasting the problem in terms of many-body coordination, we explore extended, translationally-invariant chains and lattices of receptor sites. This not only brings clarity to behaviours associated with simpler motifs, but also enables us to distil core principles for the rational design of concentration-dependent kinetics in synthetic soft-systems, which centre on the notion of geometric frustration. In doing so, we also reveal the possibility of other tunable spatio-temporal features, such as correlation lengths, mean-squared displacements and percolation-like transitions.
△ Less
Submitted 18 January, 2022;
originally announced January 2022.
-
Combating small molecule aggregation with machine learning
Authors:
Kuan Lee,
Ann Yang,
Yen-Chu Lin,
Daniel Reker,
Goncalo J. L. Bernardes,
Tiago Rodrigues
Abstract:
Biological screens are plagued by false positive hits resulting from aggregation. Thus, methods to triage small colloidally aggregating molecules (SCAMs) are in high demand. Herein, we disclose a bespoke machine-learning tool to confidently and intelligibly flag such entities. Our data demonstrate an unprecedented utility of machine learning for predicting SCAMs, achieving 80% of correct predictio…
▽ More
Biological screens are plagued by false positive hits resulting from aggregation. Thus, methods to triage small colloidally aggregating molecules (SCAMs) are in high demand. Herein, we disclose a bespoke machine-learning tool to confidently and intelligibly flag such entities. Our data demonstrate an unprecedented utility of machine learning for predicting SCAMs, achieving 80% of correct predictions in a challenging out-of-sample validation. The tool outperformed a panel of expert chemists, who correctly predicted 61 +/- 7% of the same test molecules in a Turing-like test. Further, the computational routine provided insight into molecular features governing aggregation that had remained hidden to expert intuition. Leveraging our tool, we quantify that up to 15-20% of ligands in publicly available chemogenomic databases have the high potential to aggregate at typical screening concentrations, imposing caution in systems biology and drug design programs. Our approach provides a means to augment human intuition, mitigate attrition and a pathway to accelerate future molecular medicine.
△ Less
Submitted 1 May, 2021;
originally announced May 2021.
-
Automated segmentation of choroidal layers from 3-dimensional macular optical coherence tomography scans
Authors:
Kyungmoo Lee,
Alexis K. Warren,
Michael D. Abramoff,
Andreas Wahle,
S. Scott Whitmore,
Ian C. Han,
John H. Fingert,
Todd E. Scheetz,
Robert F. Mullins,
Milan Sonka,
Elliott H. Sohn
Abstract:
Background: Changes in choroidal thickness are associated with various ocular diseases and the choroid can be imaged using spectral-domain optical coherence tomography (SDOCT) and enhanced depth imaging OCT (EDIOCT). New Method: Eighty macular SDOCT volumes from 80 patients were obtained using the Zeiss Cirrus machine. Eleven additional control subjects had two Cirrus scans done in one visit along…
▽ More
Background: Changes in choroidal thickness are associated with various ocular diseases and the choroid can be imaged using spectral-domain optical coherence tomography (SDOCT) and enhanced depth imaging OCT (EDIOCT). New Method: Eighty macular SDOCT volumes from 80 patients were obtained using the Zeiss Cirrus machine. Eleven additional control subjects had two Cirrus scans done in one visit along with EDIOCT using the Heidelberg Spectralis machine. To automatically segment choroidal layers from the OCT volumes, our graph-theoretic approach was utilized. The segmentation results were compared with reference standards from two graders, and the accuracy of automated segmentation was calculated using unsigned to signed border positioning thickness errors and Dice similarity coefficient (DSC). The repeatability and reproducibility of our choroidal thicknesses were determined by intraclass correlation coefficient (ICC), coefficient of variation (CV), and repeatability coefficient (RC). Results: The mean unsigned to signed border positioning errors for the choroidal inner and outer surfaces are 3.39plusminus1.26microns (mean plusminus SD) to minus1.52 plusminus 1.63microns and 16.09 plusminus 6.21microns to 4.73 plusminus 9.53microns, respectively. The mean unsigned to signed choroidal thickness errors are 16.54 plusminus 6.47microns to 6.25 plusminus 9.91microns, and the mean DSC is 0.949 plusminus 0.025. The ICC (95% CI), CV, RC values are 0.991 (0.977 to 0.997), 2.48%, 3.15microns for the repeatability and 0.991 (0.977 to 0.997), 2.49%, 0.53microns for the reproducibility studies, respectively. Comparison with Existing Method(s): The proposed method outperformed our previous method using choroidal vessel segmentation and inter-grader variability. Conclusions: This automated segmentation method can reliably measure choroidal thickness using different OCT platforms.
△ Less
Submitted 10 March, 2021;
originally announced March 2021.
-
The FRET-based structural dynamics challenge -- community contributions to consistent and open science practices
Authors:
Eitan Lerner,
Benjamin Ambrose,
Anders Barth,
Victoria Birkedal,
Scott C. Blanchard,
Richard Borner,
Thorben Cordes,
Timothy D. Craggs,
Taekjip Ha,
Gilad Haran,
Thorsten Hugel,
Antonino Ingargiola,
Achillefs Kapanidis,
Don C. Lamb,
Ted Laurence,
Nam ki Lee,
Edward A. Lemke,
Emmanuel Margeat,
Jens Michaelis,
Xavier Michalet,
Daniel Nettels,
Thomas-Otavio Peulen,
Benjamin Schuler,
Claus A. M. Seidel,
Hamid So-leimaninejad
, et al. (1 additional authors not shown)
Abstract:
Single-molecule Förster resonance energy transfer (smFRET) has become a mainstream technique for probing biomolecular structural dynamics. The rapid and wide adoption of the technique by an ever-increasing number of groups has generated many improvements and variations in the technique itself, in methods for sample preparation and characterization, in analysis of the data from such experiments, an…
▽ More
Single-molecule Förster resonance energy transfer (smFRET) has become a mainstream technique for probing biomolecular structural dynamics. The rapid and wide adoption of the technique by an ever-increasing number of groups has generated many improvements and variations in the technique itself, in methods for sample preparation and characterization, in analysis of the data from such experiments, and in analysis codes and algorithms. Recently, several labs that employ smFRET have joined forces to try to bring the smFRET community together in adopting a consensus on how to perform experiments and analyze results for achieving quantitative structural information. These recent efforts include multi-lab blind-tests to assess the accuracy and precision of smFRET between different labs using different procedures, the formal assembly of the FRET community and development of smFRET procedures to be considered for entries in the wwPDB. Here we delve into the different approaches and viewpoints in the field. This position paper describes the current "state-of-the field", points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and a list of resources that are openly available. To make further progress, we strongly encourage 'open science' practices. We hope that this position paper will provide a roadmap for newcomers to the field, as well as a reference for seasoned practitioners.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
Automated Measurements of Key Morphological Features of Human Embryos for IVF
Authors:
Brian D. Leahy,
Won-Dong Jang,
Helen Y. Yang,
Robbert Struyven,
Donglai Wei,
Zhe Sun,
Kylie R. Lee,
Charlotte Royston,
Liz Cam,
Yael Kalma,
Foad Azem,
Dalit Ben-Yosef,
Hanspeter Pfister,
Daniel Needleman
Abstract:
A major challenge in clinical In-Vitro Fertilization (IVF) is selecting the highest quality embryo to transfer to the patient in the hopes of achieving a pregnancy. Time-lapse microscopy provides clinicians with a wealth of information for selecting embryos. However, the resulting movies of embryos are currently analyzed manually, which is time consuming and subjective. Here, we automate feature e…
▽ More
A major challenge in clinical In-Vitro Fertilization (IVF) is selecting the highest quality embryo to transfer to the patient in the hopes of achieving a pregnancy. Time-lapse microscopy provides clinicians with a wealth of information for selecting embryos. However, the resulting movies of embryos are currently analyzed manually, which is time consuming and subjective. Here, we automate feature extraction of time-lapse microscopy of human embryos with a machine-learning pipeline of five convolutional neural networks (CNNs). Our pipeline consists of (1) semantic segmentation of the regions of the embryo, (2) regression predictions of fragment severity, (3) classification of the developmental stage, and object instance segmentation of (4) cells and (5) pronuclei. Our approach greatly speeds up the measurement of quantitative, biologically relevant features that may aid in embryo selection.
△ Less
Submitted 20 July, 2020; v1 submitted 29 May, 2020;
originally announced June 2020.
-
InteractionNet: Modeling and Explaining of Noncovalent Protein-Ligand Interactions with Noncovalent Graph Neural Network and Layer-Wise Relevance Propagation
Authors:
Hyeoncheol Cho,
Eok Kyun Lee,
Insung S. Choi
Abstract:
Expanding the scope of graph-based, deep-learning models to noncovalent protein-ligand interactions has earned increasing attention in structure-based drug design. Modeling the protein-ligand interactions with graph neural networks (GNNs) has experienced difficulties in the conversion of protein-ligand complex structures into the graph representation and left questions regarding whether the traine…
▽ More
Expanding the scope of graph-based, deep-learning models to noncovalent protein-ligand interactions has earned increasing attention in structure-based drug design. Modeling the protein-ligand interactions with graph neural networks (GNNs) has experienced difficulties in the conversion of protein-ligand complex structures into the graph representation and left questions regarding whether the trained models properly learn the appropriate noncovalent interactions. Here, we proposed a GNN architecture, denoted as InteractionNet, which learns two separated molecular graphs, being covalent and noncovalent, through distinct convolution layers. We also analyzed the InteractionNet model with an explainability technique, i.e., layer-wise relevance propagation, for examination of the chemical relevance of the model's predictions. Separation of the covalent and noncovalent convolutional steps made it possible to evaluate the contribution of each step independently and analyze the graph-building strategy for noncovalent interactions. We applied InteractionNet to the prediction of protein-ligand binding affinity and showed that our model successfully predicted the noncovalent interactions in both performance and relevance in chemical interpretation.
△ Less
Submitted 12 May, 2020;
originally announced May 2020.
-
Estimation of Infection Rate and Prediction of Initial Infected Individuals of COVID-19
Authors:
Seo Yoon Chae,
Kyoung-Eun Lee,
Hyun Min Lee,
Nam Jun,
Quang Ahn Le,
Biseko Juma Mafwele,
Tae Ho Lee,
Doo Hwan Kim,
Jae Woo Lee
Abstract:
We consider the pandemic spreading of COVID-19 for some selected countries after the outbreak of the coronavirus in Wuhan City, China. We estimated the infection rate and the initial infected individuals of COVID-19 by using the officially reported data at the early stage of the epidemic for the susceptible (S), infectable (I), quarantined (Q), and the cofirmed recovered (Rk) population model, so…
▽ More
We consider the pandemic spreading of COVID-19 for some selected countries after the outbreak of the coronavirus in Wuhan City, China. We estimated the infection rate and the initial infected individuals of COVID-19 by using the officially reported data at the early stage of the epidemic for the susceptible (S), infectable (I), quarantined (Q), and the cofirmed recovered (Rk) population model, so called SIQRk model. In the reported data we know the quarantined cases and the recovered cases. We can not know the recovered cases from the asymptomatic cases. In the SIQRk model we can estimated the model parameters and the initial infecting cases (confirmed ans asymtomatic cases) from the data fits. We obtained the infection rate in the range between 0.233 and 0.462, the basic reproduction number Ro in the range between 1.8 and 3.5, and the initial number of infected individuals in the range betwee 10 and 8409 for some selected countries. By using fitting parameters we estimated the maximum time of the infection for Germany when the government are performing the quarantine policy. The disease is undergoing to the calm state about six months after first patients were identified.
△ Less
Submitted 27 April, 2020;
originally announced April 2020.
-
Higher Criticism Tuned Regression For Weak And Sparse Signals
Authors:
Tao Jiang,
Stephanie J. London,
Mi Kyeong Lee,
Josyf C. Mychaleckyj,
Alison A. Motsinger-Reif
Abstract:
Here we propose a novel searching scheme for a tuning parameter in high-dimensional penalized regression methods to address variable selection and modeling when sample sizes are limited compared to the data dimensions. Our method is motivated by high-throughput biological data such as genome-wide association studies (GWAS) and epigenome-wide association studies (EWAS). We propose a new estimate of…
▽ More
Here we propose a novel searching scheme for a tuning parameter in high-dimensional penalized regression methods to address variable selection and modeling when sample sizes are limited compared to the data dimensions. Our method is motivated by high-throughput biological data such as genome-wide association studies (GWAS) and epigenome-wide association studies (EWAS). We propose a new estimate of the regularization parameter $λ$ in penalized regression methods based on an estimated lower bound of the proportion of false null hypotheses with confidence $(1 - α)$. The bound is estimated by applying the empirical null distribution of the higher criticism statistic, a second-level significance test constructed by dependent $p$-values using a multi-split regression and aggregation method. A tuning parameter estimate in penalized regression, $λ$, corresponds with the lower bound of the proportion of false null hypotheses. Different penalized regression methods with varied signal sparsity and strength are compared in the multi-split method setting. We demonstrate the performance of our method using both simulation experiments and the applications of real data on (1) lipid-trait genetics from the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial and (2) epigenetic analysis evaluating smoking's influence in differential methylation in the Agricultural Lung Health Study. The proposed algorithm is included in the HCTR package, available at https://cran.r-project.org/web/packages/HCTR/index.html.
△ Less
Submitted 10 February, 2020; v1 submitted 31 January, 2020;
originally announced February 2020.
-
Chunkflow: Distributed Hybrid Cloud Processing of Large 3D Images by Convolutional Nets
Authors:
Jingpeng Wu,
William M. Silversmith,
Kisuk Lee,
H. Sebastian Seung
Abstract:
It is now common to process volumetric biomedical images using 3D Convolutional Networks (ConvNets). This can be challenging for the teravoxel and even petavoxel images that are being acquired today by light or electron microscopy. Here we introduce chunkflow, a software framework for distributing ConvNet processing over local and cloud GPUs and CPUs. The image volume is divided into overlapping c…
▽ More
It is now common to process volumetric biomedical images using 3D Convolutional Networks (ConvNets). This can be challenging for the teravoxel and even petavoxel images that are being acquired today by light or electron microscopy. Here we introduce chunkflow, a software framework for distributing ConvNet processing over local and cloud GPUs and CPUs. The image volume is divided into overlapping chunks, each chunk is processed by a ConvNet, and the results are blended together to yield the output image. The frontend submits ConvNet tasks to a cloud queue. The tasks are executed by local and cloud GPUs and CPUs. Thanks to the fault-tolerant architecture of Chunkflow, cost can be greatly reduced by utilizing cheap unstable cloud instances. Chunkflow currently supports PyTorch for GPUs and PZnet for CPUs. To illustrate its usage, a large 3D brain image from serial section electron microscopy was processed by a 3D ConvNet with a U-Net style architecture. Chunkflow provides some chunk operations for general use, and the operations can be composed flexibly in a command line interface.
△ Less
Submitted 2 May, 2019; v1 submitted 23 April, 2019;
originally announced April 2019.
-
Reconstructing neuronal anatomy from whole-brain images
Authors:
James Gornet,
Kannan Umadevi Venkataraju,
Arun Narasimhan,
Nicholas Turner,
Kisuk Lee,
H. Sebastian Seung,
Pavel Osten,
Uygar Sümbül
Abstract:
Reconstructing multiple molecularly defined neurons from individual brains and across multiple brain regions can reveal organizational principles of the nervous system. However, high resolution imaging of the whole brain is a technically challenging and slow process. Recently, oblique light sheet microscopy has emerged as a rapid imaging method that can provide whole brain fluorescence microscopy…
▽ More
Reconstructing multiple molecularly defined neurons from individual brains and across multiple brain regions can reveal organizational principles of the nervous system. However, high resolution imaging of the whole brain is a technically challenging and slow process. Recently, oblique light sheet microscopy has emerged as a rapid imaging method that can provide whole brain fluorescence microscopy at a voxel size of 0.4 by 0.4 by 2.5 cubic microns. On the other hand, complex image artifacts due to whole-brain coverage produce apparent discontinuities in neuronal arbors. Here, we present connectivity-preserving methods and data augmentation strategies for supervised learning of neuroanatomy from light microscopy using neural networks. We quantify the merit of our approach by implementing an end-to-end automated tracing pipeline. Lastly, we demonstrate a scalable, distributed implementation that can reconstruct the large datasets that sub-micron whole-brain images produce.
△ Less
Submitted 17 March, 2019;
originally announced March 2019.
-
Dynamic Influence Networks for Rule-based Models
Authors:
Angus G. Forbes,
Andrew Burks,
Kristine Lee,
Xing Li,
Pierre Boutillier,
Jean Krivine,
Walter Fontana
Abstract:
We introduce the Dynamic Influence Network (DIN), a novel visual analytics technique for representing and analyzing rule-based models of protein-protein interaction networks. Rule-based modeling has proved instrumental in developing biological models that are concise, comprehensible, easily extensible, and that mitigate the combinatorial complexity of multi-state and multi-component biological mol…
▽ More
We introduce the Dynamic Influence Network (DIN), a novel visual analytics technique for representing and analyzing rule-based models of protein-protein interaction networks. Rule-based modeling has proved instrumental in developing biological models that are concise, comprehensible, easily extensible, and that mitigate the combinatorial complexity of multi-state and multi-component biological molecules. Our technique visualizes the dynamics of these rules as they evolve over time. Using the data produced by KaSim, an open source stochastic simulator of rule-based models written in the Kappa language, DINs provide a node-link diagram that represents the influence that each rule has on the other rules. That is, rather than representing individual biological components or types, we instead represent the rules about them (as nodes) and the current influence of these rules (as links). Using our interactive DIN-Viz software tool, researchers are able to query this dynamic network to find meaningful patterns about biological processes, and to identify salient aspects of complex rule-based models. To evaluate the effectiveness of our approach, we investigate a simulation of a circadian clock model that illustrates the oscillatory behavior of the KaiC protein phosphorylation cycle.
△ Less
Submitted 2 November, 2017;
originally announced November 2017.
-
Precision and accuracy of single-molecule FRET measurements - a worldwide benchmark study
Authors:
Björn Hellenkamp,
Sonja Schmid,
Olga Doroshenko,
Oleg Opanasyuk,
Ralf Kühnemuth,
Soheila Rezaei Adariani,
Anders Barth,
Victoria Birkedal,
Mark E. Bowen,
Hongtao Chen,
Thorben Cordes,
Tobias Eilert,
Carel Fijen,
Markus Götz,
Giorgos Gouridis,
Enrico Gratton,
Taekjip Ha,
Christian A. Hanke,
Andreas Hartmann,
Jelle Hendrix,
Lasse L. Hildebrandt,
Johannes Hohlbein,
Christian G. Hübner,
Eleni Kallis,
Achillefs N. Kapanidis
, et al. (28 additional authors not shown)
Abstract:
Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards ensuring both the reproducibility and accuracy of measuring FRET efficiencies are currently lacking. Here we report the results of a worldwide, comparative, blind study, in whic…
▽ More
Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards ensuring both the reproducibility and accuracy of measuring FRET efficiencies are currently lacking. Here we report the results of a worldwide, comparative, blind study, in which 20 labs determined the FRET efficiencies of several dye-labeled DNA duplexes. Using a unified and straightforward method, we show that FRET efficiencies can be obtained with a standard deviation between $Δ$E = +-0.02 and +-0.05. We further suggest an experimental and computational procedure for converting FRET efficiencies into accurate distances. We discuss potential uncertainties in the experiment and the modelling. Our extensive quantitative assessment of intensity-based smFRET measurements and correction procedures serve as an essential step towards validation of distance networks with the ultimate aim to archive reliable structural models of biomolecular systems obtained by smFRET-based hybrid methods.
△ Less
Submitted 29 December, 2017; v1 submitted 10 October, 2017;
originally announced October 2017.
-
A neuronal network model of interictal and recurrent ictal activity
Authors:
Marinho A. Lopes,
KyoungEun Lee,
Alexander V. Goltsev
Abstract:
We propose a neuronal network model which undergoes a saddle-node bifurcation on an invariant circle as the mechanism of the transition from the interictal to the ictal (seizure) state. In the vicinity of this transition, the model captures important dynamical features of both interictal and ictal states. We study the nature of interictal spikes and early warnings of the transition predicted by th…
▽ More
We propose a neuronal network model which undergoes a saddle-node bifurcation on an invariant circle as the mechanism of the transition from the interictal to the ictal (seizure) state. In the vicinity of this transition, the model captures important dynamical features of both interictal and ictal states. We study the nature of interictal spikes and early warnings of the transition predicted by this model. We further demonstrate that recurrent seizures emerge due to the interaction between two networks.
△ Less
Submitted 11 September, 2017;
originally announced September 2017.
-
Characterization of cortical motor function and imagery-related cortical activity: Potential application for prehabilitation
Authors:
Milena Korostenskaja,
Christoph Kapeller,
Ki H Lee,
Christoph Guger,
James Baumgartner,
Eduardo M. Castillo
Abstract:
To minimize functional morbidity associated with brain surgery, new preventive approaches (also referred to as "prehabilitation") by using motor-imagery-based computer interfaces (MI-BCIs) can be utilized. To achieve successful MI-BCI performance for prehabilitation purposes, the characteristics of an electrocorticographic (ECoG) signal that is associated with overt motor function ("real movement"…
▽ More
To minimize functional morbidity associated with brain surgery, new preventive approaches (also referred to as "prehabilitation") by using motor-imagery-based computer interfaces (MI-BCIs) can be utilized. To achieve successful MI-BCI performance for prehabilitation purposes, the characteristics of an electrocorticographic (ECoG) signal that is associated with overt motor function ("real movement" - RM) versus covert motor function ("motor imagery" - MI) need to be determined. In our current study, 5 patients with pharmacoresistant epilepsy (2 males, average age 25 years, SD 15), undergoing evaluation for epilepsy surgery participated in both RM and MI tasks. Although the RM- and MI- related ECoG changes had some common features, they also differed in a number of ways, such as location, frequency ranges, signal synchronization and desynchronization. These similarities and differences are discussed in a view of other neuroimaging studies, including magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI). We emphasize the need for inclusion of a broad spectrum of frequencies in ECoG analysis, when RM- and MI- related activities are concerned.
△ Less
Submitted 1 August, 2017;
originally announced August 2017.
-
Automatic Response Assessment in Regions of Language Cortex in Epilepsy Patients Using ECoG-based Functional Mapping and Machine Learning
Authors:
Harish RaviPrakash,
Milena Korostenskaja,
Eduardo Castillo,
Ki Lee,
James Baumgartner,
Ulas Bagci
Abstract:
Accurate localization of brain regions responsible for language and cognitive functions in Epilepsy patients should be carefully determined prior to surgery. Electrocorticography (ECoG)-based Real Time Functional Mapping (RTFM) has been shown to be a safer alternative to the electrical cortical stimulation mapping (ESM), which is currently the clinical/gold standard. Conventional methods for analy…
▽ More
Accurate localization of brain regions responsible for language and cognitive functions in Epilepsy patients should be carefully determined prior to surgery. Electrocorticography (ECoG)-based Real Time Functional Mapping (RTFM) has been shown to be a safer alternative to the electrical cortical stimulation mapping (ESM), which is currently the clinical/gold standard. Conventional methods for analyzing RTFM signals are based on statistical comparison of signal power at certain frequency bands. Compared to gold standard (ESM), they have limited accuracies when assessing channel responses.
In this study, we address the accuracy limitation of the current RTFM signal estimation methods by analyzing the full frequency spectrum of the signal and replacing signal power estimation methods with machine learning algorithms, specifically random forest (RF), as a proof of concept. We train RF with power spectral density of the time-series RTFM signal in supervised learning framework where ground truth labels are obtained from the ESM. Results obtained from RTFM of six adult patients in a strictly controlled experimental setup reveal the state of the art detection accuracy of $\approx 78\%$ for the language comprehension task, an improvement of $23\%$ over the conventional RTFM estimation method. To the best of our knowledge, this is the first study exploring the use of machine learning approaches for determining RTFM signal characteristics, and using the whole-frequency band for better region localization. Our results demonstrate the feasibility of machine learning based RTFM signal analysis method over the full spectrum to be a clinical routine in the near future.
△ Less
Submitted 6 August, 2017; v1 submitted 26 May, 2017;
originally announced June 2017.
-
Universal geometric constraints during epithelial jamming
Authors:
Lior Atia,
Dapeng Bi,
Yasha Sharma,
Jennifer A. Mitchel,
Bomi Gweon,
Stephan Koehler,
Stephen J. DeCamp,
Bo Lan,
Rebecca Hirsch,
Adrian F. Pegoraro,
Kyu Ha Lee,
Jacqueline Starr,
David A. Weitz,
Adam C. Martin,
Jin-Ah Park,
James P. Butler,
Jeffrey J. Fredberg
Abstract:
As an injury heals, an embryo develops, or a carcinoma spreads, epithelial cells systematically change their shape. In each of these processes cell shape is studied extensively, whereas variation of shape from cell-to-cell is dismissed most often as biological noise. But where do cell shape and variation of cell shape come from? Here we report that cell shape and shape variation are mutually const…
▽ More
As an injury heals, an embryo develops, or a carcinoma spreads, epithelial cells systematically change their shape. In each of these processes cell shape is studied extensively, whereas variation of shape from cell-to-cell is dismissed most often as biological noise. But where do cell shape and variation of cell shape come from? Here we report that cell shape and shape variation are mutually constrained through a relationship that is purely geometrical. That relationship is shown to govern maturation of the pseudostratified bronchial epithelial layer cultured from both non-asthmatic and asthmatic donors as well as formation of the ventral furrow in the epithelial monolayer of the Drosophila embryo in vivo. Across these and other vastly different epithelial systems, cell shape variation collapses to a family of distributions that is common to all and potentially universal. That distribution, in turn, is accounted for quantitatively by a mechanistic theory of cell-cell interaction showing that cell shape becomes progressively less elongated and less variable as the layer becomes progressively more jammed. These findings thus uncover a connection between jamming and geometry that is generic -spanning jammed living and inert systems alike- and demonstrate that proximity of the cell layer to the jammed state is the principal determinant of the most primitive features of epithelial cell shape and shape variation.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Evaporation-driven convective flows in suspensions of non-motile bacteria
Authors:
Jocelyn Dunstan,
Kyoung J. Lee,
Simon F. Park,
Yongyun Hwang,
Raymond E. Goldstein
Abstract:
We report a novel form of convection in suspensions of the bioluminiscent marine bacterium $Photobacterium~phosphoreum$. Suspensions of these bacteria placed in a chamber open to the air create persistent luminiscent plumes most easily visible when observed in the dark. These flows are strikingly similar to the classical bioconvection pattern of aerotactic swimming bacteria, which create an unstab…
▽ More
We report a novel form of convection in suspensions of the bioluminiscent marine bacterium $Photobacterium~phosphoreum$. Suspensions of these bacteria placed in a chamber open to the air create persistent luminiscent plumes most easily visible when observed in the dark. These flows are strikingly similar to the classical bioconvection pattern of aerotactic swimming bacteria, which create an unstable stratification by swimming upwards to an air-water interface, but they are a puzzle since the strain of $P.~phosphoreum$ used does not express flagella and therefore cannot swim. Systematic experimentation with suspensions of microspheres reveals that these flow patterns are driven not by the bacteria but by the accumulation of salt at the air-water interface due to evaporation of the culture medium; even at room temperature and humidity, and physiologically relevant salt concentrations, the rate of water evaporation is sufficient to drive convection patterns. A mathematical model is developed to understand the mechanism of plume formation, and linear stability analysis as well as numerical simulations were carried out to support the conclusions. While evaporation-driven convection has not been discussed extensively in the context of biological systems, these results suggest that the phenomenon may be relevant in other systems, particularly those using microorganisms of limited motility.
△ Less
Submitted 8 April, 2017;
originally announced April 2017.
-
A Fidelity-embedded Regularization Method for Robust Electrical Impedance Tomography
Authors:
Kyounghun Lee,
Eung Je Woo,
Jin Keun Seo
Abstract:
Electrical impedance tomography (EIT) provides functional images of an electrical conductivity distribution inside the human body. Since the 1980s, many potential clinical applications have arisen using inexpensive portable EIT devices. EIT acquires multiple trans-impedance measurements across the body from an array of surface electrodes around a chosen imaging slice. The conductivity image recons…
▽ More
Electrical impedance tomography (EIT) provides functional images of an electrical conductivity distribution inside the human body. Since the 1980s, many potential clinical applications have arisen using inexpensive portable EIT devices. EIT acquires multiple trans-impedance measurements across the body from an array of surface electrodes around a chosen imaging slice. The conductivity image reconstruction from the measured data is a fundamentally ill-posed inverse problem notoriously vulnerable to measurement noise and artifacts. Most available methods invert the ill-conditioned sensitivity or Jacobian matrix using a regularized least-squares data-fitting technique. Their performances rely on the regularization parameter, which controls the trade-off between fidelity and robustness. For clinical applications of EIT, it would be desirable to develop a method achieving consistent performance over various uncertain data, regardless of the choice of the regularization parameter. Based on the analysis of the structure of the Jacobian matrix, we propose a fidelity-embedded regularization (FER) method and a motion artifact removal filter. Incorporating the Jacobian matrix in the regularization process, the new FER method with the motion artifact removal filter offers stable reconstructions of high-fidelity images from noisy data by taking a very large regularization parameter value. The proposed method showed practical merits in experimental studies of chest EIT imaging.
△ Less
Submitted 6 March, 2017;
originally announced March 2017.
-
Multicellular self-organization of P. aeruginosa due to interactions with secreted trails
Authors:
Anatolij Gelimson,
Kun Zhao,
Calvin K. Lee,
W. Till Kranz,
Gerard C. L. Wong,
Ramin Golestanian
Abstract:
Guided movement in response to slowly diffusing polymeric trails provides a unique mechanism for self-organization of some microorganisms. To elucidate how this signaling route leads to microcolony formation, we experimentally probe the trajectory and orientation of Pseudomonas aeruginosa that propel themselves on a surface using type IV pili motility appendages, which preferentially attach to dep…
▽ More
Guided movement in response to slowly diffusing polymeric trails provides a unique mechanism for self-organization of some microorganisms. To elucidate how this signaling route leads to microcolony formation, we experimentally probe the trajectory and orientation of Pseudomonas aeruginosa that propel themselves on a surface using type IV pili motility appendages, which preferentially attach to deposited exopolysaccharides. We construct a stochastic model by analyzing single-bacterium trajectories, and show that the resulting theoretical prediction for the many-body behavior of the bacteria is in quantitative agreement with our experimental characterization of how cells explore the surface via a power law strategy.
△ Less
Submitted 19 September, 2016; v1 submitted 22 July, 2016;
originally announced July 2016.
-
Optical characterization of red blood cells from individuals with sickle cell trait and disease in Tanzania using quantitative phase imaging
Authors:
JaeHwang Jung,
Lucas E. Matemba,
KyeoReh Lee,
Paul E. Kazyoba,
Jonghee Yoon,
Julius J. Massaga,
Kyoohyun Kim,
Dong-Jin Kim,
YongKeun Park
Abstract:
Sickle cell disease (SCD) is common across Sub-Saharan Africa. However, the investigation of SCD in this area has been significantly limited mainly due to the lack of research facilities and skilled personnel. Here, we present optical measurements of individual red blood cells (RBCs) from healthy individuals and individuals with SCD and sickle cell trait in Tanzania using the quantitative phase im…
▽ More
Sickle cell disease (SCD) is common across Sub-Saharan Africa. However, the investigation of SCD in this area has been significantly limited mainly due to the lack of research facilities and skilled personnel. Here, we present optical measurements of individual red blood cells (RBCs) from healthy individuals and individuals with SCD and sickle cell trait in Tanzania using the quantitative phase imaging technique. By employing a quantitative phase imaging unit (QPIU), an existing microscope in a clinic is transformed into a powerful quantitative phase microscope providing measurements on the morphological, biochemical, and biomechanical properties of individual cells. The present approach will open up new opportunities for cost-effective investigation and diagnosis of several diseases in low resource environments.
△ Less
Submitted 22 April, 2016;
originally announced April 2016.
-
A simple and rapid method for detecting living microorganisms in food using laser speckle decorrelation
Authors:
Jonghee Yoon,
KyeoReh Lee,
YongKeun Park
Abstract:
Measuring microorganisms in food products is a critical issue for food safety and human health. Although various approaches for detecting low-levels of microorganisms in food have been developed, they require high-cost, complex equipment, invasive procedures, and skilled technicians which limit their widespread use in the food industry. Here, we present a simple, non-destructive, non-contact, and…
▽ More
Measuring microorganisms in food products is a critical issue for food safety and human health. Although various approaches for detecting low-levels of microorganisms in food have been developed, they require high-cost, complex equipment, invasive procedures, and skilled technicians which limit their widespread use in the food industry. Here, we present a simple, non-destructive, non-contact, and rapid optical method for measuring living microorganisms in meat products using laser speckle decorrelation. By simply measuring dynamic speckle intensity patterns reflected from samples and analyzing the temporal correlation time, the presence of living microorganisms can be non-invasively detected with high sensitivity. We present proof-of-principle demonstrations for detecting E. coli and B. cereus in chicken breast tissues.
△ Less
Submitted 17 March, 2016;
originally announced March 2016.
-
Magnetic fields facilitate DNA-mediated charge transport
Authors:
Jiun Ru Wong,
Kee Jin Lee,
Jian-Jun Shu,
Fangwei Shao
Abstract:
Exaggerate radical-induced DNA damage under magnetic fields is of great concerns to medical biosafety and to bio-molecular device based upon DNA electronic conductivity. In this report, the effect of applying an external magnetic field (MF) on DNA-mediated charge transport (CT) was investigated by studying guanine oxidation by a kinetics trap (8CPG) via photoirradiation of anthraquinone (AQ) in th…
▽ More
Exaggerate radical-induced DNA damage under magnetic fields is of great concerns to medical biosafety and to bio-molecular device based upon DNA electronic conductivity. In this report, the effect of applying an external magnetic field (MF) on DNA-mediated charge transport (CT) was investigated by studying guanine oxidation by a kinetics trap (8CPG) via photoirradiation of anthraquinone (AQ) in the presence of an external MF. Positive enhancement in CT efficiencies was observed in both the proximal and distal 8CPG after applying a static MF of 300 mT. MF assisted CT has shown sensitivities to magnetic field strength, duplex structures, and the integrity of base pair stacking. MF effects on spin evolution of charge injection upon AQ irradiation and alignment of base pairs to CT-active conformation during radical propagation were proposed to be the two major factors that MF attributed to facilitate DNA-mediated CT. Herein, our results suggested that the electronic conductivity of duplex DNA can be enhanced by applying an external MF. MF effects on DNA-mediated CT may offer a new avenue for designing DNA-based electronic device, and unraveled MF effects on redox and radical relevant biological processes.
△ Less
Submitted 14 August, 2015;
originally announced August 2015.
-
Noise-enhanced nonlinear response and the role of modular structure for signal detection in neuronal networks
Authors:
M. A. Lopes,
K. -E. Lee,
A. V. Goltsev,
J. F. F. Mendes
Abstract:
We find that sensory noise delivered together with a weak periodic signal not only enhances nonlinear response of neuronal networks, but also improves the synchronization of the response to the signal. We reveal this phenomenon in neuronal networks that are in a dynamical state near a saddle-node bifurcation corresponding to appearance of sustained network oscillations. In this state, even a weak…
▽ More
We find that sensory noise delivered together with a weak periodic signal not only enhances nonlinear response of neuronal networks, but also improves the synchronization of the response to the signal. We reveal this phenomenon in neuronal networks that are in a dynamical state near a saddle-node bifurcation corresponding to appearance of sustained network oscillations. In this state, even a weak periodic signal can evoke sharp nonlinear oscillations of neuronal activity. These sharp network oscillations have a deterministic form and amplitude determined by nonlinear dynamical equations. The signal-to-noise ratio reaches a maximum at an optimum level of sensory noise, manifesting stochastic resonance (SR) at the population level. We demonstrate SR by use of simulations and numerical integration of rate equations in a cortical model with stochastic neurons. Using this model, we mimic the experiments of Gluckman et al. [B. J. Gluckman et al, Phys. Rev. Lett., v. 77, 4098 (1996)] that have given evidence of SR in mammalian brain. We also study neuronal networks in which neurons are grouped in modules and every module works in the regime of SR. We find that even a few modules can strongly enhance the reliability of signal detection in comparison with the case when a modular organization is absent.
△ Less
Submitted 26 February, 2014;
originally announced February 2014.
-
Critical phenomena and noise-induced phase transitions in neuronal networks
Authors:
K. -E. Lee,
M. A. Lopes,
J. F. F. Mendes,
A. V. Goltsev
Abstract:
We study numerically and analytically first- and second-order phase transitions in neuronal networks stimulated by shot noise (a flow of random spikes bombarding neurons). Using an exactly solvable cortical model of neuronal networks on classical random networks, we find critical phenomena accompanying the transitions and their dependence on the shot noise intensity. We show that a pattern of spon…
▽ More
We study numerically and analytically first- and second-order phase transitions in neuronal networks stimulated by shot noise (a flow of random spikes bombarding neurons). Using an exactly solvable cortical model of neuronal networks on classical random networks, we find critical phenomena accompanying the transitions and their dependence on the shot noise intensity. We show that a pattern of spontaneous neuronal activity near a critical point of a phase transition is a characteristic property that can be used to identify the bifurcation mechanism of the transition. We demonstrate that bursts and avalanches are precursors of a first-order phase transition, paroxysmal-like spikes of activity precede a second-order phase transition caused by a saddle-node bifurcation, while irregular spindle oscillations represent spontaneous activity near a second-order phase transition caused by a supercritical Hopf bifurcation. Our most interesting result is the observation of the paroxysmal-like spikes. We show that a paroxysmal-like spike is a single nonlinear event that appears instantly from a low background activity with a rapid onset, reaches a large amplitude, and ends up with an abrupt return to lower activity. These spikes are similar to single paroxysmal spikes and sharp waves observed in EEG measurements. Our analysis shows that above the saddle-node bifurcation, sustained network oscillations appear with a large amplitude but a small frequency in contrast to network oscillations near the Hopf bifurcation that have a small amplitude but a large frequency. We discuss an amazing similarity between excitability of the cortical model stimulated by shot noise and excitability of the Morris-Lecar neuron stimulated by an applied current.
△ Less
Submitted 15 October, 2013;
originally announced October 2013.
-
Comment on "Nature of the Epidemic Threshold for the Susceptible-Infected-Susceptible Dynamics in Networks" [arXiv:1305.4819]
Authors:
Hyun Keun Lee,
Pyoung-Seop Shim,
Jae Dong Noh
Abstract:
Recently, Boguñá {\it et. al.} [Phys. Rev. Lett. {\bf 111}, 068701 (2013), arXiv:1305.4819] claimed that the epidemic threshold of the susceptible-infected-susceptible (SIS) model is zero on random scale-free (SF) networks with the small-world (SW) property. They drew the conclusion by taking into account a long-range reinfection mechanism. In this Comment, we will show that such an effect is too…
▽ More
Recently, Boguñá {\it et. al.} [Phys. Rev. Lett. {\bf 111}, 068701 (2013), arXiv:1305.4819] claimed that the epidemic threshold of the susceptible-infected-susceptible (SIS) model is zero on random scale-free (SF) networks with the small-world (SW) property. They drew the conclusion by taking into account a long-range reinfection mechanism. In this Comment, we will show that such an effect is too weak to guarantee an endemic phase with a {\em finite} fraction of infected nodes. Consequently, the epidemic threshold can be finite in random SF networks.
△ Less
Submitted 5 September, 2013;
originally announced September 2013.
-
Photoreceptors for a light biotransducer: a comparative study of the electrical responses of two (type-1)-opsins
Authors:
E Alfinito,
J Pousset,
L Reggiani,
K Lee
Abstract:
The increasing interest in photoactivated proteins as natural replacement of standard inorganic materials in photocells drives to the compared analysis of bacteriorhodopsin and proteorhodopsin, two widely diffused proteins belonging to the family of \textit{type-1} opsins. These proteins share similar behaviours but exhibit relevant differences in the sequential chain of the amino acids constituti…
▽ More
The increasing interest in photoactivated proteins as natural replacement of standard inorganic materials in photocells drives to the compared analysis of bacteriorhodopsin and proteorhodopsin, two widely diffused proteins belonging to the family of \textit{type-1} opsins. These proteins share similar behaviours but exhibit relevant differences in the sequential chain of the amino acids constituting their tertiary structure. The use of an impedance network analogue to model the protein main features provides a microscopic interpretation of a set of experiments on their photoconductance properties. In particular, this model links the protein electrical responses to the tertiary structure and to the interactions among neighbouring amino acids. The same model is also used to predict the small-signal response in terms of the Nyquist plot. Interesting enough, these rhodopsins are found to behave like a wide gap semiconductor with intrinsic conductivities of the order of $10^{-7}$ S/cm.
△ Less
Submitted 25 June, 2013;
originally announced June 2013.
-
Noise-induced phase transitions in neuronal networks
Authors:
K. -E. Lee,
M. A. Lopes,
A. V. Goltsev
Abstract:
Using an exactly solvable cortical model of a neuronal network, we show that, by increasing the intensity of shot noise (flow of random spikes bombarding neurons), the network undergoes first- and second-order non-equilibrium phase transitions. We study the nature of the transitions, bursts and avalanches of neuronal activity. Saddle-node and supercritical Hopf bifurcations are the mechanisms of e…
▽ More
Using an exactly solvable cortical model of a neuronal network, we show that, by increasing the intensity of shot noise (flow of random spikes bombarding neurons), the network undergoes first- and second-order non-equilibrium phase transitions. We study the nature of the transitions, bursts and avalanches of neuronal activity. Saddle-node and supercritical Hopf bifurcations are the mechanisms of emergence of sustained network oscillations. We show that the network stimulated by shot noise behaves similar to the Morris-Lecar model of a biological neuron stimulated by an applied current.
△ Less
Submitted 11 April, 2013;
originally announced April 2013.
-
Neural networks with dynamical synapses: from mixed-mode oscillations and spindles to chaos
Authors:
K. -E. Lee,
A. V. Goltsev,
M. A. Lopes,
J. F. F. Mendes
Abstract:
Understanding of short-term synaptic depression (STSD) and other forms of synaptic plasticity is a topical problem in neuroscience. Here we study the role of STSD in the formation of complex patterns of brain rhythms. We use a cortical circuit model of neural networks composed of irregular spiking excitatory and inhibitory neurons having type 1 and 2 excitability and stochastic dynamics. In the mo…
▽ More
Understanding of short-term synaptic depression (STSD) and other forms of synaptic plasticity is a topical problem in neuroscience. Here we study the role of STSD in the formation of complex patterns of brain rhythms. We use a cortical circuit model of neural networks composed of irregular spiking excitatory and inhibitory neurons having type 1 and 2 excitability and stochastic dynamics. In the model, neurons form a sparsely connected network and their spontaneous activity is driven by random spikes representing synaptic noise. Using simulations and analytical calculations, we found that if the STSD is absent, the neural network shows either asynchronous behavior or regular network oscillations depending on the noise level. In networks with STSD, changing parameters of synaptic plasticity and the noise level, we observed transitions to complex patters of collective activity: mixed-mode and spindle oscillations, bursts of collective activity, and chaotic behaviour. Interestingly, these patterns are stable in a certain range of the parameters and separated by critical boundaries. Thus, the parameters of synaptic plasticity can play a role of control parameters or switchers between different network states. However, changes of the parameters caused by a disease may lead to dramatic impairment of ongoing neural activity. We analyze the chaotic neural activity by use of the 0-1 test for chaos (Gottwald, G. & Melbourne, I., 2004) and show that it has a collective nature.
△ Less
Submitted 24 November, 2012;
originally announced November 2012.
-
Critical and resonance phenomena in neural networks
Authors:
A. V. Goltsev,
M. A. Lopes,
K. -E. Lee,
J. F. F. Mendes
Abstract:
Brain rhythms contribute to every aspect of brain function. Here, we study critical and resonance phenomena that precede the emergence of brain rhythms. Using an analytical approach and simulations of a cortical circuit model of neural networks with stochastic neurons in the presence of noise, we show that spontaneous appearance of network oscillations occurs as a dynamical (non-equilibrium) phase…
▽ More
Brain rhythms contribute to every aspect of brain function. Here, we study critical and resonance phenomena that precede the emergence of brain rhythms. Using an analytical approach and simulations of a cortical circuit model of neural networks with stochastic neurons in the presence of noise, we show that spontaneous appearance of network oscillations occurs as a dynamical (non-equilibrium) phase transition at a critical point determined by the noise level, network structure, the balance between excitatory and inhibitory neurons, and other parameters. We find that the relaxation time of neural activity to a steady state, response to periodic stimuli at the frequency of the oscillations, amplitude of damped oscillations, and stochastic fluctuations of neural activity are dramatically increased when approaching the critical point of the transition.
△ Less
Submitted 24 November, 2012;
originally announced November 2012.
-
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.
-
Modeling the formation of in vitro filopodia
Authors:
K. -C. Lee,
A. Gopinathan,
J. M. Schwarz
Abstract:
Filopodia are bundles of actin filaments that extend out ahead of the leading edge of a crawling cell to probe its upcoming environment. {\it In vitro} experiments [D. Vignjevic {\it et al.}, J. Cell Biol. {\bf 160}, 951 (2003)] have determined the minimal ingredients required for the formation of filopodia from the dendritic-like morphology of the leading edge. We model these experiments using…
▽ More
Filopodia are bundles of actin filaments that extend out ahead of the leading edge of a crawling cell to probe its upcoming environment. {\it In vitro} experiments [D. Vignjevic {\it et al.}, J. Cell Biol. {\bf 160}, 951 (2003)] have determined the minimal ingredients required for the formation of filopodia from the dendritic-like morphology of the leading edge. We model these experiments using kinetic aggregation equations for the density of growing bundle tips. In mean field, we determine the bundle size distribution to be broad for bundle sizes smaller than a characteristic bundle size above which the distribution decays exponentially. Two-dimensional simulations incorporating both bundling and cross-linking measure a bundle size distribution that agrees qualitatively with mean field. The simulations also demonstrate a nonmonotonicity in the radial extent of the dendritic region as a function of capping protein concentration, as was observed in experiments, due to the interplay between percolation and the ratcheting of growing filaments off a spherical obstacle.
△ Less
Submitted 13 May, 2010; v1 submitted 14 September, 2009;
originally announced September 2009.
-
New Proposed Mechanism of Actin-Polymerization-Driven Motility
Authors:
Kun-Chun Lee,
Andrea J. Liu
Abstract:
We present the first numerical simulation of actin-driven propulsion by elastic filaments. Specifically, we use a Brownian dynamics formulation of the dendritic nucleation model of actin-driven propulsion. We show that the model leads to a self-assembled network that exerts forces on a disk and pushes it with an average speed. This simulation approach is the first to observe a speed that varies…
▽ More
We present the first numerical simulation of actin-driven propulsion by elastic filaments. Specifically, we use a Brownian dynamics formulation of the dendritic nucleation model of actin-driven propulsion. We show that the model leads to a self-assembled network that exerts forces on a disk and pushes it with an average speed. This simulation approach is the first to observe a speed that varies non-monotonically with the concentration of branching proteins (Arp2/3), capping protein and depolymerization rate (ADF), in accord with experimental observations. Our results suggest a new interpretation of the origin of motility that can be tested readily by experiment.
△ Less
Submitted 6 June, 2008; v1 submitted 4 June, 2008;
originally announced June 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.
-
Branching, Capping, and Severing in Dynamic Actin Structures
Authors:
Ajay Gopinathan,
Kun-Chun Lee,
J. M. Schwarz,
Andrea J. Liu
Abstract:
Branched actin networks at the leading edge of a crawling cell evolve via protein-regulated processes such as polymerization, depolymerization, capping, branching, and severing. A formulation of these processes is presented and analyzed to study steady-state network morphology. In bulk, we identify several scaling regimes in severing and branching protein concentrations and find that the couplin…
▽ More
Branched actin networks at the leading edge of a crawling cell evolve via protein-regulated processes such as polymerization, depolymerization, capping, branching, and severing. A formulation of these processes is presented and analyzed to study steady-state network morphology. In bulk, we identify several scaling regimes in severing and branching protein concentrations and find that the coupling between severing and branching is optimally exploited for conditions {\it in vivo}. Near the leading edge, we find qualitative agreement with the {\it in vivo} morphology.
△ Less
Submitted 15 March, 2007;
originally announced March 2007.
-
Emergence of Chaotic Itinerancy in Simple Ecological Systems
Authors:
Pan-Jun Kim,
Tae-Wook Ko,
Hawoong Jeong,
Kyoung J. Lee,
Seung Kee Han
Abstract:
Chaotic itinerancy is a universal dynamical concept that describes itinerant motion among many different ordered states through chaotic transition in dynamical systems. Unlike the expectation of the prevalence of chaotic itinerancy in high-dimensional systems, we identify chaotic itinerant behavior from a relatively simple ecological system, which consists only of two coupled consumer-resource p…
▽ More
Chaotic itinerancy is a universal dynamical concept that describes itinerant motion among many different ordered states through chaotic transition in dynamical systems. Unlike the expectation of the prevalence of chaotic itinerancy in high-dimensional systems, we identify chaotic itinerant behavior from a relatively simple ecological system, which consists only of two coupled consumer-resource pairs. The system exhibits chaotic bursting activity, in which the explosion and the shrinkage of the population alternate indefinitely, while the explosion of one pair co-occurs with the shrinkage of the other pair. We analyze successfully the bursting activity in the framework of chaotic itinerancy, and find that large duration times of bursts tend to cluster in time, allowing the effective burst prognosis. We also investigate the control schemes on the bursting activity, and demonstrate that invoking the competitive rise of the consumer in one pair can even elongate the burst of the other pair rather than shorten it.
△ Less
Submitted 16 December, 2007; v1 submitted 3 December, 2006;
originally announced December 2006.
-
The statistical mechanics of complex signaling networks : nerve growth factor signaling
Authors:
Kevin S. Brown,
Colin C. Hill,
Guillermo A. Calero,
Kelvin H. Lee,
James P. Sethna,
Richard A. Cerione
Abstract:
It is becoming increasingly appreciated that the signal transduction systems used by eukaryotic cells to achieve a variety of essential responses represent highly complex networks rather than simple linear pathways. While significant effort is being made to experimentally measure the rate constants for individual steps in these signaling networks, many of the parameters required to describe the…
▽ More
It is becoming increasingly appreciated that the signal transduction systems used by eukaryotic cells to achieve a variety of essential responses represent highly complex networks rather than simple linear pathways. While significant effort is being made to experimentally measure the rate constants for individual steps in these signaling networks, many of the parameters required to describe the behavior of these systems remain unknown, or at best, estimates. With these goals and caveats in mind, we use methods of statistical mechanics to extract useful predictions for complex cellular signaling networks. To establish the usefulness of our approach, we have applied our methods towards modeling the nerve growth factor (NGF)-induced differentiation of neuronal cells. Using our approach, we are able to extract predictions that are highly specific and accurate, thereby enabling us to predict the influence of specific signaling modules in determining the integrated cellular response to the two growth factors. We show that extracting biologically relevant predictions from complex signaling models appears to be possible even in the absence of measurements of all the individual rate constants. Our methods also raise some interesting insights into the design and possible evolution of cellular systems, highlighting an inherent property of these systems wherein particular ''soft'' combinations of parameters can be varied over wide ranges without impacting the final output and demonstrating that a few ''stiff'' parameter combinations center around the paramount regulatory steps of the network. We refer to this property -- which is distinct from robustness -- as ''sloppiness.''
△ Less
Submitted 22 June, 2004;
originally announced June 2004.
-
Mean field approaches to the totally asymmetric exclusion process with quenched disorder and large particles
Authors:
Leah B. Shaw,
James P. Sethna,
Kelvin H. Lee
Abstract:
The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles (ribosomes) have spatial extent, covering more than one lattice site. Realistic, nonuniform gene sequences lead to quenched disorder in the particle hopping rates. We study the totally asymmetric exclusion process with large particles and quenched disorder via several mean…
▽ More
The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles (ribosomes) have spatial extent, covering more than one lattice site. Realistic, nonuniform gene sequences lead to quenched disorder in the particle hopping rates. We study the totally asymmetric exclusion process with large particles and quenched disorder via several mean field approaches and compare the mean field results with Monte Carlo simulations. Mean field equations obtained from the literature are found to be reasonably effective in describing this system. A numerical technique is developed for computing the particle current rapidly. The mean field approach is extended to include two-point correlations between adjacent sites. The two-point results are found to match Monte Carlo simulations more closely.
△ Less
Submitted 19 March, 2004;
originally announced March 2004.
-
The totally asymmetric exclusion process with extended objects, a model for protein synthesis
Authors:
Leah B. Shaw,
R. K. P. Zia,
Kelvin H. Lee
Abstract:
The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles have spatial extent, covering more than one lattice site. We expand the well studied Totally Asymmetric Exclusion Process (TASEP), in which particles typically cover a single lattice site, to include cases with extended objects. Exact solutions can be determined for a unif…
▽ More
The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles have spatial extent, covering more than one lattice site. We expand the well studied Totally Asymmetric Exclusion Process (TASEP), in which particles typically cover a single lattice site, to include cases with extended objects. Exact solutions can be determined for a uniform closed system. We analyze the uniform open system through two approaches. First, a continuum limit produces a modified diffusion equation for particle density profiles. Second, an extremal principle based on domain wall theory accurately predicts the phase diagram and currents in each phase. Finally, we briefly consider approximate approaches to a non-uniform open system with quenched disorder in the particle hopping rates and compare these approaches with Monte Carlo simulations.
△ Less
Submitted 16 September, 2003; v1 submitted 6 February, 2003;
originally announced February 2003.
-
Competing Patterns of Signaling Activity in Dictyostelium discoideum
Authors:
Kyoung J. Lee,
Edward C. Cox,
Raymond E. Goldstein
Abstract:
Quantitative experiments are described on spatio-temporal patterns of coherent chemical signaling activity in populations of {\it Dictyostelium discoideum} amoebae. We observe competition between spontaneously firing centers and rotating spiral waves that depends strongly on the overall cell density. At low densities, no complete spirals appear and chemotactic aggregation is driven by periodic c…
▽ More
Quantitative experiments are described on spatio-temporal patterns of coherent chemical signaling activity in populations of {\it Dictyostelium discoideum} amoebae. We observe competition between spontaneously firing centers and rotating spiral waves that depends strongly on the overall cell density. At low densities, no complete spirals appear and chemotactic aggregation is driven by periodic concentric waves, whereas at high densities the firing centers seen at early times nucleate and are apparently entrained by spiral waves whose cores ultimately serve as aggregation centers. Possible mechanisms for these observations are discussed.
△ Less
Submitted 28 January, 1996;
originally announced January 1996.