Size and sequence and the volume change of protein folding

Analysis of the effect of cations on protein conformational stability using solid-state nanopores

The conformation of proteins is closely related to their biological functions, and it is affected by many factors, including the type of cations in solution. However, it is difficult to detect the conformational changes of a protein in situ. As a single-molecule sensing technology, nanopores can convert molecular structural information into analyzable current signals within a reasonable time range. Herein, we detect and analyze the effects of two different types of monovalent cations (Na+ and Li+) on a model protein bovine serum albumin (BSA) conformation using SiNx nanopores with different diameters. The quantitative analysis results show that the excluded volume of BSA in LiCl salt solutions is larger than the value in NaCl solution, indicating that Li+ is more prone to unfolding the proteins and making them unstable. This study demonstrated that nanopores enable the in situ detection of the structure of proteins at the single-molecule level and provide a new approach for the quantitative analysis of proteins.

Does a Similar 3D Structure Mean a Similar Folding Pathway? The Presence of a C-Terminal α-Helical Extension in the 3D Structure of MAX60 Drastically Changes the Folding Pathway Described for Other MAX-Effectors from Magnaporthe oryzae

Does a similar 3D structure mean a similar folding pathway? This question is particularly meaningful when it concerns proteins sharing a similar 3D structure, but low sequence identity or homology. MAX effectors secreted by the phytopathogenic fungus Magnaporthe oryzae present such characteristics. They share a common 3D structure, a ß-sandwich with the same topology for all the family members, but an extremely low sequence identity/homology. In a previous study, we have investigated the folding of two MAX effectors, AVR-Pia and AVR-Pib, using High-Hydrostatic-Pressure NMR and found that they display a similar folding pathway, with a common folding intermediate. In the present work, we used a similar strategy to investigate the folding conformational landscape of another MAX effector, MAX60, and found a very different folding intermediate. Our analysis strongly supports that the presence of a C-terminal α-helical extension in the 3D structure of MAX60 could be responsible for its different folding pathway.

Active Pd Catalyst for the Selective Synthesis of Methylated Amines with Methanol

Selective N-methylation of amines with methanol is an important reaction in the synthesis of high-value-added fine chemicals, including dyes, surfactants, pharmaceuticals, agrochemicals, and materials. However, N-methylated amines possess higher reactivities and are prone to further transform into N,N-dimethylated amines. Therefore, it is still a challenge to controllably regulate the selectivity of N-methylation using heterogeneous catalysts without the use of base. Herein, we developed a series of Pd/Zn(Al)O catalysts with abundant basic sites, and the selectivity of N-methylation was controlled by a heterogeneous Pd/Zn(Al)O catalyst with a Zn/Al ratio of 10 and a Pd loading of 0.4 wt % in the pressure of H₂. The experimental results showed that the appropriate basic properties of the catalyst were beneficial to form the desired N-methylated amine. The low loading of Pd in the catalyst was highly dispersed on the support, providing sufficient active sites. These were attributed to the Zn vacancies formed by Al-doped Zn, which were beneficial to form the highly active and stable Pd sites. Furthermore, a series of amines and nitrobenzenes with different functional groups were well tolerated for the selective synthesis of N-methylated amines in the absence of base.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Heat and cold denaturation of yeast frataxin: The effect of pressure

Yfh1 is a yeast protein with the peculiar characteristic to undergo, in the absence of salt, cold denaturation at temperatures above the water freezing point. This feature makes the protein particularly interesting for studies aiming at understanding the rules that determine protein fold stability. Here, we present the phase diagram of Yfh1 unfolding as a function of pressure (0.1-500 MPa) and temperature 278-313 K (5-40°C) both in the absence and in the presence of stabilizers using Trp fluorescence as a monitor. The protein showed a remarkable sensitivity to pressure: at 293 K, pressures around 10 MPa are sufficient to cause 50% of unfolding. Higher pressures were required for the unfolding of the protein in the presence of stabilizers. The phase diagram on the pressure-temperature plane together with a critical comparison between our results and those found in the literature allowed us to draw conclusions on the mechanism of the unfolding process under different environmental conditions.

