Folding Landscape Research Articles

G-quadruplex DNA: A Longer Story.

G-quadruplexes (G4s) are distinctive four-stranded DNA or RNA structures found within cells that are thought to play functional roles in gene regulation and transcription, translation, recombination, and DNA damage/repair. While G4 structures can be uni-, bi-, or tetramolecular with respect to strands, folded unimolecular conformations are most significant in vivo. Unimolecular G4 can potentially form in sequences with runs of guanines interspersed with what will become loops in the folded structure: 5'GxLyGxLyGxLyGx, where x is typically 2-4 and y is highly variable. Such sequences are highly conserved and specifically located in genomes. In the folded structure, guanines from each run combine to form planar tetrads with four hydrogen-bonded guanine bases; these tetrads stack on one another to produce four strand segments aligned in specific parallel or antiparallel orientations, connected by the loop sequences. Three types of loops (lateral, diagonal, or "propeller") have been identified. The stacked tetrads form a central cavity that features strong coordination sites for monovalent cations that stabilize the G4 structure, with potassium or sodium preferred. A single monomeric G4 typically forms from a sequence containing roughly 20-30 nucleotides. Such short sequences have been the primary focus of X-ray crystallographic or NMR studies that have produced high-resolution structures of a variety of monomeric G4 conformations. These structures are often used as the basis for drug design efforts to modulate G4 function.We believe that the focus on monomeric G4 structures formed by such short sequences is perhaps myopic. Such short sequences for structural studies are often arbitrarily selected and removed from their native genomic sequence context, and then are often changed from their native sequences by base substitutions or deletions intended to optimize the formation of a homogeneous G4 conformation. We believe instead that G-quadruplexes prefer company and that in a longer natural sequence context multiple adjacent G4 units can form to combine into more complex multimeric G4 structures with richer topographies than simple monomeric forms. Bioinformatic searches of the human genome show that longer sequences with the potential for forming multiple G4 units are common. Telomeric DNA, for example, has a single-stranded overhang of hundreds of nucleotides with the requisite repetitive sequence with the potential for formation of multiple G4s. Numerous extended promoter sequences have similar potentials for multimeric G4 formation. X-ray crystallography and NMR methods are challenged by these longer sequences (>30 nt), so other tools are needed to explore the possible multimeric G4 landscape. We have implemented an integrated structural biology approach to address this challenge. This approach integrates experimental biophysical results with atomic-level molecular modeling and molecular dynamics simulations that provide quantitatively testable model structures. In every long sequence we have studied so far, we found that multimeric G4 structures readily form, with a surprising diversity of structures dependent on the exact native sequence used. In some cases, stable hairpin duplexes form along with G4 units to provide an even richer landscape. This Account provides an overview of our approach and recent progress and provides a new perspective on the G-quadruplex folding landscape.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype–phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

Protein folding problem: enigma, paradox, solution.

The ability of protein chains to spontaneously form their three-dimensional structures is a long-standing mystery in molecular biology. The most conceptual aspect of this mystery is how the protein chain can find its native, "working" spatial structure (which, for not too big protein chains, corresponds to the global free energy minimum) in a biologically reasonable time, without exhaustive enumeration of all possible conformations, which would take billions of years. This is the so-called "Levinthal's paradox." In this review, we discuss the key ideas and discoveries leading to the current understanding of protein folding kinetics, including folding landscapes and funnels, free energy barriers at the folding/unfolding pathways, and the solution of Levinthal's paradox. A special role here is played by the "all-or-none" phase transition occurring at protein folding and unfolding and by the point of thermodynamic (and kinetic) equilibrium between the "native" and the "unfolded" phases of the protein chain (where the theory obtains the simplest form). The modern theory provides an understanding of key features of protein folding and, in good agreement with experiments, it (i) outlines the chain length-dependent range of protein folding times, (ii) predicts the observed maximal size of "foldable" proteins and domains. Besides, it predicts the maximal size of proteins and domains that fold under solely thermodynamic (rather than kinetic) control. Complementarily, a theoretical analysis of the number of possible protein folding patterns, performed at the level of formation and assembly of secondary structures, correctly outlines the upper limit of protein folding times.

