Pressure pushes tRNALys3 into excited conformational states.

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

The microprotein Nrs1 rewires the G1/S transcriptional machinery during nitrogen limitation in budding yeast.

Commitment to cell division at the end of G1 phase, termed Start in the budding yeast Saccharomyces cerevisiae, is strongly influenced by nutrient availability. To identify new dominant activators of Start that might operate under different nutrient conditions, we screened a genome-wide ORF overexpression library for genes that bypass a Start arrest caused by absence of the G1 cyclin Cln3 and the transcriptional activator Bck2. We recovered a hypothetical gene YLR053c, renamed NRS1 for Nitrogen-Responsive Start regulator 1, which encodes a poorly characterized 108 amino acid microprotein. Endogenous Nrs1 was nuclear-localized, restricted to poor nitrogen conditions, induced upon TORC1 inhibition, and cell cycle-regulated with a peak at Start. NRS1 interacted genetically with SWI4 and SWI6, which encode subunits of the main G1/S transcription factor complex SBF. Correspondingly, Nrs1 physically interacted with Swi4 and Swi6 and was localized to G1/S promoter DNA. Nrs1 exhibited inherent transactivation activity, and fusion of Nrs1 to the SBF inhibitor Whi5 was sufficient to suppress other Start defects. Nrs1 appears to be a recently evolved microprotein that rewires the G1/S transcriptional machinery under poor nitrogen conditions.

Hidden intermediates in Mango III RNA aptamer folding revealed by pressure perturbation.

Fluorescent RNA aptamers have the potential to enable routine quantitation and localization of RNA molecules and serve as models for understanding biologically active aptamers. In recent years, several fluorescent aptamers have been selected and modified to improve their properties, revealing that small changes to the RNA or the ligands can modify significantly their fluorescent properties. Although structural biology approaches have revealed the bound, ground state of several fluorescent aptamers, characterization of low-abundance, excited states in these systems is crucial to understanding their folding pathways. Here we use pressure as an alternative variable to probe the suboptimal states of the Mango III aptamer with both fluorescence and NMR spectroscopy approaches. At moderate KCl concentrations, increasing pressure disrupted the G-quadruplex structure of the Mango III RNA and led to an intermediate with lower fluorescence. These observations indicate the existence of suboptimal RNA structural states that still bind the TO1-biotin fluorophore and moderately enhance fluorescence. At higher KCl concentration as well, the intermediate fluorescence state was populated at high pressure, but the G-quadruplex remained stable at high pressure, supporting the notion of parallel folding and/or binding pathways. These results demonstrate the usefulness of pressure for characterizing RNA folding intermediates.

Structural evolution of the ancient enzyme, dissimilatory sulfite reductase.

Dissimilatory sulfite reductase is an ancient enzyme that has linked the global sulfur and carbon biogeochemical cycles since at least 3.47 Gya. While much has been learned about the phylogenetic distribution and diversity of DsrAB across environmental gradients, far less is known about the structural changes that occurred to maintain DsrAB function as the enzyme accompanied diversification of sulfate/sulfite reducing organisms (SRO) into new environments. Analyses of available crystal structures of DsrAB from Archaeoglobus fulgidus and Desulfovibrio vulgaris, representing early and late evolving lineages, respectively, show that certain features of DsrAB are structurally conserved, including active siro-heme binding motifs. Whether such structural features are conserved among DsrAB recovered from varied environments, including hot spring environments that host representatives of the earliest evolving SRO lineage (e.g., MV2-Eury), is not known. To begin to overcome these gaps in our understanding of the evolution of DsrAB, structural models from MV2.Eury were generated and evolutionary sequence co-variance analyses were conducted on a curated DsrAB database. Phylogenetically diverse DsrAB harbor many conserved functional residues including those that ligate active siro-heme(s). However, evolutionary co-variance analysis of monomeric DsrAB subunits revealed several False Positive Evolutionary Couplings (FPEC) that correspond to residues that have co-evolved despite being too spatially distant in the monomeric structure to allow for direct contact. One set of FPECs corresponds to residues that form a structural path between the two active siro-heme moieties across the interface between heterodimers, suggesting the potential for allostery or electron transfer within the enzyme complex. Other FPECs correspond to structural loops and gaps that may have been selected to stabilize enzyme function in different environments. These structural bioinformatics results suggest that DsrAB has maintained allosteric communication pathways between subunits as SRO diversified into new environments. The observations outlined here provide a framework for future biochemical and structural analyses of DsrAB to examine potential allosteric control of this enzyme.

