Showing 1–50 of 56 results for author: Roy, S

A Review of Link Prediction Applications in Network Biology

Authors: Ahmad F. Al Musawi, Satyaki Roy, Preetam Ghosh

Abstract: In the domain of network biology, the interactions among heterogeneous genomic and molecular entities are represented through networks. Link prediction (LP) methodologies are instrumental in inferring missing or prospective associations within these biological networks. In this review, we systematically dissect the attributes of local, centrality, and embedding-based LP approaches, applied to stat… ▽ More In the domain of network biology, the interactions among heterogeneous genomic and molecular entities are represented through networks. Link prediction (LP) methodologies are instrumental in inferring missing or prospective associations within these biological networks. In this review, we systematically dissect the attributes of local, centrality, and embedding-based LP approaches, applied to static and dynamic biological networks. We undertake an examination of the current applications of LP metrics for predicting links between diseases, genes, proteins, RNA, microbiomes, drugs, and neurons. We carry out comprehensive performance evaluations on established biological network datasets to show the practical applications of standard LP models. Moreover, we compare the similarity in prediction trends among the models and the specific network attributes that contribute to effective link prediction, before underscoring the role of LP in addressing the formidable challenges prevalent in biological systems, ranging from noise, bias, and data sparseness to interpretability. We conclude the review with an exploration of the essential characteristics expected from future LP models, poised to advance our comprehension of the intricate interactions governing biological systems. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Current and future directions in network biology

Authors: Marinka Zitnik, Michelle M. Li, Aydin Wells, Kimberly Glass, Deisy Morselli Gysi, Arjun Krishnan, T. M. Murali, Predrag Radivojac, Sushmita Roy, Anaïs Baudot, Serdar Bozdag, Danny Z. Chen, Lenore Cowen, Kapil Devkota, Anthony Gitter, Sara Gosline, Pengfei Gu, Pietro H. Guzzi, Heng Huang, Meng Jiang, Ziynet Nesibe Kesimoglu, Mehmet Koyuturk, Jian Ma, Alexander R. Pico, Nataša Pržulj , et al. (12 additional authors not shown)

Abstract: Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These challenges stem from various fa… ▽ More Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These challenges stem from various factors, notably the growing complexity and volume of data together with the increased diversity of data types describing different tiers of biological organization. We discuss prevailing research directions in network biology and highlight areas of inference and comparison of biological networks, multimodal data integration and heterogeneous networks, higher-order network analysis, machine learning on networks, and network-based personalized medicine. Following the overview of recent breakthroughs across these five areas, we offer a perspective on the future directions of network biology. Additionally, we offer insights into scientific communities, educational initiatives, and the importance of fostering diversity within the field. This paper establishes a roadmap for an immediate and long-term vision for network biology. △ Less

Submitted 11 June, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 52 pages, 6 figures, 1 table

Interplay of degeneracy and non-degeneracy in fluctuations propagation in coherent feed-forward loop motif

Authors: Tuhin Subhra Roy, Mintu Nandi, Pinaki Chaudhury, Sudip Chattopadhyay, Suman K Banik

Abstract: We present a stochastic framework to decipher fluctuations propagation in classes of coherent feed-forward loops. The systematic contribution of the direct (one-step) and indirect (two-step) pathways is considered to quantify fluctuations of the output node. We also consider both additive and multiplicative integration mechanisms of the two parallel pathways (one-step and two-step). Analytical exp… ▽ More We present a stochastic framework to decipher fluctuations propagation in classes of coherent feed-forward loops. The systematic contribution of the direct (one-step) and indirect (two-step) pathways is considered to quantify fluctuations of the output node. We also consider both additive and multiplicative integration mechanisms of the two parallel pathways (one-step and two-step). Analytical expression of the output node's coefficient of variation shows contributions of intrinsic, one-step, two-step, and cross-interaction in closed form. We observe a diverse range of degeneracy and non-degeneracy in each of the decomposed fluctuations term and their contribution to the overall output fluctuations of each coherent feed-forward loop motif. Analysis of output fluctuations reveals a maximal level of fluctuations of the coherent feed-forward loop motif of type 1. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 20 pages, 4 figures

Journal ref: J. Stat. Mech. 2023 (2023) 093502

Image background assessment as a novel technique for insect microhabitat identification

Authors: Sesa Singha Roy, Reid Tingley, Alan Dorin

Abstract: The effects of climate change, urbanisation and agriculture are changing the way insects occupy habitats. Some species may utilise anthropogenic microhabitat features for their existence, either because they prefer them to natural features, or because of no choice. Other species are dependent on natural microhabitats. Identifying and analysing these insects' use of natural and anthropogenic microh… ▽ More The effects of climate change, urbanisation and agriculture are changing the way insects occupy habitats. Some species may utilise anthropogenic microhabitat features for their existence, either because they prefer them to natural features, or because of no choice. Other species are dependent on natural microhabitats. Identifying and analysing these insects' use of natural and anthropogenic microhabitats is important to assess their responses to a changing environment, for improving pollination and managing invasive pests. Traditional studies of insect microhabitat use can now be supplemented by machine learning-based insect image analysis. Typically, research has focused on automatic insect classification, but valuable data in image backgrounds has been ignored. In this research, we analysed the image backgrounds available on the ALA database to determine their microhabitats. We analysed the microhabitats of three insect species common across Australia: Drone flies, European honeybees and European wasps. Image backgrounds were classified as natural or anthropogenic microhabitats using computer vision and machine learning tools benchmarked against a manual classification algorithm. We found flies and honeybees in natural microhabitats, confirming their need for natural havens within cities. Wasps were commonly seen in anthropogenic microhabitats. Results show these insects are well adapted to survive in cities. Management of this invasive pest requires a thoughtful reduction of their access to human-provided resources. The assessment of insect image backgrounds is instructive to document the use of microhabitats by insects. The method offers insight that is increasingly vital for biodiversity management as urbanisation continues to encroach on natural ecosystems and we must consciously provide resources within built environments to maintain insect biodiversity and manage invasive pests. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: Submitted in Ecological Informatics journal, first review completed, 19 pages, 10 figures

An anatomy-based V1 model: Extraction of Low-level Features, Reduction of distortion and a V1-inspired SOM

Authors: Suvam Roy, Nikhil Ranjan Pal

Abstract: We present a model of the primary visual cortex V1, guided by anatomical experiments. Unlike most machine learning systems our goal is not to maximize accuracy but to realize a system more aligned to biological systems. Our model consists of the V1 layers 4, 2/3, and 5, with inter-layer connections between them in accordance with the anatomy. We further include the orientation selectivity of the V… ▽ More We present a model of the primary visual cortex V1, guided by anatomical experiments. Unlike most machine learning systems our goal is not to maximize accuracy but to realize a system more aligned to biological systems. Our model consists of the V1 layers 4, 2/3, and 5, with inter-layer connections between them in accordance with the anatomy. We further include the orientation selectivity of the V1 neurons and lateral influences in each layer. Our V1 model, when applied to the BSDS500 ground truth images (indicating LGN contour detection before V1), can extract low-level features from the images and perform a significant amount of distortion reduction. As a follow-up to our V1 model, we propose a V1-inspired self-organizing map algorithm (V1-SOM), where the weight update of each neuron gets influenced by its neighbors. V1-SOM can tolerate noisy inputs as well as noise in the weight updates better than SOM and shows a similar level of performance when trained with high dimensional data such as the MNIST dataset. Finally, when we applied V1 processing to the MNIST dataset to extract low-level features and trained V1-SOM with the modified MNIST dataset, the quantization error was significantly reduced. Our results support the hypothesis that the ventral stream performs gradual untangling of input spaces. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Deep learning for reconstructing protein structures from cryo-EM density maps: recent advances and future directions

Authors: Nabin Giri, Raj S. Roy, Jianlin Cheng

Abstract: Cryo-Electron Microscopy (cryo-EM) has emerged as a key technology to determine the structure of proteins, particularly large protein complexes and assemblies in recent years. A key challenge in cryo-EM data analysis is to automatically reconstruct accurate protein structures from cryo-EM density maps. In this review, we briefly overview various deep learning methods for building protein structure… ▽ More Cryo-Electron Microscopy (cryo-EM) has emerged as a key technology to determine the structure of proteins, particularly large protein complexes and assemblies in recent years. A key challenge in cryo-EM data analysis is to automatically reconstruct accurate protein structures from cryo-EM density maps. In this review, we briefly overview various deep learning methods for building protein structures from cryo-EM density maps, analyze their impact, and discuss the challenges of preparing high-quality data sets for training deep learning models. Looking into the future, more advanced deep learning models of effectively integrating cryo-EM data with other sources of complementary data such as protein sequences and AlphaFold-predicted structures need to be developed to further advance the field. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Journal ref: Current Opinion in Structural Biology Volume 79, April 2023, 102536

DRLComplex: Reconstruction of protein quaternary structures using deep reinforcement learning

Authors: Elham Soltanikazemi, Raj S. Roy, Farhan Quadir, Nabin Giri, Alex Morehead, Jianlin Cheng

