From the Cover: Solar anti-icing surface with enhanced condensate self-removing at extreme environmental conditions

The inhibition of condensation freezing under extreme conditions (i.e., ultra-low temperature and high humidity) remains a daunting challenge in the field of anti-icing. As water vapor easily condensates or desublimates and melted water refreezes instantly, these cause significant performance decrease of most anti-icing surfaces at such extreme conditions. Herein, inspired by wheat leaves, an effective condensate self-removing solar anti-icing/frosting surface (CR-SAS) is fabricated using ultrafast pulsed laser deposition technology, which exhibits synergistic effects of enhanced condensate self-removal and efficient solar anti-icing. The superblack CR-SAS displays superior anti-reflection and photothermal conversion performance, benefiting from the light trapping effect in the micro/nano hierarchical structures and the thermoplasmonic effect of the iron oxide nanoparticles. Meanwhile, the CR-SAS displays superhydrophobicity to condensed water, which can be instantly shed off from the surface before freezing through self-propelled droplet jumping, thus leading to a continuously refreshed dry area available for sunlight absorption and photothermal conversion. Under one-sun illumination, the CR-SAS can be maintained ice free even under an ambient environment of −50 °C ultra-low temperature and extremely high humidity (ice supersaturation degree of ∼260). The excellent environmental versatility, mechanical durability, and material adaptability make CR-SAS a promising anti-icing candidate for broad practical applications even in harsh environments.

Condensation freezing/frosting on solid surfaces causes severe economic and safety issues. Thus, highly efficient anti-icing/frosting approaches are vital in many engineering applications, ranging from air conditioners to power transmission systems (1–4). Tremendous efforts have been invested into designing active and passive anti-icing surfaces. Active anti-icing strategies including mechanical, chemical, and thermal methods often consume high amounts of energy and require specific facilities for deicing, which limit their practical applications (5, 6). Passive anti-icing surfaces involve strategies to delay and inhibit ice formation, such as hydrated surfaces, lubricant-infused surfaces, bioinspired anti-freezing surfaces, and superhydrophobic surfaces (SHSs) (3, 7–13). Although these passive icephobic materials offer numerous advantages to prevent ice accretion, each comes with its own limitations (14).At extreme environmental conditions where condensation and desublimation are strongly promoted, an optimal passive anti-icing surface should immediately remove condensed water droplets and leave no water for freezing. To prepare such a kind of surface, lots of efforts have been made to study the abilities of SHSs for removing condensed water. However, regular SHSs are incapable of removing microscale condensed water droplets under humid conditions due to the wetting of surface micronanostructures, which leads to the formation of highly pinned Wenzel droplets and the loss of superhydrophobicity (15–17). Tremendous investigations have shown that SHSs with specially designed structures can retain superhydrophobicity to condensed water droplets, and coalesced droplets can be spontaneously removed from the surface via self-propelled jumping (18–20). The self-propelled jumping of condensed droplets is driven by the released surface energy during droplet coalescence after overcoming solid–liquid adhesion (21–25). However, these surfaces inevitably lose their water repellency at low temperatures (i.e., < −15 °C) because of freezing (9, 21, 26). Thus, it is highly desirable to design new anti-icing surfaces that can maintain capability for removing the condensed water at extreme environmental conditions, which is highly important for anti-icing applications in many scenarios (i.e., aircrafts flying through clouds, wind turbines operating in winter, and power transmission facilities working in extremely cold and humid regions) (27).Recently, intensive research efforts have been dedicated to solar anti-icing/frosting surfaces (SASs), which can absorb sunlight efficiently and convert solar energy to heat, thereby delaying or preventing ice formation (28–30). Because of its effective utilization of clean and renewable solar energy, SASs are environmentally friendly and energy efficient. Notably, a number of photothermal conversion materials including carbon materials, conjugated polymers, two-dimensional nanostructural materials, and metallic particles have also been applied as SASs (28, 31–34). Although remarkable improvements have been made, some challenges have yet to be tackled. For instance, most of the reported SASs cannot remove the condensed water effectively, particularly in cold and humid conditions (35, 36). The accumulation of the condensed water would significantly enhance reflectance leading to reduced photothermal efficiency (37) and decreased temperature. As a result, frost formation through the freezing of condensed water (condensation frosting) will prevent sunlight from reaching the SAS, resulting in the complete loss of its photothermal capability: as light harvesting becomes less efficient, the temperature of the SAS decreases, resulting in a vicious compounding cycle of condensation freezing. Moreover, the contaminants on the SAS can also inhibit sunlight absorption (38). Therefore, it is desirable to develop highly photothermal–efficient SASs with the abilities of removing condensed water to maintain high temperature and self-cleaning for avoiding blockage of sunlight by contaminants.Herein, we present a proof-of-concept SAS with synergistically binary effects of enhanced condensate self-removing and efficient solar anti-icing. We fabricated hierarchically structured materials using ultrafast pulsed laser deposition (PLD) technology. Low-effective refractive index and multilayered iron oxide structures endow the material with broadband ultralow reflectance and high solar-to-heat conversation efficiency. The hierarchically structured SHS demonstrated the capability of removing condensates before freezing via self-propelled droplet jumping induced by droplet coalescence and evaporation flux under heating. The sustained shedding of condensed water droplets resulted in a continuously refreshed, dry and clean area available for efficient sunlight absorption and photothermal conversion; the temperature of the condensate self-removing solar anti-icing/frosting surface (CR-SAS) could be constantly maintained above the freezing temperature, which in turn ensured high-performance water repellency. With such synergistic mutual-benefitting effects of condensate self-removal and photothermal heating, under one-sun illumination, the CR-SAS remained ice-free even at an ambient temperature (T_a) of −50 °C and ultra-high humidity with a supersaturation degree (SSD) of ∼260, demonstrating superior anti-icing performances under extremely harsh operating conditions.     

Chemokine CCL5 promotes robust optic nerve regeneration and mediates many of the effects of CNTF gene therapy

Ciliary neurotrophic factor (CNTF) is a leading therapeutic candidate for several ocular diseases and induces optic nerve regeneration in animal models. Paradoxically, however, although CNTF gene therapy promotes extensive regeneration, recombinant CNTF (rCNTF) has little effect. Because intraocular viral vectors induce inflammation, and because CNTF is an immune modulator, we investigated whether CNTF gene therapy acts indirectly through other immune mediators. The beneficial effects of CNTF gene therapy remained unchanged after deleting CNTF receptor alpha (CNTFRα) in retinal ganglion cells (RGCs), the projection neurons of the retina, but were diminished by depleting neutrophils or by genetically suppressing monocyte infiltration. CNTF gene therapy increased expression of C-C motif chemokine ligand 5 (CCL5) in immune cells and retinal glia, and recombinant CCL5 induced extensive axon regeneration. Conversely, CRISPR-mediated knockdown of the cognate receptor (CCR5) in RGCs or treating wild-type mice with a CCR5 antagonist repressed the effects of CNTF gene therapy. Thus, CCL5 is a previously unrecognized, potent activator of optic nerve regeneration and mediates many of the effects of CNTF gene therapy.

Like most pathways in the mature central nervous system (CNS), the optic nerve cannot regenerate once damaged due in part to cell-extrinsic suppressors of axon growth (1, 2) and the low intrinsic growth capacity of adult retinal ganglion cells (RGCs), the projection neurons of the eye (3–5). Consequently, traumatic or ischemic optic nerve injury or degenerative diseases such as glaucoma lead to irreversible visual losses. Experimentally, some degree of regeneration can be induced by intraocular inflammation or growth factors expressed by inflammatory cells (6–10), altering the cell-intrinsic growth potential of RGCs (3–5), enhancing physiological activity (11, 12), chelating free zinc (13, 14), and other manipulations (15–19). However, the extent of regeneration achieved to date remains modest, underlining the need for more effective therapies.Ciliary neurotrophic factor (CNTF) is a leading therapeutic candidate for glaucoma and other ocular diseases (20–23). Activation of the downstream signal transduction cascade requires CNTF to bind to CNTF receptor-α (CNTFRα) (24), which leads to recruitment of glycoprotein 130 (gp130) and leukemia inhibitory factor receptor-β (LIFRβ) to form a tripartite receptor complex (25). CNTFRα anchors to the plasma membrane through a glycosylphosphatidylinositol linkage (26) and can be released and become soluble through phospholipase C-mediated cleavage (27). CNTF has been reported to activate STAT3 phosphorylation in retinal neurons, including RGCs, and to promote survival, but it is unknown whether these effects are mediated by direct action of CNTF on RGCs via CNTFRα (28). Our previous studies showed that CNTF promotes axon outgrowth from neonate RGCs in culture (29) but fails to do so in cultured mature RGCs (8) or in vivo (6). Although some studies report that recombinant CNTF (rCNTF) can promote optic nerve regeneration (20, 30, 31), others find little or no effect unless SOCS3 (suppressor of cytokine signaling-3), an inhibitor of the Jak-STAT pathway, is deleted in RGCs (5, 6, 32). In contrast, multiple studies show that adeno-associated virus (AAV)-mediated expression of CNTF in RGCs induces strong regeneration (33–40). The basis for the discrepant effects of rCNTF and CNTF gene therapy is unknown but is of considerable interest in view of the many promising clinical and preclinical outcomes obtained with CNTF to date.Because intravitreal virus injections induce inflammation (41), we investigated the possibility that CNTF, a known immune modulator (42–44), might act by elevating expression of other immune-derived factors. We report here that the beneficial effects of CNTF gene therapy in fact require immune system activation and elevation of C-C motif chemokine ligand 5 (CCL5). Depletion of neutrophils, global knockout (KO) or RGC-selective deletion of the CCL5 receptor CCR5, or a CCR5 antagonist all suppress the effects of CNTF gene therapy, whereas recombinant CCL5 (rCCL5) promotes axon regeneration and increases RGC survival. These studies point to CCL5 as a potent monotherapy for optic nerve regeneration and to the possibility that other applications of CNTF and other forms of gene therapy might similarly act indirectly through other factors.     

Cholinergic suppression of hippocampal sharp-wave ripples impairs working memory

Learning and memory are assumed to be supported by mechanisms that involve cholinergic transmission and hippocampal theta. Using G protein–coupled receptor-activation–based acetylcholine sensor (GRAB_ACh3.0) with a fiber-photometric fluorescence readout in mice, we found that cholinergic signaling in the hippocampus increased in parallel with theta/gamma power during walking and REM sleep, while ACh3.0 signal reached a minimum during hippocampal sharp-wave ripples (SPW-R). Unexpectedly, memory performance was impaired in a hippocampus-dependent spontaneous alternation task by selective optogenetic stimulation of medial septal cholinergic neurons when the stimulation was applied in the delay area but not in the central (choice) arm of the maze. Parallel with the decreased performance, optogenetic stimulation decreased the incidence of SPW-Rs. These findings suggest that septo–hippocampal interactions play a task-phase–dependent dual role in the maintenance of memory performance, including not only theta mechanisms but also SPW-Rs.