Protein unfolded states populated at high and ambient pressure are similarly compact.

The relationship between the dimensions of pressure-unfolded states of proteins compared with those at ambient pressure is controversial; resolving this issue is related directly to the mechanisms of pressure denaturation. Moreover, a significant pressure dependence of the compactness of unfolded states would complicate the interpretation of folding parameters from pressure perturbation and make comparison to those obtained using alternative perturbation approaches difficult. Here, we determined the compactness of the pressure-unfolded state of a small, cooperatively folding model protein, CTL9-I98A, as a function of temperature. This protein undergoes both thermal unfolding and cold denaturation, and the temperature dependence of the compactness at atmospheric pressure is known. High-pressure small angle x-ray scattering studies, yielding the radius of gyration and high-pressure diffusion ordered spectroscopy NMR experiments, yielding the hydrodynamic radius were carried out as a function of temperature at 250 MPa, a pressure at which the protein is unfolded. The radius of gyration values obtained at any given temperature at 250 MPa were similar to those reported previously at ambient pressure, and the trends with temperature are similar as well, although the pressure-unfolded state appears to undergo more pronounced expansion at high temperature than the unfolded state at atmospheric pressure. At 250 MPa, the compaction of the unfolded chain was maximal between 25 and 30°C, and the chain expanded upon both cooling and heating. These results reveal that the pressure-unfolded state of this protein is very similar to that observed at ambient pressure, demonstrating that pressure perturbation represents a powerful approach for observing the unfolded states of proteins under otherwise near-native conditions.

The consequences of cavity creation on the folding landscape of a repeat protein depend upon context.

The effect of introducing internal cavities on protein native structure and global stability has been well documented, but the consequences of these packing defects on folding free-energy landscapes have received less attention. We investigated the effects of cavity creation on the folding landscape of the leucine-rich repeat protein pp32 by high-pressure (HP) and urea-dependent NMR and high-pressure small-angle X-ray scattering (HPSAXS). Despite a modest global energetic perturbation, cavity creation in the N-terminal capping motif (N-cap) resulted in very strong deviation from two-state unfolding behavior. In contrast, introduction of a cavity in the most stable, C-terminal half of pp32 led to highly concerted unfolding, presumably because the decrease in stability by the mutations attenuated the N- to C-terminal stability gradient present in WT pp32. Interestingly, enlarging the central cavity of the protein led to the population under pressure of a distinct intermediate in which the N-cap and repeats 1-4 were nearly completely unfolded, while the fifth repeat and the C-terminal capping motif remained fully folded. Thus, despite modest effects on global stability, introducing internal cavities can have starkly distinct repercussions on the conformational landscape of a protein, depending on their structural and energetic context.

Quantitative High-Resolution Imaging of Live Microbial Cells at High Hydrostatic Pressure.

The majority of the Earth's microbial biomass exists in the deep biosphere, in the deep ocean, and within the Earth's crust. Although other physical parameters in these environments, such as temperature or pH, can differ substantially, they are all under high pressures. Beyond emerging genomic information, little is known about the molecular mechanisms underlying the ability of these organisms to survive and grow at pressures that can reach over 1000-fold the pressure on the Earth's surface. The mechanisms of pressure adaptation are also important in food safety, with the increasing use of high-pressure food processing. Advanced imaging represents an important tool for exploring microbial adaptation and response to environmental changes. Here, we describe implementation of a high-pressure sample chamber with a two-photon scanning microscope system, allowing for the first time, to our knowledge, quantitative high-resolution two-photon imaging at 100 MPa of living microbes from all three kingdoms of life. We adapted this setup for fluorescence lifetime imaging microscopy with phasor analysis (FLIM/Phasor) and investigated metabolic responses to pressure of live cells from mesophilic yeast and bacterial strains, as well as the piezophilic archaeon Archaeoglobus fulgidus. We also monitored by fluorescence intensity fluctuation-based methods (scanning number and brightness and raster scanning imaging correlation spectroscopy) the effect of pressure on the chromosome-associated protein HU and on the ParB partition protein in Escherichia coli, revealing partially reversible dissociation of ParB foci and concomitant nucleoid condensation. These results provide a proof of principle that quantitative, high-resolution imaging of live microbial cells can be carried out at pressures equivalent to those in the deepest ocean trenches.