Abstract: Predicted inter-chain residue-residue contacts can be used to build the quaternary structure of protein complexes from scratch. However, only a small number of methods have been developed to reconstruct protein quaternary structures using predicted inter-chain contacts. Here, we present an agent-based self-learning method based on deep reinforcement learning (DRLComplex) to build protein complex s… ▽ More Predicted inter-chain residue-residue contacts can be used to build the quaternary structure of protein complexes from scratch. However, only a small number of methods have been developed to reconstruct protein quaternary structures using predicted inter-chain contacts. Here, we present an agent-based self-learning method based on deep reinforcement learning (DRLComplex) to build protein complex structures using inter-chain contacts as distance constraints. We rigorously tested DRLComplex on two standard datasets of homodimeric and heterodimeric protein complexes (i.e., the CASP-CAPRI homodimer and Std_32 heterodimer datasets) using both true and predicted interchain contacts as inputs. Utilizing true contacts as input, DRLComplex achieved high average TM-scores of 0.9895 and 0.9881 and a low average interface RMSD (I_RMSD) of 0.2197 and 0.92 on the two datasets, respectively. When predicted contacts are used, the method achieves TM-scores of 0.73 and 0.76 for homodimers and heterodimers, respectively. Our experiments find that the accuracy of reconstructed quaternary structures depends on the accuracy of the contact predictions. Compared to other optimization methods for reconstructing quaternary structures from inter-chain contacts, DRLComplex performs similar to an advanced gradient descent method and better than a Markov Chain Monte Carlo simulation method and a simulated annealing-based method, validating the effectiveness of DRLComplex for quaternary reconstruction of protein complexes. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: 20 pages, 8 figures, 12 tables. Under review

ACM Class: I.2.1; J.3

Diversifying the Genomic Data Science Research Community

Authors: The Genomic Data Science Community Network, Rosa Alcazar, Maria Alvarez, Rachel Arnold, Mentewab Ayalew, Lyle G. Best, Michael C. Campbell, Kamal Chowdhury, Katherine E. L. Cox, Christina Daulton, Youping Deng, Carla Easter, Karla Fuller, Shazia Tabassum Hakim, Ava M. Hoffman, Natalie Kucher, Andrew Lee, Joslynn Lee, Jeffrey T. Leek, Robert Meller, Loyda B. Méndez, Miguel P. Méndez-González, Stephen Mosher, Michele Nishiguchi, Siddharth Pratap , et al. (13 additional authors not shown)

