Protein compressibility, dynamics, and pressure

Differential Responses in the Core, Active Site and Peripheral Regions of Cytochrome c Peroxidase to Extreme Pressure and Temperature

In consideration of life in extreme environments, the effects of hydrostatic pressure on proteins at the atomic level have drawn substantial interest. Large deviations of temperature and pressure from ambient conditions can shift the free energy landscape of proteins to reveal otherwise lowly populated structural states and even promote unfolding. We report the crystal structure of the heme-containing peroxidase, cytochrome c peroxidase (CcP) at 1.5 and 3.0 kbar and make comparisons to structures determined at 1.0 bar and cryo-temperatures (100 K). Pressure produces anisotropic changes in CcP, but compressibility plateaus after 1.5 kbar. CcP responds to pressure with volume declines at the periphery of the protein where B-factors are relatively high but maintains nearly intransient core structure, hydrogen bonding interactions and active site channels. Changes in active-site solvation and heme ligation reveal pressure sensitivity to protein-ligand interactions and a potential docking site for the substrate peroxide. Compression at the surface affects neither alternate side-chain conformers nor B-factors. Thus, packing in the core, which resembles a crystalline solid, limits motion and protects the active site, whereas looser packing at the surface preserves side-chain dynamics. These data demonstrate that conformational dynamics and packing densities are not fully correlated in proteins and that encapsulation of cofactors by the polypeptide can provide a precisely structured environment resistant to change across a wide range of physical conditions.

Human cytochrome C natural variants: Studying the membrane binding properties of G41S and Y48H by fluorescence energy transfer and molecular dynamics

Cytochrome C (cyt C), the protein involved in oxidative phosphorylation, plays several other crucial roles necessary for both cell life and death. Studying natural variants of cyt C offers the possibility to better characterize the structure-to-function relationship that modulates the different activities of this protein. Naturally mutations in human cyt C (G41S and Y48H) occur in the protein central Ω-loop and cause thrombocytopenia 4. In this study, we have investigated the binding of such variants and of wild type (wt) cyt C to synthetic cardiolipin-containing vesicles. The mutants have a lower propensity in membrane binding, displaying higher dissociation constants with respect to the wt protein. Compressibility measurements reveal that both variants are more flexible than the wt, suggesting that the native central Ω-loop is important for the interaction with membranes. Such hypothesis is supported by molecular dynamics simulations. A minimal distance analysis indicates that in the presence of cardiolipin the central Ω-loop of the mutants is no more in contact with the membrane, as it happens instead in the case of wt cyt C. Such finding might provide a hint for the reduced membrane binding capacity of the variants and their enhanced peroxidase activity in vivo.

Phase separation of YAP-MAML2 differentially regulates the transcriptome

Phase separation (PS) drives the formation of biomolecular condensates that are emerging biological structures involved in diverse cellular processes. Recent studies have unveiled PS-induced formation of several transcriptional factor (TF) condensates that are transcriptionally active, but how strongly PS promotes gene activation remains unclear. Here, we show that the oncogenic TF fusion Yes-associated protein 1-Mastermind like transcriptional coactivator 2 (YAP-MAML2) undergoes PS and forms liquid-like condensates that bear the hallmarks of transcriptional activity. Furthermore, we examined the contribution of PS to YAP-MAML2-mediated gene expression by developing a chemogenetic tool that dissolves TF condensates, allowing us to compare phase-separated and non-phase-separated conditions at identical YAP-MAML2 protein levels. We found that a small fraction of YAP-MAML2-regulated genes is further affected by PS, which include the canonical YAP target genes CTGF and CYR61, and other oncogenes. On the other hand, majority of YAP-MAML2-regulated genes are not affected by PS, highlighting that transcription can be activated effectively by diffuse complexes of TFs with the transcriptional machinery. Our work opens new directions in understanding the role of PS in selective modulation of gene expression, suggesting differential roles of PS in biological processes.

Inline small-angle X-ray scattering-coupled chromatography under extreme hydrostatic pressure

