Analysis of dim-light responses in rod and cone photoreceptors with altered calcium kinetics.

Rod and cone photoreceptors in the retina of vertebrates are the primary sensory neurons underlying vision. They convert light into an electrical current using a signal transduction pathway that depends on Ca[Formula: see text] feedback. It is known that manipulating the Ca[Formula: see text] kinetics affects the response shape and the photoreceptor sensitivity, but a precise quantification of these effects remains unclear. We have approached this task in mouse retina by combining numerical simulations with mathematical analysis. We consider a parsimonious phototransduction model that incorporates negative Ca[Formula: see text] feedback onto the synthesis of cyclic GMP, and fast buffering reactions to alter the Ca[Formula: see text] kinetics. We derive analytic results for the photoreceptor functioning in sufficiently dim light conditions depending on the photoreceptor type. We exploit these results to obtain conceptual and quantitative insight into how response waveform and amplitude depend on the underlying biophysical processes and the Ca[Formula: see text] feedback. With a low amount of buffering, the Ca[Formula: see text] concentration changes in proportion to the current, and responses to flashes of light are monophasic. With more buffering, the change in the Ca[Formula: see text] concentration becomes delayed with respect to the current, which gives rise to a damped oscillation and a biphasic waveform. This shows that biphasic responses are not necessarily a manifestation of slow buffering reactions. We obtain analytic approximations for the peak flash amplitude as a function of the light intensity, which shows how the photoreceptor sensitivity depends on the biophysical parameters. Finally, we study how changing the extracellular Ca[Formula: see text] concentration affects the response.

Ca2+-activated Cl- current ensures robust and reliable signal amplification in vertebrate olfactory receptor neurons.

Activation of most primary sensory neurons results in transduction currents that are carried by cations. One notable exception is the vertebrate olfactory receptor neuron (ORN), where the transduction current is carried largely by the anion [Formula: see text] However, it remains unclear why ORNs use an anionic current for signal amplification. We have sought to provide clarification on this topic by studying the so far neglected dynamics of [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] in the small space of olfactory cilia during an odorant response. Using computational modeling and simulations we compared the outcomes of signal amplification based on either [Formula: see text] or [Formula: see text] currents. We found that amplification produced by [Formula: see text] influx instead of a [Formula: see text] efflux is problematic for several reasons: First, the [Formula: see text] current amplitude varies greatly, depending on mucosal ion concentration changes. Second, a [Formula: see text] current leads to a large increase in the ciliary [Formula: see text] concentration during an odorant response. This increase inhibits and even reverses [Formula: see text] clearance by [Formula: see text] exchange, which is essential for response termination. Finally, a [Formula: see text] current increases the ciliary osmotic pressure, which could cause swelling to damage the cilia. By contrast, a transduction pathway based on [Formula: see text] efflux circumvents these problems and renders the odorant response robust and reliable.

Analysis of waveform and amplitude of mouse rod and cone flash responses.

KEY POINTS: Most vertebrate eyes have rod and cone photoreceptors, which use a signal transduction pathway consisting of many biological processes to transform light into an electrical response. We dissect and quantify the contribution of each of these processes to the photoreceptor light response by using a novel method of analysis that provides an analytical solution for the entire time course of the dim-flash light response. We find that the shape of the light response is exclusively controlled by deactivation parameters. Activation parameters scale this shape and alter the response amplitude. We show that the rising phase of the response depends on Ca2+ feedback, and we identify the deactivation parameters that control the recovery phase of the response. We devise new methods to extract values for deactivation and activation parameters from a separate analysis of response shape and response amplitude. ABSTRACT: Vertebrate eyes have rod and cone photoreceptors, which use a complex transduction pathway comprising many biological processes to transform the absorption of light into an electrical response. A fundamental question in sensory transduction is how these processes contribute to the response. To study this question, we use a well-accepted phototransduction model, which we analyse with a novel method based on the log transform of the current. We derive an analytical solution that describes the entire time course of the photoreceptor response to dim flashes of light. We use this solution to dissect and quantify the contribution of each process to the response. We find that the entire dim-flash response is proportional to the flash intensity. By normalizing responses to unit amplitude, we define a waveform that is independent of the light intensity and characterizes the invariant shape of dim-flash responses. We show that this waveform is exclusively determined by deactivation rates; activation rates only scale the waveform and affect the amplitude. This analysis corrects a previous assumption that the rising phase is determined entirely by activation rates. We further show that the rising phase depends on Ca2+ feedback to the cyclase, contrary to current belief. We identify the deactivation rates that control the recovery phase of the response, and we devise new methods to extract activation and deactivation rates from an analysis of response shape and response amplitude. In summary, we provide a comprehensive understanding of how the various transduction processes produce the cellular response.

