Effect of pressure on individual hydrogen bonds in proteins. Basic pancreatic trypsin inhibitor

Including the Ensemble of Unstructured Conformations in the Analysis of Protein's Native State by High-Pressure NMR Spectroscopy

The analysis of pressure induced changes in the chemical shift of proteins allows statements on structural fluctuations proteins exhibit at ambient pressure. The inherent issue of separating general pressure effects from structural related effects on the pressure dependence of chemical shifts has so far been addressed by considering the characteristics of random coil peptides on increasing pressure. In this work, chemically and pressure denatured states of the cold shock protein B from Bacillus subtilis (BsCspB) have been assigned in 2D 1H-15N HSQC NMR spectra and their dependence on increasing hydrostatic pressure has been evaluated. The pressure denatured polypeptide chain has been used to separate general from structural related effects on 1H and 15N chemical shifts of native BsCspB and the implications on the interpretation of pressure induced changes in the chemical shift regarding the structure of BsCspB are discussed. It has been found that the ensemble of unstructured conformations of BsCspB shows different responses to increasing pressure than random coil peptides do. Thus, the approach used for considering the general effects that arise when hydrostatic pressure increases changes the structural conclusions that are drawn from high pressure NMR spectroscopic experiments that rely on the analysis of chemical shifts.

Correlation between protein conformations and water structure and thermodynamics at high pressure: A molecular dynamics study of the Bovine Pancreatic Trypsin Inhibitor (BPTI) protein

Pressure-induced perturbation of a protein structure leading to its folding-unfolding mechanism is an important yet not fully understood phenomenon. The key point here is the role of water and its coupling with protein conformations as a function of pressure. In the current work, using extensive molecular dynamics simulation at 298 K, we systematically examine the coupling between protein conformations and water structures of pressures of 0.001, 5, 10, 15, 20 kbar, starting from (partially) unfolded structures of the protein Bovine Pancreatic Trypsin Inhibitor (BPTI). We also calculate localized thermodynamics at those pressures as a function of protein-water distance. Our findings show that both protein-specific and generic effects of pressure are operating. In particular, we found that (1) the amount of increase in water density near the protein depends on the protein structural heterogeneity; (2) the intra-protein hydrogen bond decreases with pressure, while the water-water hydrogen bond per water in the first solvation shell (FSS) increases; protein-water hydrogen bonds also found to increase with pressure, (3) with pressure hydrogen bonds of waters in the FSS getting twisted; and (4) water's tetrahedrality in the FSS decreases with pressure, but it is dependent on the local environment. Thermodynamically, at higher pressure, the structural perturbation of BPTI is due to the pressure-volume work, while the entropy decreases with the increase of pressure due to the higher translational and rotational rigidity of waters in the FSS. The local and subtle effects of pressure, found in this work, are likely to be typical of pressure-induced protein structure perturbation.

The residual structure of acid-denatured β2 -microglobulin is relevant to an ordered fibril morphology

β₂ -Microglobulin (β2m) forms amyloid fibrils in vitro under acidic conditions. Under these conditions, the residual structure of acid-denatured β2m is relevant to seeding and fibril extension processes. Disulfide (SS) bond-oxidized β2m has been shown to form rigid, ordered fibrils, whereas SS bond-reduced β2m forms curvy, less-ordered fibrils. These findings suggest that the presence of an SS bond affects the residual structure of the monomer, which subsequently influences the fibril morphology. To clarify this process, we herein performed NMR experiments. The results obtained revealed that oxidized β2m contained a residual structure throughout the molecule, including the N- and C-termini, whereas the residual structure of the reduced form was localized and other regions had a random coil structure. The range of the residual structure in the oxidized form was wider than that of the fibril core. These results indicate that acid-denatured β2m has variable conformations. Most conformations in the ensemble cannot participate in fibril formation because their core residues are hidden by residual structures. However, when hydrophobic residues are exposed, polypeptides competently form an ordered fibril. This conformational selection phase may be needed for the ordered assembly of amyloid fibrils.

Behavior of B- and Z-DNA Crystals under High Hydrostatic Pressure

Single crystals of B-DNA and Z-DNA oligomers were analyzed under high hydrostatic pressure and their behavior was compared to the A-DNA crystals already known. The amplitude of the base compression, when compared to the A-form of DNA (0.13 Å/GPa), was higher for the Z-DNA (0.32 Å/GPa) and was the highest for the B-DNA (0.42 Å/GPa). The B-DNA crystal degraded rapidly around 400–500 MPa, while the Z-structure was more resistant, up to 1.2 GPa.

