Structure and dynamics of a pentameric KCTD5/CUL3/Gßγ E3 ubiquitin ligase complex.

Heterotrimeric G proteins can be regulated by posttranslational modifications, including ubiquitylation. KCTD5, a pentameric substrate receptor protein consisting of an N-terminal BTB domain and a C-terminal domain, engages CUL3 to form the central scaffold of a cullin-RING E3 ligase complex (CRL3KCTD5) that ubiquitylates Gßγ and reduces Gßγ protein levels in cells. The cryo-EM structure of a 5:5:5 KCTD5/CUL3NTD/Gß1γ2 assembly reveals a highly dynamic complex with rotations of over 60° between the KCTD5BTB/CUL3NTD and KCTD5CTD/Gßγ moieties of the structure. CRL3KCTD5 engages the E3 ligase ARIH1 to ubiquitylate Gßγ in an E3-E3 superassembly, and extension of the structure to include full-length CUL3 with RBX1 and an ARIH1~ubiquitin conjugate reveals that some conformational states position the ARIH1~ubiquitin thioester bond to within 10 Å of lysine-23 of Gß and likely represent priming complexes. Most previously described CRL/substrate structures have consisted of monovalent complexes and have involved flexible peptide substrates. The structure of the KCTD5/CUL3NTD/Gßγ complex shows that the oligomerization of a substrate receptor can generate a polyvalent E3 ligase complex and that the internal dynamics of the substrate receptor can position a structured target for ubiquitylation in a CRL3 complex.

SARS-CoV-2 impairs interferon production via NSP2-induced repression of mRNA translation.

Viruses evade the innate immune response by suppressing the production or activity of cytokines such as type I interferons (IFNs). Here we report the discovery of a mechanism by which the SARS-CoV-2 virus coopts an intrinsic cellular machinery to suppress the production of the key immunostimulatory cytokine IFN-ß. We reveal that the SARS-CoV-2 encoded nonstructural protein 2 (NSP2) directly interacts with the cellular GIGYF2 protein. This interaction enhances the binding of GIGYF2 to the mRNA cap-binding protein 4EHP, thereby repressing the translation of the Ifnb1 mRNA. Depletion of GIGYF2 or 4EHP significantly enhances IFN-ß production, which inhibits SARS-CoV-2 replication. Our findings reveal a target for rescuing the antiviral innate immune response to SARS-CoV-2 and other RNA viruses.

Atomic model for the membrane-embedded VO motor of a eukaryotic V-ATPase.

Vacuolar-type ATPases (V-ATPases) are ATP-powered proton pumps involved in processes such as endocytosis, lysosomal degradation, secondary transport, TOR signalling, and osteoclast and kidney function. ATP hydrolysis in the soluble catalytic V1 region drives proton translocation through the membrane-embedded VO region via rotation of a rotor subcomplex. Variability in the structure of the intact enzyme has prevented construction of an atomic model for the membrane-embedded motor of any rotary ATPase. We induced dissociation and auto-inhibition of the V1 and VO regions of the V-ATPase by starving the yeast Saccharomyces cerevisiae, allowing us to obtain a ~3.9-Šresolution electron cryomicroscopy map of the VO complex and build atomic models for the majority of its subunits. The analysis reveals the structures of subunits ac8c'c″de and a protein that we identify and propose to be a new subunit (subunit f). A large cavity between subunit a and the c-ring creates a cytoplasmic half-channel for protons. The c-ring has an asymmetric distribution of proton-carrying Glu residues, with the Glu residue of subunit c″ interacting with Arg735 of subunit a. The structure suggests sequential protonation and deprotonation of the c-ring, with ATP-hydrolysis-driven rotation causing protonation of a Glu residue at the cytoplasmic half-channel and subsequent deprotonation of a Glu residue at a luminal half-channel.

Structural Characterization of Endogenous Tuberous Sclerosis Protein Complex Revealed Potential Polymeric Assembly.