Site-specific Heterogeneity of Multi-domain Human Serum Albumin and its Origin: A Red Edge Excitation Shift Study†.

Conformational heterogeneity is a defining characteristic of a protein and is vital in understanding its function and folding landscape. In the present work, we interrogated the presence of conformational heterogeneity in multi-domain human serum albumin in a domain-specific manner using red edge excitation shift (REES) in its native state and also monitored its variation along the unfolding transition. We also looked into the origin of such conformational heterogeneity by varying the solution viscosity. We observed (1) even in the native state, the heterogeneity and dynamics of the side chain exhibit varied behaviors depending on which domain of the multi-domain human serum albumin (HSA) is being examined. (2) When the protein is in the unfolded state, the extent of REES is rendered unimportant since there is a greater quantity of free water present, in addition to the disruption of the protein's structure. (3) While the rigid protein matrix provides the rigidity of domain-I and domain-III, the rigidity of domain-II is provided by water molecules, which indicates that the role of water molecules in providing the rigidity is significant. Overall, our results provide direct evidence of the rigidity and alternate side chain packing arrangement of protein core that varies domain-wise in multi-domain HSA.

Volumetric Properties of the Transition State Ensemble for Protein Folding.

Understanding how high hydrostatic pressure affects biomacromolecular interaction is important for deciphering the molecular mechanisms by which organisms adapt to live at the bottom of the ocean. The relative effect of hydrostatic pressure on the rates of folding/unfolding reactions is defined by the volumetric properties of the transition state ensemble relative to the folded and unfolded states. All-atom structure-based molecular dynamics simulations combined with quantitative computational protocol to compute volumes from three-dimensional coordinates allow volumetric mapping of protein folding landscape. This, is turn, provides qualitative understanding of the effects of hydrostatic pressure on energy landscape of proteins. The computational results for six different proteins are directly benchmark against experimental data and show an excellent agreement. Both experiments and computation show that the transition-state ensemble volume appears to be in-between the folded and unfolded state volumes, and thus the hydrostatic pressure accelerates protein unfolding.

Synergism in the Molecular Crowding of Ligand-Induced Riboswitch Folding: Kinetic/Thermodynamic Insights from Single-Molecule Spectroscopy.

Conformational dynamics in riboswitches involves ligand binding and folding of RNA, each of which can be influenced by excluded volume effects under "crowded" in vivo cellular conditions and thus incompletely characterized by in vitro studies under dilute buffer conditions. In this work, temperature-dependent single-molecule fluorescence resonance energy transfer (FRET) spectroscopy is used to characterize the thermodynamics of (i) cognate ligand and (ii) molecular crowders (PEG, polyethylene glycol) on folding of the B. subtilis LysC lysine riboswitch. With the help of detailed kinetic analysis, we isolate and study the effects of PEG on lysine binding and riboswitch folding steps individually, from which we find that PEG crowding facilitates riboswitch folding primarily via a surprising increase in affinity for the cognate ligand. This is furthermore confirmed by temperature-dependent studies, which reveal that PEG crowding is not purely entropic and instead significantly impacts both enthalpic and entropic contributions to the free energy landscape for folding. The results indicate that PEG molecular crowding/stabilization of the lysine riboswitch is more mechanistically complex and requires extension beyond the conventional picture of purely repulsive solvent-solute steric interactions arising from excluded volume and entropy. Instead, the current experimental FRET data support an alternative multistep mechanism, whereby PEG first entropically crowds the unfolded riboswitch into a "pre-folded" conformation, which in turn greatly increases the ligand binding affinity and thereby enhances the overall equilibrium for riboswitch folding.

Signal Peptide-rheostat Dynamics Delay Secretory Preprotein Folding