Conformational Properties Relevant to the Amyloidogenicity of β2-Microglobulin Analyzed Using Pressure- and Salt-Dependent Chemical Shift Data

β₂-Microglobulin (β₂m) is associated with dialysis-related amyloidosis. In vitro experiments have shown that β₂m forms amyloid fibrils at acidic pHs in the presence of moderate concentrations of salt. Previous studies suggested that acid-denatured β₂m has a hydrophobic residual structure, and the exposure of the hydrophobic residues enhances the association with seeds or other β₂m monomers. However, the nature of the residual structure relevant to its amyloidogenicity remains to be investigated. To understand the structural properties of acid-denatured β₂m and the role of salt, we investigated pressure- and salt concentration-dependent conformational changes by nuclear magnetic resonance spectroscopy and other methods. Here, pressure was utilized to characterize the conformers existing in a conformational equilibrium at ambient pressure. The obtained pressure- and salt concentration-dependent chemical shift data were simultaneously subjected to principal component analysis to characterize individual conformational change events. Unexpectedly, the addition of salt induced an expansion of the β₂m molecule, which likely resulted from the exclusion of the N-terminal region from the hydrophobic cluster region. The dissected chemical shift patterns for the salt-induced conformational change and other experimental data indicated that this conformational change caused a rigidification in the intrinsic hydrophobic cluster, leading to the observed amyloidogenicity.

Entropic allostery dominates the phosphorylation-dependent regulation of Syk tyrosine kinase release from immunoreceptor tyrosine-based activation motifs

Spleen tyrosine kinase (Syk) is an essential player in immune signaling through its ability to couple multiple classes of membrane immunoreceptors to intracellular signaling pathways. Ligand binding leads to the recruitment of Syk to a phosphorylated cytoplasmic region of the receptors called ITAM. Syk binds to ITAM with high-affinity (nanomolar Kd ) via its tandem pair of SH2 domains. The affinity between Syk and ITAM is allosterically regulated by phosphorylation at Y130 in a linker connecting the tandem SH2 domains; when Y130 is phosphorylated, the binding affinity decreases (micromolar Kd ). Previous equilibrium binding studies attribute the increase in the binding free energy to an intra-molecular binding (isomerization) step of the tandem SH2 and ITAM, but a physical basis for the increased free energy is unknown. Here, we provide evidence that Y130 phosphorylation imposes an entropy penalty to isomerization, but surprisingly, has negligible effect on the SH2 binding interactions with ITAM and thus on the binding enthalpy. An analysis of NMR chemical shift differences characterized conformational effects of ITAM binding, and binding thermodynamics were measured from isothermal titration calorimetry. Together the data support a previously unknown mechanism for the basis of regulating protein-protein interactions through protein phosphorylation. The decreased affinity for Syk association with immune receptor ITAMs by Y130 phosphorylation is an allosteric mechanism driven by an increased entropy penalty, likely contributed by conformational disorder in the SH2-SH2 inter-domain structure, while SH2-ITAM binding contacts are not affected, and binding enthalpy is unchanged.

Hydrogen bonds in the vicinity of the special pair of the bacterial reaction center probed by hydrostatic high-pressure absorption spectroscopy

Using the native bacteriochlorophyll a pigment cofactors as local probes, we investigated the response to external hydrostatic high pressure of reaction center membrane protein complexes from the photosynthetic bacterium Rhodobacter sphaeroides. Wild-type and engineered complexes were used with a varied number (0, 1 or 2) of hydrogen bonds that bind the reaction center primary donor bacteriochlorophyll cofactors to the surrounding protein scaffold. A pressure-induced breakage of hydrogen bonds was established for both detergent-purified and membrane-embedded reaction centers, but at rather different pressures: between 0.2 and 0.3GPa and at about 0.55GPa, respectively. The free energy change associated with the rupture of the single hydrogen bond present in wild-type reaction centers was estimated to be equal to 13-14kJ/mol. In the mutant with two symmetrical hydrogen bonds (FM197H) a single cooperative rupture of the two bonds was observed corresponding to an about twice stronger bond, rather than a sequential rupture of two individual bonds.

Pressure dependence of side chain 13C chemical shifts in model peptides Ac-Gly-Gly-Xxx-Ala-NH2

