Analysis of noisy transient signals based on Gaussian process regression.

Dynamic systems such as cells or tissues generate, either spontaneously or in response to stimuli, transient signals that carry information about the system. Characterization of recorded transients is often hampered by a low signal-to-noise ratio (SNR). Reduction of the noise by filtering has limited use due to partial signal distortion. Occasionally, transients can be approximated by a mathematical function, but such a function may not hold correctly if recording conditions change. We introduce here the model-independent approximation method for general noisy transient signals based on the Gaussian process regression. The method was implemented in the software TransientAnalyzer, which detects transients in a record, finds their best approximation by the Gaussian process, constructs a surrogate spline function, and estimates specified signal parameters. The method and software were tested on a cellular model of the calcium concentration transient corrupted by various SNR levels and recorded at a low sampling frequency. Statistical analysis of the model data sets provided the error of estimation <7.5% and the coefficient of variation of estimates <17% for peak SNR = 5. The performance of Gaussian process regression on signals of diverse experimental origin was even better than fitting by a function. The software and its description are available on GitHub.

Impaired Binding to Junctophilin-2 and Nanostructural Alteration in CPVT Mutation.

[Figure: see text].

In silico simulations reveal that RYR distribution affects the dynamics of calcium release in cardiac myocytes.

The dyads of cardiac myocytes contain ryanodine receptors (RYRs) that generate calcium sparks upon activation. To test how geometric factors of RYR distribution contribute to the formation of calcium sparks, which cannot be addressed experimentally, we performed in silico simulations on a large set of models of calcium release sites (CRSs). Our models covered the observed range of RYR number, density, and spatial arrangement. The calcium release function of CRSs was modeled by RYR openings, with an open probability dependent on concentrations of free Ca2+ and Mg2+ ions, in a rapidly buffered system, with a constant open RYR calcium current. We found that simulations of spontaneous sparks by repeatedly opening one of the RYRs in a CRS produced three different types of calcium release events (CREs) in any of the models. Transformation of simulated CREs into fluorescence signals yielded calcium sparks with characteristics close to the observed ones. CRE occurrence varied broadly with the spatial distribution of RYRs in the CRS but did not consistently correlate with RYR number, surface density, or calcium current. However, it correlated with RYR coupling strength, defined as the weighted product of RYR vicinity and calcium current, so that CRE characteristics of all models followed the same state-response function. This finding revealed the synergy between structure and function of CRSs in shaping dyad function. Lastly, rearrangements of RYRs simulating hypothetical experiments on splitting and compaction of a dyad revealed an increased propensity to generate spontaneous sparks and an overall increase in calcium release in smaller and more compact dyads, thus underlying the importance and physiological role of RYR arrangement in cardiac myocytes.

Magnesium Ions Moderate Calcium-Induced Calcium Release in Cardiac Calcium Release Sites by Binding to Ryanodine Receptor Activation and Inhibition Sites.

Ryanodine receptor channels at calcium release sites of cardiac myocytes operate on the principle of calcium-induced calcium release. In vitro experiments revealed competition of Ca2+ and Mg2+ in the activation of ryanodine receptors (RyRs) as well as inhibition of RyRs by Mg2+. The impact of RyR modulation by Mg2+ on calcium release is not well understood due to the technical limitations of in situ experiments. We turned instead to an in silico model of a calcium release site (CRS), based on a homotetrameric model of RyR gating with kinetic parameters determined from in vitro measurements. We inspected changes in the activity of the CRS model in response to a random opening of one of 20 realistically distributed RyRs, arising from Ca2+/Mg2+ interactions at RyR channels. Calcium release events (CREs) were simulated at a range of Mg2+-binding parameters at near-physiological Mg2+ and ATP concentrations. Facilitation of Mg2+ binding to the RyR activation site inhibited the formation of sparks and slowed down their activation. Impeding Mg-binding to the RyR activation site enhanced spark formation and speeded up their activation. Varying Mg2+ binding to the RyR inhibition site also dramatically affected calcium release events. Facilitation of Mg2+ binding to the RyR inhibition site reduced the amplitude, relative occurrence, and the time-to-end of sparks, and vice versa. The characteristics of CREs correlated dose-dependently with the effective coupling strength between RyRs, defined as a function of RyR vicinity, single-channel calcium current, and Mg-binding parameters of the RyR channels. These findings postulate the role of Mg2+ in calcium release as a negative modulator of the coupling strength among RyRs in a CRS, translating to damping of the positive feedback of the calcium-induced calcium-release mechanism.