Sec secretory proteins are distinguished from cytoplasmic ones by N-terminal signal peptides with multiple roles during post-translational translocation. They contribute to preprotein targeting to the translocase by slowing down folding, binding receptors and triggering secretion. While signal peptides get cleaved after translocation, mature domains traffic further and/or fold into functional states. How signal peptides delay folding temporarily, to keep mature domains translocation-competent, remains unclear. We previously reported that the foldon landscape of the periplasmic prolyl-peptidyl isomerase is altered by its signal peptide and mature domain features. Here, we reveal that the dynamics of signal peptides and mature domains crosstalk. This involves the signal peptide’s hydrophobic helical core, the short unstructured connector to the mature domain and the flexible rheostat at the mature domain N-terminus. Through this cis mechanism the signal peptide delays the formation of early initial foldons thus altering their hierarchy and delaying mature domain folding. We propose that sequence elements outside a protein’s native core exploit their structural dynamics to influence the folding landscape.

Emerging accessibility patterns in long telomeric overhangs

We present single-molecule experimental and computational modeling studies investigating the accessibility of human telomeric overhangs of physiologically relevant lengths. We studied 25 different overhangs that contain 4-28 repeats of GGGTTA (G-Tract) sequence and accommodate one to seven tandem G-quadruplex (GQ) structures. Using the FRET-PAINT method, we probed the distribution of accessible sites via a short imager strand, which is complementary to a G-Tract and transiently binds to available sites. We report accessibility patterns that periodically change with overhang length and interpret these patterns in terms of the underlying folding landscape and folding frustration. Overhangs that have [4n]G-Tracts, (12, 16, 20…) demonstrate the broadest accessibility patterns where the peptide nucleic acid probe accesses G-Tracts throughout the overhang. On the other hand, constructs with [4n+2]G-Tracts, (14, 18, 22…) have narrower patterns where the neighborhood of the junction between single- and double-stranded telomeres is most accessible. We interpret these results as the folding frustration being higher in [4n]G-Tract constructs compared to [4n+2]G-Tract constructs. We also developed a computational model that tests the consistency of different folding stabilities and cooperativities between neighboring GQs with the observed accessibility patterns. Our experimental and computational studies suggest the neighborhood of the junction between single- and double-stranded telomeres is least stable and most accessible, which is significant as this is a potential site where the connection between POT1/TPP1 (bound to single-stranded telomere) and other shelterin proteins (localized on double-stranded telomere) is established.

The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches.

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

Comparing the folding landscapes of evolutionarily divergent procaspase-3.

All caspases evolved from a common ancestor and subsequently developed into two general classes, inflammatory or apoptotic caspases. The caspase-hemoglobinase fold has been conserved throughout nearly one billion years of evolution and is utilized for both the monomeric and dimeric subfamilies of apoptotic caspases, called initiator and effector caspases, respectively. We compared the folding and assembly of procaspase-3b from zebrafish to that of human effector procaspases in order to examine the conservation of the folding landscape. Urea-induced equilibrium folding/unfolding of procaspase-3b showed a minimum three-state folding pathway, where the native dimer isomerizes to a partially folded dimeric intermediate, which then unfolds. A partially folded monomeric intermediate observed in the folding landscape of human procaspase-3 is not well-populated in zebrafish procaspase-3b. By comparing effector caspases from different species, we show that the effector procaspase dimer undergoes a pH-dependent conformational change, and that the conformational species in the folding landscape exhibit similar free energies. Together, the data show that the landscape for the caspase-hemoglobinase fold is conserved, yet it provides flexibility for species-specific stabilization or destabilization of folding intermediates resulting in changes in stability. The common pH-dependent conformational change in the native dimer, which yields an enzymatically inactive species, may provide an additional, albeit reversible, mechanism for controlling caspase activity in the cell.

Using ultraviolet absorption spectroscopy to study nanoswitches based on non-canonical DNA structures