The neurotransmitter acetylcholine is thought to be critical for hippocampus-dependent declarative memories (1, 2). Reduction in cholinergic neurotransmission, either in Alzheimer’s disease or in experiments with cholinergic antagonists, such as scopolamine, impairs memory function (3–8). Acetylcholine may bring about its beneficial effects on memory encoding by enhancing theta rhythm oscillations, decreasing recurrent excitation, and increasing synaptic plasticity (9–11). Conversely, drugs which activate cholinergic receptors enhance learning and, therefore, are a neuropharmacological target for the treatment of memory deficits in Alzheimer’s disease (5, 12, 13).The contribution of cholinergic mechanisms in the acquisition of long-term memories and the role of the hippocampal–entorhinal–cortical interactions are well supported by experimental data (5, 12, 13). In addition, working memory or “short-term” memory is also supported by the hippocampal–entorhinal–prefrontal cortex (14–16). Working memory in humans is postulated to be a conscious process to “keep things in mind” transiently (16). In rodents, matching to sample task, spontaneous alternation between reward locations, and the radial maze task have been suggested to function as a homolog of working memory [“working memory like” (17)].Cholinergic activity is a critical requirement for working memory (18, 19) and for sustaining theta oscillations (10, 20–22). In support of this contention, theta–gamma coupling and gamma power are significantly higher in the choice arm of the maze, compared with those in the side arms where working memory is no longer needed for correct performance (23–26). It has long been hypothesized that working memory is maintained by persistent firing of neurons, which keep the presented items in a transient store in the prefrontal cortex and hippocampal–entorhinal system (27–31), although the exact mechanisms are debated (32–37). An alternative hypothesis holds that items of working memory are stored in theta-nested gamma cycles (38). Common in these models of working memory is the need for an active, cholinergic system–dependent mechanism (39–41). However, in spontaneous alternation tasks, the animals are not moving continuously during the delay, and theta oscillations are not sustained either. During the immobility epochs, theta is replaced by intermittent sharp-wave ripples (SPW-R), yet memory performance does not deteriorate. On the contrary, artificial blockade of SPW-Rs can impair memory performance (42, 43), and prolongation of SPW-Rs improves performance (44). Under the cholinergic hypothesis of working memory, such a result is unexpected.To address the relationship between cholinergic/theta versus SPW-R mechanism in spontaneous alternation, we used a G protein–coupled receptor-activation–based acetylcholine sensor (GRAB_ACh3.0) (45) to monitor acetylcholine (ACh) activity during memory performance in mice. In addition, we optogenetically enhanced cholinergic tone, which suppresses SPW-Rs by a different mechanism than electrically or optogenetically induced silencing of neurons in the hippocampus (43, 44). We show that cholinergic signaling in the hippocampus increases in parallel with theta power/score during walking and rapid eye movement (REM) sleep and reaches a transient minimum during SPW-Rs. Selective optogenetic stimulation of medial septal cholinergic neurons decreased the incidence of SPW-Rs during non-REM sleep (46–48), as well as during the delay epoch of a working memory task and impaired memory performance. These findings demonstrate that memory performance is supported by complementary theta and SPW-R mechanisms.     

Tropospheric carbonyl sulfide mass balance based on direct measurements of sulfur isotopes

Robust estimates for the rates and trends in terrestrial gross primary production (GPP; plant CO₂ uptake) are needed. Carbonyl sulfide (COS) is the major long-lived sulfur-bearing gas in the atmosphere and a promising proxy for GPP. Large uncertainties in estimating the relative magnitude of the COS sources and sinks limit this approach. Sulfur isotope measurements (³⁴S/³²S; δ³⁴S) have been suggested as a useful tool to constrain COS sources. Yet such measurements are currently scarce for the atmosphere and absent for the marine source and the plant sink, which are two main fluxes. Here we present sulfur isotopes measurements of marine and atmospheric COS, and of plant-uptake fractionation experiments. These measurements resulted in a complete data-based tropospheric COS isotopic mass balance, which allows improved partition of the sources. We found an isotopic (δ³⁴S ± SE) value of 13.9 ± 0.1‰ for the troposphere, with an isotopic seasonal cycle driven by plant uptake. This seasonality agrees with a fractionation of −1.9 ± 0.3‰ which we measured in plant-chamber experiments. Air samples with strong anthropogenic influence indicated an anthropogenic COS isotopic value of 8 ± 1‰. Samples of seawater-equilibrated-air indicate that the marine COS source has an isotopic value of 14.7 ± 1‰. Using our data-based mass balance, we constrained the relative contribution of the two main tropospheric COS sources resulting in 40 ± 17% for the anthropogenic source and 60 ± 20% for the oceanic source. This constraint is important for a better understanding of the global COS budget and its improved use for GPP determination.

The Earth system is going through rapid changes as the climate warms and CO₂ level rises. This rise in CO₂ is mitigated by plant uptake; hence, it is important to estimate global and regional photosynthesis rates and trends (1). Yet, robust tools for investigating these processes at a large scale are scarce (2). Recent studies suggest that carbonyl sulfide (COS) could provide an improved constraint on terrestrial photosynthesis (gross primary production, GPP) (2–12). COS is the major long-lived sulfur-bearing gas in the atmosphere and the main supplier of sulfur to the stratospheric sulfate aerosol layer (13), which exerts a cooling effect on the Earth’s surface and regulates stratospheric ozone chemistry (14).During terrestrial photosynthesis, COS diffuses into leaf stomata and is consumed by photosynthetic enzymes in a similar manner to CO₂ (3–5). Contrary to CO₂, COS undergoes rapid and irreversible hydrolysis mainly by the enzyme carbonic-anhydrase (6, 7). Thus, COS can be used as a proxy for the one-way flux of CO₂ removal from the atmosphere by terrestrial photosynthesis (2, 8–11). However, the large uncertainties in estimating the COS sources weaken this approach (10–12, 15). Tropospheric COS has two main sources: the oceans and anthropogenic emissions, and one main sink–terrestrial plant uptake (8, 10–13). Smaller sources include biomass burning, soil emissions, wetlands, volcanoes, and smaller sinks include OH destruction, stratospheric destruction, and soil uptake (12). The largest source of COS to the atmosphere is the ocean, both as direct COS emission, and as indirect carbon disulfide (CS₂) and dimethylsulfide (DMS) emissions that are rapidly oxidized to COS (10, 16–20). Recent studies suggest oceanic COS emissions are in the range of 200–4,000 GgS/y (19–22). The second major COS source is the anthropogenic source, which is dominated by indirect emissions derived from CS₂ oxidation, mainly from the use of CS₂ as an industrial solvent. Direct emissions of COS are mainly derived from coal and fuel combustion (17, 23, 24). Recent studies suggest that anthropogenic emissions are in the range of 150–585 GgS/y (23, 24). The terrestrial plant uptake is estimated to be in the range of 400–1,360 GgS/y (11). Measurements of sulfur isotope ratios (δ³⁴S) in COS may be used to track COS sources and thus reduce the uncertainties in their flux estimations (15, 25–27). However, the isotopic mass balance approach works best if the COS end members are directly measured and have a significantly different isotopic signature. Previous δ³⁴S measurements of atmospheric COS are scarce and there have been no direct measurements of two important components: the δ³⁴S of oceanic COS emissions, and the isotopic fractionation of COS during plant uptake (15, 25–27). In contrast to previous studies that used assessments for these isotopic values, our aim was to directly measure the isotopic values of these missing components, and to determine the tropospheric COS δ³⁴S variability over space and time.     

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae (Effrenium voratum), which swim in circles at the air–liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Active matter exists over a wide range of spatial and temporal scales (1–6) from animal groups (7, 8) to robot swarms (9–11), to cell colonies and tissues (12–16), to cytoskeletal extracts (17–20), to man-made microswimmers (21–25). Constituent particles in active matter systems are driven out of thermal equilibrium at the individual level; they interact to develop a wealth of intriguing collective phenomena, including clustering (13, 22, 24), flocking (11, 26), swarming (12, 13), spontaneous flow (14, 20), and giant density fluctuations (10, 11). Many of these observed phenomena have been successfully described by particle-based or continuum models (1–6), which highlight the important roles of both individual motility and interparticle interactions in determining system dynamics.Current active matter research focuses primarily on linearly swimming particles which have a symmetric body and self-propel along one of the symmetry axes. However, a perfect alignment between the propulsion direction and body axis is rarely found in reality. Deviation from such a perfect alignment leads to a persistent curvature in the microswimmer trajectories; examples of such circle microswimmers include anisotropic artificial micromotors (27, 28), self-propelled nematic droplets (29, 30), magnetotactic bacteria and Janus particles in rotating external fields (31, 32), Janus particle in viscoelastic medium (33), and sperm and bacteria near interfaces (34, 35). Chiral motility of circle microswimmers, as predicted by theoretical and numerical investigations, can lead to a range of interesting collective phenomena in circular microswimmers, including vortex structures (36, 37), localization in traps (38), enhanced flocking (39), and hyperuniform states (40). However, experimental verifications of these predictions are limited (32, 35), a situation mainly due to the scarcity of suitable experimental systems.Here we address this challenge by investigating marine algae Effrenium voratum (41, 42). At air–liquid interfaces, E. voratum cells swim in circles via two eukaryotic flagella: a transverse flagellum encircling the cellular anteroposterior axis and a longitudinal one running posteriorly. Over a wide range of densities, circling E. voratum cells self-organize into disordered hyperuniform states with suppressed density fluctuations at large length scales. Hyperuniformity (43, 44) has been considered as a new form of material order which leads to novel functionalities (45–49); it has been observed in many systems, including avian photoreceptor patterns (50), amorphous ices (51), amorphous silica (52), ultracold atoms (53), soft matter systems (54–61), and stochastic models (62–64). Our work demonstrates the existence of hyperuniformity in active matter and shows that hydrodynamic interactions can be used to construct hyperuniform states.     

CXCR4/CXCL12-mediated entrapment of axons at the injury site compromises optic nerve regeneration