A kinetic analysis of mouse rod and cone photoreceptor responses.

KEY POINTS: Most vertebrate eyes have rods for dim-light vision and cones for brighter light and higher temporal sensitivity. Rods evolved from cone-like precursors through expression of different transduction genes or the same genes at different expression levels, but we do not know which molecular differences were most important. We approached this problem by analysing rod and cone responses with the same model but with different values for model parameters. We showed that, in addition to outer-segment volume, the most important differences between rods and cones are: (1) decreased transduction gain, reflecting smaller amplification in the G-protein cascade; (2) a faster rate of turnover of the second messenger cGMP in darkness; and (3) an accelerated rate of decay of the effector enzyme phosphodiesterase and perhaps also of activated visual pigment. We believe our analysis has identified the principal alterations during evolution responsible for the duplex retina. ABSTRACT: Most vertebrates have rod and cone photoreceptors, which differ in their sensitivity and response kinetics. We know that rods evolved from cone-like precursors through the expression of different transduction genes or the same genes at different levels, but we do not know which molecular differences were most important. We have approached this problem in mouse retina by analysing the kinetic differences between rod flash responses and recent voltage-clamp recordings of cone flash responses, using a model incorporating the principal features of photoreceptor transduction. We apply a novel method of analysis using the log-transform of the current, and we ask which of the model's dynamic parameters need be changed to transform the flash response of a rod into that of a cone. The most important changes are a decrease in the gain of the response, reflecting a reduction in amplification of the transduction cascade; an increase in the rate of turnover of cGMP in darkness; and an increase in the rate of decay of activated phosphodiesterase, with perhaps also an increase in the rate of decay of light-activated visual pigment. Although we cannot exclude other differences, and in particular alterations in the Ca2+ economy of the photoreceptors, we believe that we have identified the kinetic parameters principally responsible for the differences in the flash responses of the two kinds of photoreceptors, which were likely during evolution to have resulted in the duplex retina.

The PDE6 mutation in the rd10 retinal degeneration mouse model causes protein mislocalization and instability and promotes cell death through increased ion influx.

The retinal degeneration model rd10 contains a missense mutation of the catalytic PDE6 ß subunit, which hydrolyzes cGMP in response to light. This model produces cell death more slowly than others caused by PDE6 loss of function, making it of particular interest for studying potential therapeutics. We used morphology, biochemistry, and single-cell physiology to examine the mechanism of rd10 degeneration. Our results show that the mutation produces no alteration of Pde6b RNA but does dramatically decrease maximal and basal PDE6 activity, apparently caused by a decrease in protein stability and transport. The enzymatic properties of the remaining mutant PDE6 appear to be nearly normal. We demonstrate that an increase in free cGMP, which would result from decreased PDE6 activity and serve to increase opening of the cGMP-gated channels and calcium influx, is an underlying cause of cell death: degeneration of rd10/Cngb1-/- double mutants is slower than the parent rd10 line. Paradoxically, degeneration in rd10/Cngb1-/- is also slower than in Cngb1-/- This rescue is correlated with a lowering of cGMP content in Cngb1-/- retinas and suggests that it may be caused by mislocalization of active PDE6. Single-cell recordings from rd10 rods show that the rates of rise and decay of the response are significantly slower; simulations indicate that these changes are primarily the result of the decrease in PDE6 concentration and rod collecting area. Together, these results provide insights into the complex mechanisms that underlie rd10-mediated retinal degeneration and a cautionary note for analysis of therapeutic interventions.

