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.     

Heading perception depends on time-varying evolution of optic flow

There is considerable support for the hypothesis that perception of heading in the presence of rotation is mediated by instantaneous optic flow. This hypothesis, however, has never been tested. We introduce a method, termed “nonvarying phase motion,” for generating a stimulus that conveys a single instantaneous optic flow field, even though the stimulus is presented for an extended period of time. In this experiment, observers viewed stimulus videos and performed a forced-choice heading discrimination task. For nonvarying phase motion, observers made large errors in heading judgments. This suggests that instantaneous optic flow is insufficient for heading perception in the presence of rotation. These errors were mostly eliminated when the velocity of phase motion was varied over time to convey the evolving sequence of optic flow fields corresponding to a particular heading. This demonstrates that heading perception in the presence of rotation relies on the time-varying evolution of optic flow. We hypothesize that the visual system accurately computes heading, despite rotation, based on optic acceleration, the temporal derivative of optic flow.

James Gibson first remarked that the instantaneous motion of points on the retina (Fig. 1A) can be formally described as a two-dimensional (2D) field of velocity vectors called the “optic flow field” (or “optic flow”) (1). Such optic flow, caused by an observer’s movement relative to the environment, conveys information about self-motion and the structure of the visual scene (1–15). When an observer translates in a given direction along a straight path, the optic flow field radiates from a point in the image with zero velocity, or singularity, called the focus of expansion (Fig. 1B). It is well known that under such conditions, one can accurately estimate one’s “heading” (i.e., instantaneous direction of translation in retinocentric coordinates) by simply locating the focus of expansion (SI Appendix). However, if there is angular rotation in addition to translation (by moving along a curved path or by a head or eye movement), the singularity in the optic flow field will be displaced such that it no longer corresponds to the true heading (Fig. 1 C and D). In this case, if one estimates heading by locating the singularity, the estimate will be biased away from the true heading. This is known as the rotation problem (14).Open in a separate windowFig. 1.Projective geometry, the rotation problem, time-varying optic flow, and the optic acceleration hypothesis. (A) Viewer-centered coordinate frame and perspective projection. Because of motion between the viewpoint and the scene, a 3D surface point traverses a path in 3D space. Under perspective projection, the 3D path of this point projects onto a 2D path in the image plane (retina), the temporal derivative of which is called image velocity. The 2D velocities associated with all visible points define a dense 2D vector field called the optic flow field. (B–D) Illustration of the rotation problem. (B) Optic flow for pure translation (1.5-m/s translation speed, 0° heading, i.e., heading in the direction of gaze). Optic flow singularity (red circle) corresponds to heading (purple circle). (C) Pure rotation, for illustrative purposes only and not corresponding to any experimental condition (2°/s rightward rotation). (D) Translation + rotation (1.5 m/s translation speed, 0° heading, 2°/s rightward rotation). Optic flow singularity (red circle) is displaced away from heading (purple circle). (E) Three frames from a video depicting movement along a circular path with the line-of-sight initially perpendicular to a single fronto-parallel plane composed of black dots. (F) Time-varying evolution of optic flow. The first optic flow field reflects image motion between the first and second frames of the video. The second optic flow field reflects image motion between the second and third frames of the video. For this special case (circular path), the optic flow field evolves (and the optic flow singularity drifts) only due to the changing depth of the environment relative to the viewpoint. (G) Illustration of the optic acceleration hypothesis. Optic acceleration is the derivative of optic flow over time (here, approximated as the difference between the second and first optic flow fields). The singularity of the optic acceleration field corresponds to the heading direction. Acceleration vectors autoscaled for visibility.Computer vision researchers and vision scientists have developed a variety of algorithms that accurately and precisely extract observer translation and rotation from optic flow, thereby solving the rotation problem. Nearly all of these rely on instantaneous optic flow (i.e., a single optic flow field) (4, 9, 16–25) with few exceptions (26–29). However, it is unknown whether these algorithms are commensurate with the neural computations underlying heading perception.The consensus of opinion in the experimental literature is that human observers can estimate heading (30, 31) from instantaneous optic flow, in the absence of additional information (5, 10, 15, 32–34). Even so, there are reports of systematic biases in heading perception (11); the visual consequences of rotation (eye, head, and body) can bias heading judgments (10, 15, 35–37), with the amount of bias typically proportional to the magnitude of rotation. Other visual factors, such as stereo cues (38, 39), depth structure (8, 10, 40–43), and field of view (FOV) (33, 42–44) can modulate the strength of these biases. Errors in heading judgments have been reported to be greater when eye (35–37, 45, 46) or head movements (37) are simulated versus when they are real, which has been taken to mean that observers require extraretinal information, although there is also evidence to the contrary (10, 15, 33, 40, 41, 44, 47–50). Regardless, to date no one has tested whether heading perception (even with these biases) is based on instantaneous optic flow or on the information available in how the optic flow field evolves over time. Some have suggested that heading estimates rely on information accumulated over time (32, 44, 51), but no one has investigated the role of time-varying optic flow without confounding it with stimulus duration (i.e., the duration of evidence accumulation).In this study, we employed an application of an image processing technique that ensured that only a single optic flow field was available to observers, even though the stimulus was presented for an extended period of time. We called this condition “nonvarying phase motion” or “nonvarying”: The phases of two component gratings comprising each stationary stimulus patch shifted over time at a constant rate, causing a percept of motion in the absence of veridical movement (52). Phase motion also eliminated other cues that may otherwise have been used for heading judgments, including image point trajectories (15, 32) and their spatial compositions (i.e., looming) (53, 54). For nonvarying phase motion, observers exhibited large biases in heading judgments in the presence of rotation. A second condition, “time-varying phase motion,” or “time-varying,” included acceleration by varying the velocity of phase motion over time to match the evolution of a sequence of optic flow fields. Doing so allowed observers to compensate for the confounding effect of rotation on optic flow, making heading perception nearly veridical. This demonstrates that heading perception in the presence of rotation relies on the time-varying evolution of optic flow.     

The high-affinity immunoglobulin receptor FcγRI potentiates HIV-1 neutralization via antibodies against the gp41 N-heptad repeat

The HIV-1 gp41 N-heptad repeat (NHR) region of the prehairpin intermediate, which is transiently exposed during HIV-1 viral membrane fusion, is a validated clinical target in humans and is inhibited by the Food and Drug Administration (FDA)-approved drug enfuvirtide. However, vaccine candidates targeting the NHR have yielded only modest neutralization activities in animals; this inhibition has been largely restricted to tier-1 viruses, which are most sensitive to neutralization by sera from HIV-1–infected individuals. Here, we show that the neutralization activity of the well-characterized NHR-targeting antibody D5 is potentiated >5,000-fold in TZM-bl cells expressing FcγRI compared with those without, resulting in neutralization of many tier-2 viruses (which are less susceptible to neutralization by sera from HIV-1–infected individuals and are the target of current antibody-based vaccine efforts). Further, antisera from guinea pigs immunized with the NHR-based vaccine candidate (ccIZN36)₃ neutralized tier-2 viruses from multiple clades in an FcγRI-dependent manner. As FcγRI is expressed on macrophages and dendritic cells, which are present at mucosal surfaces and are implicated in the early establishment of HIV-1 infection following sexual transmission, these results may be important in the development of a prophylactic HIV-1 vaccine.