Molecular Dynamics Calculations of Partial Molar Volumes of Amino Acids in Aqueous Solutions

Partial molar volumes of amino acids in their zwitterionic and molecular forms have been calculated using molecular dynamics simulations of these systems in aqueous solutions. Calculations performed with the TIP4P, SPC (rigid and flexible), SPC/E, and polarizable water models show that the choice of water model can have a significant impact on the calculated volumes. The effect of treatment of long-range electrostatic interactions on the calculated results was also investigated. Volumes obtained in simulations with a properly chosen water model fit well the experimental data for both molecular and zwitterionic forms of amino acids.

Pressure Unfolding of Proteins: New Insights into the Role of Bound Water

High pressures can be detrimental for protein stability, resulting in unfolding and loss of function. This phenomenon occurs because the unfolding transition is accompanied by a decrease in volume, which is typically attributed to the elimination of cavities that are present within the native state as a result of packing defects. We present a novel computational approach that enables the study of pressure unfolding in atomistically detailed protein models in implicit solvent. We include the effect of pressure using a transfer free energy term that allows us to decouple the effect of protein residues and bound water molecules on the volume change upon unfolding. We discuss molecular dynamics simulations results using this protocol for two model proteins, Trp-cage and staphylococcal nuclease (SNase). We find that the volume reduction of bound water is the key energetic term that drives protein denaturation under the effect of pressure, for both Trp-cage and SNase. However, we note differences in unfolding mechanisms between the smaller Trp-cage and the larger SNase protein. Indeed, the unfolding of SNase, but not Trp-cage, is seen to be further accompanied by a reduction in the volume of internal cavities. Our results indicate that, for small peptides, like Trp-cage, pressure denaturation is driven by the increase in solvent accessibility upon unfolding, and the subsequent increase in the number of bound water molecules. For larger proteins, like SNase, the cavities within the native fold act as weak spots, determining the overall resistance to pressure denaturation. Our simulations display a striking agreement with the pressure-unfolding profile experimentally obtained for SNase and represent a promising approach for a computationally efficient and accurate exploration of pressure-induced denaturation of proteins.

The Right-Handed Parallel β-Helix Topology of Erwinia chrysanthemi Pectin Methylesterase Is Intimately Associated with Both Sequential Folding and Resistance to High Pressure

The complex topologies of large multi-domain globular proteins make the study of their folding and assembly particularly demanding. It is often characterized by complex kinetics and undesired side reactions, such as aggregation. The structural simplicity of tandem-repeat proteins, which are characterized by the repetition of a basic structural motif and are stabilized exclusively by sequentially localized contacts, has provided opportunities for dissecting their folding landscapes. In this study, we focus on the Erwinia chrysanthemi pectin methylesterase (342 residues), an all-β pectinolytic enzyme with a right-handed parallel β-helix structure. Chemicals and pressure were chosen as denaturants and a variety of optical techniques were used in conjunction with stopped-flow equipment to investigate the folding mechanism of the enzyme at 25 °C. Under equilibrium conditions, both chemical- and pressure-induced unfolding show two-state transitions, with average conformational stability (ΔG° = 35 ± 5 kJ·mol⁻¹) but exceptionally high resistance to pressure (P_m = 800 ± 7 MPa). Stopped-flow kinetic experiments revealed a very rapid (τ < 1 ms) hydrophobic collapse accompanied by the formation of an extended secondary structure but did not reveal stable tertiary contacts. This is followed by three distinct cooperative phases and the significant population of two intermediate species. The kinetics followed by intrinsic fluorescence shows a lag phase, strongly indicating that these intermediates are productive species on a sequential folding pathway, for which we propose a plausible model. These combined data demonstrate that even a large repeat protein can fold in a highly cooperative manner.

Comparative Assessment of NMR Probes for the Experimental Description of Protein Folding Pathways with High-Pressure NMR

