Anti-CRISPR Phages Cooperate to Overcome CRISPR-Cas Immunity.

Some phages encode anti-CRISPR (acr) genes, which antagonize bacterial CRISPR-Cas immune systems by binding components of its machinery, but it is less clear how deployment of these acr genes impacts phage replication and epidemiology. Here, we demonstrate that bacteria with CRISPR-Cas resistance are still partially immune to Acr-encoding phage. As a consequence, Acr-phages often need to cooperate in order to overcome CRISPR resistance, with a first phage blocking the host CRISPR-Cas immune system to allow a second Acr-phage to successfully replicate. This cooperation leads to epidemiological tipping points in which the initial density of Acr-phage tips the balance from phage extinction to a phage epidemic. Furthermore, both higher levels of CRISPR-Cas immunity and weaker Acr activities shift the tipping points toward higher initial phage densities. Collectively, these data help elucidate how interactions between phage-encoded immune suppressors and the CRISPR systems they target shape bacteria-phage population dynamics.

Six distinct NFκB signaling codons convey discrete information to distinguish stimuli and enable appropriate macrophage responses.

Macrophages initiate inflammatory responses via the transcription factor NFκB. The temporal pattern of NFκB activity determines which genes are expressed and thus, the type of response that ensues. Here, we examined how information about the stimulus is encoded in the dynamics of NFκB activity. We generated an mVenus-RelA reporter mouse line to enable high-throughput live-cell analysis of primary macrophages responding to host- and pathogen-derived stimuli. An information-theoretic workflow identified six dynamical features-termed signaling codons-that convey stimulus information to the nucleus. In particular, oscillatory trajectories were a hallmark of responses to cytokine but not pathogen-derived stimuli. Single-cell imaging and RNA sequencing of macrophages from a mouse model of Sjögren's syndrome revealed inappropriate responses to stimuli, suggestive of confusion of two NFκB signaling codons. Thus, the dynamics of NFκB signaling classify immune threats through six signaling codons, and signal confusion based on defective codon deployment may underlie the etiology of some inflammatory diseases.

Differential Expression of the Transcription Factor GATA3 Specifies Lineage and Functions of Innate Lymphoid Cells.

Lymphoid tissue inducer (LTi) cells are regarded as a subset of innate lymphoid cells (ILCs). However, these cells are not derived from the ILC common progenitor, which generates other ILC subsets and is defined by the expression of the transcription factor PLZF. Here, we examined transcription factor(s) determining the fate of LTi progenitors versus non-LTi ILC progenitors. Conditional deletion of Gata3 resulted in the loss of PLZF+ non-LTi progenitors but not the LTi progenitors that expressed the transcription factor RORγt. Consistently, PLZF+ non-LTi progenitors expressed high amounts of GATA3, whereas GATA3 expression was low in RORγt+ LTi progenitors. The generation of both progenitors required the transcriptional regulator Id2, which defines the common helper-like innate lymphoid progenitor (ChILP), but not cytokine signaling. Nevertheless, low GATA3 expression was necessary for the generation of functionally mature LTi cells. Thus, differential expression of GATA3 determines the fates and functions of distinct ILC progenitors.

Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

A dynamical systems treatment of transcriptomic trajectories in hematopoiesis.

Inspired by Waddington's illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein exogenous signaling dynamics couple to the enormously complex signaling and transcriptional machinery of a cell to elicit qualitative transitions in its collective state. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the variations of a myriad of genes. Here, we present a mathematical formalism for rigorously evaluating, from a dynamical systems perspective, whether scRNA-seq trajectories display statistical signatures consistent with bifurcations and, as a case study, pinpoint regions of multistability along the neutrophil branch of hematopoeitic differentiation. Additionally, we leverage the geometric features of linear instability to identify the low-dimensional phase plane in gene expression space within which the multistability unfolds, highlighting novel genetic players that are crucial for neutrophil differentiation. Broadly, we show that a dynamical systems treatment of scRNA-seq data provides mechanistic insights into the high-dimensional processes of cellular differentiation, taking a step toward systematic construction of mathematical models for transcriptomic dynamics.

