<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0" xmlns:media="http://search.yahoo.com/mrss/" xmlns:ynews="http://news.yahoo.com/rss/">
    <channel>
        <title>Nova Reader - Subject</title>
        <link>https://www.novareader.co</link>
        <description>Default RSS Feed</description>
        <language>en-us</language>
        <copyright>Newgen KnowledgeWorks</copyright>
        <item>
            <title><![CDATA[MyD88 oligomer size functions as a physical threshold to trigger IL1R Myddosome signaling]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766073337115-3174892f-db38-4f75-8669-c47cdab3dd7b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202012071</link>
            <description><![CDATA[<p class="para" id="N65540">Deliz-Aguirre et al. examine the molecular dynamics of interleukin 1 receptor signaling. Receptor sensing of IL-1 triggers the oligomerization of MyD88 and formation of the Myddosome signaling complex. The formation of a MyD88 oligomer of a requisite size serves as a threshold to active downstream signaling.</p><p class="para" id="N65539">A recurring feature of innate immune receptor signaling is the self-assembly of signaling proteins into oligomeric complexes. The Myddosome is an oligomeric complex that is required to transmit inflammatory signals from TLR/IL1Rs and consists of MyD88 and IRAK family kinases. However, the molecular basis for how Myddosome proteins self-assemble and regulate intracellular signaling remains poorly understood. Here, we developed a novel assay to analyze the spatiotemporal dynamics of IL1R and Myddosome signaling in live cells. We found that MyD88 oligomerization is inducible and initially reversible. Moreover, the formation of larger, stable oligomers consisting of more than four MyD88s triggers the sequential recruitment of IRAK4 and IRAK1. Notably, genetic knockout of IRAK4 enhanced MyD88 oligomerization, indicating that IRAK4 controls MyD88 oligomer size and growth. MyD88 oligomer size thus functions as a physical threshold to trigger downstream signaling. These results provide a mechanistic basis for how protein oligomerization might function in cell signaling pathways.</p>]]></description>
            <pubDate><![CDATA[2021-05-06T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Changes in ionizing radiation dose rate affect cell cycle progression in adipose derived stem cells]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766066580084-ddd60386-5e6b-41e6-adba-d1619f0f4652/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0250160</link>
            <description><![CDATA[<p class="para" id="N65539">Biomedical use of radiation is utilized in effective diagnostic and treatment tools, yet can introduce risks to healthy tissues. High energy photons used for diagnostic purposes have high penetration depth and can discriminate multiple tissues based on attenuation properties of different materials. Likewise, the ability to deposit energy at various targets within tumors make the use of photons effective treatment for cancer. Radiation focused on a tumor will deposit energy when it interacts with a biological structure (e.g. DNA), which will result in cell kill should repair capacity of the tissue be overwhelmed. Likewise, damage to normal, non-cancerous tissues is a consequence of radiation that can lead to acute or late, chronic toxicity profiles. Adipose derived stem cells (ADSCs) are mesenchymal stem cells that have been proven to have similar characteristics to bone marrow derived stem cells, except that they are much easier to obtain. Within the body, ADSCs act as immunomodulators and assist with the maintenance and repair of tissues. They have been shown to have excellent differentiation capability, making them an extremely viable option for stem cell therapies and regenerative medicine applications. Due to the tissue ADSCs are derived from, they are highly likely to be affected by radiation therapy, especially when treating tumors localized to structures with relatively high ADSC content (eg., breast cancer). For this reason, the purpose behind this research is to better understand how ADSCs are affected by doses of radiation comparable to a single fraction of radiation therapy. We also measured the response of ADSCs to exposure at different dose rates to determine if there is a significant difference in the response of ADSCs to radiation therapy relevant doses of ionizing radiation. Our findings indicate that ADSCs exposed to Cesium (Cs 137)-gamma rays at a moderate dose of 2Gy and either a low dose rate (1.40Gy/min) or a high dose rate (7.31Gy/min) slow proliferation rate, and with cell cycle arrest in some populations. These responses ADSCs were not as marked as previously measured in other stem cell types. In addition, our results indicate that differences in dose rate in the Gy/min range typically utilized in small animal or cell irradiation platforms have a minimal effect on the function of ADSCs. The potential ADSCs have in the space of regenerative medicine makes them an ideal candidate for study with ionizing radiation, as they are one of the main cell types to promote tissue healing.</p>]]></description>
            <pubDate><![CDATA[2021-04-27T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Changes in ionizing radiation dose rate affect cell cycle progression in adipose derived stem cells]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766066580084-ddd60386-5e6b-41e6-adba-d1619f0f4652/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0250160</link>
            <description><![CDATA[<p class="para" id="N65539">Biomedical use of radiation is utilized in effective diagnostic and treatment tools, yet can introduce risks to healthy tissues. High energy photons used for diagnostic purposes have high penetration depth and can discriminate multiple tissues based on attenuation properties of different materials. Likewise, the ability to deposit energy at various targets within tumors make the use of photons effective treatment for cancer. Radiation focused on a tumor will deposit energy when it interacts with a biological structure (e.g. DNA), which will result in cell kill should repair capacity of the tissue be overwhelmed. Likewise, damage to normal, non-cancerous tissues is a consequence of radiation that can lead to acute or late, chronic toxicity profiles. Adipose derived stem cells (ADSCs) are mesenchymal stem cells that have been proven to have similar characteristics to bone marrow derived stem cells, except that they are much easier to obtain. Within the body, ADSCs act as immunomodulators and assist with the maintenance and repair of tissues. They have been shown to have excellent differentiation capability, making them an extremely viable option for stem cell therapies and regenerative medicine applications. Due to the tissue ADSCs are derived from, they are highly likely to be affected by radiation therapy, especially when treating tumors localized to structures with relatively high ADSC content (eg., breast cancer). For this reason, the purpose behind this research is to better understand how ADSCs are affected by doses of radiation comparable to a single fraction of radiation therapy. We also measured the response of ADSCs to exposure at different dose rates to determine if there is a significant difference in the response of ADSCs to radiation therapy relevant doses of ionizing radiation. Our findings indicate that ADSCs exposed to Cesium (Cs 137)-gamma rays at a moderate dose of 2Gy and either a low dose rate (1.40Gy/min) or a high dose rate (7.31Gy/min) slow proliferation rate, and with cell cycle arrest in some populations. These responses ADSCs were not as marked as previously measured in other stem cell types. In addition, our results indicate that differences in dose rate in the Gy/min range typically utilized in small animal or cell irradiation platforms have a minimal effect on the function of ADSCs. The potential ADSCs have in the space of regenerative medicine makes them an ideal candidate for study with ionizing radiation, as they are one of the main cell types to promote tissue healing.</p>]]></description>
            <pubDate><![CDATA[2021-04-27T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Three corrections for overshoot effect improved the dose for step-and-shoot intensity-modulated radiation therapy]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766062404239-e3b04cb2-ad55-400f-a05e-d79662a26ac2/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0250243</link>
            <description><![CDATA[<p class="para" id="N65539">We measured the overshoot effect in a linac and reduced it using block correction, reverse-sequence correction, and index correction. A StarTrack detector was used on a Varian iX. Five segments, 1 × 10 cm<sup>2</sup> in area, were designed; the centers were at −4, −2, 0, 2, and 4 cm on the <i>x</i> axis for measuring the overshoot effect on a 10 × 10 cm<sup>2</sup> collimator setting. Block correction was applied to two segments. The first was on the new first segment at −6 cm, and the other was on the new last segment at 6 cm. Both two new segments were obtained from the 10 × 10 cm<sup>2</sup> collimator setting. The order of segments was reversed in reverse-sequence correction. Reverse-sequence correction averages the dose at every segment after two irradiations. When we used MLC Shaper, index correction reduced the first segment’s index (cumulative radiation occupation) by 60% and increased the last segment’s radiation occupation by 60% in a new MLC.log file. As for relative dose, the first segment had an overdose of 52.4% and the last segment had an underdose of 48.6%, when irradiated at 1 MU at 600 MU/min. The relative doses at the first segment, irradiated at 1 MU, after block correction, reverse-sequence correction, and index correction were applied decreased from 152.5% to 95.1%, 104.8%, and 100.1%, respectively. The relative doses at the last segment, irradiated at 1 MU, after block correction, reverse-sequence correction, and index correction were applied increased from 48.6% to 97.3%, 91.1%, and 95.9%, respectively. The overshoot effect depended on the speed of irradiation. High irradiation speeds resulted in notable overdosing and underdosing at the first and last segments, respectively. The three corrections mitigated the overshoot effect on dose. To save time and effort, the MLC.log file should be edited with a program in the future.</p>]]></description>
            <pubDate><![CDATA[2021-04-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Three corrections for overshoot effect improved the dose for step-and-shoot intensity-modulated radiation therapy]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766062404239-e3b04cb2-ad55-400f-a05e-d79662a26ac2/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0250243</link>
            <description><![CDATA[<p class="para" id="N65539">We measured the overshoot effect in a linac and reduced it using block correction, reverse-sequence correction, and index correction. A StarTrack detector was used on a Varian iX. Five segments, 1 × 10 cm<sup>2</sup> in area, were designed; the centers were at −4, −2, 0, 2, and 4 cm on the <i>x</i> axis for measuring the overshoot effect on a 10 × 10 cm<sup>2</sup> collimator setting. Block correction was applied to two segments. The first was on the new first segment at −6 cm, and the other was on the new last segment at 6 cm. Both two new segments were obtained from the 10 × 10 cm<sup>2</sup> collimator setting. The order of segments was reversed in reverse-sequence correction. Reverse-sequence correction averages the dose at every segment after two irradiations. When we used MLC Shaper, index correction reduced the first segment’s index (cumulative radiation occupation) by 60% and increased the last segment’s radiation occupation by 60% in a new MLC.log file. As for relative dose, the first segment had an overdose of 52.4% and the last segment had an underdose of 48.6%, when irradiated at 1 MU at 600 MU/min. The relative doses at the first segment, irradiated at 1 MU, after block correction, reverse-sequence correction, and index correction were applied decreased from 152.5% to 95.1%, 104.8%, and 100.1%, respectively. The relative doses at the last segment, irradiated at 1 MU, after block correction, reverse-sequence correction, and index correction were applied increased from 48.6% to 97.3%, 91.1%, and 95.9%, respectively. The overshoot effect depended on the speed of irradiation. High irradiation speeds resulted in notable overdosing and underdosing at the first and last segments, respectively. The three corrections mitigated the overshoot effect on dose. To save time and effort, the MLC.log file should be edited with a program in the future.</p>]]></description>
            <pubDate><![CDATA[2021-04-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Aurora B switches relative strength of kinetochore–microtubule attachment modes for error correction]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766061051908-b8954566-b95a-4ed0-abfb-5f6bc6d9c8cd/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202011117</link>
            <description><![CDATA[<p class="para" id="N65540">Doodhi et al. reconstitute the yeast kinetochore–microtubule interface in vitro and compare the relative strength of distinct kinetochore–microtubule interaction modes. Their results suggest how Aurora B kinase promotes the exchange of kinetochore–microtubule interactions to eliminate aberrant interactions and establish chromosome biorientation.</p><p class="para" id="N65539">To establish chromosome biorientation, aberrant kinetochore–microtubule interaction must be resolved (error correction) by Aurora B kinase. Aurora B differentially regulates kinetochore attachment to the microtubule plus end and its lateral side (end-on and lateral attachment, respectively). However, it is still unclear how kinetochore–microtubule interactions are exchanged during error correction. Here, we reconstituted the budding yeast kinetochore–microtubule interface in vitro by attaching the Ndc80 complexes to nanobeads. These Ndc80C nanobeads recapitulated in vitro the lateral and end-on attachments of authentic kinetochores on dynamic microtubules loaded with the Dam1 complex. This in vitro assay enabled the direct comparison of lateral and end-on attachment strength and showed that Dam1 phosphorylation by Aurora B makes the end-on attachment weaker than the lateral attachment. Similar reconstitutions with purified kinetochore particles were used for comparison. We suggest the Dam1 phosphorylation weakens interaction with the Ndc80 complex, disrupts the end-on attachment, and promotes the exchange to a new lateral attachment, leading to error correction.</p>]]></description>
            <pubDate><![CDATA[2021-04-14T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[<i>De novo</i> assembly of the freshwater prawn <i>Macrobrachium carcinus</i> brain transcriptome for identification of potential targets for antibody development]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766046003708-30b922a8-d6fa-4634-9680-c606fac5444c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0249801</link>
            <description><![CDATA[<p class="para" id="N65539">Crustaceans are major constituents of aquatic ecosystems and, as such, changes in their behavior and the structure and function of their bodies can serve as indicators of alterations in their immediate environment, such as those associated with climate change and anthropogenic contamination. We have used bioinformatics and a <i>de novo</i> transcriptome assembly approach to identify potential targets for developing specific antibodies to serve as nervous system function markers for freshwater prawns of the <i>Macrobrachium</i> spp. Total RNA was extracted from brain ganglia of <i>Macrobrachium carcinus</i> freshwater prawns and Illumina Next Generation Sequencing was performed using an Eel Pond mRNA Seq Protocol to construct a <i>de novo</i> transcriptome. Sequencing yielded 97,202,662 sequences: 47,630,546 paired and 1,941,570 singletons. Assembly with Trinity resulted in 197,898 assembled contigs from which 30,576 were annotated: 9,600 by orthology, 17,197 by homology, and 3,779 by transcript families. We looked for glutamate receptors contigs, due to their main role in crustacean excitatory neurotransmission, and found 138 contigs related to ionotropic receptors, 32 related to metabotropic receptors, and 18 to unidentified receptors. After performing multiple sequence alignments within different biological organisms and antigenicity analysis, we were able to develop antibodies for prawn AMPA ionotropic glutamate receptor 1, metabotropic glutamate receptor 1 and 4, and ionotropic NMDA glutamate receptor subunit 2B, with the expectation that the availability of these antibodies will help broaden knowledge regarding the underlying structural and functional mechanisms involved in prawn behavioral responses to environmental impacts. The <i>Macrobrachium carcinus</i> brain transcriptome can be an important tool for examining changes in many other nervous system molecules as a function of developmental stages, or in response to particular conditions or treatments.</p>]]></description>
            <pubDate><![CDATA[2021-04-09T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[<i>De novo</i> assembly of the freshwater prawn <i>Macrobrachium carcinus</i> brain transcriptome for identification of potential targets for antibody development]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766046003708-30b922a8-d6fa-4634-9680-c606fac5444c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0249801</link>
            <description><![CDATA[<p class="para" id="N65539">Crustaceans are major constituents of aquatic ecosystems and, as such, changes in their behavior and the structure and function of their bodies can serve as indicators of alterations in their immediate environment, such as those associated with climate change and anthropogenic contamination. We have used bioinformatics and a <i>de novo</i> transcriptome assembly approach to identify potential targets for developing specific antibodies to serve as nervous system function markers for freshwater prawns of the <i>Macrobrachium</i> spp. Total RNA was extracted from brain ganglia of <i>Macrobrachium carcinus</i> freshwater prawns and Illumina Next Generation Sequencing was performed using an Eel Pond mRNA Seq Protocol to construct a <i>de novo</i> transcriptome. Sequencing yielded 97,202,662 sequences: 47,630,546 paired and 1,941,570 singletons. Assembly with Trinity resulted in 197,898 assembled contigs from which 30,576 were annotated: 9,600 by orthology, 17,197 by homology, and 3,779 by transcript families. We looked for glutamate receptors contigs, due to their main role in crustacean excitatory neurotransmission, and found 138 contigs related to ionotropic receptors, 32 related to metabotropic receptors, and 18 to unidentified receptors. After performing multiple sequence alignments within different biological organisms and antigenicity analysis, we were able to develop antibodies for prawn AMPA ionotropic glutamate receptor 1, metabotropic glutamate receptor 1 and 4, and ionotropic NMDA glutamate receptor subunit 2B, with the expectation that the availability of these antibodies will help broaden knowledge regarding the underlying structural and functional mechanisms involved in prawn behavioral responses to environmental impacts. The <i>Macrobrachium carcinus</i> brain transcriptome can be an important tool for examining changes in many other nervous system molecules as a function of developmental stages, or in response to particular conditions or treatments.</p>]]></description>
            <pubDate><![CDATA[2021-04-09T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[QuASeR: Quantum Accelerated de novo DNA sequence reconstruction]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766030925925-c748756c-5fa4-4c3c-83f3-03998b50e99a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0249850</link>
            <description><![CDATA[<p class="para" id="N65539">In this article, we present QuASeR, a reference-free DNA sequence reconstruction implementation via de novo assembly on both gate-based and quantum annealing platforms. This is the first time this important application in bioinformatics is modeled using quantum computation. Each one of the four steps of the implementation (TSP, QUBO, Hamiltonians and QAOA) is explained with a proof-of-concept example to target both the genomics research community and quantum application developers in a self-contained manner. The implementation and results on executing the algorithm from a set of DNA reads to a reconstructed sequence, on a gate-based quantum simulator, the D-Wave quantum annealing simulator and hardware are detailed. We also highlight the limitations of current classical simulation and available quantum hardware systems. The implementation is open-source and can be found on https://github.com/QE-Lab/QuASeR.</p>]]></description>
            <pubDate><![CDATA[2021-04-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[QuASeR: Quantum Accelerated de novo DNA sequence reconstruction]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766030925925-c748756c-5fa4-4c3c-83f3-03998b50e99a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0249850</link>
            <description><![CDATA[<p class="para" id="N65539">In this article, we present QuASeR, a reference-free DNA sequence reconstruction implementation via de novo assembly on both gate-based and quantum annealing platforms. This is the first time this important application in bioinformatics is modeled using quantum computation. Each one of the four steps of the implementation (TSP, QUBO, Hamiltonians and QAOA) is explained with a proof-of-concept example to target both the genomics research community and quantum application developers in a self-contained manner. The implementation and results on executing the algorithm from a set of DNA reads to a reconstructed sequence, on a gate-based quantum simulator, the D-Wave quantum annealing simulator and hardware are detailed. We also highlight the limitations of current classical simulation and available quantum hardware systems. The implementation is open-source and can be found on https://github.com/QE-Lab/QuASeR.</p>]]></description>
            <pubDate><![CDATA[2021-04-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[New methods for quantifying rapidity of action potential onset differentiate neuron types]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766023575827-f92d8809-9c91-4dd3-a901-ef511f18c0a0/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0247242</link>
            <description><![CDATA[<p class="para" id="N65539">Two new methods for quantifying the rapidity of action potential onset have lower relative standard deviations and better distinguish neuron cell types than current methods. Action potentials (APs) in most central mammalian neurons exhibit sharp onset dynamics. The main views explaining such an abrupt onset differ. Some studies suggest sharp onsets reflect cooperative sodium channels activation, while others suggest they reflect AP backpropagation from the axon initial segment. However, AP onset rapidity is defined subjectively in these studies, often using the slope at an arbitrary value on the phase plot. Thus, we proposed more systematic methods using the membrane potential’s second-time derivative (V¨m<div class="imageVideo"><img src="" alt=""/></div>) peak width. Here, the AP rapidity was measured for four different cortical and hippocampal neuron types using four quantification methods: the inverse of full-width at the half maximum of the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak (IFWd<sup>2</sup>), the inverse of half-width at the half maximum of the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak (IHWd<sup>2</sup>), the phase plot slope, and the error ratio method. The IFWd<sup>2</sup> and IHWd<sup>2</sup> methods show the smallest variation among neurons of the same type. Furthermore, the AP rapidity, using the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak width methods, significantly differentiates between different types of neurons, indicating that AP rapidity can be used to classify neuron types. The AP rapidity measured using the IFWd<sup>2</sup> method was able to differentiate between all four neuron types analyzed. Therefore, the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak width methods provide another sensitive tool to investigate the mechanisms impacting the AP onset dynamics.</p>]]></description>
            <pubDate><![CDATA[2021-04-08T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[New methods for quantifying rapidity of action potential onset differentiate neuron types]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1766023575827-f92d8809-9c91-4dd3-a901-ef511f18c0a0/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0247242</link>
            <description><![CDATA[<p class="para" id="N65539">Two new methods for quantifying the rapidity of action potential onset have lower relative standard deviations and better distinguish neuron cell types than current methods. Action potentials (APs) in most central mammalian neurons exhibit sharp onset dynamics. The main views explaining such an abrupt onset differ. Some studies suggest sharp onsets reflect cooperative sodium channels activation, while others suggest they reflect AP backpropagation from the axon initial segment. However, AP onset rapidity is defined subjectively in these studies, often using the slope at an arbitrary value on the phase plot. Thus, we proposed more systematic methods using the membrane potential’s second-time derivative (V¨m<div class="imageVideo"><img src="" alt=""/></div>) peak width. Here, the AP rapidity was measured for four different cortical and hippocampal neuron types using four quantification methods: the inverse of full-width at the half maximum of the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak (IFWd<sup>2</sup>), the inverse of half-width at the half maximum of the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak (IHWd<sup>2</sup>), the phase plot slope, and the error ratio method. The IFWd<sup>2</sup> and IHWd<sup>2</sup> methods show the smallest variation among neurons of the same type. Furthermore, the AP rapidity, using the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak width methods, significantly differentiates between different types of neurons, indicating that AP rapidity can be used to classify neuron types. The AP rapidity measured using the IFWd<sup>2</sup> method was able to differentiate between all four neuron types analyzed. Therefore, the V¨m<div class="imageVideo"><img src="" alt=""/></div> peak width methods provide another sensitive tool to investigate the mechanisms impacting the AP onset dynamics.</p>]]></description>
            <pubDate><![CDATA[2021-04-08T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Transfection methods for high-throughput cellular assays of voltage-gated calcium and sodium channels involved in pain]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765992999812-a61c6504-beae-4297-b649-310b2017053d/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0243645</link>
            <description><![CDATA[<p class="para" id="N65539">Chemical transfection is broadly used to transiently transfect mammalian cells, although often associated with cellular stress and membrane instability, which imposes challenges for most cellular assays, including high-throughput (HT) assays. In the current study, we compared the effectiveness of calcium phosphate, FuGENE and Lipofectamine 3000 to transiently express two key voltage-gated ion channels critical in pain pathways, Ca<sub>V</sub>2.2 and Na<sub>V</sub>1.7. The expression and function of these channels were validated using two HT platforms, the Fluorescence Imaging Plate Reader FLIPR<sup>Tetra</sup> and the automated patch clamp QPatch 16X. We found that all transfection methods tested demonstrated similar effectiveness when applied to FLIPR<sup>Tetra</sup> assays. Lipofectamine 3000-mediated transfection produced the largest peak currents for automated patch clamp QPatch assays. However, the FuGENE-mediated transfection was the most effective for QPatch assays as indicated by the superior number of cells displaying GΩ seal formation in whole-cell patch clamp configuration, medium to large peak currents, and higher rates of accomplished assays for both Ca<sub>V</sub>2.2 and Na<sub>V</sub>1.7 channels. Our findings can facilitate the development of HT automated patch clamp assays for the discovery and characterization of novel analgesics and modulators of pain pathways, as well as assisting studies examining the pharmacology of mutated channels.</p>]]></description>
            <pubDate><![CDATA[2021-03-05T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Transfection methods for high-throughput cellular assays of voltage-gated calcium and sodium channels involved in pain]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765992999812-a61c6504-beae-4297-b649-310b2017053d/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0243645</link>
            <description><![CDATA[<p class="para" id="N65539">Chemical transfection is broadly used to transiently transfect mammalian cells, although often associated with cellular stress and membrane instability, which imposes challenges for most cellular assays, including high-throughput (HT) assays. In the current study, we compared the effectiveness of calcium phosphate, FuGENE and Lipofectamine 3000 to transiently express two key voltage-gated ion channels critical in pain pathways, Ca<sub>V</sub>2.2 and Na<sub>V</sub>1.7. The expression and function of these channels were validated using two HT platforms, the Fluorescence Imaging Plate Reader FLIPR<sup>Tetra</sup> and the automated patch clamp QPatch 16X. We found that all transfection methods tested demonstrated similar effectiveness when applied to FLIPR<sup>Tetra</sup> assays. Lipofectamine 3000-mediated transfection produced the largest peak currents for automated patch clamp QPatch assays. However, the FuGENE-mediated transfection was the most effective for QPatch assays as indicated by the superior number of cells displaying GΩ seal formation in whole-cell patch clamp configuration, medium to large peak currents, and higher rates of accomplished assays for both Ca<sub>V</sub>2.2 and Na<sub>V</sub>1.7 channels. Our findings can facilitate the development of HT automated patch clamp assays for the discovery and characterization of novel analgesics and modulators of pain pathways, as well as assisting studies examining the pharmacology of mutated channels.</p>]]></description>
            <pubDate><![CDATA[2021-03-05T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Nonlinear effects of intrinsic dynamics on temporal encoding in a model of avian auditory cortex]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765991295309-e948191a-f5b9-41b8-b651-5f4482e37115/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008768</link>
            <description><![CDATA[<p class="para" id="N65539">Neurons exhibit diverse intrinsic dynamics, which govern how they integrate synaptic inputs to produce spikes. Intrinsic dynamics are often plastic during development and learning, but the effects of these changes on stimulus encoding properties are not well known. To examine this relationship, we simulated auditory responses to zebra finch song using a linear-dynamical cascade model, which combines a linear spectrotemporal receptive field with a dynamical, conductance-based neuron model, then used generalized linear models to estimate encoding properties from the resulting spike trains. We focused on the effects of a low-threshold potassium current (K<sub>LT</sub>) that is present in a subset of cells in the zebra finch caudal mesopallium and is affected by early auditory experience. We found that K<sub>LT</sub> affects both spike adaptation and the temporal filtering properties of the receptive field. The direction of the effects depended on the temporal modulation tuning of the linear (input) stage of the cascade model, indicating a strongly nonlinear relationship. These results suggest that small changes in intrinsic dynamics in tandem with differences in synaptic connectivity can have dramatic effects on the tuning of auditory neurons.</p><p class="para" id="N65542">Experience-dependent developmental plasticity involves changes not only to synaptic connections, but to voltage-gated currents as well. Using biophysical models, it is straightforward to predict the effects of this intrinsic plasticity on the firing patterns of individual neurons, but it remains difficult to understand the consequences for sensory coding. We investigated this in the context of the zebra finch auditory cortex, where early exposure to a complex acoustic environment causes increased expression of a low-threshold potassium current. We simulated responses to song using a detailed biophysical model and then characterized encoding properties using generalized linear models. This analysis revealed that this potassium current has strong, nonlinear effects on how the model encodes the song’s temporal structure, and that the sign of these effects depend on the temporal tuning of the synaptic inputs. This nonlinearity gives intrinsic plasticity broad scope as a mechanism for developmental learning in the auditory system.</p>]]></description>
            <pubDate><![CDATA[2021-02-22T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Nonlinear effects of intrinsic dynamics on temporal encoding in a model of avian auditory cortex]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765991295309-e948191a-f5b9-41b8-b651-5f4482e37115/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008768</link>
            <description><![CDATA[<p class="para" id="N65539">Neurons exhibit diverse intrinsic dynamics, which govern how they integrate synaptic inputs to produce spikes. Intrinsic dynamics are often plastic during development and learning, but the effects of these changes on stimulus encoding properties are not well known. To examine this relationship, we simulated auditory responses to zebra finch song using a linear-dynamical cascade model, which combines a linear spectrotemporal receptive field with a dynamical, conductance-based neuron model, then used generalized linear models to estimate encoding properties from the resulting spike trains. We focused on the effects of a low-threshold potassium current (K<sub>LT</sub>) that is present in a subset of cells in the zebra finch caudal mesopallium and is affected by early auditory experience. We found that K<sub>LT</sub> affects both spike adaptation and the temporal filtering properties of the receptive field. The direction of the effects depended on the temporal modulation tuning of the linear (input) stage of the cascade model, indicating a strongly nonlinear relationship. These results suggest that small changes in intrinsic dynamics in tandem with differences in synaptic connectivity can have dramatic effects on the tuning of auditory neurons.</p><p class="para" id="N65542">Experience-dependent developmental plasticity involves changes not only to synaptic connections, but to voltage-gated currents as well. Using biophysical models, it is straightforward to predict the effects of this intrinsic plasticity on the firing patterns of individual neurons, but it remains difficult to understand the consequences for sensory coding. We investigated this in the context of the zebra finch auditory cortex, where early exposure to a complex acoustic environment causes increased expression of a low-threshold potassium current. We simulated responses to song using a detailed biophysical model and then characterized encoding properties using generalized linear models. This analysis revealed that this potassium current has strong, nonlinear effects on how the model encodes the song’s temporal structure, and that the sign of these effects depend on the temporal tuning of the synaptic inputs. This nonlinearity gives intrinsic plasticity broad scope as a mechanism for developmental learning in the auditory system.</p>]]></description>
            <pubDate><![CDATA[2021-02-22T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Up-down biphasic volume response of human red blood cells to PIEZO1 activation during capillary transits]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765981824364-19e7f0ef-a369-4ef2-81f0-294ac70b1dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008706</link>
            <description><![CDATA[<p class="para" id="N65539">In this paper we apply a novel JAVA version of a model on the homeostasis of human red blood cells (RBCs) to investigate the changes RBCs experience during single capillary transits. In the companion paper we apply a model extension to investigate the changes in RBC homeostasis over the approximately 200000 capillary transits during the ~120 days lifespan of the cells. These are topics inaccessible to direct experimentation but rendered mature for a computational modelling approach by the large body of recent and early experimental results which robustly constrain the range of parameter values and model outcomes, offering a unique opportunity for an in depth study of the mechanisms involved. Capillary transit times vary between 0.5 and 1.5s during which the red blood cells squeeze and deform in the capillary stream transiently opening stress-gated PIEZO1 channels allowing ion gradient dissipation and creating minuscule quantal changes in RBC ion contents and volume. Widely accepted views, based on the effects of experimental shear stress on human RBCs, suggested that quantal changes generated during capillary transits add up over time to develop the documented changes in RBC density and composition during their long circulatory lifespan, the quantal hypothesis. Applying the new red cell model (RCM) we investigated here the changes in homeostatic variables that may be expected during single capillary transits resulting from transient PIEZO1 channel activation. The predicted quantal volume changes were infinitesimal in magnitude, biphasic in nature, and essentially irreversible within inter-transit periods. A sub-second transient PIEZO1 activation triggered a sharp swelling peak followed by a much slower recovery period towards lower-than-baseline volumes. The peak response was caused by net CaCl<sub>2</sub> and fluid gain via PIEZO1 channels driven by the steep electrochemical inward Ca<sup>2+</sup> gradient. The ensuing dehydration followed a complex time-course with sequential, but partially overlapping contributions by KCl loss via Ca<sup>2+</sup>-activated Gardos channels, restorative Ca<sup>2+</sup> extrusion by the plasma membrane calcium pump, and chloride efflux by the Jacobs-Steward mechanism. The change in relative cell volume predicted for single capillary transits was around 10<sup>−5</sup>, an infinitesimal volume change incompatible with a functional role in capillary flow. The biphasic response predicted by the RCM appears to conform to the quantal hypothesis, but whether its cumulative effects could account for the documented changes in density during RBC senescence required an investigation of the effects of myriad transits over the full four months circulatory lifespan of the cells, the subject of the next paper.</p><p class="para" id="N65542">Each human red blood cell traverses narrow capillaries about 1000 to 2000 times per day. Cell deformability is essential for smooth transits and this requires tight control of cell volume well below spherical maxima. The pump-leak ion-flux balance controlling red cell volume involves sodium-potassium and calcium pumps and a constellation of passive membrane transporters. PIEZO1, a mechanosensitive ion channel plays a central role, activating transiently during capillary transits. The main question, not accessible to direct experimentation, is what happens to red cell volume and composition <i>in vivo</i> during each capillary transit. We attempted answering this question by applying novel extensions to a tried and tested model of red cell homeostasis that encodes the current knowledge on the subject. The results revealed an unexpected up-down biphasic volume response characterized by a sharp initial surge caused by CaCl<sub>2</sub> influx via PIEZO1, associated with osmotic fluid gain, followed by a slow reversal, a prediction in conflict with prevailing views. The infinitesimal amplitude of the predicted volume changes precludes them playing a significant role in flow dynamics. In the following paper we explore how the biphasic response helps explain the age-related changes red cells experience during their long circulatory lifespan.</p>]]></description>
            <pubDate><![CDATA[2021-03-03T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Up-down biphasic volume response of human red blood cells to PIEZO1 activation during capillary transits]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765981824364-19e7f0ef-a369-4ef2-81f0-294ac70b1dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008706</link>
            <description><![CDATA[<p class="para" id="N65539">In this paper we apply a novel JAVA version of a model on the homeostasis of human red blood cells (RBCs) to investigate the changes RBCs experience during single capillary transits. In the companion paper we apply a model extension to investigate the changes in RBC homeostasis over the approximately 200000 capillary transits during the ~120 days lifespan of the cells. These are topics inaccessible to direct experimentation but rendered mature for a computational modelling approach by the large body of recent and early experimental results which robustly constrain the range of parameter values and model outcomes, offering a unique opportunity for an in depth study of the mechanisms involved. Capillary transit times vary between 0.5 and 1.5s during which the red blood cells squeeze and deform in the capillary stream transiently opening stress-gated PIEZO1 channels allowing ion gradient dissipation and creating minuscule quantal changes in RBC ion contents and volume. Widely accepted views, based on the effects of experimental shear stress on human RBCs, suggested that quantal changes generated during capillary transits add up over time to develop the documented changes in RBC density and composition during their long circulatory lifespan, the quantal hypothesis. Applying the new red cell model (RCM) we investigated here the changes in homeostatic variables that may be expected during single capillary transits resulting from transient PIEZO1 channel activation. The predicted quantal volume changes were infinitesimal in magnitude, biphasic in nature, and essentially irreversible within inter-transit periods. A sub-second transient PIEZO1 activation triggered a sharp swelling peak followed by a much slower recovery period towards lower-than-baseline volumes. The peak response was caused by net CaCl<sub>2</sub> and fluid gain via PIEZO1 channels driven by the steep electrochemical inward Ca<sup>2+</sup> gradient. The ensuing dehydration followed a complex time-course with sequential, but partially overlapping contributions by KCl loss via Ca<sup>2+</sup>-activated Gardos channels, restorative Ca<sup>2+</sup> extrusion by the plasma membrane calcium pump, and chloride efflux by the Jacobs-Steward mechanism. The change in relative cell volume predicted for single capillary transits was around 10<sup>−5</sup>, an infinitesimal volume change incompatible with a functional role in capillary flow. The biphasic response predicted by the RCM appears to conform to the quantal hypothesis, but whether its cumulative effects could account for the documented changes in density during RBC senescence required an investigation of the effects of myriad transits over the full four months circulatory lifespan of the cells, the subject of the next paper.</p><p class="para" id="N65542">Each human red blood cell traverses narrow capillaries about 1000 to 2000 times per day. Cell deformability is essential for smooth transits and this requires tight control of cell volume well below spherical maxima. The pump-leak ion-flux balance controlling red cell volume involves sodium-potassium and calcium pumps and a constellation of passive membrane transporters. PIEZO1, a mechanosensitive ion channel plays a central role, activating transiently during capillary transits. The main question, not accessible to direct experimentation, is what happens to red cell volume and composition <i>in vivo</i> during each capillary transit. We attempted answering this question by applying novel extensions to a tried and tested model of red cell homeostasis that encodes the current knowledge on the subject. The results revealed an unexpected up-down biphasic volume response characterized by a sharp initial surge caused by CaCl<sub>2</sub> influx via PIEZO1, associated with osmotic fluid gain, followed by a slow reversal, a prediction in conflict with prevailing views. The infinitesimal amplitude of the predicted volume changes precludes them playing a significant role in flow dynamics. In the following paper we explore how the biphasic response helps explain the age-related changes red cells experience during their long circulatory lifespan.</p>]]></description>
            <pubDate><![CDATA[2021-03-03T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Pathogenic mutations in the kinesin-3 motor KIF1A diminish force generation and movement through allosteric mechanisms]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765978712883-f54306c7-6047-42c8-a026-88e89954ac1e/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202004227</link>
            <description><![CDATA[<p class="para" id="N65540">Kinesin-3 motors are fast and superprocessive, but their force generation properties remain unclear. The authors show that KIF1A detaches under low opposing forces but rapidly reattaches to continue motility. Rapid reattachment depends on the class-specific K-loop, whereas mutations linked to neurodevelopmental disorders impair force generation and motility.</p><p class="para" id="N65539">The kinesin-3 motor KIF1A functions in neurons, where its fast and superprocessive motility facilitates long-distance transport, but little is known about its force-generating properties. Using optical tweezers, we demonstrate that KIF1A stalls at an opposing load of ~3 pN but more frequently detaches at lower forces. KIF1A rapidly reattaches to the microtubule to resume motion due to its class-specific K-loop, resulting in a unique clustering of force generation events. To test the importance of neck linker docking in KIF1A force generation, we introduced mutations linked to human neurodevelopmental disorders. Molecular dynamics simulations predict that V8M and Y89D mutations impair neck linker docking. Indeed, both mutations dramatically reduce the force generation of KIF1A but not the motor’s ability to rapidly reattach to the microtubule. Although both mutations relieve autoinhibition of the full-length motor, the mutant motors display decreased velocities, run lengths, and landing rates and delayed cargo transport in cells. These results advance our understanding of how mutations in KIF1A can manifest in disease.</p>]]></description>
            <pubDate><![CDATA[2021-01-26T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[G1/S transcription factors assemble in increasing numbers of discrete clusters through G1 phase]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765973549167-2f9a8f54-f13f-4c7c-8c54-8ffa7c71ed03/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202003041</link>
            <description><![CDATA[<p class="para" id="N65540">Super-resolution imaging shows that G1/S transcription factors in budding yeast are spatially organized in small clusters that increase in number, not molecular content, as the transcription factor copy number increases with cell growth in G1. This spatial and temporal organization of the G1/S regulon might help coordinate the Start transition.</p><p class="para" id="N65539">In budding yeast, the transcription factors SBF and MBF activate a large program of gene expression in late G1 phase that underlies commitment to cell division, termed Start. SBF/MBF are limiting with respect to target promoters in small G1 phase cells and accumulate as cells grow, raising the questions of how SBF/MBF are dynamically distributed across the G1/S regulon and how this impacts the Start transition. Super-resolution Photo-Activatable Localization Microscopy (PALM) mapping of the static positions of SBF/MBF subunits in fixed cells revealed each transcription factor was organized into discrete clusters containing approximately eight copies regardless of cell size and that the total number of clusters increased as cells grew through G1 phase. Stochastic modeling using reasonable biophysical parameters recapitulated growth-dependent SBF/MBF clustering and predicted TF dynamics that were confirmed in live cell PALM experiments. This spatio-temporal organization of SBF/MBF may help coordinate activation of G1/S regulon and the Start transition.</p>]]></description>
            <pubDate><![CDATA[2020-08-03T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765969412596-6fc3ec41-2b10-4f9d-b7ac-16bc7756e923/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008089</link>
            <description><![CDATA[<p class="para" id="N65539">Short QT (SQT) syndrome is a genetic cardiac disorder characterized by an abbreviated QT interval of the patient’s electrocardiogram. The syndrome is associated with increased risk of arrhythmia and sudden cardiac death and can arise from a number of ion channel mutations. Cardiomyocytes derived from induced pluripotent stem cells generated from SQT patients (SQT hiPSC-CMs) provide promising platforms for testing pharmacological treatments directly in human cardiac cells exhibiting mutations specific for the syndrome. However, a difficulty is posed by the relative immaturity of hiPSC-CMs, with the possibility that drug effects observed in SQT hiPSC-CMs could be very different from the corresponding drug effect <i>in vivo</i>. In this paper, we apply a multistep computational procedure for translating measured drug effects from these cells to human QT response. This process first detects drug effects on individual ion channels based on measurements of SQT hiPSC-CMs and then uses these results to estimate the drug effects on ventricular action potentials and QT intervals of adult SQT patients. We find that the procedure is able to identify IC<sub>50</sub> values in line with measured values for the four drugs quinidine, ivabradine, ajmaline and mexiletine. In addition, the predicted effect of quinidine on the adult QT interval is in good agreement with measured effects of quinidine for adult patients. Consequently, the computational procedure appears to be a useful tool for helping predicting adult drug responses from pure <i>in vitro</i> measurements of patient derived cell lines.</p><p class="para" id="N65542">A number of cardiac disorders originate from genetic mutations affecting the function of ion channels populating the membrane of cardiomyocytes. One example is short QT syndrome, associated with increased risk of arrhythmias and sudden cardiac death. Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) provide a promising platform for testing potential pharmacological treatments for such disorders, as human cardiomyocytes exhibiting specific mutations can be generated and exposed to drugs <i>in vitro</i>. However, the electrophysiological properties of hiPSC-CMs differ significantly from those of adult native cardiomyocytes. Therefore, drug effects observed for hiPSC-CMs could possibly be very different from corresponding drug effects for adult cells <i>in vivo</i>. In this study, we apply a computational framework for translating drug effects observed for hiPSC-CMs derived from a short QT patient to drug effects for adult short QT cardiomyocytes. For one of the considered drugs, the effect on adult QT intervals has been measured and these measurements turn out to be in good agreement with the response estimated by the computational procedure. Thus, the computational framework shows promise for being a useful tool for predicting adult drug responses from measurements of hiPSC-CMs, allowing earlier identification of compounds to accurately treat cardiac diseases.</p>]]></description>
            <pubDate><![CDATA[2021-02-16T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765969412596-6fc3ec41-2b10-4f9d-b7ac-16bc7756e923/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008089</link>
            <description><![CDATA[<p class="para" id="N65539">Short QT (SQT) syndrome is a genetic cardiac disorder characterized by an abbreviated QT interval of the patient’s electrocardiogram. The syndrome is associated with increased risk of arrhythmia and sudden cardiac death and can arise from a number of ion channel mutations. Cardiomyocytes derived from induced pluripotent stem cells generated from SQT patients (SQT hiPSC-CMs) provide promising platforms for testing pharmacological treatments directly in human cardiac cells exhibiting mutations specific for the syndrome. However, a difficulty is posed by the relative immaturity of hiPSC-CMs, with the possibility that drug effects observed in SQT hiPSC-CMs could be very different from the corresponding drug effect <i>in vivo</i>. In this paper, we apply a multistep computational procedure for translating measured drug effects from these cells to human QT response. This process first detects drug effects on individual ion channels based on measurements of SQT hiPSC-CMs and then uses these results to estimate the drug effects on ventricular action potentials and QT intervals of adult SQT patients. We find that the procedure is able to identify IC<sub>50</sub> values in line with measured values for the four drugs quinidine, ivabradine, ajmaline and mexiletine. In addition, the predicted effect of quinidine on the adult QT interval is in good agreement with measured effects of quinidine for adult patients. Consequently, the computational procedure appears to be a useful tool for helping predicting adult drug responses from pure <i>in vitro</i> measurements of patient derived cell lines.</p><p class="para" id="N65542">A number of cardiac disorders originate from genetic mutations affecting the function of ion channels populating the membrane of cardiomyocytes. One example is short QT syndrome, associated with increased risk of arrhythmias and sudden cardiac death. Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) provide a promising platform for testing potential pharmacological treatments for such disorders, as human cardiomyocytes exhibiting specific mutations can be generated and exposed to drugs <i>in vitro</i>. However, the electrophysiological properties of hiPSC-CMs differ significantly from those of adult native cardiomyocytes. Therefore, drug effects observed for hiPSC-CMs could possibly be very different from corresponding drug effects for adult cells <i>in vivo</i>. In this study, we apply a computational framework for translating drug effects observed for hiPSC-CMs derived from a short QT patient to drug effects for adult short QT cardiomyocytes. For one of the considered drugs, the effect on adult QT intervals has been measured and these measurements turn out to be in good agreement with the response estimated by the computational procedure. Thus, the computational framework shows promise for being a useful tool for predicting adult drug responses from measurements of hiPSC-CMs, allowing earlier identification of compounds to accurately treat cardiac diseases.</p>]]></description>
            <pubDate><![CDATA[2021-02-16T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Arrhythmogenic effects of ultra-long and bistable cardiac action potentials]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765969395992-c9fbf4e9-06d6-4693-912e-e774f44db7c9/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008683</link>
            <description><![CDATA[<p class="para" id="N65539">Contemporary accounts of the initiation of cardiac arrhythmias typically rely on after-depolarizations as the trigger for reentrant activity. The after-depolarizations are usually triggered by calcium entry or spontaneous release within the cells of the myocardium or the conduction system. Here we propose an alternative mechanism whereby arrhythmias are triggered autonomously by cardiac cells that fail to repolarize after a normal heartbeat. We investigated the proposal by representing the heart as an excitable medium of FitzHugh-Nagumo cells where a proportion of cells were capable of remaining depolarized indefinitely. As such, those cells exhibit bistable membrane dynamics. We found that heterogeneous media can tolerate a surprisingly large number of bistable cells and still support normal rhythmic activity. Yet there is a critical limit beyond which the medium is persistently arrhythmogenic. Numerical analysis revealed that the critical threshold for arrhythmogenesis depends on both the strength of the coupling between cells and the extent to which the abnormal cells resist repolarization. Moreover, arrhythmogenesis was found to emerge preferentially at tissue boundaries where cells naturally have fewer neighbors to influence their behavior. These findings may explain why atrial fibrillation typically originates from tissue boundaries such as the cuff of the pulmonary vein.</p><p class="para" id="N65542">Cardiac fibrillation is a medical condition where normal heart function is compromised as electrical activity becomes disordered. How fibrillation arises spontaneously is not fully understood. It is generally thought to be triggered by premature depolarization of the cardiac action potential in one or more cells. Those premature beats, known as after-depolarizations, subsequently initiate a self-sustaining rotor in the otherwise normal heart tissue. In this study, we propose an alternative mechanism whereby arrhythmias are initiated by cardiac cells that fail to repolarize of their own accord but still operate normally when embedded in functional heart tissue. We find that such cells can act as focal ectopic sources under appropriate conditions of inter-cellular coupling. Moreover, those cells are more prone to initiating arrhythmia when they are located on natural tissue boundaries. This may explain why atrial fibrillation typically originates from the site where the pulmonary vein attaches to the wall of the heart.</p>]]></description>
            <pubDate><![CDATA[2021-02-16T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Arrhythmogenic effects of ultra-long and bistable cardiac action potentials]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765969395992-c9fbf4e9-06d6-4693-912e-e774f44db7c9/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008683</link>
            <description><![CDATA[<p class="para" id="N65539">Contemporary accounts of the initiation of cardiac arrhythmias typically rely on after-depolarizations as the trigger for reentrant activity. The after-depolarizations are usually triggered by calcium entry or spontaneous release within the cells of the myocardium or the conduction system. Here we propose an alternative mechanism whereby arrhythmias are triggered autonomously by cardiac cells that fail to repolarize after a normal heartbeat. We investigated the proposal by representing the heart as an excitable medium of FitzHugh-Nagumo cells where a proportion of cells were capable of remaining depolarized indefinitely. As such, those cells exhibit bistable membrane dynamics. We found that heterogeneous media can tolerate a surprisingly large number of bistable cells and still support normal rhythmic activity. Yet there is a critical limit beyond which the medium is persistently arrhythmogenic. Numerical analysis revealed that the critical threshold for arrhythmogenesis depends on both the strength of the coupling between cells and the extent to which the abnormal cells resist repolarization. Moreover, arrhythmogenesis was found to emerge preferentially at tissue boundaries where cells naturally have fewer neighbors to influence their behavior. These findings may explain why atrial fibrillation typically originates from tissue boundaries such as the cuff of the pulmonary vein.</p><p class="para" id="N65542">Cardiac fibrillation is a medical condition where normal heart function is compromised as electrical activity becomes disordered. How fibrillation arises spontaneously is not fully understood. It is generally thought to be triggered by premature depolarization of the cardiac action potential in one or more cells. Those premature beats, known as after-depolarizations, subsequently initiate a self-sustaining rotor in the otherwise normal heart tissue. In this study, we propose an alternative mechanism whereby arrhythmias are initiated by cardiac cells that fail to repolarize of their own accord but still operate normally when embedded in functional heart tissue. We find that such cells can act as focal ectopic sources under appropriate conditions of inter-cellular coupling. Moreover, those cells are more prone to initiating arrhythmia when they are located on natural tissue boundaries. This may explain why atrial fibrillation typically originates from the site where the pulmonary vein attaches to the wall of the heart.</p>]]></description>
            <pubDate><![CDATA[2021-02-16T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Ephaptic coupling in white matter fibre bundles modulates axonal transmission delays]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765922619812-ffff8ffa-4b87-4a37-b4fb-90fe4e867e4a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007858</link>
            <description><![CDATA[<p class="para" id="N65539">Axonal connections are widely regarded as faithful transmitters of neuronal signals with fixed delays. The reasoning behind this is that extracellular potentials caused by spikes travelling along axons are too small to have an effect on other axons. Here we devise a computational framework that allows us to study the effect of extracellular potentials generated by spike volleys in axonal fibre bundles on axonal transmission delays. We demonstrate that, although the extracellular potentials generated by single spikes are of the order of microvolts, the collective extracellular potential generated by spike volleys can reach several millivolts. As a consequence, the resulting depolarisation of the axonal membranes increases the velocity of spikes, and therefore reduces axonal delays between brain areas. Driving a neural mass model with such spike volleys, we further demonstrate that only ephaptic coupling can explain the reduction of stimulus latencies with increased stimulus intensities, as observed in many psychological experiments.</p><p class="para" id="N65542">Axonal fibre bundles that connect distant cortical areas contain millions of densely packed axons. When synchronous spike volleys travel through such fibre bundles, the extracellular potential within the bundles is perturbed. We use computer simulations to examine the magnitude and shape of this perturbation, and demonstrate that it is sufficiently strong to affect axonal transmission speeds. Since most spikes within a spike volley are positioned in an area where the extracellular potential is negative (relative to a distant reference), the resulting depolarisation of the axonal membranes accelerates the spike volley on average. This finding is in contrast to previous studies of ephaptic coupling effects between axons, where ephaptic coupling was found to slow down spike propagation. Our finding has consequences for information transmission and synchronisation between cortical areas.</p>]]></description>
            <pubDate><![CDATA[2021-02-08T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Ephaptic coupling in white matter fibre bundles modulates axonal transmission delays]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765922619812-ffff8ffa-4b87-4a37-b4fb-90fe4e867e4a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007858</link>
            <description><![CDATA[<p class="para" id="N65539">Axonal connections are widely regarded as faithful transmitters of neuronal signals with fixed delays. The reasoning behind this is that extracellular potentials caused by spikes travelling along axons are too small to have an effect on other axons. Here we devise a computational framework that allows us to study the effect of extracellular potentials generated by spike volleys in axonal fibre bundles on axonal transmission delays. We demonstrate that, although the extracellular potentials generated by single spikes are of the order of microvolts, the collective extracellular potential generated by spike volleys can reach several millivolts. As a consequence, the resulting depolarisation of the axonal membranes increases the velocity of spikes, and therefore reduces axonal delays between brain areas. Driving a neural mass model with such spike volleys, we further demonstrate that only ephaptic coupling can explain the reduction of stimulus latencies with increased stimulus intensities, as observed in many psychological experiments.</p><p class="para" id="N65542">Axonal fibre bundles that connect distant cortical areas contain millions of densely packed axons. When synchronous spike volleys travel through such fibre bundles, the extracellular potential within the bundles is perturbed. We use computer simulations to examine the magnitude and shape of this perturbation, and demonstrate that it is sufficiently strong to affect axonal transmission speeds. Since most spikes within a spike volley are positioned in an area where the extracellular potential is negative (relative to a distant reference), the resulting depolarisation of the axonal membranes accelerates the spike volley on average. This finding is in contrast to previous studies of ephaptic coupling effects between axons, where ephaptic coupling was found to slow down spike propagation. Our finding has consequences for information transmission and synchronisation between cortical areas.</p>]]></description>
            <pubDate><![CDATA[2021-02-08T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Candidate drugs for preventive treatment of unruptured intracranial aneurysms: A cross-sectional study]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765882314305-afdcc119-c6c2-4060-bfeb-9c5e8420757a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246865</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Background and purpose</h3><p class="para" id="N65543">Establishment of drug therapy to prevent rupture of unruptured intracranial aneurysms (IAs) is needed. Previous human and animal studies have gradually clarified candidate drugs for preventive treatment of IA rupture. However, because most of these candidates belong to classes of drugs frequently co-administered to prevent cardiovascular diseases, epidemiological studies evaluating these drugs simultaneously should be performed. Furthermore, because drugs included in the same class may have different effects in terms of disease prevention, drug-by-drug assessments are important for planning intervention trials.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Materials and methods</h3><p class="para" id="N65549">We performed a cross-sectional study enrolling patients diagnosed with IAs between July 2011 and June 2019 at our institution. Patients were divided into ruptured or unruptured groups. The drugs investigated were selected according to evidence suggested by either human or animal studies. Univariate and multivariate logistic regression analyses were performed to assess the association of drug treatment with rupture status. We also performed drug-by-drug assessments of the association, including dose-response relationships, with rupture status.</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Results</h3><p class="para" id="N65555">In total, 310 patients with ruptured and 887 patients with unruptured IAs were included. Multivariate analysis revealed an inverse association of statins (odds ratio (OR), 0.54; 95% confidence interval (CI) 0.38–0.77), calcium channel blockers (OR, 0.41; 95% CI 0.30–0.58), and angiotensin II receptor blockers (ARBs) (OR, 0.67; 95% CI 0.48–0.93) with ruptured IAs. Moreover, inverse dose-response relationships with rupture status were observed for pitavastatin and rosuvastatin among statins, benidipine, cilnidipine, and amlodipine among calcium channel blockers, and valsartan, azilsartan, candesartan, and olmesartan among ARBs. Only non-aspirin non-steroidal anti-inflammatory drugs were positively associated with ruptured IAs (OR, 3.24; 95% CI 1.71–6.13).</p></div><div class="section" id="sec004"><h3 class="BHead" id="nov000-4">Conclusions</h3><p class="para" id="N65561">The present analysis suggests that several types of statins, calcium channel blockers, and ARBs are candidate drugs for preventive treatment of unruptured IAs.</p></div>]]></description>
            <pubDate><![CDATA[2021-02-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Candidate drugs for preventive treatment of unruptured intracranial aneurysms: A cross-sectional study]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765882314305-afdcc119-c6c2-4060-bfeb-9c5e8420757a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246865</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Background and purpose</h3><p class="para" id="N65543">Establishment of drug therapy to prevent rupture of unruptured intracranial aneurysms (IAs) is needed. Previous human and animal studies have gradually clarified candidate drugs for preventive treatment of IA rupture. However, because most of these candidates belong to classes of drugs frequently co-administered to prevent cardiovascular diseases, epidemiological studies evaluating these drugs simultaneously should be performed. Furthermore, because drugs included in the same class may have different effects in terms of disease prevention, drug-by-drug assessments are important for planning intervention trials.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Materials and methods</h3><p class="para" id="N65549">We performed a cross-sectional study enrolling patients diagnosed with IAs between July 2011 and June 2019 at our institution. Patients were divided into ruptured or unruptured groups. The drugs investigated were selected according to evidence suggested by either human or animal studies. Univariate and multivariate logistic regression analyses were performed to assess the association of drug treatment with rupture status. We also performed drug-by-drug assessments of the association, including dose-response relationships, with rupture status.</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Results</h3><p class="para" id="N65555">In total, 310 patients with ruptured and 887 patients with unruptured IAs were included. Multivariate analysis revealed an inverse association of statins (odds ratio (OR), 0.54; 95% confidence interval (CI) 0.38–0.77), calcium channel blockers (OR, 0.41; 95% CI 0.30–0.58), and angiotensin II receptor blockers (ARBs) (OR, 0.67; 95% CI 0.48–0.93) with ruptured IAs. Moreover, inverse dose-response relationships with rupture status were observed for pitavastatin and rosuvastatin among statins, benidipine, cilnidipine, and amlodipine among calcium channel blockers, and valsartan, azilsartan, candesartan, and olmesartan among ARBs. Only non-aspirin non-steroidal anti-inflammatory drugs were positively associated with ruptured IAs (OR, 3.24; 95% CI 1.71–6.13).</p></div><div class="section" id="sec004"><h3 class="BHead" id="nov000-4">Conclusions</h3><p class="para" id="N65561">The present analysis suggests that several types of statins, calcium channel blockers, and ARBs are candidate drugs for preventive treatment of unruptured IAs.</p></div>]]></description>
            <pubDate><![CDATA[2021-02-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Development of a time-resolved mirrorless scintillation detector]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765882155276-96947cd6-4a3c-48f4-9715-f4f0a75bca11/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246742</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Purpose</h3><p class="para" id="N65543">We developed a compact and lightweight time-resolved mirrorless scintillation detector (TRMLSD) employing image processing techniques and a convolutional neural network (CNN) for high-resolution two-dimensional (2D) dosimetry.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Methods</h3><p class="para" id="N65549">The TRMLSD comprises a camera and an inorganic scintillator plate without a mirror. The camera was installed at a certain angle from the horizontal plane to collect scintillation from the scintillator plate. The geometric distortion due to the absence of a mirror and camera lens was corrected using a projective transform. Variations in brightness due to the distance between the image sensor and each point on the scintillator plate and the inhomogeneity of the material constituting the scintillator were corrected using a 20.0 × 20.0 cm<sup>2</sup> radiation field. Hot pixels were removed using a frame-based noise-reduction technique. Finally, a CNN-based 2D dose distribution deconvolution model was applied to compensate for the dose error in the penumbra region and a lack of backscatter. The linearity, reproducibility, dose rate dependency, and dose profile were tested for a 6 MV X-ray beam to verify dosimeter characteristics. Gamma analysis was performed for two simple and 10 clinical intensity-modulated radiation therapy (IMRT) plans.</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Results</h3><p class="para" id="N65558">The dose linearity with brightness ranging from 0.0 cGy to 200.0 cGy was 0.9998 (R-squared value), and the root-mean-square error value was 1.010. For five consecutive measurements, the reproducibility was within 3% error, and the dose rate dependency was within 1%. The depth dose distribution and lateral dose profile coincided with the ionization chamber data with a 1% mean error. In 2D dosimetry for IMRT plans, the mean gamma passing rates with a 3%/3 mm gamma criterion for the two simple and ten clinical IMRT plans were 96.77% and 95.75%, respectively.</p></div><div class="section" id="sec004"><h3 class="BHead" id="nov000-4">Conclusion</h3><p class="para" id="N65564">The verified accuracy and time-resolved characteristics of the dosimeter may be useful for the quality assurance of machines and patient-specific quality assurance for clinical step-and-shoot IMRT plans.</p></div>]]></description>
            <pubDate><![CDATA[2021-02-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Development of a time-resolved mirrorless scintillation detector]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765882155276-96947cd6-4a3c-48f4-9715-f4f0a75bca11/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246742</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Purpose</h3><p class="para" id="N65543">We developed a compact and lightweight time-resolved mirrorless scintillation detector (TRMLSD) employing image processing techniques and a convolutional neural network (CNN) for high-resolution two-dimensional (2D) dosimetry.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Methods</h3><p class="para" id="N65549">The TRMLSD comprises a camera and an inorganic scintillator plate without a mirror. The camera was installed at a certain angle from the horizontal plane to collect scintillation from the scintillator plate. The geometric distortion due to the absence of a mirror and camera lens was corrected using a projective transform. Variations in brightness due to the distance between the image sensor and each point on the scintillator plate and the inhomogeneity of the material constituting the scintillator were corrected using a 20.0 × 20.0 cm<sup>2</sup> radiation field. Hot pixels were removed using a frame-based noise-reduction technique. Finally, a CNN-based 2D dose distribution deconvolution model was applied to compensate for the dose error in the penumbra region and a lack of backscatter. The linearity, reproducibility, dose rate dependency, and dose profile were tested for a 6 MV X-ray beam to verify dosimeter characteristics. Gamma analysis was performed for two simple and 10 clinical intensity-modulated radiation therapy (IMRT) plans.</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Results</h3><p class="para" id="N65558">The dose linearity with brightness ranging from 0.0 cGy to 200.0 cGy was 0.9998 (R-squared value), and the root-mean-square error value was 1.010. For five consecutive measurements, the reproducibility was within 3% error, and the dose rate dependency was within 1%. The depth dose distribution and lateral dose profile coincided with the ionization chamber data with a 1% mean error. In 2D dosimetry for IMRT plans, the mean gamma passing rates with a 3%/3 mm gamma criterion for the two simple and ten clinical IMRT plans were 96.77% and 95.75%, respectively.</p></div><div class="section" id="sec004"><h3 class="BHead" id="nov000-4">Conclusion</h3><p class="para" id="N65564">The verified accuracy and time-resolved characteristics of the dosimeter may be useful for the quality assurance of machines and patient-specific quality assurance for clinical step-and-shoot IMRT plans.</p></div>]]></description>
            <pubDate><![CDATA[2021-02-12T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[MASSpy: Building, simulating, and visualizing dynamic biological models in Python using mass action kinetics]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765872535383-b8832290-9398-4314-917d-a2593da181f8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008208</link>
            <description><![CDATA[<p class="para" id="N65539">Mathematical models of metabolic networks utilize simulation to study system-level mechanisms and functions. Various approaches have been used to model the steady state behavior of metabolic networks using genome-scale reconstructions, but formulating dynamic models from such reconstructions continues to be a key challenge. Here, we present the Mass Action Stoichiometric Simulation Python (MASSpy) package, an open-source computational framework for dynamic modeling of metabolism. MASSpy utilizes mass action kinetics and detailed chemical mechanisms to build dynamic models of complex biological processes. MASSpy adds dynamic modeling tools to the COnstraint-Based Reconstruction and Analysis Python (COBRApy) package to provide an unified framework for constraint-based and kinetic modeling of metabolic networks. MASSpy supports high-performance dynamic simulation through its implementation of libRoadRunner: the Systems Biology Markup Language (SBML) simulation engine. Three examples are provided to demonstrate how to use MASSpy: (1) a validation of the MASSpy modeling tool through dynamic simulation of detailed mechanisms of enzyme regulation; (2) a feature demonstration using a workflow for generating ensemble of kinetic models using Monte Carlo sampling to approximate missing numerical values of parameters and to quantify biological uncertainty, and (3) a case study in which MASSpy is utilized to overcome issues that arise when integrating experimental data with the computation of functional states of detailed biological mechanisms. MASSpy represents a powerful tool to address challenges that arise in dynamic modeling of metabolic networks, both at small and large scales.</p><p class="para" id="N65542">Genome-scale reconstructions of metabolism appeared shortly after the first genome sequences became available. Constraint-based models are widely used to compute steady state properties of such reconstructions, but the attainment of dynamic models has remained elusive. We thus developed the MASSpy software package, a framework that enables the construction, simulation, and visualization of dynamic metabolic models. MASSpy is based on the mass action kinetics for each elementary step in an enzymatic reaction mechanism. MASSpy seamlessly unites existing software packages within its framework to provide the user with various modeling tools in one package. MASSpy integrates community standards to facilitate the exchange of models, giving modelers the freedom to use the software for different aspects of their own modeling workflows. Furthermore, MASSpy contains methods for generating and simulating ensembles of models, and for explicitly accounting for biological uncertainty. MASSpy has already demonstrated success in a classroom setting. We anticipate that the suite of modeling tools incorporated into MASSpy will enhance the ability of the modeling community to construct and interrogate complex dynamic models of metabolism.</p>]]></description>
            <pubDate><![CDATA[2021-01-28T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[MASSpy: Building, simulating, and visualizing dynamic biological models in Python using mass action kinetics]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765872535383-b8832290-9398-4314-917d-a2593da181f8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008208</link>
            <description><![CDATA[<p class="para" id="N65539">Mathematical models of metabolic networks utilize simulation to study system-level mechanisms and functions. Various approaches have been used to model the steady state behavior of metabolic networks using genome-scale reconstructions, but formulating dynamic models from such reconstructions continues to be a key challenge. Here, we present the Mass Action Stoichiometric Simulation Python (MASSpy) package, an open-source computational framework for dynamic modeling of metabolism. MASSpy utilizes mass action kinetics and detailed chemical mechanisms to build dynamic models of complex biological processes. MASSpy adds dynamic modeling tools to the COnstraint-Based Reconstruction and Analysis Python (COBRApy) package to provide an unified framework for constraint-based and kinetic modeling of metabolic networks. MASSpy supports high-performance dynamic simulation through its implementation of libRoadRunner: the Systems Biology Markup Language (SBML) simulation engine. Three examples are provided to demonstrate how to use MASSpy: (1) a validation of the MASSpy modeling tool through dynamic simulation of detailed mechanisms of enzyme regulation; (2) a feature demonstration using a workflow for generating ensemble of kinetic models using Monte Carlo sampling to approximate missing numerical values of parameters and to quantify biological uncertainty, and (3) a case study in which MASSpy is utilized to overcome issues that arise when integrating experimental data with the computation of functional states of detailed biological mechanisms. MASSpy represents a powerful tool to address challenges that arise in dynamic modeling of metabolic networks, both at small and large scales.</p><p class="para" id="N65542">Genome-scale reconstructions of metabolism appeared shortly after the first genome sequences became available. Constraint-based models are widely used to compute steady state properties of such reconstructions, but the attainment of dynamic models has remained elusive. We thus developed the MASSpy software package, a framework that enables the construction, simulation, and visualization of dynamic metabolic models. MASSpy is based on the mass action kinetics for each elementary step in an enzymatic reaction mechanism. MASSpy seamlessly unites existing software packages within its framework to provide the user with various modeling tools in one package. MASSpy integrates community standards to facilitate the exchange of models, giving modelers the freedom to use the software for different aspects of their own modeling workflows. Furthermore, MASSpy contains methods for generating and simulating ensembles of models, and for explicitly accounting for biological uncertainty. MASSpy has already demonstrated success in a classroom setting. We anticipate that the suite of modeling tools incorporated into MASSpy will enhance the ability of the modeling community to construct and interrogate complex dynamic models of metabolism.</p>]]></description>
            <pubDate><![CDATA[2021-01-28T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Unraveling ChR2-driven stochastic Ca<sup>2+</sup> dynamics in astrocytes: A call for new interventional paradigms]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765872492780-e3d8ad69-0734-4acf-99e7-6dd6a24065b9/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008648</link>
            <description><![CDATA[<p class="para" id="N65539">Optogenetic targeting of astrocytes provides a robust experimental model to differentially induce Ca<sup>2+</sup> signals in astrocytes <i>in vivo</i>. However, a systematic study quantifying the response of optogenetically modified astrocytes to light is yet to be performed. Here, we propose a novel stochastic model of Ca<sup>2+</sup> dynamics in astrocytes that incorporates a light sensitive component—channelrhodopsin 2 (ChR2). Utilizing this model, we investigated the effect of different light stimulation paradigms on cells expressing select variants of ChR2 (wild type, ChETA, and ChRET/TC). Results predict that depending on paradigm specification, astrocytes might undergo drastic changes in their basal Ca<sup>2+</sup> level and spiking probability. Furthermore, we performed a global sensitivity analysis to assess the effect of variation in parameters pertinent to the shape of the ChR2 photocurrent on astrocytic Ca<sup>2+</sup> dynamics. Results suggest that directing variants towards the first open state of the ChR2 photocycle (o<sub>1</sub>) enhances spiking activity in astrocytes during optical stimulation. Evaluation of the effect of Ca<sup>2+</sup> buffering and coupling coefficient in a network of ChR2-expressing astrocytes demonstrated basal level elevations in the stimulated region and propagation of calcium activity to unstimulated cells. Buffering reduced the diffusion range of Ca<sup>2+</sup> within the network, thereby limiting propagation and influencing the activity of astrocytes. Collectively, the framework presented in this study provides valuable information for the selection of light stimulation paradigms that elicit desired astrocytic activity using existing ChR2 constructs, as well as aids in the engineering of future application-oriented optogenetic variants.</p><p class="para" id="N65542">Optogenetics is a novel technique involving the targeted delivery of light-sensitive ion channels like Channelrhodopsin-2 (ChR2) to cells. Recently, this technique has been expanded to non-neuronal cell types, e.g., astrocytes. Optogenetic control of astrocytes <i>in vivo</i> can aid in further understanding the intricacies of their debated roles in the brain. Here, using a mathematical model, we evaluate the effect of short and long-term light stimulation of ChR2-expressing astrocytes on their Ca<sup>2+</sup> spiking activity and basal level. We further investigate how ChR2 gating dynamics, buffering, and coupling coefficient of Ca<sup>2+</sup> influence astrocytic activity in a single cell and a network. Results of our study can serve as a tool for experimental design and engineering of new application-oriented optogenetic constructs targeting astrocytes.</p>]]></description>
            <pubDate><![CDATA[2021-02-10T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Unraveling ChR2-driven stochastic Ca<sup>2+</sup> dynamics in astrocytes: A call for new interventional paradigms]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765872492780-e3d8ad69-0734-4acf-99e7-6dd6a24065b9/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008648</link>
            <description><![CDATA[<p class="para" id="N65539">Optogenetic targeting of astrocytes provides a robust experimental model to differentially induce Ca<sup>2+</sup> signals in astrocytes <i>in vivo</i>. However, a systematic study quantifying the response of optogenetically modified astrocytes to light is yet to be performed. Here, we propose a novel stochastic model of Ca<sup>2+</sup> dynamics in astrocytes that incorporates a light sensitive component—channelrhodopsin 2 (ChR2). Utilizing this model, we investigated the effect of different light stimulation paradigms on cells expressing select variants of ChR2 (wild type, ChETA, and ChRET/TC). Results predict that depending on paradigm specification, astrocytes might undergo drastic changes in their basal Ca<sup>2+</sup> level and spiking probability. Furthermore, we performed a global sensitivity analysis to assess the effect of variation in parameters pertinent to the shape of the ChR2 photocurrent on astrocytic Ca<sup>2+</sup> dynamics. Results suggest that directing variants towards the first open state of the ChR2 photocycle (o<sub>1</sub>) enhances spiking activity in astrocytes during optical stimulation. Evaluation of the effect of Ca<sup>2+</sup> buffering and coupling coefficient in a network of ChR2-expressing astrocytes demonstrated basal level elevations in the stimulated region and propagation of calcium activity to unstimulated cells. Buffering reduced the diffusion range of Ca<sup>2+</sup> within the network, thereby limiting propagation and influencing the activity of astrocytes. Collectively, the framework presented in this study provides valuable information for the selection of light stimulation paradigms that elicit desired astrocytic activity using existing ChR2 constructs, as well as aids in the engineering of future application-oriented optogenetic variants.</p><p class="para" id="N65542">Optogenetics is a novel technique involving the targeted delivery of light-sensitive ion channels like Channelrhodopsin-2 (ChR2) to cells. Recently, this technique has been expanded to non-neuronal cell types, e.g., astrocytes. Optogenetic control of astrocytes <i>in vivo</i> can aid in further understanding the intricacies of their debated roles in the brain. Here, using a mathematical model, we evaluate the effect of short and long-term light stimulation of ChR2-expressing astrocytes on their Ca<sup>2+</sup> spiking activity and basal level. We further investigate how ChR2 gating dynamics, buffering, and coupling coefficient of Ca<sup>2+</sup> influence astrocytic activity in a single cell and a network. Results of our study can serve as a tool for experimental design and engineering of new application-oriented optogenetic constructs targeting astrocytes.</p>]]></description>
            <pubDate><![CDATA[2021-02-10T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[GTP-dependent formation of straight tubulin oligomers leads to microtubule nucleation]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765863535038-e7325ca6-8fca-4a1d-9de1-c13ea94d11e8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202007033</link>
            <description><![CDATA[<p class="para" id="N65540">Ayukawa, Iwata, Imai, et al. visualize the early intermediates in the pathway of spontaneous nucleation of microtubules by using rapid flush negative stain electron microscopy. This study demonstrates that the formation of straight tubulin oligomers of critical size is essential for nucleation.</p><p class="para" id="N65539">Nucleation of microtubules (MTs) is essential for cellular activities, but its mechanism is unknown because of the difficulty involved in capturing rare stochastic events in the early stage of polymerization. Here, combining rapid flush negative stain electron microscopy (EM) and kinetic analysis, we demonstrate that the formation of straight oligomers of critical size is essential for nucleation. Both GDP and GTP tubulin form single-stranded oligomers with a broad range of curvatures, but upon nucleation, the curvature distribution of GTP oligomers is shifted to produce a minor population of straight oligomers. With tubulin having the Y222F mutation in the β subunit, the proportion of straight oligomers increases and nucleation accelerates. Our results support a model in which GTP binding generates a minor population of straight oligomers compatible with lateral association and further growth to MTs. This study suggests that cellular factors involved in nucleation promote it via stabilization of straight oligomers.</p>]]></description>
            <pubDate><![CDATA[2021-02-05T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Contrasting mechanisms for hidden hearing loss: Synaptopathy vs myelin defects]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765853472603-42c9ea5b-f553-43d8-8a5d-7ffb10861d8c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008499</link>
            <description><![CDATA[<p class="para" id="N65539">Hidden hearing loss (HHL) is an auditory neuropathy characterized by normal hearing thresholds but reduced amplitudes of the sound-evoked auditory nerve compound action potential (CAP). In animal models, HHL can be caused by moderate noise exposure or aging, which induces loss of inner hair cell (IHC) synapses. In contrast, recent evidence has shown that transient loss of cochlear Schwann cells also causes permanent auditory deficits in mice with similarities to HHL. Histological analysis of the cochlea after auditory nerve remyelination showed a permanent disruption of the myelination patterns at the heminode of type I spiral ganglion neuron (SGN) peripheral terminals, suggesting that this defect could be contributing to HHL. To shed light on the mechanisms of different HHL scenarios observed in animals and to test their impact on type I SGN activity, we constructed a reduced biophysical model for a population of SGN peripheral axons whose activity is driven by a well-accepted model of cochlear sound processing. We found that the amplitudes of simulated sound-evoked SGN CAPs are lower and have greater latencies when heminodes are disorganized, i.e. they occur at different distances from the hair cell rather than at the same distance as in the normal cochlea. These results confirm that disruption of heminode positions causes desynchronization of SGN spikes leading to a loss of temporal resolution and reduction of the sound-evoked SGN CAP. Another mechanism resulting in HHL is loss of IHC synapses, i.e., synaptopathy. For comparison, we simulated synaptopathy by removing high threshold IHC-SGN synapses and found that the amplitude of simulated sound-evoked SGN CAPs decreases while latencies remain unchanged, as has been observed in noise exposed animals. Thus, model results illuminate diverse disruptions caused by synaptopathy and demyelination on neural activity in auditory processing that contribute to HHL as observed in animal models and that can contribute to perceptual deficits induced by nerve damage in humans.</p><p class="para" id="N65542">Hidden hearing loss (HHL) is an auditory disorder caused by noise exposure, aging or peripheral neuropathy. It is named ‘hidden’ because it does not affect hearing thresholds, but deficits are observed in sound-evoked auditory nerve activity. Studies on animal models suggest two possible pathogenic mechanisms for HHL: (1) loss of synapses between inner hair cells and auditory nerve fibers, and (2) disruption of auditory-nerve myelin. In this study, we constructed a computational model of sound-evoked auditory nerve fiber activity to understand how each of these mechanisms affects nerve compound action potentials. We show that disruption of auditory-nerve myelin desynchronizes sound-evoked auditory neuron spiking, decreasing the amplitude and increasing the latency of the compound action potential. In addition, elongation of the initial axon segment may cause spike generation failure leading to decreased spiking probability. In contrast, synapse loss only decreases the probability of firing, thus reducing the compound action potential amplitude without disturbing its latency. This model, which accurately represents <i>in vivo</i> findings on HHL in the first stages of auditory processing, can be useful to analyze the consequences of synaptopathy and myelinopathy on downstream sound processing and perceptual deficits.</p>]]></description>
            <pubDate><![CDATA[2021-01-22T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Contrasting mechanisms for hidden hearing loss: Synaptopathy vs myelin defects]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765853472603-42c9ea5b-f553-43d8-8a5d-7ffb10861d8c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008499</link>
            <description><![CDATA[<p class="para" id="N65539">Hidden hearing loss (HHL) is an auditory neuropathy characterized by normal hearing thresholds but reduced amplitudes of the sound-evoked auditory nerve compound action potential (CAP). In animal models, HHL can be caused by moderate noise exposure or aging, which induces loss of inner hair cell (IHC) synapses. In contrast, recent evidence has shown that transient loss of cochlear Schwann cells also causes permanent auditory deficits in mice with similarities to HHL. Histological analysis of the cochlea after auditory nerve remyelination showed a permanent disruption of the myelination patterns at the heminode of type I spiral ganglion neuron (SGN) peripheral terminals, suggesting that this defect could be contributing to HHL. To shed light on the mechanisms of different HHL scenarios observed in animals and to test their impact on type I SGN activity, we constructed a reduced biophysical model for a population of SGN peripheral axons whose activity is driven by a well-accepted model of cochlear sound processing. We found that the amplitudes of simulated sound-evoked SGN CAPs are lower and have greater latencies when heminodes are disorganized, i.e. they occur at different distances from the hair cell rather than at the same distance as in the normal cochlea. These results confirm that disruption of heminode positions causes desynchronization of SGN spikes leading to a loss of temporal resolution and reduction of the sound-evoked SGN CAP. Another mechanism resulting in HHL is loss of IHC synapses, i.e., synaptopathy. For comparison, we simulated synaptopathy by removing high threshold IHC-SGN synapses and found that the amplitude of simulated sound-evoked SGN CAPs decreases while latencies remain unchanged, as has been observed in noise exposed animals. Thus, model results illuminate diverse disruptions caused by synaptopathy and demyelination on neural activity in auditory processing that contribute to HHL as observed in animal models and that can contribute to perceptual deficits induced by nerve damage in humans.</p><p class="para" id="N65542">Hidden hearing loss (HHL) is an auditory disorder caused by noise exposure, aging or peripheral neuropathy. It is named ‘hidden’ because it does not affect hearing thresholds, but deficits are observed in sound-evoked auditory nerve activity. Studies on animal models suggest two possible pathogenic mechanisms for HHL: (1) loss of synapses between inner hair cells and auditory nerve fibers, and (2) disruption of auditory-nerve myelin. In this study, we constructed a computational model of sound-evoked auditory nerve fiber activity to understand how each of these mechanisms affects nerve compound action potentials. We show that disruption of auditory-nerve myelin desynchronizes sound-evoked auditory neuron spiking, decreasing the amplitude and increasing the latency of the compound action potential. In addition, elongation of the initial axon segment may cause spike generation failure leading to decreased spiking probability. In contrast, synapse loss only decreases the probability of firing, thus reducing the compound action potential amplitude without disturbing its latency. This model, which accurately represents <i>in vivo</i> findings on HHL in the first stages of auditory processing, can be useful to analyze the consequences of synaptopathy and myelinopathy on downstream sound processing and perceptual deficits.</p>]]></description>
            <pubDate><![CDATA[2021-01-22T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Individual kinetochore-fibers locally dissipate force to maintain robust mammalian spindle structure]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765850287115-09ba7ec2-6de4-4494-a4ca-578b060de55e/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.201911090</link>
            <description><![CDATA[<p class="para" id="N65540">To segregate chromosomes, the mammalian spindle must generate and respond to force. How it does so remains poorly understood. Pulling on the spindle using microneedles, Long et al. show that it can locally dissipate sustained force by regulating microtubule dynamics and breakage, thereby preserving global spindle structure.</p><p class="para" id="N65539">At cell division, the mammalian kinetochore binds many spindle microtubules that make up the kinetochore-fiber. To segregate chromosomes, the kinetochore-fiber must be dynamic and generate and respond to force. Yet, how it remodels under force remains poorly understood. Kinetochore-fibers cannot be reconstituted in vitro, and exerting controlled forces in vivo remains challenging. Here, we use microneedles to pull on mammalian kinetochore-fibers and probe how sustained force regulates their dynamics and structure. We show that force lengthens kinetochore-fibers by persistently favoring plus-end polymerization, not by increasing polymerization rate. We demonstrate that force suppresses depolymerization at both plus and minus ends, rather than sliding microtubules within the kinetochore-fiber. Finally, we observe that kinetochore-fibers break but do not detach from kinetochores or poles. Together, this work suggests an engineering principle for spindle structural homeostasis: different physical mechanisms of local force dissipation by the k-fiber limit force transmission to preserve robust spindle structure. These findings may inform how other dynamic, force-generating cellular machines achieve mechanical robustness.</p><p class="para" id="N65542"><div class="section" id="GA"><div class="img"><div class="imgeVideo"><div class="img-fullscreenIcon" onClick="javascript:showImageContent('GA');"><img src="/public/images/journalImg/fullscreen.png"/></div><div class="imageVideo"><img src="/dataresources/secured/content-1765850287115-09ba7ec2-6de4-4494-a4ca-578b060de55e/assets/JCB_201911090_GA.jpg" alt=""/></div></div></div></div></p>]]></description>
            <pubDate><![CDATA[2020-05-20T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Single stranded DNA annealing is a conserved activity of telomere resolvases]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765838053381-adc1ea71-2675-42cc-8825-dbd287a873c1/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246212</link>
            <description><![CDATA[<p class="para" id="N65539">Bacterial species of the genera <i>Agrobacterium</i> and <i>Borrelia</i> possess chromosomes terminated by hairpin telomeres. Replication produces dimeric replication intermediates fused via replicated telomere junctions. A specialized class of enzymes, referred to as telomere resolvases, promotes the resolution of the replicated intermediate into linear monomers terminated by hairpin telomeres. Telomere resolution is catalyzed via DNA cleavage and rejoining events mechanistically similar to those promoted by topoisomerase-IB and tyrosine recombinase enzymes. Examination of the borrelial telomere resolvase, ResT, revealed unanticipated multifunctionality; aside from its expected telomere resolution activity ResT possessed a singled-stranded DNA (ssDNA) annealing activity that extended to both naked ssDNA and ssDNA complexed with its cognate single-stranded DNA binding protein (SSB). At present, the role this DNA annealing activity plays <i>in vivo</i> remains unknown. We have demonstrated here that single-stranded DNA annealing is also a conserved property of the agrobacterial telomere resolvase, TelA. This activity in TelA similarly extends to both naked ssDNA and ssDNA bound by its cognate SSB. TelA’s annealing activity was shown to stem from the N-terminal domain; removal of this domain abolished annealing without affecting telomere resolution. Further, independent expression of the N-terminal domain of TelA produced a functional annealing protein. We suggest that the apparent conservation of annealing activity in two telomere resolvases, from distantly related bacterial species, implies a role for this activity in hairpin telomere metabolism. Our demonstration of the separation of the telomere resolution and annealing activities of TelA provides a platform for future experiments aimed at identifying the role DNA annealing performs <i>in vivo</i>.</p>]]></description>
            <pubDate><![CDATA[2021-02-04T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Single stranded DNA annealing is a conserved activity of telomere resolvases]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765838053381-adc1ea71-2675-42cc-8825-dbd287a873c1/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0246212</link>
            <description><![CDATA[<p class="para" id="N65539">Bacterial species of the genera <i>Agrobacterium</i> and <i>Borrelia</i> possess chromosomes terminated by hairpin telomeres. Replication produces dimeric replication intermediates fused via replicated telomere junctions. A specialized class of enzymes, referred to as telomere resolvases, promotes the resolution of the replicated intermediate into linear monomers terminated by hairpin telomeres. Telomere resolution is catalyzed via DNA cleavage and rejoining events mechanistically similar to those promoted by topoisomerase-IB and tyrosine recombinase enzymes. Examination of the borrelial telomere resolvase, ResT, revealed unanticipated multifunctionality; aside from its expected telomere resolution activity ResT possessed a singled-stranded DNA (ssDNA) annealing activity that extended to both naked ssDNA and ssDNA complexed with its cognate single-stranded DNA binding protein (SSB). At present, the role this DNA annealing activity plays <i>in vivo</i> remains unknown. We have demonstrated here that single-stranded DNA annealing is also a conserved property of the agrobacterial telomere resolvase, TelA. This activity in TelA similarly extends to both naked ssDNA and ssDNA bound by its cognate SSB. TelA’s annealing activity was shown to stem from the N-terminal domain; removal of this domain abolished annealing without affecting telomere resolution. Further, independent expression of the N-terminal domain of TelA produced a functional annealing protein. We suggest that the apparent conservation of annealing activity in two telomere resolvases, from distantly related bacterial species, implies a role for this activity in hairpin telomere metabolism. Our demonstration of the separation of the telomere resolution and annealing activities of TelA provides a platform for future experiments aimed at identifying the role DNA annealing performs <i>in vivo</i>.</p>]]></description>
            <pubDate><![CDATA[2021-02-04T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[PGE2 upregulates the Na<sup>+</sup>/K<sup>+</sup> ATPase in HepG2 cells via EP4 receptors and intracellular calcium]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765837333795-1f82cdc5-7cc3-4b8e-bf02-8bd836ab8f6c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245400</link>
            <description><![CDATA[<p class="para" id="N65539">The Na<sup>+</sup>/K<sup>+</sup> ATPase is a key regulator of the hepatocytes ionic homeostasis, which when altered may lead to many liver disorders. We demonstrated recently, a significant stimulation of the Na<sup>+</sup>/K<sup>+</sup> ATPase in HepG2 cells treated with the S1P analogue FTY 720P, that was mediated through PGE2. The mechanism by which the prostaglandin exerts its effect was not investigated, and is the focus of this work. The type of receptors involved was determined using pharmacological inhibitors, while western blot analysis, fluorescence imaging of GFP-tagged Na<sup>+</sup>/K<sup>+</sup> ATPase, and time-lapse imaging on live cells were used to detect changes in membrane abundance of the Na<sup>+</sup>/K<sup>+</sup> ATPase. The activity of the ATPase was assayed by measuring the amount of inorganic phosphate liberated in the presence and absence of ouabain. The enhanced activity of the ATPase was not observed when EP4 receptors were blocked but still appeared in presence inhibitors of EP1, EP2 and EP3 receptors. The involvement of EP4 was confirmed by the stimulation observed with EP4 agonist. The stimulatory effect of PGE2 did not appear in presence of Rp-cAMP, an inhibitor of PKA, and was imitated by db-cAMP, a PKA activator. Chelating intracellular calcium with BAPTA-AM abrogated the effect of db-cAMP as well as that of PGE2, but PGE2 treatment in a calcium-free PBS medium did not, suggesting an involvement of intracellular calcium, that was confirmed by the results obtained with 2-APB treatment. Live cell imaging showed movement of GFP–Na+/K+ ATPase-positive vesicles to the membrane and increased abundance of the ATPase at the membrane after PGE2 treatment. It was concluded that PGE2 acts via EP4, PKA, and intracellular calcium.</p>]]></description>
            <pubDate><![CDATA[2021-01-14T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[PGE2 upregulates the Na<sup>+</sup>/K<sup>+</sup> ATPase in HepG2 cells via EP4 receptors and intracellular calcium]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765837333795-1f82cdc5-7cc3-4b8e-bf02-8bd836ab8f6c/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245400</link>
            <description><![CDATA[<p class="para" id="N65539">The Na<sup>+</sup>/K<sup>+</sup> ATPase is a key regulator of the hepatocytes ionic homeostasis, which when altered may lead to many liver disorders. We demonstrated recently, a significant stimulation of the Na<sup>+</sup>/K<sup>+</sup> ATPase in HepG2 cells treated with the S1P analogue FTY 720P, that was mediated through PGE2. The mechanism by which the prostaglandin exerts its effect was not investigated, and is the focus of this work. The type of receptors involved was determined using pharmacological inhibitors, while western blot analysis, fluorescence imaging of GFP-tagged Na<sup>+</sup>/K<sup>+</sup> ATPase, and time-lapse imaging on live cells were used to detect changes in membrane abundance of the Na<sup>+</sup>/K<sup>+</sup> ATPase. The activity of the ATPase was assayed by measuring the amount of inorganic phosphate liberated in the presence and absence of ouabain. The enhanced activity of the ATPase was not observed when EP4 receptors were blocked but still appeared in presence inhibitors of EP1, EP2 and EP3 receptors. The involvement of EP4 was confirmed by the stimulation observed with EP4 agonist. The stimulatory effect of PGE2 did not appear in presence of Rp-cAMP, an inhibitor of PKA, and was imitated by db-cAMP, a PKA activator. Chelating intracellular calcium with BAPTA-AM abrogated the effect of db-cAMP as well as that of PGE2, but PGE2 treatment in a calcium-free PBS medium did not, suggesting an involvement of intracellular calcium, that was confirmed by the results obtained with 2-APB treatment. Live cell imaging showed movement of GFP–Na+/K+ ATPase-positive vesicles to the membrane and increased abundance of the ATPase at the membrane after PGE2 treatment. It was concluded that PGE2 acts via EP4, PKA, and intracellular calcium.</p>]]></description>
            <pubDate><![CDATA[2021-01-14T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Cerebellar Golgi cell models predict dendritic processing and mechanisms of synaptic plasticity]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765837079740-7613dcd2-7aa3-49cc-8577-5cdbcffa376b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007937</link>
            <description><![CDATA[<p class="para" id="N65539">The Golgi cells are the main inhibitory interneurons of the cerebellar granular layer. Although recent works have highlighted the complexity of their dendritic organization and synaptic inputs, the mechanisms through which these neurons integrate complex input patterns remained unknown. Here we have used 8 detailed morphological reconstructions to develop multicompartmental models of Golgi cells, in which Na, Ca, and K channels were distributed along dendrites, soma, axonal initial segment and axon. The models faithfully reproduced a rich pattern of electrophysiological and pharmacological properties and predicted the operating mechanisms of these neurons. Basal dendrites turned out to be more tightly electrically coupled to the axon initial segment than apical dendrites. During synaptic transmission, parallel fibers caused slow Ca-dependent depolarizations in apical dendrites that boosted the axon initial segment encoder and Na-spike backpropagation into basal dendrites, while inhibitory synapses effectively shunted backpropagating currents. This oriented dendritic processing set up a coincidence detector controlling voltage-dependent NMDA receptor unblock in basal dendrites, which, by regulating local calcium influx, may provide the basis for spike-timing dependent plasticity anticipated by theory.</p><p class="para" id="N65542">The Golgi cells are the main inhibitory interneurons of the cerebellum granular layer and play a fundamental role in controlling cerebellar processing. However, it was unclear how spikes are processed in the dendrites by specific sets of ionic channels and how they might contribute to integrate synaptic inputs and plasticity. Here we have developed detailed multicompartmental models of Golgi cells that faithfully reproduced a large set of experimental findings and revealed the nature of signal interchange between dendrites and axo-somatic compartments. A main prediction of the models is that synaptic activation of apical dendrites can effectively trigger spike generation in the axonal initial segment followed by rapid spike backpropagation into basal dendrites. Here, incoming mossy fiber inputs and backpropagating spikes regulate the voltage-dependent unblock of NMDA channels and the induction of spike timing-dependent plasticity (STDP). STDP, which was predicted by theory, may therefore be controlled by contextual information provided by parallel fibers and integrated in apical dendrites, supporting the view that spike timing is fundamental to control synaptic plasticity at the cerebellar input stage.</p>]]></description>
            <pubDate><![CDATA[2020-12-30T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Cerebellar Golgi cell models predict dendritic processing and mechanisms of synaptic plasticity]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765837079740-7613dcd2-7aa3-49cc-8577-5cdbcffa376b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007937</link>
            <description><![CDATA[<p class="para" id="N65539">The Golgi cells are the main inhibitory interneurons of the cerebellar granular layer. Although recent works have highlighted the complexity of their dendritic organization and synaptic inputs, the mechanisms through which these neurons integrate complex input patterns remained unknown. Here we have used 8 detailed morphological reconstructions to develop multicompartmental models of Golgi cells, in which Na, Ca, and K channels were distributed along dendrites, soma, axonal initial segment and axon. The models faithfully reproduced a rich pattern of electrophysiological and pharmacological properties and predicted the operating mechanisms of these neurons. Basal dendrites turned out to be more tightly electrically coupled to the axon initial segment than apical dendrites. During synaptic transmission, parallel fibers caused slow Ca-dependent depolarizations in apical dendrites that boosted the axon initial segment encoder and Na-spike backpropagation into basal dendrites, while inhibitory synapses effectively shunted backpropagating currents. This oriented dendritic processing set up a coincidence detector controlling voltage-dependent NMDA receptor unblock in basal dendrites, which, by regulating local calcium influx, may provide the basis for spike-timing dependent plasticity anticipated by theory.</p><p class="para" id="N65542">The Golgi cells are the main inhibitory interneurons of the cerebellum granular layer and play a fundamental role in controlling cerebellar processing. However, it was unclear how spikes are processed in the dendrites by specific sets of ionic channels and how they might contribute to integrate synaptic inputs and plasticity. Here we have developed detailed multicompartmental models of Golgi cells that faithfully reproduced a large set of experimental findings and revealed the nature of signal interchange between dendrites and axo-somatic compartments. A main prediction of the models is that synaptic activation of apical dendrites can effectively trigger spike generation in the axonal initial segment followed by rapid spike backpropagation into basal dendrites. Here, incoming mossy fiber inputs and backpropagating spikes regulate the voltage-dependent unblock of NMDA channels and the induction of spike timing-dependent plasticity (STDP). STDP, which was predicted by theory, may therefore be controlled by contextual information provided by parallel fibers and integrated in apical dendrites, supporting the view that spike timing is fundamental to control synaptic plasticity at the cerebellar input stage.</p>]]></description>
            <pubDate><![CDATA[2020-12-30T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Minimal biophysical model of combined antibiotic action]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765836611622-62530cb0-96b8-4f03-a0a5-3c76f86a109b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008529</link>
            <description><![CDATA[<p class="para" id="N65539">Phenomenological relations such as Ohm’s or Fourier’s law have a venerable history in physics but are still scarce in biology. This situation restrains predictive theory. Here, we build on bacterial “growth laws,” which capture physiological feedback between translation and cell growth, to construct a minimal biophysical model for the combined action of ribosome-targeting antibiotics. Our model predicts drug interactions like antagonism or synergy solely from responses to individual drugs. We provide analytical results for limiting cases, which agree well with numerical results. We systematically refine the model by including direct physical interactions of different antibiotics on the ribosome. In a limiting case, our model provides a mechanistic underpinning for recent predictions of higher-order interactions that were derived using entropy maximization. We further refine the model to include the effects of antibiotics that mimic starvation and the presence of resistance genes. We describe the impact of a starvation-mimicking antibiotic on drug interactions analytically and verify it experimentally. Our extended model suggests a change in the type of drug interaction that depends on the strength of resistance, which challenges established rescaling paradigms. We experimentally show that the presence of unregulated resistance genes can lead to altered drug interaction, which agrees with the prediction of the model. While minimal, the model is readily adaptable and opens the door to predicting interactions of second and higher-order in a broad range of biological systems.</p><p class="para" id="N65542">Applying multiple antibiotics simultaneously can boost treatment effectiveness and aid against rampant antibiotic resistance. Because of the impractically large number of possible combinations of drugs, those that are effective are found by trial and error. Hence, a predictive theory to characterize drug cocktails would be of enormous value. Recently identified phenomenological laws ease the construction of predictive models of bacterial growth. Here, we build a model of the effects of antibiotic combinations on bacteria and show that it makes reliable predictions for experimental outcomes. Our model takes responses to individual drugs as inputs and predicts their combined effect. This output determines the type of drug interaction, which can range from antagonistic (the combined effect is weaker) to synergistic (the combined effect is stronger). We broaden the model by including the direct physical interaction on the target, drug resistance genes that alter the drug interaction, and drugs that mimic poor growth environments by choking the supply of growth-essential components, which we test experimentally. Our results prove how biophysical models that use empirical laws can predict responses to drug combinations. Importantly, such models can successfully predict mechanisms underlying interactions of drug combinations. This approach is extensible to combinations of more than two drugs and diverse biological systems.</p>]]></description>
            <pubDate><![CDATA[2021-01-07T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Minimal biophysical model of combined antibiotic action]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765836611622-62530cb0-96b8-4f03-a0a5-3c76f86a109b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008529</link>
            <description><![CDATA[<p class="para" id="N65539">Phenomenological relations such as Ohm’s or Fourier’s law have a venerable history in physics but are still scarce in biology. This situation restrains predictive theory. Here, we build on bacterial “growth laws,” which capture physiological feedback between translation and cell growth, to construct a minimal biophysical model for the combined action of ribosome-targeting antibiotics. Our model predicts drug interactions like antagonism or synergy solely from responses to individual drugs. We provide analytical results for limiting cases, which agree well with numerical results. We systematically refine the model by including direct physical interactions of different antibiotics on the ribosome. In a limiting case, our model provides a mechanistic underpinning for recent predictions of higher-order interactions that were derived using entropy maximization. We further refine the model to include the effects of antibiotics that mimic starvation and the presence of resistance genes. We describe the impact of a starvation-mimicking antibiotic on drug interactions analytically and verify it experimentally. Our extended model suggests a change in the type of drug interaction that depends on the strength of resistance, which challenges established rescaling paradigms. We experimentally show that the presence of unregulated resistance genes can lead to altered drug interaction, which agrees with the prediction of the model. While minimal, the model is readily adaptable and opens the door to predicting interactions of second and higher-order in a broad range of biological systems.</p><p class="para" id="N65542">Applying multiple antibiotics simultaneously can boost treatment effectiveness and aid against rampant antibiotic resistance. Because of the impractically large number of possible combinations of drugs, those that are effective are found by trial and error. Hence, a predictive theory to characterize drug cocktails would be of enormous value. Recently identified phenomenological laws ease the construction of predictive models of bacterial growth. Here, we build a model of the effects of antibiotic combinations on bacteria and show that it makes reliable predictions for experimental outcomes. Our model takes responses to individual drugs as inputs and predicts their combined effect. This output determines the type of drug interaction, which can range from antagonistic (the combined effect is weaker) to synergistic (the combined effect is stronger). We broaden the model by including the direct physical interaction on the target, drug resistance genes that alter the drug interaction, and drugs that mimic poor growth environments by choking the supply of growth-essential components, which we test experimentally. Our results prove how biophysical models that use empirical laws can predict responses to drug combinations. Importantly, such models can successfully predict mechanisms underlying interactions of drug combinations. This approach is extensible to combinations of more than two drugs and diverse biological systems.</p>]]></description>
            <pubDate><![CDATA[2021-01-07T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Liquid–liquid phase separation in autophagy]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765833672689-96c39090-ae35-450f-afbb-5136412ed93b/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1083/jcb.202004062</link>
            <description><![CDATA[<p class="para" id="N65540">Zhang and colleagues summarize phase separation and transition in the assembly of autophagosome formation sites and triage of protein condensates for degradation.</p><p class="para" id="N65539">Liquid–liquid phase separation (LLPS) compartmentalizes and concentrates biomacromolecules into distinct condensates. Liquid-like condensates can transition into gel and solid states, which are essential for fulfilling their different functions. LLPS plays important roles in multiple steps of autophagy, mediating the assembly of autophagosome formation sites, acting as an unconventional modulator of TORC1-mediated autophagy regulation, and triaging protein cargos for degradation. Gel-like, but not solid, protein condensates can trigger formation of surrounding autophagosomal membranes. Stress and pathological conditions cause aberrant phase separation and transition of condensates, which can evade surveillance by the autophagy machinery. Understanding the mechanisms underlying phase separation and transition will provide potential therapeutic targets for protein aggregation diseases.</p>]]></description>
            <pubDate><![CDATA[2020-06-30T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Aberrant regulation of a poison exon caused by a non-coding variant in a mouse model of <i>Scn1a</i>-associated epileptic encephalopathy]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765821971021-c9e4e2fe-42bb-4639-a761-7a9f05ab5dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pgen.1009195</link>
            <description><![CDATA[<p class="para" id="N65539">Dravet syndrome (DS) is a developmental and epileptic encephalopathy that results from mutations in the Na<sub>v</sub>1.1 sodium channel encoded by <i>SCN1A</i>. Most known DS-causing mutations are in coding regions of <i>SCN1A</i>, but we recently identified several disease-associated <i>SCN1A</i> mutations in intron 20 that are within or near to a cryptic and evolutionarily conserved “poison” exon, 20N, whose inclusion is predicted to lead to transcript degradation. However, it is not clear how these intron 20 variants alter <i>SCN1A</i> expression or DS pathophysiology in an organismal context, nor is it clear how exon 20N is regulated in a tissue-specific and developmental context. We address those questions here by generating an animal model of our index case, NM_006920.4(SCN1A):c.3969+2451G&gt;C, using gene editing to create the orthologous mutation in laboratory mice. <i>Scn1a</i> heterozygous knock-in (+/<i>KI</i>) mice exhibited an ~50% reduction in brain <i>Scn1a</i> mRNA and Na<sub>v</sub>1.1 protein levels, together with characteristics observed in other DS mouse models, including premature mortality, seizures, and hyperactivity. In brain tissue from adult <i>Scn1a</i> +/+ animals, quantitative RT-PCR assays indicated that ~1% of <i>Scn1a</i> mRNA included exon 20N, while brain tissue from <i>Scn1a +/KI</i> mice exhibited an ~5-fold increase in the extent of exon 20N inclusion. We investigated the extent of exon 20N inclusion in brain during normal fetal development in RNA-seq data and discovered that levels of inclusion were ~70% at E14.5, declining progressively to ~10% postnatally. A similar pattern exists for the homologous sodium channel Na<sub>v</sub>1.6, encoded by <i>Scn8a</i>. For both genes, there is an inverse relationship between the level of functional transcript and the extent of poison exon inclusion. Taken together, our findings suggest that poison exon usage by <i>Scn1a</i> and <i>Scn8a</i> is a strategy to regulate channel expression during normal brain development, and that mutations recapitulating a fetal-like pattern of splicing cause reduced channel expression and epileptic encephalopathy.</p><p class="para" id="N65542">Dravet syndrome (DS) is a neurological disorder affecting approximately 1:15,700 Americans that causes generalized epilepsy and associated complications. While most patients have a mutation in the <i>SCN1A</i> gene that encodes the Na<sub>v</sub>1.1 voltage-gated sodium channel, about 20% do not have a mutation identified by exome or targeted sequencing. Recently, we identified variants in intron 20, a noncoding region of <i>SCN1A</i>, in some DS patients. We hypothesized that these variants alter <i>SCN1A</i> transcript processing, decrease Na<sub>v</sub>1.1 function, and lead to DS pathophysiology via inclusion of exon 20N, a “poison” exon that leads to a premature stop codon. In this study, we generated a knock-in mouse model, <i>Scn1a+/KI</i>, of one of these variants, which resides in a genomic region that is extremely conserved across vertebrate species. We found that <i>Scn1a+/KI</i> mice have reduced levels of <i>Scn1a</i> transcript and Na<sub>v</sub>1.1 protein and develop DS-related phenotypes. We find that transcripts from brains of <i>Scn1a+/KI</i> mice show elevated rates of <i>Scn1a</i> exon 20N inclusion. We also explored the relationship between exon 20N inclusion and <i>Scn1a</i> expression during development, and found that, during brain development when <i>Scn1a</i> expression is low, exon 20N inclusion is high; postnatally, as <i>Scn1a</i> expression increases, there is a corresponding decrease in exon 20N usage. Expression of another voltage-gated sodium channel transcript, <i>Scn8a</i> (Na<sub>v</sub>1.6), was similarly regulated. Together, these data demonstrate that poison exon inclusion is a conserved mechanism to control sodium channel expression in the brain, and that an intronic mutation that disrupts the normal developmental regulation of poison exon inclusion leads to reduced Na<sub>v</sub>1.1 and generalized epilepsy.</p>]]></description>
            <pubDate><![CDATA[2021-01-07T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Aberrant regulation of a poison exon caused by a non-coding variant in a mouse model of <i>Scn1a</i>-associated epileptic encephalopathy]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765821971021-c9e4e2fe-42bb-4639-a761-7a9f05ab5dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pgen.1009195</link>
            <description><![CDATA[<p class="para" id="N65539">Dravet syndrome (DS) is a developmental and epileptic encephalopathy that results from mutations in the Na<sub>v</sub>1.1 sodium channel encoded by <i>SCN1A</i>. Most known DS-causing mutations are in coding regions of <i>SCN1A</i>, but we recently identified several disease-associated <i>SCN1A</i> mutations in intron 20 that are within or near to a cryptic and evolutionarily conserved “poison” exon, 20N, whose inclusion is predicted to lead to transcript degradation. However, it is not clear how these intron 20 variants alter <i>SCN1A</i> expression or DS pathophysiology in an organismal context, nor is it clear how exon 20N is regulated in a tissue-specific and developmental context. We address those questions here by generating an animal model of our index case, NM_006920.4(SCN1A):c.3969+2451G&gt;C, using gene editing to create the orthologous mutation in laboratory mice. <i>Scn1a</i> heterozygous knock-in (+/<i>KI</i>) mice exhibited an ~50% reduction in brain <i>Scn1a</i> mRNA and Na<sub>v</sub>1.1 protein levels, together with characteristics observed in other DS mouse models, including premature mortality, seizures, and hyperactivity. In brain tissue from adult <i>Scn1a</i> +/+ animals, quantitative RT-PCR assays indicated that ~1% of <i>Scn1a</i> mRNA included exon 20N, while brain tissue from <i>Scn1a +/KI</i> mice exhibited an ~5-fold increase in the extent of exon 20N inclusion. We investigated the extent of exon 20N inclusion in brain during normal fetal development in RNA-seq data and discovered that levels of inclusion were ~70% at E14.5, declining progressively to ~10% postnatally. A similar pattern exists for the homologous sodium channel Na<sub>v</sub>1.6, encoded by <i>Scn8a</i>. For both genes, there is an inverse relationship between the level of functional transcript and the extent of poison exon inclusion. Taken together, our findings suggest that poison exon usage by <i>Scn1a</i> and <i>Scn8a</i> is a strategy to regulate channel expression during normal brain development, and that mutations recapitulating a fetal-like pattern of splicing cause reduced channel expression and epileptic encephalopathy.</p><p class="para" id="N65542">Dravet syndrome (DS) is a neurological disorder affecting approximately 1:15,700 Americans that causes generalized epilepsy and associated complications. While most patients have a mutation in the <i>SCN1A</i> gene that encodes the Na<sub>v</sub>1.1 voltage-gated sodium channel, about 20% do not have a mutation identified by exome or targeted sequencing. Recently, we identified variants in intron 20, a noncoding region of <i>SCN1A</i>, in some DS patients. We hypothesized that these variants alter <i>SCN1A</i> transcript processing, decrease Na<sub>v</sub>1.1 function, and lead to DS pathophysiology via inclusion of exon 20N, a “poison” exon that leads to a premature stop codon. In this study, we generated a knock-in mouse model, <i>Scn1a+/KI</i>, of one of these variants, which resides in a genomic region that is extremely conserved across vertebrate species. We found that <i>Scn1a+/KI</i> mice have reduced levels of <i>Scn1a</i> transcript and Na<sub>v</sub>1.1 protein and develop DS-related phenotypes. We find that transcripts from brains of <i>Scn1a+/KI</i> mice show elevated rates of <i>Scn1a</i> exon 20N inclusion. We also explored the relationship between exon 20N inclusion and <i>Scn1a</i> expression during development, and found that, during brain development when <i>Scn1a</i> expression is low, exon 20N inclusion is high; postnatally, as <i>Scn1a</i> expression increases, there is a corresponding decrease in exon 20N usage. Expression of another voltage-gated sodium channel transcript, <i>Scn8a</i> (Na<sub>v</sub>1.6), was similarly regulated. Together, these data demonstrate that poison exon inclusion is a conserved mechanism to control sodium channel expression in the brain, and that an intronic mutation that disrupts the normal developmental regulation of poison exon inclusion leads to reduced Na<sub>v</sub>1.1 and generalized epilepsy.</p>]]></description>
            <pubDate><![CDATA[2021-01-07T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Simulation code for estimating external gamma-ray doses from a radioactive plume and contaminated ground using a local-scale atmospheric dispersion model]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765799748295-fa7ce1f1-d1cd-4a62-9e79-84d9d9e6d3bb/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245932</link>
            <description><![CDATA[<p class="para" id="N65539">In this study, we developed a simulation code powered by lattice dose-response functions (hereinafter SIBYL), which helps in the quick and accurate estimation of external gamma-ray doses emitted from a radioactive plume and contaminated ground. SIBYL couples with atmospheric dispersion models and calculates gamma-ray dose distributions inside a target area based on a map of activity concentrations using pre-evaluated dose-response functions. Moreover, SIBYL considers radiation shielding due to obstructions such as buildings. To examine the reliability of SIBYL, we investigated five typical cases for steady-state and unsteady-state plume dispersions by coupling the Gaussian plume model and the local-scale high-resolution atmospheric dispersion model using large eddy simulation. The results of this coupled model were compared with those of full Monte Carlo simulations using the particle and heavy-ion transport code system (PHITS). The dose-distribution maps calculated using SIBYL differed by up to 10% from those calculated using PHITS in most target locations. The exceptions were locations far from the radioactive contamination and those behind the intricate structures of building arrays. In addition, SIBYL’s computation time using 96 parallel processing elements was several tens of minutes even for the most computationally expensive tasks of this study. The computation using SIBYL was approximately 100 times faster than the same calculation using PHITS under the same computation conditions. From the results of the case studies, we concluded that SIBYL can estimate a ground-level dose-distribution map within one hour with accuracy that is comparable to that of the full Monte Carlo simulation.</p>]]></description>
            <pubDate><![CDATA[2021-01-25T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Simulation code for estimating external gamma-ray doses from a radioactive plume and contaminated ground using a local-scale atmospheric dispersion model]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765799748295-fa7ce1f1-d1cd-4a62-9e79-84d9d9e6d3bb/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245932</link>
            <description><![CDATA[<p class="para" id="N65539">In this study, we developed a simulation code powered by lattice dose-response functions (hereinafter SIBYL), which helps in the quick and accurate estimation of external gamma-ray doses emitted from a radioactive plume and contaminated ground. SIBYL couples with atmospheric dispersion models and calculates gamma-ray dose distributions inside a target area based on a map of activity concentrations using pre-evaluated dose-response functions. Moreover, SIBYL considers radiation shielding due to obstructions such as buildings. To examine the reliability of SIBYL, we investigated five typical cases for steady-state and unsteady-state plume dispersions by coupling the Gaussian plume model and the local-scale high-resolution atmospheric dispersion model using large eddy simulation. The results of this coupled model were compared with those of full Monte Carlo simulations using the particle and heavy-ion transport code system (PHITS). The dose-distribution maps calculated using SIBYL differed by up to 10% from those calculated using PHITS in most target locations. The exceptions were locations far from the radioactive contamination and those behind the intricate structures of building arrays. In addition, SIBYL’s computation time using 96 parallel processing elements was several tens of minutes even for the most computationally expensive tasks of this study. The computation using SIBYL was approximately 100 times faster than the same calculation using PHITS under the same computation conditions. From the results of the case studies, we concluded that SIBYL can estimate a ground-level dose-distribution map within one hour with accuracy that is comparable to that of the full Monte Carlo simulation.</p>]]></description>
            <pubDate><![CDATA[2021-01-25T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Risk factors of heart failure among patients with hypertension attending a tertiary hospital in Ibadan, Nigeria: The RISK-HHF case-control study]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765799100330-8329fb7b-aa83-4fef-9ff0-b36be9d83786/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245734</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Aim</h3><p class="para" id="N65543">Hypertension is the leading cause of heart failure (HF) in sub-Saharan Africa. Preventive public health approach to reduce the scourge of HF must seek to understand the risk factors of HF in at-risk populations. The aim of this study was to characterize the risk factors of HF among patients with hypertension attending a cardiology clinic.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Methods and results</h3><p class="para" id="N65549">One hundred and one (101) case-control age- and sex-matched pairs were recruited. The study population were adults with a clinical diagnosis of hypertensive HF (cases) and individuals with systemic hypertension without HF. They were interviewed and evaluated for cardiovascular risk factors. Associations between variables were tested with chi square test, Fisher’s exact test and independent sample t test as appropriate. Logistic regression modelling was used to determine the independent risk factors of hypertensive HF (HHF) in the study population while ‘punafcc’ package in stata12 was used to calculate the population attributable fraction (PAF) of the risk factors. Suboptimal medication adherence was the strongest adverse risk factor of HHF (medium adherence aOR: 3.53, 95%CI: 1.35–9.25; low adherence aOR: 9.44, 95%CI: 3.41–26.10) with a PAF of 67% followed by dipstick proteinuria (aOR: 4.22, 95%CI: 1.62–11.02; PAF: 34%) and alcohol consumption/day per 10grams (aOR: 1.23, 95%CI: 1.02–1.49; PAF: 22%). The protective risk factors of HHF were use of calcium channel blockers (aOR 0.25, 95%CI: 0.11–0.59; PAF: 59%), then daily fruits and vegetable consumption (aOR 0.41, 95%CI: 0.17–1.01; PAF: 46%), and eGFR (aOR 0.98, 95%CI: 0.96–0.99; PAF: 5.3%).</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Conclusions</h3><p class="para" id="N65555">The risk factors of HHF are amenable to lifestyle and dietary changes. Public health interventions and preventive cardiovascular care to improve medication adherence, promote fruit and vegetable consumption and reduce alcohol consumption among patients with hypertension are recommended. Renoprotection has utility in the prevention of HF among hypertensives.</p></div>]]></description>
            <pubDate><![CDATA[2021-01-25T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Risk factors of heart failure among patients with hypertension attending a tertiary hospital in Ibadan, Nigeria: The RISK-HHF case-control study]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765799100330-8329fb7b-aa83-4fef-9ff0-b36be9d83786/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245734</link>
            <description><![CDATA[<div class="section" id="sec001"><h3 class="BHead" id="nov000-1">Aim</h3><p class="para" id="N65543">Hypertension is the leading cause of heart failure (HF) in sub-Saharan Africa. Preventive public health approach to reduce the scourge of HF must seek to understand the risk factors of HF in at-risk populations. The aim of this study was to characterize the risk factors of HF among patients with hypertension attending a cardiology clinic.</p></div><div class="section" id="sec002"><h3 class="BHead" id="nov000-2">Methods and results</h3><p class="para" id="N65549">One hundred and one (101) case-control age- and sex-matched pairs were recruited. The study population were adults with a clinical diagnosis of hypertensive HF (cases) and individuals with systemic hypertension without HF. They were interviewed and evaluated for cardiovascular risk factors. Associations between variables were tested with chi square test, Fisher’s exact test and independent sample t test as appropriate. Logistic regression modelling was used to determine the independent risk factors of hypertensive HF (HHF) in the study population while ‘punafcc’ package in stata12 was used to calculate the population attributable fraction (PAF) of the risk factors. Suboptimal medication adherence was the strongest adverse risk factor of HHF (medium adherence aOR: 3.53, 95%CI: 1.35–9.25; low adherence aOR: 9.44, 95%CI: 3.41–26.10) with a PAF of 67% followed by dipstick proteinuria (aOR: 4.22, 95%CI: 1.62–11.02; PAF: 34%) and alcohol consumption/day per 10grams (aOR: 1.23, 95%CI: 1.02–1.49; PAF: 22%). The protective risk factors of HHF were use of calcium channel blockers (aOR 0.25, 95%CI: 0.11–0.59; PAF: 59%), then daily fruits and vegetable consumption (aOR 0.41, 95%CI: 0.17–1.01; PAF: 46%), and eGFR (aOR 0.98, 95%CI: 0.96–0.99; PAF: 5.3%).</p></div><div class="section" id="sec003"><h3 class="BHead" id="nov000-3">Conclusions</h3><p class="para" id="N65555">The risk factors of HHF are amenable to lifestyle and dietary changes. Public health interventions and preventive cardiovascular care to improve medication adherence, promote fruit and vegetable consumption and reduce alcohol consumption among patients with hypertension are recommended. Renoprotection has utility in the prevention of HF among hypertensives.</p></div>]]></description>
            <pubDate><![CDATA[2021-01-25T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Devil in the details: Mechanistic variations impact information transfer across models of transcriptional cascades]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765793407307-a14956f2-f7a8-45d0-9f64-433bdca6e931/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245094</link>
            <description><![CDATA[<p class="para" id="N65539">The transcriptional network determines a cell’s internal state by regulating protein expression in response to changes in the local environment. Due to the interconnected nature of this network, information encoded in the abundance of various proteins will often propagate across chains of noisy intermediate signaling events. The data-processing inequality (DPI) leads us to expect that this intracellular game of “telephone” should degrade this type of signal, with longer chains losing successively more information to noise. However, a previous modeling effort predicted that because the steps of these signaling cascades do not truly represent independent stages of data processing, the limits of the DPI could seemingly be surpassed, and the amount of transmitted information could actually <i>increase</i> with chain length. What that work did not examine was whether this regime of growing information transmission was attainable by a signaling system constrained by the mechanistic details of more complex protein-binding kinetics. Here we address this knowledge gap through the lens of information theory by examining a model that explicitly accounts for the binding of each transcription factor to DNA. We analyze this model by comparing stochastic simulations of the fully nonlinear kinetics to simulations constrained by the linear response approximations that displayed a regime of growing information. Our simulations show that even when molecular binding is considered, there remains a regime wherein the transmitted information can grow with cascade length, but ends after a critical number of links determined by the kinetic parameter values. This inflection point marks where correlations decay in response to an oversaturation of binding sites, screening informative transcription factor fluctuations from further propagation down the chain where they eventually become indistinguishable from the surrounding levels of noise.</p>]]></description>
            <pubDate><![CDATA[2021-01-13T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Devil in the details: Mechanistic variations impact information transfer across models of transcriptional cascades]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765793407307-a14956f2-f7a8-45d0-9f64-433bdca6e931/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0245094</link>
            <description><![CDATA[<p class="para" id="N65539">The transcriptional network determines a cell’s internal state by regulating protein expression in response to changes in the local environment. Due to the interconnected nature of this network, information encoded in the abundance of various proteins will often propagate across chains of noisy intermediate signaling events. The data-processing inequality (DPI) leads us to expect that this intracellular game of “telephone” should degrade this type of signal, with longer chains losing successively more information to noise. However, a previous modeling effort predicted that because the steps of these signaling cascades do not truly represent independent stages of data processing, the limits of the DPI could seemingly be surpassed, and the amount of transmitted information could actually <i>increase</i> with chain length. What that work did not examine was whether this regime of growing information transmission was attainable by a signaling system constrained by the mechanistic details of more complex protein-binding kinetics. Here we address this knowledge gap through the lens of information theory by examining a model that explicitly accounts for the binding of each transcription factor to DNA. We analyze this model by comparing stochastic simulations of the fully nonlinear kinetics to simulations constrained by the linear response approximations that displayed a regime of growing information. Our simulations show that even when molecular binding is considered, there remains a regime wherein the transmitted information can grow with cascade length, but ends after a critical number of links determined by the kinetic parameter values. This inflection point marks where correlations decay in response to an oversaturation of binding sites, screening informative transcription factor fluctuations from further propagation down the chain where they eventually become indistinguishable from the surrounding levels of noise.</p>]]></description>
            <pubDate><![CDATA[2021-01-13T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Synaptic polarity and sign-balance prediction using gene expression data in the <i>Caenorhabditis elegans</i> chemical synapse neuronal connectome network]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765761881628-599fb090-a8a9-465e-992b-0c8b9b124ce4/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007974</link>
            <description><![CDATA[<p class="para" id="N65539">Graph theoretical analyses of nervous systems usually omit the aspect of connection polarity, due to data insufficiency. The chemical synapse network of <i>Caenorhabditis elegans</i> is a well-reconstructed directed network, but the signs of its connections are yet to be elucidated. Here, we present the gene expression-based sign prediction of the ionotropic chemical synapse connectome of <i>C</i>. <i>elegans</i> (3,638 connections and 20,589 synapses total), incorporating available presynaptic neurotransmitter and postsynaptic receptor gene expression data for three major neurotransmitter systems. We made predictions for more than two-thirds of these chemical synapses and observed an excitatory-inhibitory (E:I) ratio close to 4:1 which was found similar to that observed in many real-world networks. Our open source tool (http://EleganSign.linkgroup.hu) is simple but efficient in predicting polarities by integrating neuronal connectome and gene expression data.</p><p class="para" id="N65542">The fundamental way neurons communicate is by activating or inhibiting each other via synapses. The balance between the two is crucial for the optimal functioning of a nervous system. However, whole-brain synaptic polarity information is unavailable for any species and experimental validation is challenging. The roundworm <i>Caenorhabditis elegans</i> possesses a fully mapped connectome with an emerging gene expression profile of its 302 neurons. Based on the consideration that the polarity of a synapse can be determined by the neurotransmitter(s) expressed in the presynaptic neuron and the receptors expressed in the postsynaptic neuron, we conceptualized and created a tool that predicts synaptic polarities based on connectivity and gene expression information. Using currently available datasets we propose for the first time that the ratio of excitatory and inhibitory synapses in a partial connectome of <i>C</i>. <i>elegans</i> is around 4 to 1 which is in line with the balance observed in many natural systems. Our method opens a way to include spatial and temporal dynamics of synaptic polarity that would add a new dimension of plasticity in the excitatory:inhibitory balance. Our tool is freely available to be used on any network accompanied by any expression atlas.</p>]]></description>
            <pubDate><![CDATA[2020-12-21T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Synaptic polarity and sign-balance prediction using gene expression data in the <i>Caenorhabditis elegans</i> chemical synapse neuronal connectome network]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765761881628-599fb090-a8a9-465e-992b-0c8b9b124ce4/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1007974</link>
            <description><![CDATA[<p class="para" id="N65539">Graph theoretical analyses of nervous systems usually omit the aspect of connection polarity, due to data insufficiency. The chemical synapse network of <i>Caenorhabditis elegans</i> is a well-reconstructed directed network, but the signs of its connections are yet to be elucidated. Here, we present the gene expression-based sign prediction of the ionotropic chemical synapse connectome of <i>C</i>. <i>elegans</i> (3,638 connections and 20,589 synapses total), incorporating available presynaptic neurotransmitter and postsynaptic receptor gene expression data for three major neurotransmitter systems. We made predictions for more than two-thirds of these chemical synapses and observed an excitatory-inhibitory (E:I) ratio close to 4:1 which was found similar to that observed in many real-world networks. Our open source tool (http://EleganSign.linkgroup.hu) is simple but efficient in predicting polarities by integrating neuronal connectome and gene expression data.</p><p class="para" id="N65542">The fundamental way neurons communicate is by activating or inhibiting each other via synapses. The balance between the two is crucial for the optimal functioning of a nervous system. However, whole-brain synaptic polarity information is unavailable for any species and experimental validation is challenging. The roundworm <i>Caenorhabditis elegans</i> possesses a fully mapped connectome with an emerging gene expression profile of its 302 neurons. Based on the consideration that the polarity of a synapse can be determined by the neurotransmitter(s) expressed in the presynaptic neuron and the receptors expressed in the postsynaptic neuron, we conceptualized and created a tool that predicts synaptic polarities based on connectivity and gene expression information. Using currently available datasets we propose for the first time that the ratio of excitatory and inhibitory synapses in a partial connectome of <i>C</i>. <i>elegans</i> is around 4 to 1 which is in line with the balance observed in many natural systems. Our method opens a way to include spatial and temporal dynamics of synaptic polarity that would add a new dimension of plasticity in the excitatory:inhibitory balance. Our tool is freely available to be used on any network accompanied by any expression atlas.</p>]]></description>
            <pubDate><![CDATA[2020-12-21T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Optically activated, customizable, excitable cells]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765745510061-ca3c61f1-e287-4b86-9aac-4679d38d1198/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0229051</link>
            <description><![CDATA[<p class="para" id="N65539">Genetically encoded fluorescent biosensors are powerful tools for studying complex signaling in the nervous system, and now both Ca<sup>2+</sup> and voltage sensors are available to study the signaling behavior of entire neural circuits. There is a pressing need for improved sensors, but improving them is challenging because testing them involves a low throughput, labor-intensive processes. Our goal was to create synthetic, excitable cells that can be activated with brief pulses of blue light and serve as a medium throughput platform for screening the next generation of sensors. In this live cell system, blue light activates an adenylyl cyclase enzyme (bPAC) that increases intracellular cAMP (Stierl M et al. 2011). In turn, the cAMP opens a cAMP-gated ion channel. This produces slow, whole-cell Ca<sup>2+</sup> transients and voltage changes. To increase the speed of these transients, we add the inwardly rectifying potassium channel Kir2.1, the bacterial voltage-gated sodium channel NAVROSD, and Connexin-43. The result is a highly reproducible, medium-throughput, live cell system that can be used to screen voltage and Ca<sup>2+</sup> sensors.</p>]]></description>
            <pubDate><![CDATA[2020-12-30T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Optically activated, customizable, excitable cells]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765745510061-ca3c61f1-e287-4b86-9aac-4679d38d1198/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0229051</link>
            <description><![CDATA[<p class="para" id="N65539">Genetically encoded fluorescent biosensors are powerful tools for studying complex signaling in the nervous system, and now both Ca<sup>2+</sup> and voltage sensors are available to study the signaling behavior of entire neural circuits. There is a pressing need for improved sensors, but improving them is challenging because testing them involves a low throughput, labor-intensive processes. Our goal was to create synthetic, excitable cells that can be activated with brief pulses of blue light and serve as a medium throughput platform for screening the next generation of sensors. In this live cell system, blue light activates an adenylyl cyclase enzyme (bPAC) that increases intracellular cAMP (Stierl M et al. 2011). In turn, the cAMP opens a cAMP-gated ion channel. This produces slow, whole-cell Ca<sup>2+</sup> transients and voltage changes. To increase the speed of these transients, we add the inwardly rectifying potassium channel Kir2.1, the bacterial voltage-gated sodium channel NAVROSD, and Connexin-43. The result is a highly reproducible, medium-throughput, live cell system that can be used to screen voltage and Ca<sup>2+</sup> sensors.</p>]]></description>
            <pubDate><![CDATA[2020-12-30T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Proteostasis is adaptive: Balancing chaperone holdases against foldases]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765745065151-4513addc-3176-4b4d-9e1c-aa2323dc63dd/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008460</link>
            <description><![CDATA[<p class="para" id="N65539">Because a cell must adapt to different stresses and growth rates, its proteostasis system must too. How do cells detect and adjust proteome folding to different conditions? Here, we explore a biophysical cost-benefit principle, namely that the cell should keep its proteome as folded as possible at the minimum possible energy cost. This can be achieved by differential expression of chaperones–balancing foldases (which accelerate folding) against holdases (which act as parking spots). The model captures changes in the foldase-holdase ratio observed both <i>within</i> organisms during aging and <i>across</i> organisms of varying metabolic rates. This work describes a simple biophysical mechanism by which cellular proteostasis adapts to meet the needs of a changing growth environment.</p><p class="para" id="N65542">Cells must maintain low levels of protein unfolding to avoid deleterious outcomes such as protein aggregation, oxidative damage, or premature degradation. The proteins responsible for this, called chaperones, come in two main varieties: ATP-consuming “foldases” that help clients fold and ATP-independent “holdases” that hold unfolded proteins until a foldase arrives. While foldases are necessary for folding, they are expensive to have in high quantities. Given that chaperones are abundant and costly, cells are under strong selective pressure to find economical combinations of foldases and holdases for maintaining low levels of unfolded protein. Yet, it is presently unclear what the ideal combination is and how it varies with growth conditions. By examining a toy model of chaperone function and minimizing the total cost of folding at different rates of protein synthesis, we find that while foldases are necessary at fast growth, holdases become increasingly effective at slow growth. Unexpectedly, total chaperone requirements were predicted to increase as synthesis slows, consistent with observations across age in worms, as well as across species with varying metabolic rates. This work thus provides a general framework for understanding the chaperone requirements of a proteome in terms of an energy minimization principle.</p>]]></description>
            <pubDate><![CDATA[2020-12-14T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Proteostasis is adaptive: Balancing chaperone holdases against foldases]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765745065151-4513addc-3176-4b4d-9e1c-aa2323dc63dd/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008460</link>
            <description><![CDATA[<p class="para" id="N65539">Because a cell must adapt to different stresses and growth rates, its proteostasis system must too. How do cells detect and adjust proteome folding to different conditions? Here, we explore a biophysical cost-benefit principle, namely that the cell should keep its proteome as folded as possible at the minimum possible energy cost. This can be achieved by differential expression of chaperones–balancing foldases (which accelerate folding) against holdases (which act as parking spots). The model captures changes in the foldase-holdase ratio observed both <i>within</i> organisms during aging and <i>across</i> organisms of varying metabolic rates. This work describes a simple biophysical mechanism by which cellular proteostasis adapts to meet the needs of a changing growth environment.</p><p class="para" id="N65542">Cells must maintain low levels of protein unfolding to avoid deleterious outcomes such as protein aggregation, oxidative damage, or premature degradation. The proteins responsible for this, called chaperones, come in two main varieties: ATP-consuming “foldases” that help clients fold and ATP-independent “holdases” that hold unfolded proteins until a foldase arrives. While foldases are necessary for folding, they are expensive to have in high quantities. Given that chaperones are abundant and costly, cells are under strong selective pressure to find economical combinations of foldases and holdases for maintaining low levels of unfolded protein. Yet, it is presently unclear what the ideal combination is and how it varies with growth conditions. By examining a toy model of chaperone function and minimizing the total cost of folding at different rates of protein synthesis, we find that while foldases are necessary at fast growth, holdases become increasingly effective at slow growth. Unexpectedly, total chaperone requirements were predicted to increase as synthesis slows, consistent with observations across age in worms, as well as across species with varying metabolic rates. This work thus provides a general framework for understanding the chaperone requirements of a proteome in terms of an energy minimization principle.</p>]]></description>
            <pubDate><![CDATA[2020-12-14T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Proton-dependent inhibition, inverted voltage activation, and slow gating of CLC-0 Chloride Channel]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765740255542-538e9fd5-cf64-4f5c-82fb-42789cd65bb3/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0240704</link>
            <description><![CDATA[<p class="para" id="N65539">CLC-0, a prototype Cl<sup>-</sup> channel in the CLC family, employs two gating mechanisms that control its ion-permeation pore: fast gating and slow gating. The negatively-charged sidechain of a pore glutamate residue, E166, is known to be the fast gate, and the swinging of this sidechain opens or closes the pore of CLC-0 on the millisecond time scale. The other gating mechanism, slow gating, operates with much slower kinetics in the range of seconds to tens or even hundreds of seconds, and it is thought to involve still-unknown conformational rearrangements. Here, we find that low intracellular pH (pH<sub>i</sub>) facilitates the closure of the CLC-0’s slow gate, thus generating current inhibition. The rate of low pH<sub>i</sub>-induced current inhibition increases with intracellular H<sup>+</sup> concentration ([H<sup>+</sup>]<sub>i</sub>)—the time constants of current inhibition by low pH<sub>i</sub> = 4.5, 5.5 and 6 are roughly 0.1, 1 and 10 sec, respectively, at room temperature. In comparison, the time constant of the slow gate closure at pH<sub>i</sub> = 7.4 at room temperature is hundreds of seconds. The inhibition by low pH<sub>i</sub> is significantly less prominent in mutants favoring the slow-gate open state (such as C212S and Y512A), further supporting the fact that intracellular H<sup>+</sup> enhances the slow-gate closure in CLC-0. A fast inhibition by low pH<sub>i</sub> causes an apparent inverted voltage-dependent activation in the wild-type CLC-0, a behavior similar to those in some channel mutants such as V490W in which only membrane hyperpolarization can open the channel. Interestingly, when V490W mutation is constructed in the background of C212S or Y512A mutation, the inverted voltage-dependent activation disappears. We propose that the slow kinetics of CLC-0’s slow-gate closure may be due to low [H<sup>+</sup>]<sub>i</sub> rather than due to the proposed large conformational change of the channel protein. Our results also suggest that the inverted voltage-dependent opening observed in some mutant channels may result from fast closure of the slow gate by the mutations.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Proton-dependent inhibition, inverted voltage activation, and slow gating of CLC-0 Chloride Channel]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765740255542-538e9fd5-cf64-4f5c-82fb-42789cd65bb3/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0240704</link>
            <description><![CDATA[<p class="para" id="N65539">CLC-0, a prototype Cl<sup>-</sup> channel in the CLC family, employs two gating mechanisms that control its ion-permeation pore: fast gating and slow gating. The negatively-charged sidechain of a pore glutamate residue, E166, is known to be the fast gate, and the swinging of this sidechain opens or closes the pore of CLC-0 on the millisecond time scale. The other gating mechanism, slow gating, operates with much slower kinetics in the range of seconds to tens or even hundreds of seconds, and it is thought to involve still-unknown conformational rearrangements. Here, we find that low intracellular pH (pH<sub>i</sub>) facilitates the closure of the CLC-0’s slow gate, thus generating current inhibition. The rate of low pH<sub>i</sub>-induced current inhibition increases with intracellular H<sup>+</sup> concentration ([H<sup>+</sup>]<sub>i</sub>)—the time constants of current inhibition by low pH<sub>i</sub> = 4.5, 5.5 and 6 are roughly 0.1, 1 and 10 sec, respectively, at room temperature. In comparison, the time constant of the slow gate closure at pH<sub>i</sub> = 7.4 at room temperature is hundreds of seconds. The inhibition by low pH<sub>i</sub> is significantly less prominent in mutants favoring the slow-gate open state (such as C212S and Y512A), further supporting the fact that intracellular H<sup>+</sup> enhances the slow-gate closure in CLC-0. A fast inhibition by low pH<sub>i</sub> causes an apparent inverted voltage-dependent activation in the wild-type CLC-0, a behavior similar to those in some channel mutants such as V490W in which only membrane hyperpolarization can open the channel. Interestingly, when V490W mutation is constructed in the background of C212S or Y512A mutation, the inverted voltage-dependent activation disappears. We propose that the slow kinetics of CLC-0’s slow-gate closure may be due to low [H<sup>+</sup>]<sub>i</sub> rather than due to the proposed large conformational change of the channel protein. Our results also suggest that the inverted voltage-dependent opening observed in some mutant channels may result from fast closure of the slow gate by the mutations.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[A Neolithic mega-tsunami event in the eastern Mediterranean: Prehistoric settlement vulnerability along the Carmel coast, Israel]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765610413160-8befcd03-1bd9-45a5-9f2c-8dafbffb6dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0243619</link>
            <description><![CDATA[<p class="para" id="N65539">Tsunami events in antiquity had a profound influence on coastal societies. Six thousand years of historical records and geological data show that tsunamis are a common phenomenon affecting the eastern Mediterranean coastline. However, the possible impact of older tsunamis on prehistoric societies has not been investigated. Here we report, based on optically stimulated luminescence chronology, the earliest documented Holocene tsunami event, between 9.91 to 9.29 ka (kilo-annum), from the eastern Mediterranean at Dor, Israel. Tsunami debris from the early Neolithic is composed of marine sand embedded within fresh-brackish wetland deposits. Global and local sea-level curves for the period, 9.91–9.29 ka, as well as surface elevation reconstructions, show that the tsunami had a run-up of at least ~16 m and traveled between 3.5 to 1.5 km inland from the palaeo-coastline. Submerged slump scars on the continental slope, 16 km west of Dor, point to the nearby “Dor-complex” as a likely cause. The near absence of Pre-Pottery Neolithic A-B archaeological sites (11.70–9.80 cal. ka) suggest these sites were removed by the tsunami, whereas younger, late Pre-Pottery Neolithic B-C (9.25–8.35 cal. ka) and later Pottery-Neolithic sites (8.25–7.80 cal. ka) indicate resettlement following the event. The large run-up of this event highlights the disruptive impact of tsunamis on past societies along the Levantine coast.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[A Neolithic mega-tsunami event in the eastern Mediterranean: Prehistoric settlement vulnerability along the Carmel coast, Israel]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765610413160-8befcd03-1bd9-45a5-9f2c-8dafbffb6dc8/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0243619</link>
            <description><![CDATA[<p class="para" id="N65539">Tsunami events in antiquity had a profound influence on coastal societies. Six thousand years of historical records and geological data show that tsunamis are a common phenomenon affecting the eastern Mediterranean coastline. However, the possible impact of older tsunamis on prehistoric societies has not been investigated. Here we report, based on optically stimulated luminescence chronology, the earliest documented Holocene tsunami event, between 9.91 to 9.29 ka (kilo-annum), from the eastern Mediterranean at Dor, Israel. Tsunami debris from the early Neolithic is composed of marine sand embedded within fresh-brackish wetland deposits. Global and local sea-level curves for the period, 9.91–9.29 ka, as well as surface elevation reconstructions, show that the tsunami had a run-up of at least ~16 m and traveled between 3.5 to 1.5 km inland from the palaeo-coastline. Submerged slump scars on the continental slope, 16 km west of Dor, point to the nearby “Dor-complex” as a likely cause. The near absence of Pre-Pottery Neolithic A-B archaeological sites (11.70–9.80 cal. ka) suggest these sites were removed by the tsunami, whereas younger, late Pre-Pottery Neolithic B-C (9.25–8.35 cal. ka) and later Pottery-Neolithic sites (8.25–7.80 cal. ka) indicate resettlement following the event. The large run-up of this event highlights the disruptive impact of tsunamis on past societies along the Levantine coast.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Tracking collective cell motion by topological data analysis]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765608618460-8a7acce1-c0d9-4d07-a5f2-58f96d02579e/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008407</link>
            <description><![CDATA[<p class="para" id="N65539">By modifying and calibrating an active vertex model to experiments, we have simulated numerically a confluent cellular monolayer spreading on an empty space and the collision of two monolayers of different cells in an antagonistic migration assay. Cells are subject to inertial forces and to active forces that try to align their velocities with those of neighboring ones. In agreement with experiments in the literature, the spreading test exhibits formation of fingers in the moving interfaces, there appear swirls in the velocity field, and the polar order parameter and the correlation and swirl lengths increase with time. Numerical simulations show that cells inside the tissue have smaller area than those at the interface, which has been observed in recent experiments. In the antagonistic migration assay, a population of fluidlike Ras cells invades a population of wild type solidlike cells having shape parameters above and below the geometric critical value, respectively. Cell mixing or segregation depends on the junction tensions between different cells. We reproduce the experimentally observed antagonistic migration assays by assuming that a fraction of cells favor mixing, the others segregation, and that these cells are randomly distributed in space. To characterize and compare the structure of interfaces between cell types or of interfaces of spreading cellular monolayers in an automatic manner, we apply topological data analysis to experimental data and to results of our numerical simulations. We use time series of data generated by numerical simulations to automatically group, track and classify the advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These techniques of topological data analysis are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant for tissue engineering, biological effects of materials, tissue and organ regeneration.</p><p class="para" id="N65542">Confluent motion of cells in tissues plays a crucial role in wound healing, tissue repair, development, morphogenesis and in numerous pathological processes such as tumor invasion and metastatic cancer. For such complex processes, controlled experiments help clarifying the roles of chemical, mechanical and biological cues. Among them, spreading of cellular tissues on an empty space and antagonistic migration assays between cancerous and normal cells are quite revealing. The interfaces between confluent cellular aggregates uncover properties thereof when a combination of modeling, numerical simulation and data analysis is used. Here we have modified an active vertex model with a dynamics that includes inertia, friction and active forces that tend to align cells based on interaction with its immediate neighborhood. Selecting appropriately junction tensions among cells and using the SAMoS software, we have succeed in simulating assays of cellular tissue spreading on an empty space and the invasion of healthy tissue by cancerous one. We have introduced topological data analysis to characterize, track and compare in an automatic manner the interfaces of the tissue both in numerical simulations and from experimental data of normal and Ras modified precancerous Human Embryonic Kidney cells. We find good agreement when normal cells are solidlike and modified cells are liquidlike according to their shape parameters. In addition, cell variability means that a fraction of randomly distributed cells favor mixing, the others segregation. Topological data analysis techniques are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant in ascertaining how the biophysical features of materials may affect tissue and organ regeneration.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Tracking collective cell motion by topological data analysis]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765608618460-8a7acce1-c0d9-4d07-a5f2-58f96d02579e/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pcbi.1008407</link>
            <description><![CDATA[<p class="para" id="N65539">By modifying and calibrating an active vertex model to experiments, we have simulated numerically a confluent cellular monolayer spreading on an empty space and the collision of two monolayers of different cells in an antagonistic migration assay. Cells are subject to inertial forces and to active forces that try to align their velocities with those of neighboring ones. In agreement with experiments in the literature, the spreading test exhibits formation of fingers in the moving interfaces, there appear swirls in the velocity field, and the polar order parameter and the correlation and swirl lengths increase with time. Numerical simulations show that cells inside the tissue have smaller area than those at the interface, which has been observed in recent experiments. In the antagonistic migration assay, a population of fluidlike Ras cells invades a population of wild type solidlike cells having shape parameters above and below the geometric critical value, respectively. Cell mixing or segregation depends on the junction tensions between different cells. We reproduce the experimentally observed antagonistic migration assays by assuming that a fraction of cells favor mixing, the others segregation, and that these cells are randomly distributed in space. To characterize and compare the structure of interfaces between cell types or of interfaces of spreading cellular monolayers in an automatic manner, we apply topological data analysis to experimental data and to results of our numerical simulations. We use time series of data generated by numerical simulations to automatically group, track and classify the advancing interfaces of cellular aggregates by means of bottleneck or Wasserstein distances of persistent homologies. These techniques of topological data analysis are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant for tissue engineering, biological effects of materials, tissue and organ regeneration.</p><p class="para" id="N65542">Confluent motion of cells in tissues plays a crucial role in wound healing, tissue repair, development, morphogenesis and in numerous pathological processes such as tumor invasion and metastatic cancer. For such complex processes, controlled experiments help clarifying the roles of chemical, mechanical and biological cues. Among them, spreading of cellular tissues on an empty space and antagonistic migration assays between cancerous and normal cells are quite revealing. The interfaces between confluent cellular aggregates uncover properties thereof when a combination of modeling, numerical simulation and data analysis is used. Here we have modified an active vertex model with a dynamics that includes inertia, friction and active forces that tend to align cells based on interaction with its immediate neighborhood. Selecting appropriately junction tensions among cells and using the SAMoS software, we have succeed in simulating assays of cellular tissue spreading on an empty space and the invasion of healthy tissue by cancerous one. We have introduced topological data analysis to characterize, track and compare in an automatic manner the interfaces of the tissue both in numerical simulations and from experimental data of normal and Ras modified precancerous Human Embryonic Kidney cells. We find good agreement when normal cells are solidlike and modified cells are liquidlike according to their shape parameters. In addition, cell variability means that a fraction of randomly distributed cells favor mixing, the others segregation. Topological data analysis techniques are scalable and could be used in studies involving large amounts of data. Besides applications to wound healing and metastatic cancer, these studies are relevant in ascertaining how the biophysical features of materials may affect tissue and organ regeneration.</p>]]></description>
            <pubDate><![CDATA[2020-12-23T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Decreased cardiac pacemaking and attenuated β-adrenergic response in TRIC-A knockout mice]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765602034859-6279c64d-9a3e-49f9-8b30-dc65a63bdf1a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0244254</link>
            <description><![CDATA[<p class="para" id="N65539">Changes in intracellular calcium levels in the sinus node modulate cardiac pacemaking (the calcium clock). Trimeric intracellular cation (TRIC) channels are counterion channels on the surface of the sarcoplasmic reticulum and compensate for calcium release from ryanodine receptors, which play a major role in calcium-induced calcium release (CICR) and the calcium clock. TRIC channels are expected to affect the calcium clock in the sinus node. However, their physiological importance in cardiac rhythm formation remains unclear. We evaluated the importance of TRIC channels on cardiac pacemaking using TRIC-A-null (TRIC-A<sup>–/–</sup>) as well as TRIC-B<sup>+/–</sup>mice. Although systolic blood pressure (SBP) was not significantly different between wild-type (WT), TRIC-B<sup>+/–</sup>, and TRIC-A<sup>–/–</sup>mice, heart rate (HR) was significantly lower in TRIC-A<sup>–/–</sup>mice than other lines. Interestingly, HR and SBP showed a positive correlation in WT and TRIC-B<sup>+/–</sup>mice, while no such correlation was observed in TRIC-A<sup>–/–</sup>mice, suggesting modification of the blood pressure regulatory system in these mice. Isoproterenol (0.3 mg/kg) increased the HR in WT mice (98.8 ± 15.1 bpm), whereas a decreased response in HR was observed in TRIC-A<sup>–/–</sup>mice (23.8 ± 5.8 bpm), suggesting decreased sympathetic responses in TRIC-A<sup>–/–</sup>mice. Electrocardiography revealed unstable R-R intervals in TRIC-A<sup>–/–</sup>mice. Furthermore, TRIC-A<sup>–/–</sup>mice sometimes showed sinus pauses, suggesting a significant role of TRIC-A channels in cardiac pacemaking. In isolated atrium contraction or action potential recording, TRIC-A<sup>–/–</sup>mice showed decreased response to a β-adrenergic sympathetic nerve agonist (isoproterenol, 100 nM), indicating decreased sympathetic responses. In summary, TRIC-A<sup>–/–</sup>mice showed decreased cardiac pacemaking in the sinus node and attenuated responses to β-adrenergic stimulation, indicating the involvement of TRIC-A channels in cardiac rhythm formation and decreased sympathetic responses.</p>]]></description>
            <pubDate><![CDATA[2020-12-21T00:00]]></pubDate>
        </item><item>
            <title><![CDATA[Decreased cardiac pacemaking and attenuated β-adrenergic response in TRIC-A knockout mice]]></title>
            <media:thumbnail url="https://storage.googleapis.com/nova-demo-unsecured-files/unsecured/content-1765602034859-6279c64d-9a3e-49f9-8b30-dc65a63bdf1a/cover.png"></media:thumbnail>
            <link>https://www.novareader.co/book/isbn/10.1371/journal.pone.0244254</link>
            <description><![CDATA[<p class="para" id="N65539">Changes in intracellular calcium levels in the sinus node modulate cardiac pacemaking (the calcium clock). Trimeric intracellular cation (TRIC) channels are counterion channels on the surface of the sarcoplasmic reticulum and compensate for calcium release from ryanodine receptors, which play a major role in calcium-induced calcium release (CICR) and the calcium clock. TRIC channels are expected to affect the calcium clock in the sinus node. However, their physiological importance in cardiac rhythm formation remains unclear. We evaluated the importance of TRIC channels on cardiac pacemaking using TRIC-A-null (TRIC-A<sup>–/–</sup>) as well as TRIC-B<sup>+/–</sup>mice. Although systolic blood pressure (SBP) was not significantly different between wild-type (WT), TRIC-B<sup>+/–</sup>, and TRIC-A<sup>–/–</sup>mice, heart rate (HR) was significantly lower in TRIC-A<sup>–/–</sup>mice than other lines. Interestingly, HR and SBP showed a positive correlation in WT and TRIC-B<sup>+/–</sup>mice, while no such correlation was observed in TRIC-A<sup>–/–</sup>mice, suggesting modification of the blood pressure regulatory system in these mice. Isoproterenol (0.3 mg/kg) increased the HR in WT mice (98.8 ± 15.1 bpm), whereas a decreased response in HR was observed in TRIC-A<sup>–/–</sup>mice (23.8 ± 5.8 bpm), suggesting decreased sympathetic responses in TRIC-A<sup>–/–</sup>mice. Electrocardiography revealed unstable R-R intervals in TRIC-A<sup>–/–</sup>mice. Furthermore, TRIC-A<sup>–/–</sup>mice sometimes showed sinus pauses, suggesting a significant role of TRIC-A channels in cardiac pacemaking. In isolated atrium contraction or action potential recording, TRIC-A<sup>–/–</sup>mice showed decreased response to a β-adrenergic sympathetic nerve agonist (isoproterenol, 100 nM), indicating decreased sympathetic responses. In summary, TRIC-A<sup>–/–</sup>mice showed decreased cardiac pacemaking in the sinus node and attenuated responses to β-adrenergic stimulation, indicating the involvement of TRIC-A channels in cardiac rhythm formation and decreased sympathetic responses.</p>]]></description>
            <pubDate><![CDATA[2020-12-21T00:00]]></pubDate>
        </item>
    </channel>
</rss>