Multidimensional NMR intrinsically provides multiple probes that can be used for deciphering the folding pathways of proteins: NH amide and CαHα groups are strategically located on the backbone of the protein, while CH₃ groups, on the side-chain of methylated residues, are involved in important stabilizing interactions in the hydrophobic core. Combined with high hydrostatic pressure, these observables provide a powerful tool to explore the conformational landscapes of proteins. In the present study, we made a comparative assessment of the NH, CαHα, and CH₃ groups for analyzing the unfolding pathway of ∆+PHS Staphylococcal Nuclease. These probes yield a similar description of the folding pathway, with virtually identical thermodynamic parameters for the unfolding reaction, despite some notable differences. Thus, if partial unfolding begins at identical pressure for these observables (especially in the case of backbone probes) and concerns similar regions of the molecule, the residues involved in contact losses are not necessarily the same. In addition, an unexpected slight shift toward higher pressure was observed in the sequence of the scenario of unfolding with CαHα when compared to amide groups.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

Pressure and Chemical Unfolding of an α-Helical Bundle Protein: The GH2 Domain of the Protein Adaptor GIPC1

When combined with NMR spectroscopy, high hydrostatic pressure is an alternative perturbation method used to destabilize globular proteins that has proven to be particularly well suited for exploring the unfolding energy landscape of small single-domain proteins. To date, investigations of the unfolding landscape of all-β or mixed-α/β protein scaffolds are well documented, whereas such data are lacking for all-α protein domains. Here we report the NMR study of the unfolding pathways of GIPC1-GH2, a small α-helical bundle domain made of four antiparallel α-helices. High-pressure perturbation was combined with NMR spectroscopy to unravel the unfolding landscape at three different temperatures. The results were compared to those obtained from classical chemical denaturation. Whatever the perturbation used, the loss of secondary and tertiary contacts within the protein scaffold is almost simultaneous. The unfolding transition appeared very cooperative when using high pressure at high temperature, as was the case for chemical denaturation, whereas it was found more progressive at low temperature, suggesting the existence of a complex folding pathway.

Molecular insights on poly(N-isopropylacrylamide) coil-to-globule transition induced by pressure

By using extensive all-atom molecular dynamics simulations of an atactic linear polymer chain, we provide microscopic insights into poly(N-isopropylacrylamide) (PNIPAM) coil-to-globule transition addressing the roles played by both temperature and pressure. We detect a coil-to-globule transition up to large pressures, showing a reentrant behavior of the critical temperature with increasing pressure in agreement with experimental observations. Furthermore, again confirming the experimental findings, we report the existence at high pressures of a new kind of globular state. It is characterized by a more structured hydration shell that is closer to PNIPAM hydrophobic domains, as compared to the globular state observed at atmospheric pressure. Our results highlight that temperature and pressure induce a PNIPAM coil-to-globule transition through different molecular mechanisms, opening the way for a systematic use of both thermodynamic variables to tune the location of the transition and the properties of the associated swollen/collapsed states.

Thermodynamic stability of hnRNP A1 low complexity domain revealed by high-pressure NMR

We have investigated the pressure- and temperature-induced conformational changes associated with the low complexity domain of hnRNP A1, an RNA-binding protein able to phase separate in response to cellular stress. Solution NMR spectra of the hnRNP A1 low-complexity domain fused with protein-G B1 domain were collected from 1 to 2500 bar and from 268 to 290 K. While the GB1 domain shows the typical pressure-induced and cold temperature-induced unfolding expected for small globular domains, the low-complexity domain of hnRNP A1 exhibits unusual pressure and temperature dependences. We observed that the low-complexity domain is pressure sensitive, undergoing a major conformational transition within the prescribed pressure range. Remarkably, this transition has the inverse temperature dependence of a typical folding-unfolding transition. Our results suggest the presence of a low-lying extended and fully solvated state(s) of the low-complexity domain that may play a role in phase separation. This study highlights the exquisite sensitivity of solution NMR spectroscopy to observe subtle conformational changes and illustrates how pressure perturbation can be used to determine the properties of metastable conformational ensembles.

