Morphological Characterization of Self-Amplifying mRNA Lipid Nanoparticles.

The mRNA technology has emerged as a rapid modality to develop vaccines during pandemic situations with the potential to protect against endemic diseases. The success of mRNA in producing an antigen is dependent on the ability to deliver mRNA to the cells using a vehicle, which typically consists of a lipid nanoparticle (LNP). Self-amplifying mRNA (SAM) is a synthetic mRNA platform that, besides encoding for the antigen of interest, includes the replication machinery for mRNA amplification in the cells. Thus, SAM can generate many antigen encoding mRNA copies and prolong expression of the antigen with lower doses than those required for conventional mRNA. This work describes the morphology of LNPs containing encapsulated SAM (SAM LNPs), with SAM being three to four times larger than conventional mRNA. We show evidence that SAM changes its conformational structure when encapsulated in LNPs, becoming more compact than the free SAM form. A characteristic "bleb" structure is observed in SAM LNPs, which consists of a lipid-rich core and an aqueous RNA-rich core, both surrounded by a DSPC-rich lipid shell. We used SANS and SAXS data to confirm that the prevalent morphology of the LNP consists of two-core compartments where components are heterogeneously distributed between the two cores and the shell. A capped cylinder core-shell model with two interior compartments was built to capture the overall morphology of the LNP. These findings provide evidence that bleb two-compartment structures can be a representative morphology in SAM LNPs and highlight the need for additional studies that elucidate the role of spherical and bleb morphologies, their mechanisms of formation, and the parameters that lead to a particular morphology for a rational design of LNPs for mRNA delivery.

Extended q-Range X-Ray Scattering Reveals High-Resolution Structural Details of Biomacromolecules in Aqueous Solutions.

We communicate a feasibility study for high-resolution structural characterization of biomacromolecules in aqueous solution from X-ray scattering experiments measured over a range of scattering vectors (q) that is approximately two orders of magnitude wider than used previously for such systems. Scattering data with such an extended q-range enables the recovery of the underlying real-space atomic pair distribution function, which facilitates structure determination. We demonstrate the potential of this method for biomacromolecules using several types of cyclodextrins (CD) as model systems. We successfully identified deviations of the tilting angles for the glycosidic units in CDs in aqueous solutions relative to their values in the crystalline forms of these molecules. Such level of structural detail is inaccessible from standard small angle scattering measurements. Our results call for further exploration of ultra-wide-angle X-ray scattering measurements for biomacromolecules.

A round-robin approach provides a detailed assessment of biomolecular small-angle scattering data reproducibility and yields consensus curves for benchmarking.

Through an expansive international effort that involved data collection on 12 small-angle X-ray scattering (SAXS) and four small-angle neutron scattering (SANS) instruments, 171 SAXS and 76 SANS measurements for five proteins (ribonuclease A, lysozyme, xylanase, urate oxidase and xylose isomerase) were acquired. From these data, the solvent-subtracted protein scattering profiles were shown to be reproducible, with the caveat that an additive constant adjustment was required to account for small errors in solvent subtraction. Further, the major features of the obtained consensus SAXS data over the q measurement range 0-1 Å-1 are consistent with theoretical prediction. The inherently lower statistical precision for SANS limited the reliably measured q-range to <0.5 Å-1, but within the limits of experimental uncertainties the major features of the consensus SANS data were also consistent with prediction for all five proteins measured in H2O and in D2O. Thus, a foundation set of consensus SAS profiles has been obtained for benchmarking scattering-profile prediction from atomic coordinates. Additionally, two sets of SAXS data measured at different facilities to q > 2.2 Å-1 showed good mutual agreement, affirming that this region has interpretable features for structural modelling. SAS measurements with inline size-exclusion chromatography (SEC) proved to be generally superior for eliminating sample heterogeneity, but with unavoidable sample dilution during column elution, while batch SAS data collected at higher concentrations and for longer times provided superior statistical precision. Careful merging of data measured using inline SEC and batch modes, or low- and high-concentration data from batch measurements, was successful in eliminating small amounts of aggregate or interparticle interference from the scattering while providing improved statistical precision overall for the benchmarking data set.

Structural and biophysical properties of farnesylated KRas interacting with the chaperone SmgGDS-558.