Endosidin 2 accelerates PIN2 endocytosis and disturbs intracellular trafficking of PIN2, PIN3, and PIN4 but not of SYT1.

We established that Endosidin2 (ES2) affected the trafficking routes of both newly synthesized and endocytic pools of PIN-FORMED2 (PIN2) in Arabidopsis root epidermal cells. PIN2 populations accumulated in separated patches, which gradually merged into large and compact ES2 aggregates (ES2As). FM4-64 endocytic tracer labeled ES2As as well. Both PIN2 pools also appeared in vacuoles. Accelerated endocytosis of PIN2, its aggregation in the cytoplasm, and redirection of PIN2 flows to vacuoles led to a substantial reduction of the abundance of this protein in the plasma membrane. Whereas PIN-FORMED3 and PIN-FORMED4 also aggregated in the cytoplasm, SYT1 was not sensitive to ES2 treatment and did not appear either in the cytoplasmic aggregates or vacuoles. Ultrastructural analysis revealed that ES2 affects the Golgi apparatus so that stacks acquired cup-shape and even circular shape surrounded by several vesicles. Abnormally shaped Golgi stacks, stack remnants, multi-lamellar structures, separated Golgi cisterna rings, tubular structures, and vesicles formed discrete clusters.

The problem of accuracy in single-channel open probability measurements.

The function of ion channels to mediate the flux of ions through membranes of living cells depends on their number, conductance, and open probability. The open probability, PO, characterizes gating of channels that is sensitive to experimental conditions and that can be determined in single-channel experiments. Individual experimental records and even whole series of single-channel activity measurements represent random samples of the stochastic gating continuous in time. The aim of this study was to understand the relationship between the accuracy (trueness and precision) of PO determination and the method of single-channel activity data collection. We used simulated single-channel experiments with variable settings of data collection for a range of open probability values. We found that at low PO, the trueness of PO determination depends on the average number of channel openings per record, while the precision of PO determination depends on the total number of channel openings in the whole dataset and on the distribution of open and closed times. We derived relationships that allow planning of single-channel experiments for the required accuracy of PO determination over a large span of open probabilities.

Structural variability of dyads relates to calcium release in rat ventricular myocytes.

Cardiac excitation-contraction coupling relies on dyads, the intracellular calcium synapses of cardiac myocytes, where the plasma membrane contacts sarcoplasmic reticulum and where electrical excitation triggers calcium release. The morphology of dyads and dynamics of local calcium release vary substantially. To better understand the correspondence between the structure and the functionality of dyads, we estimated incidences of structurally different dyads and of kinetically different calcium release sites and tested their responsiveness to experimental myocardial injury in left ventricular myocytes of rats. According to the structure of dyads estimated in random electron microscopic images of myocardial tissue, the dyads were sorted into 'compact' or 'loose' types. The calcium release fluxes, triggered at local calcium release sites in patch-clamped ventricular myocytes and recorded by laser scanning confocal fluorescence microscopy, were decomposed into 'early' and 'late' components. ANOVA tests revealed very high correlation between the relative amplitudes of early and late calcium release flux components and the relative occurrences of compact and loose dyads in the control and in the injured myocardium. This finding ascertained the relationship between the structure of dyads and the functionality of calcium release sites and the responsiveness of calcium release sites to physical load in cardiac myocytes.

Automatic assessment of the cardiomyocyte development stages from confocal microscopy images using deep convolutional networks.