Volume and compressibility differences between protein conformations revealed by high-pressure NMR

Proteins often interconvert between different conformations in ways critical to their function. Although manipulating such equilibria for biophysical study is often challenging, the application of pressure is a potential route to achieve such control by favoring the population of lower volume states. Here, we use this feature to study the interconversion of ARNT PAS-B Y456T, which undergoes a dramatic +3 slip in the β-strand register as it switches between two stably folded conformations. Using high-pressure biomolecular NMR approaches, we obtained the first, to our knowledge, quantitative data testing two key hypotheses of this process: the slipped conformation is both smaller and less compressible than the wild-type equivalent, and the interconversion proceeds through a chiefly unfolded intermediate state. Data collected in steady-state pressure and time-resolved pressure-jump modes, including observed pressure-dependent changes in the populations of the two conformers and increased rate of interconversion between conformers, support both hypotheses. Our work exemplifies how these approaches, which can be generally applied to protein conformational switches, can provide unique information that is not easily accessible through other techniques.

Enhanced accessibility and hydrophobicity of amyloidogenic intermediates of the β2-microglobulin D76N mutant revealed by high-pressure experiments

β2-Microglobulin (β2m) is the causative protein of dialysis-related amyloidosis. Its unfolding mainly proceeds along the pathway of N_C →U_C ⇄ U_T, whereas refolding follows the U_T → I_T (→N_T) →N_C pathway, in which N, I, and U are the native, intermediate, and unfolded states, respectively, with the Pro32 peptidyl-prolyl bond in cis or trans conformation as indicated by the subscript. It is noted that the I_T state is a putative amyloidogenic precursor state. Several aggregation-prone variants of β2m have been reported to date. One of these variants is D76N β2m, which is a naturally occurring amyloidogenic mutant. To elucidate the molecular mechanisms contributing to the enhanced amyloidogenicity of the mutant, we investigated the equilibrium and kinetic transitions of pressure-induced folding/unfolding equilibria in the wild type and D76N mutant by monitoring intrinsic tryptophan and 1-anilino-8-naphthalene sulfonate fluorescence. An analysis of kinetic data revealed that the different folding/unfolding behaviors of the wild type and D76N mutant were due to differences in the activation energy between the unfolded and the intermediate states as well as stability of the native state, leading to more rapid accumulation of I_T state for D76N in the refolding process. In addition, the I_T state was found to assume more hydrophobic nature. These changes induced the enhanced amyloidogenicity of the D76N mutant and the distinct pathogenic symptoms of patients. Our results suggest that the stabilization of the native state will be an effective approach for suppressing amyloid fibril formation of this mutant.

Combining High-Pressure Perturbation with NMR Spectroscopy for a Structural and Dynamical Characterization of Protein Folding Pathways

High-hydrostatic pressure is an alternative perturbation method that can be used to destabilize globular proteins. Generally perfectly reversible, pressure exerts local effects on regions or domains of a protein containing internal voids, contrary to heat or chemical denaturant that destabilize protein structures uniformly. When combined with NMR spectroscopy, high pressure (HP) allows one to monitor at a residue-level resolution the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. The use of HP-NMR has long been hampered by technical difficulties. Owing to the recent development of commercially available high-pressure sample cells, HP-NMR experiments can now be routinely performed. This review summarizes recent advances of HP-NMR techniques for the characterization at a quasi-atomic resolution of the protein folding energy landscape.

