Bertrand Garcia-Moreno

Experimental measurement of the effective dielectric in the hydrophobic core of a protein

The dielectric inside a protein is a key physical determinant of the magnitude of electrostatic interactions in proteins. We have measured this dielectric phenomenologically, in terms of the dielectric that needs to be used with the Born equation in order to reproduce the observed pKa shifts induced by burial of an ionizable group in the hydrophobic core of a protein. Mutants of staphylococcal nuclease with a buried lysine residue at position 66 were engineered for this purpose. The pKa values of buried lysines were measured by difference potentiometry. The extent of coupling between the pKa and the global stability of the protein was evaluated by measuring pKa values in hyperstable forms of nuclease engineered to be 3.3 or 6.5 kcal mol-1 more stable than the wild type. The crystallographic structure of one mutant was determined to describe the environment of the buried lysine. The dielectrics that were measured range from 10 to 12. Published pKa values of buried ionizable residues in other proteins were analyzed in a similar fashion and the dielectrics obtained from these values are consistent with the ones measured in nuclease. These results argue strongly against the prevalent use of dielectrics of 4 or lower to describe the dielectric effect inside a protein in structure-based calculations of electrostatic energies with continuum dielectric models.

Molecular determinants of the pK(a) values of Asp and Glu residues in staphylococcal nuclease.

Prior computational studies of the acid‐unfolding behavior of staphylococcal nuclease (SNase) suggest that the pKa values of its carboxylic groups are difficult to reproduce with electrostatics calculations with continuum methods. To examine the molecular determinants of the pKa values of carboxylic groups in SNase, the pKa values of all 20 Asp and Glu residues were measured with multidimensional and multinuclear NMR spectroscopy in an acid insensitive variant of SNase. The crystal structure of the protein was obtained to describe the microenvironments of the carboxylic groups. Fourteen Asp and Glu residues titrate with relatively normal pKa values that are depressed by less than 1.1 units relative to the normal pKa of Asp and Glu in water. Only six residues have pKa values shifted by more than 1.5 units. Asp‐21 has an unusually high pKa of 6.5, which is probably the result of interactions with other carboxylic groups at the active site. The most perturbed pKa values appear to be governed by hydrogen bonding and not by Coulomb interactions. The pKa values calculated with standard continuum electrostatics methods applied to static structures are more depressed than the measured values because Coulomb effects are exaggerated in the calculations. The problems persist even when the protein is treated with the dielectric constant of water. This can be interpreted to imply that structural relaxation is an important determinant of the pKa values; however, no major pH‐sensitive conformational reorganization of the backbone was detected using NMR spectroscopy. Proteins 2009.

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10‐ns‐long molecular dynamics (MD) simulations was performed with wild‐type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean‐square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior. Proteins 2005.

Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease.

His121 and His124 are embedded in a network of polar and ionizable groups on the surface of staphylococcal nuclease. To examine how membership in a network affects the electrostatic properties of ionizable groups, the tautomeric state and the pK(a) values of these histidines were measured with NMR spectroscopy in the wild-type nuclease and in 13 variants designed to disrupt the network. In the background protein, His121 and His124 titrate with pK(a) values of 5.2 and 5.6, respectively. In the variants, where the network was disrupted, the pK(a) values range from 4.03 to 6.46 for His121, and 5.04 to 5.99 for His124. The largest decrease in a pK(a) was observed when the favorable Coulomb interaction between His121 and Glu75 was eliminated; the largest increase was observed when Tyr91 or Tyr93 was substituted with Ala or Phe. In all variants, the dominant tautomeric state at neutral pH was the N(epsilon2) state. At one level the network behaves as a rigid unit that does not readily reorganize when disrupted: crystal structures of the E75A or E75Q variants show that even when the pivotal Glu75 is removed, the overall configuration of the network was unaffected. On the other hand, a few key hydrogen bonds appear to govern the conformation of the network, and when these bonds are disrupted the network reorganizes. Coulomb interactions within the network report an effective dielectric constant of 20, whereas a dielectric constant of 80 is more consistent with the magnitude of medium to long-range Coulomb interactions in this protein. The data demonstrate that when structures are treated as static, rigid bodies, structure-based pK(a) calculations with continuum electrostatics method are not useful to treat ionizable groups in cases where pK(a) values are governed by short-range polar and Coulomb interactions.

Towards a Quantitative Understanding of Protein Hydration and Volumetric Properties

Herein, we probe by pressure perturbation calorimetry (PPC) the coefficient of thermal expansion, the volumetric and the hydration properties of variants of a hyperstable variant of staphylococcal nuclease (SNase), Delta+PHS. The temperature-dependent volumetric properties of the folded and unfolded states of the wild-type protein are calculated with previously published data. The present PPC results are used to interpret the volume diagram and expansivity at a molecular level. We conclude that the expansivity of the unfolded state is, to a first approximation, temperature independent, while that of the folded state decreases with increasing temperature. Our data suggest that at low temperature the defining contribution to DeltaV comes mainly from excluded volume differences and DeltaV for unfolding is negative. In contrast, at high temperatures, differential solvation due to the increased exposed surface area of the unfolded state and, in particular, its larger thermal volume linked to the increased conformational dynamics of the unfolded state ensemble takes over and DeltaV for unfolding eventually becomes positive.

Adaptations of proteins to cellular and subcellular pH

Bioinformatics-based searches for correlations between subcellular localization and pI or charge distribution of proteins have failed to detect meaningful correlations. Recent work published in BMC Biology finds that a physicochemical metric of charge distribution correlates better with subcellular pH than does pI. See research article http://www.biomedcentral.com/1741-7007/7/69