Diverse branching forms regulated by a core auxin transport mechanism in plants.

Diverse branching forms have evolved multiple times across the tree of life to facilitate resource acquisition and exchange with the environment. In the vascular plant group, the ancestral pattern of branching involves dichotomy of a parent shoot apex to form two new daughter apices. The molecular basis of axillary branching in Arabidopsis is well understood, but few regulators of dichotomous branching are known. Through analyses of dichotomous branching in the lycophyte, Selaginella kraussiana, we identify PIN-mediated auxin transport as an ancestral branch regulator of vascular plants. We show that short-range auxin transport out of the apices promotes dichotomy and that branch dominance is globally coordinated by long-range auxin transport. Uniquely in Selaginella, angle meristems initiate at each dichotomy, and these can develop into rhizophores or branching angle shoots. We show that long-range auxin transport and a transitory drop in PIN expression are involved in angle shoot development. We conclude that PIN-mediated auxin transport is an ancestral mechanism for vascular plant branching that was independently recruited into Selaginella angle shoot development and seed plant axillary branching during evolution.

Bifurcation instructed design of multistate machines.

We propose a design paradigm for multistate machines where transitions from one state to another are organized by bifurcations of multiple equilibria of the energy landscape describing the collective interactions of the machine components. This design paradigm is attractive since, near bifurcations, small variations in a few control parameters can result in large changes to the system's state providing an emergent lever mechanism. Further, the topological configuration of transitions between states near such bifurcations ensures robust operation, making the machine less sensitive to fabrication errors and noise. To design such machines, we develop and implement a new efficient algorithm that searches for interactions between the machine components that give rise to energy landscapes with these bifurcation structures. We demonstrate a proof of concept for this approach by designing magnetoelastic machines whose motions are primarily guided by their magnetic energy landscapes and show that by operating near bifurcations we can achieve multiple transition pathways between states. This proof of concept demonstration illustrates the power of this approach, which could be especially useful for soft robotics and at the microscale where typical macroscale designs are difficult to implement.

The Rapid-reaction Kinetics of an Electron-Bifurcating Flavoprotein, the crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase EtfAB:bcd.

We have investigated the kinetic behavior of the electron-bifurcating crotonyl-CoA-dependent NADH: ferredoxin oxidoreductase EtfAB:bcd from Megasphaera elsdenii. The overall behavior of the complex in both the reductive and oxidative half-reactions is consistent with that previously determined for the individual EtfAB and bcd components. This includes an uncrossing of the half-potentials of the bifurcating flavin of the EtfAB component in the course of ferredoxin-reducing catalysis, ionization of the bcd flavin semiquinone and the appearance of a charge transfer complex upon binding of the high potential acceptor crotonyl-CoA. The observed rapid-reaction rates of ferredoxin reduction are independent of [NADH], [crotonyl-CoA] and [ferredoxin], with an observed rate of ∼0.2 s-1, consistent with the observed steady-state kinetics. In enzyme-monitored turnover experiments, an approach to steady-state where the complex's flavins become reduced but no ferredoxin is generated is followed by a steady-state phase characterized by extensive ferredoxin reduction but little change in overall levels of flavin reduction. The approach to steady-state phase can be eliminated by prior reduction of the complex, in which case there is no lag in the onset of ferredoxin reduction; this is consistent with the et FAD needing to be reduced to the level of the (anionic) semiquinone for bifurcation and concomitant ferredoxin reduction to occur. Single-turnover experiments support this conclusion, with accumulation of the anionic semiquinone of the et FAD apparently required to prime the system for subsequent bifurcation and ferredoxin reduction.

Functionally redundant formate dehydrogenases enable formate-dependent growth in Methanococcus maripaludis.