Membrane fusion between HIV-1 and host cells is mediated by the viral envelope glycoprotein (Env), a trimer consisting of the gp120 and gp41 subunits. Upon interaction with cellular receptors, Env undergoes a dramatic conformational change and forms the prehairpin intermediate (PHI) (1–3), in which the fusion peptide region at the amino terminus of gp41 inserts into the cell membrane. In the PHI, the N-heptad repeat (NHR) region of gp41 is exposed and forms a stable, three-stranded α-helical coiled coil. Subsequently, the PHI resolves when the NHR and the C-heptad repeat (CHR) regions of gp41 associate to form a trimer-of-hairpins structure that brings the viral and cell membranes into proximity, facilitating membrane fusion (Fig. 1).Open in a separate windowFig. 1.HIV-1 membrane fusion. The surface protein of the HIV-1 envelope is composed of the gp120 and gp41 subunits. After Env binds to cell-surface receptors, gp41 inserts into the host cell membrane and undergoes a conformational change to form the prehairpin intermediate. The N-heptad repeat (orange) region of gp41 is exposed in the PHI and forms a three-stranded coiled coil. To complete viral fusion, the PHI resolves to a trimer-of-hairpins structure in which the C-heptad repeat (blue) adopts a helical conformation and binds the NHR region. Fusion inhibitors such as enfuvirtide bind the NHR, preventing viral fusion by inhibiting formation of the trimer of hairpins (1–3). The membrane-proximal external region (red) is located adjacent to the transmembrane (TM) region of gp41.The NHR region of the PHI is a validated therapeutic target in humans: the Food and Drug Administration (FDA)-approved drug enfuvirtide binds the NHR and inhibits viral entry into cells (4, 5). Various versions of the three-stranded coiled coil formed by the NHR have been created and used as vaccine candidates in animals (6–10). The neutralization potencies of these antisera, as well as those of anti-NHR monoclonal antibodies (mAbs) (11–15), are modest and mostly limited to HIV-1 isolates that are highly sensitive to antibody-mediated neutralization [commonly referred to as tier-1 viruses (16)]. These results have led to skepticism about the PHI as a vaccine target.Earlier studies showed that the neutralization activities of mAbs that bound another region of gp41, the membrane-proximal external region (MPER) (Fig. 1), were enhanced as much as 5,000-fold in cells expressing FcγRI (CD64) (17, 18), an integral membrane protein that binds the Fc portion of immunoglobulin G (IgG) molecules with high (nanomolar) affinity (19, 20). This effect was not attributed to phagocytosis and occurred when the cells were preincubated with antibody and washed before adding virus (17, 18). Since the MPER is a partially cryptic epitope that is not fully exposed until after Env engages with cellular receptors (21, 22), these results suggest that by binding the Fc region, FcγRI provides a local concentration advantage for MPER mAbs at the cell surface that enhances viral neutralization (17, 18). While not expressed on T cells, FcγRI is expressed on macrophages and dendritic cells (23), which are present at mucosal surfaces and are implicated in sexual HIV-1 transmission and the early establishment of HIV-1 infection (22–34).Here we investigated whether FcγRI expression also potentiates the neutralizing activity of antibodies targeting the NHR, since that region, like the MPER, is preferentially exposed during viral fusion. We found that D5, a well-characterized anti-NHR mAb (11, 12), inhibits HIV-1 infection ∼5,000-fold more potently in TZM-bl cells expressing FcγRI (TZM-bl/FcγRI cells) than in TZM-bl cells that do not. Further, while antisera from guinea pigs immunized with (ccIZN36)₃, an NHR-based vaccine candidate (7), displayed weak neutralizing activity in TZM-bl cells, they exhibited enhanced neutralization in TZM-bl/FcγRI cells, including against some tier-2 HIV-1 isolates that are more resistant to antibody-mediated neutralization (16) and that serve as benchmarks for antibody-based vaccine efforts. These results indicate that FcγRI can play an important role in neutralization by antibodies that target the PHI. Since these receptors are expressed on cells prevalent at mucosal surfaces thought to be important for sexual HIV-1 transmission, our results motivate vaccine strategies that harness this potentiating effect.     

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.     

Structure of Ycf1p reveals the transmembrane domain TMD0 and the regulatory region of ABCC transporters

ATP binding cassette (ABC) proteins typically function in active transport of solutes across membranes. The ABC core structure is composed of two transmembrane domains (TMD1 and TMD2) and two cytosolic nucleotide binding domains (NBD1 and NBD2). Some members of the C-subfamily of ABC (ABCC) proteins, including human multidrug resistance proteins (MRPs), also possess an N-terminal transmembrane domain (TMD0) that contains five transmembrane α-helices and is connected to the ABC core by the L0 linker. While TMD0 was resolved in SUR1, the atypical ABCC protein that is part of the hetero-octameric ATP-sensitive K⁺ channel, little is known about the structure of TMD0 in monomeric ABC transporters. Here, we present the structure of yeast cadmium factor 1 protein (Ycf1p), a homolog of human MRP1, determined by electron cryo-microscopy (cryo-EM). A comparison of Ycf1p, SUR1, and a structure of MRP1 that showed TMD0 at low resolution demonstrates that TMD0 can adopt different orientations relative to the ABC core, including a ∼145° rotation between Ycf1p and SUR1. The cryo-EM map also reveals that segments of the regulatory (R) region, which links NBD1 to TMD2 and was poorly resolved in earlier ABCC structures, interacts with the L0 linker, NBD1, and TMD2. These interactions, combined with fluorescence quenching experiments of isolated NBD1 with and without the R region, suggest how posttranslational modifications of the R region modulate ABC protein activity. Mapping known mutations from MRP2 and MRP6 onto the Ycf1p structure explains how mutations involving TMD0 and the R region of these proteins lead to disease.