Tuberous sclerosis protein complex (pTSC) nucleates a proteinaceous signaling hub that integrates information about the internal and external energy status of the cell in the regulation of growth and energy consumption. Biochemical and cryo-electron microscopy studies of recombinant pTSC have revealed its structure and stoichiometry and hinted at the possibility that the complex may form large oligomers. Here, we have partially purified endogenous pTSC from fasted mammalian brains of rat and pig by leveraging a recombinant antigen binding fragment (Fab) specific for the TSC2 subunit of pTSC. We demonstrate Fab-dependent purification of pTSC from membrane-solubilized fractions of the brain homogenates. Negative stain electron microscopy of the samples purified from pig brain demonstrates rod-shaped protein particles with a width of 10 nm, a variable length as small as 40 nm, and a high degree of conformational flexibility. Larger filaments are evident with a similar 10 nm width and a ≤1 µm length in linear and weblike organizations prepared from pig brain. Immunogold labeling experiments demonstrate linear aggregates of pTSC purified from mammalian brains. These observations suggest polymerization of endogenous pTSC into filamentous superstructures.

Oncogenic and RASopathy-associated K-RAS mutations relieve membrane-dependent occlusion of the effector-binding site.

K-RAS4B (Kirsten rat sarcoma viral oncogene homolog 4B) is a prenylated, membrane-associated GTPase protein that is a critical switch for the propagation of growth factor signaling pathways to diverse effector proteins, including rapidly accelerated fibrosarcoma (RAF) kinases and RAS-related protein guanine nucleotide dissociation stimulator (RALGDS) proteins. Gain-of-function KRAS mutations occur frequently in human cancers and predict poor clinical outcome, whereas germ-line mutations are associated with developmental syndromes. However, it is not known how these mutations affect K-RAS association with biological membranes or whether this impacts signal transduction. Here, we used solution NMR studies of K-RAS4B tethered to nanodiscs to investigate lipid bilayer-anchored K-RAS4B and its interactions with effector protein RAS-binding domains (RBDs). Unexpectedly, we found that the effector-binding region of activated K-RAS4B is occluded by interaction with the membrane in one of the NMR-observable, and thus highly populated, conformational states. Binding of the RAF isoform ARAF and RALGDS RBDs induced marked reorientation of K-RAS4B from the occluded state to RBD-specific effector-bound states. Importantly, we found that two Noonan syndrome-associated mutations, K5N and D153V, which do not affect the GTPase cycle, relieve the occluded orientation by directly altering the electrostatics of two membrane interaction surfaces. Similarly, the most frequent KRAS oncogenic mutation G12D also drives K-RAS4B toward an exposed configuration. Further, the D153V and G12D mutations increase the rate of association of ARAF-RBD with lipid bilayer-tethered K-RAS4B. We revealed a mechanism of K-RAS4B autoinhibition by membrane sequestration of its effector-binding site, which can be disrupted by disease-associated mutations. Stabilizing the autoinhibitory interactions between K-RAS4B and the membrane could be an attractive target for anticancer drug discovery.

Biochemical Classification of Disease-associated Mutants of RAS-like Protein Expressed in Many Tissues (RIT1).

RAS-like protein expressed in many tissues 1 (RIT1) is a disease-associated RAS subfamily small guanosine triphosphatase (GTPase). Recent studies revealed that germ-line and somatic RIT1 mutations can cause Noonan syndrome (NS), and drive proliferation of lung adenocarcinomas, respectively, akin to RAS mutations in these diseases. However, the locations of these RIT1 mutations differ significantly from those found in RAS, and do not affect the three mutational "hot spots" of RAS. Moreover, few studies have characterized the GTPase cycle of RIT1 and its disease-associated mutants. Here we developed a real-time NMR-based GTPase assay for RIT1 and investigated the effect of disease-associated mutations on GTPase cycle. RIT1 exhibits an intrinsic GTP hydrolysis rate similar to that of H-RAS, but its intrinsic nucleotide exchange rate is ∼4-fold faster, likely as a result of divergent residues near the nucleotide binding site. All of the disease-associated mutations investigated increased the GTP-loaded, activated state of RIT1 in vitro, but they could be classified into two groups with different intrinsic GTPase properties. The S35T, A57G, and Y89H mutants exhibited more rapid nucleotide exchange, whereas F82V and T83P impaired GTP hydrolysis. A RAS-binding domain pulldown assay indicated that RIT1 A57G and Y89H were highly activated in HEK293T cells, whereas T83P and F82V exhibited more modest activation. All five mutations are associated with NS, whereas two (A57G and F82V) have also been identified in urinary tract cancers and myeloid malignancies. Characterization of the effects on the GTPase cycle of RIT1 disease-associated mutations should enable better understanding of their role in disease processes.