Abstract: Over the last 20 years, there has been an explosion of genomic data collected for disease association, functional analyses, and other large-scale discoveries. At the same time, there have been revolutions in cloud computing that enable computational and data science research, while making data accessible to anyone with a web browser and an internet connection. However, students at institutions wit… ▽ More Over the last 20 years, there has been an explosion of genomic data collected for disease association, functional analyses, and other large-scale discoveries. At the same time, there have been revolutions in cloud computing that enable computational and data science research, while making data accessible to anyone with a web browser and an internet connection. However, students at institutions with limited resources have received relatively little exposure to curricula or professional development opportunities that lead to careers in genomic data science. To broaden participation in genomics research, the scientific community needs to support students, faculty, and administrators at Underserved Institutions (UIs) including Community Colleges, Historically Black Colleges and Universities, Hispanic-Serving Institutions, and Tribal Colleges and Universities in taking advantage of these tools in local educational and research programs. We have formed the Genomic Data Science Community Network (http://www.gdscn.org/) to identify opportunities and support broadening access to cloud-enabled genomic data science. Here, we provide a summary of the priorities for faculty members at UIs, as well as administrators, funders, and R1 researchers to consider as we create a more diverse genomic data science community. △ Less

Submitted 9 June, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 42 pages, 3 figures

Information transmission in a two-step cascade: Interplay of activation and repression

Authors: Tuhin Subhra Roy, Mintu Nandi, Ayan Biswas, Pinaki Chaudhury, Suman K Banik

Abstract: We present an information-theoretic formalism to study signal transduction in four architectural variants of a model two-step cascade with increasing input population. Our results categorize these four types into two classes depending upon the effect played out by activation and repression on mutual information, net synergy, and signal-to-noise ratio. Within the Gaussian framework and using the li… ▽ More We present an information-theoretic formalism to study signal transduction in four architectural variants of a model two-step cascade with increasing input population. Our results categorize these four types into two classes depending upon the effect played out by activation and repression on mutual information, net synergy, and signal-to-noise ratio. Within the Gaussian framework and using the linear noise approximation, we derive the analytic expressions for these metrics to establish their underlying relationships in terms of the biochemical parameters. We also verify our approximations through stochastic simulations. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 10 pages, 4 figures

A Fokker-Planck feedback control framework for optimal personalized therapies in colon cancer-induced angiogenesis

Authors: Souvik Roy, Zui Pan, Suvra Pal

Abstract: In this paper, a new framework for obtaining personalized optimal treatment strategies in colon cancer-induced angiogenesis is presented. The dynamics of colon cancer is given by a Itó stochastic process, which helps in modeling the randomness present in the system. The stochastic dynamics is then represented by the Fokker-Planck (FP) partial differential equation (PDE) that governs the evolution… ▽ More In this paper, a new framework for obtaining personalized optimal treatment strategies in colon cancer-induced angiogenesis is presented. The dynamics of colon cancer is given by a Itó stochastic process, which helps in modeling the randomness present in the system. The stochastic dynamics is then represented by the Fokker-Planck (FP) partial differential equation (PDE) that governs the evolution of the associated probability density function. The optimal therapies are obtained using a three step procedure. First, a finite dimensional FP-constrained optimization problem is formulated that takes input individual noisy patient data, and is solved to obtain the unknown parameters corresponding to the individual tumor characteristics. Next, a sensitivity analysis of the optimal parameter set is used to determine the parameters to be controlled, thus, helping in assessing the types of treatment therapies. Finally, a feedback FP control problem is solved to determine the optimal combination therapies. Numerical results with the combination drug, comprising of Bevacizumab and Capecitabine, demonstrate the efficiency of the proposed framework. △ Less

Submitted 13 November, 2021; v1 submitted 31 August, 2021; originally announced September 2021.

MSC Class: 35R30; 49J20; 49K20; 62D99; 65M08; 82C31

Optimal Background Correction in Double Quantum Coherence Electron Spin Resonance Spectroscopy for Accurate Data Analysis

Authors: Aritro Sinha Roy

Abstract: Electron spin resonance (ESR) pulsed dipolar spectroscopy (PDS) is used in protein 3D structure determination. However, the accuracy of the signal analysis depends heavily on the background correction process. In this work, we derive the functional forms of double quantum coherence (DQC) ESR signal in typical frozen samples of micro-molar concentration, quantifying both the intramolecular and the… ▽ More Electron spin resonance (ESR) pulsed dipolar spectroscopy (PDS) is used in protein 3D structure determination. However, the accuracy of the signal analysis depends heavily on the background correction process. In this work, we derive the functional forms of double quantum coherence (DQC) ESR signal in typical frozen samples of micro-molar concentration, quantifying both the intramolecular and the background contributions. △ Less

Submitted 29 June, 2021; originally announced June 2021.

An Automated Global Method for Extraction of Distance Distributions from Electron Spin Resonance Pulsed Dipolar Signals

Authors: Aritro Sinha Roy

Abstract: Electron spin resonance (ESR) pulsed dipolar spectroscopy (PDS) is used effectively in measuring nano-meter range distances for protein structure prediction. The current global approach in extracting the distance distribution from time domain PDS signal has multiple limitations. We present a parameter free global method, which is more efficacious and less sensitive to signal noise compared to the… ▽ More Electron spin resonance (ESR) pulsed dipolar spectroscopy (PDS) is used effectively in measuring nano-meter range distances for protein structure prediction. The current global approach in extracting the distance distribution from time domain PDS signal has multiple limitations. We present a parameter free global method, which is more efficacious and less sensitive to signal noise compared to the current method. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 4 pages, 3 figures

Analyzing Host-Viral Interactome of SARS-CoV-2 for Identifying Vulnerable Host Proteins during COVID-19 Pathogenesis

Authors: Jayanta Kumar Das, Swarup Roy, Pietro Hiram Guzzi

Abstract: The development of therapeutic targets for COVID-19 treatment is based on the understanding of the molecular mechanism of pathogenesis. The identification of genes and proteins involved in the infection mechanism is the key to shed out light into the complex molecular mechanisms. The combined effort of many laboratories distributed throughout the world has produced the accumulation of both protein… ▽ More The development of therapeutic targets for COVID-19 treatment is based on the understanding of the molecular mechanism of pathogenesis. The identification of genes and proteins involved in the infection mechanism is the key to shed out light into the complex molecular mechanisms. The combined effort of many laboratories distributed throughout the world has produced the accumulation of both protein and genetic interactions. In this work we integrate these available results and we obtain an host protein-protein interaction network composed by 1432 human proteins. We calculate network centrality measures to identify key proteins. Then we perform functional enrichment of central proteins. We observed that the identified proteins are mostly associated with several crucial pathways, including cellular process, signalling transduction, neurodegenerative disease. Finally, we focused on proteins involved in causing disease in the human respiratory tract. We conclude that COVID19 is a complex disease, and we highlighted many potential therapeutic targets including RBX1, HSPA5, ITCH, RAB7A, RAB5A, RAB8A, PSMC5, CAPZB, CANX, IGF2R, HSPA1A, which are central and also associated with multiple diseases △ Less

Submitted 5 February, 2021; originally announced February 2021.

Scale-free networks may not necessarily witness cooperation

Authors: Deep Nath, Saptarshi Sinha, Soumen Roy

Abstract: Networks with a scale-free degree distribution are widely thought to promote cooperation in various games. Herein, by studying the well-known prisoner's dilemma game, we demonstrate that this need not necessarily be true. For the very same degree sequence and degree distribution, we present a variety of possible behaviour. We reassess the perceived importance of hubs in a network towards the maint… ▽ More Networks with a scale-free degree distribution are widely thought to promote cooperation in various games. Herein, by studying the well-known prisoner's dilemma game, we demonstrate that this need not necessarily be true. For the very same degree sequence and degree distribution, we present a variety of possible behaviour. We reassess the perceived importance of hubs in a network towards the maintenance of cooperation. We also reevaluate the dependence of cooperation on network clustering and assortativity. △ Less

Submitted 19 June, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: Accepted to Europhysics Letters (EPL)

Journal ref: Europhysics Letters, 134, 60002 (2021)

A Generalized Epidemiological Model for COVID-19 with Dynamic and Asymptomatic Population

Authors: Anirban Ghatak, Shivshanker Singh Patel, Soham Bonnerjee, Subhrajyoty Roy

Abstract: In this paper, we develop an extension of standard epidemiological models, suitable for COVID-19. This extension incorporates the transmission due to pre-symptomatic or asymptomatic carriers of the virus. Furthermore, this model also captures the spread of the disease due to the movement of people to/from different administrative boundaries within a country. The model describes the probabilistic r… ▽ More In this paper, we develop an extension of standard epidemiological models, suitable for COVID-19. This extension incorporates the transmission due to pre-symptomatic or asymptomatic carriers of the virus. Furthermore, this model also captures the spread of the disease due to the movement of people to/from different administrative boundaries within a country. The model describes the probabilistic rise in the number of confirmed cases due to the concomitant effects of (incipient) human transmission and multiple compartments. The associated parameters in the model can help architect the public health policy and operational management of the pandemic. For instance, this model demonstrates that increasing the testing for symptomatic patients does not have any major effect on the progression of the pandemic, but testing rate of the asymptomatic population has an extremely crucial role to play. The model is executed using the data obtained for the state of Chhattisgarh in the Republic of India. The model is shown to have significantly better predictive capability than the other epidemiological models. This model can be readily applied to any administrative boundary (state or country). Moreover, this model can be applied for any other epidemic as well. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Feasibility Assessment of an Optically Powered Digital Retinal Prosthesis Architecture for Retinal Ganglion Cell Stimulation

Authors: William Lemaire, Maher Benhouria, Konin Koua, Wei Tong, Gabriel Martin-Hardy, Melanie Stamp, Kumaravelu Ganesan, Louis-Philippe Gauthier, Marwan Besrour, Arman Ahnood, David John Garrett, Sébastien Roy, Michael Ibbotson, Steven Prawer, Réjean Fontaine

Abstract: Clinical trials previously demonstrated the notable capacity to elicit visual percepts in blind patients affected with retinal diseases by electrically stimulating the remaining neurons on the retina. However, these implants restored very limited visual acuity and required transcutaneous cables traversing the eyeball, leading to reduced reliability and complex surgery with high postoperative infec… ▽ More Clinical trials previously demonstrated the notable capacity to elicit visual percepts in blind patients affected with retinal diseases by electrically stimulating the remaining neurons on the retina. However, these implants restored very limited visual acuity and required transcutaneous cables traversing the eyeball, leading to reduced reliability and complex surgery with high postoperative infection risks. To overcome the limitations imposed by cables, a retinal implant architecture in which near-infrared illumination carries both power and data through the pupil to a digital stimulation controller is presented. A high efficiency multi-junction photovoltaic cell transduces the optical power to a CMOS stimulator capable of delivering flexible interleaved sequential stimulation through a diamond microelectrode array. To demonstrate the capacity to elicit a neural response with this approach while complying with the optical irradiance limit at the pupil, fluorescence imaging with a calcium indicator is used on a degenerate rat retina. The power delivered by the laser at the permissible irradiance of 4 mW/mm2 at 850 nm is shown to be sufficient to both power the stimulator ASIC and elicit a response in retinal ganglion cells (RGCs), with the ability to generate of up to 35 000 pulses per second at the average stimulation threshold. This confirms the feasibility of generating a response in RGCs with an infrared-powered digital architecture capable of delivering complex sequential stimulation patterns at high repetition rates, albeit with some limitations. △ Less

Submitted 13 October, 2023; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: 11 pages, 13 figures

Mass testing and proactiveness affect epidemic spreading

Authors: Saptarshi Sinha, Deep Nath, Soumen Roy

Abstract: The detection and management of diseases become quite complicated when pathogens contain asymptomatic phenotypes amongst their ranks, as evident during the recent COVID-19 pandemic. Spreading of diseases has been studied extensively under the paradigm of Susceptible - Infected - Recovered - Deceased (SIRD) dynamics. Various game-theoretic approaches have also addressed disease spread, many of whic… ▽ More The detection and management of diseases become quite complicated when pathogens contain asymptomatic phenotypes amongst their ranks, as evident during the recent COVID-19 pandemic. Spreading of diseases has been studied extensively under the paradigm of Susceptible - Infected - Recovered - Deceased (SIRD) dynamics. Various game-theoretic approaches have also addressed disease spread, many of which consider S, I, R, and D as strategies rather than as states. Remarkably, most studies from the above approaches do not account for the distinction between the symptomatic or asymptomatic aspect of the disease. It is well-known that precautionary measures like washing hands, wearing masks and social distancing significantly mitigate the spread of many contagious diseases. Herein, we consider the adoption of such precautions as strategies and treat S, I, R, and D as states. We also attempt to capture the differences in epidemic spreading arising from symptomatic and asymptomatic diseases on various network topologies. Through extensive computer simulations, we examine that the cost of maintaining precautionary measures as well as the extent of mass testing in a population affects the final fraction of socially responsible individuals. We observe that the lack of mass testing could potentially lead to a pandemic in case of asymptomatic diseases. Network topology also seems to play an important role. We further observe that the final fraction of proactive individuals depends on the initial fraction of both infected as well as proactive individuals. Additionally, edge density can significantly influence the overall outcome. Our findings are in broad agreement with the lessons learnt from the ongoing COVID-19 pandemic. △ Less

Submitted 14 May, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: To appear in the special issue on "Emerging trends in modeling and analysis of complex systems" of Journal of the Indian Institute of Science (J IISc) (2021)

Journal ref: Journal of the Indian Institute of Science, 101, 371-380 (2021)

An India-specific Compartmental Model for Covid-19: Projections and Intervention Strategies by Incorporating Geographical, Infrastructural and Response Heterogeneity

Authors: Sanit Gupta, Sahil Shah, Sumit Chaturvedi, Pranav Thakkar, Parvinder Solanki, Soham Dibyachintan, Sandeepan Roy, M. B. Sushma, Adwait Godbole, Noufal Jaseem, Pradumn Kumar, Sucheta Ravikanti, Aritra Das, Giridhara R. Babu, Tarun Bhatnagar, Avijit Maji, Mithun K. Mitra, Sai Vinjanampathy

Abstract: We present a compartmental meta-population model for the spread of Covid-19 in India. Our model simulates populations at a district or state level using an epidemiological model that is appropriate to Covid-19. Different districts are connected by a transportation matrix developed using available census data. We introduce uncertainties in the testing rates into the model that takes into account th… ▽ More We present a compartmental meta-population model for the spread of Covid-19 in India. Our model simulates populations at a district or state level using an epidemiological model that is appropriate to Covid-19. Different districts are connected by a transportation matrix developed using available census data. We introduce uncertainties in the testing rates into the model that takes into account the disparate responses of the different states to the epidemic and also factors in the state of the public healthcare system. Our model allows us to generate qualitative projections of Covid-19 spread in India, and further allows us to investigate the effects of different proposed interventions. By building in heterogeneity at geographical and infrastructural levels and in local responses, our model aims to capture some of the complexity of epidemiological modeling appropriate to a diverse country such as India. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: code and updates available at https://github.com/vinjanampathy/Heterogenous_COVID19_India_Model

Topology dependent payoffs can lead to escape from prisoner's dilemma

Authors: Saptarshi Sinha, Deep Nath, Soumen Roy

Abstract: The maintenance of cooperation in the presence of spatial restrictions has been studied extensively. It is well-established that the underlying graph topology can significantly influence the outcome of games on graphs. Maintenance of cooperation could be difficult, especially in the absence of spatial restrictions. The evolution of cooperation would naturally depend on payoffs. However, payoffs ar… ▽ More The maintenance of cooperation in the presence of spatial restrictions has been studied extensively. It is well-established that the underlying graph topology can significantly influence the outcome of games on graphs. Maintenance of cooperation could be difficult, especially in the absence of spatial restrictions. The evolution of cooperation would naturally depend on payoffs. However, payoffs are generally considered to be invariant in a given game. A natural yet unexplored question is whether the topology of the underlying structures on which the games are played, possesses no role whatsoever in the determination of payoffs. Herein, we introduce the notion of cooperator graphs and defector graphs as well as a new form of game payoff, which is weakly dependent on the underlying network topology. These concepts are inspired by the well-known microbial phenomenon of quorum sensing. We demonstrate that even with such a weak dependence, the fundamental game dynamics and indeed the very nature of the game may be altered. Such changes in the nature of a game have been well-reported in theoretical and experimental studies. △ Less

Submitted 12 March, 2021; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: Replaced by the version accepted to European Physical Journal B

Journal ref: The European Physical Journal B, 94, 80 (2021)

Analyzing COVID-19 pandemic with a new growth model for population ecology

Authors: Deeptak Biswas, Sulagna Roy

Abstract: We have proposed a new form of growth rate for population ecology. Generally, the growth rate is dependent on the size of the population at that particular epoch. We have introduced an alternative time-dependent form of growth rate. This form satisfies essential conditions to represent population growth and can be an alternative form for growth models to analyze population dynamics. We have employ… ▽ More We have proposed a new form of growth rate for population ecology. Generally, the growth rate is dependent on the size of the population at that particular epoch. We have introduced an alternative time-dependent form of growth rate. This form satisfies essential conditions to represent population growth and can be an alternative form for growth models to analyze population dynamics. We have employed the generalized Richards model as a guideline to compare our results. Further, we have applied our model in the case of epidemics. To check the efficacy of our model, we have verified the 2003 SARS data. This model has estimated the final epidemic size with good accuracy. Thereafter, we intend to describe the present COVID-2019 pandemic. We have performed our analysis with data for Italy, Spain, and Germany. Following, we have tried to predict the number of COVID-19 cases and the turning point for the USA, UK, and India. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Interpretable brain age prediction using linear latent variable models of functional connectivity

Authors: Ricardo Pio Monti, Alex Gibberd, Sandipan Roy, Matt Nunes, Romy Lorenz, Robert Leech, Takeshi Ogawa, Motoaki Kawanabe, Aapo Hyvarinen

Abstract: Neuroimaging-driven prediction of brain age, defined as the predicted biological age of a subject using only brain imaging data, is an exciting avenue of research. In this work we seek to build models of brain age based on functional connectivity while prioritizing model interpretability and understanding. This way, the models serve to both provide accurate estimates of brain age as well as allow… ▽ More Neuroimaging-driven prediction of brain age, defined as the predicted biological age of a subject using only brain imaging data, is an exciting avenue of research. In this work we seek to build models of brain age based on functional connectivity while prioritizing model interpretability and understanding. This way, the models serve to both provide accurate estimates of brain age as well as allow us to investigate changes in functional connectivity which occur during the ageing process. The methods proposed in this work consist of a two-step procedure: first, linear latent variable models, such as PCA and its extensions, are employed to learn reproducible functional connectivity networks present across a cohort of subjects. The activity within each network is subsequently employed as a feature in a linear regression model to predict brain age. The proposed framework is employed on the data from the CamCAN repository and the inferred brain age models are further demonstrated to generalize using data from two open-access repositories: the Human Connectome Project and the ATR Wide-Age-Range. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 21 pages, 11 figures

SeizureNet: Multi-Spectral Deep Feature Learning for Seizure Type Classification

Authors: Umar Asif, Subhrajit Roy, Jianbin Tang, Stefan Harrer

Abstract: Automatic classification of epileptic seizure types in electroencephalograms (EEGs) data can enable more precise diagnosis and efficient management of the disease. This task is challenging due to factors such as low signal-to-noise ratios, signal artefacts, high variance in seizure semiology among epileptic patients, and limited availability of clinical data. To overcome these challenges, in this… ▽ More Automatic classification of epileptic seizure types in electroencephalograms (EEGs) data can enable more precise diagnosis and efficient management of the disease. This task is challenging due to factors such as low signal-to-noise ratios, signal artefacts, high variance in seizure semiology among epileptic patients, and limited availability of clinical data. To overcome these challenges, in this paper, we present SeizureNet, a deep learning framework which learns multi-spectral feature embeddings using an ensemble architecture for cross-patient seizure type classification. We used the recently released TUH EEG Seizure Corpus (V1.4.0 and V1.5.2) to evaluate the performance of SeizureNet. Experiments show that SeizureNet can reach a weighted F1 score of up to 0.94 for seizure-wise cross validation and 0.59 for patient-wise cross validation for scalp EEG based multi-class seizure type classification. We also show that the high-level feature embeddings learnt by SeizureNet considerably improve the accuracy of smaller networks through knowledge distillation for applications with low-memory constraints. △ Less

Submitted 29 September, 2020; v1 submitted 7 March, 2019; originally announced March 2019.

Seizure Type Classification using EEG signals and Machine Learning: Setting a benchmark

Authors: Subhrajit Roy, Umar Asif, Jianbin Tang, Stefan Harrer

Abstract: Accurate classification of seizure types plays a crucial role in the treatment and disease management of epileptic patients. Epileptic seizure types not only impact the choice of drugs but also the range of activities a patient can safely engage in. With recent advances being made towards artificial intelligence enabled automatic seizure detection, the next frontier is the automatic classification… ▽ More Accurate classification of seizure types plays a crucial role in the treatment and disease management of epileptic patients. Epileptic seizure types not only impact the choice of drugs but also the range of activities a patient can safely engage in. With recent advances being made towards artificial intelligence enabled automatic seizure detection, the next frontier is the automatic classification of seizure types. On that note, in this paper, we explore the application of machine learning algorithms for multi-class seizure type classification. We used the recently released TUH EEG seizure corpus (V1.4.0 and V1.5.2) and conducted a thorough search space exploration to evaluate the performance of a combination of various pre-processing techniques, machine learning algorithms, and corresponding hyperparameters on this task. We show that our algorithms can reach a weighted $F1$ score of up to 0.901 for seizure-wise cross validation and 0.561 for patient-wise cross validation thereby setting a benchmark for scalp EEG based multi-class seizure type classification. △ Less

Submitted 11 August, 2020; v1 submitted 3 February, 2019; originally announced February 2019.

Comments: 5 pages, 2 figure, 4 table

Neural (EEG) Response during Creation and Appreciation: A Novel Study with Hindustani Raga Music

Authors: Archi Banerjee, Shankha Sanyal, Souparno Roy, Sourya Sengupta, Sayan Biswas, Sayan Nag, Ranjan Sengupta, Dipak Ghosh

Abstract: What happens inside the performers brain when he is performing and composing a particular raga. Are there some specific regions in brain which are activated when an artist is creating or imaging a raga in his brain. Do the regions remain the same when the artist is listening to the same raga sung by him. These are the questions that perplexed neuroscientists for a long time. In this study we striv… ▽ More What happens inside the performers brain when he is performing and composing a particular raga. Are there some specific regions in brain which are activated when an artist is creating or imaging a raga in his brain. Do the regions remain the same when the artist is listening to the same raga sung by him. These are the questions that perplexed neuroscientists for a long time. In this study we strive to answer these questions by using latest state of the art techniques to assess brain response. An EEG experiment was conducted for two eminent performers of Indian classical music, when they mentally created the imagery of a raga Jay Jayanti in their mind, as well as when they listened to the same raga. The beauty of Hindustani music lies in the fact that the musician is himself the composer and recreates the imagery of the raga in his mind while performing, hence the scope of creative improvisations are immense. The alpha and theta frequency rhythms were segregated from each of the time series data and analyzed using robust non MFDXA technique to quantitatively assess the degree of cross-correlation of each EEG frequency rhythm in different combination of electrodes from frontal, occipital and temporal lobes. A strong response was found in the occipital and fronto occipital region during mental improvisation of the raga, which is an interesting revelation of this study. Strong retentive features were obtained in regard to both alpha and theta rhythms in musical listening in the fronto temporal and occipital temporal region while the features were almost absent in the thinking part. Further, other specific regions have been identified separately for the two separate conditions in which the correlations among the different lobes were the strongest. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 14 pages, 5 figures

Gestalt Phenomenon in Music? A Neurocognitive Physics Study with EEG

Authors: Shankha Sanyal, Archi Banerjee, Souparno Roy, Sourya Sengupta, Sayan Biswas, Sayan Nag, Ranjan Sengupta, Dipak Ghosh

Abstract: The term gestalt has been widely used in the field of psychology which defined the perception of human mind to group any object not in part but as a unified whole. Music in general is polytonic i.e. a combination of a number of pure tones (frequencies) mixed together in a manner that sounds harmonius. The study of human brain response due to different frequency groups of acoustic signal can give u… ▽ More The term gestalt has been widely used in the field of psychology which defined the perception of human mind to group any object not in part but as a unified whole. Music in general is polytonic i.e. a combination of a number of pure tones (frequencies) mixed together in a manner that sounds harmonius. The study of human brain response due to different frequency groups of acoustic signal can give us an excellent insight regarding the neural and functional architecture of brain functions. In this work we have tried to analyze the effect of different frequency bands of music on the various frequency rhythms of human brain obtained from EEG data of 5 participants. Four (4) widely popular Rabindrasangeet clips were subjected to Wavelet Transform method for extracting five resonant frequency bands from the original music signal. These resonant frequency bands were presented to the subjects as auditory stimulus and EEG signals recorded simultaneously in 19 different locations of the brain. The recorded EEG signals were noise cleaned and subjected to Multifractal Detrended Fluctuation Analysis (MFDFA) technique on the alpha, theta and gamma frequency range. Thus, we obtained the complexity values (in the form of multifractal spectral width) in alpha, theta and gamma EEG rhythms corresponding to different frequency bands of music. We obtain frequency specific arousal based response in different lobes of brain as well as in specific EEG bands corresponding to musical stimuli. This revelation can be of immense importance when it comes to the field of cognitive music therapy. △ Less

Submitted 19 March, 2017; originally announced March 2017.

Comments: 14 Pages, 5 Figures, Presented in International Conference on Creativity and Cognition in Art and Design, NIMHANS, Bangalore; 19-21 January, 2017

New Empirical Evidence on Disjunction Effect and Cultural Dependence

Authors: Indranil Mukhopadhyay, Nithin Nagaraj, Sisir Roy

Abstract: We perform new experiment using almost the same sample size considered by Tversky and Shafir to test the validity of classical probability theory in decision making. The results clearly indicate that the disjunction effect depends also on culture and more specifically on gender (females rather than males). We did more statistical analysis rather that putting the actual values done by previous auth… ▽ More We perform new experiment using almost the same sample size considered by Tversky and Shafir to test the validity of classical probability theory in decision making. The results clearly indicate that the disjunction effect depends also on culture and more specifically on gender (females rather than males). We did more statistical analysis rather that putting the actual values done by previous authors. We propose different kind of disjunction effect i.e. strong and weak based on our statistical analysis. △ Less

Submitted 1 March, 2017; originally announced March 2017.

Comments: 14 pages

Ambiguity invokes Creativity : looking through Quantum physics

Authors: Souparno Roy, Archi Banerjee, Ranjan Sengupta, Dipak Ghosh

Abstract: Creativity, defined as the tendency to generate or recognize new ideas or alternatives and to make connections between seemingly unrelated phenomena, is too vast a horizon to be summed up in such a simple sentence. The extreme abstractness of creativity makes it harder to quantify in its entirety. Yet, a lot of efforts have been made both by psychologists and neurobiologists to identify its signat… ▽ More Creativity, defined as the tendency to generate or recognize new ideas or alternatives and to make connections between seemingly unrelated phenomena, is too vast a horizon to be summed up in such a simple sentence. The extreme abstractness of creativity makes it harder to quantify in its entirety. Yet, a lot of efforts have been made both by psychologists and neurobiologists to identify its signature. A general conformity is expressed in the Free association theory, i.e. the more freely a persons conceptual nodes are connected, the more divergent thinker (also, creative) he or she is. Also, tolerance of ambiguity is found to be related to divergent thinking. In this study, we approach the problem of creativity from a theoretical physics standpoint. Theoretically, for the initial conceptual state, the next jump to any other node is equally probable and non-deterministic. Repeated intervention of external stimulus (analogous to a measurement) is responsible for such jumps. And to study such a non-deterministic system with continuous measurements, Quantum theory has been proven the most successful, time and again. We suggest that this collection of nodes form a system which is likely to be governed by quantum physics and specify the transformations which could help explain the conceptual jump between states. Our argument, from the point of view of physics is that the initial evolution of the creative process is identical, person or field independent. To answer the next obvious question about individual creativity, we hypothesize that the quantum system, under continuous measurements (in the form of external stimuli) evolves with chaotic dynamics, hence separating a painter from a musician. Possible experimental methodology of these effects has also been suggested using ambiguous figures. △ Less

Submitted 23 February, 2017; originally announced February 2017.

Comments: 11 pages, 2 figures

Non-classicality in mental states : an experimental study with ambiguous audio (music) stimuli

Authors: Souparno Roy, Ranjan Sengupta, Tarit Guhathakurata, Dipak Ghosh

Abstract: This paper attempts to address the question that whether the present physical or mathematical theories are sufficient to understand the complexities of human brain when it interacts with the external environment in the form of an auditory stimulus.There have been efforts reporting that the introduction of ambiguity in visual stimuli causes effects which can't be explained classically.In this paper… ▽ More This paper attempts to address the question that whether the present physical or mathematical theories are sufficient to understand the complexities of human brain when it interacts with the external environment in the form of an auditory stimulus.There have been efforts reporting that the introduction of ambiguity in visual stimuli causes effects which can't be explained classically.In this paper,it is investigated whether ambiguity in auditory stimuli can introduce any non-classical effects in human brain.Simple experiments were performed on normal subjects where they listened to an ambiguous auditory signal and responded to a question with 'yes' or 'no'.The outcome of the test showed that the classical formula of total probability does not hold true in this case.Results were interesting and indicate that there is a definite non-classicality in mental states in perception of ambiguous audio stimuli. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: Conference submission;comments are welcome

Temporal divergence in cropping pattern and its implications on geospatial drought assessment

Authors: C. S. Murthy, M. V. R. Sesha Sai, M. Naresh Kumar, P. S. Roy

Abstract: Time series data on cropping pattern at disaggregated level were analysed and its implications on geospatial drought assessment were demonstrated. An index of Cropping Pattern Dissimilarity (CP-DI) between a pair of years, developed in this study, proved that the cropping pattern of a year has a higher degree of similarity with that of recent past years only and tends to be dissimilar with longer… ▽ More Time series data on cropping pattern at disaggregated level were analysed and its implications on geospatial drought assessment were demonstrated. An index of Cropping Pattern Dissimilarity (CP-DI) between a pair of years, developed in this study, proved that the cropping pattern of a year has a higher degree of similarity with that of recent past years only and tends to be dissimilar with longer time difference. The temporal divergence in cropping pattern has direct implications on geospatial approach of drought assessment, in which, time series NDVI data are compared for drought interpretation. It was found that, seasonal NDVI profiles of drought year and normal year did not show any anomaly when the cropping patterns were dissimilar and two normal years having dissimilar cropping pattern showed different NDVI profiles. Therefore, it is suggested that such temporal comparisons of NDVI are better restricted to recent past years to achieve more objective interpretation. △ Less

Submitted 26 September, 2015; originally announced September 2015.

Comments: 9 pages, 14 figures

Journal ref: Geocarto International 10/2009, 24(5):377-395

A new hybrid spectral similarity measure for discrimination of Vigna species

Authors: M. Naresh Kumar, M. V. R Seshasai, K. S Vara Prasad, V. Kamala, K. V Ramana, R. S. Dwivedi, P. S. Roy

Abstract: The reflectance spectrum of the species in a hyperspectral data can be modelled as an n-dimensional vector. The spectral angle mapper computes the angle between the vectors which is used to discriminate the species. The spectral information divergence models the data as a probability distribution so that the spectral variability between the bands can be extracted using the stochastic measures. The… ▽ More The reflectance spectrum of the species in a hyperspectral data can be modelled as an n-dimensional vector. The spectral angle mapper computes the angle between the vectors which is used to discriminate the species. The spectral information divergence models the data as a probability distribution so that the spectral variability between the bands can be extracted using the stochastic measures. The hybrid approach of spectral angle mapper and spectral information divergence is found to be better discriminator than spectral angle mapper or spectral information divergence alone. The spectral correlation angle is computed as a cosine of the angle of the Pearsonian correlation coefficient between the vectors. The spectral correlation angle is a better measure than the spectral angle mapper as it considers only standardized values of the vectors rather than the absolute values of the vector. In the present paper a new hybrid measure is proposed which is based on the spectral correlation angle and the spectral information divergence. The proposed method has been compared with the hybrid approach of spectral information divergence and spectral angle mapper for discrimination of crops belonging to Vigna species using measures like relative spectral discriminatory power, relative discriminatory probability and relative discriminatory entropy in different spectral regions. Experimental results using the laboratory spectra show that the proposed method gives higher relative discriminatory power in 400nm-700nm spectral region. △ Less

Submitted 18 September, 2015; originally announced September 2015.

Journal ref: International Journal of Remote Sensing, 32(14),4041-4053, 2011

Using complex networks towards information retrieval and diagnostics in multidimensional imaging

Authors: Soumya Jyoti Banerjee, Mohammad Azharuddin, Debanjan Sen, Smruti Savale, Himadri Datta, Anjan Kr Dasgupta, Soumen Roy

Abstract: We present a fresh and broad yet simple approach towards information retrieval in general and diagnostics in particular by applying the theory of complex networks on multidimensional, dynamic images. We demonstrate a successful use of our method with the time series generated from high content thermal imaging videos of patients suffering from the aqueous deficient dry eye (ADDE) disease. Remarkabl… ▽ More We present a fresh and broad yet simple approach towards information retrieval in general and diagnostics in particular by applying the theory of complex networks on multidimensional, dynamic images. We demonstrate a successful use of our method with the time series generated from high content thermal imaging videos of patients suffering from the aqueous deficient dry eye (ADDE) disease. Remarkably, network analyses of thermal imaging time series of contact lens users and patients upon whom Laser-Assisted in situ Keratomileusis (Lasik) surgery has been conducted, exhibit pronounced similarity with results obtained from ADDE patients. We also propose a general framework for the transformation of multidimensional images to networks for futuristic biometry. Our approach is general and scalable to other fluctuation-based devices where network parameters derived from fluctuations, act as effective discriminators and diagnostic markers. △ Less

Submitted 1 December, 2015; v1 submitted 8 June, 2015; originally announced June 2015.

Comments: Replaced by published version. Detailed Supplementary Information on journal website

Journal ref: Scientific Reports, 5: 17271 (2015)

Control of human immune response function by T-cell population fluctuation and relaxation dynamics

Authors: Susmita Roy, Biman Bagchi

Abstract: Clinical studies have indicated that in malignant surveillances fluctuations in the population of certain effector T-cell repertoire become suppressed. Motivated by such observations and in an attempt to quantify adaptive human response to pathogens, we define an immune response function (IMRF) in terms of mean square fluctuations of T-cell concentrations. We employ a recently developed kinetic mo… ▽ More Clinical studies have indicated that in malignant surveillances fluctuations in the population of certain effector T-cell repertoire become suppressed. Motivated by such observations and in an attempt to quantify adaptive human response to pathogens, we define an immune response function (IMRF) in terms of mean square fluctuations of T-cell concentrations. We employ a recently developed kinetic model of T-cell regulation that contains the essential immunosuppressive effects of vitamin-D. We employ Gillespie algorithm to make the first study of fluctuations along the stochastic trajectories. This fluctuation-based IMRF can differentiate responses of different individuals after pathogenic incursion both under healthy and disease conditions. We find that relative fluctuations in T-cells (and hence IMRF) are different in strongly regulated (malignant prone) and weakly regulated (autoimmune prone) regions. The cross-over from one steady state (weakly regulated) to the other (strongly regulated) is accompanied by a divergence-like growth in the fluctuation of both the effector and regulatory T-cell concentration over a wide range of pathogenic stimulation, displaying a dynamical phase transition like behavior. The growth in fluctuation in this desired immune response regime is found to arise from an intermittent fluctuation between regulatory and effector T-cells that results in a bimodal distribution of population of each, indicating bistability. The signature of intermittent behavior is further confirmed by calculating the power spectrum of the corresponding fluctuation of time correlation function. The calculated time correlation functions of fluctuations show that the slow fluctuation causes the bistabilty in healthy state. Thus, in diseases diagnosis process, such steady state response parameters can provide immense information which might become helpful to define an immune status. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 6 Figures. arXiv admin note: text overlap with arXiv:1404.5111

Slow poisoning and destruction of networks: Edge proximity and its implications for biological and infrastructure networks

Authors: Soumya Jyoti Banerjee, Saptarshi Sinha, Soumen Roy

Abstract: We propose a network metric, edge proximity, ${\cal P}_e$, which demonstrates the importance of specific edges in a network, hitherto not captured by existing network metrics. The effects of removing edges with high ${\cal P}_e$ might initially seem inconspicuous but are eventually shown to be very harmful for networks. Compared to existing strategies, the removal of edges by ${\cal P}_e$ leads to… ▽ More We propose a network metric, edge proximity, ${\cal P}_e$, which demonstrates the importance of specific edges in a network, hitherto not captured by existing network metrics. The effects of removing edges with high ${\cal P}_e$ might initially seem inconspicuous but are eventually shown to be very harmful for networks. Compared to existing strategies, the removal of edges by ${\cal P}_e$ leads to a remarkable increase in the diameter and average shortest path length in undirected real and random networks till the first disconnection and well beyond. ${\cal P}_e$ can be consistently used to rupture the network into two nearly equal parts, thus presenting a very potent strategy to greatly harm a network. Targeting by ${\cal P}_e$ causes notable efficiency loss in U.S. and European power grid networks. ${\cal P}_e$ identifies proteins with essential cellular functions in protein-protein interaction networks. It pinpoints regulatory neural connections and important portions of the neural and brain networks, respectively. Energy flow interactions identified by ${\cal P}_e$ form the backbone of long food web chains. Finally, we scrutinize the potential of ${\cal P}_e$ in edge controllability dynamics of directed networks. △ Less

Submitted 11 February, 2015; v1 submitted 29 September, 2014; originally announced September 2014.

Comments: 5 Pages, 5 figures, Revtex-4.1

Journal ref: Phys. Rev. E 91, 022807 (2015)

Hydrophobic hydration driven self-assembly of Curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

Authors: Milan Kumar Hazra, Susmita Roy, Biman Bagchi

Abstract: As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of… ▽ More As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting strongly with hydrophilic groups at the surface of curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. △ Less

Submitted 20 June, 2014; originally announced June 2014.

Comments: 31 pages, 9 figures

A stochastic chemical dynamic approach to correlate autoimmunity and optimal vitamin-D range

Authors: Susmita Roy, Krishna Shrinivas, Biman Bagchi

Abstract: Motivated by several recent experimental observations that vitamin-D could interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells) to promote and to regulate different stages of immune response, we developed a coarse grained kinetic model in an attempt to quantify the role of vitamin-D in immunomodulatory responses. Our kinetic model, developed using the ideas of chemical ne… ▽ More Motivated by several recent experimental observations that vitamin-D could interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells) to promote and to regulate different stages of immune response, we developed a coarse grained kinetic model in an attempt to quantify the role of vitamin-D in immunomodulatory responses. Our kinetic model, developed using the ideas of chemical network theory, leads to a system of nine coupled equations that we solve both by direct and by stochastic (Gillespie) methods. Both the analyses consistently provide detail information on the dependence of immune response to the variation of critical rate parameters. We find that although vitamin-D plays a negligible role in the initial immune response, it exerts a profound influence in the long term, especially in helping the system to achieve a new, stable steady state. The study explores the role of vitamin-D in preserving an observed bistability in the phase diagram (spanned by system parameters) of immune regulation, thus allowing the response to tolerate a wide range of pathogenic stimulation which could help in resisting autoimmune diseases. We also study how vitamin-D affects the time dependent population of dendritic cells that connect between innate and adaptive immune responses. Variations in dose dependent response in anti-inflammatory and pro-inflammatory T-cell populations to vitamin-D correlate well with recent experimental results. Our kinetic model allows for an estimation of the range of optimum level of vitamin-D required for smooth functioning of the immune system and for control of both hyper-regulation and inflammation. Most importantly, the present study reveals that an overdose or toxic level of vitamin-D or any steroid analogue could give rise to too large a tolerant response, leading to an inefficacy in adaptive immune function. △ Less