Effect of Internal Cavities on Folding Rates and Routes Revealed by Real-Time Pressure-Jump NMR Spectroscopy

The time required to fold proteins usually increases significantly under conditions of high pressure. Taking advantage of this general property of proteins, we combined P-jump experiments with NMR spectroscopy to examine in detail the folding reaction of staphylococcal nuclease (SNase) and of some of its cavity-containing variants. The nearly 100 observables that could be measured simultaneously collectively describe the kinetics of folding as a function of pressure and denaturant concentration with exquisite site-specific resolution. SNase variants with cavities in the central core of the protein exhibit a highly heterogeneous transition-state ensemble (TSE) with a smaller solvent-excluded void volume than the TSE of the parent SNase. This heterogeneous TSE experiences Hammond behavior, becoming more native-like (higher molar volume) with increasing denaturant concentration. In contrast, the TSE of the L125A variant, which has a cavity at the secondary core, is only slightly different from that of the parent SNase. Because pressure acts mainly to eliminate solvent-excluded voids, which are heterogeneously distributed throughout structures, it perturbs the protein more selectively than chemical denaturants, thereby facilitating the characterization of intermediates and the consequences of packing on folding mechanisms. Besides demonstrating how internal cavities can affect the routes and rates of folding of a protein, this study illustrates how the combination of P-jump and NMR spectroscopy can yield detailed mechanistic insight into protein folding reactions with exquisite site-specific temporal information.

Direct Evidence for Deprotonation of a Lysine Side Chain Buried in the Hydrophobic Core of a Protein

We report direct evidence for deprotonation of a lysine side chain buried in the hydrophobic core of a protein, demonstrating heteronuclear 1H-15N NMR data on the Lys-66 side chain amine (Nzeta) group in the delta-PHS/V66K variant of staphylococcal nuclease. Previous crystallographic study has shown that the Lys-66 Nzeta group is completely buried in the hydrophobic core. On the basis of double and triple resonance experiments, we found that the 1Hzeta and 15Nzeta chemical shifts at pH 8.0 and 6 degrees C for the buried lysine are 0.81 and 23.3 ppm, respectively, which are too abnormal to correspond to the protonated (NH3+) state. Further investigations using a model system suggested that the abnormal 1H and 15N chemical shifts represent the deprotonated (NH2) state of the Lys-66 Nzeta group. More straightforward evidence for the deprotonation was obtained with 2D F1-1H-coupled 1H-15N heteronuclear correlation experiments. Observed 15N multiplets clearly indicated that the spin system for the Lys-66 Nzeta group is AX2 (NH2) rather than AX3 (NH3+). Interestingly, although the amine group is buried in the hydrophobic core, the hydrogen exchange between water and the Lys-66 Nzeta group was found to be relatively rapid (93 s(-1) at -1 degrees C), which suggests the presence of a dynamic process such as local unfolding or water penetration. The partial self-decoupling effect on 15Nzeta multiplets due to the rapid hydrogen exchange is also discussed.

Remodeling of the Folding Free Energy Landscape of Staphylococcal Nuclease by Cavity-Creating Mutations

The folding of staphylococcal nuclease (SNase) is known to proceed via a major intermediate in which the central OB subdomain is folded and the C-terminal helical subdomain is disordered. To identify the structural and energetic determinants of this folding free energy landscape, we have examined in detail, using high-pressure NMR, the consequences of cavity creating mutations in each of the two subdomains of an ultrastable SNase, Δ+PHS. The stabilizing mutations of Δ+PHS enhanced the population of the major folding intermediate. Cavity creation in two different regions of the Δ+PHS reference protein, despite equivalent effects on global stability, had very distinct consequences on the complexity of the folding free energy landscape. The L125A substitution in the C-terminal helix of Δ+PHS slightly suppressed the major intermediate and promoted an additional excited state involving disorder in the N-terminus, but otherwise decreased landscape heterogeneity with respect to the Δ+PHS background protein. The I92A substitution, located in the hydrophobic OB-fold core, had a much more profound effect, resulting in a significant increase in the number of intermediate states and implicating the entire protein structure. Denaturant (GuHCl) had very subtle and specific effects on the landscape, suppressing some states and favoring others, depending upon the mutational context. These results demonstrate that disrupting interactions in a region of the protein with highly cooperative, unfrustrated folding has very profound effects on the roughness of the folding landscape, whereas the effects are less pronounced for an energetically equivalent substitution in an already frustrated region.

High-Pressure SAXS Study of Folded and Unfolded Ensembles of Proteins

A structural interpretation of the thermodynamic stability of proteins requires an understanding of the structural properties of the unfolded state. High-pressure small-angle x-ray scattering was used to measure the effects of temperature, pressure, denaturants, and stabilizing osmolytes on the radii of gyration of folded and unfolded state ensembles of staphylococcal nuclease. A set of variants with the internal Val-66 replaced with Ala, Tyr, or Arg was used to examine how changes in the volume and polarity of an internal microcavity affect the dimensions of the native state and the pressure sensitivity of the ensemble. The unfolded state ensembles achieved for these proteins with high pressure were more compact than those achieved at high temperature, and were all very sensitive to the presence of urea and glycerol. Substitutions at the hydrophobic core detectably altered the conformation of the protein, even in the folded state. The introduction of a charged residue, such as Arg, inside the hydrophobic interior of a protein could dramatically alter the structural properties, even those of the unfolded state. The data suggest that a charge at an internal position can interfere with the formation of transient hydrophobic clusters in the unfolded state, and ensure that the pressure-unfolded form of a protein occupies the maximum volume possible. Only at high temperatures does the radius of gyration of the unfolded state ensemble approach the value for a statistical random coil.