As continuing discoveries highlight the surprising abundance and resilience of deep ocean and subsurface microbial life, the effects of extreme hydrostatic pressure on biological structure and function have attracted renewed interest. Biological small-angle X-ray scattering (BioSAXS) is a widely used method of obtaining structural information from biomolecules in solution under a wide range of solution conditions. Due to its ability to reduce radiation damage, remove aggregates, and separate monodisperse components from complex mixtures, size-exclusion chromatography-coupled SAXS (SEC-SAXS) is now the dominant form of BioSAXS at many synchrotron beamlines. While BioSAXS can currently be performed with some difficulty under pressure with non-flowing samples, it has not been clear how, or even if, continuously flowing SEC-SAXS, with its fragile media-packed columns, might work in an extreme high-pressure environment. Here we show, for the first time, that reproducible chromatographic separations coupled directly to high-pressure BioSAXS can be achieved at pressures up to at least 100 MPa and that pressure-induced changes in folding and oligomeric state and other properties can be observed. The apparatus described here functions at a range of temperatures (0°C-50°C), expanding opportunities for understanding biomolecular rules of life in deep ocean and subsurface environments.

Unravelling the Adaptation Mechanisms to High Pressure in Proteins

Life is thought to have appeared in the depth of the sea under high hydrostatic pressure. Nowadays, it is known that the deep biosphere hosts a myriad of life forms thriving under high-pressure conditions. However, the evolutionary mechanisms leading to their adaptation are still not known. Here, we show the molecular bases of these mechanisms through a joint structural and dynamical study of two orthologous proteins. We observed that pressure adaptation involves the decoupling of protein–water dynamics and the elimination of cavities in the protein core. This is achieved by rearranging the charged residues on the protein surface and using bulkier hydrophobic residues in the core. These findings will be the starting point in the search for a complete genomic model explaining high-pressure adaptation.

Volumetric Properties of Four-Stranded DNA Structures

Four-stranded non-canonical DNA structures including G-quadruplexes and i-motifs have been found in the genome and are thought to be involved in regulation of biological function. These structures have been implicated in telomere biology, genomic instability, and regulation of transcription and translation events. To gain an understanding of the molecular determinants underlying the biological role of four-stranded DNA structures, their biophysical properties have been extensively studied. The limited libraries on volume, expansibility, and compressibility accumulated to date have begun to provide insights into the molecular origins of helix-to-coil and helix-to-helix conformational transitions involving four-stranded DNA structures. In this article, we review the recent progress in volumetric investigations of G-quadruplexes and i-motifs, emphasizing how such data can be used to characterize intra-and intermolecular interactions, including solvation. We describe how volumetric data can be interpreted at the molecular level to yield a better understanding of the role that solute-solvent interactions play in modulating the stability and recognition events of nucleic acids. Taken together, volumetric studies facilitate unveiling the molecular determinants of biological events involving biopolymers, including G-quadruplexes and i-motifs, by providing one more piece to the thermodynamic puzzle describing the energetics of cellular processes in vitro and, by extension, in vivo.

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.

Effect of X-ray free-electron laser-induced shockwaves on haemoglobin microcrystals delivered in a liquid jet

X-ray free-electron lasers (XFELs) enable obtaining novel insights in structural biology. The recently available MHz repetition rate XFELs allow full data sets to be collected in shorter time and can also decrease sample consumption. However, the microsecond spacing of MHz XFEL pulses raises new challenges, including possible sample damage induced by shock waves that are launched by preceding pulses in the sample-carrying jet. We explored this matter with an X-ray-pump/X-ray-probe experiment employing haemoglobin microcrystals transported via a liquid jet into the XFEL beam. Diffraction data were collected using a shock-wave-free single-pulse scheme as well as the dual-pulse pump-probe scheme. The latter, relative to the former, reveals significant degradation of crystal hit rate, diffraction resolution and data quality. Crystal structures extracted from the two data sets also differ. Since our pump-probe attributes were chosen to emulate EuXFEL operation at its 4.5 MHz maximum pulse rate, this prompts concern about such data collection.

High-pressure small-angle X-ray scattering cell for biological solutions and soft materials

