-
Insights into Protein Unfolding under pH, Temperature, and Shear using Molecular Dynamics Simulations
Authors:
Yinhao Jia,
Clare Cocker,
Janani Sampath
Abstract:
Protein biologics hold immense potential in therapeutic applications, but their ephemeral nature has hindered their widespread application. The effects of different stressors on protein folding have long been studied, but whether these stressors induce protein unfolding through different pathways remains unclear. In this work, we conduct all-atom molecular dynamics simulations to investigate the u…
▽ More
Protein biologics hold immense potential in therapeutic applications, but their ephemeral nature has hindered their widespread application. The effects of different stressors on protein folding have long been studied, but whether these stressors induce protein unfolding through different pathways remains unclear. In this work, we conduct all-atom molecular dynamics simulations to investigate the unfolding of bovine serum albumin (BSA) under three distinct external stressors: high temperature, acidic pH, and shear stress. Our findings reveal that each stressor induces unique unfolding patterns in BSA, indicating stressor-specific unfolding pathways. Detailed structural analysis showed that high temperature significantly disrupts the protein's secondary structure, while acidic pH causes notable alterations in the tertiary structure, leading to domain separation and an extended shape. Shear stress initially perturbs the tertiary structure, initiating structural rearrangements followed by a loss of secondary structure. These distinct unfolding behaviors suggest that different stabilization strategies are required to enhance protein stability under various denaturation conditions. Insights from these unfolding studies can inform the design of materials, especially polymers, aimed at improving protein stability.
△ Less
Submitted 20 August, 2024;
originally announced August 2024.
-
SIU: A Million-Scale Structural Small Molecule-Protein Interaction Dataset for Unbiased Bioactivity Prediction
Authors:
Yanwen Huang,
Bowen Gao,
Yinjun Jia,
Hongbo Ma,
Wei-Ying Ma,
Ya-Qin Zhang,
Yanyan Lan
Abstract:
Small molecules play a pivotal role in modern medicine, and scrutinizing their interactions with protein targets is essential for the discovery and development of novel, life-saving therapeutics. The term "bioactivity" encompasses various biological effects resulting from these interactions, including both binding and functional responses. The magnitude of bioactivity dictates the therapeutic or t…
▽ More
Small molecules play a pivotal role in modern medicine, and scrutinizing their interactions with protein targets is essential for the discovery and development of novel, life-saving therapeutics. The term "bioactivity" encompasses various biological effects resulting from these interactions, including both binding and functional responses. The magnitude of bioactivity dictates the therapeutic or toxic pharmacological outcomes of small molecules, rendering accurate bioactivity prediction crucial for the development of safe and effective drugs. However, existing structural datasets of small molecule-protein interactions are often limited in scale and lack systematically organized bioactivity labels, thereby impeding our understanding of these interactions and precise bioactivity prediction. In this study, we introduce a comprehensive dataset of small molecule-protein interactions, consisting of over a million binding structures, each annotated with real biological activity labels. This dataset is designed to facilitate unbiased bioactivity prediction. We evaluated several classical models on this dataset, and the results demonstrate that the task of unbiased bioactivity prediction is challenging yet essential.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Multi-level Interaction Modeling for Protein Mutational Effect Prediction
Authors:
Yuanle Mo,
Xin Hong,
Bowen Gao,
Yinjun Jia,
Yanyan Lan
Abstract:
Protein-protein interactions are central mediators in many biological processes. Accurately predicting the effects of mutations on interactions is crucial for guiding the modulation of these interactions, thereby playing a significant role in therapeutic development and drug discovery. Mutations generally affect interactions hierarchically across three levels: mutated residues exhibit different si…
▽ More
Protein-protein interactions are central mediators in many biological processes. Accurately predicting the effects of mutations on interactions is crucial for guiding the modulation of these interactions, thereby playing a significant role in therapeutic development and drug discovery. Mutations generally affect interactions hierarchically across three levels: mutated residues exhibit different sidechain conformations, which lead to changes in the backbone conformation, eventually affecting the binding affinity between proteins. However, existing methods typically focus only on sidechain-level interaction modeling, resulting in suboptimal predictions. In this work, we propose a self-supervised multi-level pre-training framework, ProMIM, to fully capture all three levels of interactions with well-designed pretraining objectives. Experiments show ProMIM outperforms all the baselines on the standard benchmark, especially on mutations where significant changes in backbone conformations may occur. In addition, leading results from zero-shot evaluations for SARS-CoV-2 mutational effect prediction and antibody optimization underscore the potential of ProMIM as a powerful next-generation tool for developing novel therapeutic approaches and new drugs.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Full-Atom Peptide Design with Geometric Latent Diffusion
Authors:
Xiangzhe Kong,
Yinjun Jia,
Wenbing Huang,
Yang Liu
Abstract:
Peptide design plays a pivotal role in therapeutics, allowing brand new possibility to leverage target binding sites that are previously undruggable. Most existing methods are either inefficient or only concerned with the target-agnostic design of 1D sequences. In this paper, we propose a generative model for full-atom \textbf{Pep}tide design with \textbf{G}eometric \textbf{LA}tent \textbf{D}iffus…
▽ More
Peptide design plays a pivotal role in therapeutics, allowing brand new possibility to leverage target binding sites that are previously undruggable. Most existing methods are either inefficient or only concerned with the target-agnostic design of 1D sequences. In this paper, we propose a generative model for full-atom \textbf{Pep}tide design with \textbf{G}eometric \textbf{LA}tent \textbf{D}iffusion (PepGLAD). We first establish a benchmark consisting of both 1D sequences and 3D structures from Protein Data Bank (PDB) and literature for systematic evaluation. We then identify two major challenges of leveraging current diffusion-based models for peptide design: the full-atom geometry and the variable binding geometry. To tackle the first challenge, PepGLAD derives a variational autoencoder that first encodes full-atom residues of variable size into fixed-dimensional latent representations, and then decodes back to the residue space after conducting the diffusion process in the latent space. For the second issue, PepGLAD explores a receptor-specific affine transformation to convert the 3D coordinates into a shared standard space, enabling better generalization ability across different binding shapes. Experimental Results show that our method not only improves diversity and binding affinity significantly in the task of sequence-structure co-design, but also excels at recovering reference structures for binding conformation generation.
△ Less
Submitted 21 May, 2024; v1 submitted 21 February, 2024;
originally announced February 2024.
-
Protein-ligand binding representation learning from fine-grained interactions
Authors:
Shikun Feng,
Minghao Li,
Yinjun Jia,
Weiying Ma,
Yanyan Lan
Abstract:
The binding between proteins and ligands plays a crucial role in the realm of drug discovery. Previous deep learning approaches have shown promising results over traditional computationally intensive methods, but resulting in poor generalization due to limited supervised data. In this paper, we propose to learn protein-ligand binding representation in a self-supervised learning manner. Different f…
▽ More
The binding between proteins and ligands plays a crucial role in the realm of drug discovery. Previous deep learning approaches have shown promising results over traditional computationally intensive methods, but resulting in poor generalization due to limited supervised data. In this paper, we propose to learn protein-ligand binding representation in a self-supervised learning manner. Different from existing pre-training approaches which treat proteins and ligands individually, we emphasize to discern the intricate binding patterns from fine-grained interactions. Specifically, this self-supervised learning problem is formulated as a prediction of the conclusive binding complex structure given a pocket and ligand with a Transformer based interaction module, which naturally emulates the binding process. To ensure the representation of rich binding information, we introduce two pre-training tasks, i.e.~atomic pairwise distance map prediction and mask ligand reconstruction, which comprehensively model the fine-grained interactions from both structure and feature space. Extensive experiments have demonstrated the superiority of our method across various binding tasks, including protein-ligand affinity prediction, virtual screening and protein-ligand docking.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
Morphological entropy encodes cellular migration strategies on multiple length scales
Authors:
Yanping Liu,
Yang Jiao,
Qihui Fan,
Xinwei Li,
Zhichao Liu,
Jun Hu,
Jianwei Shuai,
Liyu Liu,
Zhangyong Li
Abstract:
Cell migration is crucial to many physiological and pathological processes. During migration, a cell adapts its morphology, including the overall morphology and nucleus morphology, in response to various cues in complex microenvironments, e.g. topotaxis and chemotaxis. Thus, cellular morphology dynamics can encode migration strategies based on which various migration mechanisms can be inferred. Ho…
▽ More
Cell migration is crucial to many physiological and pathological processes. During migration, a cell adapts its morphology, including the overall morphology and nucleus morphology, in response to various cues in complex microenvironments, e.g. topotaxis and chemotaxis. Thus, cellular morphology dynamics can encode migration strategies based on which various migration mechanisms can be inferred. However, how to decipher cell migration mechanisms encoded in the morphology dynamics remains a challenging problem. Here we introduce a novel universal metric, namely cell morphological entropy (CME), by combining parametric morphological analysis with Shannon entropy. The utility of CME, which accurately quantifies the complex cellular morphology on multiple length scales through the deviation from the perfect circular shape, is demonstrated using a variety of normal and tumorous cell lines in distinct in vitro microenvironments. Our results reveal that 1) the effects of geometric constraints on cell nucleus, 2) the emerging interplays of MCF-10A cells migrating on collagen gel, and 3) the critical transition of tumor spheroid from proliferation to invasion. The analysis indicates that the CME offers a physically interpretable and efficient tool to quantify morphology on multiple length scales in real-time, which provides more insights into cell migration, and further contributing to the understanding of the diverse behavioral modes as well as collective cell motility in more complex microenvironment.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
Machine Learning-guided Lipid Nanoparticle Design for mRNA Delivery
Authors:
Daisy Yi Ding,
Yuhui Zhang,
Yuan Jia,
Jiuzhi Sun
Abstract:
While RNA technologies hold immense therapeutic potential in a range of applications from vaccination to gene editing, the broad implementation of these technologies is hindered by the challenge of delivering these agents effectively. Lipid nanoparticles have emerged as one of the most widely used delivery agents, but their design optimization relies on laborious and costly experimental methods. W…
▽ More
While RNA technologies hold immense therapeutic potential in a range of applications from vaccination to gene editing, the broad implementation of these technologies is hindered by the challenge of delivering these agents effectively. Lipid nanoparticles have emerged as one of the most widely used delivery agents, but their design optimization relies on laborious and costly experimental methods. We propose to in silico optimize LNP design with machine learning models. On a curated dataset of 622 LNPs from published studies, we demonstrate the effectiveness of our model in predicting the transfection efficiency of unseen LNPs, with the multilayer perceptron achieving a classification accuracy of 98% on the test set. Our work represents a pioneering effort in combining ML and LNP design, offering significant potential for improving screening efficiency by computationally prioritizing LNP candidates for experimental validation and accelerating the development of effective mRNA delivery systems.
△ Less
Submitted 28 August, 2023; v1 submitted 2 August, 2023;
originally announced August 2023.
-
Perfused and Nonperfused Microaneurysms Identified and Characterized by Structural and Angiographic OCT
Authors:
Min Gao,
Tristan T. Hormel,
Yukun Guo,
Kotaro Tsuboi,
Christina J. Flaxel,
David Huang,
Thomas S. Hwang,
Yali Jia
Abstract:
Purpose: Microaneurysms (MAs) have distinct, oval-shaped, hyperreflective walls on structural OCT, and inconsistent flow signal in the lumen with OCT angiography (OCTA). Their relationship to regional macular edema in diabetic retinopathy (DR) has not been quantitatively explored. Participants: A total of 99 participants, including 23 with mild, NPDR, 25 with moderate NPDR, 34 with severe NPDR, an…
▽ More
Purpose: Microaneurysms (MAs) have distinct, oval-shaped, hyperreflective walls on structural OCT, and inconsistent flow signal in the lumen with OCT angiography (OCTA). Their relationship to regional macular edema in diabetic retinopathy (DR) has not been quantitatively explored. Participants: A total of 99 participants, including 23 with mild, NPDR, 25 with moderate NPDR, 34 with severe NPDR, and 17 with proliferative DR. Methods: We obtained 3 x 3-mm scans with a commercial device (Solix, Visionix/Optovue) in 99 patients with DR. Trained graders manually identified MAs and their location relative to the anatomic layers from cross-sectional OCT. Microaneurysms were first classified as perfused if flow signal was present in the OCTA channel. Then, perfused MAs were further classified into fully and partially perfused MAs based on the flow characteristics in en face OCTA. The presence of retinal fluid based on OCT near MAs was compared between perfused and nonperfused types. We also compared OCT-based MA detection to fundus photography (FP)- and fluorescein angiography (FA)-based detection. Results: We identified 308 MAs (166 fully perfused, 88 partially perfused, 54 nonperfused) in 42 eyes using OCT and OCTA. Nearly half of the MAs identified in this study straddle the inner nuclear layer and outer plexiform layer. Compared with partially perfused and nonperfused MAs, fully perfused MAs were more likely to be associated with local retinal fluid. The associated fluid volumes were larger with fully perfused MAs compared with other types. OCT/OCTA detected all MAs found on FP. Although not all MAs seen with FA were identified with OCT, some MAs seen with OCT were not visible with FA or FP. Conclusions: OCT-identified MAs with colocalized flow on OCTA are more likely to be associated with DME than those without flow.
△ Less
Submitted 9 October, 2023; v1 submitted 23 March, 2023;
originally announced March 2023.
-
Evaluating Evasion Strategies in Zebrafish Larvae
Authors:
Yusheng Jiao,
Brendan Colvert,
Yi Man,
Matthew J. McHenry,
Eva Kanso
Abstract:
An effective evasion strategy allows prey to survive encounters with predators. Prey are generally thought to escape in a direction that is either random or serves to maximize the minimum distance from the predator. Here we introduce a comprehensive approach to determine the most likely evasion strategy among multiple hypotheses and the role of biomechanical constraints on the escape response of p…
▽ More
An effective evasion strategy allows prey to survive encounters with predators. Prey are generally thought to escape in a direction that is either random or serves to maximize the minimum distance from the predator. Here we introduce a comprehensive approach to determine the most likely evasion strategy among multiple hypotheses and the role of biomechanical constraints on the escape response of prey fish. Through a consideration of six strategies with sensorimotor noise and previous kinematic measurements, our analysis shows that zebrafish larvae generally escape in a direction orthogonal to the predator's heading. By sensing only the predator's heading, this orthogonal strategy maximizes the distance from fast-moving predators, and, when operating within the biomechanical constraints of the escape response, it provides the best predictions of prey behavior among all alternatives. This work demonstrates a framework for resolving the strategic basis of evastion in predator-prey interactions, which could be applied to a broad diversity of animals.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
Computational interference microscopy enabled by deep learning
Authors:
Yuheng Jiao,
Yuchen R. He,
Mikhail E. Kandel,
Xiaojun Liu,
Wenlong Lu,
Gabriel Popescu
Abstract:
Quantitative phase imaging (QPI) has been widely applied in characterizing cells and tissues. Spatial light interference microscopy (SLIM) is a highly sensitive QPI method, due to its partially coherent illumination and common path interferometry geometry. However, its acquisition rate is limited because of the four-frame phase-shifting scheme. On the other hand, off-axis methods like diffraction…
▽ More
Quantitative phase imaging (QPI) has been widely applied in characterizing cells and tissues. Spatial light interference microscopy (SLIM) is a highly sensitive QPI method, due to its partially coherent illumination and common path interferometry geometry. However, its acquisition rate is limited because of the four-frame phase-shifting scheme. On the other hand, off-axis methods like diffraction phase microscopy (DPM), allows for single-shot QPI. However, the laser-based DPM system is plagued by spatial noise due to speckles and multiple reflections. In a parallel development, deep learning was proven valuable in the field of bioimaging, especially due to its ability to translate one form of contrast into another. Here, we propose using deep learning to produce synthetic, SLIM-quality, high-sensitivity phase maps from DPM, single-shot images as input. We used an inverted microscope with its two ports connected to the DPM and SLIM modules, such that we have access to the two types of images on the same field of view. We constructed a deep learning model based on U-net and trained on over 1,000 pairs of DPM and SLIM images. The model learned to remove the speckles in laser DPM and overcame the background phase noise in both the test set and new data. Furthermore, we implemented the neural network inference into the live acquisition software, which now allows a DPM user to observe in real-time an extremely low-noise phase image. We demonstrated this principle of computational interference microscopy (CIM) imaging using blood smears, as they contain both erythrocytes and leukocytes, in static and dynamic conditions.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
Learning to swim in potential flow
Authors:
Yusheng Jiao,
Feng Ling,
Sina Heydari,
Nicolas Heess,
Josh Merel,
Eva Kanso
Abstract:
Fish swim by undulating their bodies. These propulsive motions require coordinated shape changes of a body that interacts with its fluid environment, but the specific shape coordination that leads to robust turning and swimming motions remains unclear. To address the problem of underwater motion planning, we propose a simple model of a three-link fish swimming in a potential flow environment and w…
▽ More
Fish swim by undulating their bodies. These propulsive motions require coordinated shape changes of a body that interacts with its fluid environment, but the specific shape coordination that leads to robust turning and swimming motions remains unclear. To address the problem of underwater motion planning, we propose a simple model of a three-link fish swimming in a potential flow environment and we use model-free reinforcement learning for shape control. We arrive at optimal shape changes for two swimming tasks: swimming in a desired direction and swimming towards a known target. This fish model belongs to a class of problems in geometric mechanics, known as driftless dynamical systems, which allow us to analyze the swimming behavior in terms of geometric phases over the shape space of the fish. These geometric methods are less intuitive in the presence of drift. Here, we use the shape space analysis as a tool for assessing, visualizing, and interpreting the control policies obtained via reinforcement learning in the absence of drift. We then examine the robustness of these policies to drift-related perturbations. Although the fish has no direct control over the drift itself, it learns to take advantage of the presence of moderate drift to reach its target.
△ Less
Submitted 7 December, 2020; v1 submitted 30 September, 2020;
originally announced September 2020.
-
Focal Loss Analysis of Nerve Fiber Layer Reflectance for Glaucoma Diagnosis
Authors:
Ou Tan,
Liang Liu,
Qisheng You,
Jie Wang,
Aiyin Chen,
Eliesa Ing,
John C. Morrison,
Yali Jia,
David Huang
Abstract:
Purpose: To evaluate nerve fiber layer (NFL) reflectance for glaucoma diagnosis. Methods: Participants were imaged with 4.5X4.5-mm volumetric disc scans using spectral-domain optical coherence tomography (OCT). The normalized NFL reflectance map was processed by an azimuthal filter to reduce directional reflectance bias due to variation of beam incidence angle. The peripapillary area of the map wa…
▽ More
Purpose: To evaluate nerve fiber layer (NFL) reflectance for glaucoma diagnosis. Methods: Participants were imaged with 4.5X4.5-mm volumetric disc scans using spectral-domain optical coherence tomography (OCT). The normalized NFL reflectance map was processed by an azimuthal filter to reduce directional reflectance bias due to variation of beam incidence angle. The peripapillary area of the map was divided into 160 superpixels. Average reflectance was the mean of superpixel reflectance. Low-reflectance superpixels were identified as those with NFL reflectance below the 5 percentile normative cutoff. Focal reflectance loss was measure by summing loss in low-reflectance superpixels. Results: Thirty-five normal, 30 pre-perimetric and 35 perimetric glaucoma participants were enrolled. Azimuthal filtering improved the repeatability of the normalized NFL reflectance, as measured by the pooled superpixel standard deviation (SD), from 0.73 to 0.57 dB (p<0.001, paired t-test) and reduced the population SD from 2.14 to 1.78 dB (p<0.001, t-test). Most glaucomatous reflectance maps showed characteristic patterns of contiguous wedge or diffuse defects. Focal NFL reflectance loss had significantly higher diagnostic sensitivity than the best NFL thickness parameter (overall, inferior, or focal loss volume): 53% v. 23% (p=0.027) in PPG eyes and 100% v. 80% (p=0.023) in PG eyes, with the specificity fixed at 99%. Conclusions: Azimuthal filtering reduces the variability of NFL reflectance measurements. Focal NFL reflectance loss has excellent glaucoma diagnostic accuracy compared to the standard NFL thickness parameters. The reflectance map may be useful for localizing NFL defects.
△ Less
Submitted 24 June, 2020;
originally announced June 2020.
-
Modeling cell migration regulated by cell-ECM micromechanical coupling
Authors:
Yu Zheng,
Hanqing Nan,
Qihui Fan,
Xiaochen Wang,
Liyu Liu,
Ruchuan Liu,
Fangfu Ye,
Bo Sun,
Yang Jiao
Abstract:
Cell migration in fibreous extracellular matrix (ECM) is crucial to many physiological and pathological processes such as tissue regeneration, immune response and cancer progression. During migration, individual cells can generate active pulling forces via actin filament contraction, which are transmitted to the ECM fibers through focal adhesion complexes, remodel the ECM, and eventually propagate…
▽ More
Cell migration in fibreous extracellular matrix (ECM) is crucial to many physiological and pathological processes such as tissue regeneration, immune response and cancer progression. During migration, individual cells can generate active pulling forces via actin filament contraction, which are transmitted to the ECM fibers through focal adhesion complexes, remodel the ECM, and eventually propagate to and can be sensed by other cells in the system. The microstructure and physical properties of the ECM can also significantly influence cell migration, e.g., via durotaxis and contact guidance. Here, we develop a computational model for cell migration regulated by cell-ECM micro-mechanical coupling. Our model explicitly takes into account a variety of cellular level processes including focal adhesion formation and disassembly, active traction force generation and cell locomotion due to actin filament contraction, transmission and propagation of tensile forces in the ECM, as well as the resulting ECM remodeling. We validate our model by accurately reproducing single-cell dynamics of MCF-10A breast cancer cells migrating on collagen gels and show that the durotaxis and contact guidance effects naturally arise as a consequence of the cell-ECM micro-mechanical interactions considered in the model. Moreover, our model predicts strongly correlated multi-cellular migration dynamics, which are resulted from the ECM-mediated mechanical coupling among the migrating cell and are subsequently verified in {\it in vitro} experiments using MCF-10A cells. Our computational model provides a robust tool to investigate emergent collective dynamics of multi-cellular systems in complex {\it in vivo} micro-environment and can be utilized to design {\it in vitro} micro-environments to guide collective behaviors and self-organization of cells.
△ Less
Submitted 16 May, 2019;
originally announced May 2019.
-
Collective dynamics and self-organization of active particles on random networks
Authors:
Hanqing Nan,
Yu Zheng,
Shaohua Chen,
Christopher Eddy,
Jianxiang Tian,
Wenxiang Xu,
Bo Sun,
Yang Jiao
Abstract:
Collective cell migration in 3D extracellular matrix (ECM) is crucial to many physiological and pathological processes. Migrating cells can generate active pulling forces via actin filament contraction, which are transmitted to the ECM fibers and lead to a dynamically evolving force network in the system. Here, we elucidate the role of such force network in regulating collective cell behaviors usi…
▽ More
Collective cell migration in 3D extracellular matrix (ECM) is crucial to many physiological and pathological processes. Migrating cells can generate active pulling forces via actin filament contraction, which are transmitted to the ECM fibers and lead to a dynamically evolving force network in the system. Here, we elucidate the role of such force network in regulating collective cell behaviors using a minimal active-particle-on-network (APN) model, in which active particles can pull the fibers and hop between neighboring nodes of the network following local durotaxis. Our model reveals a dynamic transition as the particle number density approaches a critical value, from an "absorbing" state containing isolated stationary small particle clusters, to a "active" state containing a single large cluster undergone constant dynamic reorganization. This reorganization is dominated by a subset of highly dynamic "radical" particles in the cluster, whose number also exhibits a transition at the same critical density. The transition is underlaid by the percolation of "influence spheres" due to the particle pulling forces. Our results suggest a robust mechanism based on ECM-mediated mechanical coupling for collective cell behaviors in 3D ECM.
△ Less
Submitted 25 April, 2019;
originally announced April 2019.
-
Modelling Hierarchical Flocking
Authors:
Yongnan Jia,
Tamas Vicsek
Abstract:
We present a general framework for modeling a wide selection of flocking scenarios under free boundary conditions. Several variants have been considered - including examples for the widely observed behavior of hierarchically interacting units. The models we have simulated correspond to classes of various realistic situations. Our primary goal was to investigate the stability of a flock in the pres…
▽ More
We present a general framework for modeling a wide selection of flocking scenarios under free boundary conditions. Several variants have been considered - including examples for the widely observed behavior of hierarchically interacting units. The models we have simulated correspond to classes of various realistic situations. Our primary goal was to investigate the stability of a flock in the presence of noise. Some of our findings are counter-intuitive in the first approximation, e.g., if the hierarchy is based purely on dominance (an uneven contribution of the neighbors to the decision about the direction of flight of a given individual) the flock is more prone to loose coherence due to perturbations even when a comparison with the standard egalitarian flock is made. Thus, we concentrated on building models based on leaderfollower relationships. And, indeed, our findings support the concept that hierarchical organization can be very efficient in important practical cases, especially if the leaderfollower interactions (corresponding to an underlying directed network of interactions) have several levels. Efficiency here is associated with remaining stable (coherent and cohesive) even in cases when collective motion is destroyed by random perturbations. The framework we present allows a the study of several further complex interactions among the members of flocking agents.
△ Less
Submitted 21 April, 2019;
originally announced April 2019.
-
Collective ECM remodeling organizes 3D collective cancer invasion
Authors:
Jihan Kim,
Yu Zheng,
Amani A. Alobaidi,
Hanqing Nan,
Jianxiang Tian,
Yang Jiao,
Bo Sun
Abstract:
Tumor metastasis, traditionally considered as random spreading of individual cancer cells, has been shown to also involve coordinated collective invasion in the tissue space. Here we demonstrate experimentally and computationally that physical interactions between cells and extracellular matrix (ECM) support coordinated dissemination of cells from tumor organoids. We find that collective remodelin…
▽ More
Tumor metastasis, traditionally considered as random spreading of individual cancer cells, has been shown to also involve coordinated collective invasion in the tissue space. Here we demonstrate experimentally and computationally that physical interactions between cells and extracellular matrix (ECM) support coordinated dissemination of cells from tumor organoids. We find that collective remodeling of the ECM fibrous structure by cell-generated forces produces reciprocal cues to bias cell motility. As a result, dissemination of cells from tumor organoids are controlled by the organoid geometry. Our results indicate that mechanics and geometry are closely coupled in multicellular processes during metastasis and morphogenesis. This study also suggests that migrating cells in 3D ECM represent a distinct class of active particle system where the collective dynamics is governed by the micromechanical remodeling of the environment rather than direct particle-particle interactions.
△ Less
Submitted 3 April, 2019; v1 submitted 8 March, 2019;
originally announced March 2019.
-
Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy
Authors:
Hang Xie,
Yang Jiao,
Qihui Fan,
Miaomiao Hai,
Jiaen Yang,
Zhijian Hu,
Yue Yang,
Jianwei Shuai,
Guo Chen,
Ruchuan Liu,
Liyu Liu
Abstract:
A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate…
▽ More
A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations reproduce the well-established observation that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.
△ Less
Submitted 7 March, 2018;
originally announced March 2018.
-
Heterogeneous Force Chains in Cellularized Biopolymer Network
Authors:
Long Liang,
Christopher Jones,
Bo Sun,
Yang Jiao
Abstract:
Biopolymer Networks play an important role in coordinating and regulating collective cellular dynamics via a number of signaling pathways. Here, we investigate the mechanical response of a model biopolymer network due to the active contraction of embedded cells. Specifically, a graph (bond-node) model derived from confocal microscopy data is used to represent the network microstructure, and cell c…
▽ More
Biopolymer Networks play an important role in coordinating and regulating collective cellular dynamics via a number of signaling pathways. Here, we investigate the mechanical response of a model biopolymer network due to the active contraction of embedded cells. Specifically, a graph (bond-node) model derived from confocal microscopy data is used to represent the network microstructure, and cell contraction is modeled by applying correlated displacements at specific nodes, representing the focal adhesion sites. A force-based stochastic relaxation method is employed to obtain force-balanced network under cell contraction. We find that the majority of the forces are carried by a small number of heterogeneous force chains emitted from the contracting cells. The force chains consist of fiber segments that either possess a high degree of alignment before cell contraction or are aligned due to the reorientation induced by cell contraction. Large fluctuations of the forces along different force chains are observed. Importantly, the decay of the forces along the force chains is significantly slower than the decay of radially averaged forces in the system. These results suggest that the fibrous nature of biopolymer network structure can support long-range force transmission and thus, long-range mechanical signaling between cells.
△ Less
Submitted 18 August, 2015; v1 submitted 9 August, 2015;
originally announced August 2015.
-
A Cellular Automaton Model for Tumor Dormancy: Emergence of a Proliferative Switch
Authors:
Duyu Chen,
Yang Jiao,
Salvatore Torquato
Abstract:
Malignant cancers that lead to fatal outcomes for patients may remain dormant for very long periods of time. Although individual mechanisms such as cellular dormancy, angiogenic dormancy and immunosurveillance have been proposed, a comprehensive understanding of cancer dormancy and the "switch" from a dormant to a proliferative state still needs to be strengthened from both a basic and clinical po…
▽ More
Malignant cancers that lead to fatal outcomes for patients may remain dormant for very long periods of time. Although individual mechanisms such as cellular dormancy, angiogenic dormancy and immunosurveillance have been proposed, a comprehensive understanding of cancer dormancy and the "switch" from a dormant to a proliferative state still needs to be strengthened from both a basic and clinical point of view. Computational modeling enables one to explore a variety of scenarios for possible but realistic microscopic dormancy mechanisms and their predicted outcomes. The aim of this paper is to devise such a predictive computational model of dormancy with an emergent "switch" behavior. Specifically, we generalize a previous cellular automaton (CA) model for proliferative growth of solid tumor that now incorporates a variety of cell-level tumor-host interactions and different mechanisms for tumor dormancy, for example the effects of the immune system. Our new CA rules induce a natural "competition" between the tumor and tumor suppression factors in the microenvironment. This competition either results in a "stalemate" for a period of time in which the tumor either eventually wins (spontaneously emerges) or is eradicated; or it leads to a situation in which the tumor is eradicated before such a "stalemate" could ever develop. We also predict that if the number of actively dividing cells within the proliferative rim of the tumor reaches a critical, yet low level, the dormant tumor has a high probability to resume rapid growth. Our findings may shed light on the fundamental understanding of cancer dormancy.
△ Less
Submitted 17 October, 2014;
originally announced October 2014.