Source author record

Padmini Rangamani

Padmini Rangamani appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2026arXiv

A finite element framework for solving coupled multiphysics problem with moving boundaries in cell biophysics

Cellular morphodynamics requires solving systems of coupled partial differential equations on moving bulk and surface domains, where advection-dominant transport, structure preservation, and severe mesh distortions make robust simulation difficult. We present a holistic finite element framework that jointly addresses these obstacles for biophysical applications by combining model-agnostic structure-preserving postprocessing, ALE-based mesh redistribution strategies driven by surface-tangential velocities, and stabilized discretization for advection-diffusion-reaction problems tailored to evolving domains. The methodology is modular and applies to advection-diffusion-reaction systems, Cahn-Hilliard phase separation, Helfrich-type geometric flows, as well as their staggered and potentially mixed-dimensional couplings. We provide a concise notation for evolving bulk and surface geometries, extend positivity-, bound-, and mass-preserving projections to moving meshes, and develop a two-step redistribution procedure that maintains element quality without remeshing. Convergence studies, manufactured solutions, and biologically motivated test cases -- including tumor-growth surrogates and phase segregation on deformable membranes -- demonstrate accuracy, stability, and versatility across the problem classes considered.

preprint2023arXiv

Bayesian Parameter Estimation for Dynamical Models in Systems Biology

Dynamical systems modeling, particularly via systems of ordinary differential equations, has been used to effectively capture the temporal behavior of different biochemical components in signal transduction networks. Despite the recent advances in experimental measurements, including sensor development and '-omics' studies that have helped populate protein-protein interaction networks in great detail, modeling in systems biology lacks systematic methods to estimate kinetic parameters and quantify associated uncertainties. This is because of multiple reasons, including sparse and noisy experimental measurements, lack of detailed molecular mechanisms underlying the reactions, and missing biochemical interactions. Additionally, the inherent nonlinearities with respect to the states and parameters associated with the system of differential equations further compound the challenges of parameter estimation. In this study, we propose a comprehensive framework for Bayesian parameter estimation and complete quantification of the effects of uncertainties in the data and models. We apply these methods to a series of signaling models of increasing mathematical complexity. Systematic analysis of these dynamical systems showed that parameter estimation depends on data sparsity, noise level, and model structure, including the existence of multiple steady states. These results highlight how focused uncertainty quantification can enrich systems biology modeling and enable additional quantitative analyses for parameter estimation.

preprint2022arXiv

Can biophysical models give insight into the synaptic changes associated with addiction?

Effective treatments that prevent or reduce drug relapse vulnerability should be developed to relieve the high burden of drug addiction to society. This will only be possible by enhancing the understanding of the molecular mechanisms underlying the neurobiology of addiction. Recent experimental data have shown that dendritic spines, small protrusions from the dendrites that receive input from excitatory neurons, from spiny neurons in the nucleus accumbens exhibit morphological changes during drug exposure and withdrawal. Moreover, these changes relate to the characteristic drug-seeking behavior of addiction. However, due to the complexity of the dendritic spines, we do not yet fully understand the processes underlying their structural changes in response to different inputs. We propose that biophysical models can enhance the current understanding of these processes by incorporating different, and sometimes, discrepant experimental data to identify the shared underlying mechanisms and generate experimentally testable hypotheses. This review aims to give an up-to-date report on biophysical models of dendritic spines, focusing on those models that describe their shape changes, which are well-known to relate to learning and memory. Moreover, it examines how these models can enhance our understanding of the effect of the drugs and the synaptic changes during disease progression.

preprint2022arXiv

Formation of protein-mediated tubes is governed by a snapthrough transition