Pressure-Temperature Analysis of the Stability of the CTL9 Domain Reveals Hidden Intermediates.

The observation of two-state unfolding for many small single-domain proteins by denaturants has led to speculation that protein sequences may have evolved to limit the population of partially folded states that could be detrimental to fitness. How such strong cooperativity arises from a multitude of individual interactions is not well understood. Here, we investigate the stability and folding cooperativity of the C-terminal domain of the ribosomal protein L9 in the pressure-temperature plane using site-specific NMR. In contrast to apparent cooperative unfolding detected with denaturant-induced and thermal-induced unfolding experiments and stopped-flow refolding studies at ambient pressure, NMR-detected pressure unfolding revealed significant deviation from two-state behavior, with a core region that was selectively destabilized by increasing temperature. Comparison of pressure-dependent NMR signals from both the folded and unfolded states revealed the population of at least one invisible excited state at atmospheric pressure. The core destabilizing cavity-creating I98A mutation apparently increased the cooperativity of the loss of folded-state peak intensity while also increasing the population of this invisible excited state present at atmospheric pressure. These observations highlight how local stability is subtly modulated by sequence to tune protein conformational landscapes and illustrate the ability of pressure- and temperature-dependent studies to reveal otherwise hidden states.

Structural mapping of fluorescently-tagged, functional nhTMEM16 scramblase in a lipid bilayer.

Most members of the TransMEMbrane protein 16 (TMEM16) family are Ca2+-regulated scramblases that facilitate the bidirectional movement of phospholipids across membranes necessary for diverse physiological processes. The nhTMEM16 scramblase (from the fungus Nectria hematococca) is a homodimer with a large cytoplasmic region and a hydrophilic, membrane-exposed groove in each monomer. The groove provides the transbilayer conduit for lipids, but the mechanism by which Ca2+ regulates it is not clear. Because fusion of large protein tags at either the N or C terminus abolishes nhTMEM16 activity, we hypothesized that its cytoplasmic portion containing both termini may regulate lipid translocation via a Ca2+-dependent conformational change. To test this hypothesis, here we used fluorescence methods to map key distances within the nhTMEM16 homodimer and between its termini and the membrane. To this end, we developed functional nhTMEM16 variants bearing an acyl carrier protein (ACP) tag at one or both of the termini. These constructs were fluorescently labeled by ACP synthase-mediated insertion of CoA-conjugated fluorophores and reconstituted into vesicles containing fluorescent lipids to obtain the distance of closest approach between the labeled tag and the membrane via FRET. Fluorescence lifetime measurements with phasor analysis were used to determine the distance between the N and C termini of partnering monomers in the nhTMEM16 homodimer. We now report that the measured distances do not vary significantly between Ca2+-replete and EGTA-treated samples, indicating that whereas the cytoplasmic portion of the protein is important for function, it does not appear to regulate scramblase activity via a detectable conformational change.

High pressure activation of the Mrr restriction endonuclease in Escherichia coli involves tetramer dissociation.

A sub-lethal hydrostatic pressure (HP) shock of ∼100 MPa elicits a RecA-dependent DNA damage (SOS) response in Escherichia coli K-12, despite the fact that pressure cannot compromise the covalent integrity of DNA. Prior screens for HP resistance identified Mrr (Methylated adenine Recognition and Restriction), a Type IV restriction endonuclease (REase), as instigator for this enigmatic HP-induced SOS response. Type IV REases tend to target modified DNA sites, and E. coli Mrr activity was previously shown to be elicited by expression of the foreign M.HhaII Type II methytransferase (MTase), as well. Here we measured the concentration and stoichiometry of a functional GFP-Mrr fusion protein using in vivo fluorescence fluctuation microscopy. Our results demonstrate that Mrr is a tetramer in unstressed cells, but shifts to a dimer after HP shock or co-expression with M.HhaII. Based on the differences in reversibility of tetramer dissociation observed for wild-type GFP-Mrr and a catalytic mutant upon HP shock compared to M.HhaII expression, we propose a model by which (i) HP triggers Mrr activity by directly pushing inactive Mrr tetramers to dissociate into active Mrr dimers, while (ii) M.HhaII triggers Mrr activity by creating high affinity target sites on the chromosome, which pull the equilibrium from inactive tetrameric Mrr toward active dimer.

