Mechanism of allosteric inhibition in the Plasmodium falciparum cGMP-dependent protein kinase

Most malaria deaths are caused by the protozoan parasite Plasmodium falciparum. Its life cycle is regulated by a cGMP-dependent protein kinase (PfPKG), whose inhibition is a promising antimalaria strategy. Allosteric kinase inhibitors, such as cGMP analogs, offer enhanced selectivity relative to competitive kinase inhibitors. However, the mechanisms underlying allosteric PfPKG inhibition are incompletely understood. Here, we show that 8-NBD-cGMP is an effective PfPKG antagonist. Using comparative NMR analyses of a key regulatory domain, PfD, in its apo, cGMP-bound, and cGMP analog–bound states, we elucidated its inhibition mechanism of action. Using NMR chemical shift analyses, molecular dynamics simulations, and site-directed mutagenesis, we show that 8-NBD-cGMP inhibits PfPKG not simply by reverting a two-state active versus inactive equilibrium, but by sampling also a distinct inactive “mixed” intermediate. Surface plasmon resonance indicates that the ability to stabilize a mixed intermediate provides a means to effectively inhibit PfPKG, without losing affinity for the cGMP analog. Our proposed model may facilitate the rational design of PfPKG-selective inhibitors for improved management of malaria.

Malaria is a serious health risk for about 40% of the world's population, with more than 200 million diagnoses and 400,000 deaths globally estimated in 2018 (1). Most malarial deaths are caused by Plasmodium falciparum with life cycles in both the human and the mosquito hosts (2,3). The regulation of such cycles relies on cyclic nucleotide signaling pathways (4). Specifically, the P. falciparum guanosine-3Ј,5Ј-cyclic monophos-phate (cGMP) 2 -dependent kinase (PfPKG) plays a major role in both replication and transmission of the Plasmodium parasite, including exflagellation and gametogenesis in the asexual blood stage, invasion of hepatocytes by sporozoites, and gliding motility of ookinetes (5)(6)(7)(8)(9). Hence, PfPKG has emerged as a promising target for development of antimalaria therapies complementary to current malaria drugs, whose effectiveness is limited by increasing resistance (2,10,11). Furthermore, PfPKG inhibition is an effective and novel strategy not only for malaria treatment, but also for malaria prophylaxis (2,6,(12)(13)(14)(15).
One strategy to inhibit PfPKG is through targeting the active site in the kinase domain. Inhibitors that target such kinase active sites are effective, but poorly selective, because the active sites are highly conserved among eukaryotic kinases (16). One approach to circumvent this is to target allosteric sites, such as the cGMP-binding sites in PfPKG.
PfPKG consists of an N-terminal regulatory region followed by a C-terminal catalytic domain (Fig. 1A), similar to mammalian PKG (17)(18)(19). In the absence of cGMP, the regulatory region of PfPKG autoinhibits the catalytic domain. Upon cGMP binding to the regulatory region, the PfPKG autoinhibition is allosterically removed, resulting in activation. The regulatory domain is composed of four cGMP-binding domains (CBDs): A, B, C, and D (Fig. 1A). CBD-C is degenerate and does not bind to cGMP (20). CBD-D (here referred to as "PfD") is adjacent to the catalytic domain (Fig. 1A) and exhibits the highest affinity and selectivity toward cGMP, with a K d of 40 nM, compared with K d Ն 1 M for CBD-A and CBD-B (14). In addition, the construct spanning only PfD and the catalytic domain displays low basal activity in the apo state and is efficiently activated by cGMP (21). Hence, PfD serves as the critical controlling unit for both PfPKG autoinhibition and activation.
PfD is composed of a contiguous ␤-subdomain and noncontiguous ␣-subdomain (Fig. 1B) (14,22). The ␤-subdomain forms a ␤-barrel critical for cGMP recognition and includes the base-binding region (BBR) and the phosphate-binding cassette (PBC) (Fig. 1, B and C), whereas the ␣-subdomain mediates inhibitory interactions with the catalytic domain and includes three N-terminal and two C-terminal helices, called the N3A motif and ␣B-␣C, respectively (Fig. 1, B and D). In the apo form, the N3A contacts the ␤-barrel, whereas the C-terminal helices point away from it ( Fig. 1B) (14). Upon cGMP binding, the N3A moves away from the ␤-subdomain, whereas the C-terminal helices are recruited to it by an arginine residue in the PBC, Arg 484 (Fig. 1B). Arg 484 bridges the PBC to two residues in the ␣C and defines a unique cGMP-capping triad (Fig. 1B), markedly distinct from human PKG, in which cGMP is capped by a Tyr in the C-terminal helix (23). Disruption of the PfPKG-specific capping triad by mutations blocks proliferation in P. falciparum (14), confirming that PfD is an excellent target for treating malaria. Additionally, Arg 528 in the ␣C plays a pivotal role in regulating PfPKG activity by controlling the ␣C conformation through a hydrogen bond with Tyr 480 , which stabilizes the active conformation ( Fig. 1D) (14,21,24).
The allosteric transition of PfD provides a simple model to rationalize the cGMP-dependent activation of the fulllength protein through a four-state allosteric thermodynamic cycle, which arises from the coupling of the PfD inactive-active autoinhibitory equilibrium and the apo-holo cGMP-binding equilibrium (Fig. 1E) (25)(26)(27). Of the four states (Fig. 1E), the holo-inactive is particularly relevant for inhibition. This state is selectively stabilized by partial agonists and/or antagonists, which decouple binding from the intramolecular allosteric pathways necessary for activation. We hypothesize that the stabilization of the holo-inactive state arise from a reversal of the two-state inactive/active conformational equilibrium (Fig. 1F) or from a multistate equilibrium sampling a distinct intermediate conformation that exhibits components of both the inactive and the active states of PfD in different regions (Fig. 1G).
To test these hypotheses and understand the mechanism of action of allosteric partial agonists and antagonists that serve as leads for new PfPKG-selective inhibitors, we use comparative NMR to analyze PfD in its apo, cGMP-bound, and cGMP analog-bound forms. The cGMP analogs investigated here include 8-NBD-cGMP, 8-pCPT-cGMP, and PET-cGMP ( Fig. 1, I-K). Our results reveal that 8-NBD-cGMP and 8-pCPT-cGMP significantly reduce the kinase activity and induce major perturbations throughout PfD by stabilizing an additional intermediate conformer that is distinct from both the native apo-inactive and holo-active conformations. This "mixed" intermediate state exhibits inhibitory potential comparable with that of the apo-inactive state, as the C-terminal lid is disengaged, but it preserves an affinity similar to that of cGMP, as the key cGMP-binding regions, including the PBC and the adjacent ␣B, are still engaged. Hence, the stabilization of an allosterically mixed intermediate emerges as a new strategy for selectively inhibiting PfPKG with high potency and efficacy.