Regenerative failure in the mammalian optic nerve is generally attributed to axotomy-induced retinal ganglion cell (RGC) death, an insufficient intrinsic regenerative capacity, and an extrinsic inhibitory environment. Here, we show that a chemoattractive CXCL12/CXCR4-dependent mechanism prevents the extension of growth-stimulated axons into the distal nerve. The chemokine CXCL12 is chemoattractive toward axonal growth cones in an inhibitory environment, and these effects are entirely abolished by the specific knockout of its receptor, CXCR4 (CXCR4^−/−), in cultured regenerating RGCs. Notably, 8% of naïve RGCs express CXCL12 and transport the chemokine along their axons in the nerve. Thus, axotomy causes its release at the injury site. However, most osteopontin-positive α-RGCs, the main neuronal population that survives optic nerve injury, express CXCR4 instead. Thus, CXCL12-mediated attraction prevents growth-stimulated axons from regenerating distally in the nerve, indicated by axons returning to the lesion site. Accordingly, specific depletion of CXCR4 in RGC reduces aberrant axonal growth and enables long-distance regeneration. Likewise, CXCL12 knockout in RGCs fully mimics these CXCR4^−/− effects. Thus, active CXCL12/CXCR4-mediated entrapment of regenerating axons to the injury site contributes to regenerative failure in the optic nerve.

Retinal ganglion cells (RGCs) convey the visual input from the eye through the optic nerve and optic tract into the brain’s target regions. As typical neurons of the central nervous system (CNS), mammalian RGCs lose most of their capability to regrow injured axons after birth (1, 2), leading to an irreversible functional loss after optic nerve damage. To date, regenerative failure has been mainly attributed to three leading causes: 1) axotomy-induced apoptosis of RGCs, 2) the low intrinsic capacity to regrow axons, and 3) the external inhibitory environment with CNS myelin and glial scar proteins (3, 4).One widely used approach to delay axotomy-induced RGC degeneration and activate the intrinsic regenerative capacity of injured axons is inflammatory stimulation (IS) in the eye induced by a lens injury, intravitreal Pam₃Cys, or zymosan injection (5–7). IS leads to the expression and release of CNTF, LIF, and IL-6 from retinal astrocytes and Müller cells (8–10), which directly interact with RGCs and activate neuroprotective/regenerative signaling such as the JAK/STAT3 pathway (8, 9, 11, 12). IS, therefore, enables moderate axon regeneration beyond the lesion site of the optic nerve. Although combinatorial strategies, together with measures overcoming the inhibitory CNS environment synergistically, further improve IS-mediated optic nerve regeneration (13–17), the overall outcome remains mostly unsatisfactory. Thus, additional unknown mechanisms besides neurodegeneration, low intrinsic capacity, and the inhibitory environment might contribute to optic nerve regeneration failure.The chemokine receptor CXCR4, a seven-transmembrane G protein–coupled receptor, is expressed in embryonic and adult neurons (18–20). We have recently shown that this receptor is also expressed in the somata and axons of adult rat RGCs (18). Next to its role as a coreceptor for HIV entry and cancer-cell migration/proliferation (21, 22), CXCR4 is reportedly involved in neurogenesis and axonal pathfinding during the embryonal development of RGCs (20, 23, 24). CXCR4 regulates different signaling pathways upon binding its ligand CXCL12 (also known as stromal cell–derived factor 1, SDF-1), which is part of the chemokine family of chemotactic cytokines in the immune system involved in the attraction of lymphocytes (25, 26). CXCL12 is also reportedly expressed by some CNS neurons, astrocytes, and microglia (19, 27–30). As the CXCR4/CXCL12 axis is highly conserved between different species (31) and involved in axonal pathfinding during embryonal development of RGCs (20, 32), we speculated that CXCR4 expression in adult RGCs might also play a role in the regenerative processes of mature axons.The current study shows that growth-stimulated axons of RGCs are actively attracted and entrapped at the lesion site of the optic nerve by a CXCL12/CXCR4-dependent mechanism. CXCL12 is expressed in a subpopulation of RGCs and axonally transported, implying its release at the injury site. A different RGC subpopulation expressed CXCR4, causing axons in the distal nerve to return to the injury site. Specific depletion of CXCR4 or CXCL12 in RGCs abolished aberrant growth. It enabled long-distance regeneration in the optic nerve, with some axons reaching the optic chiasm 3 wk after injury. Thus, active CXCL12/CXCR4-mediated entrapment markedly compromises axon extension into the distal optic nerve and contributes to regenerative failure in the optic nerve.     

Identifying hydrophobic protein patches to inform protein interaction interfaces

Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.

Protein–protein interactions play a crucial role in numerous biological processes, ranging from signal transduction and immune response to protein aggregation and phase behavior (1–3). Consequently, being able to understand, predict, and modulate protein interactions has important implications for understanding cellular processes and mitigating the progression of disease (4, 5). A necessary first step toward this ambitious goal is uncovering the interfaces through which proteins interact (6–8). In principle, identifying hydrophobic protein regions, which interact weakly with water, should be a promising strategy for uncovering protein interaction interfaces (9, 10). Indeed, the release of weakly interacting hydration waters from hydrophobic regions can drive protein interactions, as well as other aqueous assemblies (11–13). However, even when the structure of a protein is available at atomistic resolution, it is challenging to identify its hydrophobic patches because they are not uniformly nonpolar, but display variations in polarity and charge at the nanoscale. Moreover, the emergent hydrophobicity of a protein patch stems from the collective response of protein hydration waters to the nanoscale chemical and topographical patterns displayed by the patch (14–20) and cannot be captured by simply counting the number of nonpolar groups in the patch, or even through more involved additive approaches, such as hydropathy scales or surface-area models (21–28).To address this challenge, we build upon seminal theoretical advances and molecular simulation studies, which have shown that near a hydrophobic surface, it is easier to disrupt surface–water interactions and form interfacial cavities (29–34). To uncover protein regions that have the weakest interactions with water, here, we employ specialized molecular simulations, wherein protein–water interactions are disrupted by systematically displacing water molecules from the protein hydration shell (35–37). By identifying the protein patches that nucleate cavities most readily in our simulations, we are then able to uncover the most hydrophobic protein regions. Interestingly, we find that both hydrophobic and hydrophilic protein patches are highly heterogeneous and contain comparable numbers of nonpolar and polar atoms. Our results thus highlight the nontrivial relationship between the chemical composition of protein patches and their emergent hydrophobicity (24–26), and further emphasize the importance of accounting for the collective solvent response in characterizing protein hydrophobicity (16). We then interrogate whether the most hydrophobic protein patches, which nucleate cavities readily, are also likely to mediate protein interactions. Across five proteins that participate in either homodimer or heterodimer formation, we find that roughly 60 to 70% of interfacial contacts and only about 10 to 20% of noncontacts nucleate cavities. Our work thus provides a versatile computational framework for characterizing hydrophobicity and uncovering interaction interfaces of not just proteins, but also of other complex amphiphilic solutes, such as cavitands, dendrimers, and patchy nanoparticles (38–41).     

From the Cover: Evidence from South Africa for a protracted end-Permian extinction on land

Earth’s largest biotic crisis occurred during the Permo–Triassic Transition (PTT). On land, this event witnessed a turnover from synapsid- to archosauromorph-dominated assemblages and a restructuring of terrestrial ecosystems. However, understanding extinction patterns has been limited by a lack of high-precision fossil occurrence data to resolve events on submillion-year timescales. We analyzed a unique database of 588 fossil tetrapod specimens from South Africa’s Karoo Basin, spanning ∼4 My, and 13 stratigraphic bin intervals averaging 300,000 y each. Using sample-standardized methods, we characterized faunal assemblage dynamics during the PTT. High regional extinction rates occurred through a protracted interval of ∼1 Ma, initially co-occurring with low origination rates. This resulted in declining diversity up to the acme of extinction near the Daptocephalus–Lystrosaurus declivis Assemblage Zone boundary. Regional origination rates increased abruptly above this boundary, co-occurring with high extinction rates to drive rapid turnover and an assemblage of short-lived species symptomatic of ecosystem instability. The “disaster taxon” Lystrosaurus shows a long-term trend of increasing abundance initiated in the latest Permian. Lystrosaurus comprised 54% of all specimens by the onset of mass extinction and 70% in the extinction aftermath. This early Lystrosaurus abundance suggests its expansion was facilitated by environmental changes rather than by ecological opportunity following the extinctions of other species as commonly assumed for disaster taxa. Our findings conservatively place the Karoo extinction interval closer in time, but not coeval with, the more rapid marine event and reveal key differences between the PTT extinctions on land and in the oceans.