Methanogens are essential for the complete remineralization of organic matter in anoxic environments. Most cultured methanogens are hydrogenotrophic, using H2 as an electron donor to reduce CO2 to CH4, but in the absence of H2 many can also use formate. Formate dehydrogenase (Fdh) is essential for formate oxidation, where it transfers electrons for the reduction of coenzyme F420 or to a flavin-based electron bifurcating reaction catalyzed by heterodisulfide reductase (Hdr), the terminal reaction of methanogenesis. Furthermore, methanogens that use formate encode at least two isoforms of Fdh in their genomes, but how these different isoforms participate in methanogenesis is unknown. Using Methanococcus maripaludis, we undertook a biochemical characterization of both Fdh isoforms involved in methanogenesis. Both Fdh1 and Fdh2 interacted with Hdr to catalyze the flavin-based electron bifurcating reaction, and both reduced F420 at similar rates. F420 reduction preceded flavin-based electron bifurcation activity for both enzymes. In a Δfdh1 mutant background, a suppressor mutation was required for Fdh2 activity. Genome sequencing revealed that this mutation resulted in the loss of a specific molybdopterin transferase (moeA), allowing for Fdh2-dependent growth, and the metal content of the proteins suggested that isoforms are dependent on either molybdenum or tungsten for activity. These data suggest that both isoforms of Fdh are functionally redundant, but their activities in vivo may be limited by gene regulation or metal availability under different growth conditions. Together these results expand our understanding of formate oxidation and the role of Fdh in methanogenesis.

Collective cell migration is spatiotemporally regulated during mammary epithelial bifurcation.

Branched epithelial networks are generated through an iterative process of elongation and bifurcation. We sought to understand bifurcation of the mammary epithelium. To visualize this process, we utilized three-dimensional (3D) organotypic culture and time-lapse confocal microscopy. We tracked cell migration during bifurcation and observed local reductions in cell speed at the nascent bifurcation cleft. This effect was proximity dependent, as individual cells approaching the cleft reduced speed, whereas cells exiting the cleft increased speed. As the cells slow down, they orient both migration and protrusions towards the nascent cleft, while cells in the adjacent branches orient towards the elongating tips. We next tested the hypothesis that TGF-ß signaling controls mammary branching by regulating cell migration. We first validated that addition of TGF-ß1 (TGFB1) protein increased cleft number, whereas inhibition of TGF-ß signaling reduced cleft number. Then, consistent with our hypothesis, we observed that pharmacological inhibition of TGF-ß1 signaling acutely decreased epithelial migration speed. Our data suggest a model for mammary epithelial bifurcation in which TGF-ß signaling regulates cell migration to determine the local sites of bifurcation and the global pattern of the tubular network.

Structures and Electron Transport Paths in the Four Families of Flavin-Based Electron Bifurcation Enzymes.

Oxidoreductases facilitating electron transfer between molecules are pivotal in metabolic pathways. Flavin-based electron bifurcation (FBEB), a recently discovered energy coupling mechanism in oxidoreductases, enables the reversible division of electron pairs into two acceptors, bridging exergonic and otherwise unfeasible endergonic reactions. This chapter explores the four distinct FBEB complex families and highlights a decade of structural insights into FBEB complexes. In this chapter, we discuss the architecture, electron transfer routes, and conformational changes across all FBEB families, revealing the structural foundation that facilitate these remarkable functions.

Superposition strategies: How and why White people say contradictory things about race.

