Near-infrared voltage-sensitive dyes based on chromene donor.

Voltage-sensitive dyes (VSDs) are used to image electrical activity in cells and tissues with submillisecond time resolution. Most of these fast sensors are constructed from push-pull chromophores whose fluorescence spectra are modulated by the electric field across the cell membrane. It was found that the substitution of naphthalene with chromene produces a 60 to 80 nm red-shift in absorption and emission spectra while maintaining fluorescence quantum efficiency and voltage sensitivity. One dye was applied to ex vivo murine heart with excitation at 730 nm, by far the longest wavelength reported in voltage imaging. This VSD resolves cardiac action potentials in single trials with 12% ΔF/F per action potential. The well-separated excitation spectra between these long-wavelength VSDs and channelrhodopsin (ChR2) enabled monitoring of action potential propagation in ChR2 hearts without any perturbation of electrical dynamics. Importantly, by employing spatially localized optogenetic manipulation, action potential dynamics can be assessed in an all-optical fashion with no artifact related to optical cross-talk between the reporter and actuator. These new environmentally sensitive chromene-based chromophores are also likely to have applications outside voltage imaging.

Optogenetic confirmation of transverse-tubular membrane excitability in intact cardiac myocytes.

T-tubules (TT) form a complex network of sarcolemmal membrane invaginations, essential for well-co-ordinated excitation-contraction coupling (ECC) and thus homogeneous mechanical activation of cardiomyocytes. ECC is initiated by rapid depolarization of the sarcolemmal membrane. Whether TT membrane depolarization is active (local generation of action potentials; AP) or passive (following depolarization of the outer cell surface sarcolemma; SS) has not been experimentally validated in cardiomyocytes. Based on the assessment of ion flux pathways needed for AP generation, we hypothesize that TT are excitable. We therefore explored TT excitability experimentally, using an all-optical approach to stimulate and record trans-membrane potential changes in TT that were structurally disconnected, and hence electrically insulated, from the SS membrane by transient osmotic shock. Our results establish that cardiomyocyte TT can generate AP. These AP show electrical features that differ substantially from those observed in SS, consistent with differences in the density of ion channels and transporters in the two different membrane domains. We propose that TT-generated AP represent a safety mechanism for TT AP propagation and ECC, which may be particularly relevant in pathophysiological settings where morpho-functional changes reduce the electrical connectivity between SS and TT membranes. KEY POINTS: Cardiomyocytes are characterized by a complex network of membrane invaginations (the T-tubular system) that propagate action potentials to the core of the cell, causing uniform excitation-contraction coupling across the cell. In the present study, we investigated whether the T-tubular system is able to generate action potentials autonomously, rather than following depolarization of the outer cell surface sarcolemma. For this purpose, we developed a fully optical platform to probe and manipulate the electrical dynamics of subcellular membrane domains. Our findings demonstrate that T-tubules are intrinsically excitable, revealing distinct characteristics of self-generated T-tubular action potentials. This active electrical capability would protect cells from voltage drops potentially occurring within the T-tubular network.

MolClustPy: a Python package to characterize multivalent biomolecular clusters.

SUMMARY: Low-affinity interactions among multivalent biomolecules may lead to the formation of molecular complexes that undergo phase transitions to become supply-limited large clusters. In stochastic simulations, such clusters display a wide range of sizes and compositions. We have developed a Python package, MolClustPy, which performs multiple stochastic simulation runs using NFsim (Network-Free stochastic simulator); MolClustPy characterizes and visualizes the distribution of cluster sizes, molecular composition, and bonds across molecular clusters. The statistical analysis offered by MolClustPy is readily applicable to other stochastic simulation software, such as SpringSaLaD and ReaDDy. AVAILABILITY AND IMPLEMENTATION: The software is implemented in Python. A detailed Jupyter notebook is provided to enable convenient running. Code, user guide, and examples are freely available at https://molclustpy.github.io/.

Beyond analytic solution: Analysis of FRAP experiments by spatial simulation of the forward problem.

Fluorescence redistribution after photobleaching is a commonly used method to understand the dynamic behavior of molecules within cells. Analytic solutions have been developed for specific, well-defined models of dynamic behavior in idealized geometries, but these solutions are inaccurate in complex geometries or when complex binding and diffusion behaviors exist. We demonstrate the use of numerical reaction-diffusion simulations using the Virtual Cell software platform to model fluorescence redistribution after photobleaching experiments. Multiple simulations employing parameter scans and varying bleaching locations and sizes can help to bracket diffusion coefficients and kinetic rate constants in complex image-based geometries. This approach is applied to problems in membrane surface diffusion as well as diffusion and binding in cytosolic volumes in complex cell geometries. In addition, we model diffusion and binding within phase-separated biomolecular condensates (liquid droplets). These are modeled as spherical low-affinity binding domains that also define a high viscosity medium for exchange of the free fluorescently labeled ligand with the external cytosol.