Non-canonical forms of DNA are attracting increasing interest for applications in nanotechnology. It is frequently convenient to characterize DNA molecules using a label-free approach such as ultraviolet absorption spectroscopy. In this paper we present the results of our investigation into the use of this technique to probe the folding of quadruplex and triplex nanoswitches. We confirmed that four G-quartets were necessary for folding at sub-mM concentrations of potassium and found that the wrong choice of sequence for the linker between G-tracts could dramatically disrupt folding, presumably due to the presence of kinetic traps in the folding landscape. In the case of the triplex nanoswitch we examined, we found that the UV spectrum showed a small change in absorbance when a triplex was formed. We anticipate that our results will be of interest to researchers seeking to design DNA nanoswitches based on quadruplexes and triplexes.

Ultrasound-induced protein restructuring and ordered aggregation to form amyloid crystals

Amyloid crystals, a form of ordered protein aggregates documented relatively recently, have not been studied as extensively as amyloid fibres. This study investigates the formation of amyloid crystals with low frequency ultrasound (20 kHz) using β-lactoglobulin, as a model protein for amyloid synthesis. Acoustic cavitation generates localised zones of intense shear, with extreme heat and pressure that could potentially drive the formation of amyloid structures at ambient bulk fluid temperatures (20 ± 1 °C). Thioflavin T fluorescence and electron microscopy showed that low-frequency ultrasound at 20 W/cm3 input power induced β-stacking to produce amyloid crystals in the mesoscopic size range, with a mean length of approximately 22 µm. FTIR spectroscopy indicated a shift towards increased intermolecular antiparallel β-sheet content. An increase in sonication time (0–60 min) and input power (4–24 W/cm3) increased the mean crystal length, but this increase was not linearly proportional to sonication time and input power due to the delayed onset of crystal growth. We propose that acoustic cavitation causes protein unfolding and aggregation and imparts energy to aggregates to cross the torsion barrier, to achieve their lowest energy state as amyloid crystals. The study contributes to a further understanding of protein chemistry relating to the energy landscape of folding and aggregation. Ultrasound presents opportunities for practical applications of amyloid structures, presenting a more adaptable and scalable approach for synthesis.Graphical abstract

Combining High-Pressure NMR and Geometrical Sampling to Obtain a Full Topological Description of Protein Folding Landscapes: Application to the Folding of Two MAX Effectors from Magnaporthe oryzae.

Despite advances in experimental and computational methods, the mechanisms by which an unstructured polypeptide chain regains its unique three-dimensional structure remains one of the main puzzling questions in biology. Single-molecule techniques, ultra-fast perturbation and detection approaches and improvement in all-atom and coarse-grained simulation methods have greatly deepened our understanding of protein folding and the effects of environmental factors on folding landscape. However, a major challenge remains the detailed characterization of the protein folding landscape. Here, we used high hydrostatic pressure 2D NMR spectroscopy to obtain high-resolution experimental structural information in a site-specific manner across the polypeptide sequence and along the folding reaction coordinate. We used this residue-specific information to constrain Cyana3 calculations, in order to obtain a topological description of the entire folding landscape. This approach was used to describe the conformers populating the folding landscape of two small globular proteins, AVR-Pia and AVR-Pib, that belong to the structurally conserved but sequence-unrelated MAX effectors superfamily. Comparing the two folding landscapes, we found that, in spite of their divergent sequences, the folding pathway of these two proteins involves a similar, inescapable, folding intermediate, even if, statistically, the routes used are different.

Ortho-halogen effects: n→π* interactions, halogen bonding, and deciphering chiral attributes in N-aryl glycine peptoid foldamers

Noncovalent interactions play a prominent role in stabilizing the conformation and provide a fertile environment for crafting new structures. Here, we report our ensuing discovery of noncovalent interactions in N-aryl glycine peptoid foldamers emanate through the interaction between the ortho-halogen with the backbone tertiary amide carbonyl. N-aryl peptoids incorporating ortho-substituted halogens can facilitate simultaneous presence of n→π* interactions (Caryl-X···C=O), and halogen bonding (Caryl-X···O=C-) as established by computational modeling, and X-ray crystallography. We observed chiral features including pyramidalization of amide carbonyl carbon as n→π* signature, and stereogenic Caryl-N chiral axis in N-aryl glycine peptoid dimers and a trimer. Inclusion and precise tuning of these noncovalent interactions in the structured peptoid oligomers would further promote the structure and dynamics of a predictable peptoid folding landscape.