Structure-guided mutation of the conserved G3-box glycine in Rheb generates a constitutively activated regulator of mammalian target of rapamycin (mTOR).

Constitutively activated variants of small GTPases, which provide valuable functional probes of their role in cellular signaling pathways, can often be generated by mutating the canonical catalytic residue (e.g. Ras Q61L) to impair GTP hydrolysis. However, this general approach is ineffective for a substantial fraction of the small GTPase family in which this residue is not conserved (e.g. Rap) or not catalytic (e.g. Rheb). Using a novel engineering approach, we have manipulated nucleotide binding through structure-guided substitutions of an ultraconserved glycine residue in the G3-box motif (DXXG). Substitution of Rheb Gly-63 with alanine impaired both intrinsic and TSC2 GTPase-activating protein (GAP)-mediated GTP hydrolysis by displacing the hydrolytic water molecule, whereas introduction of a bulkier valine side chain selectively blocked GTP binding by steric occlusion of the γ-phosphate. Rheb G63A stimulated phosphorylation of the mTORC1 substrate p70S6 kinase more strongly than wild-type, thus offering a new tool for mammalian target of rapamycin (mTOR) signaling.

Direct structural analysis of a single acyl carrier protein domain in fatty acid synthase from the fungus Saccharomyces cerevisiae.

Acyl carrier protein (ACP) is the work horse of polyketide (PKS) and fatty acid synthases (FAS) and acts as a substrate shuttling domain in these mega enzymes. In fungi, FAS forms a 2.6 MDa symmetric assembly with six identical copies of FAS1 and FAS2 polypeptides. However, ACP spatial distribution is not restricted by symmetry owing to the long and flexible loops that tether the shuttling domain to its corresponding FAS2 polypeptide. This symmetry breaking has hampered experimental investigation of substrate shuttling route in fungal FAS. Here, we develop a protein engineering and expression method to isolate asymmetric fungal FAS proteins containing odd numbers of ACP domains. Electron cryomicroscopy (cryoEM) observation of the engineered complex reveals a non-uniform distribution of the substrate shuttling domain relative to its corresponding FAS2 polypeptide at 2.9 Å resolution. This work lays the methodological foundation for experimental study of ACP shuttling route in fungi.

The HisRS-like domain of GCN2 is a pseudoenzyme that can bind uncharged tRNA.

GCN2 is a stress response kinase that phosphorylates the translation initiation factor eIF2α to inhibit general protein synthesis when activated by uncharged tRNA and stalled ribosomes. The presence of a HisRS-like domain in GCN2, normally associated with tRNA aminoacylation, led to the hypothesis that eIF2α kinase activity is regulated by the direct binding of this domain to uncharged tRNA. Here we solved the structure of the HisRS-like domain in the context of full-length GCN2 by cryoEM. Structure and function analysis shows the HisRS-like domain of GCN2 has lost histidine and ATP binding but retains tRNA binding abilities. Hydrogen deuterium exchange mass spectrometry, site-directed mutagenesis and computational docking experiments support a tRNA binding model that is partially shifted from that employed by bona fide HisRS enzymes. These results demonstrate that the HisRS-like domain of GCN2 is a pseudoenzyme and advance our understanding of GCN2 regulation and function.

Membrane-dependent modulation of the mTOR activator Rheb: NMR observations of a GTPase tethered to a lipid-bilayer nanodisc.

Like most Ras superfamily proteins, the GTPase domain of Ras homologue enriched in brain (Rheb) is tethered to cellular membranes through a prenylated cysteine in a flexible C-terminal region; however, little is known about how Rheb or other GTPases interact with the membrane or how this environment may affect their GTPase functions. We used NMR methods to characterize Rheb tethered to nanodiscs, monodisperse protein-encapsulated lipid bilayers with a diameter of 10 nm. Membrane conjugation markedly reduced the rate of intrinsic nucleotide exchange, while GTP hydrolysis was unchanged. NMR measurements revealed that the GTPase domain interacts transiently with the surface of the bilayer in two distinct preferred orientations, which are determined by the bound nucleotide. We propose models of membrane-dependent signal regulation by Rheb that shed light on previously unexplained in vivo properties of this GTPase. The study presented provides a general approach for direct experimental investigation of membrane-dependent properties of other Ras-superfamily GTPases.