Plasma membrane tubes are ubiquitous in cellular membranes and in the membranes of intracellular organelles. They play crucial roles in trafficking, ion transport, and cellular motility. The formation of plasma membrane tubes can be due to localized forces acting on the membrane or by curvature-induced by membrane-bound proteins. Here, we present a mathematical framework to model cylindrical tubular protrusions formed by proteins that induce anisotropic spontaneous curvature. Our analysis revealed that the tube radius depends on an effective tension that includes contributions from the bare membrane tension and the protein-induced curvature. We also found that the length of the tube undergoes an abrupt transition from a short, dome-shaped membrane to a long cylinder and this transition is characteristic of a snapthrough instability. Finally, we show that the snapthrough instability depends on the different parameters including coat area, bending modulus, and extent of protein-induced curvature. Our findings have implications for tube formation due to BAR-domain proteins in processes such as endocytosis, t-tubule formation in myocytes, and cristae formation in mitochondria.

preprint2021arXiv

Mem3DG: Modeling Membrane Mechanochemical Dynamics in 3D using Discrete Differential Geometry

Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared to approximation-based methods (e.g., finite element method), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, $\texttt{Mem3DG}$, for modeling 3D mechanochemical dynamics of membranes based on Discrete Differential Geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, $\texttt{Mem3DG}$ is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.

preprint2020arXiv

Stability analysis of a signaling circuit with dual species of GTPase switches

GTPases are molecular switches that regulate a wide range of cellular processes, such as organelle biogenesis, position, shape, and signal transduction. These enzymes operate by toggling between an active ("ON") guanosine triphosphate (GTP)-bound state and an inactive ("OFF") guanosine diphosphate (GDP)-bound state; such a toggle is regulated by GEFs (guanine nucleotide exchange factors) and GAPs (GTPase activating proteins). Here we dissect a network motif between monomeric (m) and trimeric (t) GTPases assembled exclusively in eukaryotic cells of multicellular organisms. To this end, we develop a system of ordinary differential equations in which these two classes of GTPases are interlinked conditional to their ON/OFF states within a motif through feedforward and feedback loops. We provide formulas for the steady states of the system and perform local stability analysis to investigate the role of the different connections between the GTPase switches. A feedforward from the active mGTPase to the GEF of the tGTPase was sufficient to provide two locally stable states: one where both active/inactive forms of the mGTPase can be interpreted as having low concentrations and the other where both m- and tGTPase have high concentrations. When a feedback loop from the GEF of the tGTPase to the GAP of the mGTPase was added to the feedforward system, two other locally stable states emerged, both having the tGTPase inactivated and being interpreted as having low active tGTPase concentrations. Finally, the addition of a second feedback loop, from the active tGTPase to the GAP of the mGTPase, gives rise to a family of steady states parametrized by the inactive tGTPase concentrations. Our findings reveal that the coupling of these two different GTPase motifs can dramatically change their steady state behaviors and shed light on how such coupling may impact information processing in eukaryotic cells.

preprint2019arXiv

3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.

preprint2019arXiv

An Open Source Mesh Generation Platform for Biophysical Modeling Using Realistic Cellular Geometries

Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this study, we outline the steps for a naïve user to approach GAMer 2, a mesh generation code written in C++ designed to convert structural datasets to realistic geometric meshes, while preserving the underlying shapes. We present two example cases, 1) mesh generation at the subcellular scale as informed by electron tomography, and 2) meshing a protein with structure from x-ray crystallography. We further demonstrate that the meshes generated by GAMer are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.

preprint2016arXiv

An introduction to linear stability analysis for deciphering spatial patterns in signaling networks

Mathematical modeling is now used commonly in the analysis of signaling networks. With advances in high resolution microscopy, the spatial location of different signaling molecules and the spatio-temporal dynamics of signaling microdomains are now widely acknowledged as key features of biochemical signal transduction. Reaction-diffusion mechanisms are commonly used to model such features, often with a heavy reliance on numerical simulations to obtain results. However, simulations are parameter dependent and may not be able to provide an understanding of the full range of the system responses. Analytical approaches on the other hand provide a framework to study the entire phase space. In this tutorial, we provide a largely analytical method for studying reaction-diffusion models and analyzing their stability properties. Using two representative biological examples, we demonstrate how this approach can guide experimental design.