Due to the centrality of race and racism in social, economic, and political life, coupled with the racially privileged position of White people, the assessment of White racial attitudes is an ongoing concern. There is a great deal of survey-based, quantitative work that demonstrates a compelling case of White attitudinal polarization-a grouping of authoritarian, racist attitudes versus another alliance of progressive, antiracist attitudes-an increasingly racialized culture war. However, other studies, largely qualitative and open-ended, demonstrate the heterogeneous, shifting, and hypocritical nature of White discourse about race. To resolve this paradox, I refrain from the assumption that White racial "attitudes" are essentially bifurcated, while I also refuse the contention that White people produce spontaneous narratives whole-cloth. Rather, I argue that with sustained attention to time, context, and triangulation, we can better understand how and why White people speak of People of Color in positive ways one moment and negative the next, marshaling both to defend, rationalize, or improve their racialized subject position. I argue that these contradictions are-à la Schrödinger's famous thought experiment-"superposition strategies." Both racist and antiracist attitudes are simultaneously alive and dead in the same individual or group. Contradictory White discourse helps maintain a sense of self-efficacy and coherent White racial identity within conflictual and politically supercharged social situations, as well as within racially unequal social structures.

An uncharacteristically low-potential flavin governs the energy landscape of electron bifurcation.

Electron bifurcation, an energy-conserving process utilized extensively throughout all domains of life, represents an elegant means of generating high-energy products from substrates with less reducing potential. The coordinated coupling of exergonic and endergonic reactions has been shown to operate over an electrochemical potential of ∼1.3 V through the activity of a unique flavin cofactor in the enzyme NADH-dependent ferredoxin-NADP+ oxidoreductase I. The inferred energy landscape has features unprecedented in biochemistry and presents novel energetic challenges, the most intriguing being a large thermodynamically uphill step for the first electron transfer of the bifurcation reaction. However, ambiguities in the energy landscape at the bifurcating site deriving from overlapping flavin spectral signatures have impeded a comprehensive understanding of the specific mechanistic contributions afforded by thermodynamic and kinetic factors. Here, we elucidate an uncharacteristically low two-electron potential of the bifurcating flavin, resolving the energetic challenge of the first bifurcation event.

Fin ray branching is defined by TRAP+ osteolytic tubules in zebrafish.

The shaping of bone structures relies on various cell types and signaling pathways. Here, we use the zebrafish bifurcating fin rays during regeneration to investigate bone patterning. We found that the regenerating fin rays form via two mineralization fronts that undergo an osteoblast-dependent fusion/stitching until the branchpoint, and that bifurcation is not simply the splitting of one unit into two. We identified tartrate-resistant acid phosphatase-positive osteolytic tubular structures at the branchpoints, hereafter named osteolytic tubules (OLTs). Chemical inhibition of their bone-resorbing activity strongly impairs ray bifurcation, indicating that OLTs counteract the stitching process. Furthermore, by testing different osteoactive compounds, we show that the position of the branchpoint depends on the balance between bone mineralization and resorption activities. Overall, these findings provide a unique perspective on fin ray formation and bifurcation, and reveal a key role for OLTs in defining the proximo-distal position of the branchpoint.

Specific CaMKIIs mediate convergent extension cell movements in early zebrafish development.

BACKGROUND: Noncanonical Wnts are morphogens that can elevate intracellular Ca2+, activate the Ca2+/calmodulin-dependent protein kinase, CaMKII, and promote cell movements during vertebrate gastrulation. RESULTS: Zebrafish express seven CaMKII genes during embryogenesis; two of these, camk2b1 and camk2g1, are necessary for convergent extension (CE) cell movements. CaMKII morphant phenotypes were observed as early as epiboly. At the 1-3 somite stage, neuroectoderm and paraxial cells remained unconverged in both morphants. Later, somites lacked their stereotypical shape and were wider, more closely spaced, and body gap angles increased. At 24hpf, somite compression and notochord undulation coincided with a shorter and broader body axis. A camk2b1 crispant was generated which phenocopied the camk2b1 morphant. The levels of cell proliferation, apoptosis and paraxial and neuroectodermal markers were unchanged in morphants. Hyperactivation of CaMKII during gastrulation by transient pharmacological intervention (thapsigargin) also caused CE defects. Mosaically expressed dominant-negative CaMKII recapitulated these phenotypes and showed significant midline bifurcation. Finally, the introduction of CaMKII partially rescued Wnt11 morphant phenotypes. CONCLUSIONS: Overall, these data support a model whereby cyclically activated CaMKII encoded from two genes enables cell migration during the process of CE.