The maximum solubility product marks the threshold for condensation of multivalent biomolecules.

Clustering of weakly interacting multivalent biomolecules underlies the formation of membraneless compartments known as condensates. As opposed to single-component (homotypic) systems, the concentration dependence of multicomponent (heterotypic) condensate formation is not well understood. We previously proposed the solubility product (SP), the product of monomer concentrations in the dilute phase, as a tool for understanding the concentration dependence of multicomponent systems. In this study, we further explore the limits of the SP concept using spatial Langevin dynamics and rule-based stochastic simulations. We show, for a variety of idealized molecular structures, how the maximum SP coincides with the onset of the phase transition, i.e., the formation of large clusters. We reveal the importance of intracluster binding in steering the free and cluster phase molecular distributions. We also show how structural features of biomolecules shape the SP profiles. The interplay of flexibility, length, and steric hindrance of linker regions controls the phase transition threshold. Remarkably, when SPs are normalized to nondimensional variables and plotted against the concentration scaled to the threshold for phase transition, the curves all coincide independent of the structural features of the binding partners. Similar coincidence is observed for the normalized clustering versus concentration plots. Overall, the principles derived from these systematic models will help guide and interpret in vitro and in vivo experiments on the biophysics of biomolecular condensates.

Optical mapping of cardiac electromechanics in beating in vivo hearts.

Optical mapping has been widely used in the study of cardiac electrophysiology in motion-arrested, ex vivo heart preparations. Recent developments in motion artifact mitigation techniques have made it possible to optically map beating ex vivo hearts, enabling the study of cardiac electromechanics using optical mapping. However, the ex vivo setting imposes limitations on optical mapping such as altered metabolic states, oversimplified mechanical loads, and the absence of neurohormonal regulation. In this study, we demonstrate optical electromechanical mapping in an in vivo heart preparation. Swine hearts were exposed via median sternotomy. Voltage-sensitive dye, either di-4-ANEQ(F)PTEA or di-5-ANEQ(F)PTEA, was injected into the left anterior descending artery. Fluorescence was excited by alternating green and amber light for excitation ratiometry. Cardiac motion during sinus and paced rhythm was tracked using a marker-based method. Motion tracking and excitation ratiometry successfully corrected most motion artifact in the membrane potential signal. Marker-based motion tracking also allowed simultaneous measurement of epicardial deformation. Reconstructed membrane potential and mechanical deformation measurements were validated using monophasic action potentials and sonomicrometry, respectively. Di-5-ANEQ(F)PTEA produced longer working time and higher signal/noise ratio than di-4-ANEQ(F)PTEA. In addition, we demonstrate potential applications of the new optical mapping system including electromechanical mapping during vagal nerve stimulation, fibrillation/defibrillation. and acute regional ischemia. In conclusion, although some technical limitations remain, optical mapping experiments that simultaneously image electrical and mechanical function can be conducted in beating, in vivo hearts.

Recent advances and current limitations of available technology to optically manipulate and observe cardiac electrophysiology.

Optogenetics, utilising light-reactive proteins to manipulate tissue activity, are a relatively novel approach in the field of cardiac electrophysiology. We here provide an overview of light-activated transmembrane channels (optogenetic actuators) currently applied in strategies to modulate cardiac activity, as well as newly developed variants yet to be implemented in the heart. In addition, we touch upon genetically encoded indicators (optogenetic sensors) and fluorescent dyes to monitor tissue activity, including cardiac transmembrane potential and ion homeostasis. The combination of the two allows for all-optical approaches to monitor and manipulate the heart without any physical contact. However, spectral congestion poses a major obstacle, arising due to the overlap of excitation/activation and emission spectra of various optogenetic proteins and/or fluorescent dyes, resulting in optical crosstalk. Therefore, optogenetic proteins and fluorescent dyes should be carefully selected to avoid optical crosstalk and consequent disruptions in readouts and/or cellular activity. We here present a novel approach to simultaneously monitor transmembrane potential and cytosolic calcium, while also performing optogenetic manipulation. For this, we used the novel voltage-sensitive dye ElectroFluor 730p and the cytosolic calcium indicator X-Rhod-1 in mouse hearts expressing channelrhodopsin-2 (ChR2). By exploiting the isosbestic point of ElectroFluor 730p and avoiding the ChR2 activation spectrum, we here introduce a novel optical imaging and manipulation approach with minimal crosstalk. Future developments in both optogenetic proteins and fluorescent dyes will allow for additional and more optimised strategies, promising a bright future for all-optical approaches in the field of cardiac electrophysiology.