The Deep Sea Osmolyte Trimethylamine N-Oxide and Macromolecular Crowders Rescue the Antiparallel Conformation of the Human Telomeric G-Quadruplex from Urea and Pressure Stress.

Organisms are thriving in the deep sea at pressures up to the 1 kbar level, which imposes severe stress on the conformational dynamics and stability of their biomolecules. The impact of osmolytes and macromolecular crowders, mimicking intracellular conditions, on the effect of pressure on the conformational dynamics of a human telomeric G-quadruplex (G4) DNA is explored in this study employing single-molecule Förster resonance energy transfer (FRET) experiments. In neat buffer, pressurization favors the parallel/hybrid state of the G4-DNA over the antiparallel conformation at ≈400 bar, finally leading to unfolding beyond 1000 bar. High-pressure NMR data support these findings. The folded topological conformers have different solvent accessible surface areas and cavity volumes, leading to different volumetric properties and hence pressure stabilities. The deep-sea osmolyte trimethylamine N-oxide (TMAO) and macromolecular crowding agents are able to effectively rescue the G4-DNA from unfolding in the whole pressure range encountered on Earth.

Putting the Piezolyte Hypothesis under Pressure.

A group of small molecules that stabilize proteins against high hydrostatic pressure has been classified as piezolytes, a subset of stabilizing cosolutes. This distinction would imply that piezolytes counteract the effects of high hydrostatic pressure through effects on the volumetric properties of the protein. The purpose of this study was to determine if cosolutes proposed to be piezolytes have an effect on the volumetric properties of proteins through direct experimental measurements of volume changes upon unfolding of model proteins lysozyme and ribonuclease A, in solutions containing varying cosolute concentrations. Solutions containing the proposed piezolytes glutamate, sarcosine, and betaine were used, as well as solutions containing the denaturants guanidinium hydrochloride and urea. Changes in thermostability were monitored using differential scanning calorimetry whereas changes in volume were monitored using pressure perturbation calorimetry. Our findings indicate that increasing stabilizing cosolute concentration increases the stability and transition temperature of the protein, but does not change the temperature dependence of volume changes upon unfolding. The results suggest that the pressure stability of a protein in solution is not directly affected by the presence of these proposed piezolytes, and so they cannot be granted this distinction.

V67L Mutation Fills an Internal Cavity To Stabilize RecA Mtu Intein.

Inteins mediate protein splicing, which has found extensive applications in protein science and biotechnology. In the Mycobacterium tuberculosis RecA mini-mini intein (ΔΔIhh), a single valine to leucine substitution at position 67 (V67L) dramatically increases intein stability and activity. However, crystal structures show that the V67L mutation causes minimal structural rearrangements, with a root-mean-square deviation of 0.2 Å between ΔΔIhh-V67 and ΔΔIhh-L67. Thus, the structural mechanisms for V67L stabilization and activation remain poorly understood. In this study, we used intrinsic tryptophan fluorescence, high-pressure nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations to probe the structural basis of V67L stabilization of the intein fold. Guanidine hydrochloride denaturation monitored by fluorescence yielded free energy changes (ΔGf°) of -4.4 and -6.9 kcal mol-1 for ΔΔIhh-V67 and ΔΔIhh-L67, respectively. High-pressure NMR showed that ΔΔIhh-L67 is more resistant to pressure-induced unfolding than ΔΔIhh-V67 is. The change in the volume of folding (ΔVf) was significantly larger for V67 (71 ± 2 mL mol-1) than for L67 (58 ± 3 mL mol-1) inteins. The measured difference in ΔVf (13 ± 3 mL mol-1) roughly corresponds to the volume of the additional methylene group for Leu, supporting the notion that the V67L mutation fills a nearby cavity to enhance intein stability. In addition, we performed MD simulations to show that V67L decreases side chain dynamics and conformational entropy at the active site. It is plausible that changes in cavities in V67L can also mediate allosteric effects to change active site dynamics and enhance intein activity.