Pressure is a fundamental thermodynamic parameter controlling the behavior of biological macromolecules. Pressure affects protein denaturation, kinetic parameters of enzymes, ligand binding, membrane permeability, ion trans-duction, expression of genetic information, viral infectivity, protein association and aggregation, and chemical processes. In many cases pressure alters the molecular shape. Small-angle X-ray scattering (SAXS) is a primary method to determine the shape and size of macromolecules. However, relatively few SAXS cells described in the literature are suitable for use at high pressures and with biological materials. Described here is a novel high-pressure SAXS sample cell that is suitable for general facility use by prioritization of ease of sample loading, temperature control, mechanical stability and X-ray background minimization. Cell operation at 14 keV is described, providing a q range of 0.01 < q < 0.7 Å-1, pressures of 0-400 MPa and an achievable temperature range of 0-80°C. The high-pressure SAXS cell has recently been commissioned on the ID7A beamline at the Cornell High Energy Synchrotron Source and is available to users on a peer-reviewed proposal basis.

Shape-preserving elastic solid models of macromolecules

Mass-spring models have been a standard approach in molecular modeling for the last few decades, such as elastic network models (ENMs) that are widely used for normal mode analysis. In this work, we present a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as ENMs in producing the equilibrium dynamics and moreover, offers some significant new features that may greatly benefit the research community. ESM is different from ENM in that it treats macromolecules as elastic solids. Our particular version of ESM presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape and thus makes normal mode computations and visualization of extremely large complexes more manageable. Secondly, as a solid model, ESM’s close link to finite element analysis renders it ideally suited for studying mechanical responses of macromolecules under external force. Lastly, we show that ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy. The complete MATLAB code of αESM is provided.

Mass-spring models have been a standard approach in classical molecular modeling where atoms are modeled as spheres with a mass and their interactions modeled as springs. The models have been extremely successful. Thinking ahead, however, as molecular systems of our interest grow more quickly in size or dimension than what our computation resources can keep up with, some adjustments in methodology are timely. This work presents a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as mass-spring models in producing the equilibrium dynamics and moreover, offers some unique features that make it suitable for much larger systems. ESM is different from ENMs in that it treats macromolecules as elastic solids. Our particular version of ESM model presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape. ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (R_h) of biomolecules in solution. For unfolded proteins, the R_h provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the R_h of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by ¹³C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that R_h values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the R_h of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Porcine parvovirus VP1/VP2 on a time series epitope mapping: exploring the effects of high hydrostatic pressure on the immune recognition of antigens

Porcine parvovirus (PPV) is a DNA virus that causes reproductive failure in gilts and sows, resulting in embryonic and fetal losses worldwide. Epitope mapping of PPV is important for developing new vaccines. In this study, we used spot synthesis analysis for epitope mapping of the capsid proteins of PPV (NADL-2 strain) and correlated the findings with predictive data from immunoinformatics. The virus was exposed to three conditions prior to inoculation in pigs: native (untreated), high hydrostatic pressure (350 MPa for 1 h) at room temperature and high hydrostatic pressure (350 MPa for 1 h) at − 18 °C, and was compared with a commercial vaccine produced using inactivated PPV. The screening of serum samples detected 44 positive spots corresponding to 20 antigenic sites. Each type of inoculated antigen elicited a distinct epitope set. In silico prediction located linear and discontinuous epitopes in B cells that coincided with several epitopes detected in spot synthesis of sera from pigs that received different preparations of inoculum. The conditions tested elicited antibodies against the VP1/VP2 antigen that differed in relation to the response time and the profile of structurally available regions that were recognized.

On empirical decomposition of volumetric data