High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Deep sea biology is known to thrive at pressures up to ≈1 kbar, which motivates fundamental biophysical studies of biomolecules under such extreme environments. In this work, the conformational equilibrium of the lysine riboswitch has been systematically investigated by single molecule FRET (smFRET) microscopy at pressures up to 1500 bar. The lysine riboswitch preferentially unfolds with increasing pressure, which signals an increase in free volume (ΔV0 > 0) upon folding of the biopolymer. Indeed, the effective lysine binding constant increases quasi-exponentially with pressure rise, which implies a significant weakening of the riboswitch-ligand interaction in a high-pressure environment. The effects of monovalent/divalent cations and osmolytes on folding are also explored to acquire additional insights into cellular mechanisms for adapting to high pressures. For example, we find that although Mg2+ greatly stabilizes folding of the lysine riboswitch (ΔΔG0 < 0), there is negligible impact on changes in free volume (ΔΔV0 ≈ 0) and thus any pressure induced denaturation effects. Conversely, osmolytes (commonly at high concentrations in deep sea marine species) such as the trimethylamine N-oxide (TMAO) significantly reduce free volumes (ΔΔV0 < 0) and thereby diminish pressure-induced denaturation. We speculate that, besides stabilizing RNA structure, enhanced levels of TMAO in cells might increase the dynamic range for competent riboswitch folding by suppressing the pressure-induced denaturation response. This in turn could offer biological advantage for vertical migration of deep-sea species, with impacts on food searching in a resource limited environment.

Hydration properties of a protein at low and high pressures: Physics of pressure denaturation

Using experimentally determined structures of ubiquitin at 1 and 3000 bar, we generate sufficiently large ensembles of model structures in the native and pressure-induced (denatured) states by means of molecular dynamics simulations with explicit water. We calculate the values of a free-energy function (FEF), which comprises the hydration free energy (HFE) and the intramolecular (conformational) energy and entropy, for the two states at 1 and 3000 bar. The HFE and the conformational entropy, respectively, are calculated using our statistical-mechanical method, which has recently been shown to be accurate, and the Boltzmann-quasi-harmonic method. The HFE is decomposed into a variety of physically insightful components. We show that the FEF of the native state is lower than that of the denatured state at 1 bar, whereas the opposite is true at 3000 bar, thus being successful in reproducing the pressure denaturation. We argue that the following two quantities of hydration play essential roles in the denaturation: the WASA-dependent term in the water-entropy loss upon cavity creation for accommodating the protein (WASA is the water-accessible surface area of the cavity) and the protein-water Lennard-Jones interaction energy. At a high pressure, the mitigation of the serious water crowding in the system is the most important, and the WASA needs to be sufficiently enlarged with the increase in the excluded-volume being kept as small as possible. The denatured structure thus induced is characterized by the water penetration into the protein interior. The pressure denaturation is accompanied by a significantly large gain of water entropy.

Monitoring Unfolding of Titin I27 Single and Bi Domain with High-Pressure NMR Spectroscopy

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the folding energy landscape. Simulations, when corroborated by experimental data yielding global information on the folding process, can provide this level of insight. Molecular dynamics (MD) has often been combined with force spectroscopy experiments to decipher the unfolding mechanism of titin immunoglobulin-like single or multidomain, the giant multimodular protein from sarcomeres, yielding information on the sequential events during titin unfolding under stretching. Here, we used high-pressure NMR to monitor the unfolding of titin I27 Ig-like single domain and tandem. Because this method brings residue-specific information on the folding process, it can provide quasiatomic details on this process without the help of MD simulations. Globally, the results of our high-pressure analysis are in agreement with previous results obtained by the combination of experimental measurements and MD simulation and/or protein engineering, although the intermediate folding state caused by the early detachment of the AB β-sheet, often reported in previous works based on MD or force spectroscopy, cannot be detected. On the other hand, the A'G parallel β-sheet of the β-sandwich has been confirmed as the Achilles heel of the three-dimensional scaffold: its disruption yields complete unfolding with very similar characteristics (free energy, unfolding volume, kinetics rate constants) for the two constructs.

Folding of the Ig-Like Domain of the Dengue Virus Envelope Protein Analyzed by High-Hydrostatic-Pressure NMR at a Residue-Level Resolution

