There is an urgent need for the development of novel and truly rapid (equal or less than 1 hour) antibiotic susceptibility testing (AST) platforms in order to provide best antimicrobial prescribing practices and to help reduce the increasing global threat of antibiotic resistance. A 785 nm surface enhanced Raman spectroscopy (SERS) based phenotypic methodology is described that results in accurate minimum inhibitory concentration (MIC) determinations for all tested strain-antibiotic pairs. The SERS-AST procedure results in accurate MICs, the key quantitative measure of in vitro drug susceptibility, in 1 hour, including a 30-minute incubation period. The method is effective for both Gram positive and negative species, and for antibiotics with different initial primary targets of antibiotic activity, and for both bactericidal and bacteriostatic antibiotics. The molecular level mechanism of this methodology is described. Bacterial SERS spectra are due to secreted purine nucleotide degradation products (principally adenine, guanine, xanthine and hypoxanthine) resulting from water washing induced bacterial stringent response and the resulting (p)ppGpp alarmone mediates nucleobase formation from unneeded tRNA and rRNA. The rewiring of metabolic responses resulting from the secondary metabolic effects of antibiotic exposure during the 30-minute incubation period accounts for the dose dependence of the SERS spectral intensities which are used to accurately yield the MIC. This is the fastest demonstrated AST method yielding MICs.
Voltage distribution in sub-cellular micro-domains such as neuronal synapses, small protrusions or dendritic spines regulates the opening and closing of ionic channels, energy production and thus cellular homeostasis and excitability. Yet how voltage changes at such a small scale in vivo remains challenging due to the experimental diffraction limit, large signal fluctuations and the still limited resolution of fast voltage indicators. Here, we study the voltage distribution in nano-compartments using a computational approach based on the Poisson-Nernst-Planck equations for the electro-diffusion motion of ions, where inward and outward fluxes are generated between channels. We report a current-voltage (I-V) logarithmic relationship generalizing Nernst law that reveals how the local membrane curvature modulates the voltage. We further find that an influx current penetrating a cellular electrolyte can lead to perturbations from tens to hundreds of nanometers deep depending on the local channels organization. Finally, we show that the neck resistance of dendritic spines can be completely shunted by the transporters located on the head boundary, facilitating ionic flow. To conclude, we propose that voltage is regulated at a subcellular level by channels organization, membrane curvature and narrow passages.
Oliver M. Drozdowski, Falko Ziebert, Ulrich S. Schwarz
Cell crawling on flat substrates is based on intracellular flows of the actin cytoskeleton that are driven by both actin polymerization at the front and myosin contractility at the back. The new experimental tool of optogenetics makes it possible to spatially control contraction and thereby possibly also cell migration. Here we analyze this situation theoretically using a one-dimensional active gel model in which the excluded volume interactions of myosin and their aggregation into minifilaments is modeled by a supercritical van der Waals fluid. This physically simple and transparent, but nonlinear and thermodynamically rigorous model predicts bistability between sessile and motile solutions. We then show that one can switch between these two states at realistic parameter ranges via optogenetic activation or inhibition of contractility, in agreement with recent experiments. We also predict the required activation strengths and initiation times.
A highly organized and densely packed lattice of molecular machinery within the sarcomeres of muscle cells powers contraction. Although many of the proteins that drive contraction have been studied extensively, the mechanical impact of fluid shearing within the lattice of molecular machinery has received minimal attention. It was recently proposed that fluid flow augments substrate transport in the sarcomere, however, this analysis used analytical models of fluid flow in the molecular machinery that could not capture its full complexity. By building a finite element model of the sarcomere, we estimate the explicit flow field, and contrast it with analytical models. Our results demonstrate that viscous drag forces on sliding filaments are surprisingly small in contrast to the forces generated by single myosin molecular motors. This model also indicates that the energetic cost of fluid flow through viscous shearing with lattice proteins is likely minimal. The model also highlights a steep velocity gradient between sliding filaments and demonstrates that the maximal radial fluid velocity occurs near the tips of the filaments. To our knowledge, this is the first computational analysis of fluid flow within the highly structured sarcomere.
Nickolay Korabel, Daniel Han, Alessandro Taloni
et al.
Trajectories of endosomes inside living eukaryotic cells are highly heterogeneous in space and time and diffuse anomalously due to a combination of viscoelasticity, caging, aggregation and active transport. Some of the trajectories display switching between persistent and anti-persistent motion while others jiggle around in one position for the whole measurement time. By splitting the ensemble of endosome trajectories into slow moving sub-diffusive and fast moving super-diffusive endosomes, we analyzed them separately. The mean squared displacements and velocity auto-correlation functions confirm the effectiveness of the splitting methods. Applying the local analysis, we show that both ensembles are characterized by a spectrum of local anomalous exponents and local generalized diffusion coefficients. Slow and fast endsomes have exponential distributions of local anomalous exponents and power law distributions of generalized diffusion coefficients. This suggests that heterogeneous fractional Brownian motion is an appropriate model for both fast and slow moving endosomes. This article is part of a Special Issue entitled: "Recent Advances In Single-Particle Tracking: Experiment and Analysis" edited by Janusz Szwabiński and Aleksander Weron.
We study the dynamics of membrane vesicle motor transport into dendritic spines, which are bulbous intracellular compartments in neurons that play a key role in transmitting signals between neurons. We consider the stochastic analog of the vesicle transport model in [Park and Fai, The Dynamics of Vesicles Driven Into Closed Constrictions by Molecular Motors. Bull. Math. Biol. 82, 141 (2020)]. The stochastic version, which may be considered as an agent-based model, relies mostly on the action of individual myosin motors to produce vesicle motion. To aid in our analysis, we coarse-grain this agent-based model using a master equation combined with a partial differential equation describing the probability of local motor positions. We confirm through convergence studies that the coarse-graining captures the essential features of bistability in velocity (observed in experiments) and waiting-time distributions to switch between steady-state velocities. Interestingly, these results allow us to reformulate the translocation problem in terms of the mean first passage time for a run-and-tumble particle moving on a finite domain with absorbing boundaries at the two ends. We conclude by presenting numerical and analytical calculations of vesicle translocation.
Shawn D. Ryan, Zachary McCarthy, Mykhailo Potomkin
Many cellular processes rely on the cell's ability to transport material to and from the nucleus. Networks consisting of many microtubules and actin filaments are key to this transport. Recently, the inhibition of intracellular transport has been implicated in neurodegenerative diseases such as Alzheimer's disease and Amyotrophic Lateral Sclerosis (ALS). Furthermore, microtubules may contain so-called defective regions where motor protein velocity is reduced due to accumulation of other motors and microtubule associated proteins. In this work, we propose a new mathematical model describing the motion of motor proteins on microtubules which incorporate a defective region. We take a mean-field approach derived from a first principle lattice model to study motor protein dynamics and density profiles. In particular, given a set of model parameters we obtain a closed-form expression for the equilibrium density profile along a given microtubule. We then verify the analytic results using mathematical analysis on the discrete model and Monte Carlo simulations. This work will contribute to the fundamental understanding of inhomogeneous microtubules providing insight into microscopic interactions that may result in the onset of neurodegenerative diseases. Our results for inhomogeneous microtubules are consistent with prior work studying the homogeneous case.
Stefan Schuster, Jan Ewald, Thomas Dandekar
et al.
In host-pathogen interactions, often the host (attacked organism) defends itself by some toxic compound and the parasite, in turn, responds by producing an enzyme that inactivates that compound. In some cases, the host can respond by producing an inhibitor of that enzyme, which can be considered as a counter-counter defence. An example is provided by cephalosporins, beta-lactamases and clavulanic acid (an inhibitor of beta-lactamases). Here, we tackle the question under which conditions it pays, during evolution, to establish a counter-counter defence rather than to intensify or widen the defence mechanisms. We establish a mathematical model describing this phenomenon, based on enzyme kinetics for competitive inhibition. We use an objective function based on Haber's rule, which says that the toxic effect is proportional to the time integral of toxin concentration. The optimal allocation of defence and counter-counter defence can be calculated in an analytical way despite the nonlinearity in the underlying differential equation. The calculation provides a threshold value for the dissociation constant of the inhibitor. Only if the inhibition constant is below that threshold, that is, in the case of strong binding of the inhibitor, it pays to have a counter-counter defence. This theoretical prediction accounts for the observation that not for all defence mechanisms, a counter-counter defence exists. Our results should be of interest for computing optimal mixtures of beta-lactam antibiotics and beta-lactamase inhibitors such as sulbactam, as well as for plant-herbivore and other molecular-ecological interactions and to fight antibiotic resistance in general.
How cells regulate the number of organelles is a fundamental question in cell biology. While decades of experimental work have uncovered four fundamental processes that regulate organelle biogenesis, namely, de novo synthesis, fission, fusion and decay, a comprehensive understanding of how these processes together control organelle abundance remains elusive. Recent fluorescence microscopy experiments allow for the counting of organelles at the single-cell level. These measurements provide information about the cell-to-cell variability in organelle abundance in addition to the mean level. Motivated by such measurements, we build upon a recent study and analyze a general stochastic model of organelle biogenesis. We compute the exact analytical expressions for the probability distribution of organelle numbers, their mean, and variance across a population of single cells. It is shown that different mechanisms of organelle biogenesis lead to distinct signatures in the distribution of organelle numbers which allows us to discriminate between these various mechanisms. By comparing our theory against published data for peroxisome abundance measurements in yeast, we show that a widely believed model of peroxisome biogenesis that involves de novo synthesis, fission, and decay is inadequate in explaining the data. Also, our theory predicts bimodality in certain limits of the model. Overall, the framework developed here can be harnessed to gain mechanistic insights into the process of organelle biogenesis.
Thiruvallur R. Gowrishankar, Julie V. Stern, Kyle C. Smith
et al.
Extraordinarily large but short electric field pulses are reported by many experiments to cause bipolar cancellation (BPC). This unusual cell response occurs if a first pulse is followed by a second pulse with opposite polarity. Possibly universal, BPC presently lacks a mechanistic explanation. Multiple versions of the "standard model" of cell electroporation (EP) fail to account for BPC. Here we show, for the first time, how an extension of the standard model can account for a key experimental observation that essentially defines BPC: the amount of a tracer that enters a cell, and how tracer influx can be decreased by the second part of a bipolar pulse. The extended model can also account for the recovery of BPC wherein the extent of BPC is diminished if the spacing between the first and second pulses is increased. Our approach is reverse engineering, meaning that we identify and introduce an additional biophysical mechanism that allows pore transport to change. We hypothesize that occluding molecules from outside the membrane enter or relocate within a pore. Significantly, the additional mechanism is fundamental and general, involving a combination of partitioning and hindrance. Molecules near the membrane can enter pores to block transport of tracer molecules while still passing small ions (+/- 1) that govern electrical behavior. Accounting for such behavior requires an extension of the standard model.
In permissive environments, E. coli can double its dry mass every 21 minutes. During this time, ribosomes, RNA polymerases, and the proteome are all doubled. Yet, the question of how to relate bacterial doubling time to other biologically relevant time scales in the growth process remains illusive, due to the complex temporal nesting pattern of these processes. In particular, the relation between the cell's doubling time and the ribosome assembly time is not known. Here we develop a model that connects growth rate to ribosome assembly time and show that the existence of a self-assembly step increases the overall growth rate, because during ribosome self-assembly existing ribosomes can start a new round of reproduction, by making a new batch of ribosomal proteins prior to the completion of the previous round. This overlapping of ribosome reproduction cycles increases growth rate beyond the serial-limit that is typically assumed to hold. Using recent data from ribosome profiling and well known measurements of the average translation rate, rigid bounds on the in-vivo ribosome self-assembly time are set, which are robust to the assumptions regarding the biological noises involved. At 21 minutes doubling time, the ribosome assembly time is found to be approximately 6 minutes --- three fold larger than the common estimate. We further use our model to explain the detrimental effect of a recently discovered ribosome assembly inhibitor drug, and predict the effect of limiting the expression of ribosome assembly chaperons on the overall growth rate.
Martin Carballo-Pacheco, Jonathan Desponds, Tatyana Gavrilchenko
et al.
Cells need to reliably sense external ligand concentrations to achieve various biological functions such as chemotaxis or signaling. The molecular recognition of ligands by surface receptors is degenerate in many systems leading to crosstalk between different receptors. Crosstalk is often thought of as a deviation from optimal specific recognition, as the binding of non-cognate ligands can interfere with the detection of the receptor's cognate ligand, possibly leading to a false triggering of a downstream signaling pathway. Here we quantify the optimal precision of sensing the concentrations of multiple ligands by a collection of promiscuous receptors. We demonstrate that crosstalk can improve precision in concentration sensing and discrimination tasks. To achieve superior precision, the additional information about ligand concentrations contained in short binding events of the non-cognate ligand should be exploited. We present a proofreading scheme to realize an approximate estimation of multiple ligand concentrations that reaches a precision close to the derived optimal bounds. Our results help rationalize the observed ubiquity of receptor crosstalk in molecular sensing.
Emma M. Keizer, Bjorn Bastian, Robert W. Smith
et al.
It is well known that the kinetics of an intracellular biochemical network is stochastic. This is due to intrinsic noise arising from the random timing of biochemical reactions in the network as well as due to extrinsic noise stemming from the interaction of unknown molecular components with the network and from the cell's changing environment. While there are many methods to study the effect of intrinsic noise on the system dynamics, few exist to study the influence of both types of noise. Here we show how one can extend the conventional linear-noise approximation to allow for the rapid evaluation of the molecule numbers statistics of a biochemical network influenced by intrinsic noise and by slow lognormally distributed extrinsic noise. The theory is applied to simple models of gene regulatory networks and its validity confirmed by comparison with exact stochastic simulations. In particular we show how extrinsic noise modifies the dependence of the variance of the molecule number fluctuations on the rate constants, the mutual information between input and output signalling molecules and the robustness of feed-forward loop motifs.
Christopher E. Miles, Sean D. Lawley, James P. Keener
We propose and analyze a mathematical model of cargo transport by non-processive molecular motors. In our model, the motors change states by random discrete events (corresponding to stepping and binding/unbinding), while the cargo position follows a stochastic differential equation (SDE) that depends on the discrete states of the motors. The resulting system for the cargo position is consequently an SDE that randomly switches according to a Markov jump process governing motor dynamics. To study this system we (1) cast the cargo position in a renewal theory framework and generalize the renewal reward theorem and (2) decompose the continuous and discrete sources of stochasticity and exploit a resulting pair of disparate timescales. With these mathematical tools, we obtain explicit formulas for experimentally measurable quantities, such as cargo velocity and run length. Analyzing these formulas then yields some predictions regarding so-called non-processive clustering, the phenomenon that a single motor cannot transport cargo, but that two or more motors can. We find that having motor stepping, binding, and unbinding rates depend on the number of bound motors, due to geometric effects, is necessary and sufficient to explain recent experimental data on non-processive motors.
Information processing at the molecular scale is limited by thermal fluctuations. This can cause undesired consequences in copying information since thermal noise can lead to errors that can compromise the functionality of the copy. For example, a high error rate during DNA duplication can lead to cell death. Given the importance of accurate copying at the molecular scale, it is fundamental to understand its thermodynamic features. In this paper, we derive a universal expression for the copy error as a function of entropy production and {\cred work dissipated by the system during wrong incorporations}. Its derivation is based on the second law of thermodynamics, hence its validity is independent of the details of the molecular machinery, be it any polymerase or artificial copying device. Using this expression, we find that information can be copied in three different regimes. In two of them, work is dissipated to either increase or decrease the error. In the third regime, the protocol extracts work while correcting errors, reminiscent of a Maxwell demon. As a case study, we apply our framework to study a copy protocol assisted by kinetic proofreading, and show that it can operate in any of these three regimes. We finally show that, for any effective proofreading scheme, error reduction is limited by the chemical driving of the proofreading reaction.
We introduce a model for microtubule mechanics containing lateral bonds between dimers in neighboring protofilaments, bending rigidity of dimers, and repulsive interactions between protofilaments modeling steric constraints to investigate the influence of mechanical forces on hydrolysis and catastrophes. We use the allosteric dimer model, where tubulin dimers are characterized by an equilibrium bending angle, which changes from $0^\circ$ to $22^\circ$ by hydrolysis of a dimer. This also affects the lateral interaction and bending energies and, thus, the mechanical equilibrium state of the microtubule. As hydrolysis gives rise to conformational changes in dimers, mechanical forces also influence the hydrolysis rates by mechanical energy changes modulating the hydrolysis rate. The interaction via the microtubule mechanics then gives rise to correlation effects in the hydrolysis dynamics, which have not been taken into account before. Assuming a dominant influence of mechanical energies on hydrolysis rates, we investigate the most probable hydrolysis pathways both for vectorial and random hydrolysis. Investigating the stability with respect to lateral bond rupture, we identify initiation configurations for catastrophes along the hydrolysis pathways and values for a lateral bond rupture force. If we allow for rupturing of lateral bonds between dimers in neighboring protofilaments above this threshold force, our model exhibits avalanche-like catastrophe events.
In the past one hundred years, deterministic rate equations have been successfully used to infer enzyme-catalysed reaction mechanisms and to estimate rate constants from reaction kinetics experiments conducted in vitro. In recent years, sophisticated experimental techniques have been developed that allow the measurement of enzyme- catalysed and other biopolymer-mediated reactions inside single cells at the single molecule level. Time course data obtained by these methods are considerably noisy because molecule numbers within cells are typically quite small. As a consequence, the interpretation and analysis of single cell data requires stochastic methods, rather than deterministic rate equations. Here we concisely review both experimental and theoretical techniques which enable single molecule analysis with particular emphasis on the major developments in the field of theoretical stochastic enzyme kinetics, from its inception in the mid-twentieth century to its modern day status. We discuss the differences between stochastic and deterministic rate equation models, how these depend on enzyme molecule numbers and substrate inflow into the reaction compartment and how estimation of rate constants from single cell data is possible using recently developed stochastic approaches.