Volumetric characterization of proteins and their recognition events has been instrumental in providing information on the role of intra- and intermolecular interactions, including hydration, in stabilizing biomolecules. The credibility of molecular models and interpretation schemes used to rationalize experimental data are essential for the validity of microscopic insights derived from volumetric results. Current empirical schemes used to interpret volumetric data suffer from a lack of theoretical and computational substantiation. In this contribution, we take advantage age of recent MD simulations of proteins in solution coupled with Voronoi-Delaunay tessellation of simulated structures that have provided an exceptional level of structural detail on the nature of protein-water interfaces. We use these structural insights to re-evaluate empirical frameworks used for interpretation of volumetric data. An important issue in this respect is the actual dividing surface between water and protein atoms that is used in volumetric studies when the solute and solvent are treated as hard spheres enclosed within their respective van der Waals surfaces. In one development, using Voronoi tessellation of MD simulated protein-water systems the dividing surface has been defined as the points equidistant from the water and protein atoms. The interstitial void volume between the solute and the dividing surface corresponds to thermal volume envisaged by Scaled Particle Theory. In this communication, we explicitly account for the contributions of thermal volume to the partial molar volume, compressibility, and expansibility of proteins and re-examine and redefine the intrinsic and hydration volumetric contributions. We discuss the implications of our results for protein transitions and association events.

Enzymes from piezophiles

The discovery of microbial communities in extreme conditions that would seem hostile to life leads to the question of how the molecules making up these microbes can maintain their structure and function. While microbes that live under extremes of temperature have been heavily studied, those that live under extremes of pressure, or "piezophiles", are now increasingly being studied because of advances in sample collection and high-pressure cells for biochemical and biophysical measurements. Here, adaptations of enzymes in piezophiles against the effects of pressure are discussed in light of recent experimental and computational studies. However, while concepts from studies of enzymes from temperature extremophiles can provide frameworks for understanding adaptations by piezophile enzymes, the effects of temperature and pressure on proteins differ in significant ways. Thus, the state of the knowledge of adaptation in piezophile enzymes is still in its infancy and many more experiments and computational studies on different enzymes from a variety of piezophiles are needed.

Experimental Investigation of Triplet Correlation Approximations for Fluid Water

Triplet correlations play a central role in our understanding of fluids and their properties. Of particular interest is the relationship between the pair and triplet correlations. Here we use a combination of Fluctuation Solution Theory and experimental pair radial distribution functions to investigate the accuracy of the Kirkwood Superposition Approximation (KSA), as given by integrals over the relevant pair and triplet correlation functions, at a series of state points for pure water using only experimental quantities. The KSA performs poorly, in agreement with a variety of other studies. Several additional approximate relationships between the pair and triplet correlations in fluids are also investigated and generally provide good agreement for the fluid thermodynamics for regions of the phase diagram where the compressibility is small. A simple power law relationship between the pair and triplet fluctuations is particularly successful for state points displaying low to moderately high compressibilities.

Quasiharmonic analysis of protein energy landscapes from pressure-temperature molecular dynamics simulations

Positional fluctuations of an atom in a protein can be described as motion in an effective local energy minimum created by the surrounding protein atoms. The dependence of atomic fluctuations on both temperature (T) and pressure (P) has been used to probe the nature of these minima, which are generally described as harmonic in experiments such as x-ray crystallography and neutron scattering. Here, a quasiharmonic analysis method is presented in which the P-T dependence of atomic fluctuations is in terms of an intrinsic isobaric thermal expansivity α_P and an intrinsic isothermal compressibility κ_T. The method is tested on previously reported mean-square displacements from P-T molecular dynamics simulations of lysozyme, which were interpreted to have a pressure-independent dynamical transition T_g near 200 K and a change in the pressure dependence near 480 MPa. Our quasiharmonic analysis of the same data shows that the P-T dependence can be described in terms of α_P and κ_T where below T_g, the temperature dependence is frozen at the T_g value. In addition, the purported transition at 480 MPa is reinterpreted as a consequence of the pressure dependence of T_g and the quasiharmonic frequencies. The former also indicates that barrier heights between substates are pressure dependent in these data. Furthermore, the insights gained from this quasiharmonic analysis, which was of the energy landscape near the native state of a protein, suggest that similar analyses of other simulations may be useful in understanding such phenomena as pressure-induced protein unfolding.

High pressure Raman study of type-I collagen