Mechanism of actin filament nucleation.

We used computational methods to analyze the mechanism of actin filament nucleation. We assumed a pathway where monomers form dimers, trimers, and tetramers that then elongate to form filaments but also considered other pathways. We aimed to identify the rate constants for these reactions that best fit experimental measurements of polymerization time courses. The analysis showed that the formation of dimers and trimers is unfavorable because the association reactions are orders of magnitude slower than estimated in previous work rather than because of rapid dissociation of dimers and trimers. The 95% confidence intervals calculated for the four rate constants spanned no more than one order of magnitude. Slow nucleation reactions are consistent with published high-resolution structures of actin filaments and molecular dynamics simulations of filament ends. One explanation for slow dimer formation, which we support with computational analysis, is that actin monomers are in a conformational equilibrium with a dominant conformation that cannot participate in the nucleation steps.

The Interplay of Structural and Cellular Biophysics Controls Clustering of Multivalent Molecules.

Dynamic molecular clusters are assembled through weak multivalent interactions and are platforms for cellular functions, especially receptor-mediated signaling. Clustering is also a prerequisite for liquid-liquid phase separation. It is not well understood, however, how molecular structure and cellular organization control clustering. Using coarse-grained kinetic Langevin dynamics, we performed computational experiments on a prototypical ternary system modeled after membrane-bound nephrin, the adaptor Nck1, and the actin nucleation promoting factor NWASP. Steady-state cluster size distributions favored stoichiometries that optimized binding (stoichiometry matching) but still were quite broad. At high concentrations, the system can be driven beyond the saturation boundary such that cluster size is limited only by the number of available molecules. This behavior would be predictive of phase separation. Domains close to binding sites sterically inhibited clustering much less than terminal domains because the latter effectively restrict access to the cluster interior. Increased flexibility of interacting molecules diminished clustering by shielding binding sites within compact conformations. Membrane association of nephrin increased the cluster size distribution in a density-dependent manner. These properties provide insights into how molecular ensembles function to localize and amplify cell signaling.

mol2sphere: spherical decomposition of multi-domain molecules for visualization and coarse grained spatial modeling.

Motivation: Proteins, especially those involved in signaling pathways are composed of functional modules connected by linker domains with varying degrees of flexibility. To understand the structure-function relationships in these macromolecules, it is helpful to visualize the geometric arrangement of domains. Furthermore, accurate spatial representation of domain structure is necessary for coarse-grain models of the multi-molecular interactions that comprise signaling pathways. Results: We introduce a new tool, mol2sphere, that transforms the atomistic structure of a macromolecule into a series of linked spheres corresponding to domains. It does this with a k-means clustering algorithm. It may be used for visualization or for coarse grain modeling and simulation. Availability and implementation: PyMOL plugin, source, and documentation.https://nmrbox.org/registry/mol2sphere. SpringSaLaD executables and documentation: http://vcell.org/ssalad, SpringSaLaD v.2 source: https://github.com/jmasison/SpringSaLaD.

Fragility of foot process morphology in kidney podocytes arises from chaotic spatial propagation of cytoskeletal instability.

Kidney podocytes' function depends on fingerlike projections (foot processes) that interdigitate with those from neighboring cells to form the glomerular filtration barrier. The integrity of the barrier depends on spatial control of dynamics of actin cytoskeleton in the foot processes. We determined how imbalances in regulation of actin cytoskeletal dynamics could result in pathological morphology. We obtained 3-D electron microscopy images of podocytes and used quantitative features to build dynamical models to investigate how regulation of actin dynamics within foot processes controls local morphology. We find that imbalances in regulation of actin bundling lead to chaotic spatial patterns that could impair the foot process morphology. Simulation results are consistent with experimental observations for cytoskeletal reconfiguration through dysregulated RhoA or Rac1, and they predict compensatory mechanisms for biochemical stability. We conclude that podocyte morphology, optimized for filtration, is intrinsically fragile, whereby local transient biochemical imbalances may lead to permanent morphological changes associated with pathophysiology.

Perspectives on Sharing Models and Related Resources in Computational Biomechanics Research.