For evaluating the pressure responses of folded as well as intrinsically unfolded proteins detectable by NMR spectroscopy the availability of data from well-defined model systems is indispensable. In this work we report the pressure dependence of ¹³C chemical shifts of the side chain atoms in the protected tetrapeptides Ac-Gly-Gly-Xxx-Ala-NH₂ (Xxx, one of the 20 canonical amino acids). Contrary to expectation the chemical shifts of a number of nuclei have a nonlinear dependence on pressure in the range from 0.1 to 200 MPa. The size of the polynomial pressure coefficients B ₁ and B ₂ is dependent on the type of atom and amino acid studied. For H^N, N and C^α the first order pressure coefficient B ₁ is also correlated to the chemical shift at atmospheric pressure. The first and second order pressure coefficients of a given type of carbon atom show significant linear correlations suggesting that the NMR observable pressure effects in the different amino acids have at least partly the same physical cause. In line with this observation the magnitude of the second order coefficients of nuclei being direct neighbors in the chemical structure also are weakly correlated. The downfield shifts of the methyl resonances suggest that gauche conformers of the side chains are not preferred with pressure. The valine and leucine methyl groups in the model peptides were assigned using stereospecifically ¹³C enriched amino acids with the pro-R carbons downfield shifted relative to the pro-S carbons.

Pressure dependence of backbone chemical shifts in the model peptides Ac-Gly-Gly-Xxx-Ala-NH2

For a better understanding of nuclear magnetic resonance (NMR) detected pressure responses of folded as well as unstructured proteins the availability of data from well-defined model systems are indispensable. In this work we report the pressure dependence of chemical shifts of the backbone atoms (1)H(α), (13)C(α) and (13)C' in the protected tetrapeptides Ac-Gly-Gly-Xxx-Ala-NH2 (Xxx one of the 20 canonical amino acids). Contrary to expectation the chemical shifts of these nuclei have a nonlinear dependence on pressure in the range from 0.1 to 200 MPa. The polynomial pressure coefficients B 1 and B 2 are dependent on the type of amino acid studied. The coefficients of a given nucleus show significant linear correlations suggesting that the NMR observable pressure effects in the different amino acids have at least partly the same physical cause. In line with this observation the magnitude of the second order coefficients of nuclei being direct neighbors in the chemical structure are also weakly correlated.

Cavities and Excited States in Proteins

Protein cavities or voids are observed as defects in atomic packing. Cavities have long been suggested to play important roles in protein dynamics and function, but the underlying origin and mechanism remains elusive. Here, recent studies about the cavities characterized by high-pressure NMR spectroscopy have been reviewed. Analysis of the pressure-dependent chemical shifts showed both linear and nonlinear response of proteins to pressure. The linear response corresponded to compression within the native ensemble, while the nonlinear response indicated the involvement of low-lying excited states that were different from the native state. The finding of non-linear pressure shifts in various proteins suggested that the existence of the low-lying excited states was common for globular proteins. However, the absolute nonlinear coefficient values varied significantly from protein to protein, and showed a good correlation with the density of cavities. Extensive studies on hen lysozyme as a model system showed that the cavity hydration and water penetration into the interior of proteins was an origin of the conformational transition to the excited states. The importance of cavities for protein function and evolution has also been explained. In addition to these "equilibrium" cavities, there are also "transient" cavities formed in the interior of the protein structure, as manifested by the ring flip motions of aromatic rings. The significance of transient cavities, reflecting an intrinsic dynamic nature within the native state, has also been discussed.

Functional Sub-states by High-pressure Macromolecular Crystallography

At the molecular level, high-pressure perturbation is of particular interest for biological studies as it allows trapping conformational substates. Moreover, within the context of high-pressure adaptation of deep-sea organisms, it allows to decipher the molecular determinants of piezophily. To provide an accurate description of structural changes produced by pressure in a macromolecular system, developments have been made to adapt macromolecular crystallography to high-pressure studies. The present chapter is an overview of results obtained so far using high-pressure macromolecular techniques, from nucleic acids to virus capsid through monomeric as well as multimeric proteins.

A hypothesis to reconcile the physical and chemical unfolding of proteins