The high pressure response of type-I collagen from bovine Achilles tendon is investigated with micro-Raman spectroscopy. Fluorinert™ and methanol-ethanol mixtures were used as pressure transmitting media (PTM) in a diamond anvil cell. The Raman spectrum of collagen is dominated by three bands centred at approximately 1450, 1660 and 2930 cm^-1 , attributed to C-H deformation, C=O stretching of the peptide bond (amide-I band) and C-H stretching modes respectively. Upon pressure increase, using Fluorinert™ as PTM, a shift towards higher frequencies of the C-H stretching and deformation peaks is observed. Contrary, the amide-I band peaks are shifted to lower frequencies with moderate pressure slopes. On the other hand, when using the alcohol mixture as PTM, the amide-I band exhibits more pronounced C=O bond softening, deduced from the shift to lower frequencies, suggesting a strengthening of the hydrogen bonds between glycine and proline residues of different collagen chains due to the presence of the polar alcohol molecules. Furthermore, some of the peaks exhibit abrupt changes in their pressure slopes at approximately 2 GPa, implying a variation in the compressibility of the collagen fibres. This could be attributed to a pitch change from 10/3 to 7/2, sliding of the tropocollagen molecules, twisting variation at the molecular level and/or elimination of the D-gaps induced by kink compression. All spectral changes are reversible upon pressure release, which indicates that denaturation has not taken place. Finally, a minor lipid phase contamination was detected in some sample spots. Its pressure response is also monitored.

Dissimilar flexibility of α and β subunits of human adult hemoglobin influences the protein dynamics and its alteration induced by allosteric effectors

The general question by what mechanism an "effector" molecule and the hemes of hemoglobin interact over widely separated intramolecular distances to change the oxygen affinity has been extensively investigated, and still has remained of central interest. In the present work we were interested in clarifying the general role of the protein matrix and its dynamics in the regulation of human adult hemoglobin (HbA). We used a spectroscopy approach that yields the compressibility (κ) of the protein matrix around the hemes of the subunits in HbA and studied how the binding of heterotropic allosteric effectors modify this parameter. κ is directly related to the variance of volume fluctuation, therefore it characterizes the molecular dynamics of the protein structure. For the experiments the heme groups either in the α or in the β subunits of HbA were replaced by fluorescent Zn-protoporphyrinIX, and series of fluorescence line narrowed spectra were measured at varied pressures. The evaluation of the spectra yielded the compressibility that showed significant dynamic asymmetry between the subunits: κ of the α subunit was 0.17±0.05/GPa, while for the β subunit it was much higher, 0.36±0.07/GPa. The heterotropic effectors, chloride ions, inositol hexaphosphate and bezafibrate did not cause significant changes in κ of the α subunits, while in the β subunits the effectors lead to a significant reduction down to 0.15±0.04/GPa. We relate our results to structural data, to results of recent functional studies and to those of molecular dynamics simulations, and find good agreements. The observed asymmetry in the flexibility suggests a distinct role of the subunits in the regulation of Hb that results in the observed changes of the oxygen binding capability.

Simulated pressure denaturation thermodynamics of ubiquitin

Simulations of protein thermodynamics are generally difficult to perform and provide limited information. It is desirable to increase the degree of detail provided by simulation and thereby the potential insight into the thermodynamic properties of proteins. In this study, we outline how to analyze simulation trajectories to decompose conformation-specific, parameter free, thermodynamically defined protein volumes into residue-based contributions. The total volumes are obtained using established methods from Fluctuation Solution Theory, while the volume decomposition is new and is performed using a simple proximity method. Native and fully extended ubiquitin are used as the test conformations. Changes in the protein volumes are then followed as a function of pressure, allowing for conformation-specific protein compressibility values to also be obtained. Residue volume and compressibility values indicate significant contributions to protein denaturation thermodynamics from nonpolar and coil residues, together with a general negative compressibility exhibited by acidic residues.

Electronic characterization of Geobacter sulfurreducens pilins in self-assembled monolayers unmasks tunnelling and hopping conduction pathways

The peptide subunit of Geobacter nanowires (pili) metal-reducing bacterium Geobacter sulfurreducens was self-assembled as a conductive monolayer. Its electronic characterized revealed tunneling and hopping regimes.