Mass extinctions are major perturbations of the biosphere resulting from a wide range of different causes including glaciations and sea level fall (1), large igneous provinces (2), and bolide impacts (3, 4). These events caused permanent changes to Earth’s ecosystems, altering the evolutionary trajectory of life (5). However, links between the broad causal factors of mass extinctions and the biological and ecological disturbances that lead to species extinctions have been difficult to characterize. This is because ecological disturbances unfold on timescales much shorter than the typical resolution of paleontological studies (6), particularly in the terrestrial record (6–8). Coarse-resolution studies have demonstrated key mass extinction phenomena including high extinction rates and lineage turnover (7, 9), changes in species richness (10), ecosystem instability (11), and the occurrence of disaster taxa (12). However, finer time resolutions are central to determining the association and relative timings of these effects, their potential causal factors, and their interrelationships. Achieving these goals represents a key advance in understanding the ecological mechanisms of mass extinctions.The end-Permian mass extinction (ca. 251.9 Ma) was Earth’s largest biotic crisis as measured by taxon last occurrences (13–15). Large outpourings from Siberian Trap volcanism (2) are the likely trigger of calamitous climatic changes, including a runaway greenhouse effect and ocean acidification, which had profound consequences for life on land and in the oceans (16–18). An estimated 81% of marine species (19) and 89% of tetrapod genera became extinct as established Permian ecosystems gave way to those of the Triassic. In the ocean, this included the complete extinction of reef-forming tabulate and rugose corals (20, 21) and significant losses in previously diverse ammonoid, brachiopod, and crinoid families (22). On land, many nonmammalian synapsids became extinct (16), and the glossopterid-dominated floras of Gondwana also disappeared (23). Stratigraphic sequences document a global “coral gap” and “coal gap” (24, 25), suggesting reef and forest ecosystems were rare or absent for up to 5 My after the event (26). Continuous fossil-bearing deposits documenting patterns of turnover across the Permian–Triassic transition (PTT) on land (27) and in the oceans (28) are geographically widespread (29, 30), including marine and continental successions that are known from China (31, 32) and India (33). Continental successions are known from Russia (34), Australia (35), Antarctica (36), and South Africa’s Karoo Basin (Fig. 1 and 37–40), the latter providing arguably the most densely sampled and taxonomically scrutinized (41–43) continental record of the PTT. The main extinction has been proposed to occur at the boundary between two biostratigraphic zones with distinctive faunal assemblages, the Daptocephalus and Lystrosaurus declivis assemblage zones (Fig. 1), which marks the traditional placement of the Permian–Triassic geologic boundary [(37) but see ref. 44]. Considerable research has attempted to understand the anatomy of the PTT in South Africa (38, 39, 45–52) and to place it in the context of biodiversity changes across southern Gondwana (53, 54) and globally (29, 31, 32, 44, 47, 55).Open in a separate windowFig. 1.Map of South Africa depicting the distribution of the four tetrapod fossil assemblage zones (Cistecephalus, Daptocephalus, Lystrosaurus declivis, Cynognathus) and our two study sites where fossils were collected in this study (sites A and B). Regional lithostratigraphy and biostratigraphy within the study interval are shown alongside isotope dilution–thermal ionization mass spectrometry dates retrieved by Rubidge et al., Botha et al., and Gastaldo et al. (37, 44, 80). The traditional (dashed red line) and associated PTB hypotheses for the Karoo Basin (37, 44) are also shown. Although traditionally associated with the PTB, the Daptocephalus–Lystrosaurus declivis Assemblage Zone boundary is defined by first appearances of co-occurring tetrapod assemblages, so its position relative to the three PTB hypotheses is unchanged. The Ripplemead member (*) has yet to be formalized by the South African Committee for Stratigraphy.Decades of research have demonstrated the richness of South Africa’s Karoo Basin fossil record, resulting in hundreds of stratigraphically well-documented tetrapod fossils across the PTT (37, 39, 56). This wealth of data has been used qualitatively to identify three extinction phases and an apparent early postextinction recovery phase (39, 45, 51). Furthermore, studies of Karoo community structure and function have elucidated the potential role of the extinction and subsequent recovery in breaking the incumbency of previously dominant clades, including synapsids (11, 57). Nevertheless, understanding patterns of faunal turnover and recovery during the PTT has been limited by the scarcity of quantitative investigations. Previous quantitative studies used coarsely sampled data (i.e., assemblage zone scale, 2 to 3 Ma time intervals) to identify low species richness immediately after the main extinction, potentially associated with multiple “boom and bust” cycles of primary productivity based on δ¹³C variation during the first 5 My of the Triassic (41, 58). However, many details of faunal dynamics in this interval remain unknown. Here, we investigate the dynamics of this major tetrapod extinction at an unprecedented time resolution (on the order of hundreds of thousands of years), using sample-standardized methods to quantify multiple aspects of regional change across the Cistecephalus, Daptocephalus, and Lystrosaurus declivis assemblage zones.     

Evolution and global charge conservation for polarization singularities emerging from non-Hermitian degeneracies

Core concepts in singular optics, especially the polarization singularities, have rapidly penetrated the surging fields of topological and non-Hermitian photonics. For open photonic structures with non-Hermitian degeneracies in particular, polarization singularities would inevitably encounter another sweeping concept of Berry phase. Several investigations have discussed, in an inexplicit way, connections between both concepts, hinting at that nonzero topological charges for far-field polarizations on a loop are inextricably linked to its nontrivial Berry phase when degeneracies are enclosed. In this work, we reexamine the seminal photonic crystal slab that supports the fundamental two-level non-Hermitian degeneracies. Regardless of the invariance of nontrivial Berry phase (concerning near-field Bloch modes defined on the momentum torus) for different loops enclosing both degeneracies, we demonstrate that the associated far polarization fields (defined on the momentum sphere) exhibit topologically inequivalent patterns that are characterized by variant topological charges, including even the trivial scenario of zero charge. Moreover, the charge carried by the Fermi arc actually is not well defined, which could be different on opposite bands. It is further revealed that for both bands, the seemingly complex evolutions of polarizations are bounded by the global charge conservation, with extra points of circular polarizations playing indispensable roles. This indicates that although not directly associated with any local charges, the invariant Berry phase is directly linked to the globally conserved charge, physical principles underlying which have all been further clarified by a two-level Hamiltonian with an extra chirality term. Our work can potentially trigger extra explorations beyond photonics connecting Berry phase and singularities.

Pioneered by Pancharatnam, Berry, Nye, and others (1–10), Berry phase and singularities have become embedded languages all across different branches of photonics. Optical Berry phase is largely manifested through either polarization evolving Pancharatnam–Berry phase or the spin-redirection Bortolotti–Rytov–Vladimirskii–Berry phase (2, 4, 5, 11–15); while optical singularities are widely observed as singularities of intensity (caustics) (6), phase (vortices) (7), or polarization (8–10). As singularities for complex vectorial waves, polarization singularities are skeletons of electromagnetic waves and are vitally important for understanding various interference effects underlying many applications (16–20).There is a superficial similarity between the aforementioned two concepts: Both the topological charge of polarization field [Hopf index of line field (21, 22)] and Berry phase are defined on a closed circuit. Despite this, it is quite unfortunate that almost no definite connections have been established between them in optics. This is fully understandable: Berry phase is defined on the Pancharatnam connection (parallel transport) that decides the phase contrast between neighboring states on the loop (3, 4); while the polarization charge reflects accumulated orientation rotations of polarization ellipses, which has no direct relevance to the overall phase of each state. This explains why in pioneering works where both concepts were present (23–27), their interplay was rarely elaborated on.Spurred by studies into bound states in the continuum, polarization singularities have gained enormous renewed interest in open periodic photonic structures, manifested in different morphologies with both fundamental and higher-order half-integer charges (28–50). Simultaneously, the significance of Berry phase has been further reinforced in surging fields of topological and non-Hermitian photonics (1, 23, 26, 51–56). In open periodic structures involving band degeneracies, Berry phase and polarization singularity would inevitably meet, which sparks the influential work on non-Hermitian degeneracy (36) and several other following studies (40, 43, 45) discussing both concepts simultaneously. Although not claimed explicitly, those works hint that nontrivial Berry phase produces nonzero polarization charges.Aiming to bridge Berry phase and polarization singularity, we reexamine the seminal photonic crystal slab (PCS) that supports elementary two-level non-Hermitian degeneracies. It is revealed that with an invariant nontrivial π Berry phase, the corresponding polarization fields on different isofrequency contours enclosing both non-Hermitian degenerate points (or equivalently exceptional points [EPs]) (26) exhibit diverse patterns characterized by different polarization charges, even including the trivial zero charge. It is further revealed that the charge carried by the Fermi arc is actually not well defined, which could be different on opposite bands. We also discover that such complexity of field evolutions is constrained by global charge conservation for both bands, with extra points of circular polarizations (C points) playing pivotal roles. This reveals the explicit connection between globally conserved charge and the invariant Berry phase, underlying which the physical mechanisms have been further clarified by a two-level Hamiltonian with an extra chirality term (25). We show that such an unexpected connection is generically manifest in various structures, despite the fact that Berry phase and polarization charge actually characterize different entities of near-field Bloch modes and their projected far polarization fields, respectively: Bloch modes are defined on the momentum torus and can be folded into the irreducible Brillouin zone; while polarization fields are defined on the momentum sphere, due to the involvement of out-of-plane wave vectors along which there is no periodicity. Our study can spur further investigations in other subjects beyond photonics to explore conceptual interconnectedness, where both the concepts of Berry phase and singularities are present.     

From the Cover: Phase space transport in the interaction between shocks and plasma turbulence

The interaction of collisionless shocks with fully developed plasma turbulence is numerically investigated. Hybrid kinetic simulations, where a turbulent jet is slammed against an oblique shock, are employed to address the role of upstream turbulence on plasma transport. A technique, using coarse graining of the Vlasov equation, is proposed, showing that the particle transport strongly depends on upstream turbulence properties, such as strength and coherency. These results might be relevant for the understanding of acceleration and heating processes in space plasmas.

A turbulent plasma wind flows from the sun and permeates the heliosphere, encountering several magnetic obstacles, leading to shocks that continuously interact with the incoming complex solar wind—a scenario that becomes a prototype for understanding many other systems characterized by the presence of shocks. Despite decades of research, the interaction of shocks with plasma turbulence and the subsequent energetic particle production still remain poorly understood (1, 2). Shocks are well-known efficient, natural particle accelerators (3) and have been modeled in a number of theories (4–8). Less understood is the interaction of shocks and turbulence that characterizes spectacular high-energy events, as in supernovae explosions propagating through the interstellar turbulent medium, as in the case of coronal mass ejections that stream through the turbulent solar wind, and as for the complex Earth’s bow shock environment. In many of the above examples, oblique shocks are known to generate coherent field-aligned beams (FABs), as observed at Earth’s bow shock (9). FABs are an important source of free energy throughout the interplanetary medium (10). Turbulence-generated coherent structures and waves might interact with the shock discontinuity, in an interplay that is likely to play a pivotal role in particle acceleration and plasma heating (11–13).Turbulence is populated by a variety of structures that can work effectively as particle “traps” and “corridors” that either hinder or enable their motion (14) and represents another crucial source of accelerated particles (15–18). An example of such an energization process has been observed in the patterns of local reconnection that develop in turbulence (19, 20). In order to understand such mechanisms, the transport properties need to be explored in the plasma phase space (21).Due to the difference between the spatial and temporal scales involved in accelerating particles, shocks and turbulence are often considered theoretically in isolation rather than together. However, fundamental studies have suggested that these are inextricably linked: Shocks are likely to propagate in turbulent media, and turbulence is responsible for changing fundamental aspects of shock transitions (22–27). Inspired by these studies, here we quantitatively explore the intimate relation between these two phenomena.     

Cognitive aging is associated with redistribution of synaptic weights in the hippocampus

Behaviors that rely on the hippocampus are particularly susceptible to chronological aging, with many aged animals (including humans) maintaining cognition at a young adult-like level, but many others the same age showing marked impairments. It is unclear whether the ability to maintain cognition over time is attributable to brain maintenance, sufficient cognitive reserve, compensatory changes in network function, or some combination thereof. While network dysfunction within the hippocampal circuit of aged, learning-impaired animals is well-documented, its neurobiological substrates remain elusive. Here we show that the synaptic architecture of hippocampal regions CA1 and CA3 is maintained in a young adult-like state in aged rats that performed comparably to their young adult counterparts in both trace eyeblink conditioning and Morris water maze learning. In contrast, among learning-impaired, but equally aged rats, we found that a redistribution of synaptic weights amplifies the influence of autoassociational connections among CA3 pyramidal neurons, yet reduces the synaptic input onto these same neurons from the dentate gyrus. Notably, synapses within hippocampal region CA1 showed no group differences regardless of cognitive ability. Taking the data together, we find the imbalanced synaptic weights within hippocampal CA3 provide a substrate that can explain the abnormal firing characteristics of both CA3 and CA1 pyramidal neurons in aged, learning-impaired rats. Furthermore, our work provides some clarity with regard to how some animals cognitively age successfully, while others’ lifespans outlast their “mindspans.”