The apo versus cGMP-bound comparative NMR analyses of PfD reveal pervasive allosteric perturbations
As a first step to probe the allosteric effects of cGMP on the PfD construct (residues 401-542) used in our NMR studies, we acquired NH HSQC spectra of apo and cGMP-bound PfD ( Fig.  2A). The cGMP-bound spectrum was assigned through tripleresonance NMR experiments and selective amino acid labeling. However, the apo construct was prone to precipitation during acquisition of 3D NMR data sets. Hence, the apo assignments were obtained by transferring the holo assignments to the apo spectrum through 2D N z -exchange (Fig. 2, B and D) and 2D difference N z -exchange experiments acquired in the presence of substoichiometric amounts of cGMP (Fig. 2, C and E). Based on these assignments, the cGMP-induced residue-specific compounded 15 N-1 H chemical shift changes (⌬CCS) were measured (Fig. 2F).
The majority of the residues exhibit significant ⌬CCS values upon cGMP binding (Fig. 2F). As expected, significant ppm shifts are observed in the PBC and BBR regions, which directly interact with cGMP, as well as the ␤2-3 loop, which is adjacent to the PBC (Fig. 2F). Furthermore, major ppm changes are also detected for several ␣-helical regions, especially the ␣A, ␣B, and ␣C helices (Fig. 2F). The ppm shifts in ␣B-␣C are consistent with the cGMP-dependent structural changes for these regions, as shown by the local root mean square deviation between the apo and cGMP-bound structures (Fig. 2G). The ␣B-␣C structural changes also affect adjacent regions, such as the ␣A helix and the ␤8 strand ( Fig.  2F), although such regions are not subject to major structural changes upon cGMP binding (Fig. 2G). Hence, the cGMP-dependent ⌬CCS values (Fig. 2F) are in full agreement with the structures of apo-and holo-PfD, providing an initial validation of our NMR assignments. Our NMR assignments were further validated by the agreement between the measured and the computed C␣ and C␤ chemical shift values based on the apo-and holo-PfD structures (Fig. 2, H and I). Once the NMR assignments were validated, we extended our analyses to evaluate the effects of cGMP analogs.