Prediction of stability changes upon mutation in an icosahedral capsid

Identifying the contributions to thermodynamic stability of capsids is of fundamental and practical importance. Here we use simulation to assess how mutations affect the stability of lumazine synthase from the hyperthermophile Aquifex aeolicus, a T = 1 icosahedral capsid; in the simulations the icosahedral symmetry of the capsid is preserved by simulating a single pentamer and imposing crystal symmetry, in effect simulating an infinite cubic lattice of icosahedral capsids. The stability is assessed by estimating the free energy of association using an empirical method previously proposed to identify biological units in crystal structures. We investigate the effect on capsid formation of seven mutations, for which it has been experimentally assessed whether they disrupt capsid formation or not. With one exception, our approach predicts the effect of the mutations on the capsid stability. The method allows the identification of interaction networks, which drive capsid assembly, and highlights the plasticity of the interfaces between subunits in the capsid. Proteins 2015; 83:1733–1741. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc

Experimental triplet and quadruplet fluctuation densities and spatial distribution function integrals for pure liquids

Fluctuation solution theory has provided an alternative view of many liquid mixture properties in terms of particle number fluctuations. The particle number fluctuations can also be related to integrals of the corresponding two body distribution functions between molecular pairs in order to provide a more physical picture of solution behavior and molecule affinities. Here, we extend this type of approach to provide expressions for higher order triplet and quadruplet fluctuations, and thereby integrals over the corresponding distribution functions, all of which can be obtained from available experimental thermodynamic data. The fluctuations and integrals are then determined using the International Association for the Properties of Water and Steam Formulation 1995 (IAPWS-95) equation of state for the liquid phase of pure water. The results indicate small, but significant, deviations from a Gaussian distribution for the molecules in this system. The pressure and temperature dependence of the fluctuations and integrals, as well as the limiting behavior as one approaches both the triple point and the critical point, are also examined.

Volume and Compressibility of Proteins

The partial specific (or molar) volume, expansibility, and compressibility of a protein are fundamental thermodynamic quantities for characterizing its structure in solution. We review the definitions, measurements, and implications of these volumetric quantities in relation to protein structural biology. The partial specific volumes under constant molality (isomolal) and chemical potential (isopotential) conditions of the cosolvent (multicomponent systems) are explained in terms of preferential solvent interactions relevant to the solubility and stability of proteins. The partial expansibility is briefly discussed in terms of the effects of temperature on protein-solvent interactions (hydration) and internal packing defects (cavities). We discuss the compressibility-structure-function relationships of proteins based on analyses of the correlations between the partial adiabatic compressibilities and the structures or functions of various globular proteins (including mutants), focusing on the roles of the internal cavities in structural fluctuations. The volume and compressibility changes associated with various conformational transitions are also discussed in terms of the changes in hydration and cavities in order to elucidate the nonnative structures and the transition mechanisms, especially those associated with pressure denaturation.

Driving Forces in Pressure-Induced Protein Transitions

The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when exposed to water in an unfolded conformation. This chapter presents a volumetric analysis of the peptide group and the 20 naturally occurring amino acid side chains in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and in the unfolded conformation modeled by low-molecular analogs of proteins. The transfer of a peptide group from the protein interior to water becomes increasingly favorable as pressure increases. This observation classifies solvation of peptide groups as a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. An inference can be drawn that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures this core, owing to the absence of structural constraints, may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced protein denaturation. Volumetric data presented here have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensembles.

Pressure-Dependent Conformation and Fluctuation in Folded Protein Molecules

Hydrostatic pressure leads to nonuniform compression of proteins. The structural change is on average only about 0.1 Å kbar(-1), and is therefore within the range of fluctuations at ambient pressure. The largest changes are around cavities and buried water molecules. Sheets distort much more than helices. Hydrogen bonds compress about 0.012 Å kbar(-1), although there is a wide range, including some hydrogen bonds that lengthen. In the presence of ligands and inhibitors, structural changes are smaller. Pressure has little effect on rapid fluctuations, but with larger scale slower motions, pressure increases the population of excited states (if they have smaller overall volume), and slows the fluctuations. In barnase, pressure is shown to be a useful way to characterise fluctuations on the timescale of microseconds, and helps to show that fluctuations in barnase are hierarchical, with the faster fluctuations providing a platform for the slower ones. The excited states populated at high pressure are probably functionally important.