KRas is a small GTPase and membrane-bound signaling protein. Newly synthesized KRas is post-translationally modified with a membrane-anchoring prenyl group. KRas chaperones are therapeutic targets in cancer due to their participation in trafficking oncogenic KRas to membranes. SmgGDS splice variants are chaperones for small GTPases with basic residues in their hypervariable domain (HVR), including KRas. SmgGDS-607 escorts pre-prenylated small GTPases, while SmgGDS-558 escorts prenylated small GTPases. We provide a structural description of farnesylated and fully processed KRas (KRas-FMe) in complex with SmgGDS-558 and define biophysical properties of this interaction. Surface plasmon resonance measurements on biomimetic model membranes quantified the thermodynamics of the interaction of SmgGDS with KRas, and small-angle x-ray scattering was used to characterize complexes of SmgGDS-558 and KRas-FMe structurally. Structural models were refined using Monte Carlo and molecular dynamics simulations. Our results indicate that SmgGDS-558 interacts with the HVR and the farnesylated C-terminus of KRas-FMe, but not its G-domain. Therefore, SmgGDS-558 interacts differently with prenylated KRas than prenylated RhoA, whose G-domain was found in close contact with SmgGDS-558 in a recent crystal structure. Using immunoprecipitation assays, we show that SmgGDS-558 binds the GTP-bound, GDP-bound, and nucleotide-free forms of farnesylated and fully processed KRas in cells, consistent with SmgGDS-558 not engaging the G-domain of KRas. We found that the dissociation constant, Kd, for KRas-FMe binding to SmgGDS-558 is comparable with that for the KRas complex with PDEδ, a well-characterized KRas chaperone that also does not interact with the KRas G-domain. These results suggest that KRas interacts in similar ways with the two chaperones SmgGDS-558 and PDEδ. Therapeutic targeting of the SmgGDS-558/KRas complex might prove as useful as targeting the PDEδ/KRas complex in KRas-driven cancers.

Conformational heterogeneity of UCAAUC RNA oligonucleotide from molecular dynamics simulations, SAXS, and NMR experiments.

We describe the conformational ensemble of the single-stranded r(UCAAUC) oligonucleotide obtained using extensive molecular dynamics (MD) simulations and Rosetta's FARFAR2 algorithm. The conformations observed in MD consist of A-form-like structures and variations thereof. These structures are not present in the pool generated using FARFAR2. By comparing with available nuclear magnetic resonance (NMR) measurements, we show that the presence of both A-form-like and other extended conformations is necessary to quantitatively explain experimental data. To further validate our results, we measure solution X-ray scattering (SAXS) data on the RNA hexamer and find that simulations result in more compact structures than observed from these experiments. The integration of simulations with NMR via a maximum entropy approach shows that small modifications to the MD ensemble lead to an improved description of the conformational ensemble. Nevertheless, we identify persisting discrepancies in matching experimental SAXS data.

Ionization and structural properties of mRNA lipid nanoparticles influence expression in intramuscular and intravascular administration.

Lipid Nanoparticles (LNPs) are used to deliver siRNA and COVID-19 mRNA vaccines. The main factor known to determine their delivery efficiency is the pKa of the LNP containing an ionizable lipid. Herein, we report a method that can predict the LNP pKa from the structure of the ionizable lipid. We used theoretical, NMR, fluorescent-dye binding, and electrophoretic mobility methods to comprehensively measure protonation of both the ionizable lipid and the formulated LNP. The pKa of the ionizable lipid was 2-3 units higher than the pKa of the LNP primarily due to proton solvation energy differences between the LNP and aqueous medium. We exploited these results to explain a wide range of delivery efficiencies in vitro and in vivo for intramuscular (IM) and intravascular (IV) administration of different ionizable lipids at escalating ionizable lipid-to-mRNA ratios in the LNP. In addition, we determined that more negatively charged LNPs exhibit higher off-target systemic expression of mRNA in the liver following IM administration. This undesirable systemic off-target expression of mRNA-LNP vaccines could be minimized through appropriate design of the ionizable lipid and LNP.

Structural Characterization and Modeling of a Respiratory Syncytial Virus Fusion Glycoprotein Nanoparticle Vaccine in Solution.

