-
Manifold Transform by Recurrent Cortical Circuit Enhances Robust Encoding of Familiar Stimuli
Authors:
Weifan Wang,
Xueyan Niu,
Tai-Sing Lee
Abstract:
A ubiquitous phenomenon observed throughout the primate hierarchical visual system is the sparsification of the neural representation of visual stimuli as a result of familiarization by repeated exposure, manifested as the sharpening of the population tuning curves and suppression of neural responses at the population level. In this work, we investigated the computational implications and circuit…
▽ More
A ubiquitous phenomenon observed throughout the primate hierarchical visual system is the sparsification of the neural representation of visual stimuli as a result of familiarization by repeated exposure, manifested as the sharpening of the population tuning curves and suppression of neural responses at the population level. In this work, we investigated the computational implications and circuit mechanisms underlying these neurophysiological observations in an early visual cortical circuit model. We found that such a recurrent neural circuit, shaped by BCM Hebbian learning, can also reproduce these phenomena. The resulting circuit became more robust against noises in encoding the familiar stimuli. Analysis of the geometry of the neural response manifold revealed that recurrent computation and familiar learning transform the response manifold and the neural dynamics, resulting in enhanced robustness against noise and better stimulus discrimination. This prediction is supported by preliminary physiological evidence. Familiarity training increases the alignment of the slow modes of network dynamics with the invariant features of the learned images. These findings revealed how these rapid plasticity mechanisms can improve contextual visual processing in even the early visual areas in the hierarchical visual system.
△ Less
Submitted 20 August, 2024;
originally announced August 2024.
-
Incremental Learning and Self-Attention Mechanisms Improve Neural System Identification
Authors:
Isaac Lin,
Tianye Wang,
Shang Gao,
Shiming Tang,
Tai Sing Lee
Abstract:
Convolutional neural networks (CNNs) have been shown to be the state-of-the-art approach for modeling the transfer functions of visual cortical neurons. Cortical neurons in the primary visual cortex are are sensitive to contextual information mediated by extensive horizontal and feedback connections. Standard CNNs can integrate global spatial image information to model such contextual modulation v…
▽ More
Convolutional neural networks (CNNs) have been shown to be the state-of-the-art approach for modeling the transfer functions of visual cortical neurons. Cortical neurons in the primary visual cortex are are sensitive to contextual information mediated by extensive horizontal and feedback connections. Standard CNNs can integrate global spatial image information to model such contextual modulation via two mechanisms: successive rounds of convolutions and a fully connected readout layer. In this paper, we find that non-local networks or self-attention (SA) mechanisms, theoretically related to context-dependent flexible gating mechanisms observed in the primary visual cortex, improve neural response predictions over parameter-matched CNNs in two key metrics: tuning curve correlation and tuning peak. We factorize networks to determine the relative contribution of each context mechanism. This reveals that information in the local receptive field is most important for modeling the overall tuning curve, but surround information is critically necessary for characterizing the tuning peak. We find that self-attention can replace subsequent spatial-integration convolutions when learned in an incremental manner, and is further enhanced in the presence of a fully connected readout layer, suggesting that the two context mechanisms are complementary. Finally, we find that learning a receptive-field-centric model with self-attention, before incrementally learning a fully connected readout, yields a more biologically realistic model in terms of center-surround contributions.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Spatial modeling algorithms for reactions and transport (SMART) in biological cells
Authors:
Emmet A. Francis,
Justin G. Laughlin,
Jørgen S. Dokken,
Henrik N. T. Finsberg,
Christopher T. Lee,
Marie E. Rognes,
Padmini Rangamani
Abstract:
Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their many functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Herein, we present our software, Spatial Modeling Algorithms for Reactions and Transport (SMART), a package that takes…
▽ More
Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their many functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Herein, we present our software, Spatial Modeling Algorithms for Reactions and Transport (SMART), a package that takes in high-level user specifications about cell signaling networks and molecular transport, and then assembles and solves the associated mathematical and computational systems. SMART uses state-of-the-art finite element analysis, via the FEniCS Project software, to efficiently and accurately resolve cell signaling events over discretized cellular and subcellular geometries. We demonstrate its application to several different biological systems, including YAP/TAZ mechanotransduction, calcium signaling in neurons and cardiomyocytes, and ATP generation in mitochondria. Throughout, we utilize experimentally-derived realistic cellular geometries represented by well-conditioned tetrahedral meshes. These scenarios demonstrate the applicability, flexibility, accuracy and efficiency of SMART across a range of temporal and spatial scales.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
A large calcium-imaging dataset reveals a systematic V4 organization for natural scenes
Authors:
Tianye Wang,
Haoxuan Yao,
Tai Sing Lee,
Jiayi Hong,
Yang Li,
Hongfei Jiang,
Ian Max Andolina,
Shiming Tang
Abstract:
The visual system evolved to process natural scenes, yet most of our understanding of the topology and function of visual cortex derives from studies using artificial stimuli. To gain deeper insights into visual processing of natural scenes, we utilized widefield calcium-imaging of primate V4 in response to many natural images, generating a large dataset of columnar-scale responses. We used this d…
▽ More
The visual system evolved to process natural scenes, yet most of our understanding of the topology and function of visual cortex derives from studies using artificial stimuli. To gain deeper insights into visual processing of natural scenes, we utilized widefield calcium-imaging of primate V4 in response to many natural images, generating a large dataset of columnar-scale responses. We used this dataset to build a digital twin of V4 via deep learning, generating a detailed topographical map of natural image preferences at each cortical position. The map revealed clustered functional domains for specific classes of natural image features. These ranged from surface-related attributes like color and texture to shape-related features such as edges, curvature, and facial features. We validated the model-predicted domains with additional widefield calcium-imaging and single-cell resolution two-photon imaging. Our study illuminates the detailed topological organization and neural codes in V4 that represent natural scenes.
△ Less
Submitted 23 July, 2023; v1 submitted 3 July, 2023;
originally announced July 2023.
-
SMART: Spatial Modeling Algorithms for Reaction and Transport
Authors:
Justin G. Laughlin,
Jørgen S. Dokken,
Henrik N. T. Finsberg,
Emmet A. Francis,
Christopher T. Lee,
Marie E. Rognes,
Padmini Rangamani
Abstract:
Recent advances in microscopy and 3D reconstruction methods have allowed for characterization of cellular morphology in unprecedented detail, including the irregular geometries of intracellular subcompartments such as membrane-bound organelles. These geometries are now compatible with predictive modeling of cellular function. Biological cells respond to stimuli through sequences of chemical reacti…
▽ More
Recent advances in microscopy and 3D reconstruction methods have allowed for characterization of cellular morphology in unprecedented detail, including the irregular geometries of intracellular subcompartments such as membrane-bound organelles. These geometries are now compatible with predictive modeling of cellular function. Biological cells respond to stimuli through sequences of chemical reactions generally referred to as cell signaling pathways. The propagation and reaction of chemical substances in cell signaling pathways can be represented by coupled nonlinear systems of reaction-transport equations. These reaction pathways include numerous chemical species that react across boundaries or interfaces (e.g., the cell membrane and membranes of organelles within the cell) and domains (e.g., the bulk cell volume and the interior of organelles). Such systems of multi-dimensional partial differential equations (PDEs) are notoriously difficult to solve because of their high dimensionality, non-linearities, strong coupling, stiffness, and potential instabilities. In this work, we describe Spatial Modeling Algorithms for Reactions and Transport (SMART), a high-performance finite-element-based simulation package for model specification and numerical simulation of spatially-varying reaction-transport processes. SMART is based on the FEniCS finite element library, provides a symbolic representation framework for specifying reaction pathways, and supports geometries in 2D and 3D including large and irregular cell geometries obtained from modern ultrastructural characterization methods.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Conceptual Framework and Documentation Standards of Cystoscopic Media Content for Artificial Intelligence
Authors:
Okyaz Eminaga,
Timothy Jiyong Lee,
Jessie Ge,
Eugene Shkolyar,
Mark Laurie,
Jin Long,
Lukas Graham Hockman,
Joseph C. Liao
Abstract:
Background: The clinical documentation of cystoscopy includes visual and textual materials. However, the secondary use of visual cystoscopic data for educational and research purposes remains limited due to inefficient data management in routine clinical practice. Methods: A conceptual framework was designed to document cystoscopy in a standardized manner with three major sections: data management…
▽ More
Background: The clinical documentation of cystoscopy includes visual and textual materials. However, the secondary use of visual cystoscopic data for educational and research purposes remains limited due to inefficient data management in routine clinical practice. Methods: A conceptual framework was designed to document cystoscopy in a standardized manner with three major sections: data management, annotation management, and utilization management. A Swiss-cheese model was proposed for quality control and root cause analyses. We defined the infrastructure required to implement the framework with respect to FAIR (findable, accessible, interoperable, re-usable) principles. We applied two scenarios exemplifying data sharing for research and educational projects to ensure the compliance with FAIR principles. Results: The framework was successfully implemented while following FAIR principles. The cystoscopy atlas produced from the framework could be presented in an educational web portal; a total of 68 full-length qualitative videos and corresponding annotation data were sharable for artificial intelligence projects covering frame classification and segmentation problems at case, lesion and frame levels. Conclusion: Our study shows that the proposed framework facilitates the storage of the visual documentation in a standardized manner and enables FAIR data for education and artificial intelligence research.
△ Less
Submitted 18 January, 2023; v1 submitted 14 January, 2023;
originally announced January 2023.
-
Acoustical Analysis of Speech Under Physical Stress in Relation to Physical Activities and Physical Literacy
Authors:
Si-Ioi Ng,
Rui-Si Ma,
Tan Lee,
Raymond Kim-Wai Sum
Abstract:
Human speech production encompasses physiological processes that naturally react to physic stress. Stress caused by physical activity (PA), e.g., running, may lead to significant changes in a person's speech. The major changes are related to the aspects of pitch level, speaking rate, pause pattern, and breathiness. The extent of change depends presumably on physical fitness and well-being of the p…
▽ More
Human speech production encompasses physiological processes that naturally react to physic stress. Stress caused by physical activity (PA), e.g., running, may lead to significant changes in a person's speech. The major changes are related to the aspects of pitch level, speaking rate, pause pattern, and breathiness. The extent of change depends presumably on physical fitness and well-being of the person, as well as intensity of PA. The general wellness of a person is further related to his/her physical literacy (PL), which refers to a holistic description of engagement in PA. This paper presents the development of a Cantonese speech database that contains audio recordings of speech before and after physical exercises of different intensity levels. The corpus design and data collection process are described. Preliminary results of acoustical analysis are presented to illustrate the impact of PA on pitch level, pitch range, speaking and articulation rate, and time duration of pauses. It is also noted that the effect of PA is correlated to some of the PA and PL measures.
△ Less
Submitted 25 April, 2022; v1 submitted 20 November, 2021;
originally announced November 2021.
-
Mem3DG: Modeling Membrane Mechanochemical Dynamics in 3D using Discrete Differential Geometry
Authors:
Cuncheng Zhu,
Christopher T. Lee,
Padmini Rangamani
Abstract:
Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared to approximation-based methods (e.g., finite element method), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimens…
▽ More
Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared to approximation-based methods (e.g., finite element method), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, $\texttt{Mem3DG}$, for modeling 3D mechanochemical dynamics of membranes based on Discrete Differential Geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, $\texttt{Mem3DG}$ is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.
△ Less
Submitted 31 October, 2021;
originally announced November 2021.
-
Exploiting Pre-Trained ASR Models for Alzheimer's Disease Recognition Through Spontaneous Speech
Authors:
Ying Qin,
Wei Liu,
Zhiyuan Peng,
Si-Ioi Ng,
Jingyu Li,
Haibo Hu,
Tan Lee
Abstract:
Alzheimer's disease (AD) is a progressive neurodegenerative disease and recently attracts extensive attention worldwide. Speech technology is considered a promising solution for the early diagnosis of AD and has been enthusiastically studied. Most recent works concentrate on the use of advanced BERT-like classifiers for AD detection. Input to these classifiers are speech transcripts produced by au…
▽ More
Alzheimer's disease (AD) is a progressive neurodegenerative disease and recently attracts extensive attention worldwide. Speech technology is considered a promising solution for the early diagnosis of AD and has been enthusiastically studied. Most recent works concentrate on the use of advanced BERT-like classifiers for AD detection. Input to these classifiers are speech transcripts produced by automatic speech recognition (ASR) models. The major challenge is that the quality of transcription could degrade significantly under complex acoustic conditions in the real world. The detection performance, in consequence, is largely limited. This paper tackles the problem via tailoring and adapting pre-trained neural-network based ASR model for the downstream AD recognition task. Only bottom layers of the ASR model are retained. A simple fully-connected neural network is added on top of the tailored ASR model for classification. The heavy BERT classifier is discarded. The resulting model is light-weight and can be fine-tuned in an end-to-end manner for AD recognition. Our proposed approach takes only raw speech as input, and no extra transcription process is required. The linguistic information of speech is implicitly encoded in the tailored ASR model and contributes to boosting the performance. Experiments show that our proposed approach outperforms the best manual transcript-based RoBERTa by an absolute margin of 4.6% in terms of accuracy. Our best-performing models achieve the accuracy of 83.2% and 78.0% in the long-audio and short-audio competition tracks of the 2021 NCMMSC Alzheimer's Disease Recognition Challenge, respectively.
△ Less
Submitted 4 October, 2021;
originally announced October 2021.
-
A study of computational and conceptual complexities of compartment and agent based models
Authors:
Prateek Kunwar,
Oleksandr Markovichenko,
Monique Chyba,
Yuriy Mileyko,
Alice Koniges,
Thomas Lee
Abstract:
The ongoing COVID-19 pandemic highlights the essential role of mathematical models in understanding the spread of the virus along with a quantifiable and science-based prediction of the impact of various mitigation measures. Numerous types of models have been employed with various levels of success. This leads to the question of what kind of a mathematical model is most appropriate for a given sit…
▽ More
The ongoing COVID-19 pandemic highlights the essential role of mathematical models in understanding the spread of the virus along with a quantifiable and science-based prediction of the impact of various mitigation measures. Numerous types of models have been employed with various levels of success. This leads to the question of what kind of a mathematical model is most appropriate for a given situation. We consider two widely used types of models: equation-based models (such as standard compartmental epidemiological models) and agent-based models. We assess their performance by modeling the spread of COVID-19 on the Hawaiian island of Oahu under different scenarios. We show that when it comes to information crucial to decision making, both models produce very similar results. At the same time, the two types of models exhibit very different characteristics when considering their computational and conceptual complexity. Consequently, we conclude that choosing the model should be mostly guided by available computational and human resources.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
Ubiquitous proximity to a critical state for collective neural activity in the CA1 region of freely moving mice
Authors:
Yi-Ling Chen,
Chun-Chung Chen,
Yu-Ying Mei,
Ning Zhou,
Dongchuan Wu,
Ting-Kuo Lee
Abstract:
Using miniscope recordings of calcium fluorescence signals in the CA1 region of the hippocampus of mice, we monitor the neural activity of hippocampal regions while the animals are freely moving in an open chamber. Using a data-driven statistical modeling approach, the statistical properties of the recorded data are mapped to spin-glass models with pairwise interactions. Considering the parameter…
▽ More
Using miniscope recordings of calcium fluorescence signals in the CA1 region of the hippocampus of mice, we monitor the neural activity of hippocampal regions while the animals are freely moving in an open chamber. Using a data-driven statistical modeling approach, the statistical properties of the recorded data are mapped to spin-glass models with pairwise interactions. Considering the parameter space of the model, the observed system is generally near a critical state between two distinct phases. The close proximity to the criticality is found to be robust against different ways of sampling and segmentation of the measured data. By independently altering the coupling distribution and the network structure of the statistical model, the network structures are found to be vital to maintain the proximity to the critical state. We further find the observed assignment of the coupling strengths makes the net coupling at each site more balanced with slight variation, which likely helps the maintenance of the critical state. Network analysis on the connectivity obtained by thresholding the coupling strengths find the connectivity of the networks to be well described by a random network model. These results are consistent across different experiments, sampling and segmentation choices in our analysis.
△ Less
Submitted 25 February, 2021;
originally announced February 2021.
-
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.
-
Estimating Fiber Orientation Distribution through Blockwise Adaptive Thresholding with Application to HCP Young Adults Data
Authors:
Seungyong Hwang,
Thomas C. M. Lee,
Debashis Paul,
Jie Peng
Abstract:
Due to recent technological advances, large brain imaging data sets can now be collected. Such data are highly complex so extraction of meaningful information from them remains challenging. Thus, there is an urgent need for statistical procedures that are computationally scalable and can provide accurate estimates that capture the neuronal structures and their functionalities. We propose a fast me…
▽ More
Due to recent technological advances, large brain imaging data sets can now be collected. Such data are highly complex so extraction of meaningful information from them remains challenging. Thus, there is an urgent need for statistical procedures that are computationally scalable and can provide accurate estimates that capture the neuronal structures and their functionalities. We propose a fast method for estimating the fiber orientation distribution(FOD) based on diffusion MRI data. This method models the observed dMRI signal at any voxel as a convolved and noisy version of the underlying FOD, and utilizes the spherical harmonics basis for representing the FOD, where the spherical harmonic coefficients are adaptively and nonlinearly shrunk by using a James-Stein type estimator. To further improve the estimation accuracy by enhancing the localized peaks of the FOD, as a second step a super-resolution sharpening process is then applied. The resulting estimated FODs can be fed to a fiber tracking algorithm to reconstruct the white matter fiber tracts. We illustrate the overall methodology using both synthetic data and data from the Human Connectome Project.
△ Less
Submitted 28 June, 2021; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Recurrent Feedback Improves Feedforward Representations in Deep Neural Networks
Authors:
Siming Yan,
Xuyang Fang,
Bowen Xiao,
Harold Rockwell,
Yimeng Zhang,
Tai Sing Lee
Abstract:
The abundant recurrent horizontal and feedback connections in the primate visual cortex are thought to play an important role in bringing global and semantic contextual information to early visual areas during perceptual inference, helping to resolve local ambiguity and fill in missing details. In this study, we find that introducing feedback loops and horizontal recurrent connections to a deep co…
▽ More
The abundant recurrent horizontal and feedback connections in the primate visual cortex are thought to play an important role in bringing global and semantic contextual information to early visual areas during perceptual inference, helping to resolve local ambiguity and fill in missing details. In this study, we find that introducing feedback loops and horizontal recurrent connections to a deep convolution neural network (VGG16) allows the network to become more robust against noise and occlusion during inference, even in the initial feedforward pass. This suggests that recurrent feedback and contextual modulation transform the feedforward representations of the network in a meaningful and interesting way. We study the population codes of neurons in the network, before and after learning with feedback, and find that learning with feedback yielded an increase in discriminability (measured by d-prime) between the different object classes in the population codes of the neurons in the feedforward path, even at the earliest layer that receives feedback. We find that recurrent feedback, by injecting top-down semantic meaning to the population activities, helps the network learn better feedforward paths to robustly map noisy image patches to the latent representations corresponding to important visual concepts of each object class, resulting in greater robustness of the network against noises and occlusion as well as better fine-grained recognition.
△ Less
Submitted 22 December, 2019;
originally announced December 2019.
-
Complexity and Diversity in Sparse Code Priors Improve Receptive Field Characterization of Macaque V1 Neurons
Authors:
Ziniu Wu,
Harold Rockwell,
Yimeng Zhang,
Shiming Tang,
Tai Sing Lee
Abstract:
System identification techniques -- projection pursuit regression models (PPRs) and convolutional neural networks (CNNs) -- provide state-of-the-art performance in predicting visual cortical neurons' responses to arbitrary input stimuli. However, the constituent kernels recovered by these methods are often noisy and lack coherent structure, making it difficult to understand the underlying componen…
▽ More
System identification techniques -- projection pursuit regression models (PPRs) and convolutional neural networks (CNNs) -- provide state-of-the-art performance in predicting visual cortical neurons' responses to arbitrary input stimuli. However, the constituent kernels recovered by these methods are often noisy and lack coherent structure, making it difficult to understand the underlying component features of a neuron's receptive field.
In this paper, we show that using a dictionary of diverse kernels with complex shapes learned from natural scenes based on efficient coding theory, as the front-end for PPRs and CNNs can improve their performance in neuronal response prediction as well as algorithmic data efficiency and convergence speed. Extensive experimental results also indicate that these sparse-code kernels provide important information on the component features of a neuron's receptive field. In addition, we find that models with the complex-shaped sparse code front-end are significantly better than models with a standard orientation-selective Gabor filter front-end for modeling V1 neurons that have been found to exhibit complex pattern selectivity. We show that the relative performance difference due to these two front-ends can be used to produce a sensitive metric for detecting complex selectivity in V1 neurons.
△ Less
Submitted 1 October, 2021; v1 submitted 19 November, 2019;
originally announced November 2019.
-
An Open Source Mesh Generation Platform for Biophysical Modeling Using Realistic Cellular Geometries
Authors:
Christopher T. Lee,
Justin G. Laughlin,
John B. Moody,
Rommie E. Amaro,
J. Andrew McCammon,
Michael J. Holst,
Padmini Rangamani
Abstract:
Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution o…
▽ More
Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this study, we outline the steps for a naïve user to approach GAMer 2, a mesh generation code written in C++ designed to convert structural datasets to realistic geometric meshes, while preserving the underlying shapes. We present two example cases, 1) mesh generation at the subcellular scale as informed by electron tomography, and 2) meshing a protein with structure from x-ray crystallography. We further demonstrate that the meshes generated by GAMer are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.
△ Less
Submitted 16 December, 2019; v1 submitted 10 September, 2019;
originally announced September 2019.
-
3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries
Authors:
Christopher T. Lee,
Justin G. Laughlin,
Nils Angliviel de La Beaumelle,
Rommie E. Amaro,
J. Andrew McCammon,
Ravi Ramamoorthi,
Michael J. Holst,
Padmini Rangamani
Abstract:
Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer sim…
▽ More
Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.
△ Less
Submitted 17 December, 2019; v1 submitted 29 January, 2019;
originally announced January 2019.
-
An Evolutionary Game Theoretic Model of Rhino Horn Devaluation
Authors:
Nikoleta E. Glynatsi,
Vincent Knight,
Tamsin E. Lee
Abstract:
Rhino populations are at a critical level due to the demand for rhino horn and the subsequent poaching. Wildlife managers attempt to secure rhinos with approaches to devalue the horn, the most common of which is dehorning. Game theory has been used to examine the interaction of poachers and wildlife managers where a manager can either `dehorn' their rhinos or leave the horn attached and poachers m…
▽ More
Rhino populations are at a critical level due to the demand for rhino horn and the subsequent poaching. Wildlife managers attempt to secure rhinos with approaches to devalue the horn, the most common of which is dehorning. Game theory has been used to examine the interaction of poachers and wildlife managers where a manager can either `dehorn' their rhinos or leave the horn attached and poachers may behave `selectively' or `indiscriminately'. The approach described in this paper builds on this previous work and investigates the interactions between the poachers. We build an evolutionary game theoretic model and determine which strategy is preferred by a poacher in various different populations of poachers. The purpose of this work is to discover whether conditions which encourage the poachers to behave selectively exist, that is, they only kill those rhinos with full horns.
The analytical results show that full devaluation of all rhinos will likely lead to indiscriminate poaching. In turn it shows that devaluing of rhinos can only be effective when implemented along with a strong disincentive framework. This paper aims to contribute to the necessary research required for informed discussion about the lively debate on legalising rhino horn trade.
△ Less
Submitted 12 October, 2018; v1 submitted 20 December, 2017;
originally announced December 2017.
-
Convolutional neural networks automate detection for tracking of submicron scale particles in 2D and 3D
Authors:
Jay M. Newby,
Alison M. Schaefer,
Phoebe T. Lee,
M. Gregory Forest,
Samuel K. Lai
Abstract:
Particle tracking is a powerful biophysical tool that requires conversion of large video files into position time series, i.e. traces of the species of interest for data analysis. Current tracking methods, based on a limited set of input parameters to identify bright objects, are ill-equipped to handle the spectrum of spatiotemporal heterogeneity and poor signal-to-noise ratios typically presented…
▽ More
Particle tracking is a powerful biophysical tool that requires conversion of large video files into position time series, i.e. traces of the species of interest for data analysis. Current tracking methods, based on a limited set of input parameters to identify bright objects, are ill-equipped to handle the spectrum of spatiotemporal heterogeneity and poor signal-to-noise ratios typically presented by submicron species in complex biological environments. Extensive user involvement is frequently necessary to optimize and execute tracking methods, which is not only inefficient but introduces user bias. To develop a fully automated tracking method, we developed a convolutional neural network for particle localization from image data, comprised of over 6,000 parameters, and employed machine learning techniques to train the network on a diverse portfolio of video conditions. The neural network tracker provides unprecedented automation and accuracy, with exceptionally low false positive and false negative rates on both 2D and 3D simulated videos and 2D experimental videos of difficult-to-track species.
△ Less
Submitted 6 October, 2018; v1 submitted 10 April, 2017;
originally announced April 2017.
-
Effect of fractional blood flow on plasma skimming in the microvasculature
Authors:
Jiho Yang,
Sung Sic Yoo,
Tae-Rin Lee
Abstract:
Although redistribution of red blood cells at bifurcated vessels is highly dependent on flow rate, it is still challenging to quantitatively express the dependency of flow rate in plasma skimming due to nonlinear cellular interactions. We suggest a plasma skimming model that can involve the effect of fractional blood flow at each bifurcation point. For validating the new model, it is compared with…
▽ More
Although redistribution of red blood cells at bifurcated vessels is highly dependent on flow rate, it is still challenging to quantitatively express the dependency of flow rate in plasma skimming due to nonlinear cellular interactions. We suggest a plasma skimming model that can involve the effect of fractional blood flow at each bifurcation point. For validating the new model, it is compared with \textit{in vivo} data at single bifurcation points, as well as microvascular network systems. In the simulation results, the exponential decay of plasma skimming parameter, $M$, along fractional flow rate shows the best performance in both cases.
△ Less
Submitted 6 April, 2017; v1 submitted 23 January, 2017;
originally announced January 2017.
-
Generalized Plasma Skimming Model for Cells and Drug Carriers in the Microvasculature
Authors:
Tae-Rin Lee,
Sung Sic Yoo,
Jiho Yang
Abstract:
In microvascular transport, where both blood and drug carriers are involved, plasma skimming has a key role on changing hematocrit level and drug carrier concentration in capillary beds after continuous vessel bifurcation in the microvasculature. While there have been numerous studies on modeling the plasma skimming of blood, previous works lacked in consideration of its interaction with drug carr…
▽ More
In microvascular transport, where both blood and drug carriers are involved, plasma skimming has a key role on changing hematocrit level and drug carrier concentration in capillary beds after continuous vessel bifurcation in the microvasculature. While there have been numerous studies on modeling the plasma skimming of blood, previous works lacked in consideration of its interaction with drug carriers. In this paper, a generalized plasma skimming model is suggested to predict the redistributions of both the cells and drug carriers at each bifurcation. In order to examine its applicability, this new model was applied on a single bifurcation system to predict the redistribution of red blood cells and drug carriers. Furthermore, this model was tested at microvascular network level under different plasma skimming conditions for predicting the concentration of drug carriers. Based on these results, the applicability of this generalized plasma skimming model is fully discussed and future works along with the model's limitations are summarized.
△ Less
Submitted 20 September, 2016; v1 submitted 16 September, 2016;
originally announced September 2016.
-
DNA-Level Splice Junction Prediction using Deep Recurrent Neural Networks
Authors:
Byunghan Lee,
Taehoon Lee,
Byunggook Na,
Sungroh Yoon
Abstract:
A eukaryotic gene consists of multiple exons (protein coding regions) and introns (non-coding regions), and a splice junction refers to the boundary between a pair of exon and intron. Precise identification of spice junctions on a gene is important for deciphering its primary structure, function, and interaction. Experimental techniques for determining exon/intron boundaries include RNA-seq, which…
▽ More
A eukaryotic gene consists of multiple exons (protein coding regions) and introns (non-coding regions), and a splice junction refers to the boundary between a pair of exon and intron. Precise identification of spice junctions on a gene is important for deciphering its primary structure, function, and interaction. Experimental techniques for determining exon/intron boundaries include RNA-seq, which is often accompanied by computational approaches. Canonical splicing signals are known, but computational junction prediction still remains challenging because of a large number of false positives and other complications. In this paper, we exploit deep recurrent neural networks (RNNs) to model DNA sequences and to detect splice junctions thereon. We test various RNN units and architectures including long short-term memory units, gated recurrent units, and recently proposed iRNN for in-depth design space exploration. According to our experimental results, the proposed approach significantly outperforms not only conventional machine learning-based methods but also a recent state-of-the-art deep belief network-based technique in terms of prediction accuracy.
△ Less
Submitted 16 December, 2015;
originally announced December 2015.
-
Known Structure, Unknown Function: An Inquiry-based Undergraduate Biochemistry Laboratory Course
Authors:
Cynthia Gray,
Carol W. Price,
Christopher T. Lee,
Alison H. Dewald,
Matthew A. Cline,
Charles E. McAnany,
Linda Columbus,
Cameron Mura
Abstract:
Undergraduate biochemistry laboratory courses often do not provide students with an authentic research experience, particularly when the express purpose of the laboratory is purely instructional. However, an instructional laboratory course that is inquiry- and research-based could simultaneously impart scientific knowledge and foster a student's research expertise and confidence. We have developed…
▽ More
Undergraduate biochemistry laboratory courses often do not provide students with an authentic research experience, particularly when the express purpose of the laboratory is purely instructional. However, an instructional laboratory course that is inquiry- and research-based could simultaneously impart scientific knowledge and foster a student's research expertise and confidence. We have developed a year-long undergraduate biochemistry laboratory curriculum wherein students determine, via experiment and computation, the function of a protein of known three-dimensional structure. The first half of the course is inquiry-based and modular in design; students learn general biochemical techniques while gaining preparation for research experiments in the second semester. Having learned standard biochemical methods in the first semester, students independently pursue their own (original) research projects in the second semester. This new curriculum has yielded an improvement in student performance and confidence as assessed by various metrics. To disseminate teaching resources to students and instructors alike, a freely accessible Biochemistry Laboratory Education resource is available at http://biochemlab.org.
△ Less
Submitted 8 July, 2015;
originally announced July 2015.
-
Global asymptotic properties for a Leslie-Gower food chain model
Authors:
Andrei Korobeinikov,
William T. Lee
Abstract:
We study global asymptotic properties of a continuous time Leslie-Gower food chain model. We construct a Lyapunov function which enables us to establish global asymptotic stability of the unique coexisting equilibrium state.
We study global asymptotic properties of a continuous time Leslie-Gower food chain model. We construct a Lyapunov function which enables us to establish global asymptotic stability of the unique coexisting equilibrium state.
△ Less
Submitted 2 April, 2011;
originally announced April 2011.