High pressure (HP) or urea is commonly used to disturb folding species. Pressure favors the reversible unfolding of proteins by causing changes in the volumetric properties of the protein-solvent system. However, no mechanistic model has fully elucidated the effects of urea on structure unfolding, even though protein-urea interactions are considered to be crucial. Here, we provide NMR spectroscopy and 3D reconstructions from X-ray scattering to develop the "push-and-pull" hypothesis, which helps to explain the initial mechanism of chemical unfolding in light of the physical events triggered by HP. In studying MpNep2 from Moniliophthora perniciosa, we tracked two cooperative units using HP-NMR as MpNep2 moved uphill in the energy landscape; this process contrasts with the overall structural unfolding that occurs upon reaching a threshold concentration of urea. At subdenaturing concentrations of urea, we were able to trap a state in which urea is preferentially bound to the protein (as determined by NMR intensities and chemical shifts); this state is still folded and not additionally exposed to solvent [fluorescence and small-angle X-ray scattering (SAXS)]. This state has a higher susceptibility to pressure denaturation (lower p1/2 and larger ΔVu); thus, urea and HP share concomitant effects of urea binding and pulling and water-inducing pushing, respectively. These observations explain the differences between the molecular mechanisms that control the physical and chemical unfolding of proteins, thus opening up new possibilities for the study of protein folding and providing an interpretation of the nature of cooperativity in the folding and unfolding processes.

High Pressure NMR Spectroscopy

The combination of high-resolution NMR spectroscopy with pressure perturbation, known as variable-pressure NMR spectroscopy or simply high pressure NMR spectroscopy, is a relatively recent accomplishment, but is a technique expanding rapidly with high promise in future. The importance of the method is that it allows, for the first time in history, a systematic means of detecting and analyzing the structures and thermodynamic stability of high-energy sub-states in proteins. High-energy sub-states have been only vaguely known so far, as normally their populations are too low to be detected by conventional spectroscopic techniques including NMR spectroscopy. By now, however, high pressure NMR spectroscopy has established unequivocally that high-energy conformers are universally present in proteins in equilibrium with their stable folded counterparts. This chapter describes briefly the techniques of high pressure NMR spectroscopy and its unique and novel aspects as a method to explore protein structure in its high-energy paradigm with illustrative examples. It is now well established that high pressure NMR spectroscopy is a method to study intrinsic fluctuations of proteins, rather than those forced by pressure, by detecting structural changes amplified by pressure. Extension of the method to other bio-macromolecular systems is considered fairly straightforward.

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.

Effect of pressure on thermal stability of g-quadruplex DNA and double-stranded DNA structures

Pressure is a thermodynamic parameter that can induce structural changes in biomolecules due to a volumetric decrease. Although most proteins are denatured by pressure over 100 MPa because they have the large cavities inside their structures, the double-stranded structure of DNA is stabilized or destabilized only marginally depending on the sequence and salt conditions. The thermal stability of the G-quadruplex DNA structure, an important non-canonical structure that likely impacts gene expression in cells, remarkably decreases with increasing pressure. Volumetric analysis revealed that human telomeric DNA changed by more than 50 cm3 mol-1 during the transition from a random coil to a quadruplex form. This value is approximately ten times larger than that for duplex DNA under similar conditions. The volumetric analysis also suggested that the formation of G-quadruplex DNA involves significant hydration changes. The presence of a cosolute such as poly(ethylene glycol) largely repressed the pressure effect on the stability of G-quadruplex due to alteration in stabilities of the interactions with hydrating water. This review discusses the importance of local perturbations of pressure on DNA structures involved in regulation of gene expression and highlights the potential for application of high-pressure chemistry in nucleic acid-based nanotechnology.

Impact of hydrostatic pressure on an intrinsically disordered protein: a high-pressure NMR study of α-synuclein

The impact of pressure on the backbone (15) N, (1) H and (13) C chemical shifts in N-terminally acetylated α-synuclein has been evaluated over a pressure range 1-2500 bar. Even while the chemical shifts fall very close to random coil values, as expected for an intrinsically disordered protein, substantial deviations in the pressure dependence of the chemical shifts are seen relative to those in short model peptides. In particular, the nonlinear pressure response of the (1) H(N) chemical shifts, which commonly is associated with the presence of low-lying "excited states", is much larger in α-synuclein than in model peptides. The linear pressure response of (1) H(N) chemical shift, commonly linked to H-bond length change, correlates well with those in short model peptides, and is found to be anticorrelated with its temperature dependence. The pressure dependence of (13) C chemical shifts shows remarkably large variations, even when accounting for residue type, and do not point to a clear shift in population between different regions of the Ramachandran map. However, a nearly universal decrease in (3) JHN-Hα by 0.22 ± 0.05 Hz suggests a slight increase in population of the polyproline II region at 2500 bar. The first six residues of N-terminally acetylated synuclein show a transient of approximately 15% population of α-helix, which slightly diminishes at 2500 bar. The backbone dynamics of the protein is not visibly affected beyond the effect of slight increase in water viscosity at 2500 bar.