Computer assisted image acquisition techniques, including confocal microscopy, require efficient tools for an automatic sorting of vast amount of generated image data. The complexity of the classification process, absence of adequate tools, and insufficient amount of reference data has made the automated processing of images challenging. Mastering of this issue would allow implementation of statistical analysis in research areas such as in research on formation of t-tubules in cardiac myocytes. We developed a system aimed at automatic assessment of cardiomyocyte development stages (SAACS). The system classifies confocal images of cardiomyocytes with fluorescent dye stained sarcolemma. We based SAACS on a densely connected convolutional network (DenseNet) topology. We created a set of labelled source images, proposed an appropriate data augmentation technique and designed a class probability graph. We showed that the DenseNet topology, in combination with the augmentation technique is suitable for the given task, and that high-resolution images are instrumental for image categorization. SAACS, in combination with the automatic high-throughput confocal imaging, will allow application of statistical analysis in the research of the tubular system development or remodelling and loss.

Calcium Signaling and Contractility in Cardiac Myocyte of Wolframin Deficient Rats.

Wolframin (Wfs1) is a membrane protein of the sarco/endoplasmic reticulum. Wfs1 mutations are responsible for the Wolfram syndrome, characterized by diabetic and neurological symptoms. Although Wfs1 is expressed in cardiac muscle, its role in this tissue is not clear. We have characterized the effect of invalidation of Wfs1 on calcium signaling-related processes in isolated ventricular myocytes of exon5-Wfs1 deficient rats (Wfs1-e5/-e5) before the onset of overt disease. Calcium transients and contraction were measured in field-stimulated isolated myocytes using confocal microscopy with calcium indicator fluo-3 AM and sarcomere length detection. Calcium currents and their calcium release-dependent inactivation were characterized in whole-cell patch-clamp experiments. At 4 months, Wfs1-e5/-e5 animals were euglycemic, and echocardiographic examination revealed fully compensated cardiac function. In field-stimulated isolated ventricular myocytes, both the amplitude and the duration of contraction of Wfs1-e5/-e5 animals were elevated relative to control Wfs1+/+ littermates. Increased contractility of myocytes resulted largely from prolonged cytosolic calcium transients. Neither the amplitude of calcium currents nor their voltage dependence of activation differed between the two groups. Calcium currents in Wfs1-e5/-e5 myocytes showed a larger extent of inactivation by short voltage prepulses applied to selectively induce calcium release-dependent inactivation of calcium current. Neither the calcium content of the sarcoplasmic reticulum, measured by application of 20 mmol/l caffeine, nor the expression of SERCA2, determined from Western blots, differed significantly in myocytes of Wfs1-e5/-e5 animals compared to control ones. These experiments point to increased duration of calcium release in ventricular myocytes of Wfs1-e5/-e5 animals. We speculate that the lack of functional wolframin might cause changes leading to upregulation of RyR2 channels resulting in prolongation of channel openings and/or a delay in termination of calcium release.

The mannose 6-phosphate/insulin-like growth factor 2 receptor mediates plasminogen-induced efferocytosis.

The plasminogen system is harnessed in a wide variety of physiological processes, such as fibrinolysis, cell migration, or efferocytosis; and accordingly, it is essential upon inflammation, tissue remodeling, wound healing, and for homeostatic maintenance in general. Previously, we identified a plasminogen receptor in the mannose 6-phosphate/insulin-like growth factor 2 receptor (M6P/IGF2R, CD222). Here, we demonstrate by means of genetic knockdown, knockout, and rescue approaches combined with functional studies that M6P/IGF2R is up-regulated on the surface of macrophages, recognizes plasminogen exposed on the surface of apoptotic cells, and mediates plasminogen-induced efferocytosis. The level of uptake of plasminogen-coated apoptotic cells inversely correlates with the TNF-α production by phagocytes indicating tissue clearance without inflammation by this mechanism. Our results reveal an up-to-now undetermined function of M6P/IGF2R in clearance of apoptotic cells, which is crucial for tissue homeostasis.