Submitted 21 April, 2014; originally announced April 2014.

Comments: arXiv admin note: substantial text overlap with arXiv:1304.7193

A Comparative Study of Protein Unfolding in Aqueous Urea and DMSO Solutions: Surface Polarity, Solvent Specificity and Sequence of Secondary Structure Melting

Authors: Susmita Roy, Biman Bagchi

Abstract: Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here we alter the solvent environment in a controlled manner by choosing two different co-solvents of water, urea and dimethyl sulphoxide (DMSO), and study unfolding of four different proteins to understand the respective sequence of meltin… ▽ More Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here we alter the solvent environment in a controlled manner by choosing two different co-solvents of water, urea and dimethyl sulphoxide (DMSO), and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha-helices and beta-sheets in these two solvents. For example, at 8M urea solution, beta-sheet parts of a protein is found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8M DMSO solution unfolds alpha helices first, followed by the separation of beta-sheets for majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solution demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and non-polar residues of a protein and the mode of molecular action of a solvent on that protein. △ Less

Submitted 20 April, 2014; originally announced April 2014.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

Authors: Kushal Bagchi, Susmita Roy

Abstract: In addition to the biologically active monomer of the protein Insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The Insulin monomer contains two distinct surfaces, namely the dimer forming surface (DFS) and the hexamer forming surface (HFS) that are specifically designed to facilitate the formation of the dimer and the hexamer, re… ▽ More In addition to the biologically active monomer of the protein Insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The Insulin monomer contains two distinct surfaces, namely the dimer forming surface (DFS) and the hexamer forming surface (HFS) that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behaviour of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of Insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly a clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inverted orientation owing to the flat curvature of hydrophobic surface and interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces and constructed free energy surfaces of these water molecules.These free energy surfaces reveal the differences between the two hydration layers. △ Less