ATP binding cassette (ABC) proteins are a large family of membrane proteins found in all kingdoms of life (1, 2). ABC proteins have a core structure composed of two transmembrane (TM) domains (TMD1 and TMD2) and two cytosolic nucleotide binding domains (NBD1 and NBD2) (Fig. 1A and SI Appendix, Fig. S1A) (3–5). Through ATP binding and hydrolysis at the NBDs, ABC proteins actively transport solutes across cell membranes, regulate activities of other proteins, or function as channels (1, 2). Thus, ABC proteins are involved in many biological processes, including lipid homeostasis, cellular metal trafficking, and antigen peptide transport. Mutations in human ABC proteins cause diseases such as Tangier disease, adenoleukodystrophy, cystic fibrosis, Dubin–Johnson syndrome, and pseudoxanthoma elasticum (PXE) (1, 2). Furthermore, the export of a wide range of cancer chemotherapeutics, antibiotics, and anti-fungal drugs by ABC transporters confers multidrug resistance to tumor cells, bacteria, and fungal pathogens, respectively (1, 2, 6, 7).Open in a separate windowFig. 1.Ycf1p structure. (A) Ycf1p domain arrangement. TMD, transmembrane domain; L0, L0 linker; NBD, nucleotide binding domain; and R, regulatory (R) region. The Pep4p proteolytic digestion site within the luminal loop 6 of TMD1 is denoted by a pink “*.” Phosphorylation sites in the L0 linker (S251) and R region (S908 and T911) are depicted with a “P” circled in red. (B) Cryo-EM density of Ycf1p with domains colored as in A. (C) Example of the atomic model for individual TM helices in TMD0 and the R region fit into the corresponding map densities. (D) Schematic ribbon diagram of Ycf1p colored as in A and B and with the proteolytic digestion site denoted by a pink “*.”Human ABC proteins are divided into seven subfamilies (A to G) based in part on the sequence of their NBDs and TMDs in the core ABC structure (1, 2). The C-subfamily is the most diverse and includes the cystic fibrosis transmembrane conductance regulator (CFTR), the sulphonylurea receptors that form regulatory subunits in ATP-sensitive K⁺ (K_ATP) channels, and the multidrug resistance proteins (MRPs). In addition to the ABC core, ABCC proteins contain an N-terminal extension that is either composed of an additional TM domain (TMD0) and L0 linker (Fig. 1A, orange and tan, respectively, and SI Appendix, Fig. S1A) or just an L0 tail (5, 8). A TMD0, but not L0 linker, is also found in some ABCB proteins (3, 5). These N-terminal extensions are involved in trafficking, endosomal recycling, protein interactions, and/or regulation of ABC activity (9–18). The existence of disease-causing mutations in TMD0 and the L0 linker of different ABCC proteins (8, 13, 18) indicates that these regions play important roles in protein function when present.High-resolution structural information for TMD0 is available only for the atypical ABCC protein SUR1 (19, 20), which is part of the large hetero-octameric K_ATP channel complex. In contrast, structures of monomeric ABC transporters showed only low-resolution density for TMD0 that was insufficient for building a full atomic model or lacked density for the domain altogether (14, 21–24). The vacuolar ABCC protein yeast cadmium factor 1 (Ycf1p) from Saccharomyces cerevisiae is a close homolog of human MRPs and a model for ABCC proteins that function as monomers. Ycf1p transports glutathione-conjugated heavy metals, such as Cd²⁺, from the cytosol into the vacuole, detoxifying the cell (25, 26). Human MRP1 can rescue Cd²⁺ transport activity in a YCF1 deletion strain (27).Like other ABCC proteins, Ycf1p contains a relatively long and mostly disordered linker that connects NBD1 and TMD2 (25, 28, 29) (Fig. 1A and SI Appendix, Fig. S1A). This linker contains stimulatory phosphorylation sites (25, 28), similar to the phospho-regulatory (R) region in the ABCC protein CFTR (30–32). Ycf1p also contains an inhibitory phosphorylation site in the L0 linker (33). However, how the R region interacts with the ABCC core and how its phosphorylation modulates protein function remain poorly understood for most ABCC proteins. Structural studies of Ycf1p presented here reveal how TMD0 and the R region exert their regulatory functions in MRP-like ABCC proteins.     

A rate threshold mechanism regulates MAPK stress signaling and survival

Cells are exposed to changes in extracellular stimulus concentration that vary as a function of rate. However, how cells integrate information conveyed from stimulation rate along with concentration remains poorly understood. Here, we examined how varying the rate of stress application alters budding yeast mitogen-activated protein kinase (MAPK) signaling and cell behavior at the single-cell level. We show that signaling depends on a rate threshold that operates in conjunction with stimulus concentration to determine the timing of MAPK signaling during rate-varying stimulus treatments. We also discovered that the stimulation rate threshold and stimulation rate-dependent cell survival are sensitive to changes in the expression levels of the Ptp2 phosphatase, but not of another phosphatase that similarly regulates osmostress signaling during switch-like treatments. Our results demonstrate that stimulation rate is a regulated determinant of cell behavior and provide a paradigm to guide the dissection of major stimulation rate dependent mechanisms in other systems.

All cells employ signal transduction pathways to respond to physiologically relevant changes in extracellular stressors, nutrient levels, hormones, morphogens, and other stimuli that vary as functions of both concentration and rate in healthy and diseased states (1–7). Switch-like “instantaneous” changes in the concentrations of stimuli in the extracellular environment have been widely used to show that the strength of signaling and overall cellular response are dependent on the stimulus concentration, which in many cases needs to exceed a certain threshold (8, 9). Previous studies have shown that the rate of stimulation can also influence signaling output in a variety of pathways (10–17) and that stimulation profiles of varying rates can be used to probe underlying signaling pathway circuitry (4, 18, 19). However, it is still not clear how cells integrate information conveyed by changes in both the stimulation rate and concentration in determining signaling output. It is also not clear if cells require stimulation gradients to exceed a certain rate in order to commence signaling.Recent investigations have demonstrated that stimulation rate can be a determining factor in signal transduction. In contrast to switch-like perturbations, which trigger a broad set of stress-response pathways, slow stimulation rates activate a specific response to the stress applied in Bacillus subtilis cells (10). Meanwhile, shallow morphogen gradient stimulation fails to activate developmental pathways in mouse myoblast cells in culture, even when concentrations sufficient for activation during pulsed treatment are delivered (12). These observations raise the possibility that stimulation profiles must exceed a set minimum rate or rate threshold to achieve signaling activation. Although such rate thresholds would help cells decide if and how to respond to dynamic changes in stimulus concentration, the possibility of signaling regulation by a rate threshold has never been directly investigated in any system. Further, no study has experimentally examined how stimulation rate requirements impact cell phenotype or how cells molecularly regulate the stimulation rate required for signaling activation. As such, the biological significance of any existing rate threshold regulation of signaling remains unknown.The budding yeast Saccharomyces cerevisiae high osmolarity glycerol (HOG) pathway provides an ideal model system for addressing these issues (Fig. 1A). The evolutionarily conserved mitogen-activated protein kinase (MAPK) Hog1 serves as the central signaling mediator of this pathway (20–22). It is well established that instantaneous increases in osmotic stress concentration induce Hog1 phosphorylation, activation, and translocation to the nucleus (18, 21, 23–30). Activated Hog1 governs the majority of the cellular osmoadaptation response that enables cells to survive (23, 31, 32). Multiple apparently redundant MAPK phosphatases dephosphorylate and inactivate Hog1, which, along with the termination of upstream signaling after adaptation, results in its return to the cytosol (Fig. 1A) (23, 25, 26, 33–39). Because of this behavior, time-lapse analysis of Hog1 nuclear enrichment in single cells has proven an excellent and sensitive way to monitor signaling responses to dynamic stimulation patterns in real time (18, 27–30, 40, 41). Further, such assays have been readily combined with traditional growth and molecular genetic approaches to link observed signaling responses with cell behavior and signaling pathway architecture (27–29).Open in a separate windowFig. 1.Hog1 signaling and cell survival are sensitive to the rate of preconditioning osmotic stress application. (A) Schematic of the budding yeast HOG response. (B) Preconditioning protection assay workflow indicating the first stress treatments to a final concentration of 0.4 M NaCl (Left), high-stress exposure (Middle), and colony formation readout (Right). (C) High-stress survival as a function of each first treatment relative to the untreated first stress condition. Bars and errors are means and SD from three biological replicates. *Statistically significant by Kolmogorov–Smirnov test (P < 0.05). NS = not significant. (D) Treatment concentration over time. (E) Treatment rate over time for quadratic and pulse treatment. The rate for the pulse is briefly infinite (blue vertical line) before it drops to 0. (F) Hog1 nuclear localization during the treatments depicted in D and E. (Inset) Localization pattern in the quadratic-treated sample. Lines represent means and shaded error represents the SD from three to four biological replicates.Here, we use systematically designed osmotic stress treatments imposed at varying rates of increase to show that a rate threshold condition regulates yeast high-stress survival and Hog1 MAPK signaling. We demonstrate that only stimulus profiles that satisfy both this rate threshold condition and a concentration threshold condition result in robust signaling. We go on to show that the protein tyrosine phosphatase Ptp2, but not the related Ptp3 phosphatase, serves as a major rate threshold regulator. By expressing PTP2 under the control of a series of different enhancer–promoter DNA constructs, we demonstrate that changes in the level of Ptp2 expression can alter the stimulation rate required for signaling induction and survival. These findings establish rate thresholds as a critical and regulated component of signaling biology akin to concentration thresholds.     