Volumetric properties underlying ligand binding in a monomeric hemoglobin: a high-pressure NMR study

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using (1)H-(15)N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Pressure Dependence of 15N Chemical Shifts in Model Peptides Ac-Gly-Gly-X-Ala-NH2

High pressure NMR spectroscopy has developed into an important tool for studying conformational equilibria of proteins in solution. We have studied the amide proton and nitrogen chemical shifts of the 20 canonical amino acids X in the random-coil model peptide Ac-Gly-Gly-X-Ala-NH₂, in a pressure range from 0.1 to 200 MPa, at a proton resonance frequency of 800 MHz. The obtained data allowed the determination of first and second order pressure coefficients with high accuracy at 283 K and pH 6.7. The mean first and second order pressure coefficients <B₁^15N> and <B₂^15N> for nitrogen are 2.91 ppm/GPa and −2.32 ppm/GPa², respectively. The corresponding values <B₁^1H> and <B₂^1H> for the amide protons are 0.52 ppm/GPa and −0.41 ppm/GPa². Residual dependent ¹J_1H15N-coupling constants are shown.

High-pressure macromolecular crystallography and NMR: status, achievements and prospects

Biomacromolecules are thermodynamic entities that exist in general as an equilibrium mixture of the basic folded state and various higher-energy substates including all functionally relevant ones. Under physiological conditions, however, the higher-energy substates are usually undetectable on spectroscopy, as their equilibrium populations are extremely low. Hydrostatic pressure gives a general solution to this problem. As proteins generally have smaller partial molar volumes in higher-energy states than in the basic folded state, pressure can shift the equilibrium toward the former substantially, and allows their direct detection and analysis with X-ray crystallography or NMR spectroscopy at elevated pressures. These techniques are now mature, and their status and selected applications are presented with future prospects.

The stability of almond β-glucosidase during combined high pressure-thermal processing: a kinetic study

The thermal and the combined high pressure-thermal inactivation kinetics of almond β-glucosidase (β-D-glucoside glucohydrolase, EC 3.2.1.21) were investigated at pressures from 0.1 to 600 MPa and temperatures ranging from 30 to 80 °C. Thermal treatments at temperatures higher than 50 °C resulted in significant inactivation with complete inactivation after 2 min of treatment at 80 °C. Both the thermal and high pressure inactivation kinetics were described well by first-order model. Application of pressure increased the inactivation kinetics of the enzyme except at moderate temperatures (50 to 70 °C) and pressures between 0.1 and 100 MPa where slight pressure stabilisation of the enzyme against thermal denaturation was observed. The activation energy for the inactivation of the enzyme at atmospheric pressure was estimated to be 216.2±8.6 kJ/mol decreasing to 55.2±3.9 kJ/mol at 600 MPa. The activation volumes were negative at all temperature conditions excluding the temperature-pressure range where slight pressure stabilisation was observed. The values of the activation volumes were estimated to be -29.6±0.6, -29.8±1.7, -20.6±3.2, -41.2±4.8, -36.5±1.8, -39.6±4.3, -31.0±4.5 and -33.8±3.9 cm3/mol at 30, 35, 40, 45, 50, 60, 65 and 70 °C, respectively, with no clear trend with temperature. The pressure-temperature dependence of the inactivation rate constants was well described by an empirical third-order polynomial model.

A delicate interplay of structure, dynamics, and thermodynamics for function: a high pressure NMR study of outer surface protein A

Outer surface protein A (OspA) is a crucial protein in the infection of Borrelia burgdorferi causing Lyme disease. We studied conformational fluctuations of OspA with high-pressure (15)N/(1)H two-dimensional NMR along with high-pressure fluorescence spectroscopy. We found evidence within folded, native OspA for rapid local fluctuations of the polypeptide backbone in the nonglobular single layer β-sheet connecting the N- and C-terminal domains with τ << ms, which may give the two domains certain independence in mobility and thermodynamic stability. Furthermore, we found that folded, native OspA is in equilibrium (τ >> ms) with a minor conformer I, which is almost fully disordered and hydrated for the entire C-terminal part of the polypeptide chain from β8 to the C-terminus. Conformer I is characterized with ΔG(0) = 32 ± 9 kJ/mol and ΔV(0) = -140 ± 40 mL/mol, populating only ∼0.001% at 40°C at 0.1 MPa, pH 5.9. Because in the folded conformer the receptor binding epitope of OspA is buried in the C-terminal domain, its transition into conformer I under in vivo conditions may be critical for the infection of B. burgdorferi. The formation and stability of the peculiar conformer I are apparently supported by a large packing defect or cavity located in the C-terminal domain.