How rods respond to single photons: Key adaptations of a G-protein cascade that enable vision at the physical limit of perception.

Rod photoreceptors are among the most sensitive light detectors in nature. They achieve their remarkable sensitivity across a wide variety of species through a number of essential adaptations: a specialized cellular geometry, a G-protein cascade with an unusually stable receptor molecule, a low-noise transduction mechanism, a nearly perfect effector enzyme, and highly evolved mechanisms of feedback control and receptor deactivation. Practically any change in protein expression, enzyme activity, or feedback control can be shown to impair photon detection, either by decreasing sensitivity or signal-to-noise ratio, or by reducing temporal resolution. Comparison of mammals to amphibians suggests that rod outer-segment morphology and the molecules and mechanism of transduction may have evolved together to optimize light sensitivity in darkness, which culminates in the extraordinary ability of these cells to respond to single photons at the ultimate limit of visual perception.

Computational and mathematical methods for morphogenetic gradient analysis, boundary formation and axonal targeting.

Morphogenesis and axonal targeting are key processes during development that depend on complex interactions at molecular, cellular and tissue level. Mathematical modeling is essential to bridge this multi-scale gap in order to understand how the emergence of large structures is controlled at molecular level by interactions between various signaling pathways. We summarize mathematical modeling and computational methods for time evolution and precision of morphogenetic gradient formation. We discuss tissue patterning and the formation of borders between regions labeled by different morphogens. Finally, we review models and algorithms that reveal the interplay between morphogenetic gradients and patterned activity for axonal pathfinding and the generation of the retinotopic map in the visual system.

Detection of single photons by toad and mouse rods.

Amphibian and mammalian rods can both detect single photons of light even though they differ greatly in physical dimensions, mammalian rods being much smaller in diameter than amphibian rods. To understand the changes in physiology and biochemistry required by such large differences in outer segment geometry, we developed a computational approach, taking into account the spatial organization of the outer segment divided into compartments, together with molecular dynamics simulations of the signaling cascade. We generated simulations of the single-photon response together with intrinsic background fluctuations in toad and mouse rods. Combining this computational approach with electrophysiological data from mouse rods, we determined key biochemical parameters. On average around one phosphodiesterase (PDE) molecule is spontaneously active per mouse compartment, similar to the value for toad, which is unexpected due to the much smaller diameter in mouse. A larger number of spontaneously active PDEs decreases dark noise, thereby improving detection of single photons; it also increases cGMP turnover, which accelerates the decay of the light response. These constraints explain the higher PDE density in mammalian compared with amphibian rods that compensates for the much smaller diameter of mammalian disks. We further find that the rate of cGMP hydrolysis by light-activated PDE is diffusion limited, which is not the case for spontaneously activated PDE. As a consequence, in the small outer segment of a mouse rod only a few activated PDEs are sufficient to generate a signal that overcomes noise, which permits a shorter lifetime of activated rhodopsin and greater temporal resolution.

Engrailed homeoprotein recruits the adenosine A1 receptor to potentiate ephrin A5 function in retinal growth cones.

Engrailed 1 and engrailed 2 homeoprotein transcription factors (collectively Engrailed) display graded expression in the chick optic tectum where they participate in retino-tectal patterning. In vitro, extracellular Engrailed guides retinal ganglion cell (RGC) axons and synergises with ephrin A5 to provoke the collapse of temporal growth cones. In vivo disruption of endogenous extracellular Engrailed leads to misrouting of RGC axons. Here we characterise the signalling pathway of extracellular Engrailed. Our results show that Engrailed/ephrin A5 synergy in growth cone collapse involves adenosine A1 receptor activation after Engrailed-dependent ATP synthesis, followed by ATP secretion and hydrolysis to adenosine. This is, to our knowledge, the first evidence for a role of the adenosine A1 receptor in axon guidance. Based on these results, together with higher expression of the adenosine A1 receptor in temporal than nasal growth cones, we propose a computational model that illustrates how the interaction between Engrailed, ephrin A5 and adenosine could increase the precision of the retinal projection map.