Geometrical manipulation of complex supramolecular tessellations by hierarchical assembly of amphiphilic platinum(II) complexes

Here we report complex supramolecular tessellations achieved by the directed self-assembly of amphiphilic platinum(II) complexes. Despite the twofold symmetry, these geometrically simple molecules exhibit complicated structural hierarchy in a columnar manner. A possible key to such an order increase is the topological transition into circular trimers, which are noncovalently interlocked by metal···metal and π–π interactions, thereby allowing for cofacial stacking in a prismatic assembly. Another key to success is to use the immiscibility of the tailored hydrophobic and hydrophilic sidechains. Their phase separation leads to the formation of columnar crystalline nanostructures homogeneously oriented on the substrate, featuring an unusual geometry analogous to a rhombitrihexagonal Archimedean tiling. Furthermore, symmetry lowering of regular motifs by design results in an orthorhombic lattice obtained by the coassembly of two different platinum(II) amphiphiles. These findings illustrate the potentials of supramolecular engineering in creating complex self-assembled architectures of soft materials.

Tessellation in two dimensions (2D) is a very old topic in geometry on how one or more shapes can be periodically arranged to fill a Euclidean plane without any gaps. Tessellation principles have been extensively applied in decorative art since the early times. In natural sciences, there has been a growing attention on creating ordered structures with increasingly complex architectures inspired by semi-regular Archimedean tilings (ATs) and quasicrystalline textures on account of their intriguing physical properties (1–5) and biological functions (6). Recent advances in this regard have been achieved in various fields of supramolecular science, including the programmable self-assembly of DNA molecules (7), coordination-driven assembly (8–10), supramolecular interfacial engineering (11–13), crystallization of organic polygons (14, 15), colloidal particle superlattices (16), and other soft-matter systems (17–20). Moreover, tessellation in 2D can overcome the topological frustration to generate complex semi- or non-regular patterns by using geometrically simple motifs. As exemplified by the self-templating assembly of spherical soft microparticles (21), a vast array of 2D micropatterns encoding non-regular tilings, such as rectangular, rhomboidal, hexagonal, and herringbone superlattices were obtained by layer-by-layer strategy at a liquid–liquid interface. Tessellation principles have also been extended to the self-assembly of giant molecules in three dimensions (3D). Superlattices with high space-group symmetry (Im3¯m, Pm3¯n, and P4₂/mnm) were reported in dendrimers and dendritic polymers by Percec and coworkers (22–24). Recently, Cheng and coworkers identified the highly ordered Frank–Kasper phases obtained from giant amphiphiles containing molecular nanoparticles (25–28). Despite such advancements made in the field of soft matter, an understanding of how structural ordering in supramolecular materials is influenced by the geometric factors of its constituent molecules has so far remained elusive.In light of these developments and the desire to explore the supramolecular systems, square-planar platinum(II) (Pt^II) polypyridine complexes may serve as an ideal candidate for model studies not only because of their intriguing spectroscopic and luminescence properties (29, 30), but also because of their propensity to form supramolecular polymers or oligomers via noncovalent Pt···Pt and π–π interactions (31–39). Although rod-shaped and lamellar structures are the most commonly observed in the self-assembly of planar Pt^II complexes (34–39), 2D-ordered nanostructures, such as the hexagonally packed columns (31, 40) and honeycomb-like networks (41–43), were recently first demonstrated by us.Herein, we report a serendipitous discovery of a C_2h-symmetric Pt^II amphiphile (Fig. 1A) that can hierarchically self-assemble into a 3D-ordered nanostructure with hexagonal geometry. Interestingly, this structurally anisotropic molecule possibly undergoes topological transition and interlocks to form its circular trimer by noncovalent Pt···Pt and π–π interactions (Fig. 1B). The resultant triangular motif is architecturally stabilized and preorganized for one-dimensional (1D) prismatic assembly (Fig. 1C). Together with the phase separation of the tailored hydrophobic and hydrophilic sidechains, an unusual and unique 3D hexagonal lattice is formed (Fig. 1D), in which the Pt centers adopt a rare rhombitrihexagonal AT-like order. Finally, the nanoarchitecture develops in a hierarchical manner on the substrate due to the homogeneous nucleation (Fig. 1E).Open in a separate windowFig. 1.Hierarchical self-assembly of Pt^II amphiphile into hexagonal ordering. (A) Space-filling (CPK) model of a C_2h-symmetric Pt^II amphiphile (1). All of the hydrogen atoms and counterions are omitted for clarity. (B) CPK representations of possible models of regular triangular, tetragonal, pentagonal, and hexagonal motifs formed with Pt···Pt and π–π stacking. These motifs possess a hydrophilic core (red) with various diameters wrapped by a hydrophobic shell comprising long alkyl chains (gray). (C) CPK representation of a 1D prismatic structure consisting of circular trimers with long-range Pt···Pt and π–π stacking. (D) CPK representation of a 3D columnar lattice constructed by the prismatic assemblies adopting a rare rhombitrihexagonal AT-like order. With the assistance of the phase separation, the hydrophobic domain serves as a discrete column associated with six prismatic neighbors. (E) Schematic representation of the nanoarchitecture with homogeneous orientation.     

A solution to a sex ratio puzzle in Melittobia wasps

The puzzling sex ratio behavior of Melittobia wasps has long posed one of the greatest questions in the field of sex allocation. Laboratory experiments have found that, in contrast to the predictions of theory and the behavior of numerous other organisms, Melittobia females do not produce fewer female-biased offspring sex ratios when more females lay eggs on a patch. We solve this puzzle by showing that, in nature, females of Melittobia australica have a sophisticated sex ratio behavior, in which their strategy also depends on whether they have dispersed from the patch where they emerged. When females have not dispersed, they lay eggs with close relatives, which keeps local mate competition high even with multiple females, and therefore, they are selected to produce consistently female-biased sex ratios. Laboratory experiments mimic these conditions. In contrast, when females disperse, they interact with nonrelatives, and thus adjust their sex ratio depending on the number of females laying eggs. Consequently, females appear to use dispersal status as an indirect cue of relatedness and whether they should adjust their sex ratio in response to the number of females laying eggs on the patch.