Aging is the biggest risk factor for Alzheimer’s disease, but many aged individuals nevertheless retain the ability to perform cognitive tasks with young adult (YA)-like competency, and are thus resilient to age-related cognitive decline and dementias (1, 2). The mechanisms of such resilience are unknown, but are thought to involve neural or cognitive reserve, brain network adaptations, or simply the ability to maintain cognitive brain circuits in a YA-like state (3–5). Much of the cellular and functional insight into the concept or risk of/resilience against age-related cognitive impairments has come from animal models of normal/nonpathological aging (6–10). Many of these studies have shown that circuit function abnormalities are associated with behavioral impairments. The cellular and structural bases for such functional aberrations, however, remain largely unknown.Two of the most well-studied cognitive domains that show susceptibility to chronological aging in both rodents and nonhuman primates are working memory and spatial/temporal memory (6–10). Importantly, these cognitive domains engage anatomically distinct neurocognitive systems, with the former relying on prefrontal/orbitofrontal cortical circuits and the latter relying on hippocampal circuitry. Interestingly, although behavioral deficits in these two domains (in the case of rat models of cognitive aging) begin to emerge, worsen, and become increasingly prevalent between 12 and 18 mo of age in most strains (reviewed in refs. 9 and 11), cognitive aging within hippocampus-dependent forms of learning and memory are relatively independent of those that engage the prefrontal/orbitofrontal cortical neural systems (6–9, 12–15).Neither the mechanisms underlying the conservation of memory function across chronological aging nor those contributing to the age-related emergence and exacerbation of memory impairments are clearly understood for either neurocognitive system. It is clear, however, that neither frank neuronal loss (16, 17) nor overall synapse loss (18) contributes to cognitive aging within the medial temporal lobe/hippocampal memory system. Rather, the intriguing idea that has emerged from work in both the hippocampal and the prefrontal/orbitofrontal cortical memory systems is that there are functional alterations in the synaptic connections in individual microcircuits embedded within these larger neuroanatomical systems (6–10, 19–31).Axospinous synapses (including those in hippocampal and cortical circuits) are characterized on the basis of the three-dimensional morphology of their postsynaptic densities (PSDs) (20, 32–34). The most-abundant axospinous synaptic subtype has a continuous, macular, disk-shaped PSD, as compared to the less-abundant perforated synaptic subtype, which has at least one discontinuity in its PSD (34). In addition to differing substantially with regard to relative frequency, perforated and nonperforated synapses also harbor major differences in size and synaptic AMPA-type and NMDA-type receptor expression levels (AMPAR and NMDAR, respectively) (34–38). There is also evidence that perforated and nonperforated synapses are differentially involved in synaptic plasticity (39–44) and in preservation of—or reductions in—memory function during chronological aging (6, 20, 45). Layered onto these general distinctions between perforated and nonperforated synapses are more specific differences in their characteristics when considered within neural circuits. For example, perforated synapses have a stronger and more consistent influence on neuronal computation within hippocampal region CA1 than their nonperforated counterparts, which nevertheless outnumber the former by a roughly 9-to-1 ratio (34, 46, 47).These and other circuit-specific differences necessitate a circuit-based approach to understanding the synaptic bases underlying the retention or loss of YA-like memory function in the aging brain. In many ways, the hippocampal system is particularly convenient for such circuit-based approaches (48, 49). Information about the internal and external world is funneled to the parahippocampal system and then relayed via the entorhinal cortex to the dentate gyrus, the first component of the so-called trisynaptic circuit in the hippocampus proper. Granule cells in the dentate gyrus then transmit their computations to hippocampal region CA3 via the mossy fibers, which form very large and anatomically distinct synapses called mossy fiber bouton–thorny excrescence synaptic complexes in the stratum lucidum (SL). CA3 pyramidal neurons then integrate information from their autoassociational connections in the stratum radiatum (SR) and stratum oriens (SO), with both direct entorhinal inputs in stratum lacunosum-moleculare (SLM) and the dentate gyrus inputs in the SL, and convey this information to hippocampal CA1 pyramidal neurons. Neurons in hippocampal CA1 then integrate this information in their basal and apical SR dendrites with direct entorhinal cortical inputs in their most distal, tufted dendrites in the SLM, and represent the first and largest extrahippocampal output from the hippocampus proper. Thus, the computations performed both within individual hippocampal subregions and between them as an interconnected neurocognitive system are complex, and involve a combination of intrinsic (i.e., membrane-bound ion channels that regulate membrane excitability) and synaptic (i.e., ligand-gated ion channels expressed at both excitatory and inhibitory synapses) influences. Additionally, age-related changes at any level of these complex circuits will have downstream consequences on the accuracy/reliability of the information being relayed to extrahippocampal regions via CA1 pyramidal neurons.Given the amount of evidence supporting a possible synaptic explanation for age-related learning and memory impairments in hippocampus-dependent forms of cognition (6–10), we combined patch-clamp physiology, serial section conventional and immunogold electron microscopy (EM), quantitative Western blot analyses, and behavioral characterization using two hippocampus-dependent forms of learning in YA (6- to 8-mo old) and aged rats (28- to 29-mo old) to examine two interconnected hippocampal regions implicated in cognitive aging: Regions CA1 and CA3. We focused on CA1 and CA3 because of their central location in the hippocampal circuit (48–50), their similar laminar dendritic structure (48–50), and their well-documented age-related changes in place field specificity and reliability (51–56).We find that the synaptic architecture and balance of synaptic weights in YA and aged, learning-unimpaired (AU) rats is remarkably similar, but that both are different in aged, learning-impaired (AI) rats. Moreover, this restructuring among “unsuccessful” cognitive agers has an intriguing specificity: It involves only AMPARs, only perforated axospinous synapses, and only hippocampal region CA3, which together shift the balance of synaptic weights that drive action potential output in CA3 pyramidal neurons maladaptively toward an overemphasis of the autoassociational synapses that interconnect CA3 pyramidal neurons.     

From the Cover: Radiation-induced neoantigens broaden the immunotherapeutic window of cancers with low mutational loads

Immunotherapies are a promising advance in cancer treatment. However, because only a subset of cancer patients benefits from these treatments it is important to find mechanisms that will broaden the responding patient population. Generally, tumors with high mutational burdens have the potential to express greater numbers of mutant neoantigens. As neoantigens can be targets of protective adaptive immunity, highly mutated tumors are more responsive to immunotherapy. Given that external beam radiation 1) is a standard-of-care cancer therapy, 2) induces expression of mutant proteins and potentially mutant neoantigens in treated cells, and 3) has been shown to synergize clinically with immune checkpoint therapy (ICT), we hypothesized that at least one mechanism of this synergy was the generation of de novo mutant neoantigen targets in irradiated cells. Herein, we use Kras^G12D x p53^−/− sarcoma cell lines (KP sarcomas) that we and others have shown to be nearly devoid of mutations, are poorly antigenic, are not controlled by ICT, and do not induce a protective antitumor memory response. However, following one in vitro dose of 4- or 9-Gy irradiation, KP sarcoma cells acquire mutational neoantigens and become sensitive to ICT in vivo in a T cell-dependent manner. We further demonstrate that some of the radiation-induced mutations generate cytotoxic CD8⁺ T cell responses, are protective in a vaccine model, and are sufficient to make the parental KP sarcoma line susceptible to ICT. These results provide a proof of concept that induction of new antigenic targets in irradiated tumor cells represents an additional mechanism explaining the clinical findings of the synergy between radiation and immunotherapy.

Immune checkpoint therapy (ICT) can lead to durable responses in subsets of cancer patients (1–8). On the basis of computational analyses, the patients who most benefit from ICT are those with cancers that have high mutational burden (9–18). For example, patients bearing tumors with high mutational burden caused by environmental exposure (such as ultraviolet-induced melanoma) or deficiencies in DNA repair (such as microsatellite instability-high colorectal cancers) tend to respond well to immunotherapy (18–26). Presumably the sensitivity of such cancers reflects the increased likelihood of formation of immunogenic, tumor-specific mutant neoantigens (27). We and others previously showed that certain tumor-specific neoantigens are major targets of natural and therapeutically induced antitumor responses in both mice and humans (28–41). Therefore, the presence of immunogenic tumor neoantigens is currently thought to contribute to tumor sensitivity to immunotherapy.However, many cancer patients do not respond to ICT, suggesting that their neoantigen burden is either of insufficient magnitude or immunogenicity to function as targets for T cell-dependent antitumor mechanisms. Indeed, there are many tumor types, such as acute myeloid leukemia, estrogen receptor-positive breast, and prostate cancers, that have limited mutational burdens and display low response rates to ICT (9, 13, 42, 43). Additionally, tumor cell clones expressing immunogenic neoantigens that develop during tumor evolution may be eliminated from tumors with high mutational burden by the process of cancer immunoediting, resulting in outgrowth of tumor cell clones with reduced immunogenicity that can then grow progressively in the presence of the unmanipulated immune system (33, 44, 45). Therefore, a process by which tumors with low neoantigen burden can acquire immunogenic mutations has the potential to expand the number of patients able to benefit from ICT.Ionizing radiation has been shown to elicit DNA damage in tumor cells, leading to an increase in overall mutational load (46–52). This damage is thought to occur primarily through generation of reactive oxygen species which induce base pair substitutions by mechanisms involving transitions, transversions, and/or faulty DNA repair (53). Multiple preclinical studies have demonstrated antitumor responses when focal radiation is combined with ICT in tumors that do not respond to ICT alone (54–60) and several clinical studies have demonstrated that human tumor patients have improved responsiveness to ICT following focal radiation (e.g., {"type":"clinical-trial","attrs":{"text":"NCT02303990","term_id":"NCT02303990"}}NCT02303990, {"type":"clinical-trial","attrs":{"text":"NCT02298946","term_id":"NCT02298946"}}NCT02298946, {"type":"clinical-trial","attrs":{"text":"NCT02383212","term_id":"NCT02383212"}}NCT02383212) (61–67). Radiation has been demonstrated to function as an in vivo tumor vaccine by inducing damage-associated molecular patterns (DAMP)-dependent immunogenic cell death (68), inducing DNA damage sensed by pattern recognition receptors (69, 70), enhancing access of immune effector cells to their cognate targets through tumor cell debulking and vasculature changes (71, 72), up-regulating major histocompatibility complex class I (MHC-I) receptors (73), up-regulating cell-surface molecules such as Fas (74), and augmenting tumor antigen cross-presentation by specific subsets of dendritic cells through up-regulation of type I interferon (IFN), which results in increased numbers and action of tumor-specific CD8⁺ T cells (75–77). However, none of these explanations take into account that following irradiation, tumor cells acquire novel mutations that may function as effective tumor neoantigens. In fact, two groups have demonstrated broadening of the T cell repertoire following radiation treatment of mouse 4T1 mammary tumors and B16F10 melanoma tumors (56, 78). Radiation-induced neoantigens may partially explain the broadening of the T cell repertoire reported during noncurative doses of irradiation.Given the above observations, we specifically explored whether one dose of in vitro irradiation could increase the immunogenicity of poorly immunogenic tumor cell lines through mechanisms involving the de novo generation of tumor-specific mutant neoantigens. For this purpose, we used a mouse Kras^G12D x p53^−/− sarcoma cell line as a model system since the R.D.S. and T.J. laboratories have previously shown that these tumor cells express a very limited number of somatic mutations, are essentially devoid of mutational neoantigens, and are nonimmunogenic and grow progressively in syngeneic wild-type (WT) mice either following treatment with control antibody or the combination of anti–PD-1/anti–CTLA-4 (34, 41). We find that treating these cell lines with noncurative doses of irradiation induces expression of somatic mutations, some of which function as neoantigens and render the sarcoma cells susceptible to ICT in vivo. These data support the concept that an additional mechanism underlying the synergy between radiation therapy and immunotherapy is that the former induces immunogenic mutations in tumors that now function as targets for the latter.     