Transient access to the protein interior: simulation versus NMR

Many proteins rely on rare structural fluctuations for their function, whereby solvent and other small molecules gain transient access to internal cavities. In magnetic relaxation dispersion (MRD) experiments, water molecules buried in such cavities are used as intrinsic probes of the intermittent protein motions that govern their exchange with external solvent. While this has allowed a detailed characterization of exchange kinetics for several proteins, little is known about the exchange mechanism. Here, we use a millisecond all-atom MD trajectory produced by Shaw et al. (Science2010, 330, 341) to characterize water exchange from the four internal hydration sites in the protein bovine pancreatic trypsin inhibitor. Using a recently developed stochastic point process approach, we compute the survival correlation function probed by MRD experiments as well as other quantities designed to validate the exchange-mediated orientational randomization (EMOR) model used to interpret the MRD data. The EMOR model is found to be quantitatively accurate, and the simulation reproduces the experimental mean survival times for all four sites with activation energy discrepancies in the range 0-3 kBT. On the other hand, the simulated hydration sites are somewhat too flexible, and the water flip barrier is underestimated by up to 6 kBT. The simulation reveals that water molecules gain access to the internal sites by a transient aqueduct mechanism, migrating as single-file water chains through transient (<5 ns) tunnels or pores. The present study illustrates the power of state-of-the-art molecular dynamics simulations in validating and extending experimental results.

Pressure-induced changes in the fluorescence behavior of red fluorescent proteins

We present an experimental study on the fluorescence behavior of the red fluorescent proteins TagRFP-S, TagRFP-T, mCherry, mOrange2, mStrawberry, and mKO as a function of pressure up to several GPa. TagRFP-S, TagRFP-T, mOrange2, and mStrawberry show an initial increase in fluorescence intensity upon application of pressure above ambient conditions. At higher pressures, the fluorescence intensity decreases dramatically for all proteins under study, probably due to denaturing of the proteins. Small blue shifts in the fluorescence spectra with increasing pressure were seen in all proteins under study, hinting at increased rigidity of the chromophore environment. In addition, mOrange2 and mStrawberry exhibit strong and abrupt changes in their fluorescence spectra at certain pressures. These changes are likely due to structural modifications of the hydrogen bonding environment of the chromophore. The strong differences in behavior between proteins with identical or very similar chromophores highlight how the chromophore environment contributes to pressure-induced behavior of the fluorescence performance.

On the role of the difference in surface tensions involved in the allosteric regulation of NHE-1 induced by low to mild osmotic pressure, membrane tension and lipid asymmetry

The sodium-proton exchanger 1 (NHE-1) is a membrane transporter that exchanges Na(+) for H(+) ion across the membrane of eukaryotic cells. It is cooperatively activated by intracellular protons, and this allosteric regulation is modulated by the biophysical properties of the plasma membrane and related lipid environment. Consequently, NHE-1 is a mechanosensitive transporter that responds to osmotic pressure, and changes in membrane composition. The purpose of this study was to develop the relationship between membrane surface tension, and the allosteric balance of a mechanosensitive transporter such as NHE-1. In eukaryotes, the asymmetric composition of membrane leaflets results in a difference in surface tensions that is involved in the creation of a reservoir of intracellular vesicles and membrane buds contributing to buffer mechanical constraints. Therefore, we took this phenomenon into account in this study and developed a set of relations between the mean surface tension, membrane asymmetry, fluid phase endocytosis and the allosteric equilibrium constant of the transporter. We then used the experimental data published on the effects of osmotic pressure and membrane modification on the NHE-1 allosteric constant to fit these equations. We show here that NHE-1 mechanosensitivity is more based on its high sensitivity towards the asymmetry between the bilayer leaflets compared to mean global membrane tension. This compliance to membrane asymmetry is physiologically relevant as with their slower transport rates than ion channels, transporters cannot respond as high pressure-high conductance fast-gating emergency valves.