Facilitated diffusion framework for transcription factor search with conformational changes.

Cellular responses often require the fast activation or repression of specific genes, which depends on transcription factors (TFs) that have to quickly find the promoters of these genes within a large genome. TFs search for their DNA promoter target by alternating between bulk diffusion and sliding along the DNA, a mechanism known as facilitated diffusion. We study a facilitated diffusion framework with switching between three search modes: a bulk mode and two sliding modes triggered by conformational changes between two protein conformations. In one conformation (search mode) the TF interacts unspecifically with the DNA backbone resulting in fast sliding. In the other conformation (recognition mode) it interacts specifically and strongly with DNA base pairs leading to slow displacement. From the bulk, a TF associates with the DNA at a random position that is correlated with the previous dissociation point, which implicitly is a function of the DNA structure. The target affinity depends on the conformation. We derive exact expressions for the mean first passage time (MFPT) to bind to the promoter and the conditional probability to bind before detaching when arriving at the promoter site. We systematically explore the parameter space and compare various search scenarios. We compare our results with experimental data for the dimeric Lac repressor search in E. coli bacteria. We find that a coiled DNA conformation is absolutely necessary for a fast MFPT. With frequent spontaneous conformational changes, a fast search time is achieved even when a TF becomes immobilized in the recognition state due to the specific bindings. We find a MFPT compatible with experimental data in presence of a specific TF-DNA interaction energy that has a Gaussian distribution with a large variance.

Dissection of a Krox20 positive feedback loop driving cell fate choices in hindbrain patterning.

Although feedback loops are essential in development, their molecular implementation and precise functions remain elusive. Using enhancer knockout in mice, we demonstrate that a direct, positive autoregulatory loop amplifies and maintains the expression of Krox20, a transcription factor governing vertebrate hindbrain segmentation. By combining quantitative data collected in the zebrafish with biophysical modelling that accounts for the intrinsic stochastic molecular dynamics, we dissect the loop at the molecular level. We find that it underpins a bistable switch that turns a transient input signal into cell fate commitment, as we observe in single cell analyses. The stochasticity of the activation process leads to a graded input-output response until saturation is reached. Consequently, the duration and strength of the input signal controls the size of the hindbrain segments by modulating the distribution between the two cell fates. Moreover, segment formation is buffered from severe variations in input level. Finally, the progressive extinction of Krox20 expression involves a destabilization of the loop by repressor molecules. These mechanisms are of general significance for cell type specification and tissue patterning.

Astroglial networks scale synaptic activity and plasticity.

Astrocytes dynamically interact with neurons to regulate synaptic transmission. Although the gap junction proteins connexin 30 (Cx30) and connexin 43 (Cx43) mediate the extensive network organization of astrocytes, their role in synaptic physiology is unknown. Here we show, by inactivating Cx30 and Cx43 genes, that astroglial networks tone down hippocampal synaptic transmission in CA1 pyramidal neurons. Gap junctional networking facilitates extracellular glutamate and potassium removal during synaptic activity through modulation of astroglial clearance rate and extracellular space volume. This regulation limits neuronal excitability, release probability, and insertion of postsynaptic AMPA receptors, silencing synapses. By controlling synaptic strength, connexins play an important role in synaptic plasticity. Altogether, these results establish connexins as critical proteins for extracellular homeostasis, important for the formation of functional synapses.

Synapsin 2a tetramerisation selectively controls the presynaptic nanoscale organisation of reserve synaptic vesicles.