preprint2016arXiv

Membrane tension is a key determinant of bud morphology in clathrin-mediated endocytosis

In clathrin-mediated endocytosis (CME), clathrin and various adaptor proteins coat a patch of the plasma membrane, which is reshaped to form a budded vesicle. Experimental studies have demonstrated that elevated membrane tension can inhibit bud formation by a clathrin coat. In this study, we investigate the impact of membrane tension on the mechanics of membrane budding by simulating clathrin coats that either grow in area or progressively induce greater curvature. At low membrane tension, progressively increasing the area of a curvature-generating coat causes the membrane to smoothly evolve from a flat to budded morphology, whereas the membrane remains essentially flat at high membrane tensions. Interestingly, at physiologically relevant, intermediate membrane tensions, the shape evolution of the membrane undergoes a snapthrough instability in which increasing coat area causes the membrane to "snap" from an open, U-shaped bud to a closed, $Ω$-shaped bud. This instability is accompanied by a large energy barrier, which could cause a developing endocytic pit to stall if the binding energy of additional coat is insufficient to overcome this barrier. Similar results were found for a coat of constant area in which the spontaneous curvature progressively increases. Additionally, a pulling force on the bud, simulating a force from actin polymerization, is sufficient to drive a transition from an open to closed bud, overcoming the energy barrier opposing this transition.

preprint2016arXiv

Shear-Induced Nitric Oxide Production by Endothelial Cells

We present a biochemical model of the wall shear stress (WSS)-induced activation of endothelial nitric oxide synthase (eNOS) in an endothelial cell (EC). The model includes three key mechanotransducers: mechanosensing ion channels, integrins and G-protein-coupled receptors. The reaction cascade consists of two interconnected parts. The first is rapid activation of calcium, which results in formation of calcium-calmodulin complexes, followed by recruitment of eNOS from caveolae. The second is phosphoryaltion of eNOS by protein kinases PKC and AKT. The model also includes a negative feedback loop due to inhibition of calcium influx into the cell by cyclic guanosine monophosphate (cGMP). In this feedback, increased nitric oxide (NO) levels cause an increase in cGMP levels, so that cGMP inhibition of calcium influx can limit NO production. The model was used to predict the dynamics of NO production by an EC subjected to a step increase of WSS from zero to a finite physiologically relevant value. Among several experimentally observed features, the model predicts a highly nonlinear, biphasic transient behavior of eNOS activation and NO production: a rapid initial activation due to the very rapid influx of calcium into the cytosol (occurring within 1 to 5 minutes) is followed by a sustained period of activation due to protein kinases.

preprint2014arXiv

ATP concentration regulates enzyme kinetics

Adenosine 5'-triphosphate (ATP) is the nearly ubiquitous "energy currency" of living organisms, and thus is a crucial participant in the majority of enzymatic reactions. The standard models in enzyme kinetics generally ignore the temporal dynamics of ATP because it is assumed to be present in large excess. However, this assumption may not hold in many situations of cellular stress where ATP concentrations may be comparable to substrate levels. Here, we demonstrate the importance of ATP concentration on the dynamics of multi-enzyme reactions by explicit consideration of ATP as a secondary substrate for an enzyme. We apply our model to the mitogen-activated protein (MAP) kinase cascade, which is involved in the regulation of a vast range of cellular activities. We show that three fundamental features of this signaling network --- (i) duration of response, (ii) signal amplification, and (iii) ultrasensitivity to stimulus concentration --- are all dependent on ATP concentration. Our results indicate that the concentration of ATP regulates the response of the MAP kinase activation network, and potentially suggests another possible mechanism for disruption of the cascade in pathogenic states.

preprint2014arXiv

Protein-induced membrane curvature changes membrane tension

Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this paper, we show that the classical elastic model of lipid membranes cannot account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, and a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the new viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvatures. This change in membrane tension can have a functional impact in many biological phenomena where proteins interact with membranes.