The respiratory syncytial virus (RSV) fusion (F) protein/polysorbate 80 (PS80) nanoparticle vaccine is the most clinically advanced vaccine for maternal immunization and protection of newborns against RSV infection. It is composed of a near-full-length RSV F glycoprotein, with an intact membrane domain, formulated into a stable nanoparticle with PS80 detergent. To understand the structural basis for the efficacy of the vaccine, a comprehensive study of its structure and hydrodynamic properties in solution was performed. Small-angle neutron scattering experiments indicate that the nanoparticle contains an average of 350 PS80 molecules, which form a cylindrical micellar core structure and five RSV F trimers that are arranged around the long axis of the PS80 core. All-atom models of full-length RSV F trimers were built from crystal structures of the soluble ectodomain and arranged around the long axis of the PS80 core, allowing for the generation of an ensemble of conformations that agree with small-angle neutron and X-ray scattering data as well as transmission electron microscopy (TEM) images. Furthermore, the hydrodynamic size of the RSV F nanoparticle was found to be modulated by the molar ratio of PS80 to protein, suggesting a mechanism for nanoparticle assembly involving addition of RSV F trimers to and growth along the long axis of the PS80 core. This study provides structural details of antigen presentation and conformation in the RSV F nanoparticle vaccine, helping to explain the induction of broad immunity and observed clinical efficacy. Small-angle scattering methods provide a general strategy to visualize surface glycoproteins from other pathogens and to structurally characterize nanoparticle vaccines.

Chemical shifts-based similarity restraints improve accuracy of RNA structures determined via NMR.

Determination of structure of RNA via NMR is complicated in large part by the lack of a precise parameterization linking the observed chemical shifts to the underlying geometric parameters. In contrast to proteins, where numerous high-resolution crystal structures serve as coordinate templates for this mapping, such models are rarely available for smaller oligonucleotides accessible via NMR, or they exhibit crystal packing and counter-ion binding artifacts that prevent their use for the chemical shifts analysis. On the other hand, NMR-determined structures of RNA often are not solved at the density of restraints required to precisely define the variable degrees of freedom. In this study we sidestep the problems of direct parameterization of the RNA chemical shifts/structure relationship and examine the effects of imposing local fragmental coordinate similarity restraints based on similarities of the experimental secondary ribose 13C/1H chemical shifts instead. The effect of such chemical shift similarity (CSS) restraints on the structural accuracy is assessed via residual dipolar coupling (RDC)-based cross-validation. Improvements in the coordinate accuracy are observed for all of the six RNA constructs considered here as test cases, which argues for routine inclusion of these terms during NMR-based oligonucleotide structure determination. Such accuracy improvements are expected to facilitate derivation of the chemical shift/structure relationships for RNA.

HIV-1 gp120-CD4-Induced Antibody Complex Elicits CD4 Binding Site-Specific Antibody Response in Mice.

Elicitation of broadly neutralizing Ab (bNAb) responses toward the conserved HIV-1 envelope (Env) CD4 binding site (CD4bs) by vaccination is an important goal for vaccine development and yet to be achieved. The outcome of previous immunogenicity studies suggests that the limited accessibility of the CD4bs and the presence of predominant nonneutralizing determinants (nND) on Env may impede the elicitation of bNAbs and their precursors by vaccination. In this study, we designed a panel of novel immunogens that 1) preferentially expose the CD4bs by selective elimination of glycosylation sites flanking the CD4bs, and 2) minimize the nND immune response by engineering fusion proteins consisting of gp120 Core and one or two CD4-induced (CD4i) mAbs for masking nND epitopes, referred to as gp120-CD4i fusion proteins. As expected, the fusion proteins possess improved antigenicity with retained affinity for VRC01-class, CD4bs-directed bNAbs and dampened affinity for nonneutralizing Abs. We immunized C57BL/6 mice with these fusion proteins and found that overall the fusion proteins elicit more focused CD4bs Ab response than prototypical gp120 Core by serological analysis. Consistently, we found that mice immunized with selected gp120-CD4i fusion proteins have higher frequencies of germinal center-activated B cells and CD4bs-directed memory B cells than those inoculated with parental immunogens. We isolated three mAbs from mice immunized with selected gp120-CD4i fusion proteins and found that their footprints on Env are similar to VRC01-class bNAbs. Thus, using gp120-CD4i fusion proteins with selective glycan deletion as immunogens could focus Ab response toward CD4bs epitope.

Structure of the cell-binding component of the Clostridium difficile binary toxin reveals a di-heptamer macromolecular assembly.