The role of computational modeling for biomechanics research and related clinical care will be increasingly prominent. The biomechanics community has been developing computational models routinely for exploration of the mechanics and mechanobiology of diverse biological structures. As a result, a large array of models, data, and discipline-specific simulation software has emerged to support endeavors in computational biomechanics. Sharing computational models and related data and simulation software has first become a utilitarian interest, and now, it is a necessity. Exchange of models, in support of knowledge exchange provided by scholarly publishing, has important implications. Specifically, model sharing can facilitate assessment of reproducibility in computational biomechanics and can provide an opportunity for repurposing and reuse, and a venue for medical training. The community's desire to investigate biological and biomechanical phenomena crossing multiple systems, scales, and physical domains, also motivates sharing of modeling resources as blending of models developed by domain experts will be a required step for comprehensive simulation studies as well as the enhancement of their rigor and reproducibility. The goal of this paper is to understand current perspectives in the biomechanics community for the sharing of computational models and related resources. Opinions on opportunities, challenges, and pathways to model sharing, particularly as part of the scholarly publishing workflow, were sought. A group of journal editors and a handful of investigators active in computational biomechanics were approached to collect short opinion pieces as a part of a larger effort of the IEEE EMBS Computational Biology and the Physiome Technical Committee to address model reproducibility through publications. A synthesis of these opinion pieces indicates that the community recognizes the necessity and usefulness of model sharing. There is a strong will to facilitate model sharing, and there are corresponding initiatives by the scientific journals. Outside the publishing enterprise, infrastructure to facilitate model sharing in biomechanics exists, and simulation software developers are interested in accommodating the community's needs for sharing of modeling resources. Encouragement for the use of standardized markups, concerns related to quality assurance, acknowledgement of increased burden, and importance of stewardship of resources are noted. In the short-term, it is advisable that the community builds upon recent strategies and experiments with new pathways for continued demonstration of model sharing, its promotion, and its utility. Nonetheless, the need for a long-term strategy to unify approaches in sharing computational models and related resources is acknowledged. Development of a sustainable platform supported by a culture of open model sharing will likely evolve through continued and inclusive discussions bringing all stakeholders at the table, e.g., by possibly establishing a consortium.

Compartmental and Spatial Rule-Based Modeling with Virtual Cell.

In rule-based modeling, molecular interactions are systematically specified in the form of reaction rules that serve as generators of reactions. This provides a way to account for all the potential molecular complexes and interactions among multivalent or multistate molecules. Recently, we introduced rule-based modeling into the Virtual Cell (VCell) modeling framework, permitting graphical specification of rules and merger of networks generated automatically (using the BioNetGen modeling engine) with hand-specified reaction networks. VCell provides a number of ordinary differential equation and stochastic numerical solvers for single-compartment simulations of the kinetic systems derived from these networks, and agent-based network-free simulation of the rules. In this work, compartmental and spatial modeling of rule-based models has been implemented within VCell. To enable rule-based deterministic and stochastic spatial simulations and network-free agent-based compartmental simulations, the BioNetGen and NFSim engines were each modified to support compartments. In the new rule-based formalism, every reactant and product pattern and every reaction rule are assigned locations. We also introduce the rule-based concept of molecular anchors. This assures that any species that has a molecule anchored to a predefined compartment will remain in this compartment. Importantly, in addition to formulation of compartmental models, this now permits VCell users to seamlessly connect reaction networks derived from rules to explicit geometries to automatically generate a system of reaction-diffusion equations. These may then be simulated using either the VCell partial differential equations deterministic solvers or the Smoldyn stochastic simulator.

Rule-based modeling with Virtual Cell.