Disruption of the creb3l1 gene causes defects in caudal fin regeneration and patterning in zebrafish Danio rerio.

BACKGROUND: The gene cAMP-Responsive Element Binding protein 3-like-1 (CREB3L1) has been implicated in bone development in mice, with CREB3L1 knock-out mice exhibiting fragile bones, and in humans, with CREB3L1 mutations linked to osteogenesis imperfecta. However, the mechanism through which Creb3l1 regulates bone development is not fully understood. RESULTS: To probe the role of Creb3l1 in organismal physiology, we used CRISPR-Cas9 genome editing to generate a Danio rerio (zebrafish) model of Creb3l1 deficiency. In contrast to mammalian phenotypes, the Creb3l1 deficient fish do not display abnormalities in osteogenesis, except for a decrease in the bifurcation pattern of caudal fin. Both, skeletal morphology and overall bone density appear normal in the mutant fish. However, the regeneration of caudal fin postamputation is significantly affected, with decreased overall regenerate and mineralized bone area. Moreover, the mutant fish exhibit a severe patterning defect during regeneration, with a significant decrease in bifurcation complexity of the fin rays and distalization of the bifurcation sites. Analysis of genes implicated in bone development showed aberrant patterning of shha and ptch2 in Creb3l1 deficient fish, linking Creb3l1 with Sonic Hedgehog signaling during fin regeneration. CONCLUSIONS: Our results uncover a novel role for Creb3l1 in regulating tissue growth and patterning during regeneration.

Modulation of SK Channels via Calcium Buffering Tunes Intrinsic Excitability of Parvalbumin Interneurons in Neuropathic Pain: A Computational and Experimental Investigation.

Parvalbumin-expressing interneurons (PVINs) play a crucial role within the dorsal horn of the spinal cord by preventing touch inputs from activating pain circuits. In both male and female mice, nerve injury decreases PVINs' output via mechanisms that are not fully understood. In this study, we show that PVINs from nerve-injured male mice change their firing pattern from tonic to adaptive. To examine the ionic mechanisms responsible for this decreased output, we used a reparametrized Hodgkin-Huxley type model of PVINs, which predicted (1) the firing pattern transition is because of an increased contribution of small conductance calcium-activated potassium (SK) channels, enabled by (2) impairment in intracellular calcium buffering systems. Analyzing the dynamics of the Hodgkin-Huxley type model further demonstrated that a generalized Hopf bifurcation differentiates the two types of state transitions observed in the transient firing of PVINs. Importantly, this predicted mechanism holds true when we embed the PVIN model within the neuronal circuit model of the spinal dorsal horn. To experimentally validate this hypothesized mechanism, we used pharmacological modulators of SK channels and demonstrated that (1) tonic firing PVINs from naive male mice become adaptive when exposed to an SK channel activator, and (2) adapting PVINs from nerve-injured male mice return to tonic firing on SK channel blockade. Our work provides important insights into the cellular mechanism underlying the decreased output of PVINs in the spinal dorsal horn after nerve injury and highlights potential pharmacological targets for new and effective treatment approaches to neuropathic pain.SIGNIFICANCE STATEMENT Parvalbumin-expressing interneurons (PVINs) exert crucial inhibitory control over Aß fiber-mediated nociceptive pathways at the spinal dorsal horn. The loss of their inhibitory tone leads to neuropathic symptoms, such as mechanical allodynia, via mechanisms that are not fully understood. This study identifies the reduced intrinsic excitability of PVINs as a potential cause for their decreased inhibitory output in nerve-injured condition. Combining computational and experimental approaches, we predict a calcium-dependent mechanism that modulates PVINs' electrical activity following nerve injury: a depletion of cytosolic calcium buffer allows for the rapid accumulation of intracellular calcium through the active membranes, which in turn potentiates SK channels and impedes spike generation. Our results therefore pinpoint SK channels as potential therapeutic targets for treating neuropathic symptoms.