High-Resolution Mapping of a Repeat Protein Folding Free Energy Landscape.

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the conformational ensemble along the entire folding reaction coordinate. Simulations can provide this level of insight for small proteins. In contrast, with the exception of hydrogen exchange, which does not monitor folding directly, experimental studies of protein folding have not yielded such structural and energetic detail. NMR can provide residue specific atomic level structural information, but its implementation in protein folding studies using chemical or temperature perturbation is problematic. Here we present a highly detailed structural and energetic map of the entire folding landscape of the leucine-rich repeat protein, pp32 (Anp32), obtained by combining pressure-dependent site-specific 1H-15N HSQC data with coarse-grained molecular dynamics simulations. The results obtained using this equilibrium approach demonstrate that the main barrier to folding of pp32 is quite broad and lies near the unfolded state, with structure apparent only in the C-terminal region. Significant deviation from two-state unfolding under pressure reveals an intermediate on the folded side of the main barrier in which the N-terminal region is disordered. A nonlinear temperature dependence of the population of this intermediate suggests a large heat capacity change associated with its formation. The combination of pressure, which favors the population of folding intermediates relative to chemical denaturants; NMR, which allows their observation; and constrained structure-based simulations yield unparalleled insight into protein folding mechanisms.

High Pressure ZZ-Exchange NMR Reveals Key Features of Protein Folding Transition States.

Understanding protein folding mechanisms and their sequence dependence requires the determination of residue-specific apparent kinetic rate constants for the folding and unfolding reactions. Conventional two-dimensional NMR, such as HSQC experiments, can provide residue-specific information for proteins. However, folding is generally too fast for such experiments. ZZ-exchange NMR spectroscopy allows determination of folding and unfolding rates on much faster time scales, yet even this regime is not fast enough for many protein folding reactions. The application of high hydrostatic pressure slows folding by orders of magnitude due to positive activation volumes for the folding reaction. We combined high pressure perturbation with ZZ-exchange spectroscopy on two autonomously folding protein domains derived from the ribosomal protein, L9. We obtained residue-specific apparent rates at 2500 bar for the N-terminal domain of L9 (NTL9), and rates at atmospheric pressure for a mutant of the C-terminal domain (CTL9) from pressure dependent ZZ-exchange measurements. Our results revealed that NTL9 folding is almost perfectly two-state, while small deviations from two-state behavior were observed for CTL9. Both domains exhibited large positive activation volumes for folding. The volumetric properties of these domains reveal that their transition states contain most of the internal solvent excluded voids that are found in the hydrophobic cores of the respective native states. These results demonstrate that by coupling it with high pressure, ZZ-exchange can be extended to investigate a large number of protein conformational transitions.

The stoichiometry of scaffold complexes in living neurons - DLC2 functions as a dimerization engine for GKAP.

Quantitative spatio-temporal characterization of protein interactions in living cells remains a major challenge facing modern biology. We have investigated in living neurons the spatial dependence of the stoichiometry of interactions between two core proteins of the N-methyl-D-aspartate (NMDA)-receptor-associated scaffolding complex, GKAP (also known as DLGAP1) and DLC2 (also known as DYNLL2), using a novel variation of fluorescence fluctuation microscopy called two-photon scanning number and brightness (sN&B). We found that dimerization of DLC2 was required for its interaction with GKAP, which, in turn, potentiated GKAP self-association. In the dendritic shaft, the DLC2-GKAP hetero-oligomeric complexes were composed mainly of two DLC2 and two GKAP monomers, whereas, in spines, the hetero-complexes were much larger, with an average of ∼16 DLC2 and ∼13 GKAP monomers. Disruption of the GKAP-DLC2 interaction strongly destabilized the oligomers, decreasing the spine-preferential localization of GKAP and inhibiting NMDA receptor activity. Hence, DLC2 serves a hub function in the control of glutamatergic transmission by ordering GKAP-containing complexes in dendritic spines. Beyond illuminating the role of DLC2-GKAP interactions in glutamatergic signaling, these data underscore the power of the sN&B approach for quantitative spatio-temporal imaging of other important protein complexes.

Recruitment, assembly, and molecular architecture of the SpoIIIE DNA pump revealed by superresolution microscopy.