Electrostatics of the protein-water interface and the dynamical transition in proteins

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.

Protein dynamics and pressure: what can high pressure tell us about protein structural flexibility?

After a brief overview of NMR and X-ray crystallography studies on protein flexibility under pressure, we summarize the effects of hydrostatic pressure on the native fold of azurin from Pseudomonas aeruginosa as inferred from the variation of the intrinsic phosphorescence lifetime and the acrylamide bimolecular quenching rate constants of the buried Trp residue. The pressure/temperature response of the globular fold and modulation of its dynamical structure is analyzed both in terms of a reduction of internal cavities and of the hydration of the polypeptide. The study of the effect of two single point cavity forming mutations, F110S and I7S, on the unfolding volume change (ΔV(0)) of azurin and on the internal dynamics of the protein fold under pressure demonstrate that the creation of an internal cavity will enhance the plasticity and lower the stability of the globular structure. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Isothermal compressibility of macromolecular crystals and macromolecules derived from high-pressure X-ray crystallography

The compressibility of several nucleic acid and globular protein crystals has been investigated by high-pressure macromolecular crystallography. Further, crystal structures at four different pressures allowed the determination of the intrinsic compressibilityversuspressure of d(GGTATACC)2and hen egg-white lysozyme. For lysozyme, the values for the intrinsic molecular compressibility at atmospheric pressure and the nonlinearity index were 0.070 GPa−1and 8.15, respectively. On the basis of two crystal structures at atmospheric and high pressure, similar, albeit less complete, information was derived for d(CGCGAATTCGCG)2and bovine erythrocyte Cu,Zn superoxide dismutase. Using these data and accurate calculations of the solvent-excluded volume, the apparent solvent compressibility in the crystalline state was determined as a function of pressure and compared with results from a simple model that assumes invariant unit-cell content, with the conclusion that solvent compressibility was abnormal for three out of the five crystals investigated. Experimental results suggest that macromolecular crystals submitted to high pressure may have a variable unit-cell mass due to solvent exchange with the surrounding pool, as already observed in other hydrated crystals such as zeolites.

Multiscale mechanics of fibrin polymer: gel stretching with protein unfolding and loss of water

Blood clots and thrombi consist primarily of a mesh of branched fibers made of the protein fibrin. We propose a molecular basis for the marked extensibility and negative compressibility of fibrin gels based on the structural and mechanical properties of clots at the network, fiber, and molecular levels. The force required to stretch a clot initially rises linearly and is accompanied by a dramatic decrease in clot volume and a peak in compressibility. These macroscopic transitions are accompanied by fiber alignment and bundling after forced protein unfolding. Constitutive models are developed to integrate observations at spatial scales that span six orders of magnitude and indicate that gel extensibility and expulsion of water are both manifestations of protein unfolding, which is not apparent in other matrix proteins such as collagen.

Detection without deflection? A hypothesis for direct sensing of sound pressure by hair cells

It is widely thought that organisms detect sound by sensing the deflection of hair-like projections, the stereocilia, at the apex of hair cells. In the case of mammals, the standard interpretation is that hair cells in the cochlea respond to deflection of stereocilia induced by motion generated by a hydrodynamic travelling wave. But in the light of persistent anomalies, an alternative hypothesis seems to have some merit: that sensing cells (in particular the outer hair cells) may, at least at low intensities, be reacting to a different stimulus - the rapid pressure wave that sweeps through the cochlear fluids at the speed of sound in water. This would explain why fast responses are sometimes seen before the peak of the travelling wave. Yet how could cells directly sense fluid pressure? Here, a model is constructed of the outer hair cell as a pressure vessel able to sense pressure variations across its cuticular pore, and this 'fontanelle' model, based on the sensing action of the basal body at this compliant spot, could explain the observed anomalies. Moreover, the fontanelle model can be applied to a wide range of other organisms, suggesting that direct pressure detection is a general mode of sensing complementary to stereociliar displacement.