Sex allocation has produced many of the greatest success stories in the study of social behaviors (1–4). Time and time again, relatively simple theory has explained variation in how individuals allocate resources to male and female reproduction. Hamilton’s local mate competition (LMC) theory predicts that when n diploid females lay eggs on a patch and the offspring mate before the females disperse, the evolutionary stable proportion of male offspring (sex ratio) is (n − 1)/2n (Fig. 1) (5). A female-biased sex ratio is favored to reduce competition between sons (brothers) for mates and to provide more mates (daughters) for those sons (6–8). Consistent with this prediction, females of >40 species produce female-biased sex ratios and reduce this female bias when multiple females lay eggs on the same patch (higher n; Fig. 1) (9). The fit of data to theory is so good that the sex ratio under LMC has been exploited as a “model trait” to study the factors that can constrain “perfect adaptation” (4, 10–13).Open in a separate windowFig. 1.LMC. The sex ratio (proportion of sons) is plotted versus the number of females laying eggs on a patch. The bright green dashed line shows the LMC theory prediction for the haplodiploid species (5, 39). A more female-biased sex ratio is favored in haplodiploids because inbreeding increases the relative relatedness of mothers to their daughters (7, 32). Females of many species adjust their offspring sex ratio as predicted by theory, such as the parasitoid Nasonia vitripennis (green diamonds) (82). In contrast, the females of several Melittobia species, such as M. australica, continue to produce extremely female-biased sex ratios, irrespective of the number of females laying eggs on a patch (blue squares) (15).In stark contrast, the sex ratio behavior of Melittobia wasps has long been seen as one of the greatest problems for the field of sex allocation (3, 4, 14–21). The life cycle of Melittobia wasps matches the assumptions of Hamilton’s LMC theory (5, 15, 19, 21). Females lay eggs in the larvae or pupae of solitary wasps and bees, and then after emergence, female offspring mate with the short-winged males, who do not disperse. However, laboratory experiments on four Melittobia species have found that females lay extremely female-biased sex ratios (1 to 5% males) and that these extremely female-biased sex ratios change little with increasing number of females laying eggs on a patch (higher n; Fig. 1) (15, 17–20, 22). A number of hypotheses to explain this lack of sex ratio adjustment have been investigated and rejected, including sex ratio distorters, sex differential mortality, asymmetrical male competition, and reciprocal cooperation (15–18, 20, 22–26).We tested whether Melittobia’s unusual sex ratio behavior can be explained by females being related to the other females laying eggs on the same patch. After mating, some females disperse to find new patches, while some may stay at the natal patch to lay eggs on previously unexploited hosts (Fig. 2). If females do not disperse, they can be related to the other females laying eggs on the same host (27–31). If females laying eggs on a host are related, this increases the extent to which relatives are competing for mates and so can favor an even more female-biased sex ratio (28, 32–35). Although most parasitoid species appear unable to directly assess relatedness, dispersal behavior could provide an indirect cue of whether females are with close relatives (36–38). Consequently, we predict that when females do not disperse and so are more likely to be with closer relatives, they should maintain extremely female-biased sex ratios, even when multiple females lay eggs on a patch (28, 35).Open in a separate windowFig. 2.Host nest and dispersal manners of Melittobia. (A) Photograph of the prepupae of the leaf-cutter bee C. sculpturalis nested in a bamboo cane and (B) a diagram showing two ways that Melittobia females find new hosts. The mothers of C. sculpturalis build nursing nests with pine resin consisting of individual cells in which their offspring develop. If Melittobia wasps parasitize a host in a cell, female offspring that mate with males inside the cell find a different host on the same patch (bamboo cane) or disperse by flying to other patches.We tested whether the sex ratio of Melittobia australica can be explained by dispersal status in a natural population. We examined how the sex ratio produced by females varies with the number of females laying eggs on a patch and whether or not they have dispersed before laying eggs. To match our data to the predictions of theory, we developed a mathematical model tailored to the unique population structure of Melittobia, where dispersal can be a cue of relatedness. We then conducted a laboratory experiment to test whether Melittobia females are able to directly access the relatedness to other females and adjust their sex ratio behavior accordingly. Our results suggest that females are adjusting their sex ratio in response to both the number of females laying eggs on a patch and their relatedness to the other females. However, relatedness is assessed indirectly by whether or not they have dispersed. Consequently, the solution to the puzzling behavior reflects a more-refined sex ratio strategy.     

Structure,self-assembly,and properties of a truncated reflectin variant

Naturally occurring and recombinant protein-based materials are frequently employed for the study of fundamental biological processes and are often leveraged for applications in areas as diverse as electronics, optics, bioengineering, medicine, and even fashion. Within this context, unique structural proteins known as reflectins have recently attracted substantial attention due to their key roles in the fascinating color-changing capabilities of cephalopods and their technological potential as biophotonic and bioelectronic materials. However, progress toward understanding reflectins has been hindered by their atypical aromatic and charged residue-enriched sequences, extreme sensitivities to subtle changes in environmental conditions, and well-known propensities for aggregation. Herein, we elucidate the structure of a reflectin variant at the molecular level, demonstrate a straightforward mechanical agitation-based methodology for controlling this variant’s hierarchical assembly, and establish a direct correlation between the protein’s structural characteristics and intrinsic optical properties. Altogether, our findings address multiple challenges associated with the development of reflectins as materials, furnish molecular-level insight into the mechanistic underpinnings of cephalopod skin cells’ color-changing functionalities, and may inform new research directions across biochemistry, cellular biology, bioengineering, and optics.