Replisome bypass of a protein-based R-loop block by Pif1

Efficient and faithful replication of the genome is essential to maintain genome stability. Replication is carried out by a multiprotein complex called the replisome, which encounters numerous obstacles to its progression. Failure to bypass these obstacles results in genome instability and may facilitate errors leading to disease. Cells use accessory helicases that help the replisome bypass difficult barriers. All eukaryotes contain the accessory helicase Pif1, which tracks in a 5′–3′ direction on single-stranded DNA and plays a role in genome maintenance processes. Here, we reveal a previously unknown role for Pif1 in replication barrier bypass. We use an in vitro reconstituted Saccharomyces cerevisiae replisome to demonstrate that Pif1 enables the replisome to bypass an inactive (i.e., dead) Cas9 (dCas9) R-loop barrier. Interestingly, dCas9 R-loops targeted to either strand are bypassed with similar efficiency. Furthermore, we employed a single-molecule fluorescence visualization technique to show that Pif1 facilitates this bypass by enabling the simultaneous removal of the dCas9 protein and the R-loop. We propose that Pif1 is a general displacement helicase for replication bypass of both R-loops and protein blocks.

Efficient and faithful replication of the genome is essential to maintain genome stability and is carried out by a multiprotein complex called the replisome (1–4). There are numerous obstacles to progression of the replisome during the process of chromosome duplication. These obstacles include RNA-DNA hybrids (R-loops), DNA secondary structures, transcribing RNA polymerases, and other tightly bound proteins (5–9). Failure to bypass these barriers may result in genome instability, which can lead to cellular abnormalities and genetic disease. Cells contain various accessory helicases that help the replisome bypass these difficult barriers (10–20). A subset of these helicases act on the opposite strand of the replicative helicase (1, 2, 14, 19).All eukaryotes contain an accessory helicase, Pif1, which tracks in a 5′–3′ direction on single-stranded DNA (ssDNA) (11–16). Pif1 is important in pathways such as Okazaki-fragment processing and break-induced repair that require the removal of DNA-binding proteins as well as potential displacement of R-loops (11–13, 21, 15–18, 22–25). Genetic studies and immunoprecipitation pull-down assays indicate that Pif1 interacts with PCNA (the DNA sliding clamp), Pol ε (the leading-strand polymerase), the MCMs (the motor subunits of the replicative helicase CMG), and RPA (the single-stranded DNA-binding protein) (15, 26, 27). Pif1 activity in break-induced repair strongly depends on its interaction with PCNA (26). These interactions with replisomal components suggest that Pif1 could interact with the replisome during replication. In Escherichia coli, the replicative helicase is the DnaB homohexamer that encircles the lagging strand and moves in a 5′–3′ direction (20). E. coli accessory helicases include the monomeric UvrD (helicase II) and Rep, which move in the 3′–5′ direction and operate on the opposite strand from the DnaB hexamer. It is known that these monomeric helicases promote the bypass of barriers during replication such as stalled RNA polymerases (5). The eukaryotic replicative helicase is the 11-subunit CMG (Cdc45, Mcm2–7, GINS) and tracks in the 3′–5′ direction, opposite to the direction of Pif1 (25, 28). Once activated by Mcm10, the MCM motor domains of CMG encircle the leading strand (29–32). We hypothesized that, similar to UvrD and Rep in E. coli, Pif1 interacts with the replisome tracking in the opposite direction to enable bypass of replication obstacles.In this report, we use an in vitro reconstituted Saccharomyces cerevisiae replisome to study the role of Pif1 in bypass of a “dead” Cas9 (dCas9), which is a Cas9 protein that is deactivated in DNA cleavage but otherwise fully functional in DNA binding. As with Cas9, dCas9 is a single-turnover enzyme that can be programmed with a guide RNA (gRNA) to target either strand. The dCas9–gRNA complex forms a roadblock consisting of an R-loop and a tightly bound protein (dCas9), a construct that is similar to a stalled RNA polymerase. This roadblock (hereafter dCas9 R-loop) arrests replisomes independent of whether the dCas9 R-loop is targeted to the leading or lagging strand (30). Besides its utility due to its programmable nature (33), the use of the dCas9 R-loop allows us to answer several mechanistic questions. For example, the ability to program the dCas9 R-loop block to any specific sequence enables us to observe whether block removal is different depending on whether the block is on the leading or lagging strand. Furthermore, the inner diameter of CMG can accommodate double-stranded DNA (dsDNA) and possibly an R-loop, but not a dCas9 protein. Using the dCas9 R-loop block allows us to determine the fate of each of its components.Here, we report that Pif1 enables the bypass of the dCas9 R-loop by the replisome. Interestingly, dCas9 R-loops targeted to either the leading or lagging strand are bypassed with similar efficiency. In addition, the PCNA clamp is not required for bypass of the block, indicating that Pif1 does not need to interact with PCNA during bypass of the block. We used a single-molecule fluorescence imaging to show that both the dCas9 and the R-loop are displaced as an intact nucleoprotein complex. We propose that Pif1 is a general displacement helicase for replication bypass of both R-loops and protein blocks.     

Environmental noise degrades hippocampus-related learning and memory

The neural mechanisms underlying the impacts of noise on nonauditory function, particularly learning and memory, remain largely unknown. Here, we demonstrate that rats exposed postnatally (between postnatal days 9 and 56) to structured noise delivered at a sound pressure level of ∼65 dB displayed significantly degraded hippocampus-related learning and memory abilities. Noise exposure also suppressed the induction of hippocampal long-term potentiation (LTP). In parallel, the total or phosphorylated levels of certain LTP-related key signaling molecules in the synapses of the hippocampus were down-regulated. However, no significant changes in stress-related processes were found for the noise-exposed rats. These results in a rodent model indicate that even moderate-level noise with little effect on stress status can substantially impair hippocampus-related learning and memory by altering the plasticity of synaptic transmission. They support the importance of more thoroughly defining the unappreciated hazards of moderately loud noise in modern human environments.

The noise pollution accompanying industrialization is a risk factor to human health. Earlier studies have extensively examined the deleterious impacts of noise in the auditory systems of both humans and animal models (1–6), showing that noise exposure either early or late in life can induce progressive hearing loss, change neural coding along the auditory pathway, and alter auditory-related perception and behavior.The auditory system, however, contains direct and indirect pathways to other systems and structures of the brain that are necessary for functional integration. For example, earlier studies found that the hippocampus, the core area of the brain associated with learning and memory processes, receives neuronal inputs from the auditory system through the lemniscal and nonlemniscal pathways (7–11). It is thus conceivable that noise-evoked activities might be transmitted via these connections to the hippocampus, thereby affecting learning and memory. Indeed, animal studies have shown that exposure to loud noise (e.g., above a sound pressure level [SPL] of 95 dB) that induces temporary or permanent shifts in the auditory threshold disrupts hippocampal histology, decreases neurogenesis in the hippocampus, and impairs hippocampus-related learning and memory abilities (12–16). In addition, epidemiological studies have demonstrated that environment noise has substantially negative effects on children’s learning outcomes and cognitive abilities (17–19). While the usual explanations for the origins of these noise-induced effects on nonauditory functions have relied on stress-related processes (15, 16, 20–23), the underlying neural mechanisms remain largely unknown.In this study, we exposed rat pups to structured noise delivered at ∼65 dB SPL for a 7-wk period. Exposure to a moderate level of modulated broad-spectrum noise more realistically models the noise environments people encountered in industrial workplaces and other modern acoustic settings (2, 4, 24–26). We then evaluated the behavioral consequences of noise exposure on hippocampus-related learning and memory for these noise-exposed rats. In addition, we explored the mechanisms underlying possible postexposure changes in learning and memory via physiological and molecular assessments of the hippocampus.     

Photosynthesis tunes quantum-mechanical mixing of electronic and vibrational states to steer exciton energy transfer

Photosynthetic species evolved to protect their light-harvesting apparatus from photoxidative damage driven by intracellular redox conditions or environmental conditions. The Fenna–Matthews–Olson (FMO) pigment–protein complex from green sulfur bacteria exhibits redox-dependent quenching behavior partially due to two internal cysteine residues. Here, we show evidence that a photosynthetic complex exploits the quantum mechanics of vibronic mixing to activate an oxidative photoprotective mechanism. We use two-dimensional electronic spectroscopy (2DES) to capture energy transfer dynamics in wild-type and cysteine-deficient FMO mutant proteins under both reducing and oxidizing conditions. Under reducing conditions, we find equal energy transfer through the exciton 4–1 and 4–2-1 pathways because the exciton 4–1 energy gap is vibronically coupled with a bacteriochlorophyll-a vibrational mode. Under oxidizing conditions, however, the resonance of the exciton 4–1 energy gap is detuned from the vibrational mode, causing excitons to preferentially steer through the indirect 4–2-1 pathway to increase the likelihood of exciton quenching. We use a Redfield model to show that the complex achieves this effect by tuning the site III energy via the redox state of its internal cysteine residues. This result shows how pigment–protein complexes exploit the quantum mechanics of vibronic coupling to steer energy transfer.