Phenol-soluble modulins PSMα3 and PSMβ2 form nanotubes that are cross-α amyloids

Phenol-soluble modulins (PSMs) are peptide-based virulence factors that play significant roles in the pathogenesis of staphylococcal strains in community-associated and hospital-associated infections. In addition to cytotoxicity, PSMs display the propensity to self-assemble into fibrillar species, which may be mediated through the formation of amphipathic conformations. Here, we analyze the self-assembly behavior of two PSMs, PSMα3 and PSMβ2, which are derived from peptides expressed by methicillin-resistant Staphylococcus aureus (MRSA), a significant human pathogen. In both cases, we observed the formation of a mixture of self-assembled species including twisted filaments, helical ribbons, and nanotubes, which can reversibly interconvert in vitro. Cryo–electron microscopy structural analysis of three PSM nanotubes, two derived from PSMα3 and one from PSMβ2, revealed that the assemblies displayed remarkably similar structures based on lateral association of cross-α amyloid protofilaments. The amphipathic helical conformations of PSMα3 and PSMβ2 enforced a bilayer arrangement within the protofilaments that defined the structures of the respective PSMα3 and PSMβ2 nanotubes. We demonstrate that, similar to amyloids based on cross-β protofilaments, cross-α amyloids derived from these PSMs display polymorphism, not only in terms of the global morphology (e.g., twisted filament, helical ribbon, and nanotube) but also with respect to the number of protofilaments within a given peptide assembly. These results suggest that the folding landscape of PSM derivatives may be more complex than originally anticipated and that the assemblies are able to sample a wide range of supramolecular structural space.

PH Effects on the Stability and Folding of Monomeric Caspases

Caspases are an ancient class of cysteinyl proteases that play an integral part in cell development and apoptosis as an evolutionarily conserved function. During apoptosis, the intracellular pH decreases from 7.4 to ~ 6.8 and cytosolic acidification affects the activation of caspase‐3 and pH‐dependent conformational changes affect the rate of auto maturation in the dimeric procaspase‐3, which in turn influences the ability of this enzyme to execute the cell. Proteins are required to fold correctly and attain their native conformation which is vital for protein function and the survival of the organism. The protein folding landscape is sculpted through evolution and there exists an evolutionary selection for functional residues by selecting residues that stabilize the interactions to attain the native conformation while repressing the selection of non‐native interactions that take place during folding. Through the course of evolution, how these proteins fold to the native conformation to contribute to the fitness of the organism while maintaining or altering the energy landscape that is guided by trade‐offs between the stability and the activity of the protein remains elusive. To understand how the energy landscapes changed over time, we resurrected an ancestral caspase (existed about 650 milllion years ago) and we examined the folding pathways of caspases from Cnidarians, specifically corals. We studied the equilibrium unfolding of caspases from O.faveolata (species susceptible to disease) and P.astreoides (species resistant to disease) as a function of pH (between pH 3 and pH 10.5) to examine the changes in the context of stability and folding. The data shows that these monomeric proteins have a high secondary structure at the native state and unfold via a three‐state process closer to the physiological pH (pH 6.5~pH8.5), via an intermediate state which is aggregation prone. However, outside of this pH range, the native protein loses its secondary structure and exists in an intermediate state which is stable. Due to this, these proteins undergo a two‐state unfolding mechanism wherein, the protein transitions from an intermediate state to the unfolded state, reversibly. Overall, the free energy change for these primitive coral caspases is ~12.5kcal/mol, which is higher than what is observed for human caspases. This higher stability might be due to the decrease in the activity and the presence of certain stabilizing residues that are not observed in the human caspases.

Conserved and divergent chaperoning effects of Hsp60/10 chaperonins on protein folding landscapes