Materials from naturally occurring and recombinant proteins are frequently employed for the study of fundamental biological processes and leveraged for applications in fields as diverse as electronics, optics, bioengineering, medicine, and fashion (1–13). Such broad utility is enabled by the numerous advantageous characteristics of protein-based materials, which include sequence modularity, controllable self-assembly, stimuli-responsiveness, straightforward processability, inherent biological compatibility, and customizable functionality (1–13). Within this context, unique structural proteins known as reflectins have recently attracted substantial attention because of their key roles in the fascinating color-changing capabilities of cephalopods, such as the squid shown in Fig. 1A, and have furthermore demonstrated their utility for unconventional biophotonic and bioelectronic technologies (11–40). For example, in vivo, Bragg stack-like ultrastructures from reflectin-based high refractive index lamellae (membrane-enclosed platelets) are responsible for the angle-dependent narrowband reflectance (iridescence) of squid iridophores, as shown in Fig. 1B (15–20). Analogously, folded membranes containing distributed reflectin-based particle arrangements within sheath cells lead to the mechanically actuated iridescence of squid chromatophore organs, as shown in Fig. 1C (15, 16, 21, 22). Moreover, in vitro, films processed from squid reflectins not only exhibit proton conductivities on par with some state-of-the-art artificial materials (23–27) but also support the growth of murine and human neural stem cells (28, 29). Additionally, morphologically variable coatings assembled from different reflectin isoforms can enable the functionality of chemically and electrically actuated color-changing devices, dynamic near-infrared camouflage platforms, and stimuli-responsive photonic architectures (27, 30–34). When considered together, these discoveries and demonstrations constitute compelling motivation for the continued exploration of reflectins as model biomaterials.Open in a separate windowFig. 1.(A) A camera image of a D. pealeii squid for which the skin contains light-reflecting cells called iridophores (bright spots) and pigmented organs called chromatophores (colored spots). Image credit: Roger T. Hanlon (photographer). (B) An illustration of an iridophore (Left), which shows internal Bragg stack-like ultrastructures from reflectin-based lamellae (i.e., membrane-enclosed platelets) (Inset). (C) An illustration of a chromatophore organ (Left), which shows arrangements of reflectin-based particles within the sheath cells (Inset). (D) The logo of the 28-residue-long N-terminal motif (RMN), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (E) The logo of the 28-residue-long internal motif (RMI), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (F) The logo of the 21-residue-long C-terminal motif (RMC), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (G) The amino acid sequence of full-length D. pealeii reflectin A1, which contains a single RMN motif (gray oval) and five RMI motifs (orange ovals). (H) An illustration of the selection of the prototypical truncated reflectin variant (denoted as RfA1TV) from full-length D. pealeii reflectin A1.Given reflectins’ demonstrated significance from both fundamental biology and applications perspectives, some research effort has been devoted to resolving their three-dimensional (3D) structures (30, 31, 35–39). For example, fibers drawn from full-length Euprymna scolopes reflectin 1a and films processed from truncated E. scolopes reflectin 1a were shown to possess secondary structural elements (i.e., α-helices or β-sheets) (30, 31). In addition, precipitated nanoparticles and drop-cast films from full-length Doryteuthis pealeii reflectin A1 have exhibited β-character, which was seemingly associated with their conserved motifs (35, 36). Moreover, nanoparticles assembled from both full-length and truncated Sepia officinalis reflectin 2 variants have demonstrated signatures consistent with β-sheet or α-helical secondary structure, albeit in the presence of surfactants (38). However, such studies were made exceedingly challenging by reflectins’ atypical primary sequences enriched in aromatic and charged residues, documented extreme sensitivities to subtle changes in environmental conditions, and well-known propensities for poorly controlled aggregation (12, 14, 15, 30–32, 34–39). Consequently, the reported efforts have all suffered from multiple drawbacks, including the need for organic solvents or denaturants, the evaluation of only polydisperse or aggregated (rather than monomeric) proteins, a lack of consensus among different experimental techniques, inadequate resolution that precluded molecular-level insight, imperfect agreement between computational predictions and experimental observations, and/or the absence of conclusive correlations between structure and optical functionality. As such, there has emerged an exciting opportunity for investigating reflectins’ molecular structures, which remain poorly understood and the subject of some debate.Herein, we elucidate the structure of a reflectin variant at the molecular level, demonstrate a robust methodology for controlling this variant’s hierarchical assembly, and establish a direct correlation between its structural characteristics and optical properties. We first rationally select a prototypical reflectin variant expected to recapitulate the behavior of its parent protein by using a bioinformatics-guided approach. We next map the conformational and energetic landscape accessible to our selected protein by means of all-atom molecular dynamics (MD) simulations. We in turn produce our truncated reflectin variant with and without isotopic labeling, develop solution conditions that maintain the protein in a monomeric state, and characterize the variant’s size and shape with small-angle X-ray scattering (SAXS). We subsequently resolve our protein’s dynamic secondary and tertiary structures and evaluate its backbone conformational fluctuations with NMR spectroscopy. Finally, we demonstrate a straightforward mechanical agitation-based approach to controlling our truncated reflectin variant’s secondary structure, hierarchical self-assembly, and bulk refractive index distribution. Overall, our findings address multiple challenges associated with the development of reflectins as materials, furnish molecular-level insight into the mechanistic underpinnings of cephalopod skin cells’ color-changing functionalities, and appear poised to inform new directions across biochemistry, cellular biology, bioengineering, and optics.     

Developmental influence on evolutionary rates and the origin of placental mammal tooth complexity

Development has often been viewed as a constraining force on morphological adaptation, but its precise influence, especially on evolutionary rates, is poorly understood. Placental mammals provide a classic example of adaptive radiation, but the debate around rate and drivers of early placental evolution remains contentious. A hallmark of early dental evolution in many placental lineages was a transition from a triangular upper molar to a more complex upper molar with a rectangular cusp pattern better specialized for crushing. To examine how development influenced this transition, we simulated dental evolution on “landscapes” built from different parameters of a computational model of tooth morphogenesis. Among the parameters examined, we find that increases in the number of enamel knots, the developmental precursors of the tooth cusps, were primarily influenced by increased self-regulation of the molecular activator (activation), whereas the pattern of knots resulted from changes in both activation and biases in tooth bud growth. In simulations, increased activation facilitated accelerated evolutionary increases in knot number, creating a lateral knot arrangement that evolved at least ten times on placental upper molars. Relatively small increases in activation, superimposed on an ancestral tritubercular molar growth pattern, could recreate key changes leading to a rectangular upper molar cusp pattern. Tinkering with tooth bud geometry varied the way cusps initiated along the posterolingual molar margin, suggesting that small spatial variations in ancestral molar growth may have influenced how placental lineages acquired a hypocone cusp. We suggest that development could have enabled relatively fast higher-level divergence of the placental molar dentition.