Probing the GTPase cycle with real-time NMR: GAP and GEF activities in cell extracts.

The Ras superfamily of small GTPases is a large family of switch-like proteins that control diverse cellular functions, and their deregulation is associated with multiple disease processes. When bound to GTP they adopt a conformation that interacts with effector proteins, whereas the GDP-bound state is generally biologically inactive. GTPase activating proteins (GAPs) promote hydrolysis of GTP, thus impeding the biological activity of GTPases, whereas guanine nucleotide exchange factors (GEFs) promote exchange of GDP for GTP and activate GTPase proteins. A number of methods have been developed to assay GTPase nucleotide hydrolysis and exchange, as well as the activity of GAPs and GEFs. The kinetics of these reactions are often studied with purified proteins and fluorescent nucleotide analogs, which have been shown to non-specifically impact hydrolysis and exchange. Most GAPs and GEFs are large multidomain proteins subject to complex regulation that is challenging to reconstitute in vitro. In cells, the activities of full-length GAPs or GEFs are typically assayed indirectly on the basis of nucleotide loading of the cognate GTPase, or by exploiting their interaction with effector proteins. Here, we describe a recently developed real-time NMR method to assay kinetics of nucleotide exchange and hydrolysis reactions by direct monitoring of nucleotide-dependent structural changes in an isotopically labeled GTPase. The unambiguous readout of this method makes it possible to precisely measure GAP and GEF activities from extracts of mammalian cells, enabling studies of their catalytic and regulatory mechanisms. We present examples of NMR-based assays of full-length GAPs and GEFs overexpressed in mammalian cells.

Structure of the immunoregulatory sialidase NEU1.

Sialic acids linked to glycoproteins and glycolipids are important mediators of cell and protein recognition events. These sugar residues are removed by neuraminidases (sialidases). Neuraminidase-1 (sialidase-1 or NEU1) is a ubiquitously expressed mammalian sialidase located in lysosomes and on the cell membrane. Because of its modulation of multiple signaling processes, it is a potential therapeutic target for cancers and immune disorders. Genetic defects in NEU1 or in its protective protein cathepsin A (PPCA, CTSA) cause the lysosomal storage diseases sialidosis and galactosialidosis. To further our understanding of this enzyme's function at the molecular level, we determined the three-dimensional structure of murine NEU1. The enzyme oligomerizes through two self-association interfaces and displays a wide substrate-binding cavity. A catalytic loop adopts an inactive conformation. We propose a mechanism of activation involving a conformational change in this loop upon binding to its protective protein. These findings may facilitate the development of selective inhibitor and agonist therapies.

Atomic model for core modifying region of human fatty acid synthase in complex with Denifanstat.

Fatty acid synthase (FASN) catalyzes the de novo synthesis of palmitate, a 16-carbon chain fatty acid that is the primary precursor of lipid metabolism and an important intracellular signaling molecule. FASN is an attractive drug target in diabetes, cancer, fatty liver diseases, and viral infections. Here, we develop an engineered full-length human FASN (hFASN) that enables isolation of the condensing and modifying regions of the protein post-translation. The engineered protein enables electron cryo-microscopy (cryoEM) structure determination of the core modifying region of hFASN to 2.7 Å resolution. Examination of the dehydratase dimer within this region reveals that unlike its close homolog, porcine FASN, the catalytic cavity is close-ended and is accessible only through one opening in the vicinity of the active site. The core modifying region exhibits two major global conformational variabilities that describe long-range bending and twisting motions of the complex in solution. Finally, we solved the structure of this region bound to an anti-cancer drug, Denifanstat (i.e., TVB-2640), demonstrating the utility of our approach as a platform for structure guided design of future hFASN small molecule inhibitors.

Identification of triazenyl indoles as inhibitors of fungal fatty acid biosynthesis with broad-spectrum activity.

Rising drug resistance among pathogenic fungi, paired with a limited antifungal arsenal, poses an increasing threat to human health. To identify antifungal compounds, we screened the RIKEN natural product depository against representative isolates of four major human fungal pathogens. This screen identified NPD6433, a triazenyl indole with broad-spectrum activity against all screening strains, as well as the filamentous mold Aspergillus fumigatus. Mechanistic studies indicated that NPD6433 targets the enoyl reductase domain of fatty acid synthase 1 (Fas1), covalently inhibiting its flavin mononucleotide-dependent NADPH-oxidation activity and arresting essential fatty acid biosynthesis. Robust Fas1 inhibition kills Candida albicans, while sublethal inhibition impairs diverse virulence traits. At well-tolerated exposures, NPD6433 extended the lifespan of nematodes infected with azole-resistant C. albicans. Overall, identification of NPD6433 provides a tool with which to explore lipid homeostasis as a therapeutic target in pathogenic fungi and reveals a mechanism by which Fas1 function can be inhibited.