Photosynthetic organisms convert solar photons into chemical energy by taking advantage of the quantum mechanical nature of their molecular systems and the chemistry of their environment (1–4). Antenna complexes, composed of one or more pigment–protein complexes, facilitate the first steps in the photosynthesis process: They absorb photons and determine which proportion of excitations to move to reaction centers, where charge separation occurs (4). In oxic environments, excitations can generate highly reactive singlet oxygen species. These pigment–protein complexes can quench excess excitations in these environments with molecular moieties such as quinones and cysteine residues (1, 5–7).The Fenna–Matthews–Olson (FMO) complex, a trimer of pigment–protein complexes found in the green sulfur bacterium Chlorobaculum tepidum (8), has emerged as a model system to study the photophysical properties of photosynthetic antenna complexes (9–19). Each subunit in the FMO complex contains eight bacteriochlorophyll-a site molecules (Protein Data Bank, ID code: 3ENI) that are coupled to form a basis of eight partially delocalized excited states called excitons (Fig. 1) (20–23). Previous experiments on FMO have observed the presence of long-lived coherences in nonlinear spectroscopic signals at both cryogenic and physiological temperatures (11, 13). The coherent signals are thought to arise from some combination of electronic (24–26), vibrational (16–18), and vibronic (27) coherences in the system (28–30). One previous study reported that the coherent signals in FMO remain unchanged upon mutagenesis of the protein, suggesting that the signals are ground state vibrational coherences (17). Others discuss the role of vibronic coupling, where electronic and nuclear degrees of freedom become coupled (29). Other dimeric model systems have demonstrated the regimes in which these vibronically coupled states produce coherent or incoherent transport and vibronic coherences (31–33). Recent spectroscopic data has suggested that vibronic coupling plays a role in driving efficient energy transfer through photosynthetic complexes (27, 31, 33, 34), but to date there is no direct experimental evidence suggesting that biological systems use vibronic coupling as part of their biological function.Open in a separate windowFig. 1.(Left) Numbered sites and sidechains of cysteines C353 and C49 in the FMO pigment–protein complex (PDB ID code: 3ENI) (20). (Right) Site densities for excitons 4, 2, and 1 in reducing conditions with the energy transfer branching ratios for the WT oxidized and reduced protein. The saturation of pigments in each exciton denotes the relative contribution number to the exciton. The C353 residue is located near excitons 4 and 2, which have most electron density along one side of the complex, and other redox-active residues such as the Trp/Tyr chain. C353 and C49 surround site III, which contains the majority of exciton 1 density. Excitons 2 and 4 are generally delocalized over sites IV, V, and VII.It has been shown that redox conditions affect excited state properties in pigment-protein complexes, yet little is known about the underlying microscopic mechanisms for these effects (1, 9). Many commonly studied light-harvesting complexes—including the FMO complex (20), light-harvesting complex 2 (LH2) (35), the PC645 phycobiliprotein (36), and the cyanobacterial antenna complex isiA (37)—contain redox-active cysteine residues in close proximity to their chromophores. As the natural low light environment of C. tepidum does not necessitate photoprotective responses to light quantity and quality, its primary photoprotective mechanism concerns its response to oxidative stress. C. tepidum is an obligate anaerobe, but the presence of many active anoxygenic genes such as sodB for superoxide dismutase and roo for rubredoxin oxygen oxidoreductase (38) suggests that it is frequently exposed to molecular oxygen (7, 39). Using time-resolved fluorescence measurements, Orf et al. demonstrated that two cysteine residues in the FMO complex, C49 and C353, quench excitons under oxidizing conditions (1), which could protect the excitation from generating reactive oxygen species (7, 40–42). In two-dimensional electronic spectroscopy (2DES) experiments, Allodi et al. showed that redox conditions in both the wild-type and C49A/C353A double-mutant proteins affect the ultrafast dynamics through the FMO complex (9, 43). The recent discovery that many proteins across the evolutionary landscape possess chains of tryptophan and tyrosine residues provides evidence that these redox-active residues may link the internal protein behavior with the chemistry of the surrounding environment (41, 43).In this paper, we present data showing that pigment–protein complexes tune the vibronic coupling of their chromophores and that the absence of this vibronic coupling activates an oxidative photoprotective mechanism. We use 2DES to show that a pair of cysteine residues in FMO, C49 and C353, can steer excitations toward quenching sites in oxic environments. The measured reaction rate constants demonstrate unusual nonmonotonic behavior. We then use a Redfield model to determine how the exciton energy transfer (EET) time constants arise from changing chlorophyll site energies and their system-bath couplings (44, 45). The analysis reveals that the cysteine residues tune the resonance between exciton 4–1 energy gap and an intramolecular chlorophyll vibration in reducing conditions to induce vibronic coupling and detune the resonance in oxidizing conditions. This redox-dependent modulation of the vibronic coupling steers excitations through different pathways in the complex to change the likelihood that they interact with exciton quenchers.     

A general model of hippocampal and dorsal striatal learning and decision making

Humans and other animals use multiple strategies for making decisions. Reinforcement-learning theory distinguishes between stimulus–response (model-free; MF) learning and deliberative (model-based; MB) planning. The spatial-navigation literature presents a parallel dichotomy between navigation strategies. In “response learning,” associated with the dorsolateral striatum (DLS), decisions are anchored to an egocentric reference frame. In “place learning,” associated with the hippocampus, decisions are anchored to an allocentric reference frame. Emerging evidence suggests that the contribution of hippocampus to place learning may also underlie its contribution to MB learning by representing relational structure in a cognitive map. Here, we introduce a computational model in which hippocampus subserves place and MB learning by learning a “successor representation” of relational structure between states; DLS implements model-free response learning by learning associations between actions and egocentric representations of landmarks; and action values from either system are weighted by the reliability of its predictions. We show that this model reproduces a range of seemingly disparate behavioral findings in spatial and nonspatial decision tasks and explains the effects of lesions to DLS and hippocampus on these tasks. Furthermore, modeling place cells as driven by boundaries explains the observation that, unlike navigation guided by landmarks, navigation guided by boundaries is robust to “blocking” by prior state–reward associations due to learned associations between place cells. Our model, originally shaped by detailed constraints in the spatial literature, successfully characterizes the hippocampal–striatal system as a general system for decision making via adaptive combination of stimulus–response learning and the use of a cognitive map.

Behavioral and neuroscientific studies suggest that animals can apply multiple strategies to the problem of maximizing future reward, referred to as the reinforcement-learning (RL) problem (1, 2). One strategy is to build a model of the environment that can be used to simulate the future to plan optimal actions (3) and the past for episodic memory (4–6). An alternative, model-free (MF) approach uses trial and error to estimate a direct mapping from the animal’s state to its expected future reward, which the agent caches and looks up at decision time (7, 8), potentially supporting procedural memory (9). This computation is thought to be carried out in the brain through prediction errors signaled by phasic dopamine responses (10). These strategies are associated with different tradeoffs (2). The model-based (MB) approach is powerful and flexible, but computationally expensive and, therefore, slow at decision time. MF methods, in contrast, enable rapid action selection, but these methods learn slowly and adapt poorly to changing environments. In addition to MF and MB methods, there are intermediate solutions that rely on learning useful representations that reduce burdens on the downstream RL process (11–13).In the spatial-memory literature, a distinction has been observed between “response learning” and “place learning” (14–16). When navigating to a previously visited location, response learning involves learning a sequence of actions, each of which depends on the preceding action or sensory cue (expressed in egocentric terms). For example, one might remember a sequence of left and right turns starting from a specific landmark. An alternative place-learning strategy involves learning a flexible internal representation of the spatial layout of the environment (expressed in allocentric terms). This “cognitive map” is thought to be supported by the hippocampal formation, where there are neurons tuned to place and heading direction (17–19). Spatial navigation using this map is flexible because it can be used with arbitrary starting locations and destinations, which need not be marked by immediate sensory cues.We posit that the distinction between place and response learning is analogous to that between MB and MF RL (20). Under this view, associative reinforcement is supported by the DLS (21, 22). Indeed, there is evidence from both rodents (23–25) and humans (26, 27) that spatial-response learning relies on the same basal ganglia structures that support MF RL. Evidence also suggests an analogy between MB reasoning and hippocampus (HPC)-based place learning (28, 29). However, this equivalence is not completely straightforward. For example, in rodents, multiple hippocampal lesion and inactivation studies failed to elicit an effect on action-outcome learning, a hallmark of MB planning (30–35). Nevertheless, there are indications that HPC might contribute to a different aspect of MB RL: namely, the representation of relational structure. Tasks that require memory of the relationships between stimuli do show dependence on HPC (36–42).Here, we formalize the perspective that hippocampal contributions to MB learning and place learning are the same, as are the dorsolateral striatal contributions to MF and response learning. In our model, HPC supports flexible behavior by representing the relational structure among different allocentric states, while dorsolateral striatum (DLS) supports associative reinforcement over egocentric sensory features. The model arbitrates between the use of these systems by weighting each system’s action values by the reliability of the system, as measured by a recent average of prediction errors, following Wan Lee et al. (43). We show that HPC and DLS maintain these roles across multiple task domains, including a range of spatial and nonspatial tasks. Our model can quantitatively explain a range of seemingly disparate findings, including the choice between place and response strategies in spatial navigation (23, 44) and choices on nonspatial multistep decision tasks (45, 46). Furthermore, it explains the puzzling finding that landmark-guided navigation is sensitive to the blocking effect, whereas boundary-guided navigation is not (27), and that these are supported by the DLS and HPC, respectively (26). Thus, different RL strategies that manage competing tradeoffs can explain a longstanding body of spatial navigation and decision-making literature under a unified model.     

ICAM-1 induced rearrangements of capsid and genome prime rhinovirus 14 for activation and uncoating