Whether developmental processes bias or constrain morphological adaptation is a long-standing question in evolutionary biology (1–4). Many of the distinctive features of a species derive from pattern formation processes that establish the position and number of anatomical structures (5). If developmental processes like pattern formation are biased toward generating only particular kinds of variation, adaptive radiations may often be directed along developmental–genetic “lines of least resistance” (2, 4, 6, 7). Generally, the evolutionary consequences of this developmental bias have been considered largely in terms of how it might influence the pattern of character evolution (e.g., refs. 1, 2, 8–10). But development could also influence evolutionary rates by controlling how much variation is accessible to natural selection in a given generation (11).For mammals, the dentition is often the only morphological system linking living and extinct species (12). Correspondingly, tooth morphology plays a crucial role in elucidating evolutionary relationships, time calibrating phylogenetic trees, and reconstructing adaptive responses to past environmental change (e.g., refs. 13–15). One of the most pervasive features of dental evolution among mammals is an increase in the complexity of the tooth occlusal surface, primarily through the addition of new tooth cusps (16, 17). These increases in tooth complexity are functionally and ecologically significant because they enable more efficient mechanical breakdown of lower-quality foods like plant leaves (18).Placental mammals are the most diverse extant mammalian group, comprising more than 6,000 living species spread across 19 extant orders, and this taxonomic diversity is reflected in their range of tooth shapes and dietary ecologies (12). Many extant placental orders, especially those with omnivorous or herbivorous ecologies (e.g., artiodactyls, proboscideans, rodents, and primates), convergently evolved a rectangular upper molar cusp pattern from a placental ancestor with a more triangular cusp pattern (19–21). This resulted from separate additions in each lineage of a novel posterolingual cusp, the "hypocone'''' [sensu (19)], to the tritubercular upper molar (Fig. 1), either through modification of a posterolingual cingulum (“true” hypocone) or another posterolingual structure, like a metaconule (pseudohypocone) (19). The fossil record suggests that many of the basic steps in the origin of this rectangular cusp pattern occurred during an enigmatic early diversification window associated with the divergence and early radiation of several placental orders (20, 21; Fig. 1). However, there remains debate about the rate and pattern of early placental divergence (22–24). On the one hand, most molecular phylogenies suggest that higher-level placental divergence occurred largely during the Late Cretaceous (25, 26), whereas other molecular phylogenies and paleontological analyses suggest more rapid divergence near the Cretaceous–Paleogene (K–Pg) boundary (21, 24, 27–29). Most studies agree that ecological opportunity created in the aftermath of the K–Pg extinction probably played an important role in ecomorphological diversification within the placental orders (30, 31). But exactly how early placentals acquired the innovations needed to capitalize on ecological opportunity remains unclear. Dental innovations, especially those which facilitated increases in tooth complexity, may have been important because they would have promoted expansion into plant-based dietary ecologies left largely vacant after the K–Pg extinction event (32).Open in a separate windowFig. 1.Placental mammal lineages separately evolved complex upper molar teeth with a rectangular cusp pattern composed of two lateral pairs of cusps from a common ancestor with a simpler, triangular cusp pattern. Many early relatives of the extant placental orders, such as Eritherium, possessed a hypocone cusp and a more rectangular primary cusp pattern. Examples of complex upper molars are the following: Proboscidea, the gomphothere Anancus; Rodentia, the wood mouse Apodemus; and Artiodactyla, the suid Nyanzachoerus.Mammalian tooth cusps form primarily during the “cap” and “bell” stage of dental development, when signaling centers called enamel knots establish the future sites of cusp formation within the inner dental epithelium (33, 34). The enamel knots secrete molecules that promote proliferation and changes in cell–cell adhesion, which facilitates invagination of the dental epithelium into an underlying layer of mesenchymal cells (34, 35). Although a range of genes are involved in tooth cusp patterning (36–38), the basic dynamics can be effectively modeled using reaction–diffusion models with just three diffusible morphogens: an activator, an inhibitor, and a growth factor (39–41). Candidate activator genes in mammalian tooth development include Bmp4, Activin A, Fgf20, and Wnt genes, whereas potential inhibitors include Shh and Sostdc, and Fgf4 and Bmp2 have been hypothesized to act as growth factors (38, 40–43). In computer models of tooth development, activator molecules up-regulated in the underlying mesenchyme stimulate differentiation of overlying epithelium into nondividing enamel knot cells. These in turn secrete molecules that inhibit further differentiation of epithelium into knot cells, while also promoting cell proliferation that creates the topographic relief of the cusp (40). Although many molecular, cellular, and physical processes have the potential to influence cusp formation, and thereby tooth complexity (35, 37), parameters that control the strength and conductance of the activator and inhibitor signals, the core components of the reaction–diffusion cusp patterning mechanism (39, 40) are likely to be especially important.Here, we integrate a previous computer model of tooth morphogenesis called ToothMaker (41), with simulations of trait evolution and data from the fossil record (Fig. 2), to examine the developmental origins of tooth complexity in placental mammals. Specifically, we ask the following: 1) What developmental processes can influence how many cusps form? 2) How might these developmental processes influence the evolution of tooth cusp number, especially rates? And 3) what developmental changes may have been important in the origins of the fourth upper molar cusp, the hypocone, in placental mammal evolution?Open in a separate windowFig. 2.Workflow for simulations of tooth complexity evolution. (A) Tooth shape is varied for five signaling and growth parameters in ToothMaker. (B) From an ancestral state, each parameter is varied in 2.5% increments up to a maximum of ± 50% of the ancestral state. (C) Tooth complexity and enamel knot (EK) pattern were quantified for each parameter combination. Tooth complexity was measured using cusp number/EK number and OPC. ToothMaker and placental upper second molars were classified into categories based on EK/cusp pattern. (D) The parameter space was populated with pattern and tooth complexity datums to build a developmental landscape. (E) Tooth complexity evolution was simulated on each developmental landscape. (F) Resulting diversity and pattern of tooth complexity was compared with placental mammal molar diversity.     

Stable metal anodes enabled by a labile organic molecule bonded to a reduced graphene oxide aerogel

Metallic anodes (lithium, sodium, and zinc) are attractive for rechargeable battery technologies but are plagued by an unfavorable metal–electrolyte interface that leads to nonuniform metal deposition and an unstable solid–electrolyte interphase (SEI). Here we report the use of electrochemically labile molecules to regulate the electrochemical interface and guide even lithium deposition and a stable SEI. The molecule, benzenesulfonyl fluoride, was bonded to the surface of a reduced graphene oxide aerogel. During metal deposition, this labile molecule not only generates a metal-coordinating benzenesulfonate anion that guides homogeneous metal deposition but also contributes lithium fluoride to the SEI to improve Li surface passivation. Consequently, high-efficiency lithium deposition with a low nucleation overpotential was achieved at a high current density of 6.0 mA cm⁻². A Li|LiCoO₂ cell had a capacity retention of 85.3% after 400 cycles, and the cell also tolerated low-temperature (−10 °C) operation without additional capacity fading. This strategy was applied to sodium and zinc anodes as well.

Rechargeable batteries based on metal anodes including lithium (Li), sodium (Na), and zinc (Zn) show great promise in achieving high energy density (1–3). Unfortunately, the electrochemical interface of the metal anodes is not favorable for metal deposition. Metal nucleation is inhomogeneous at the surface, leading to the growth of metal dendrites (4–7) and the formation of an unstable solid–electrolyte interphase (SEI) that is incapable of protecting metals from the side reactions with the electrolyte (8–12).Substantial efforts have been devoted to stabilizing the interface of metal anodes, especially for Li metal. These include the design of artificial protective layers (13–17), alternative electrolytes (18–24), and sacrificial additives (25–30) to stabilize the metal–electrolyte interface, the development of mechanically robust coatings (31–34) to block Li dendrite growth, and the use of structured scaffolds to host dendrite-free Li deposition by reducing local current densities (35–43). However, the performance of metal anodes remains poor under high-current or low-temperature conditions. This is because the inhomogeneous Li nucleation and unstable SEI problems have not been well addressed, and these problems at the interface are even exacerbated under critical operating conditions, especially high-current densities and low temperatures (5, 6, 44).Toward this end, we report a simple molecular approach for regulating the electrochemical interface of metal anodes, which enables even Li deposition and stable SEI formation in a conventional electrolyte. This was realized by bonding a labile organic molecule, benzenesulfonyl fluoride (BSF), to a reduced graphene oxide (rGO) aerogel surface as the Li anode host (Fig. 1A). During Li deposition, BSF molecules electrochemically decompose at the interface and generate benzenesulfonate anions bonded to the rGO aerogel (Fig. 1B). The conjugated anions have a strong binding affinity for Li, serving as lithiophilic sites on the rGO surface to synergistically induce homogeneous Li nucleation of Li on the rGO surface. At the same time, BSF molecules contribute LiF to the SEI layer, which facilitates Li surface passivation (Fig. 1C). As a result, high-efficiency (99.2%) Li deposition was achieved at a Li deposition amount of 6.0 mAh cm⁻² and a current density of 6.0 mA cm⁻²; the barrier to Li nucleation was markedly reduced, as evidenced by the low nucleation overpotentials at high-current density (6.0 mA cm⁻²) or at a low temperature (−10 °C). A 400-cycle life with a capacity retention of 83.6% was achieved for a Li|LiCoO₂ (LCO) cell in a conventional carbonate electrolyte. Moreover, with the organic molecule-tuned interface, the Li|LCO cell can be stably cycled at a low operating temperature (−10 °C). This approach was applied to Na and Zn metal anodes as well.Open in a separate windowFig. 1.Illustration of a stable interface for Li deposition using a labile organic molecule, benzenesulfonyl fluoride (BSF). (A) Covalently bonded BSF on the rGO aerogel surface. (B) In situ generation of a lithiophilic conjugated anion (benzenesulfonate) and LiF on the surface during Li deposition. (C) Li nucleation preferentially occurs at the conjugated anion sites owing to the strong Li binding affinity, which leads to uniform Li deposition. In addition, the LiF that is formed is in the SEI layer and passivates the Li surface.     