Targeting Clostridium difficile infection is challenging because treatment options are limited, and high recurrence rates are common. One reason for this is that hypervirulent C. difficile strains often have a binary toxin termed the C. difficile toxin, in addition to the enterotoxins TsdA and TsdB. The C. difficile toxin has an enzymatic component, termed CDTa, and a pore-forming or delivery subunit termed CDTb. CDTb was characterized here using a combination of single-particle cryoelectron microscopy, X-ray crystallography, NMR, and other biophysical methods. In the absence of CDTa, 2 di-heptamer structures for activated CDTb (1.0 MDa) were solved at atomic resolution, including a symmetric (SymCDTb; 3.14 Å) and an asymmetric form (AsymCDTb; 2.84 Å). Roles played by 2 receptor-binding domains of activated CDTb were of particular interest since the receptor-binding domain 1 lacks sequence homology to any other known toxin, and the receptor-binding domain 2 is completely absent in other well-studied heptameric toxins (i.e., anthrax). For AsymCDTb, a Ca2+ binding site was discovered in the first receptor-binding domain that is important for its stability, and the second receptor-binding domain was found to be critical for host cell toxicity and the di-heptamer fold for both forms of activated CDTb. Together, these studies represent a starting point for developing structure-based drug-design strategies to target the most severe strains of C. difficile.

Accuracy of MD solvent models in RNA structure refinement assessed via liquid-crystal NMR and spin relaxation data.

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA.

Maximizing accuracy of RNA structure in refinement against residual dipolar couplings.

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

Comment on "Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water".

Riback et al (Reports, 13 October 2017, p. 238) used small-angle x-ray scattering (SAXS) experiments to infer a degree of compaction for unfolded proteins in water versus chemical denaturant that is highly consistent with the results from Förster resonance energy transfer (FRET) experiments. There is thus no "contradiction" between the two methods, nor evidence to support their claim that commonly used FRET fluorophores cause protein compaction.

PAGE4 and Conformational Switching: Insights from Molecular Dynamics Simulations and Implications for Prostate Cancer.

Prostate-associated gene 4 (PAGE4) is an intrinsically disordered protein implicated in prostate cancer. Thestress-response kinase homeodomain-interacting protein kinase 1 (HIPK1) phosphorylates two residues in PAGE4, serine 9 and threonine 51. Phosphorylation of these two residues facilitates the interaction of PAGE4 with activator protein-1 (AP-1) transcription factor complex to potentiate AP-1's activity. In contrast, hyperphosphorylation of PAGE4 by CDC-like kinase 2 (CLK2) attenuates this interaction with AP-1. Small-angleX-ray scattering and single-molecule fluorescence resonance energy transfer measurements have shown that PAGE4 expands upon hyperphosphorylation and that this expansion is localized to its N-terminal half. To understand the interactions underlying this structural transition, we performed molecular dynamics simulations using Atomistic AWSEM, a multi-scale molecular model that combines atomistic and coarse-grained simulation approaches. Our simulations show that electrostatic interactions drive transient formation of an N-terminal loop, the destabilization of which accounts for the dramatic change in size upon hyperphosphorylation. Phosphorylation also changes the preference of secondary structure formation of the PAGE4 ensemble, which leads to a transition between states that display different degrees of disorder. Finally, we construct a mechanism-based mathematical model that allows us to capture the interactions ofdifferent phosphoforms of PAGE4 with AP-1 and its downstream target, the androgen receptor (AR)-a key therapeutic target in prostate cancer. Our model predicts intracellular oscillatory dynamics of HIPK1-PAGE4, CLK2-PAGE4, and AR activity, indicating phenotypic heterogeneity in an isogenic cell population. Thus, conformational switching of PAGE4 may potentially affect the efficiency of therapeutically targeting AR activity.

A trapped human PPM1A-phosphopeptide complex reveals structural features critical for regulation of PPM protein phosphatase activity.

Metal-dependent protein phosphatases (PPM) are evolutionarily unrelated to other serine/threonine protein phosphatases and are characterized by their requirement for supplementation with millimolar concentrations of Mg2+ or Mn2+ ions for activity in vitro The crystal structure of human PPM1A (also known as PP2Cα), the first PPM structure determined, displays two tightly bound Mn2+ ions in the active site and a small subdomain, termed the Flap, located adjacent to the active site. Some recent crystal structures of bacterial or plant PPM phosphatases have disclosed two tightly bound metal ions and an additional third metal ion in the active site. Here, the crystal structure of the catalytic domain of human PPM1A, PPM1Acat, complexed with a cyclic phosphopeptide, c(MpSIpYVA), a cyclized variant of the activation loop of p38 MAPK (a physiological substrate of PPM1A), revealed three metal ions in the active site. The PPM1Acat D146E-c(MpSIpYVA) complex confirmed the presence of the anticipated third metal ion in the active site of metazoan PPM phosphatases. Biophysical and computational methods suggested that complex formation results in a slightly more compact solution conformation through reduced conformational flexibility of the Flap subdomain. We also observed that the position of the substrate in the active site allows solvent access to the labile third metal-binding site. Enzyme kinetics of PPM1Acat toward a phosphopeptide substrate supported a random-order, bi-substrate mechanism, with substantial interaction between the bound substrate and the labile metal ion. This work illuminates the structural and thermodynamic basis of an innate mechanism regulating the activity of PPM phosphatases.