Most rhinoviruses, which are the leading cause of the common cold, utilize intercellular adhesion molecule-1 (ICAM-1) as a receptor to infect cells. To release their genomes, rhinoviruses convert to activated particles that contain pores in the capsid, lack minor capsid protein VP4, and have an altered genome organization. The binding of rhinoviruses to ICAM-1 promotes virus activation; however, the molecular details of the process remain unknown. Here, we present the structures of virion of rhinovirus 14 and its complex with ICAM-1 determined to resolutions of 2.6 and 2.4 Å, respectively. The cryo-electron microscopy reconstruction of rhinovirus 14 virions contains the resolved density of octanucleotide segments from the RNA genome that interact with VP2 subunits. We show that the binding of ICAM-1 to rhinovirus 14 is required to prime the virus for activation and genome release at acidic pH. Formation of the rhinovirus 14–ICAM-1 complex induces conformational changes to the rhinovirus 14 capsid, including translocation of the C termini of VP4 subunits, which become poised for release through pores that open in the capsids of activated particles. VP4 subunits with altered conformation block the RNA–VP2 interactions and expose patches of positively charged residues. The conformational changes to the capsid induce the redistribution of the virus genome by altering the capsid–RNA interactions. The restructuring of the rhinovirus 14 capsid and genome prepares the virions for conversion to activated particles. The high-resolution structure of rhinovirus 14 in complex with ICAM-1 explains how the binding of uncoating receptors enables enterovirus genome release.

Human rhinoviruses are the cause of more than half of common colds (1). Medical visits and missed days of school and work cost tens of billions of US dollars annually (2, 3). There is currently no cure for rhinovirus infections, and the available treatments are only symptomatic. Rhinoviruses belong to the family Picornaviridae, genus Enterovirus, and are classified into species A, B, and C (4). Rhinoviruses A and B can belong to either “major” or “minor” groups, based on their utilization of intercellular adhesion molecule-1 (ICAM-1) or low-density lipoprotein receptor for cell entry (5–7). Type C rhinoviruses use CDHR3 as a receptor (8). Rhinovirus 14 belongs to the species rhinovirus B and uses ICAM-1 as a receptor. Receptors recognized by rhinoviruses and other enteroviruses can be divided into two groups based on their function in the infection process (9). Attachment receptors such as DAF, PSGL1, KREMEN1, CDHR3, and sialic acid enable the binding and endocytosis of virus particles into cells (10–13). In contrast, uncoating receptors including ICAM-1, CD155, CAR, and SCARB2 enable virus cell entry but also promote genome release from virus particles (5, 14–16).Virions of rhinoviruses are nonenveloped and have icosahedral capsids (17). Genomes of rhinoviruses are 7,000 to 9,000 nucleotide-long single-stranded positive-sense RNA molecules (1, 17). The rhinovirus genome encodes a single polyprotein that is co- and posttranslationally cleaved into functional protein subunits. Capsid proteins VP1, VP3, and VP0, originating from one polyprotein, form a protomer, 60 of which assemble into a pseudo-T = 3 icosahedral capsid. To render the virions mature and infectious, VP0 subunits are cleaved into VP2 and VP4 (18, 19). VP1 subunits form pentamers around fivefold symmetry axes, whereas subunits VP2 and VP3 form heterohexamers centered on threefold symmetry axes. The major capsid proteins VP1 through 3 have a jelly roll β-sandwich fold formed by two β-sheets, each containing four antiparallel β-strands, which are conventionally named B to I (20–22). The two β-sheets contain the strands BIDG and CHEF, respectively. The C termini of the capsid proteins are located at the virion surface, whereas the N termini mediate interactions between the capsid proteins and the RNA genome on the inner surface of the capsid. VP4 subunits are attached to the inner face of the capsid formed by the major capsid proteins. The surfaces of rhinovirus virions are characterized by circular depressions called canyons, which are centered around fivefold symmetry axes of the capsids (21).The VP1 subunits of most rhinoviruses, but not those of rhinovirus 14, contain hydrophobic pockets, which are filled by molecules called pocket factors (17, 21, 23, 24). It has been speculated that pocket factors are fatty acids or lipids (25). The pockets are positioned immediately below the canyons. The exposure of rhinoviruses to acidic pH induces expulsion of the pocket factors, which leads to the formation of activated particles and genome release (17, 26–32). The activated particles are characterized by capsid expansion, a reduction in interpentamer contacts, the release of VP4 subunits, externalization of N termini of VP1 subunits, and changes in the distribution of RNA genomes (17, 26–29, 33, 34). Artificial hydrophobic compounds that bind to VP1 pockets with high affinity inhibit infection by rhinoviruses (35, 36).ICAM-1 is an endothelial- and leukocyte-associated protein that stabilizes cell–cell interactions and facilitates the movement of leukocytes through endothelia (37). ICAM-1 can be divided into an extracellular amino-terminal part composed of five immunoglobulin domains, a single transmembrane helix, and a 29-residue–long carboxyl-terminal cytoplasmic domain. The immunoglobulin domains are characterized by a specific fold that consists of seven to eight β-strands, which form two antiparallel β-sheets in a sandwich arrangement (38–40). The immunoglobulin domains of ICAM-1 are stabilized by disulfide bonds and glycosylation (38–41). The connections between the immunoglobulin domains are formed by flexible linkers that enable bending of the extracellular part of ICAM-1. For example, the angle between domains 1 and 2 differs by 8° between molecules in distinct crystal forms (38, 42). As a virus receptor, ICAM-1 enables the virus particles to sequester at the cell surface and mediates their endocytosis (5). The structures of complexes of rhinoviruses 3, 14, and 16, and CVA21 with ICAM-1 have been determined to resolutions of 9 to 28 Å (42–46). It was shown that ICAM-1 molecules bind into the canyons at the rhinovirus surface, approximately between fivefold and twofold symmetry axes (42–46). ICAM-1 interacts with residues from all three major capsid proteins. It has been speculated that the binding of ICAM-1 triggers the transition of virions of rhinovirus 14 to activated particles and initiates genome release (45, 47). However, the limited resolution of the previous studies prevented characterization of the corresponding molecular mechanism.Here, we present the cryo-electron microscopy (cryo-EM) reconstruction of the rhinovirus 14 virion, which contains resolved density of octanucleotide segments of the RNA genome that interact with VP2 subunits. Furthermore, we show that the binding of ICAM-1 to rhinovirus 14 induces changes in its capsid and genome, which are required for subsequent virus activation and genome release at acidic pH.     

Improved bounds on entropy production in living systems

Living systems maintain or increase local order by working against the second law of thermodynamics. Thermodynamic consistency is restored as they consume free energy, thereby increasing the net entropy of their environment. Recently introduced estimators for the entropy production rate have provided major insights into the efficiency of important cellular processes. In experiments, however, many degrees of freedom typically remain hidden to the observer, and, in these cases, existing methods are not optimal. Here, by reformulating the problem within an optimization framework, we are able to infer improved bounds on the rate of entropy production from partial measurements of biological systems. Our approach yields provably optimal estimates given certain measurable transition statistics. In contrast to prevailing methods, the improved estimator reveals nonzero entropy production rates even when nonequilibrium processes appear time symmetric and therefore may pretend to obey detailed balance. We demonstrate the broad applicability of this framework by providing improved bounds on the energy consumption rates in a diverse range of biological systems including bacterial flagella motors, growing microtubules, and calcium oscillations within human embryonic kidney cells.

Thermodynamic laws place fundamental limits on the efficiency and fitness of living systems (1, 2). To maintain cellular order and perform essential biological functions such as sensing (3–6), signaling (7), replication (8, 9) or locomotion (10), organisms consume energy and dissipate heat. In doing so, they increase the entropy of their environment (2), in agreement with the second law of thermodynamics (11). Obtaining reliable estimates for the energy consumption and entropy production in living matter holds the key to understanding the physical boundaries (12–14) that constrain the range of theoretically and practically possible biological processes (3). Recent experimental (6, 15, 16) and theoretical (17–20) advances in the imaging and modeling of cellular and subcellular dynamics have provided groundbreaking insights into the thermodynamic efficiency of molecular motors (14, 21), biochemical signaling (16, 22, 23) and reaction (24) networks, and replication (9) and adaption (25) phenomena. Despite such major progress, however, it is also known that the currently available entropy production estimators (26, 27) can fail under experimentally relevant conditions, especially when only a small set of observables is experimentally accessible or nonequilibrium transport currents (28–30) vanish.To help overcome these limitations, we introduce here a generic optimization framework that can produce significantly improved bounds on the entropy production in living systems. We will prove that these bounds are optimal given certain measurable statistics. From a practical perspective, our method requires observations of only a few coarse-grained state variables of an otherwise hidden Markovian network. We demonstrate the practical usefulness by determining improved entropy production bounds for bacterial flagella motors (10, 31), growing microtubules (32, 33), and calcium oscillations (7, 34) in human embryonic kidney cells.Generally, entropy production rates can be estimated from the time series of stochastic observables (35). Thermal equilibrium systems obey the principle of detailed balance, which means that every forward trajectory is as likely to be observed as its time reversed counterpart, neutralizing the arrow of time (36). By contrast, living organisms operate far from equilibrium, which means that the balance between forward and reversed trajectories is broken and net fluxes may arise (1, 37–39). When all microscopic details of a nonequilibrium system are known, one can measure the rate of entropy production by comparing the likelihoods of forward and reversed trajectories in sufficiently large data samples (35, 36). However, in most if not all biophysical experiments, many degrees of freedom remain hidden to the observer, demanding methods (28, 40, 41) that do not require complete knowledge of the system. A powerful alternative is provided by thermodynamic uncertainty relations (TURs), which use the mean and variance of steady-state currents to bound entropy production rates (18, 19, 26, 42–48). Although highly useful when currents can be measured (44–47), or when the system can be externally manipulated (40, 49), these methods give, by construction, trivial zero bounds for current-free nonequilibrium systems, such as driven one-dimensional (1D) nonperiodic oscillators. In the absence of currents, potential asymmetries in the forward and reverse trajectories can still be exploited to bound the entropy production rate (29, 30, 50), but to our knowledge no existing method is capable of producing nonzero bounds when forward and reverse trajectories are statistically identical. Moreover, even though previous bounds can become tight in some cases (51), optimal entropy production estimators for nonequilibrium systems are in general unknown.To obtain bounds that are provably optimal under reasonable conditions on the available data, we reformulate the problem here within an optimization framework. Formally, we consider any steady-state Markovian dynamics for which only coarse-grained variables are observable, where these observables may appear non-Markovian. We then search over all possible underlying Markovian systems to identify the one which minimizes the entropy production rate while obeying the observed statistics. More specifically, our algorithmic implementation leverages information about successive transitions, allowing us to discover nonzero bounds on entropy production even when the coarse-grained statistics appear time symmetric. We demonstrate this for both synthetic test data and experimental data (52) for flagella motors. Subsequently, we consider the entropy production of microtubules (33), which slowly grow before rapidly shrinking in steady state, to show how refined coarse graining in space and time leads to improved bounds. The final application to calcium oscillations in human embryonic kidney cells (34) illustrates how external stimulation with drugs can increase entropy production.