Rapid-reaction kinetics of the bifurcating NAD+-dependent NADPH:ferredoxin oxidoreductase NfnI from Pyrococcus furiosus.

We have investigated the kinetics of NAD+-dependent NADPH:ferredoxin oxidoreductase (NfnI), a bifurcating transhydrogenase that takes two electron pairs from NADPH to reduce two ferredoxins and one NAD+ through successive bifurcation events. NADPH reduction takes place at the bifurcating FAD of NfnI's large subunit, with high-potential electrons transferred to the [2Fe-2S] cluster and S-FADH of the small subunit, ultimately on to NAD+; low-potential electrons are transferred to two [4Fe-4S] clusters of the large subunit and on to ferredoxin. Reduction of NfnI by NADPH goes to completion only at higher pH, with a limiting kred of 36 ± 1.6 s-1 and apparent KdNADPH of 5 ± 1.2 µM. Reduction of one of the [4Fe-4S] clusters of NfnI occurs within a second, indicating that in the absence of NAD+, the system can bifurcate and generate low-potential electrons without NAD+. When enzyme is reduced by NADPH in the absence of NAD+ but the presence of ferredoxin, up to three equivalents of ferredoxin become reduced, although the reaction is considerably slower than seen during steady-state turnover. Bifurcation appears to be limited by transfer of the first, high-potential electron into the high-potential pathway. Ferredoxin reduction without NAD+ demonstrates that electron bifurcation is an intrinsic property of the bifurcating FAD and is not dependent on the simultaneous presence of NAD+ and ferredoxin. The tight coupling between NAD+ and ferredoxin reduction observed under multiple-turnover conditions is instead simply due to the need to remove reducing equivalents from the high-potential electron pathway under multiple-turnover conditions.

Rapid-reaction kinetics of the butyryl-CoA dehydrogenase component of the electron-bifurcating crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase from Megasphaera elsdenii.

We have investigated the equilibrium properties and rapid-reaction kinetics of the isolated butyryl-CoA dehydrogenase (bcd) component of the electron-bifurcating crotonyl-CoA-dependent NADH:ferredoxin oxidoreductase (EtfAB-bcd) from Megasphaera elsdenii. We find that a neutral FADH• semiquinone accumulates transiently during both reduction with sodium dithionite and with NADH in the presence of catalytic concentrations of EtfAB. In both cases full reduction of bcd to the hydroquinone is eventually observed, but the accumulation of FADH• indicates that a substantial portion of reduction occurs in sequential one-electron processes rather than a single two-electron event. In rapid-reaction experiments following the reaction of reduced bcd with crotonyl-CoA and oxidized bcd with butyryl-CoA, long-wavelength-absorbing intermediates are observed that are assigned to bcdred:crotonyl-CoA and bcdox:butyryl-CoA charge-transfer complexes, demonstrating their kinetic competence in the course of the reaction. In the presence of crotonyl-CoA there is an accumulation of semiquinone that is unequivocally the anionic FAD•- rather than the neutral FADH• seen in the absence of substrate, indicating that binding of substrate/product results in ionization of the bcd semiquinone. In addition to fully characterizing the rapid-reaction kinetics of both the oxidative and reductive half-reactions, our results demonstrate that one-electron processes play an important role in the reduction of bcd in EtfAB-bcd.

Noncovalent interactions that tune the reactivities of the flavins in bifurcating electron transferring flavoprotein.

Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize noncovalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin (ETflavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin (Bfflavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of ETflavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of ETflavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement, and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the Bfflavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the Bfflavin O4. This would stabilize the anionic Bfflavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.