UNLABELLED: Rule-based modeling is invaluable when the number of possible species and reactions in a model become too large to allow convenient manual specification. The popular rule-based software tools BioNetGen and NFSim provide powerful modeling and simulation capabilities at the cost of learning a complex scripting language which is used to specify these models. Here, we introduce a modeling tool that combines new graphical rule-based model specification with existing simulation engines in a seamless way within the familiar Virtual Cell (VCell) modeling environment. A mathematical model can be built integrating explicit reaction networks with reaction rules. In addition to offering a large choice of ODE and stochastic solvers, a model can be simulated using a network free approach through the NFSim simulation engine. AVAILABILITY AND IMPLEMENTATION: Available as VCell (versions 6.0 and later) at the Virtual Cell web site (http://vcell.org/). The application installs and runs on all major platforms and does not require registration for use on the user's computer. Tutorials are available at the Virtual Cell website and Help is provided within the software. Source code is available at Sourceforge. CONTACT: vcell_support@uchc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Defects in T-tubular electrical activity underlie local alterations of calcium release in heart failure.

Action potentials (APs), via the transverse axial tubular system (TATS), synchronously trigger uniform Ca(2+) release throughout the cardiomyocyte. In heart failure (HF), TATS structural remodeling occurs, leading to asynchronous Ca(2+) release across the myocyte and contributing to contractile dysfunction. In cardiomyocytes from failing rat hearts, we previously documented the presence of TATS elements which failed to propagate AP and displayed spontaneous electrical activity; the consequence for Ca(2+) release remained, however, unsolved. Here, we develop an imaging method to simultaneously assess TATS electrical activity and local Ca(2+) release. In HF cardiomyocytes, sites where T-tubules fail to conduct AP show a slower and reduced local Ca(2+) transient compared with regions with electrically coupled elements. It is concluded that TATS electrical remodeling is a major determinant of altered kinetics, amplitude, and homogeneity of Ca(2+) release in HF. Moreover, spontaneous depolarization events occurring in failing T-tubules can trigger local Ca(2+) release, resulting in Ca(2+) sparks. The occurrence of tubule-driven depolarizations and Ca(2+) sparks may contribute to the arrhythmic burden in heart failure.

SpringSaLaD: A Spatial, Particle-Based Biochemical Simulation Platform with Excluded Volume.

We introduce Springs, Sites, and Langevin Dynamics (SpringSaLaD), a comprehensive software platform for spatial, stochastic, particle-based modeling of biochemical systems. SpringSaLaD models biomolecules in a coarse-grained manner as a group of linked spherical sites with excluded volume. This mesoscopic approach bridges the gap between highly detailed molecular dynamics simulations and the various methods used to study network kinetics and diffusion at the cellular level. SpringSaLaD is a standalone tool that supports model building, simulation, visualization, and data analysis, all through a user-friendly graphical user interface that should make it more accessible than tools built into more comprehensive molecular dynamics infrastructures. Importantly, for bimolecular reactions we derive an exact expression relating the macroscopic on-rate to the various microscopic parameters with the inclusion of excluded volume; this makes SpringSaLaD more accurate than other tools, which rely on approximate relationships between these parameters.

T-Tubular Electrical Defects Contribute to Blunted ß-Adrenergic Response in Heart Failure.

Alterations of the ß-adrenergic signalling, structural remodelling, and electrical failure of T-tubules are hallmarks of heart failure (HF). Here, we assess the effect of ß-adrenoceptor activation on local Ca(2+) release in electrically coupled and uncoupled T-tubules in ventricular myocytes from HF rats. We employ an ultrafast random access multi-photon (RAMP) microscope to simultaneously record action potentials and Ca(2+) transients from multiple T-tubules in ventricular cardiomyocytes from a HF rat model of coronary ligation compared to sham-operated rats as a control. We confirmed that ß-adrenergic stimulation increases the frequency of Ca(2+) sparks, reduces Ca(2+) transient variability, and hastens the decay of Ca(2+) transients: all these effects are similarly exerted by ß-adrenergic stimulation in control and HF cardiomyocytes. Conversely, ß-adrenergic stimulation in HF cells accelerates a Ca(2+) rise exclusively in the proximity of T-tubules that regularly conduct the action potential. The delayed Ca(2+) rise found at T-tubules that fail to conduct the action potential is instead not affected by ß-adrenergic signalling. Taken together, these findings indicate that HF cells globally respond to ß-adrenergic stimulation, except at T-tubules that fail to conduct action potentials, where the blunted effect of the ß-adrenergic signalling may be directly caused by the lack of electrical activity.

Pathway Commons at virtual cell: use of pathway data for mathematical modeling.

UNLABELLED: Pathway Commons is a resource permitting simultaneous queries of multiple pathway databases. However, there is no standard mechanism for using these data (stored in BioPAX format) to annotate and build quantitative mathematical models. Therefore, we developed a new module within the virtual cell modeling and simulation software. It provides pathway data retrieval and visualization and enables automatic creation of executable network models directly from qualitative connections between pathway nodes. AVAILABILITY AND IMPLEMENTATION: Available at Virtual Cell (http://vcell.org/). Application runs on all major platforms and does not require registration for use on the user's computer. Tutorials and video are available at user guide page.

Proteo-lipobeads for the oriented encapsulation of membrane proteins.

As a surrogate of live cells, proteo-lipobeads are presented, encapsulating functional membrane proteins in a strict orientation into a lipid bilayer. Assays can be performed just as on live cells, for example using fluorescence measurements. As a proof of concept, we have demonstrated proton transport through cytochrome c oxidase.