Neurotransmitter release relies on the regulated fusion of synaptic vesicles (SVs) that are tightly packed within the presynaptic bouton of neurons. The mechanism by which SVs are clustered at the presynapse, while preserving their ability to dynamically recycle to support neuronal communication, remains unknown. Synapsin 2a (Syn2a) tetramerization has been suggested as a potential clustering mechanism. Here, we used Dual-pulse sub-diffractional Tracking of Internalised Molecules (DsdTIM) to simultaneously track single SVs from the recycling and the reserve pools, in live hippocampal neurons. The reserve pool displays a lower presynaptic mobility compared to the recycling pool and is also present in the axons. Triple knockout of Synapsin 1-3 genes (SynTKO) increased the mobility of reserve pool SVs. Re-expression of wild-type Syn2a (Syn2aWT), but not the tetramerization-deficient mutant K337Q (Syn2aK337Q), fully rescued these effects. Single-particle tracking revealed that Syn2aK337QmEos3.1 exhibited altered activity-dependent presynaptic translocation and nanoclustering. Therefore, Syn2a tetramerization controls its own presynaptic nanoclustering and thereby contributes to the dynamic immobilisation of the SV reserve pool.

Data-driven multiscale dynamical framework to control a pandemic evolution with non-pharmaceutical interventions.

Before the availability of vaccines, many countries have resorted multiple times to drastic social restrictions to prevent saturation of their health care system, and to regain control over an otherwise exponentially increasing COVID-19 pandemic. With the advent of data-sharing, computational approaches are key to efficiently control a pandemic with non-pharmaceutical interventions (NPIs). Here we develop a data-driven computational framework based on a time discrete and age-stratified compartmental model to control a pandemic evolution inside and outside hospitals in a constantly changing environment with NPIs. Besides the calendrical time, we introduce a second time-scale for the infection history, which allows for non-exponential transition probabilities. We develop inference methods and feedback procedures to successively recalibrate model parameters as new data becomes available. As a showcase, we calibrate the framework to study the pandemic evolution inside and outside hospitals in France until February 2021. We combine national hospitalization statistics from governmental websites with clinical data from a single hospital to calibrate hospitalization parameters. We infer changes in social contact matrices as a function of NPIs from positive testing and new hospitalization data. We use simulations to infer hidden pandemic properties such as the fraction of infected population, the hospitalisation probability, or the infection fatality ratio. We show how reproduction numbers and herd immunity levels depend on the underlying social dynamics.

Syntaxin1A lateral diffusion reveals transient and local SNARE interactions.

At the synapse, vesicles stably dock at the active zone. However, in cellular membranes, proteins undergo a diffusive motion. It is not known how the motion of membrane proteins involved in vesicle exocytosis is compatible with both vesicle docking and the dynamic remodeling of the plasma membrane imposed by cycles of exocytosis and endocytosis. To address this question, we studied the motion of the presynaptic membrane protein syntaxin1A at both the population and single-molecule levels in primary cultures of rat spinal cord neurons. Syntaxin1A was rapidly exchanged between synaptic and extrasynaptic regions. Changes in syntaxin1A mobility were associated with interactions related to the formation of the exocytotic complex. Finally, we propose a reaction-diffusion model reconciling the observed diffusive properties of syntaxin at the population level and at the molecular level. This work allows us to describe the diffusive behavior and kinetics of interactions between syntaxin1A and its partners that lead to its transient stabilization at the synapse.

Gated narrow escape time for molecular signaling.

The mean time for a diffusing ligand to activate a target protein located on the surface of a microdomain can regulate cellular signaling. When the ligand switches between various states induced by chemical interactions or conformational changes, while target activation occurs in only one state, this activation time is affected. We investigate this dynamics using new equations for the sojourn times spent in each state. For two states, we obtain exact solutions in dimension one, and asymptotic ones confirmed by Brownian simulations in dimension 3. We find that the activation time is quite sensitive to changes of the switching rates, which can be used to modulate signaling. Interestingly, our analysis reveals that activation can be fast although the ligand spends most of the time "hidden" in the nonactivating state. Finally, we obtain a new formula for the narrow escape time in the presence of switching.