The projection analysis of NMR chemical shifts reveals extended EPAC autoinhibition determinants.

EPAC is a cAMP-dependent guanine nucleotide exchange factor that serves as a prototypical molecular switch for the regulation of essential cellular processes. Although EPAC activation by cAMP has been extensively investigated, the mechanism of EPAC autoinhibition is still not fully understood. The steric clash between the side chains of two conserved residues, L273 and F300 in EPAC1, has been previously shown to oppose the inactive-to-active conformational transition in the absence of cAMP. However, it has also been hypothesized that autoinhibition is assisted by entropic losses caused by quenching of dynamics that occurs if the inactive-to-active transition takes place in the absence of cAMP. Here, we test this hypothesis through the comparative NMR analysis of several EPAC1 mutants that target different allosteric sites of the cAMP-binding domain (CBD). Using what to our knowledge is a novel projection analysis of NMR chemical shifts to probe the effect of the mutations on the autoinhibition equilibrium of the CBD, we find that whenever the apo/active state is stabilized relative to the apo/inactive state, dynamics are consistently quenched in a conserved loop (ß2-ß3) and helix (α5) of the CBD. Overall, our results point to the presence of conserved and nondegenerate determinants of CBD autoinhibition that extends beyond the originally proposed L273/F300 residue pair, suggesting that complete activation necessitates the simultaneous suppression of multiple autoinhibitory mechanisms, which in turn confers added specificity for the cAMP allosteric effector.

The Trisubstituted Isoxazole MMV688766 Exerts Broad-Spectrum Activity against Drug-Resistant Fungal Pathogens through Inhibition of Lipid Homeostasis.

Candida species are among the most prevalent causes of systemic fungal infection, posing a growing threat to public health. While Candida albicans is the most common etiological agent of systemic candidiasis, the frequency of infections caused by non-albicans Candida species is rising. Among these is Candida auris, which has emerged as a particular concern. Since its initial discovery in 2009, it has been identified worldwide and exhibits resistance to all three principal antifungal classes. Here, we endeavored to identify compounds with novel bioactivity against C. auris from the Medicines for Malaria Venture's Pathogen Box library. Of the five hits identified, the trisubstituted isoxazole MMV688766 emerged as the only compound displaying potent fungicidal activity against C. auris, as well as other evolutionarily divergent fungal pathogens. Chemogenomic profiling, as well as subsequent metabolomic and phenotypic analyses, revealed that MMV688766 disrupts cellular lipid homeostasis, driving a decrease in levels of early sphingolipid intermediates and fatty acids and a concomitant increase in lysophospholipids. Experimental evolution to further probe MMV688766's mode of action in the model fungus Saccharomyces cerevisiae revealed that loss of function of the transcriptional regulator HAL9 confers resistance to MMV688766, in part through the upregulation of the lipid-binding chaperone HSP12, a response that appears to assist in tolerating MMV688766-induced stress. The novel mode of action we have uncovered for MMV688766 against drug-resistant fungal pathogens highlights the broad utility of targeting lipid homeostasis to disrupt fungal growth and how screening structurally-diverse chemical libraries can provide new insights into resistance-conferring stress responses of fungi. IMPORTANCE As widespread antimicrobial resistance threatens to propel the world into a postantibiotic era, there is a pressing need to identify mechanistically distinct antimicrobial agents. This is of particular concern when considering the limited arsenal of drugs available to treat fungal infections, coupled with the emergence of highly drug-resistant fungal pathogens, including Candida auris. In this work, we demonstrate that existing libraries of drug-like chemical matter can be rich resources for antifungal molecular scaffolds. We discovered that the small molecule MMV688766, from the Pathogen Box library, displays previously undescribed broad-spectrum fungicidal activity through perturbation of lipid homeostasis. Characterization of the mode of action of MMV688766 provided new insight into the protective mechanisms fungi use to cope with the disruption of lipid homeostasis. Our findings highlight that elucidating the genetic circuitry required to survive in the presence of cellular stress offers powerful insights into the biological pathways that govern this important phenotype.