Submitted 21 December, 2013; originally announced December 2013.

Comments: 34 pages, 10 figures

Solvent Sensitivity of Protein Unfolding: Study of Chicken Villin Headpiece Subdomain in Water-Ethanol and Water-DMSO Mixtures

Authors: Rikhia Ghosh, Susmita Roy, Biman Bagchi

Abstract: In the present work we study and compare unfolding of a small protein, chicken villin headpiece (HP-36) in two different aqueous binary mixtures, namely water-ethanol (EtOH) and water-dimethyl sulphoxide (DMSO). In both the binary mixtures, HP-36 is found to unfold (fully or partially, depending on the mixture) under ambient conditions, that otherwise requires temperature as high as ~600 K to dena… ▽ More In the present work we study and compare unfolding of a small protein, chicken villin headpiece (HP-36) in two different aqueous binary mixtures, namely water-ethanol (EtOH) and water-dimethyl sulphoxide (DMSO). In both the binary mixtures, HP-36 is found to unfold (fully or partially, depending on the mixture) under ambient conditions, that otherwise requires temperature as high as ~600 K to denature in pure aqueous solvent. In all the cases, first step of unfolding is found to be similar, i.e. separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely Phe-7, Phe-11 and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. Subsequent unfolding steps follow different paths in different chemical environments. As both water-DMSO and water-ethanol show composition dependent anomalies, so do the details of unfolding dynamics. With an increase of co-solvent concentration different partially unfolded intermediates are found to be formed in both the cases. This is reflected in a remarkable non-monotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states is particularly attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol and methyl groups of DMSO. While in DMSO the protein gradually attains a completely unfolded state at xDMSO=0.30, unfolding in water-ethanol appears to be more complex and sensitive to solvent composition. △ Less