Prediction of nearest neighbor effects on backbone torsion angles and NMR scalar coupling constants in disordered proteins.

Using fine-tuned hydrogen bonding criteria, a library of coiled peptide fragments has been generated from a large set of high-resolution protein X-ray structures. This library is shown to be an improved representation of ϕ/ψ torsion angles seen in intrinsically disordered proteins (IDPs). The ϕ/ψ torsion angle distribution of the library, on average, provides good agreement with experimentally observed chemical shifts and 3 JHN-Hα coupling constants for a set of five disordered proteins. Inspection of the coil library confirms that nearest-neighbor effects significantly impact the ϕ/ψ distribution of residues in the coil state. Importantly, 3 JHN-Hα coupling constants derived from the nearest-neighbor modulated backbone ϕ distribution in the coil library show improved agreement to experimental values, thereby providing a better way to predict 3 JHN-Hα coupling constants for IDPs, and for identifying locations that deviate from fully random behavior.

Hybrid Applications of Solution Scattering to Aid Structural Biology.

Biomolecular applications of solution X-ray and neutron scattering (SAXS and SANS, respectively) started in late 1960s - early 1970s but were relatively limited in their ability to provide a detailed structural picture and lagged behind what became the two primary methods of experimental structural biology - X-ray crystallography and NMR. However, improvements in both data analysis and instrumentation led to an explosive growth in the number of studies that used small-angle scattering (SAS) for investigation of macromolecular structure, often in combination with other biophysical techniques. Such hybrid applications are nowadays quickly becoming a norm whenever scattering data are used for two reasons. First, it is generally accepted that SAS data on their own cannot lead to a uniquely defined high-resolution structural model, creating a need for supplementing them with information from complementary techniques. Second, solution scattering data are frequently applied in situations when a method such NMR or X-ray crystallography cannot provide a satisfactory structural picture, which makes these additional restraints highly desirable. Maturation of the hybrid bio-SAS approaches brings to light new questions including completeness of the conformational space sampling, model validation, and data compatibility.

Phosphorylation-induced conformational dynamics in an intrinsically disordered protein and potential role in phenotypic heterogeneity.

Intrinsically disordered proteins (IDPs) that lack a unique 3D structure and comprise a large fraction of the human proteome play important roles in numerous cellular functions. Prostate-Associated Gene 4 (PAGE4) is an IDP that acts as a potentiator of the Activator Protein-1 (AP-1) transcription factor. Homeodomain-Interacting Protein Kinase 1 (HIPK1) phosphorylates PAGE4 at S9 and T51, but only T51 is critical for its activity. Here, we identify a second kinase, CDC-Like Kinase 2 (CLK2), which acts on PAGE4 and hyperphosphorylates it at multiple S/T residues, including S9 and T51. We demonstrate that HIPK1 is expressed in both androgen-dependent and androgen-independent prostate cancer (PCa) cells, whereas CLK2 and PAGE4 are expressed only in androgen-dependent cells. Cell-based studies indicate that PAGE4 interaction with the two kinases leads to opposing functions. HIPK1-phosphorylated PAGE4 (HIPK1-PAGE4) potentiates c-Jun, whereas CLK2-phosphorylated PAGE4 (CLK2-PAGE4) attenuates c-Jun activity. Consistent with the cellular data, biophysical measurements (small-angle X-ray scattering, single-molecule fluorescence resonance energy transfer, and NMR) indicate that HIPK1-PAGE4 exhibits a relatively compact conformational ensemble that binds AP-1, whereas CLK2-PAGE4 is more expanded and resembles a random coil with diminished affinity for AP-1. Taken together, the results suggest that the phosphorylation-induced conformational dynamics of PAGE4 may play a role in modulating changes between PCa cell phenotypes. A mathematical model based on our experimental data demonstrates how differential phosphorylation of PAGE4 can lead to transitions between androgen-dependent and androgen-independent phenotypes by altering the AP-1/androgen receptor regulatory circuit in PCa cells.

Consistent View of Polypeptide Chain Expansion in Chemical Denaturants from Multiple Experimental Methods.

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment.

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.