Dengue fever is a mosquito-borne endemic disease in tropical and subtropical regions, causing a significant public health problem in Southeast Asia. Domain III (ED3) of the viral envelope protein contains the two dominant putative epitopes and part of the heparin sulfate receptor binding region that drives the dengue virus (DENV)’s fusion with the host cell. Here, we used high-hydrostatic-pressure nuclear magnetic resonance (HHP-NMR) to obtain residue-specific information on the folding process of domain III from serotype 4 dengue virus (DEN4-ED3), which adopts the classical three-dimensional (3D) ß-sandwich structure known as the Ig-like fold. Interestingly, the folding pathway of DEN4-ED3 shares similarities with that of the Titin I27 module, which also adopts an Ig-like fold, but is functionally unrelated to ED3. For both proteins, the unfolding process starts by the disruption of the N- and C-terminal strands on one edge of the ß-sandwich, yielding a folding intermediate stable over a substantial pressure range (from 600 to 1000 bar). In contrast to this similarity, pressure-jump kinetics indicated that the folding transition state is considerably more hydrated in DEN4-ED3 than in Titin I27.

Thermodynamic properties of aqueous osmolyte solutions at high-pressure conditions

Living organisms can be encountered in nature under extreme conditions. At the seabed, pressure may reach 1000 bar. Yet microorganisms can be found that still function under these conditions. On the one hand, it is known that high pressure even has a positive effect on piezophile enzymes increasing their activity. On the other hand, such microorganisms might contain up to very high concentrations of osmolytes that counteract osmotic stress. To better understand high-pressure influences on biochemical systems, fundamental knowledge about pressure effects on thermodynamic properties of such osmolytes is important. However, literature data is scarce and experiments at high-pressure conditions are challenging. Hence, new high-pressure density data of aqueous osmolyte solutions were measured in this work at temperatures between 298.15 K and 318.15 K and at osmolyte concentrations up to 3 mol/kg water. Further, the thermodynamic model PC-SAFT has been applied recently to successfully model vapor pressures of water and density of water up to 10 kbar [M. Knierbein et al., Density variations of TMAO solutions in the kilobar range: experiments, PC-SAFT predictions, and molecular dynamics simulations, Biophysical chemistry, (2019)]. This allowed accurately predicting effects of temperature and osmolyte concentration on thermodynamic properties (especially mixture densities) up to very high pressures. Common osmolytes (trimethylamine-N-oxide, urea, ectoine, glycerol, glycine) as well as the dipeptides acetyl-N-methylglycine amide, acetyl-N-methylalanine amide, and acetyl-N-methylleucine amide were under investigation.

Molecular-Level Profiling of Human Serum Transferrin Protein through Assessment of Nanopore-Based Electrical and Chemical Responsiveness

In this study, we investigated the voltage and pH responsiveness of human serum transferrin (hSTf) protein using silicon nitride (Si _xN _y) nanopores. The Fe(III)-rich holo form of hSTf was dominant when pH > pI, while the Fe(III)-free apo form was dominant when pH < pI. The translocations of hSTf were purely in an electrophoretic sense, thus depended on its pI and the solution pH. With increasing voltage, voltage driven protein unfolding became prominent which was indicated by the trends associated with change in conductance, due to hSTf translocation, and in the excluded electrolyte volume. Additionally, analysis of the translocation events of the pure apo form of hSTf showed a clear difference in the event population compared to that of the holo form. The results obtained demonstrate the successful application of nanopore devices to distinguish between the holo and apo forms of hSTf in a mixture and to analyze its folding and unfolding phenomenon over a range of pH and applied voltages.

The role of hydrophobicity in the cold denaturation of proteins under high pressure: A study on apomyoglobin

An exciting debate arises when microscopic mechanisms involved in the denaturation of proteins at high pressures are explained. In particular, the issue emerges when the hydrophobic effect is invoked, given that hydrophobicity cannot elucidate by itself the volume changes measured during protein unfolding. In this work, we study by the use of molecular dynamics simulations and essential dynamics analysis the relation between the solvation dynamics, volume, and water structure when apomyoglobin is subjected to a hydrostatic pressure regime. Accordingly, the mechanism of cold denaturation of proteins under high-pressure can be related to the disruption of the hydrogen-bond network of water favoring the coexistence of two states, low-density and high-density water, which directly implies in the formation of a molten globule once the threshold of 200 MPa has been overcome.