Submitted 10 July, 2013; originally announced July 2013.

Turtle: Identifying frequent k-mers with cache-efficient algorithms

Authors: Rajat Shuvro Roy, Debashish Bhattacharya, Alexander Schliep

Abstract: Counting the frequencies of k-mers in read libraries is often a first step in the analysis of high-throughput sequencing experiments. Infrequent k-mers are assumed to be a result of sequencing errors. The frequent k-mers constitute a reduced but error-free representation of the experiment, which can inform read error correction or serve as the input to de novo assembly methods. Ideally, the memory… ▽ More Counting the frequencies of k-mers in read libraries is often a first step in the analysis of high-throughput sequencing experiments. Infrequent k-mers are assumed to be a result of sequencing errors. The frequent k-mers constitute a reduced but error-free representation of the experiment, which can inform read error correction or serve as the input to de novo assembly methods. Ideally, the memory requirement for counting should be linear in the number of frequent k-mers and not in the, typically much larger, total number of k-mers in the read library. We present a novel method that balances time, space and accuracy requirements to efficiently extract frequent k-mers even for high coverage libraries and large genomes such as human. Our method is designed to minimize cache-misses in a cache-efficient manner by using a Pattern-blocked Bloom filter to remove infrequent k-mers from consideration in combination with a novel sort-and-compact scheme, instead of a Hash, for the actual counting. While this increases theoretical complexity, the savings in cache misses reduce the empirical running times. A variant can resort to a counting Bloom filter for even larger savings in memory at the expense of false negatives in addition to the false positives common to all Bloom filter based approaches. A comparison to the state-of-the-art shows reduced memory requirements and running times. Note that we also provide the first competitive method to count k-mers up to size 64. △ Less