STUDY ON COMPLEXES OF TRYPSIN AND ITS INHIBITORS BY MEANS OF ATOM-BOND ELECTRONEGATIVITY EQUALIZATION METHOD FUSED INTO MOLECULAR MECHANICS (ABEEM/MM)

Trypsin is one of the most important enzymes and plays important roles in the regulation of biological processes. The newly developed trypsin–inhibitor interaction potential model based on atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) was employed to study complexes of trypsin and its inhibitors. Some structural properties, including root-mean-square deviations (RMSD) of bond length, bond angle and key dihedral, and coordinated RMS Shifts of atoms and hydrogen bond, were studied using ABEEM/MM method and compared with OPLS-AA force-field. At the same time, comparative study on the charges of the hydrogen atoms at the tail of the ligand was also discussed in order to investigate the effect of hydrogen bond between trypsin and its ligand. This work demonstrates that ABEEM/MM model can well reproduce good structures for these trypsin–inhibitor complexes with rather small RMSD of bond length, bond angles, key dihedrals, and RMS Shifts of atomic coordinates with respect to the experimental crystal structures, compared with the results from the OPLS-AA method. The charges obtained by ABEEM/MM model are in good accordance with those from an ab initio calculation. Moreover, both the polarization and the salt-bridge effects have been taken into account. It is shown that ABEEM charges can properly describe the electrostatic interactions between the protein trypsin and its inhibitors.

NMR structures of apo L. casei dihydrofolate reductase and its complexes with trimethoprim and NADPH: contributions to positive cooperative binding from ligand-induced refolding, conformational changes, and interligand hydrophobic interactions

In order to examine the origins of the large positive cooperativity (ΔG(0)(coop) = -2.9 kcal mol(-1)) of trimethoprim (TMP) binding to a bacterial dihydrofolate reductase (DHFR) in the presence of NADPH, we have determined and compared NMR solution structures of L. casei apo DHFR and its binary and ternary complexes with TMP and NADPH and made complementary thermodynamic measurements. The DHFR structures are generally very similar except for the A-B loop region and part of helix B (residues 15-31) which could not be directly detected for L. casei apo DHFR because of line broadening from exchange between folded and unfolded forms. Thermodynamic and NMR measurements suggested that a significant contribution to the cooperativity comes from refolding of apo DHFR on binding the first ligand (up to -0.95 kcals mol(-1) if 80% of A-B loop requires refolding). Comparisons of Cα-Cα distance differences and domain rotation angles between apo DHFR and its complexes indicated that generally similar conformational changes involving domain movements accompany formation of the binary complexes with either TMP or NADPH and that the binary structures are approaching that of the ternary complex as would be expected for positive cooperativity. These favorable ligand-induced structural changes upon binding the first ligand will also contribute significantly to the cooperative binding. A further substantial contribution to cooperative binding results from the proximity of the bound ligands in the ternary complex: this reduces the solvent accessible area of the ligand and provides a favorable entropic hydrophobic contribution (up to -1.4 kcal mol(-1)).

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.

Application of high pressure laser flash photolysis in studies on selected hemoprotein reactions

This article focuses on the application of high pressure laser flash photolysis for studies on selected hemoprotein reactions with the objective to establish details of the underlying reaction mechanisms. In this context, particular attention is given to the reactions of small molecules such as dioxygen, carbon monoxide, and nitric oxide with selected hemoproteins (hemoglobin, myoglobin, neuroglobin and cytochrome P450(cam)), as well as to photo-induced electron transfer reactions occurring in hemoproteins (particularly in various types of cytochromes). Mechanistic conclusions based on the interpretation of the obtained activation volumes are discussed in this account.

Combined chemical shift changes and amino acid specific chemical shift mapping of protein-protein interactions