cGMP analog versus cGMP chemical shift differences report on inhibitory conformational changes of PfD
We focused on three cGMP base analogs with distinct base modifications: 8-NBD-, 8-pCPT-, and PET-cGMP ( Fig. 1, I-K), as they act as PfPKG activators to different degrees. Whereas PET-cGMP acts as an agonist, as shown by kinase assays (Fig. 1L), 8-pCPT-and 8-NBD-cGMP elicit progressively increasing potential as antagonists. 8-NBD-cGMP is the most potent inhibitor of PfPKG (Fig. S1), followed by 8-pCPT-cGMP. This range of kinase responses (Fig. 1L) makes this group of cGMP analogs an excellent toolset to investigate the mechanisms underlying the allosteric inhibition of PfPKG. For this purpose, we first measured the cGMP analog versus cGMP differences in compounded chemical shifts (Fig. 3A). The overall distributions of such ppm differences (Fig. 3B) rank in full agreement with the inhibitory potencies of the cGMP analogs (Fig. 1L). Furthermore, the EDITORS' PICK: Mechanism of PfPKG allosteric inhibition pervasiveness of the effects caused by the cGMP modification (Fig. 3, A and C-E) also reflects the extent of inhibition. The perturbations resulting from the replacement of cGMP with either 8-NBD-cGMP or 8-pCPT-cGMP extend throughout the PfD structure, whereas the least inhibitory analog, PET-cGMP, perturbs only a relatively limited region near the cGMP-binding pocket (Fig. 3, A and C-E).
Overall, the data of Fig. 3 suggest that chemical shifts are ideally suited to investigate the dynamical allosteric transitions underlying PfPKG inhibition. This notion is fully consistent

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
with the fast exchange between the inactive and the active states of CBDs (27)(28)(29)(30)(31)(32)(33), whereby the observed NMR peak positions are population-weighted averages of the pure inactive and active ppm values. Hence, chemical shifts report on the position of the dynamic inactive-active CBD conformational equilibrium. Chemical shift projection analysis (CHESPA) is an effective means to evaluate the perturbation of such dynamic equilibria caused by ligand modifications (28,29,31,32). In the CHESPA analysis, the NMR peak positions of the cGMP analog-bound sample are evaluated relative to a reference vector connecting the peaks of the apo-inactive and the cGMPbound active samples (Fig. 4, A and D). Specifically, the angle () arising from the cGMP analog-bound versus cGMP-bound vector relative to the reference vector and the normalized projections onto the reference vector (i.e. fractional shift, X) are quantified through CHESPA (Fig. 4, A and D). The cos() values report on the direction of the cGMP analog-induced shifts, whereas the X values gauge the extent of such shift toward inhibition or activation.

CHESPA of the cGMP analog-bound PfDs reveals the presence of a third "mixed" conformer with an engaged pre-lid but a disengaged lid
The CHESPA of 8-NBD-cGMP shows that the majority of residues exhibit negative X values (Fig. 4B), suggestive of most residues shifting toward the inactive state upon binding of 8-NBD-cGMP. However, the magnitudes of the X values are highly residue-dependent, pointing to deviations from the pure inactiveactive two-state equilibrium, which predicts similar X values for residues outside the binding site. Another characteristic of a twostate model is the co-linearity of the perturbation and reference vectors (i.e. ͉cos()͉ ϳ 1) for residues outside the binding site. However, the CHESPA of 8-NBD-cGMP reveals that several residues not in the binding site exhibit ͉cos()͉ values significantly less than 1 (Fig. 4C), suggesting deviations from the two-state model and pointing to the presence of additional state(s) sampled by the PfD: 8-NBD-cGMP complex. Similar trends are observed for the 8-pCPT-cGMP-bound sample (Fig. 4, E and F).  (54). I, similar to H but for the C␣ chemical shifts.

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
The presence of multiple conformational equilibria is further supported by the width of the X value distributions (Fig. 4G). The X distributions measured for 8-NBD-cGMP-and 8-pCPT-cGMP-bound PfD are significantly wider compared with PET-cGMP ( Fig. 4G), which is nearly a full agonist and samples mainly the active state (Fig. S2). These observations suggest that the PfD:8-NBD-cGMP and PfD:8-pCPT-cGMP complexes sample a multistate ensemble, including not only the inactive and active conformers represented by the apo and cGMP-bound crystal structures, but also an additional state or states.
To simplify the analysis, we focused on the C-terminal helices that are directly linked to the catalytic domain and are one of the allosteric elements most critical in the inhibition and activation of the kinase (14,21,24). The C-terminal helices span two regions with clearly distinct average X values (͗X͘) in the PfD:8-NBD-cGMP complex: residues 515-530, denoted as the pre-lid with ͗X͘ ϭ Ϫ0.48, and residues 534 -535, denoted as the lid with ͗X͘ ϭ Ϫ0.88 (Fig. 4, B, H, and I). This means that for the PfD:8-NBD-cGMP complex, the pre-lid only shifts about halfway toward the inactive state, whereas the lid region shifts almost completely to the inactive state. A similar trend is observed for the PfD:8-pCPT-cGMP complex with average X values of Ϫ0.30 and Ϫ0.64 in the pre-lid and the lid regions, respectively.
One of the simplest models to explain both the enzymatic and the NMR data in Figs. 1L and 4, respectively, is a three-state equilibrium mechanism involving an inactive, an active, and an intermediate state (Fig. 5A). The inactive state is similar to the apo-PfD structure, where both the pre-lid and the lid are disengaged. The active state is similar to the cGMP-bound structure, where both the pre-lid and the lid are engaged. The intermediate state exhibits a mixed character, as the pre-lid is engaged similarly to the active state, but the lid is disengaged, similarly to the inactive state (Fig. 5A). Next, we checked whether it is possible to utilize this three-state model to back-calculate the relative kinase activities.