Lessons from pressure denaturation of proteins

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

Pressure and cosolvent modulation of the catalytic activity of amyloid fibrils

We report on the effects of pressure and cosolvents on the catalytic activity of a designed amyloid fibril by applying a high-pressure stopped-flow methodology with rapid spectroscopic detection. FTIR spectroscopic data revealed a remarkable pressure and temperature stability of the fibrillar catalyst. The activity is further enhanced by osmolytes and macromolecular crowding.

Pressure Effect on Zero-Field Splitting Parameter of Hemin: Model Case of Hemoproteins under Pressure

We experimentally studied the pressure dependence of the zero-field splitting (ZFS) parameter of hemin (iron(III) protoporphyrin IX chloride), which is a model complex of hemoproteins, via high-frequency and high-field electron paramagnetic resonance (HFEPR) under pressure. Owing to the large ZFS, the pressure effect on the electronic structure of iron-porphyrin complexes has not yet been explored using EPR. Therefore, we systematically studied this effect using our newly developed sub-terahertz EPR spectroscopy system in the frequency range of 80-515 GHz, under magnetic fields up to 10 T and pressure up to 2 GPa. We observed a systematic shift of the resonance fields of hemin upon pressure application, from which the axial component of the ZFS parameter was found to increase from D = 6.9 to 7.9 cm^-1 at 2 GPa. In contrast to the previous methods used to study proteins under pressure, which mainly focused on conformational changes, our HFEPR technique can obtain more microscopic insights into the electronic structures of metal ions under pressure. In this sense, our technique provides novel opportunities to study the pressure effects on biofunctional active centers of versatile metalloproteins.

Theoretical analysis of high-pressure effects on conformational equilibria

Along with temperature, pressure is the most important physical parameter determining the thermodynamic properties and reactivity of chemical systems. In this work, we discuss the effects of high pressure on conformational properties of organic molecules and propose an approach toward calculation of conformational volume changes based on molecular dynamics simulations. The results agree well with the experimental data. Furthermore, we demonstrate that pressure can be used as an instrument for fine-tuning of molecular conformations and to propel a properly constructed molecular rotor possessing a suitable combination of energy and volume profiles.

The effects of glycine, TMAO and osmolyte mixtures on the pressure dependent enzymatic activity of α-chymotrypsin

Different natural osmolytes modulate the pressure dependent enzyme kinetics in different ways.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

Excess Volumes from the Pressure Derivative of the Excess Chemical Potential: Testing Simple Models for Cavity Formation in Water