Alterations in the health of hibernating bats under pathogen pressure.

In underground hibernacula temperate northern hemisphere bats are exposed to Pseudogymnoascus destructans, the fungal agent of white-nose syndrome. While pathological and epidemiological data suggest that Palearctic bats tolerate this infection, we lack knowledge about bat health under pathogen pressure. Here we report blood profiles, along with body mass index (BMI), infection intensity and hibernation temperature, in greater mouse-eared bats (Myotis myotis). We sampled three European hibernacula that differ in geomorphology and microclimatic conditions. Skin lesion counts differed between contralateral wings of a bat, suggesting variable exposure to the fungus. Analysis of blood parameters suggests a threshold of ca. 300 skin lesions on both wings, combined with poor hibernation conditions, may distinguish healthy bats from those with homeostatic disruption. Physiological effects manifested as mild metabolic acidosis, decreased glucose and peripheral blood eosinophilia which were strongly locality-dependent. Hibernating bats displaying blood homeostasis disruption had 2 °C lower body surface temperatures. A shallow BMI loss slope with increasing pathogen load suggested a high degree of infection tolerance. European greater mouse-eared bats generally survive P. destructans invasion, despite some health deterioration at higher infection intensities (dependant on hibernation conditions). Conservation measures should minimise additional stressors to conserve constrained body reserves of bats during hibernation.

Hibernation temperature-dependent Pseudogymnoascus destructans infection intensity in Palearctic bats.

White-nose syndrome (WNS) is a fungal disease caused by Pseudogymnoascus destructans that is devastating to Nearctic bat populations but tolerated by Palearctic bats. Temperature is a factor known to be important for fungal growth and bat choice of hibernation. Here we investigated the effect of temperature on the pathogenic fungal growth in the wild across the Palearctic. We modelled body surface temperature of bats with respect to fungal infection intensity and disease severity and were able to relate this to the mean annual surface temperature at the site. Bats that hibernated at lower temperatures had less fungal growth and fewer skin lesions on their wings. Contrary to expectation derived from laboratory P. destructans culture experiments, natural infection intensity peaked between 5 and 6°C and decreased at warmer hibernating temperature. We made predictive maps based on bat species distributions, temperature and infection intensity and disease severity data to determine not only where P. destructans will be found but also where the infection will be invasive to bats across the Palearctic. Together these data highlight the mechanistic model of the interplay between environmental and biological factors, which determine progression in a wildlife disease.

Calcium release-dependent inactivation precedes formation of the tubular system in developing rat cardiac myocytes.

Developing cardiac myocytes undergo substantial structural and functional changes transforming the mechanism of excitation-contraction coupling from the embryonic form, based on calcium influx through sarcolemmal DHPR calcium channels, to the adult form, relying on local calcium release through RYR calcium channels of sarcoplasmic reticulum stimulated by calcium influx. We characterized day-by-day the postnatal development of the structure of sarcolemma, using techniques of confocal fluorescence microscopy, and the development of the calcium current, measured by the whole-cell patch-clamp in isolated rat ventricular myocytes. We characterized the appearance and expansion of the t-tubule system and compared it with the appearance and progress of the calcium current inactivation induced by the release of calcium ions from sarcoplasmic reticulum as structural and functional measures of direct DHPR-RYR interaction. The release-dependent inactivation of calcium current preceded the development of the t-tubular system by several days, indicating formation of the first DHPR-RYR couplons at the surface sarcolemma and their later spreading close to contractile myofibrils with the growing t-tubules. Large variability of both of the measured parameters among individual myocytes indicates uneven maturation of myocytes within the growing myocardium.

The N-Terminal Region of the Ryanodine Receptor Affects Channel Activation.