The CHESPA-based three-state model of PfD inhibition explains the experimental kinase activities
Using the average X values for the pre-lid and the lid regions (Fig. 4, B and E) and  Table S1). Similarly, the three state populations can be calculated for the 8-pCPT-cGMP-and PET-cGMP-bound PfD complexes based on the respective CHESPA data (Fig. 4E,  Fig. S2, and Table S1). Using these populations and assuming that the inactive and intermediate states are fully inhibitory,

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
given the lid disengagement (14), while the active state is noninhibitory, the relative kinase activity at ligand saturation induced by these cGMP analogs was predicted. The predicted activation values correlate well with the experimental kinase activity values (Fig. 5B).
The ability of our three-state model to recapitulate the enzymatic data has three key implications. First, the measurements on the PfD construct are relevant for the longer PfPKG construct used in the kinase assays to probe the PfPKG function. Second, the NMR chemical shift projection analysis and the resulting fractional shift values are suitable for predicting the populations of different functional states. Third, the mixed intermediate state, where only the lid region is disengaged, is inhibitory to an extent similar to the inactive state, further corroborating the pivotal role of the lid region in activation. Hence, we further probed the inhibitory mixed intermediate stabilized by the 8-NBD-cGMP analog.

8-NBD-cGMP binds to PfD in a syn conformation and perturbs interactions critical for activation
As a first step toward investigating the mixed intermediate, we examined the base orientation of 8-NBD-cGMP as bound to PfD. Typically, the base of cyclic nucleotides is oriented either in syn or anti relative to the ribose ring (34 -38). To determine whether the guanine base of 8-NBD-cGMP is syn or anti, we could not rely on NOE-based assessments, as this cGMP analog lacks nonexchangeable protons in the base. Hence, we focused on the chemical shifts of the C␥ 2 methyl of Ile 465 , which in the PfD:cGMP complex is in direct contact with the guanine base, and of the C␥ methyl of Thr 493 , whose hydroxyl interacts with the amino group of cGMP (Fig. 1C). The chemical shifts of both Ile 465 C␥ 2 and Thr 493 C␥ show only marginal changes when cGMP is replaced by 8-NBD-cGMP (Fig. 6, A and B), suggesting that the 8-NBD-cGMP guanine base binds in a syn conformation, similar to cGMP. This conclusion was confirmed also by

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
the 3D map of 8-NBD-cGMP versus cGMP NH CCS changes for residues in the ␤-barrel (Fig. 6C, red surface). Because residues in the ␤-barrel do not undergo significant structural changes upon cNMP binding, the perturbations sensed by the ⌬CCS report mainly on direct binding contacts (33,(39)(40)(41). The ⌬CCS map is fully consistent with 8-NBD-cGMP binding in syn rather than anti guanine orientation (Fig. 6C, residues within the red surface). Similarly, we confirmed through the ␤-barrel ⌬CCS map that a syn orientation is preserved also for the guanine base of PET-cGMP (Fig. 6D, blue surface). For this analog, the substitution is at the 1-and 2-rather than the 8-position (Fig. 1, I and K), causing the most significant CCS changes to cluster at the opposite side of the ␤-barrel for PET-cGMP compared with 8-NBD-cGMP, when both analogs maintain a syn guanine orientation (Fig. 6, C and D, blue versus red  surfaces).
Although 8-NBD-cGMP preserves a syn guanine orientation when bound to PfD, similar to cGMP, the introduction of the 8-NBD moiety has a profound effect on the conformational ensemble accessed by PfD, as shown by the NMR-based threestate model of Fig. 5A. Of the three PfD conformers in Fig. 5A (Fig.  5A), relying on the cGMP-bound coordinates (PDB code 4OGF) (14) for residues 401-532 and on the apo coordinates (PDB code 5DYK) (24) for residues 534 -542, which span the lid region. The hybrid model provided the initial coordinates for a 1-s MD trajectory, and Fig. 6 (E and F) shows representative structures selected based on the cluster analysis of such trajectory (Table S2).
In agreement with the data of Fig. 6 (A-C), Fig. 6E shows that the PfD-bound 8-NBD-cGMP preserves a syn guanine orientation. Fig. 6F illustrates that the lid region of the intermediate state is disengaged, whereas the pre-lid is engaged. Fig. 6G shows that when the intermediate conformation is overlaid with the cGMP-bound structure (PDB code 4OFG), the Arg 528 side chain of 4OFG clashes with the NBD moiety of 8-NBD-cGMP (Fig. 6G, gray sticks). The steric clash explains why in the simulated intermediate structure, the Arg 528 side chain moves away from the 8-NBD-cGMP and rotates outward, breaking the interaction with Tyr 480 (Fig. 6G, red sticks). Considering that the Tyr 480 -Arg 528 interaction is critical for activation (14,21), the perturbation of Arg 528 caused by 8-NBD is consistent with the lack of kinase activity observed upon replacement of cGMP with 8-NBD-cGMP.
This perturbation revealed by the MD simulations is consistent with the significant 8-NBD-cGMP versus cGMP chemical shift changes occurring at or near Tyr 480 (Fig. 3A). The Tyr 480 NH HSQC cross-peak of the cGMP:PfD complex moves completely toward the apo position when cGMP is replaced by 8-NBD-cGMP (Fig. 6H), confirming that the 8-NBD moiety perturbs the Tyr 480 -Arg 528 interaction. To further corroborate this result, we generated the R528K mutant and compared the R528K versus WT chemical shifts for both cGMP-and 8-NBD-cGMP-bound samples (Fig. 6I). For the cGMP:PfD complex, the largest R528K versus WT ⌬CCS values cluster primarily close to Tyr 480 (Fig. 6I), as expected based on the Arg 528 -Tyr 480 interaction observed in the cGMP:PfD structure (14). However, when cGMP is replaced by 8-NBD-cGMP, the R528K versus WT CCS changes in the vicinity of Tyr 480 are markedly suppressed (Fig. 6I), confirming that in the 8-NBD-cGMP:PfD complex, the Arg 528 -Tyr 480 interaction is significantly weakened compared with the cGMP:PfD structure.
Similar analyses were carried out to probe the effect of 8-NBD-cGMP on the capping triad interactions between Arg 484 , Gln 532 , and Asp 533 , which are also critical for activation (14). We found that the Arg 484 -Asp 533 interaction is perturbed in the simulated intermediate structure of the 8-NBD-cGMP: PfD complex (Fig. S3A). To confirm such perturbation, we measured the CCS changes between WT and the R484A mutant for  (Table  S1). Error bars for the experimental kinase activity represent the standard deviation of multiple measurements, and for the predicted kinase activity represent the errors propagated from the errors on the populations, which are determined from the standard deviation of the fractional shift (X) values of the residues used to calculate the average X values.

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
both the cGMP-and 8-NBD-cGMP-bound samples (Fig.  S3B). Fig. S3B reveals that the R484A versus WT CCS changes in the C-terminal helices, especially in the vicinity of Gln 532 -Asp 533 , are suppressed when cGMP is replaced with 8-NBD-cGMP, confirming that 8-NBD-cGMP perturbs the capping interaction between Arg 484 and the C-terminal helices. These results are in full agreement with the MD simulations, suggesting that the 8-NBD substituent destabilizes also the capping triad.
Hence, we conclude that upon introduction of the NBD group at the 8-guanine position, the syn base orientation is pre-served but the Tyr 480 -Arg 528 contact and the capping triad interactions are significantly perturbed, thus compromising the full engagement of the C-helix and explaining why 8-NBD-cGMP functions as an effective allosteric antagonist of the PfPKG kinase.

Sampling a mixed intermediate state is beneficial for maximizing kinase inhibition while minimizing affinity losses
The three-state model with the mixed intermediate state featuring an engaged pre-lid and a disengaged lid (Fig. 5A) pro-

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
vides an excellent strategy for PfPKG kinase inhibition without significantly compromising affinity. Reconciling high inhibitory efficacy and potency would be more challenging for a simple two-state model. In this case, the disengagement of the lid necessary for inhibition also implies the concurrent disengagement of other PfD elements necessary for high-affinity binding, such as the PBC in its entirety. On the contrary, in the mixed intermediate state (Fig. 5A), which exhibits an inhibitory effect similar to that of the inactive state, it maintains the pre-lid and the adjacent PBC engaged. For example, whereas the capping carboxylate of Asp 533 is displaced away from Arg 484 in the PBC, the side chain of this arginine in the mixed intermediate remains in an orientation similar to the active conformation (Fig. S3A), thus minimizing perturbations of the interactions with the cyclic nucleotide. Hence, based on the three-state model of Fig. 5A, we hypothesize that the 8-NBD-cGMP preserves an affinity comparable with the endogenous allosteric effector, cGMP.
Our hypothesis is confirmed through surface plasmon resonance (SPR) measurements (Fig. 7) showing that the affinities of 8-NBD-cGMP and cGMP are comparable (i.e. K d ϭ 59 Ϯ 4 nM versus 51 Ϯ 7 nM; Table S3), despite their markedly different efficacies in terms of kinase activation ( Fig. 1L and Fig. S1A). This conservation of affinity is especially important, because 8-NBD-cGMP acts as a competitive antagonist (Fig. S1B). Hence, 8-NBD-cGMP serves as an excellent example for how sampling a mixed intermediate state enables significant inhibition without major affinity losses by selectively disengaging only structural units strictly necessary for kinase activation. Similar inhibition strategies are anticipated to inform the design of future PfPKG allosteric inhibitors, which may serve as anti-malaria drug leads.

Expression and purification
The PfPKG CBD-D construct (residues 401-542) (i.e. PfD) within a His-tagged expression vector pQTEV (14) was expressed in the Escherichia coli strain BL21 (DE3). The transformed E. coli was grown in a minimal medium enriched with 15 NH 4 Cl and [ 13 C]glucose and supplemented with D-biotin, thiamine HCl, MgSO 4 , and CaCl 2 . The cells were grown at 37°C to an optical density of 0.7 (at ϭ 600 nm), induced with 0.5 mM isopropyl 1-thio-␤-D-galactopyranoside, and then incubated for an additional 18 h at 18°C. Cells were harvested, resuspended in the lysis buffer (50 mM Tris, pH 7.4, 500 mM NaCl, 1 mM BME), and lysed using a cell disruptor. The cell debris was removed through centrifugation, and then the supernatant was loaded onto an Ni 2ϩ -Sepharose resin. The cell lysate was allowed to flow through, and the column was rinsed with 50 mM Tris, pH 7.4, 500 mM NaCl, 1 mM BME, and 20 mM imidazole. After thorough rinsing, the protein was eluted with 50 mM Tris, pH 7.4, 500 mM NaCl, 1 mM BME, and 300 mM imidazole. The eluted protein was then dialyzed in 50 mM Tris, pH 7.4, 100 mM NaCl, 1 mM BME with tobacco etch virus protease for 24 h at 4°C. The cleaved products were loaded on the Ni 2ϩ -Sepharose resin. The protein without the His-tag was obtained through collecting the flow-through and was further purified with a HiLoad 16/60 Superdex 75 gel filtration column (GE Healthcare) equilibrated with NMR buffer (50 mM Tris, pH 7, 100 mM NaCl, 10 mM MgCl 2 , 1 mM DTT). The His-tagged PfPKG 401-853 construct including the catalytic domain was expressed in E. coli TP2000 and purified using Ni-NTA and size-exclusion chromatography, as described previously (21).

NMR spectroscopy
NMR experiments were performed on a Bruker Avance 700-MHz spectrometer equipped with a 5-mm TCI cryoprobe or a Bruker Avance 850-MHz spectrometer equipped with a tripleresonance TXI probe. All experiments were acquired at 306 K in NMR buffer (50 mM Tris, pH 7, 100 mM NaCl, 10 mM MgCl 2 , 1 mM DTT) with 5% D 2 O. All NMR spectra were processed with NMRPipe (42) and analyzed using Sparky (43). Three-dimensional triple-resonance NMR experiments (i.e. HNCACB, CBCA(CO)NH, HNCA, HN(CO)CA, HNCO, HCCH-TOCSY, and (H)CC(CO)NH) were acquired and used to assign the twodimensional 1 H-15 N HSQC and 1 H-13 C HSQC NMR spectra of cGMP-bound PfD. Because the apo sample was prone to precipitation, 2D N z -exchange, 2D difference N z -exchange, and 2D C z -exchange experiments with mixing time of 250 ms were used to transfer the chemical shift assignments from the cGMP-bound PfD to the apo-PfD spectra. The cGMP analogbound spectra were assigned through spectral comparison with  Table S3. A Langmuir 1:1 model was applied for fitting.

EDITORS' PICK: Mechanism of PfPKG allosteric inhibition
cGMP-bound spectra, if no ambiguities were present, and through 2D N z exchange and 2D difference N z exchange with the apo state. Chemical shift assignments of R528K:cGMP and R528K:8-NBD-cGMP were obtained through spectral comparison with WT spectra if no ambiguities were present.

Chemical shift analyses
All spectra utilized for chemical shift analyses were acquired with 100 M PfD and saturating amounts of cGMP or cGMP analogs (2 mM) or without any cNMP addition for apo spectra. For accurate chemical shift comparisons, an internal reference (i.e. 15 N-labeled N-acetylglycine) was used to align the spectra in Sparky, and the resulting chemical shift values were used to compute the compounded chemical shift differences and perform CHESPA following protocols described previously (29,30,34,44,45). The fractional shift (X) values and cos() values were computed for residues with cGMP analog-bound versus cGMP-bound CCS difference Ͼ0.05 ppm. The average ͗X͘ values for the pre-lid and the lid were calculated using residues that were assigned in all five samples (i.e. apo, cGMP-bound, 8-NBD-cGMP-bound, 8-pCPT-cGMP-bound, and PET-cGMP-bound PfD), had X values between 0 and Ϫ1, and had ⌬CCS (cGMP analog-bound versus cGMP-bound) values greater than 0.05 ppm, ensuring reliable X and cos() values. The distribution of X values for cGMP analogs was computed using the residues that exhibited ͉cos()͉ Ͼ 0.75 to minimize biases from nearest-neighbor effects.

Site-directed mutagenesis
The R528K mutation was created through site-directed mutagenesis using a KOD Hot Start Master mix (Novagen) and primer pairs (synthesized by IDT): forward, CTGGCCCACTT-GGAAGAGAAGATTAAGATGCAGGATAC; reverse, GTAT-CCTGCATCTTAATCTTCTCTTCCAAGTGGGCCAG. The reaction protocol provided by the KOD Hot Start kit was followed. Template DNA was degraded with DpnI for 2 h at 37°C, and the reaction product was used for transformation of E. coli top10 cells. The transformed colony was grown in lysogeny broth or Luria-Bertani medium medium for 18 h at 37°C. Plasmids were isolated from the E. coli using the GeneJET Plasmid Miniprep Kit (Thermo Fisher Scientific). The purified plasmids were sequenced to verify the mutation. The expression and purification of PfD R528K were carried out using the same protocol as for the WT PfD.

Kinase assays
Specific kinase activities were determined using an enzymecoupled spectrophotometric assay (46). In a quartz cuvette, purified PfPKG 401-853 was mixed with the reaction buffer (100 mM MOPS, pH 7.0, 10 mM MgCl 2 , 5 mM ␤-mercaptoethanol, 1 mM ATP, 1 mM phosphoenolpyruvate, 230 M NADH, 15 units/ml lactate dehydrogenase, 8.4 units/ml pyruvate kinase). The reactions were started by adding 1 mM PKStide (GRTGRRNSI; GeneCust, Luxembourg) and monitored in the absence and presence of a 100 M concentration of each cyclic nucleotide (cGMP, 8-pCPT-cGMP, and PET-cGMP) for at least 1 min. Specific kinase activity was calculated according to the Lambert-Beer law.
Kinase activity induced by 8-NBD-cGMP was measured by microfluidic mobility shift assay using a Caliper DeskTop Profiler (PerkinElmer Life Sciences). Here, 100 M cGMP or 8-NBD-cGMP was mixed with assay buffer (20 mM MOPS (pH 7), 300 mM NaCl, 10 mM MgCl 2 , 1 mM ATP, 1 mM DTT, 990 M PKStide, 10 M FITC-PKStide (FITC-GRTGRRNSI; GeneCust, Luxembourg, 0.05% L-31) and a fixed protein concentration and continued as described previously (21). The assay was also conducted with multiple concentrations of 8-NBD-cGMP or cGMP, ranging from 0.1 nM to 300 M. For the competition assay, the protein was preincubated with 350 nM cGMP, varying concentrations of 8-NBD-cGMP were added to the reaction, and substrate conversion was monitored over 2 h. The combined kinase inhibition data were further analyzed and normalized using GraphPad Prism 8.3.0 (GraphPad Software).

SPR
Direct binding data were obtained via SPR using a Biacore T200 (GE Healthcare Life Sciences). Briefly, 200 -300 RU of Ni 2ϩ (5 mM NiCl 2 in running buffer) were immobilized on a poly-NTA-derivatized sensor chip (NIHC 1000M; XanTec Bioanalytics GmbH). Subsequently, His-tagged PfPKG 401-853 in running buffer (20 mM HEPES, 150 mM NaCl, 50 M EDTA, 0.01% P20, pH 7.3) was captured at a flow rate of 10 l/min to an immobilization level of 500 -2500 RU. Next, association was measured by injection of a dilution series of the respective cyclic nucleotide for 45-75 s. Then running buffer was injected for 60 -180 s to initiate dissociation. Regeneration of the sensor chip surface was achieved by consecutively injecting 0.5 M EDTA (pH 8) for 600 s, followed by injections of 3 M guanidinium HCl, 0.5% SDS, and 3 M guanidinium HCl for 240 s each. Unless stated otherwise, all kinetic measurements were performed at 25°C in running buffer and a flow rate of 50 l/min. Data were evaluated with the Biacore T200 Evaluation Software 3.0 using a global fit and assuming a 1:1 Langmuir binding model. Further data processing was performed using GraphPad Prism 8.3.0 (GraphPad Software).

Molecular dynamic simulations: initial model preparation
The structure of 8-NBD-cGMP in a syn conformation was generated using RDKit (47). A model of the mixed intermediate state sampled by the PfD:8-NBD-cGMP complex was then built using residues 401-533 from the active PfD:cGMP crystal structure (PDB code 4OFG) (14) and lid region residues 534 -542 from the inactive apo-PfD crystal structure (PDB code 5DYK) (24). The 8-NBD-cGMP ligand was docked into the resulting mixed structure of PfD using AutoDock4 with default parameters generated by AutoDockTools (48). Residues 463-494, where the PBC and BBR regions are located, were assigned as the binding pocket for AutoDock. The structure in which the ribose-phosphate and guanine moieties of the docked ligand were best aligned with the position of cGMP in the cGMPbound active structure (PDB code 4OFG) (14) was selected as the initial coordinates for MD simulations.

Molecular dynamics simulation protocol
All simulations were performed using the AMBER 16 software PMEMD.CUDA (49) on the Shared Hierarchical Aca-EDITORS' PICK: Mechanism of PfPKG allosteric inhibition demic Research Computing Network (SHARCNET). The Amber ff99SBnmr force field was employed for the protein, whereas RESP charge and the general Amber force field were applied to generate the parameters for 8-NBD-cGMP using the Gaussian09 program (50) and the antechamber module implemented in AMBER 16, as described previously (51). The TIP3P water model was applied in the form of a rectangular MD solvent box surrounding the PfD:8-NBD-cGMP structure, with a minimum distance of 12 Å between the PfD:8-NBD-cGMP structure and the box edge. To mimic the experimental conditions, the charged residues and N/C termini were set to the protonation states expected at pH 7, following a protocol similar to that implemented previously (52), and Na ϩ and Cl Ϫ ions were added to mimic a 100 mM salt concentration.
The simulations started with an energy minimization performed with harmonic restraints on the protein backbone (with a force constant of 5 kcal/mol⅐Å 2 ). The system was then gradually heated from 0 to 100 K at constant volume and then heated from 100 to 306 K in the NPT ensemble, and all of the heating process employed harmonic restraints on protein backbone (with a force constant of 3 kcal/mol⅐Å 2 ) for a total of 1 ns followed by an equilibration period of 1.2 ns with main-chain atoms restrained with a force constant of 1 kcal/mol⅐Å 2 . Then a 1-s production MD trajectory was generated in the NPT ensemble at 306 K and 1 atm with a weak-coupling algorithm, saving structures every 10 ps for subsequent analysis. A cluster analysis for 20,000 frames recorded from the MD trajectory (i.e. the structures recorded at 50-ps intervals) was performed, with ␤-core alignment, using the CPPTRAJ protocol (53).

Data availability
The NMR chemical shift assignment for cGMP-bound PfD has been deposited in the Biological Magnetic Resonance Data Bank (BMRB) with the code 50203. All remaining data are contained within the article.