Narrow escape time to a structured target located on the boundary of a microdomain.

The forward binding rate of chemical reactions is the reciprocal of the mean time for a Brownian molecule to hit its molecular target. When the target is embedded in the surface of a microdomain, this time is known as the narrow escape time, and it has been computed for various geometries. However, for large targets that extend from the surface far into the cytosol the classical computations do not apply and new ones are needed. In this work we generalize the narrow escape time formula to a three-dimensional spine attached to the boundary. We treat in detail the case of an ellipsoidal spine and validate our analysis by Brownian simulations. Finally, we compute the narrow escape time when the spine is uniformly covered by small partially absorbing disks and estimate the homogenized trapping rate of such a patchy surface.

Diffusion in narrow domains and application to phototransduction.

The mean time for a Brownian particle to find a small target inside a narrow domain is a key parameter for many chemical reactions occurring in cellular microstructures. Although current estimations are given for a large class of domains, they cannot be used for narrow domains often encountered in cellular biology, such as the synaptic cleft, narrow compartments in the outer segment of vertebrate photoreceptors, or neuron-glia contact. We compute here the mean time for a Brownian particle to hit a small target placed on the surface of a narrow cylinder. We then use this result to estimate the rate constant of cyclic-GMP (cGMP) hydrolysis by the activated enzyme phosphodiesterase (PDE) in the narrow microdomains that build up the outer segment of a rod photoreceptor. By controlling the cGMP concentration, PDE activity is at the basis of the early photoresponse chemical reaction cascade. Our approach allows us to compute the cGMP rate constant as a function of biophysical parameters.

Asymmetric ephaptic inhibition between compartmentalized olfactory receptor neurons.

In the Drosophila antenna, different subtypes of olfactory receptor neurons (ORNs) housed in the same sensory hair (sensillum) can inhibit each other non-synaptically. However, the mechanisms underlying this underexplored form of lateral inhibition remain unclear. Here we use recordings from pairs of sensilla impaled by the same tungsten electrode to demonstrate that direct electrical ("ephaptic") interactions mediate lateral inhibition between ORNs. Intriguingly, within individual sensilla, we find that ephaptic lateral inhibition is asymmetric such that one ORN exerts greater influence onto its neighbor. Serial block-face scanning electron microscopy of genetically identified ORNs and circuit modeling indicate that asymmetric lateral inhibition reflects a surprisingly simple mechanism: the physically larger ORN in a pair corresponds to the dominant neuron in ephaptic interactions. Thus, morphometric differences between compartmentalized ORNs account for highly specialized inhibitory interactions that govern information processing at the earliest stages of olfactory coding.

The dynamics of phosphodiesterase activation in rods and cones.

Phototransduction starts with the activation of a rhodopsin (respectively, coneopsin) molecule, located in the outer segment of rod (respectively, cone) photoreceptors. The subsequent amplification pathway proceeds via the G-protein transducin to the activation of phosphodiesterase (PDE), a G-protein coupled effector enzyme. In this article, we study the dynamics of PDE activation by constructing a Markov model that is based on the underlying chemical reactions including multiple rhodopsin phosphorylations. We derive explicit equations for the mean and the variance of activated PDE. Our analysis reveals that a low rhodopsin lifetime variance is neither necessary nor sufficient to achieve reliable PDE activation. The numerical simulations show that during the rising phase the variability of PDE activation is much lower compared to the recovery phase, and this property depends crucially on the transducin activation rates. Furthermore, we find that the dynamics of the activation process greatly differs depending on whether rhodopsin or PDE deactivation limits the recovery of the photoresponse. Finally, our simulations for cones show that only very few PDEs are activated by an excited photopigment, which might explain why in S-cones no single photon response can be observed.