Mutations in the cardiac ryanodine receptor (RyR2), the ion channel responsible for release of calcium ions from intracellular stores into cytoplasm, are the cause of several inherited cardiac arrhythmias. At the molecular level, disease symptoms can be mimicked by domain peptides from mutation-prone regions of RyR2 that bind to RyR2 and activate it. Here we show that the domain peptide DPcpvtN2, corresponding to the central helix of the N-terminal region of RyR2, activates the RyR2 channel. Structural modeling of interaction between DPcpvtN2 and the N-terminal region of RyR2 in the closed and open conformation provided three plausible structures of the complex. Only one of them could explain the dependence of RyR2 activity on concentration of DPcpvtN2. The structure of the complex was at odds with the previously proposed "domain switch" mechanism of competition between domain peptides and ryanodine receptor domains. Likewise, in structural models of the N-terminal region, the conformational changes induced by DPcpvtN2 binding were different from those induced by mutation of central helix amino acids. The activating effect of DPcpvtN2 binding and of mutations in the central helix could be explained by their similar effect on the transition energy between the closed and open conformation of RyR2.

RyR2R420Q catecholaminergic polymorphic ventricular tachycardia mutation induces bradycardia by disturbing the coupled clock pacemaker mechanism.

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a lethal genetic arrhythmia that manifests syncope or sudden death in children and young adults under stress conditions. CPVT patients often present bradycardia and sino-atrial node (SAN) dysfunction. However, the mechanism remains unclear. We analyzed SAN function in two CPVT families and in a novel knock-in (KI) mouse model carrying the RyR2R420Q mutation. Humans and KI mice presented slower resting heart rate. Accordingly, the rate of spontaneous intracellular Ca2+ ([Ca2+]i) transients was slower in KI mouse SAN preparations than in WT, without any significant alteration in the "funny" current (If ). The L-type Ca2+ current was reduced in KI SAN cells in a [Ca2+]i-dependent way, suggesting that bradycardia was due to disrupted crosstalk between the "voltage" and "Ca2+" clock, and the mechanisms of pacemaking was induced by aberrant spontaneous RyR2- dependent Ca2+ release. This finding was consistent with a higher Ca2+ leak during diastolic periods produced by long-lasting Ca2+ sparks in KI SAN cells. Our results uncover a mechanism for the CPVT-causing RyR2 N-terminal mutation R420Q, and they highlight the fact that enhancing the Ca2+ clock may slow the heart rhythm by disturbing the coupling between Ca2+ and voltage clocks.

Ectoparasites may serve as vectors for the white-nose syndrome fungus.

BACKGROUND: Vertebrate ectoparasites frequently play a role in transmission of infectious agents. Pseudogymnoascus destructans is a psychrophilic fungus known to cause white-nose syndrome (WNS), an emerging infectious disease of bats. It is transmitted with direct contact between bats or with contaminated environment. The aim of this study was to examine wing mites from the family Spinturnicidae parasitizing hibernating bats for the presence of P. destructans propagules as another possible transmission route. METHODS: Wing mites collected from 33 bats at four hibernation sites in the Czech Republic were inspected for the presence and load of pathogen's DNA using quantitative PCR. Simultaneously, wing damage of inspected bats caused by WNS was quantified using ultraviolet light (UV) transillumination and the relationship between fungal load on wing mites and intensity of infection was subjected to correlation analysis. RESULTS: All samples of wing mites were positive for the presence of DNA of P. destructans, indicating a high probability of their role in the transmission of the pathogen's propagules between bats. CONCLUSIONS: Mechanical transport of adhesive P. destructans spores and mycelium fragments on the body of spinturnicid mites is highly feasible. The specialised lifestyle of mites, i.e., living on bat wing membranes, the sites most typically affected by fungal growth, enables pathogen transport. Moreover, P. destructans metabolic traits suggest an ability to grow and sporulate on a range of organic substrates, including insects, which supports the possibility of growth on bat ectoparasites, at least in periods when bats roost in cold environments and enter torpor. In addition to transport of fungal propagules, mites may facilitate entry of fungal hyphae into the epidermis through injuries caused by biting.