The GroEL/ES chaperonin cavity surface charge properties, especially the negative charges, play an important role in its capacity to assist intracavity protein folding. Remarkably, the larger fraction of GroEL/ES negative charges are not conserved among different bacterial species, resulting in a large variation in negative-charge density in the GroEL/ES cavity across prokaryotes. Intriguingly, eukaryotic GroEL/ES homologs have the lowest negative-charge density in the chaperonin cavity. This prompted us to investigate if GroEL’s chaperoning mechanism changed during evolution. Using a model in vivo GroEL/ES substrate, we show that the ability of GroEL/ES to buffer entropic traps in the folding pathway of its substrate was partially dependent upon the negative-charge density inside its cavity. While this activity of GroEL/ES was found to be essential for Escherichia coli, it has been perfected in some organisms and diminished in others. However, irrespective of their charges, all the tested homologs retained their ability to regulate polypeptide chain collapse and remove enthalpic traps from folding pathways. The ability of these GroEL/ES homologs to buffer mutational variations in a model substrate correlated with their negative-charge density. Thus, Hsp60/10 chaperonins in different organisms may have changed to accommodate a different spectrum of mutations on their substrates.

Direct observation of chaperone-modulated talin mechanics with single-molecule resolution

Talin as a critical focal adhesion mechanosensor exhibits force-dependent folding dynamics and concurrent interactions. Being a cytoplasmic protein, talin also might interact with several cytosolic chaperones; however, the roles of chaperones in talin mechanics remain elusive. To address this question, we investigated the force response of a mechanically stable talin domain with a set of well-known unfoldase (DnaJ, DnaK) and foldase (DnaKJE, DsbA) chaperones, using single-molecule magnetic tweezers. Our findings demonstrate that chaperones could affect adhesion proteins’ stability by changing their folding mechanics; while unfoldases reduce their unfolding force from ~11 pN to ~6 pN, foldase shifts it upto ~15 pN. Since talin is mechanically synced within 2 pN force ranges, these changes are significant in cellular conditions. Furthermore, we determined that chaperones directly reshape the energy landscape of talin: unfoldases decrease the unfolding barrier height from 26.8 to 21.7 kBT, while foldases increase it to 33.5 kBT. We reconciled our observations with eukaryotic Hsp70 and Hsp40 and observed their similar function of decreasing the talin unfolding barrier. Quantitative mapping of this chaperone-induced talin folding landscape directly illustrates that chaperones perturb the adhesion protein stability under physiological force, thereby, influencing their force-dependent interactions and adhesion dynamics.

Ligands with polyfluorophenyl moieties promote a local structural rearrangement in the Spinach2 and Broccoli aptamers that increases ligand affinities.

The interaction of nucleic acids with their molecular targets often involves structural reorganization that may traverse a complex folding landscape. With the more recent recognition that many RNAs, both coding and noncoding, may regulate cellular activities by interacting with target molecules, it becomes increasingly important to understand how nucleic acids interact with their targets and how drugs might be developed that can influence critical folding transitions. We have extensively investigated the interaction of the Spinach2 and Broccoli aptamers with a library of small molecule ligands modified by various extensions from the imido nitrogen of DFHBI [(Z)-5-(3,5-difluoro-4-hydroxybenzylidene)-2,3-dimethyl-3,5-dihydro-4H-imidazol-4-one] that reach out from the Spinach2 ligand binding pocket. Studies of the interaction of these compounds with the aptamers revealed that polyfluorophenyl-modified ligands initiate a slow change in aptamer affinity that takes an extended time (half-life of ∼40 min) to achieve. The change in affinity appears to involve an initial disruption of the entrance to the ligand binding pocket followed by a gradual transition to a more defined structure for which the most likely driving force is an interaction of the gateway adenine with a nearby 2'OH group. These results suggest that polyfluorophenyl modifications might increase the ability of small molecule drugs to disrupt local structure and promote RNA remodeling.

The endoplasmic reticulum proteostasis network profoundly shapes the protein sequence space accessible to HIV envelope.

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, at least in theory, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate the levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV-1) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody- and drug-escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s induction. Our data indicate that this phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly, and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is induced, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.