Protein-protein interactions are often studied by chemical shift mapping using solution NMR spectroscopy. When heteronuclear data are available the interaction interface is usually predicted by combining the chemical shift changes of different nuclei to a single quantity, the combined chemical shift perturbation Deltadelta comb In this paper different procedures (published and non-published) to calculate Deltadelta comb are examined that include a variety of different functional forms and weighting factors for each nucleus. The predictive power of all shift mapping methods depends on the magnitude of the overlap of the chemical shift distributions of interacting and non-interacting residues and the cut-off criterion used. In general, the quality of the prediction on the basis of chemical shift changes alone is rather unsatisfactory but the combination of chemical shift changes on the basis of the Hamming or the Euclidian distance can improve the result. The corrected standard deviation to zero of the combined chemical shift changes can provide a reasonable cut-off criterion. As we show combined chemical shifts can also be applied for a more reliable quantitative evaluation of titration data.

Complexity of aromatic ring-flip motions in proteins: Y97 ring dynamics in cytochrome c observed by cross-relaxation suppressed exchange NMR spectroscopy

Dynamics of large-amplitude conformational motions in proteins are complex and less understood, although these processes are intimately associated with structure, folding, stability, and function of proteins. Here, we use a large set of spectra obtained by cross-relaxation suppressed exchange NMR spectroscopy (EXSY) to study the 180 degrees flipping motion of the Y97 ring of horse ferricytochrome c as a function of near-physiological temperature in the 288-308 K range. With rising temperature, the ring-flip rate constant makes a continuous transition from Arrhenius to anti-Arrhenius behavior through a narrow Arrhenius-like zone. This behavior is seen not only for the native state of the protein, but also for native-like states generated by adding subdenaturing amounts of guanidine deuterochloride (GdnDCl). Moderately destabilizing concentrations of the denaturant (1.5 M GdnDCl) completely removes the Arrhenius-like feature from the temperature window employed. The Arrhenius to anti-Arrhenius transition can be explained by the heat capacity model where temperature strengthens ground state interactions, perhaps hydrophobic in nature. The effect of the denaturant may appear to arise from direct protein-denaturant interactions that are structure-stabilizing under subdenaturing conditions. The temperature distribution of rate constants under different stability conditions also suggests that the prefactor in Arrhenius-like relations is temperature dependent. Although the use of the transition state theory (TST) offers several challenges associated with data interpretation, the present results and a consideration of others published earlier provide evidence for complexity of ring-flip dynamics in proteins.

Adaptation of the base-paired double-helix molecular architecture to extreme pressure

The behaviour of the d(GGTATACC) oligonucleotide has been investigated by X-ray crystallography at 295 K in the range from ambient pressure to 2 GPa (approximately 20,000 atm). Four 3D-structures of the A-DNA form (at ambient pressure, 0.55, 1.09 and 1.39 GPa) were refined at 1.60 or 1.65 A resolution. In addition to the diffraction pattern of the A-form, the broad meridional streaks previously explained by occluded B-DNA octamers within the channels of the crystalline A-form matrix were observed up to at least 2 GPa. This work highlights an important property of nucleic acids, their capability to withstand very high pressures, while keeping in such conditions a nearly invariant geometry of base pairs that store and carry genetic information. The double-helix base-paired architecture behaves as a molecular spring, which makes it especially adapted to very harsh conditions. These features may have contributed to the emergence of a RNA World at prebiotic stage.

Species-specific differences in the intermediate states of human and Syrian hamster prion protein detected by high pressure NMR spectroscopy

Human (huPrP) and Syrian hamster (ShaPrP) prion proteins have barriers for mutual infectivity, although they fold into almost an identical structure. The pressure responses of huPrP and ShaPrP characterized by high pressure NMR spectroscopy show differences in their excited states, as monitored by pressure-induced chemical shifts and intensity changes of individual residues in the (15)N/(1)H HSQC spectra. Both proteins fluctuate rapidly between two well folded (native) conformations N(1) and N(2) and less frequently between N and the excited states I(1) and I(2) with local disorder that may present structural intermediates on the way to PrP(Sc). These four structural states can be observed in the hamster and human PrP. At ambient pressure, less than 5 molecules of 10,000 are in the intermediate state I(2). From the structural point of view, the different states are mutually different, particularly in positions strategically important for generating species barriers for infection. The results point to the notion that excited state conformers are important for infection and that their structural differences may crucially determine species barriers for infection.

Experiments with proteins at low temperature: What do we learn on properties in their functional state?