Submitted 8 May, 2013; originally announced May 2013.

Vitamin D sensitivity to the immune responses and autoimmunity: A chemical network model study

Authors: Susmita Roy, Krishna Shrinivas, Biman Bagchi

Abstract: Although Vitamin D is believed to be involved in a large number of immune responses our understanding of these processes at cellular level remained at infancy. We develop and solve a coarse grained kinetic network model to quantify the effects of variation of vitamin D on human immunity. The system of equations accounts for known inter-relation between active and inactive vitamin D, antigen presen… ▽ More Although Vitamin D is believed to be involved in a large number of immune responses our understanding of these processes at cellular level remained at infancy. We develop and solve a coarse grained kinetic network model to quantify the effects of variation of vitamin D on human immunity. The system of equations accounts for known inter-relation between active and inactive vitamin D, antigen presenting cells, effector T cells, regulatory T cells and pathogen. Both time dependent and steady state solutions are obtained. The time dependent solution of the system of equations reveals that the immune response is rather strongly regulated in presence of vitamin D. We found quantitatively that lower than optimum levels of concentration of active vitamin D correspond to weak regulation where, once a pathogen/antigen enters the body, the nature of the immune response would be less regulatory and hence more aggressive, or inflammatory. The steady state solution of our model shows that vitamin D enhances the tolerance level of immune system, thereby increasing resistance to autoimmune diseases. Our model and accompanied numerical analyses reveal another important aspect of immunity: While extremely low levels of vitamin D could lead to increased risk of autoimmune responses, an overdose (toxic) level would give rise to too large a tolerant response, leading to the increased risk of tumors and cancerous cell growth. △ Less

Submitted 26 April, 2013; originally announced April 2013.

Comments: 44 pages, 9 figures

An Efficient Biological Sequence Compression Technique Using LUT And Repeat In The Sequence

Authors: Subhankar Roy, Sunirmal Khatua, Sudipta Roy, Samir K. Bandyopadhyay

Abstract: Data compression plays an important role to deal with high volumes of DNA sequences in the field of Bioinformatics. Again data compression techniques directly affect the alignment of DNA sequences. So the time needed to decompress a compressed sequence has to be given equal priorities as with compression ratio. This article contains first introduction then a brief review of different biological se… ▽ More Data compression plays an important role to deal with high volumes of DNA sequences in the field of Bioinformatics. Again data compression techniques directly affect the alignment of DNA sequences. So the time needed to decompress a compressed sequence has to be given equal priorities as with compression ratio. This article contains first introduction then a brief review of different biological sequence compression after that my proposed work then our two improved Biological sequence compression algorithms after that result followed by conclusion and discussion, future scope and finally references. These algorithms gain a very good compression factor with higher saving percentage and less time for compression and decompression than the previous Biological Sequence compression algorithms. Keywords: Hash map table, Tandem repeats, compression factor, compression time, saving percentage, compression, decompression process. △ Less

Submitted 11 November, 2012; v1 submitted 26 September, 2012; originally announced September 2012.

Comments: 9 pages, 3 figures, 5 tables

Key to Network Controllability

Authors: Soumya Jyoti Banerjee, Soumen Roy

Abstract: Liu et al recently proposed a minimum number of driver nodes needed to obtain full structural controllability over a directed network. Driver nodes are unmatched nodes, from which there are directed paths to all matched nodes. Their most important assertion is that a system's controllability is to a great extent encoded by the underlying network's degree distribution, $P(k_{in}, k_{out})$. Is the… ▽ More Liu et al recently proposed a minimum number of driver nodes needed to obtain full structural controllability over a directed network. Driver nodes are unmatched nodes, from which there are directed paths to all matched nodes. Their most important assertion is that a system's controllability is to a great extent encoded by the underlying network's degree distribution, $P(k_{in}, k_{out})$. Is the controllability of a network decided almost completely by the immediate neighbourhood of a node, while, even slightly distant nodes play no role at all? Motivated by the above question, in this communication, we argue that an effective understanding of controllability in directed networks can be reached using distance based measures of closeness centrality and betweenness centrality and may not require the knowledge of local connectivity measures like in-degree and out-degree. △ Less

Submitted 17 September, 2012; originally announced September 2012.

Comments: Arising from YY Liu, JJ Slotine and AL Barabasi, Nature 473, 167 - 173 (2011)

Systems biology beyond degree, hubs and scale-free networks: the case for multiple metrics in complex networks

Authors: Soumen Roy

Abstract: Modeling and topological analysis of networks in biological and other complex systems, must venture beyond the limited consideration of very few network metrics like degree, betweenness or assortativity. A proper identification of informative and redundant entities from many different metrics, using recently demonstrated techniques, is essential. A holistic comparison of networks and growth models… ▽ More Modeling and topological analysis of networks in biological and other complex systems, must venture beyond the limited consideration of very few network metrics like degree, betweenness or assortativity. A proper identification of informative and redundant entities from many different metrics, using recently demonstrated techniques, is essential. A holistic comparison of networks and growth models is best achieved only with the use of such methods. △ Less

Submitted 3 June, 2012; v1 submitted 10 May, 2012; originally announced May 2012.

Comments: To appear in Systems and Synthetic Biology (Springer)

Journal ref: Systems and Synthetic Biology (2012) 6:31-34

An Information Theoretical Analysis of Kinase Activated Phosphorylation Dephosphorylation Cycle

Authors: Hong Qian, Sumit Roy

Abstract: Signal transduction, the information processing mechanism in biological cells, is carried out by a network of biochemical reactions. The dynamics of driven biochemical reactions can be studied in terms of nonequilibrium statistical physics. Such systems may also be studied in terms of Shannon's information theory. We combine these two perspectives in this study of the basic units (modules) of cell… ▽ More Signal transduction, the information processing mechanism in biological cells, is carried out by a network of biochemical reactions. The dynamics of driven biochemical reactions can be studied in terms of nonequilibrium statistical physics. Such systems may also be studied in terms of Shannon's information theory. We combine these two perspectives in this study of the basic units (modules) of cellular signaling: the phosphorylation dephosphorylation cycle (PdPC) and the guanosine triphosphatase (GTPase). We show that the channel capacity is zero if and only if the free energy expenditure of biochemical system is zero. In fact, a positive correlation between the channel capacity and free energy expenditure is observed. In terms of the information theory, a linear signaling cascade consisting of multiple steps of PdPC can function as a distributed "multistage code". With increasing number of steps in the cascade, the system trades channel capacity with the code complexity. Our analysis shows that while a static code can be molecular structural based; a biochemical communication channel has to have energy expenditure. △ Less