ATP-fuelled molecular motors are responsible for rapid and specific transfer of double-stranded DNA during several fundamental processes, such as cell division, sporulation, bacterial conjugation, and viral DNA transport. A dramatic example of intercompartmental DNA transfer occurs during sporulation in Bacillus subtilis, in which two-thirds of a chromosome is transported across a division septum by the SpoIIIE ATPase. Here, we use photo-activated localization microscopy, structured illumination microscopy, and fluorescence fluctuation microscopy to investigate the mechanism of recruitment and assembly of the SpoIIIE pump and the molecular architecture of the DNA translocation complex. We find that SpoIIIE assembles into ∼45 nm complexes that are recruited to nascent sites of septation, and are subsequently escorted by the constriction machinery to the center of sporulation and division septa. SpoIIIE complexes contain 47±20 SpoIIIE molecules, a majority of which are assembled into hexamers. Finally, we show that directional DNA translocation leads to the establishment of a compartment-specific, asymmetric complex that exports DNA. Our data are inconsistent with the notion that SpoIIIE forms paired DNA conducting channels across fused membranes. Rather, our results support a model in which DNA translocation occurs through an aqueous DNA-conducting pore that could be structurally maintained by the divisional machinery, with SpoIIIE acting as a checkpoint preventing membrane fusion until completion of chromosome segregation. Our findings and proposed mechanism, and our unique combination of innovating methodologies, are relevant to the understanding of bacterial cell division, and may illuminate the mechanisms of other complex machineries involved in DNA conjugation and protein transport across membranes.

Why and How Does Pressure Unfold Proteins?

This year, 2014, marks the 100th anniversary of the first publication reporting the denaturation of proteins by high hydrostatic pressure (Bridgman 1914). Since that time a large and recently increasing number of studies of pressure effects on protein stability have been published, yet the mechanism for the action of pressure on proteins remains subject to considerable debate. This review will present an overview from this author's perspective of where this debate stands today.

Exploring the Protein Folding Pathway with High-Pressure NMR: Steady-State and Kinetics Studies.

Defining the physical-chemical determinants of protein folding and stability, under normal and pathological conditions has constituted a major subfield in biophysical chemistry for over 50 years. Although a great deal of progress has been made in recent years towards this goal, a number of important questions remain. These include characterizing the structural, thermodynamic and dynamic properties of the barriers between conformational states on the protein energy landscape, understanding the sequence dependence of folding cooperativity, defining more clearly the role of solvation in controlling protein stability and dynamics and probing the high energy thermodynamic states in the native state basin and their role in misfolding and aggregation. Fundamental to the elucidation of these questions is a complete thermodynamic parameterization of protein folding determinants. In this chapter, we describe the use of high-pressure coupled to Nuclear Magnetic Resonance (NMR) spectroscopy to reveal unprecedented details on the folding energy landscape of proteins.

Evolutionarily Conserved Pattern of Interactions in a Protein Revealed by Local Thermal Expansion Properties.

The way in which the network of intramolecular interactions determines the cooperative folding and conformational dynamics of a protein remains poorly understood. High-pressure NMR spectroscopy is uniquely suited to examine this problem because it combines the site-specific resolution of the NMR experiments with the local character of pressure perturbations. Here we report on the temperature dependence of the site-specific volumetric properties of various forms of staphylococcal nuclease (SNase), including three variants with engineered internal cavities, as measured with high-pressure NMR spectroscopy. The strong temperature dependence of pressure-induced unfolding arises from poorly understood differences in thermal expansion between the folded and unfolded states. A significant inverse correlation was observed between the global thermal expansion of the folded proteins and the number of strong intramolecular hydrogen bonds, as determined by the temperature coefficient of the backbone amide chemical shifts. Comparison of the identity of these strong H-bonds with the co-evolution of pairs of residues in the SNase protein family suggests that the architecture of the interactions detected in the NMR experiments could be linked to a functional aspect of the protein. Moreover, the temperature dependence of the residue-specific volume changes of unfolding yielded residue-specific differences in expansivity and revealed how mutations impact intramolecular interaction patterns. These results show that intramolecular interactions in the folded states of proteins impose constraints against thermal expansion and that, hence, knowledge of site-specific thermal expansivity offers insight into the patterns of strong intramolecular interactions and other local determinants of protein stability, cooperativity, and potentially also of function.