The authors compared the spectral response of Zn-substituted horseradish peroxidase in a glycerol/water solvent to hydrostatic pressure at 2 K and ambient temperature. The low temperature experiments clearly demonstrate the presence of at least three different conformations with drastically different elastic properties. However, the main conformation, which determines the fluorescence spectrum at ambient temperature, did not show any significant difference between low and high temperature and pressure. The authors conclude that the local compressibility of the heme pocket of the protein depends only very weakly on temperature.

Protein under pressure: Molecular dynamics simulation of the arc repressor

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.

Conformational fluctuations of proteins revealed by variable pressure NMR

With the high-resolution variable-pressure NMR spectroscopy, one can study conformational fluctuations of proteins in a much wider conformational space than hitherto explored by NMR and other spectroscopic techniques. This is because a protein in solution generally exists as a dynamic mixture of conformers mutually differing in partial molar volume, and pressure can select the population of a conformer according to its relative volume. In this review, we describe how variable-pressure NMR can be used to probe conformational fluctuations of proteins in a wide conformational space from the folded to the fully unfolded structures, with actual examples. Furthermore, the newly emerging technique "NMR snapshots" expresses amply fluctuating protein structures as changes in atomic coordinates. Finally, the concept of conformational fluctuation is extended to include intermolecular association leading to amyloidosis.

Effect of pressure on the conformation of proteins. A molecular dynamics simulation of lysozyme

The effect of pressure on the structure and mobility of lysozyme was studied by molecular dynamics computer simulation at 1 and 3 kbar (1 atm = 1.01325 bar = 101.325 kPa). The results have good agreement with the available experimental data, allowing the analysis of other features of the effect of pressure on the protein solution. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. From the analysis of secondary structure along the trajectories it is observed that the conformation under pressure is more stable, suggesting that pressure acts as a 'conformer selector' on the protein. The difference in solvent-accessed surface (SAS) with pressure shows a clear inversion of the hydrophilic/hydrophobic SAS ratio, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.

Insights into the role of hydration in protein structure and stability obtained through hydrostatic pressure studies

A thorough understanding of protein structure and stability requires that we elucidate the molecular basis for the effects of both temperature and pressure on protein conformational transitions. While temperature effects are relatively well understood and the change in heat capacity upon unfolding has been reasonably well parameterized, the state of understanding of pressure effects is much less advanced. Ultimately, a quantitative parameterization of the volume changes (at the basis of pressure effects) accompanying protein conformational transitions will be required. The present report introduces a qualitative hypothesis based on available model compound data for the molecular basis of volume change upon protein unfolding and its dependence on temperature.

Infrequent cavity-forming fluctuations in HPr from Staphylococcus carnosus revealed by pressure- and temperature-dependent tyrosine ring flips

Infrequent structural fluctuations of a globular protein is seldom detected and studied in detail. One tyrosine ring of HPr from Staphylococcus carnosus, an 88-residue phosphocarrier protein with no disulfide bonds, undergoes a very slow ring flip, the pressure and temperature dependence of which is studied in detail using the on-line cell high-pressure nuclear magnetic resonance technique in the pressure range from 3 MPa to 200 MPa and in the temperature range from 257 K to 313 K. The ring of Tyr6 is buried sandwiched between a beta-sheet and alpha-helices (the water-accessible area is less than 0.26 nm2), its hydroxyl proton being involved in an internal hydrogen bond. The ring flip rates 10(1)-10(5) s(-1) were determined from the line shape analysis of H(delta1, delta2) and H(epsilon1,epsilon2) of Tyr6, giving an activation volume DeltaV++ of 0.044 +/- 0.008 nm3 (27 mL mol(-1)), an activation enthalpy DeltaH++ of 89 +/- 10 kJ mol(-1), and an activation entropy DeltaS++ of 16 +/- 2 JK(-1) mol(-1). The DeltaV++) and DeltaH++ values for HPr found previously for Tyr and Phe ring flips of BPTI and cytochrome c fall within the range of DeltaV(double dagger) of 28 to 51 mL mol(-1) and DeltaH++ of 71 to 155 kJ mol(-1). The fairly common DeltaV++ and DeltaH++ values are considered to represent the extra space or cavity required for the ring flip and the extra energy required to create a cavity, respectively, in the core part of a globular protein. Nearly complete cold denaturation was found to take place at 200 MPa and 257 K independently from the ring reorientation process.