Excess volumes and excess compressibilities for hard spheres in water were computed by pressure derivatives of the excess chemical potential, which is equivalent to the work of cavity formation. This is relevant to the application of continuum solvation methods at various pressures. The excess chemical potential was modeled within phenomenological expressions for curved surfaces plus a pressure-volume term, for which two approaches were adopted, differing for the radius of the spherical volume. This implies a different dependence on pressure of parameters. In all cases, in the surface term, for the pressure derivative of parameters of the curvature function, use was made of the previously proposed expressions for the first two moments obtained from the density and radial distribution of oxygens in liquid water. Only for the parameter which has the dimension of surface tension (γ̃) was explicit dependence on pressure considered and results are affected by the specific polynomial used. In agreement with what inferred from simulation results obtained for cavities in TIP4P water, negative and positive adsorptions at the contact radius were extrapolated for a very large cavity at 1 and 8000 atm, respectively. The expressions here employed for the excess chemical potential predict the zero value of asymptotic adsorption to be at a pressure between 500 and 800 atm, which can be compared to results from the revised scaled particle theory. In the same range, for a nanometer-sized cavity, a change of behavior occurs regarding the ratio between the excess Helmholtz free energy and the product between pressure and excess volume.

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 ¹H-¹⁵N 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.

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 (ΔG_f°) 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 (ΔV_f) was significantly larger for V67 (71 ± 2 mL mol^-1) than for L67 (58 ± 3 mL mol^-1) inteins. The measured difference in ΔV_f (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-pressure NMR techniques for the study of protein dynamics, folding and aggregation

High-pressure is a well-known perturbation method used to destabilize globular proteins and dissociate protein complexes or aggregates. The heterogeneity of the response to pressure offers a unique opportunity to dissect the thermodynamic contributions to protein stability. In addition, pressure perturbation is generally reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility of monitoring at an atomic resolution the structural transitions occurring upon unfolding and determining the kinetic properties of the process. The recent development of commercially available high-pressure sample cells greatly increased the potential applications for high-pressure NMR experiments that can now be routinely performed. This review summarizes the recent applications and future directions of high-pressure NMR techniques for the characterization of protein conformational fluctuations, protein folding and the stability of protein complexes and aggregates.

Native State Volume Fluctuations in Proteins as a Mechanism for Dynamic Allostery

Allostery enables tight regulation of protein function in the cellular environment. Although existing models of allostery are firmly rooted in the current structure-function paradigm, the mechanistic basis for allostery in the absence of structural change remains unclear. In this study, we show that a typical globular protein is able to undergo significant changes in volume under native conditions while exhibiting no additional changes in protein structure. These native state volume fluctuations were found to correlate with changes in internal motions that were previously recognized as a source of allosteric entropy. This finding offers a novel mechanistic basis for allostery in the absence of canonical structural change. The unexpected observation that function can be derived from expanded, low density protein states has broad implications for our understanding of allostery and suggests that the general concept of the native state be expanded to allow for more variable physical dimensions with looser packing.

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.

Equilibria of oligomeric proteins under high pressure - A theoretical description

High pressure methods have become a useful tool for studying protein structure and stability. Using them, various physico-chemical processes including protein unfolding, aggregation, oligomer dissociation or enzyme-activity decrease were studied on many different proteins. Oligomeric protein dissociation is a process that can perfectly utilize the potential of high-pressure techniques, as the high pressure shifts the equilibria to higher concentrations making them better observable by spectroscopic methods. This can be especially useful when the oligomeric form is highly stable at atmospheric pressure. These applications may be, however, hindered by less intensive experimental response as well as interference of the oligomerization equilibria with unfolding or aggregation of the subunits, but also by more complex theoretical description. In this study we develop mathematical models describing different kinds of oligomerization equilibria, both closed (equilibrium of monomer and the highest possible oligomer without any intermediates) and consecutive. Closed homooligomer equilibria are discussed for any oligomerization degree, while the more complex heterooligomer equilibria and the consecutive equilibria in both homo- and heterooligomers are taken into account only for dimers and trimers. In all the cases, fractions of all the relevant forms are evaluated as functions of pressure and concentration. Significant points (inflection points and extremes) of the resulting transition curves, that can be determined experimentally, are evaluated as functions of pressure and/or concentration. These functions can be further used in order to evaluate the thermodynamic parameters of the system, i.e. atmospheric-pressure equilibrium constants and volume changes of the individual steps of the oligomer-dissociation processes.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.

The energetics of a three-state protein folding system probed by high-pressure relaxation dispersion NMR spectroscopy

The energetic and volumetric properties of a three-state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure-dependent (15) N-relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT ) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol(-1) and 40 μL mol(-1)  bar(-1) , respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on-pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter.

The unfolding effects on the protein hydration shell and partial molar volume: a computational study

In this paper we apply the computational analysis recently proposed by our group to characterize the solvation properties of a native protein in aqueous solution, and to four model aqueous solutions of globular proteins in their unfolded states thus characterizing the protein unfolded state hydration shell and quantitatively evaluating the protein unfolded state partial molar volumes.