White-nose syndrome without borders: Pseudogymnoascus destructans infection tolerated in Europe and Palearctic Asia but not in North America.

A striking feature of white-nose syndrome, a fungal infection of hibernating bats, is the difference in infection outcome between North America and Europe. Here we show high WNS prevalence both in Europe and on the West Siberian Plain in Asia. Palearctic bat communities tolerate similar fungal loads of Pseudogymnoascus destructans infection as their Nearctic counterparts and histopathology indicates equal focal skin tissue invasiveness pathognomonic for WNS lesions. Fungal load positively correlates with disease intensity and it reaches highest values at intermediate latitudes. Prevalence and fungal load dynamics in Palearctic bats remained persistent and high between 2012 and 2014. Dominant haplotypes of five genes are widespread in North America, Europe and Asia, expanding the source region of white-nose syndrome to non-European hibernacula. Our data provides evidence for both endemicity and tolerance to this persistent virulent fungus in the Palearctic, suggesting that host-pathogen interaction equilibrium has been established.

Ryanodine receptor gating controls generation of diastolic calcium waves in cardiac myocytes.

The role of cardiac ryanodine receptor (RyR) gating in the initiation and propagation of calcium waves was investigated using a mathematical model comprising a stochastic description of RyR gating and a deterministic description of calcium diffusion and sequestration. We used a one-dimensional array of equidistantly spaced RyR clusters, representing the confocal scanning line, to simulate the formation of calcium sparks. Our model provided an excellent description of the calcium dependence of the frequency of diastolic calcium sparks and of the increased tendency for the production of calcium waves after a decrease in cytosolic calcium buffering. We developed a hypothesis relating changes in the propensity to form calcium waves to changes of RyR gating and tested it by simulation. With a realistic RyR gating model, increased ability of RyR to be activated by Ca2+ strongly increased the propensity for generation of calcium waves at low (0.05-0.1-µM) calcium concentrations but only slightly at high (0.2-0.4-µM) calcium concentrations. Changes in RyR gating altered calcium wave formation by changing the calcium sensitivity of spontaneous calcium spark activation and/or the average number of open RyRs in spontaneous calcium sparks. Gating changes that did not affect RyR activation by Ca2+ had only a weak effect on the propensity to form calcium waves, even if they strongly increased calcium spark frequency. Calcium waves induced by modulating the properties of the RyR activation site could be suppressed by inhibiting the spontaneous opening of the RyR. These data can explain the increased tendency for production of calcium waves under conditions when RyR gating is altered in cardiac diseases.

Structural insights into the human RyR2 N-terminal region involved in cardiac arrhythmias.

Human ryanodine receptor 2 (hRyR2) mediates calcium release from the sarcoplasmic reticulum, enabling cardiomyocyte contraction. The N-terminal region of hRyR2 (amino acids 1-606) is the target of >30 arrhythmogenic mutations and contains a binding site for phosphoprotein phosphatase 1. Here, the solution and crystal structures determined under near-physiological conditions, as well as a homology model of the hRyR2 N-terminal region, are presented. The N-terminus is held together by a unique network of interactions among its three domains, A, B and C, in which the central helix (amino acids 410-437) plays a prominent stabilizing role. Importantly, the anion-binding site reported for the mouse RyR2 N-terminal region is notably absent from the human RyR2. The structure concurs with the differential stability of arrhythmogenic mutations in the central helix (R420W, I419F and I419F/R420W) which are owing to disparities in the propensity of mutated residues to form energetically favourable or unfavourable contacts. In solution, the N-terminus adopts a globular shape with a prominent tail that is likely to involve residues 545-606, which are unresolved in the crystal structure. Docking the N-terminal domains into cryo-electron microscopy maps of the closed and open RyR1 conformations reveals C(α) atom movements of up to 8 Šupon channel gating, and predicts the location of the leucine-isoleucine zipper segment and the interaction site for spinophilin and phosphoprotein phosphatase 1 on the RyR surface.