Type III secretion system effector proteins are mechanically labile

Multiple gram-negative bacteria encode type III secretion systems (T3SS) that allow them to inject effector proteins directly into host cells to facilitate colonization. To be secreted, effector proteins must be at least partially unfolded to pass through the narrow needle-like channel (diameter <2 nm) of the T3SS. Fusion of effector proteins to tightly packed proteins—such as GFP, ubiquitin, or dihydrofolate reductase (DHFR)—impairs secretion and results in obstruction of the T3SS. Prior observation that unfolding can become rate-limiting for secretion has led to the model that T3SS effector proteins have low thermodynamic stability, facilitating their secretion. Here, we first show that the unfolding free energy (ΔGunfold0) of two Salmonella effector proteins, SptP and SopE2, are 6.9 and 6.0 kcal/mol, respectively, typical for globular proteins and similar to published ΔGunfold0 for GFP, ubiquitin, and DHFR. Next, we mechanically unfolded individual SptP and SopE2 molecules by atomic force microscopy (AFM)-based force spectroscopy. SptP and SopE2 unfolded at low force (F_unfold ≤ 17 pN at 100 nm/s), making them among the most mechanically labile proteins studied to date by AFM. Moreover, their mechanical compliance is large, as measured by the distance to the transition state (Δx^‡ = 1.6 and 1.5 nm for SptP and SopE2, respectively). In contrast, prior measurements of GFP, ubiquitin, and DHFR show them to be mechanically robust (F_unfold > 80 pN) and brittle (Δx^‡ < 0.4 nm). These results suggest that effector protein unfolding by T3SS is a mechanical process and that mechanical lability facilitates efficient effector protein secretion.

Type III secretion systems (T3SS) are large nanomachines utilized by both pathogenic and symbiotic bacteria to inject effector proteins directly into the cytoplasm of host cells (1–3). Once delivered, effector proteins facilitate host cell colonization through a variety of mechanisms (4–7), including down-regulation of the host immune response (8) and rearrangement of the cytoskeleton (9, 10). The T3SS apparatus, known as the injectisome, is a syringe-like structure with a hollow needle that spans the inner and outer bacterial membranes, the extracellular space, and the host membrane, enabling proteins to pass directly from bacteria to host cells (Fig. 1A) (2). Specialized bacterial chaperones often bind the N-terminal 50 to 100 amino acids (aa) of the effector proteins, known as the chaperone binding domain, and help maintain the effector N-terminal domain in an extended conformation. C-terminal to the chaperone binding domain, effector proteins contain one or more globular domains, which adopt their folded conformations even when in complex with their cognate chaperone (4, 11, 12). The effector proteins, or their chaperone complexes, are recognized by the base of the injectisome prior to secretion (13). At its narrowest point, the injectisome needle’s inner diameter is less than 2 nm (14–16). As a result, effector proteins must be mostly unfolded to be secreted (17–20). Secretion is thus thought to proceed by a “threading-the-needle mechanism,” where the N-terminal extended domain is released from the chaperone and fed to the injectisome, followed by unfolding of the C-terminal effector domain (21).Open in a separate windowFig. 1.Thermodynamic stability of T3SS effector proteins SptP_CD and SopE2_CD. (A) Schematic depiction of protein transport through the T3SS showing effector proteins, which are at least partially folded in the bacterial cytoplasm. Such effector proteins interact with an associated unfoldase to passage through the T3SS, which has an inner channel with a diameter <2 nm. Once inside the host cytoplasm, effector proteins refold to carry out their function. (B) Crystal structures of SptP_CD (Protein Data Bank [PDB] ID code 1G4U) and SopE2_CD (PDB ID code 1R9K). (C) Ellipticity from CD at λ = 222 nm plotted as a function of urea concentrations for SptP_CD (orange) and SopE2_CD (green). A fit of the data with Eq. 1 yielded the free energy of unfolding ΔGunfold0 for SptP_CD (6.9 ± 0.2 kcal/mol [mean ± fit error]) and SopE2_CD (6.0 ± 0.2 kcal/mol [mean ± fit error]). Data points are the result of at least three independent measurements. Error bars represent SD.Before proteins are secreted through the T3SS, they interact with a hexameric ATPase at the base of the T3SS that is capable of mediating chaperone release from effector proteins and effector-protein unfolding (15, 22). Indeed, most in vivo unfolding is catalyzed by unfoldases that work from one end of the substrate protein in stark contrast to the global effects of temperature, pH, or chemical denaturants. The most common examples of targeted protein unfolding are catalyzed by ATPases of the AAA(+) family that mechanically unfold their substrates (23, 24). For example, the AAA(+) ATPase ClpX forms a ring-shaped hexamer that mechanically pulls its substrates through its narrow central pore to unfold them (25). These are powerful unfoldases that can unfold even tightly packed proteins such as GFP, ubiquitin, and dihydrofolate reductase (DHFR) (23, 24, 26, 27). However, the T3SS ATPase does not belong to the AAA(+) family of ATPases. Instead, it is structurally similar to the catalytic β-subunit of the F₁F₀ ATP synthase, a rotary motor that normally couples proton gradient dissipation to ATP synthesis but can also run in reverse and hydrolyze ATP to do work (15, 28–30). The T3SS ATPase is not as powerful an unfoldase as the AAA(+) family, as fusions of effector proteins with GFP, ubiquitin, or DHFR stall in the injectisome and are poorly secreted (20, 22, 31, 32). These observations have led to the current model that T3SS effector proteins have low thermodynamic stability to facilitate their secretion (22, 31–33).While thermodynamic stability is the most common metric of protein stability, mechanical stability is a distinct metric that quantifies how easily a protein unfolds under force (F_unfold). Mechanical stability is typically measured by pulling across the N and C termini of single molecules via force spectroscopy using optical tweezers (34, 35) or an atomic force microscope (AFM) (36). Early force spectroscopy studies showed that thermodynamic stability does not correlate with mechanical stability (37–41). For example, titin’s I28 domain requires ∼20% more force to unfold than titin’s I27 domain [I85 and I91, respectively, in the new nomenclature (42)], despite I27 having approximately twofold higher thermodynamic stability (43). Importantly, AFM studies have shown that GFP (44), ubiquitin (45), and DHFR (46) are mechanically robust, requiring high forces to unfold despite their typical thermodynamic stabilities. These three proteins each stall the T3SS; thus, mechanical stability may be the physical determinant to proteins being secreted by the T3SS, rather than thermodynamic stability.Here, we determine the thermodynamic and mechanical stabilities of SptP and SopE2, two effector proteins from Salmonella enterica. These effectors are ideal candidates for this study as they have known crystal structures (10, 47), have characterized in vivo secretion kinetics (48), and represent effector proteins of different size and structure (Fig. 1B). We show that the catalytic domains of SptP and SopE2 have unremarkable thermodynamic stabilities, similar to many other previously characterized proteins, including GFP, ubiquitin, and DHFR. Conversely, our AFM-based force spectroscopy measurements demonstrate that SptP and SopE2 are among the most mechanically labile proteins studied to date by AFM. These two T3SS effector proteins are therefore mechanically labile while being thermodynamically stable, supporting the hypothesis that it is mechanical stability, not thermodynamic stability, that predicts efficient protein secretion by the T3SS.