Real-time NMR study of three small GTPases reveals that fluorescent 2'(3')-O-(N-methylanthraniloyl)-tagged nucleotides alter hydrolysis and exchange kinetics.

The Ras family of small GTPases control diverse signaling pathways through a conserved "switch" mechanism, which is turned on by binding of GTP and turned off by GTP hydrolysis to GDP. Full understanding of GTPase switch functions requires reliable, quantitative assays for nucleotide binding and hydrolysis. Fluorescently labeled guanine nucleotides, such as 2'(3')-O-(N-methylanthraniloyl) (mant)-substituted GTP and GDP analogs, have been widely used to investigate the molecular properties of small GTPases, including Ras and Rho. Using a recently developed NMR method, we show that the kinetics of nucleotide hydrolysis and exchange by three small GTPases, alone and in the presence of their cognate GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors, are affected by the presence of the fluorescent mant moiety. Intrinsic hydrolysis of mantGTP by Ras homolog enriched in brain (Rheb) is approximately 10 times faster than that of GTP, whereas it is 3.4 times slower with RhoA. On the other hand, the mant tag inhibits TSC2GAP-catalyzed GTP hydrolysis by Rheb but promotes p120 RasGAP-catalyzed GTP hydrolysis by H-Ras. Guanine nucleotide exchange factor-catalyzed nucleotide exchange for both H-Ras and RhoA was inhibited by mant-substituted nucleotides, and the degree of inhibition depends highly on the GTPase and whether the assay measures association of mantGTP with, or dissociation of mantGDP from the GTPase. These results indicate that the mant moiety has significant and unpredictable effects on GTPase reaction kinetics and underscore the importance of validating its use in each assay.

Structure of the murine lysosomal multienzyme complex core.

The enzymes ß-galactosidase (GLB1) and neuraminidase 1 (NEU1; sialidase 1) participate in the degradation of glycoproteins and glycolipids in the lysosome. To remain active and stable, they associate with PPCA [protective protein cathepsin A (CTSA)] into a high-molecular weight lysosomal multienzyme complex (LMC), of which several forms exist. Genetic defects in these three proteins cause the lysosomal storage diseases GM1-gangliosidosis/mucopolysaccharidosis IV type B, sialidosis, and galactosialidosis, respectively. To better understand the interactions between these enzymes, we determined the three-dimensional structure of the murine LMC core. This 0.8-MDa complex is composed of three GLB1 dimers and three CTSA dimers, adopting a triangular architecture maintained through six copies of a unique GLB1-CTSA polar interface. Mutations in this contact surface that occur in GM1-gangliosidosis prevent formation of the LMC in vitro. These findings may facilitate development of therapies for lysosomal storage disorders.

Steric occlusion regulates proximal interactions of acyl carrier protein domain in fungal fatty acid synthase.

The acyl carrier protein (ACP) domain shuttles substrates and reaction intermediates in type I fungal fatty acid synthases via transient protein-protein interactions. Here, using electron cryo-microscopy (cryoEM), we report the structure of a fungal FAS stalled at the dehydration reaction, which precedes the final enoyl reduction in the fatty acid biosynthesis cycle. This conformation revealed multiple contact sites between ACP and the dehydratase (DH) and enoyl reductase (ER) domains that occluded the ACP binding to the adjacent ER domain. Our data suggests a minimal path from the DH to the ER reaction site that requires minute changes in the coordinates of the structured N- and C- termini of the ACP domain.

Label-free assay for thermodynamic analysis of protein-ligand interactions: a multivariate strategy for allosteric ligand screening.

The binding of allosteric ligands to protein can induce changes to the holoprotein conformation, stability, and activity that have an impact on unfolding dynamics. Herein we report a label-free strategy for determining the dissociation constant of protein-ligand interactions over a wide dynamic range (>10(4), Kd from nano- to millimolar) using capillary electrophoresis that overcomes the constraints of an ideal two-state protein unfolding model. Multivariate analysis of thermodynamic parameters associated with holoprotein unfolding and ligand binding is demonstrated for the classification of cyclic nucleotide analogues that function as allosteric modulators of regulatory proteins, such as the exchange protein directly activated by cAMP (EPAC).