Introduction {#s1}
============

CRISPR systems are widespread in archaea and bacteria, providing adaptive immunity against invading mobile genetic elements (MGE) ([@bib39]; [@bib29]). Type III CRISPR systems ([Figure 1](#fig1){ref-type="fig"}) are multi-functional effector proteins that have specialised in the detection of foreign RNA ([@bib41]; [@bib45]). The large subunit, Cas10, harbours two enzyme active sites that are activated by target RNA binding: a DNA-cleaving HD nuclease domain ([@bib38]; [@bib11]; [@bib12]; [@bib21]) and a cyclase domain for cyclic oligoadenylate (cOA) synthesis ([@bib22]; [@bib33]; [@bib36]). The third enzymatic activity of type III systems is situated in the Cas7 subunit of the complex, which cleaves bound RNA targets and in turn regulates Cas10 enzymatic activities ([@bib40]; [@bib36]; [@bib20]; [@bib32]). The cyclase domain polymerises ATP into cOA species consisting of between 3 and 6 AMP subunits (denoted cA~3~, cA~4~ etc.), in varying proportions ([@bib22]; [@bib33]; [@bib36]; [@bib15]; [@bib32]). cOA second messengers activate CRISPR ancillary nucleases of the Csx1/Csm6, Can1 (CRISPR associated nuclease 1) and NucC families, which drive the immune response against MGEs ([@bib22]; [@bib33]; [@bib36]; [@bib15]; [@bib25]; [@bib30]). To date, cA~4~ appears to be the most widely used signalling molecule by type III CRISPR systems ([@bib15]). The ribonuclease activity of Csx1/Csm6 is crucial for the clearance of MGEs ([@bib16]; [@bib13]; [@bib15]), particularly when viral genes are transcribed late in infection, at low levels or mutated ([@bib16]; [@bib18]; [@bib35]).

![Cartoon of type III CRISPR cyclic nucleotide signalling and defence in *Sulfolobus solfataricus*.\
The Cas10 subunit of the type III CRISPR complex synthesises cyclic tetra-adenylate (cA~4~) when viral RNA transcripts are detected. Target RNA cleavage shuts-off cA~4~ synthesis. cA~4~ binds to the CARF (CRISPR associated Rossmann Fold) domain of CRISPR ancillary nuclease Csx1 and allosterically activates its HEPN (Higher Eukaryotes and Prokaryotes Nucleotide binding) domain, which degrades RNA non-specifically within the cell. Extant cA~4~ is degraded slowly by CRISPR ring nucleases (Crn1 family) which likely facilitate cell recovery after clearing the virus. Viral anti-CRISPR ring nucleases (AcrIII-1 family) degrade cA~4~ rapidly to stop activation of ancillary defence enzymes such as Csx1 and supress antiviral immunity.](elife-55852-fig1){#fig1}

In our previous study, we demonstrated that the type III-D system from *Sulfolobus solfataricus* synthesises predominantly cA~4,~ which activates the CRISPR ancillary ribonuclease Csx1. We examined the first regulatory step in cOA synthesis in detail and demonstrated that target RNA cleavage and dissociation from the complex shut-off cOA synthesis ([@bib36]). Since CRISPR ancillary nucleases degrade nucleic acids non-specifically, cellular as well as viral targets are destroyed. Collateral cleavage of self-transcripts by a Csm6 enzyme has previously been shown to result in cell growth arrest ([@bib35]). Therefore, in addition to regulating the synthesis of cOA, cells need a mechanism to remove extant cOA if they are to return to normal growth. To solve this problem, *S. solfataricus* encodes CRISPR-associated ring nuclease 1 (Crn1) family enzymes ([@bib2]). Crn1 enzymes slowly degrade cA~4~ to yield di-adenylate products incapable of activating Csx1. In other species Csm6 proteins have evolved catalytic CARF domains capable of degrading cA~4~, thereby acting as their own 'off-switches' to their RNase activity ([@bib3]; [@bib17]). Unsurprisingly, archaeal viruses and bacteriophage have co-opted this regulatory strategy in order to subvert type III CRISPR defence. Many archaeal viruses and bacteriophage encode a ring nuclease anti-CRISPR (AcrIII-1), unrelated to Crn1, which neutralises the type III response by rapidly degrading cA~4~ to prevent ancillary nuclease activation ([@bib4]).

It is clear that the cA~4~ antiviral second messenger is at the centre of a network of interactions that are crucial for effective type III CRISPR defence against MGE. Here, we show that detection of even a single molecule of invading RNA has the potential to generate a large signal amplification via synthesis of cA~4~ that in turn activates the non-specific degradative ribonuclease Csx1. We explore how a cellular ring nuclease can return the cell to a basal state and how viruses can subvert the system. By quantifying and modelling the equilibria and reactions that take place in the arena of type III CRISPR defence, we build a comprehensive model of this dynamic, life or death process.

Results {#s2}
=======

While the control of cOA synthesis by target RNA binding and cleavage is now understood reasonably well, the full implications of cOA generation in a virally-infected cell are not. This requires a detailed knowledge of the levels of cOA produced, consequences for antiviral defence enzymes and the effects of cOA degrading enzymes from cellular and viral sources. These were the aims of our study.

Signal amplification on cA~4~ production {#s2-1}
----------------------------------------

We first investigated the extent of signal amplification that occurs in a cell from detection of a single viral RNA and generation of the cA~4~ second messenger. Using the *S. solfataricus* type III-D CRISPR effector, we varied the concentration of target RNA and quantified the resultant cA~4~ production. As previously observed ([@bib36]), increasing the target RNA concentration resulted in increased cA~4~ production ([Figure 2](#fig2){ref-type="fig"}). Quantification of the concentration of cA~4~ generated was accomplished by using α-^32^P-ATP and quantification of products using a phosphorimager in comparison to standards ([Figure 2---figure supplement 1](#fig2s1){ref-type="fig"}). We observed that approximately 1000 molecules (980 ± 24) of cA~4~ were generated per molecule of RNA, over a range of 10--100 nM target RNA ([Figure 2](#fig2){ref-type="fig"}). When a poorly-cleavable target RNA species containing phosphorothioates was used as the substrate, the amount of cA~4~ generated increased approximately 3-fold (3100 ± 750, [Figure 2](#fig2){ref-type="fig"}), confirming the important role of RNA cleavage for deactivation of the cyclase domain ([@bib36]; [@bib32]). Additionally, analysis of previously published data enabled us to determine the catalytic rate constant (*k*) of cA~4~ synthesis by the *S. solfataricus* type III effector as 0.04 ± 0.01 min^−1^ at 70°C ([Figure 2---figure supplement 2](#fig2s2){ref-type="fig"}).

![Approximately 1000 molecules cA~4~ are made per molecule of RNA target.\
(**A**) Upper panel shows phosphorimages of thin-layer chromatography of cyclic tetra-adenylate (cA~4~) made by *S. solfataricus* (Sso) Csm complex (470 nM carrying the CRISPR RNA A26) across a range of RNA target concentrations (0.1, 1, 10, 25,100 nM) complementary to the A26 CRISPR RNA at 70°C. Lower panel shows cA~4~ synthesised with a cleavage resistant (phosphorothioate) form of the RNA target. (**B**) Bar graph of the concentration of cA~4~ generated with increasing cleavable and cleavage-resistant RNA target generated by quantifying the densiometric signals from A, with an α-^32^P-ATP standard curve ([Figure 2---figure supplement 1](#fig2s1){ref-type="fig"}). Error bars indicate the standard deviation of the mean of three technical replicates, with individual data points shown as clear circles. No data are shown for 1 nM cleavable RNA target as cA~4~ generated was below detection limits. (**C**) Bar chart quantifying the number of molecules of cA~4~ generated per molecule of cleavable or cleavage resistant target RNA across a range of RNA target concentrations. On average SsoCsm synthesised 980 ± 24 and 3100 ± 750 molecules of cA~4~ per molecule of cleavable and cleavage resistant target RNA, respectively. (**D**) Cartoon depicting the cellular implications of \~1000 molecules of cA~4~ generated per molecule of RNA target, which in *S. solfataricus* would equate to \~6 µM cA~4~ within the cell.\
Figure 2---source data 1.Excel spreadsheet with raw data.](elife-55852-fig2){#fig2}

Given that *S. solfataricus* cells are cocci with a diameter of approximately 0.7 µm, the volume of an average cell can be calculated as approximately 0.3 fL (by comparison, *E. coli* has a cell volume of 1 fL \[[@bib23]\]). Using Avogadro's number, 1000 molecules equates to an intracellular concentration of 6 µM cA~4~ in *S. solfataricus*. Thus, detection of one viral RNA in the cell could result in the synthesis of 6 µM cA~4~, 10 RNAs -- 60 µM, etc. The upper limits of cA~4~ generation could be defined by the number of viral target RNAs present, the number of type III effectors carrying a crRNA matching that target, or even conceivably the amount of ATP available for cA~4~ generation.

Kinetic parameters of the Csx1 ribonuclease {#s2-2}
-------------------------------------------

The cA~4~ second messenger binds to CARF family proteins to elicit an immune response. To understand the concentration of cA~4~ required to activate an antiviral response, we determined the dissociation constant of the major ancillary ribonuclease Csx1 for the cA~4~ activator. Using radioactive cA~4~, we titrated an increasing concentration of Csx1 protein and subjected the mixture to native gel electrophoresis ([Figure 3A,B](#fig3){ref-type="fig"}). cA~4~ was bound by Csx1 with a dissociation constant (*K*~D~) of 130 ± 20 nM. Thus, even one viral target RNA detected by the type III CRISPR system should generate enough cA~4~ (6 µM) to fully activate the Csx1 ribonuclease for defence. We proceeded to estimate the binding affinity of a ribonuclease-deficient Csx1 variant for its RNA target, yielding an apparent dissociation constant of approximately 5 µM ([Figure 3C](#fig3){ref-type="fig"}). Finally, we measured the initial velocity of RNA degradation at a variety of RNA concentrations under steady state conditions and determined the multiple-turnover kinetic constant (*k*~cat~) for cA~4~-activated RNA cleavage by Csx1 as 0.44 ± 0.03 min^−1^ at 70°C ([Figure 4](#fig4){ref-type="fig"}).

![Csx1 binds cA~4~ with high affinity and RNA with relatively low affinity.\
(**A**) Phosphorimage of native gel electrophoresis visualising cA~4~ (20 nM) binding by Csx1 (concentrations as indicated in the figure). The other bands visible are due to unreacted ATP and other linear nucleotide products. (**B**) Plot of fraction of cA~4~ bound by Csx1. Error bars indicate the standard deviation of the mean of four technical replicates and the data were fitted to a quadratic equation incorporating a term for non-specific binding. The inset plot is a magnification of cA~4~ binding between 0.01 and 5 µM Csx1 dimer concentrations. (**C**) Phosphorimage of native gel electrophoresis visualising A1 substrate RNA binding by Csx1 H345N protein dimer in the absence (left hand-side) or presence (right hand-side) of unlabelled cA~4~ (20 µM). The image shown is representative of three technical replicates. Control c -- RNA alone.\
Figure 3---source data 1.Csx1 binding to cA~4~ and RNA.](elife-55852-fig3){#fig3}

![Degradation of RNA by Csx1.\
Analysis of multiple-turnover, steady state kinetics of A1 RNA cleavage by Csx1 (125 nM dimer) at 70°C. The data were fitted to the Michaelis-Menten equation and error bars show the standard deviation of the mean of three technical replicates.\
Figure 4---source data 1.Excel spreadsheet with data for kinetics of Csx1.](elife-55852-fig4){#fig4}

Kinetic and equilibrium constants of the ring nucleases Crn1 and AcrIII-1 {#s2-3}
-------------------------------------------------------------------------

We have previously established that Crn1 cleaved cA~4~ at a rate of 0.089 ± 0.003 min^−1^ at 50°C, while AcrIII-1 cleaved cA~4~ at a rate of 5.4 ± 0.38 min^−1^, about 60-fold faster ([@bib4]). The difference in reaction rates probably reflects the different roles of the two enzymes, with Crn1 working in conjunction with the type III CRISPR defence and AcrIII-1 opposing it. To quantify the interaction between ring nucleases and cA~4~, we titrated radioactively labelled cA~4~ with either Crn1 or AcrIII-1 and visualised cA~4~ binding by phosphorimaging following native gel electrophoresis. Crn1 bound cA~4~ with an apparent *K~D~* of \~50 nM, while the inactive H47A variant of AcrIII-1 bound cA~4~ with an apparent *K*~D~ of \~25 nM ([Figure 5](#fig5){ref-type="fig"}). Thus, both ring nucleases bound cA~4~ three to four-fold more tightly than Csx1.

![Crn1 and AcrIII-1 bind cA~4~ with high affinity.\
Phosphorimages of native gel electrophoresis visualising radiolabelled cyclic oligoadenylate (cOA) binding by (**A**) Crn1 (**B**) and catalytically inactive AcrIII-1 (SIRV1 gp29 H47A). Crn1 binds cA~4~ (10 nM) with an apparent dissociation constant (*K*~D~) of \~50 nM, whereas AcrIII-1 binds cA~4~ with an apparent *K*~D~ of \~25 nM. The images shown are representative of three technical replicates. Control c -- cOA alone. The other bands near the bottom of the gel are caused by unreacted ATP and other linear products.](elife-55852-fig5){#fig5}

Kinetic modelling of the antiviral signalling pathway and its regulation by cA~4~ degrading enzymes {#s2-4}
---------------------------------------------------------------------------------------------------

We entered the experimentally determined kinetic and equilibria parameters into the KinTek Global Kinetic Explorer software package and generated a model to simulate RNA degradation by Csx1 and the effects of ring nucleases over time ([Figure 6A](#fig6){ref-type="fig"}). Enzyme concentrations were set initially at 1 µM, based on published studies of transcript levels ([@bib34]; [@bib44]), but were varied during modelling to assess the influence of enzyme concentration on RNA cleavage. We standardised the reaction temperature at 70°C, close to the growth temperature of *S. solfataricus*. This necessitated an estimation of the rate constants of Crn1 and AcrIII-1, which were measured at lower temperatures, based on a 2-fold increase in activity for each 10°C rise in temperature, in line with our previous studies of these enzymes ([@bib2]; [@bib4]). To model the generation of cA~4~ by the Csm:target RNA complex, we set the concentration of that complex at 6, 60 or 600 µM (equivalent to low, medium and high levels of infection), allowing the generation of the equivalent concentration of cA~4~ with the measured rate constant of 0.04 min^−1^ ([Figure 6A](#fig6){ref-type="fig"}). Henceforth, this will be referred to as 6, 60 or 600 µM cA~4~ for simplicity.

![Modelling of *S. solfataricus* antiviral signalling.\
(**A**) Schematic showing kinetic and equilibrium parameters inserted into the KinTek Global Kinetic Explorer software for modelling the type III CRISPR defence illustrated in [Figure 1](#fig1){ref-type="fig"}. Parameters have been determined in this study with the following exceptions: Crn1 rate constant of 0.46 min^−1^ at 70°C estimated from rate of 0.23 min^−1^ measured at 60°C and AcrIII-1 rate constant of 21.6 min^−1^ at 70°C estimated from rate of 5.4 min^−1^ measured at 50°C ([@bib4]). The parameter 'target-Csm-ATP' was set at 6, 60 or 600 µM in simulations. Underscores connecting two variables indicate their relationship in a complex. cA~4~, cyclic tetra-adenylate; Crn1, CRISPR ring nuclease 1; AcrIII-1, viral ring nuclease anti-CRISPR SIRV1 gp29; Csx1, CRISPR ancillary ribonuclease; A~2~ \>P, di-adenylate containing 2',3' cyclic phosphate (product of cA~4~ cleavage). (**B**) Free cA~4~ (600 µM, blue; 60 µM, red; 6 µM, green) concentration and (**C**), RNA cleavage in presence of 1 µM Crn1 and 1 µM Csx1. (**D**) and (**E**) show equivalent plots in the presence of 1 µM AcrIII-1. Insets where present show a magnified view of the start of each plot. (**F**) 3D plot visualising concentration of RNA (1000 µM at start) cleaved by Csx1 in response to 60 µM cA~4~ made by Csm complex, 1 µM Crn1 and varying amounts of AcrIII-1 across a range of doubling endpoints. (**G**) 3D plot visualising concentration of RNA (1000 µM at start) cleaved by Csx1 in response to 60 µM cA~4~ made by Csm complex, and varying concentrations of Csx1, Crn1 and AcrIII-1 enzymes.\
Figure 6---source data 1.Raw data from modelling.\
Figure 6---source data 2.Extracted data from modelling.\
Figure 6---source code 1.KinTek model code for kinetic modelling in [Figure 6](#fig6){ref-type="fig"}.](elife-55852-fig6){#fig6}

In the absence of any ring nuclease activity, 1 µM Csx1 was fully activated, degrading all 1000 µM RNA provided in the simulation within 2400 min, regardless of the simulated level of infection ([Figure 6---figure supplement 1](#fig6s1){ref-type="fig"}). The lack of titration of Csx1 activity observed here is due to its high affinity for cA~4~, resulting in full activation even at very low simulated infection levels. We observed a rapid initial increase in the levels of cA~4~ due to synthesis by the activated Csm complex ([Figure 6B](#fig6){ref-type="fig"}). In the simulations, 1 µM Crn1 degraded 60 µM cA~4~ within 120 min, but took over 1200 min to degrade 600 µM cA~4~ ([Figure 6B](#fig6){ref-type="fig"}). These differences were reflected in the levels of Csx1-mediated RNA degradation, which were much higher for 600 µM cA~4~ than for the lower concentrations ([Figure 6C](#fig6){ref-type="fig"}). Thus, the kinetic modelling demonstrated that addition of a ring nuclease activity allows the cell to respond differently to varied levels of infection, and therefore cA~4~ concentration, resulting in a range of RNA degradation levels.

We next introduced the AcrIII-1 enzyme to the simulation, to model the effect of the Acr on the host defence system. Strikingly, when AcrIII-1 was present at 1 µM, even 600 µM cA~4~ was degraded rapidly (the bulk within 10 min). As a result, Csx1 activity and therefore RNA degradation were strongly supressed ([Figure 6D,E](#fig6){ref-type="fig"}). The long tail of cA~4~ observed in these simulations represents cA~4~ bound to the Csx1 enzyme and thus unavailable for degradation by AcrIII-1, reflected in the persistence of active Csx1-cA~4~ complex in the simulations over an extended period ([Figure 6---figure supplement 2](#fig6s2){ref-type="fig"}). Notably, at the highest simulated infection level of 600 µM cA~4~, Csx1 persists in a 100% activated form for over 1200 min in the absence of AcrIII-1, but is reduced to 50% activated within 150 min in its presence, with similar effects seen at lower cA~4~ concentrations ([Figure 6---figure supplement 2](#fig6s2){ref-type="fig"}).

The concentration of AcrIII-1 within *S. solfataricus* cells during infection is not known. We therefore varied AcrIII-1 concentration in the model (with 60 µM cA~4~, 1 µM Csx1 and Crn1) and simulated RNA cleavage ([Figure 6F](#fig6){ref-type="fig"}). AcrIII-1 concentrations as low as 100 nM reduced RNA cleavage by about 30%, and higher levels reduced RNA cleavage by Csx1 markedly ([Figure 6F](#fig6){ref-type="fig"}). This relationship held in simulations with 600 µM cA~4~ ([Figure 6F](#fig6){ref-type="fig"}, [Figure 6---figure supplement 3](#fig6s3){ref-type="fig"}). During infection, a positive correlation between AcrIII-1 concentration and viral transcript levels would be required for continued escape from Csx1-mediated CRISPR defence -- a reasonable assumption. Finally, we simulated the effects of increasing Csx1 on RNA cleavage in our infection model. Increasing Csx1 concentration 10-fold led to a significantly more robust immune response (reflected in high levels of RNA degradation) ([Figure 6G](#fig6){ref-type="fig"}) that could not be fully overcome by increasing the concentrations of Crn1 or AcrIII-1 in proportion (compare for example 10,10,10 with 1,1,1 µM in [Figure 6G](#fig6){ref-type="fig"}). Hence, increasing the expression levels of defence nucleases such as Csx1 in response to viral infection could provide an effective means for cells to combat viruses armed with Acrs such as AcrIII-1, but perhaps with the cost of increased toxicity.

Discussion {#s3}
==========

Signal amplification in type III CRISPR defence {#s3-1}
-----------------------------------------------

In this study, we used biochemical data to build a kinetic model of the type III CRISPR antiviral signalling pathway within *S. solfataricus* cells and examined the capacity of CRISPR and anti-CRISPR ring nucleases for its regulation. Quantification of cA~4~ generated by the SsoCsm complex in vitro revealed that \~ 1000 molecules of cA~4~ are made per RNA target, amounting to a concentration of 6 µM in the cell if replicated in vivo. This large degree of signal amplification would ensure that detection of 1 RNA target could generate sufficient amounts of cA~4~ to fully activate the ribonuclease effector protein Csx1, which has a dissociation constant for cA~4~ of 0.13 µM. Given the large signal amplification observed here, it seems likely that some means of cOA degradation, either via self-limiting ribonucleases ([@bib3]; [@bib17]) or dedicated ring nucleases ([@bib2]), will be essential for type III CRISPR systems to provide immunity rather than elicit abortive infection. Indeed, growth arrest has been observed for cOA activated Csm6 during bacteriophage infection ([@bib35]). This life or death decision in response to genotoxic stress has also been observed in *S. islandicus*, which becomes dormant upon viral infection and subsequently dies if virus remains in culture ([@bib5]). In recent years, diverse CRISPR systems have been implicated in abortive infection or cell dormancy. The Type I-F CRISPR system of *Pectobacterium atrosepticum* was found to provide population protection by aborting infection when infected by virulent phage ([@bib42]). Likewise, the in-trans collateral RNA cleavage of *Listeria seeligeri* Cas13a resulted in cell dormancy, providing herd immunity to the bacterial population ([@bib31]). Similarly, in ecological contexts, it is possible that different multiplicities of viral infection illicit different outcomes from the type III CRISPR response that benefit either the individual cell or the population.

Cellular and viral ring nucleases reset the system in fundamentally different ways {#s3-2}
----------------------------------------------------------------------------------

Biochemical comparison of Crn1 and AcrIII-1 revealed that both enzymes bind cA~4~ with dissociation constants around 40 nM, around 3-fold tighter than observed for Csx1. However, Crn1 is a much slower enzyme. Kinetic modelling of the antiviral signalling pathway confirms that Crn1 is effective only at low levels of viral gene expression, where it has the potential to neutralise the toxicity associated with cA~4-~activated ribonucleases to offer a route for cell recovery without abrogating immunity. In contrast, the much faster reaction kinetics of the anti-CRISPR ring nuclease means it can rapidly deactivate Csx1 and immunosuppress cells even under very high RNA target (and thus cA~4~) levels.

Our modelling suggests that the rapid turnover of cA~4~ by AcrIII-1 over a wide concentration range greatly limits RNA cleavage by deactivating defence enzymes. Therefore, the deployment of AcrIII-1 upon viral infection may not only promote viral propagation but also safeguard cellular integrity until viral release by lysis. Recent studies have uncovered that sequentially infecting phage evade CRISPR defences by exploiting the immunosuppression achieved by Acr enzymes from failed infections ([@bib6]; [@bib24]). Further, these immunosuppressed cells have been shown to be susceptible to Acr-negative phage infections, highlighting the complex ecological consequences of supressing CRISPR immunity ([@bib8]). In *Sulfolobus* Turreted Icosahedral virus (STIV), the AcrIII-1 gene *B116* is expressed early in the viral life cycle ([@bib34]). Therefore, AcrIII-1 accumulation in the cell, possibly from early expression by unsuccessful viruses may, as our models demonstrate, favour the success of latter viral infections. Type III CRISPR systems also conditionally tolerate prophage ([@bib14]), and unsurprisingly, AcrIII-1 is found in a number of prophages and integrative and conjugative elements. In these cases, constitutively expressed AcrIII-1 may further immunocompromise cells, and sensitise them to infection by phage otherwise eradicated by type III CRISPR defence. In the ongoing virus-host conflict, while increasing Csx1 concentration may allow better immunity when faced with AcrIII-1, upregulating AcrIII-1 expression in cells will undoubtedly offer viruses an avenue for counter offence.

It should be noted that the type III CRISPR locus of *S. solfataricus* contains a number of CARF domain proteins and their contribution to immunity has not yet been studied. In particular, the CARF-family putative transcription factor Csa3 appears to be involved in transcriptional regulation of CRISPR loci, including the adaptation and type I-A effector genes, when activated by cA~4~ ([@bib27]; [@bib28]). These observations suggest that the cA~4~ signal may transcend type III CRISPR defence in some cell types by activating multiple defence systems. However, by degrading the second messenger, AcrIII-1 has the potential to neutralise all of these.

Cyclic nucleotides in prokaryotic defence systems {#s3-3}
-------------------------------------------------

Cyclic nucleotide-based defence systems are emerging as powerful cellular sentinels against parasitic elements in prokaryotes. Mirroring the role of cyclic GMP-AMP synthase (cGAS) in eukaryotic defence against viruses as part of the cGAS-STING pathway, bacterial cGAS enzymes have recently been discovered that abort infection by activating phospholipases through cGAMP signaling ([@bib9]). Termed the cyclic-oligonucleotide-based antiphage signaling system (CBASS), a large number of additional cOA sensing effector proteins associated with CBASS loci remain uncharacterised, highlighting great diversity in the cellular arsenal used for defence ([@bib7]; [@bib9]). Furthermore, diverse nucleotide cyclases have been identified that generate a range of cyclic nucleotides including cUMP-AMP, c-di-UMP and cAAG, which are also likely to function in novel antiviral signal transduction pathways ([@bib43]). Type III systems also generate cyclic tri-adenylate (cA~3~) and cyclic penta-adenylate (cA~5~) molecules. Whereas no signalling role has yet been ascribed to cA~5~, cA~3~ has been demonstrated to activate a family of DNases termed NucC which abort infection by degrading the host genome prior to completion of the phage replication cycle ([@bib25]).

The balance between immunity, abortive infection and successful pathogen replication is likely to be governed by enzymes that synthesise and degrade these cyclic nucleotide second messengers. Just as prokaryotes with type III CRISPR require a means to degrade cOA in appropriate circumstances, eukaryotic cells have enzymes that degrade cGAMP to regulate cGAS-STING mediated immunity ([@bib26]). Likewise, while prokaryotic viruses utilise AcrIII-1 to rapidly degrade cA~4~, eukaryotic poxviruses utilise Poxins to subvert host immunity by destroying cGAMP ([@bib10]), and pathogenic Group B *Streptococci* degrade host c-di-AMP using the CndP enzyme to circumvent innate immunity ([@bib1]). The rate of discovery of new defence pathways and cyclic nucleotide signals is breath-taking. Analysis of the dynamic interplay between enzymes that leads to fluctuations in the levels of these second messengers is therefore of crucial importance if we are to achieve an understanding of these processes.

Materials and methods {#s4}
=====================

  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  Reagent type\                                       Designation                    Source or reference                                             Identifiers          Additional\
  (species) or resource                                                                                                                                                   information
  --------------------------------------------------- ------------------------------ --------------------------------------------------------------- -------------------- ------------------------------
  Gene (*Sulfolobus solfataricus*)                    Csm complex (eight subunits)   PMID:[24119402](https://www.ncbi.nlm.nih.gov/pubmed/24119402)                        virus expression construct

  Gene\                                               Csx1                           PMID:[29963983](https://www.ncbi.nlm.nih.gov/pubmed/29963983)   UniProtKB - Q97YD5   plasmid expression construct
  (*Sulfolobus solfataricus*)                                                                                                                                             

  Gene (*Sulfolobus solfataricus*)                    Crn1                           PMID:[30232454](https://www.ncbi.nlm.nih.gov/pubmed/30232454)   UniProtKB - Q7LYJ6   plasmid expression construct

  Gene (*Sulfolobus islandicus rod-shaped virus 1*)   AcrIII-1                       PMID:[31942067](https://www.ncbi.nlm.nih.gov/pubmed/31942067)   UniProtKB - Q8QL27   plasmid expression construct

  Software, algorithm                                 KinTek\                        PMID:[19897109](https://www.ncbi.nlm.nih.gov/pubmed/19897109)                        model constructed\
                                                      Kinetic Explorer                                                                                                    for this paper
  ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Cyclic oligoadenylate (cOA) synthesis and visualisation {#s4-1}
-------------------------------------------------------

Cyclic tetra-adenylate (cA~4~) made per RNA target (0.01, 0.1, 1, 10, 25 or 50 nM) was investigated in a 20 μl reaction volume incubating A26 RNA target or A26 phosphorothioate RNA target ([Table 1](#table1){ref-type="table"}) with 13.5 μg *Sulfolobus solfataricus* (Sso)Csm complex (\~470 nM carrying A26 CRISPR RNA) in Csx1 buffer containing 20 mM MES pH 5.5, 100 mM K-glutamate, 1 mM DTT and 3 units SUPERase•In Inhibitor supplemented with 1 mM ATP, 5 nM α-^32^P-ATP and 2 mM MgCl~2~ at 70°C for 2 hr. All samples were deproteinised by phenol-chloroform extraction (Ambion) followed by chloroform (Sigma-Aldrich) extraction prior to separating the cOA products by thin-layer chromatography (TLC). TLC was carried out as previously described ([@bib37]). In brief, 1 μl of radiolabelled cOA product was spotted 1 cm from the bottom of a 20 × 20 cm silica gel TLC plate (Supelco Sigma-Aldrich). The TLC plate was placed in a sealed glass chamber pre-warmed at 37°C containing 0.5 cm of a running buffer composed of 30% H~2~O, 70% ethanol and 0.2 M ammonium bicarbonate, pH 9.2. After TLC the plate was air dried and sample migration visualised by phosphor imaging. For analysis, densiometric signals corresponding to cA~4~ was quantified as previously described ([@bib37]).

###### Oligonucleotides.

CRISPR RNA A26 is shown 3' to 5'. Phosphorothioate linkages are indicated with an asterisk. Regions complementary to CRISPR RNA A26 are italicized.

  --------------------------------- --------------------------------------------------------------
  Crispr rna a26                    3'-GCAACAATTCTTGCTGCAACAATCTTCAACCCATACCAGAAAGUUA
  Name                              Sequence (5'−3')
  Target RNA A26                    AGGGUCGUUGUUAAGAACGACGUUGUUAGAAGUUGGGUAUGGUGGAGA
  Phosphorothioate target RNA A26   AGGGUCGUUGUUAAGAACGACGUUGU\*U\*A\*GAAGUUGGGU\*A\*U\*GGUGGAGA
  A1 substrate RNA                  AGGGUAUUAUUUGUUUGUUUCUUCUAAACUAUAAGCUAGUUCUGGAGA
  --------------------------------- --------------------------------------------------------------

Generation of α-^32^P-ATP standard curves {#s4-2}
-----------------------------------------

cA~4~ synthesis was visualised by incorporation of 5 nM α-^32^P-ATP added together with 0.5 mM ATP at the start of the reaction. Therefore, to calculate the concentration of ATP used for cA~4~ synthesis, α-^32^P-ATP standard curves were generated in duplicate, starting with 5 nM α-^32^P-ATP within a 20 μl volume to represent the densiometric signal corresponding to the complete conversion of 0.5 mM ATP into cOA. Serial two-fold dilutions of 5 nM α-^32^P-ATP and 0.5 mM ATP starting from a 20 μl volume were made and 1 μl of each dilution was spotted on a silica plate and phosphorimaged alongside TLC separating cOA made with varying RNA target concentrations. After phosphorimaging, the densiometric signals of the serial dilutions were quantified, averaged and plotted against ATP concentration starting from 0.5 mM and halving with each two-fold dilution. A line of best fit was then drawn. The concentration of ATP used to synthesise cA~4~ was calculated by entering the densiometric signal of the cA~4~ product into to equation of the line of best fit for the α-^32^P-ATP standard curve. The concentration of cA~4~ generated was derived by dividing the concentration of ATP incorporated by four to account for polymerisation of four ATP molecules to generate one molecule of cA~4~. Finally, the molecules of cA~4~ made per RNA was calculated by dividing the cA~4~ concentration generated by the concentration A26 RNA target used for cOA synthesis.

Calculation determining the concentration of cA~4~ made when one RNA target is detected within a *S. solfataricus* cell of ≈ 0.8 µm (0.6--1.0 µm) diameter.$$Volume\ \left( V \right) = \frac{4}{3}\pi r^{3}\,\,{and}\,\, r = \frac{1}{2}d$$$$r = \frac{1}{2}\  \times \ 0.8\,{\mu m}$$$$r = 0.4\,\mu m$$$$V = \frac{4}{3}\pi \times (0.4\,\mu m)^{3}$$$$V = 0.268{\,\mu m}^{3}$$$$1\,{\mu m}^{3} = 1{fL}$$$$0.268{\,\mu m}^{3} = 0.268{\, fL} = 2.68 \times 10^{- 13}{\, mL}$$$$1\,\, mole\,\, of\,\, RNA\,\, = \,\, 6.022 \times 10^{23}\,\, molecules\,\, of\,\, RNA$$$$1\,\, molecule\,\, of\,\, RNA\,\, = 1 \div \,\, 6.022 \times 10^{23}\,\, = 1.661 \times 10^{- 24}\,\, moles\,\, of\,\, RNA$$

As \~1000 molecules of cA~4~ is made per 1 molecule of RNA$$1.661\  \times 10^{- 24\ }moles\  \times \ 1000\  = 1.661\  \times 10^{- 21}\ mole{s\,\, of\,\, cA}_{4}$$

Concentration (M) = moles / Volume (L)$$1.661\  \times 10^{- 21}\ moles\  \div \ 2.68 \times 10^{- 16}\, L = 6.20 \times 10^{- 6}\,{M\, or\, 6.20\,\mu M\,\, cA_{4}}$$

Electrophoretic mobility shift assays to determine cA~4~ equilibrium binding constants {#s4-3}
--------------------------------------------------------------------------------------

\~20 nM radioactively-labelled cA~4~ generated using the SsoCsm was incubated with increasing concentrations of Csx1 (0.01, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80, 1,0, 2.0, 4.0, 8.0, 10.0, 20.0 μM protein dimer) in buffer containing 20 mM Tris-HCl pH 7.5, 150 mM NaCl, 2 mM MgCl~2~ supplemented with 2 µM Ultrapure Bovine Serum Albumin (Invitrogen) for 10 min at 25°C. A reaction volume equivalent of 20% (v/v) glycerol was then added prior to loading the samples on a 15% polyacrylamide, 1 X TBE gel. Electrophoresis was carried out at 28°C and 250 V. Gels were phosphor imaged overnight at −80°C. For investigating RNA binding, 50 nM 5'-end radiolabelled and gel purified A1 RNA was incubated with Csx1 variant H345N (0.01, 0.10, 1.0, 5.0, 10.0, 20.0 μM protein dimer) in the presence or absence of 20 µM cA~4~ for 15 min at 40°C. To examine cA~4~ binding by Crn1,\~10 nM radiolabelled SsoCsm cA~4~ was incubated with Sso2081 (0.01, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80, 1,0, 2.0, 4.0, 8.0, 10.0, 20.0 μM protein dimer) on ice for 15 min before gel electrophoresis as described above but at 300V and at 4°C. cA~4~ binding by AcrIII-1 was examined by incubating \~10 nM radiolabelled SsoCsm cA~4~ with SIRV1 gp49 H47A (0.001, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08, 0.10, 1.0, 10.0, 20.0 µM protein dimer) for 10 min at 25°C before gel electrophoresis at 30°C as described above. For analysis densiometric signal corresponding to cA~4~ bound protein was quantified. The densiometric count corresponding to cA~4~ bound to 20 µM Csx1 dimer was used to represent 100% binding and densiometric counts from other lanes were normalised to this value within each replicate. Error of the 100% bound (20 µM Csx1 dimer) densiometric count was derived by calculating the area adjusted count for each replicate and then the standard deviation of their mean, reporting the standard deviation as a fraction of the mean set as 100% bound. The data were fitted to a quadratic equation with an adjust for nonspecific binding (Y = (ymax-ymin)\*((x + \[ligand\] + Kd) - ((x+\[ligand\]+Kd)\^2--4\*x\*\[ligand\])\^0.5)/(2\*\[ligand\]) +a\*x+b) using GraphPad Prism 8.

Multiple turnover kinetics of RNA cleavage by Csx1 {#s4-4}
--------------------------------------------------

Multiple turnover kinetic experiments were carried out by incubating Csx1 (0.125 µM dimer) with radiolabelled (50 nM) and unlabelled RNA A1 to a final concentration of 1.25, 1.5, 2.0, 5.0, 10.0, or 20.0 μM RNA in Csx1 buffer at 70°C. Control reactions with no protein and with protein and RNA in the absence of cA~4~ were included. 10 μl reaction aliquots were quenched by adding to phenol-chloroform and vortexing and the different time-points at which reactions were quenched for each RNA concentration is indicated in data transparency. Deproteinised products were run on a 7 M urea, 20% acrylamide, 1 X TBE gel at 45°C as previously described ([@bib37]), and phosphorimaged overnight at −80°C. Each experiment was carried out in triplicate. Fraction of RNA cut was plotted against time and a linear regression was carried out on the linear portion of each plot to determine the initial rate of RNA cleavage at each RNA concentration. To determine the parameters *k*~cat~ and *K*~M~, the average initial rate of RNA cleavage was plotted against RNA concentration and fitted to the Michaelis-Menten equation using GraphPad Prism 8.

Modelling antiviral signalling and its control by ring nucleases {#s4-5}
----------------------------------------------------------------

Modelling was carried out using the KinTek Explorer eight software package ([@bib19]), which is available from (<https://kintekcorp.com/software>). Experiments were modelled and simulated using kinetic and equilibrium paramters detemined experimentally as described in [Figure 6A](#fig6){ref-type="fig"}. The following steps were inserted to generate the model:

1.  target_Csm_ATP = target_Csm + cA4 (irreversible)

2.  cA4 + Csx1 = cA4_Csx1

3.  cA4_Csx1 + RNA = cA4_Csx1_RNA

4.  cA4_Csx1_RNA = cA4_Csx1_cleavedRNA (irreversible)

5.  cA4_Csx1_cleavedRNA = cA4_Csx1 + cleavedRNA

6.  cA4 + Crn1 = cA4_Crn1

7.  cA4_Crn1 = A2 + Crn1 (irreversible)

8.  cA4 + Vrn = cA4_Vrn cA4_Vrn = A2 + Vrn (irreversible)

Simulations were carried out varying target_Csm_ATP concentration (6, 60 and 600 µM) while Csx1, Crn1 (Sso2081) and AcrIII-1 concentrations were fixed at 1 µM dimer, or varied depending on the simulation, with total substrate RNA in the cell fixed at 1000 µM.

Funding Information
===================

This paper was supported by the following grants:

-   http://dx.doi.org/10.13039/501100000268Biotechnology and Biological Sciences Research Council BB/S000313/1 to Malcolm F White.

-   http://dx.doi.org/10.13039/100004440Wellcome 210486/Z/18/Z to Clarissa M Czekster.

This work was supported by a grant from the Biotechnology and Biological Sciences Research Council (Grant REF BB/S000313/1 to MFW) and the Wellcome Trust (Grant 210486/Z/18/Z to CMC)

Additional information {#s5}
======================

No competing interests declared.

Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing.

Methodology, Writing - review and editing.

Methodology, Writing - review and editing.

Methodology, Writing - review and editing.

Visualization, Methodology, Writing - review and editing.

Conceptualization, Formal analysis, Supervision, Funding acquisition, Writing - original draft, Project administration, Writing - review and editing.

Additional files {#s6}
================

Data availability {#s7}
=================

All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 2, 3, 4 and 6.

10.7554/eLife.55852.sa1

Decision letter

Wiedenheft

Blake

Reviewing Editor

Montana State University

United States

Lawrence

Martin

Reviewer

Sorek

Rotem

Reviewer

Weizmann Institute of Science

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Thank you for submitting your article \"The dynamic interplay of host and viral enzymes in type III CRISPR-mediated cyclic nucleotide signalling\" for consideration by *eLife*. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Gisela Storz as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Philip Kranzusch (Reviewer \#1); Martin Lawrence (Reviewer \#2); Rotem Sorek (Reviewer \#3).

The reviewers universally agree that this work is nearly ready for publication. No additional experiments are requested.

Summary:

It is becoming increasingly clear that cyclic nucleotides play a critical role in the host response to a viral infection and that viruses have evolved mechanisms to intercept these signaling molecules. In this work, the authors focused on the synthesis of cyclic-oligo-adenylate (cA~4~) produced by the type-III CRISPR system in *Sulfolobus solfataricus*, the immune proteins that are activated by cA~4~, and viral anti-CRISPR proteins that degrade cA~4~. Biochemical data is used to build a kinetic model, which helps explains how cyclases and ring nucleases govern cyclic nucleotide levels and determine infection outcomes in virus-host conflict.

Essential revisions:

1\) Rates of cA~4~ synthesis by the *S. solfataricus* Csm complex have been previously measured (Rouillon et al., 2018; Rouillon et al., 2019). These rates should be incorporated into the mathematical model to understand how the relative rates of cA~4~ synthesis and cA~4~ degradation may influence effector function dynamics.

2\) The current model does not take into account previously observations from the lab demonstrating that cA~4~ synthesis wanes after \~10 min if no new activating RNA is delivered to the system. Please discuss this in the context of your model.

10.7554/eLife.55852.sa2

Author response

> Essential revisions:
>
> 1\) Rates of cA~4~ synthesis by the *S. solfataricus* Csm complex have been previously measured (Rouillon et al., 2018; Rouillon et al., 2019). These rates should be incorporated into the mathematical model to understand how the relative rates of cA~4~ synthesis and cA~4~ degradation may influence effector function dynamics.

We agree that this would strengthen the model. We have determined the rate of cA~4~ synthesis (Figure 2---figure supplement 2) from previously published data (Rouillon et al., 2018) and have incorporated this into the model. We have also adjusted the kinetic parameters of our model to the reaction temperature of 70 °C, which better represents infection kinetics at the growth temperature of *S. solfataricus*. We have also increased the cellular RNA concentration to 1000 µM, from 100 µM in the previous model, which we believe better represents cellular conditions and means that RNA is never fully depleted in the simulations. For comparison and transparency, we include a supplementary version of Figure 6 with simulations starting at fixed cA~4~ concentrations (as in the original submission) carried with the current model as per our initial submission. See Author response images 1 and 2. Ultimately, the changes we have made, including incorporation of the rate of cA~4~ synthesis do not radically alter any of our previous conclusions.

![](elife-55852-resp-fig1){#respfig1}

![](elife-55852-resp-fig2){#respfig2}

> 2\) The current model does not take into account previously observations from the lab demonstrating that cA~4~ synthesis wanes after \~10 min if no new activating RNA is delivered to the system. Please discuss this in the context of your model.

This query relates to point 1, and we have addressed this point in the new model, which now incorporates the time dependent synthesis of cA~4~, dependent on the concentration of target RNA.

[^1]: Department of Molecular Sensory Systems, Center of Advanced European Studies and Research (Stiftung Caesar), Bonn, Germany.