Submitted 20 December, 2011; originally announced December 2011.

Comments: 17 pages, 7 figures

Journal ref: IEEE Transactions on NanoBioscience, Volume 11, Issue 3, pp. 289 - 295 (2012)

SLIQ: Simple Linear Inequalities for Efficient Contig Scaffolding

Authors: Rajat S. Roy, Kevin C. Chen, Anirvan M. Sengupta, Alexander Schliep

Abstract: Scaffolding is an important subproblem in "de novo" genome assembly in which mate pair data are used to construct a linear sequence of contigs separated by gaps. Here we present SLIQ, a set of simple linear inequalities derived from the geometry of contigs on the line that can be used to predict the relative positions and orientations of contigs from individual mate pair reads and thus produce a c… ▽ More Scaffolding is an important subproblem in "de novo" genome assembly in which mate pair data are used to construct a linear sequence of contigs separated by gaps. Here we present SLIQ, a set of simple linear inequalities derived from the geometry of contigs on the line that can be used to predict the relative positions and orientations of contigs from individual mate pair reads and thus produce a contig digraph. The SLIQ inequalities can also filter out unreliable mate pairs and can be used as a preprocessing step for any scaffolding algorithm. We tested the SLIQ inequalities on five real data sets ranging in complexity from simple bacterial genomes to complex mammalian genomes and compared the results to the majority voting procedure used by many other scaffolding algorithms. SLIQ predicted the relative positions and orientations of the contigs with high accuracy in all cases and gave more accurate position predictions than majority voting for complex genomes, in particular the human genome. Finally, we present a simple scaffolding algorithm that produces linear scaffolds given a contig digraph. We show that our algorithm is very efficient compared to other scaffolding algorithms while maintaining high accuracy in predicting both contig positions and orientations for real data sets. △ Less

Submitted 9 November, 2011; v1 submitted 6 November, 2011; originally announced November 2011.

Comments: 16 pages, 6 figures, 7 tables

Relative velocity of sliding of microtubules by the action of Kinesin-5

Authors: Sthitadhi Roy

Abstract: Kinesin-5, also known as Eg5 in vertebrates is a processive motor with 4 heads, which moves on filamentous tracks called microtubules. The basic function of Kinesin-5 is to slide apart two anti-parallel microtubules by simultaneously walking on both the microtubules. We develop an analytical expression for the steady-state relative velocity of this sliding in terms of the rates of attachments and… ▽ More Kinesin-5, also known as Eg5 in vertebrates is a processive motor with 4 heads, which moves on filamentous tracks called microtubules. The basic function of Kinesin-5 is to slide apart two anti-parallel microtubules by simultaneously walking on both the microtubules. We develop an analytical expression for the steady-state relative velocity of this sliding in terms of the rates of attachments and detachments of motor heads with the ATPase sites on the microtubules. We first analyse the motion of one pair of motor heads on one microtubule and then couple it to the motion of the other pair of motor heads of the same motor on the second microtubule to get the relative velocity of sliding. △ Less

Submitted 10 March, 2011; originally announced March 2011.

Comments: 6 pages Revtex, 5 eps figures

Strong associations between microbe phenotypes and their network architecture

Authors: Soumen Roy, Vladimir Filkov

Abstract: Understanding the dependence and interplay between architecture and function in biological networks has great relevance to disease progression, biological fabrication and biological systems in general. We propose methods to assess the association of various microbe characteristics and phenotypes with the topology of their networks. We adopt an automated approach to characterize metabolic network… ▽ More Understanding the dependence and interplay between architecture and function in biological networks has great relevance to disease progression, biological fabrication and biological systems in general. We propose methods to assess the association of various microbe characteristics and phenotypes with the topology of their networks. We adopt an automated approach to characterize metabolic networks of 32 microbial species using 11 topological metrics from complex networks. Clustering allows us to extract the indispensable, independent and informative metrics. Using hierarchical linear modeling, we identify relevant subgroups of these metrics and establish that they associate with microbial phenotypes surprisingly well. This work can serve as a stepping stone to cataloging biologically relevant topological properties of networks and towards better modeling of phenotypes. The methods we use can also be applied to networks from other disciplines. △ Less

Submitted 7 October, 2009; v1 submitted 7 January, 2009; originally announced January 2009.

Comments: Replaced by the version scheduled to appear in Phys. Rev. E (Rapid Comm.)

Journal ref: Phys. Rev. E 80, 040902(R) (2009) (Rapid Comm.)

Modeling and verifying a broad array of network properties

Authors: Vladimir Filkov, Zachary M. Saul, Soumen Roy, Raissa M. D'Souza, Premkumar T. Devanbu

Abstract: Motivated by widely observed examples in nature, society and software, where groups of already related nodes arrive together and attach to an existing network, we consider network growth via sequential attachment of linked node groups, or graphlets. We analyze the simplest case, attachment of the three node V-graphlet, where, with probability alpha, we attach a peripheral node of the graphlet, a… ▽ More Motivated by widely observed examples in nature, society and software, where groups of already related nodes arrive together and attach to an existing network, we consider network growth via sequential attachment of linked node groups, or graphlets. We analyze the simplest case, attachment of the three node V-graphlet, where, with probability alpha, we attach a peripheral node of the graphlet, and with probability (1-alpha), we attach the central node. Our analytical results and simulations show that tuning alpha produces a wide range in degree distribution and degree assortativity, achieving assortativity values that capture a diverse set of many real-world systems. We introduce a fifteen-dimensional attribute vector derived from seven well-known network properties, which enables comprehensive comparison between any two networks. Principal Component Analysis (PCA) of this attribute vector space shows a significantly larger coverage potential of real-world network properties by a simple extension of the above model when compared against a classic model of network growth. △ Less

Submitted 26 March, 2009; v1 submitted 12 May, 2008; originally announced May 2008.

Comments: To appear in Europhysics Letters

Journal ref: Europhysics Letters, 86 (2009) 28003

Non Markovian Noise mediated through Anamolous Diffusion within Ion Channels

Authors: Sisir Roy, Indranil Mitra, Rodolfo Llinas

Abstract: It is quite clear from a wide range of experiments that gating phenomena of ion channels is inherently stochastic. It has been discussed using BD simulations in a recent paper that memory effects in ion transport is negligible, unless the barrier height is high. In this brief report we like to state using Differential Stochastic Methods (DSM's) that the Markovian property of exponential dwell ti… ▽ More It is quite clear from a wide range of experiments that gating phenomena of ion channels is inherently stochastic. It has been discussed using BD simulations in a recent paper that memory effects in ion transport is negligible, unless the barrier height is high. In this brief report we like to state using Differential Stochastic Methods (DSM's) that the Markovian property of exponential dwell times do indeed give rise to a high barrier, which in turn indicates that memory effects need not be ignored. We have thus constructed a Generalized Langevin Equation which contains a combination of Non Markovian at different time scales & Markovian processes and develop an algorithm to describe the scheme of events. We see that the oscillatory function behaviour with exponential decay is obtained in the Markovian limit and two distinct time scales corresponding to the processes of diffusion & drift may be obtained from preliminary simulation results. We propose that the results need much more inspection and it will be worthwhile to reproduce using MD simulations. The most important idea which we like to propose in this paper is that the rise of time scales and memory effects may be inherently related to the differential behaviour of shear viscosity in the cytoplasm & extracellular matrix. △ Less

Submitted 23 January, 2008; originally announced January 2008.

Comments: 10 pages, 6 figures

Co-operativity in neurons and the role of noise in brain

Authors: Indranil Mitra, Sisir Roy, Gary Hastings

Abstract: In view of some recent results in case of the dopaminergic neurons exhibiting long range correlations in VTA of the limbic brain we are interested to find out whether any stochastic nonlinear response may be reproducible in the nano scales usimg the results of quantum mechanics. We have developed a scheme to investigate this situation in this paper by taking into consideration the Schrodinger eq… ▽ More In view of some recent results in case of the dopaminergic neurons exhibiting long range correlations in VTA of the limbic brain we are interested to find out whether any stochastic nonlinear response may be reproducible in the nano scales usimg the results of quantum mechanics. We have developed a scheme to investigate this situation in this paper by taking into consideration the Schrodinger equation (SE) in an arbitrary manifold with a metric, which is in some sense a special case of the heat kernel equation. The special case of this heat kernel equation is the diffusion equation, which may reproduce some key phenomena of the neural activities. We make a dual equivalent circuit model of SE and incorporate non commutativity and noise inside the circuit scheme. The behaviour of the circuit elements with interesting limits are investigated. The most bizarre part is the long range response of the model by dint of the Central Limit Theorem, which is responsible for coherent behaviour of a large assembly of neurons. △ Less

Submitted 15 August, 2007; originally announced August 2007.

Comments: 12 Pages, Updated Version, SPIE 2007 Conference