Early Strides in NMR Dynamics Measurements.

The study of protein dynamics using the measurement of relaxation times by NMR was based on a set of studies in the mid-20th century that outlined theories and methods. However, the complexity of protein NMR was such that these simple experiments were not practical for application to proteins. The advent of techniques in the 1980s for isotopic labeling of proteins meant that pulse sequences could now be applied in multidimensional NMR experiments to enable per-residue information about the local relaxation times. One of the earliest advances was published in Biochemistry in 1989. The paper "Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease" by Lewis Kay, Dennis Torchia, and Ad Bax delineated a set of pulse sequences that are used with minor modifications even today. This paper, with others from a limited number of other laboratories, forms the basis for the experimental determination of the backbone dynamics of proteins. The biological insights obtained from such measurements have only increased in the past 30 years. Sometimes, the best and perhaps only way to advance a field is an advancement in the technical capabilities that allows new perspectives to be reached.

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins.

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.

Folding stability and cooperativity of the three forms of 1-110 residues fragment of staphylococcal nuclease.

Folding stability and cooperativity of the three forms of 1-110 residues fragment of staphylococcal nuclease (SNase110) have been studied by various biophysical and NMR methods. Samples of G-88W- and V-66W-mutant SNase110, namely G-88W110 and V-66W110, in aqueous solution and SNase110 in 2.0 M TMAO are adopted in this study. The unfolding transitions and folded conformations of the three SNase fragments were detected by far- and near-ultraviolet circular dichroism and intrinsic tryptophan fluorescence measurements. The tertiary structures and internal motions of the fragments were determined by NMR spectroscopy. Both G-88W and V-66W single mutations as well as a small organic osmolyte (Trimethylamine N-oxide, TMAO) can fold the fragment into a native-like conformation. However, the tertiary structures of the three fragments exhibit different degrees of folding stability and compactness. G-88W110 adopts a relatively rigid structure representing a most stable native-like beta-subdomain conformation of the three fragments. V-66W110- and TMAO-stabilized SNase110 produce less compact structures having a less stable "beta-barrel" structural region. The different folding status accounts for the different backbone dynamic and urea-unfolding transition features of the three fragments. The G-20I/G-29I-mutant variants of the three fragments have provided the evidence that the folding status is correlated closely to the packing of the beta-strands in the beta-barrel of the fragments. The native-like beta-barrel structural region acts as a nonlocal nucleus for folding the fragment. The tertiary folding of the three fragments is initiated by formation of the local nucleation sites at two beta-turn regions, I-18-D-21 and Y-27-Q-30, and developed by the formation of a nonlocal nucleation site at the beta-barrel region. The formation of beta-barrel and overall structure is concerted, but the level of cooperativity is different for the three 1-110 residues SNase fragments.

High apparent dielectric constant inside a protein reflects structural reorganization coupled to the ionization of an internal Asp.

The dielectric properties of proteins are poorly understood and difficult to describe quantitatively. This limits the accuracy of methods for structure-based calculation of electrostatic energies and pK(a) values. The pK(a) values of many internal groups report apparent protein dielectric constants of 10 or higher. These values are substantially higher than the dielectric constants of 2-4 measured experimentally with dry proteins. The structural origins of these high apparent dielectric constants are not well understood. Here we report on structural and equilibrium thermodynamic studies of the effects of pH on the V66D variant of staphylococcal nuclease. In a crystal structure of this protein the neutral side chain of Asp-66 is buried in the hydrophobic core of the protein and hydrated by internal water molecules. Asp-66 titrates with a pK(a) value near 9. A decrease in the far UV-CD signal was observed, concomitant with ionization of this aspartic acid, and consistent with the loss of 1.5 turns of alpha-helix. These data suggest that the protein dielectric constant needed to reproduce the pK(a) value of Asp-66 with continuum electrostatics calculations is high because the dielectric constant has to capture, implicitly, the energetic consequences of the structural reorganization that are not treated explicitly in continuum calculations with static structures.

Pressure-dependent transition in protein dynamics at about revealed by molecular dynamics simulation.

Molecular dynamics simulations of a crystalline protein, Staphylococcal nuclease, over the pressure range 1 bar to 15 kbar reveal a qualitative change in the internal protein motions at approximately 4 kbar. This change involves the existence of two linear regimes in the mean-square displacement for internal protein motion, (P) with a twofold decrease in the slope for P>4 kbar. The major effect of pressure on the dynamics is a loss, with increasing pressure of large amplitude, collective protein modes below 2 THz effective frequency, accompanied by restriction of large-scale solvent translational motion.

Packaging of up to 240 subunits of a 17 kDa nuclease into the interior of recombinant hepatitis B virus capsids.

The icosahedral nucleocapsid of hepatitis B virus (HBV) consists of multiple subunits of a single 183 amino acids (aa) core protein encasing the viral genome. However, recombinant core protein alone also forms capsid-like particles. We have recently shown that a 238 aa protein centrally inserted into the core protein can be displayed on the particle surface. Here we demonstrate that replacement of the C-terminal basic domain by the 17 kDa Staphylococcus aureus nuclease also yields particles but that in these the foreign domains are located in the interior. The packaged nuclease is enzymatically active, and the chimeric protein forms mosaic particles with the wild-type core protein. Hence the HBV capsid is useful as a molecular platform which, dependent on the fusion site, allows foreign protein domains to either be packaged into or be exposed on the exterior of the particle. These results are of relevance for the use of the HBV capsid as a vaccine carrier, and as a target for antiviral therapy.

Structure-based calculation of the equilibrium folding pathway of proteins. Correlation with hydrogen exchange protection factors.

A new statistical thermodynamic formalism has been developed in order to describe the equilibrium folding pathway of proteins. The resulting formalism allows calculation of the probabilities that individual amino acid residues will be in a native or native-like conformation for any given degree of folding of the protein molecule. The residue probabilities are defined by the probability distribution of conformational states and can be used to calculate experimental quantities like native-state, hydrogen exchange protection factors. A combinatorial algorithm aimed at generating a large ensemble of conformational states (10(4) to 10(6)) using the native structure as a template has been developed. The Gibbs energy and corresponding probability of each conformational state is estimated by using a previously developed structural parametrization of the energetics. The approach has been applied to five different proteins: hen egg-white lysozyme, equine lysozyme, bovine pancreatic trypsin inhibitor, staphylococcal nuclease and turkey ovomucoid third domain. The validity of the approach has been tested by comparing predicted and experimental hydrogen exchange protection factors. It is shown that for the above proteins 76%, 73%, 74%, 78% and 81% of all observed protection factors are predicted correctly. Furthermore, on average, the magnitude of the predicted protection factors, expressed as apparent free energies per residue deviate less than 1 kcal/mol from those obtained experimentally. These results represent the first attempt at predicting both the location and magnitude of hydrogen exchange protection factors from the high-resolution structure of a protein. The good agreement between experimental and predicted values has permitted a close examination of the nature of the equilibrium folding intermediates existing under conditions of maximal stability of the native state.

NMR docking of the competitive inhibitor thymidine 3',5'-diphosphate into the X-ray structure of staphylococcal nuclease.

In the X-ray structure of the ternary staphylococcal nuclease-Ca(2+)-3',5'-pdTp complex, the conformation of the bound inhibitor 3',5'-pdTp is distorted by Lys-70* and Lys-71* from an adjacent molecule of the enzyme in the crystal lattice (Loll, P. J. and Lattman, E. E. Proteins 5:183-201, 1989; Serpersu, E. H., Hibler, D. W., Gerlt, J. A., and Mildvan, A. S. Biochemistry 28:1539-1548, 1989). Since this interaction does not occur in solution, the NMR docking procedure has been used to correct this problem. Based on 8 Co(2+)-nucleus distances measured by paramagnetic effects on T1, and 9 measured and 45 lower limit interproton distances determined by 1D and 2D NOE studies of the ternary Ca2+ complex, the conformation of enzyme-bound 3',5'-pdTp is high-anti (chi = 58 +/- 10 degrees) with a C2' endo/O1' endo sugar pucker (delta = 143 +/- 2 degrees), (-) synclinal about the C3'-O3' bond (epsilon = 273 +/- 4 degrees), trans, gauche about the C4'-C5' bond (gamma = 301 +/- 29 degrees) and either (-) or (+) clinal about the C5'-O5' bond (beta = 92 +/- 8 degrees or 274 +/- 3 degrees). The structure of 3',5'-pdTp in the crystalline complex differs due to rotations about the C4'-C5' bond (gamma = 186 +/- 12 degrees, gauche, trans) and the C5'-O5' bond [beta = 136 +/- 10 degrees, (+) anticlinal]. The undistorted conformation of enzyme-bound metal-3',5'-pdTp determined by NMR was docked into the X-ray structure of the enzyme, using 19 intermolecular NOEs from ring proton resonances of Tyr-85, Tyr-113, and Tyr-115 to proton resonances of the inhibitor. van der Waals overlaps were then removed by energy minimization. Subsequent molecular dynamics and energy minimization produced no significant changes, indicating the structure to be in a global rather than in a local minimum. While the metal-coordinated 5'-phosphate of the NMR-docked structure of 3',5'-pdTp overlaps with that in the X-ray structure, and similarly receives bifunctional hydrogen bonds from both Arg-35 and Arg-87, the thymine, deoxyribose, and 3'-phosphate are significantly displaced from their positions in the X-ray structure, with the 3'-phosphate receiving hydrogen bonds from Lys-49 rather than from Lys-84 and Tyr-85. The repositioned thymine ring permits hydrogen bonding to the phenolic hydroxyl of Tyr-115.(ABSTRACT TRUNCATED AT 400 WORDS)

Mutational tests of the NMR-docked structure of the staphylococcal nuclease-metal-3',5'-pdTp complex.

In the X-ray structure of the staphylococcal nuclease-Ca(2+)-3',5'-pdTp complex, the conformation of the inhibitor 3',5'-pdTp is distorted by Lys-70* and Lys-71* from an adjacent molecule of staphylococcal nuclease (Loll, P.J., Lattman, E.E. Proteins 5:183-201, 1989). In order to correct this crystal packing problem, the solution conformation of enzyme-bound 3',5'-pdTp in the staphylococcal nuclease-metal-pdTp complex determined by NMR methods was docked into the X-ray structure of the enzyme [Weber, D.J., Serpersu, E.H., Gittis, A.G., Lattman, E.E., Mildvan, A.S. (preceding paper)]. In the NMR-docked structure, the 5'-phosphate of 3',5'-pdTp overlaps with that in the X-ray structure. However, the 3'-phosphate accepts a hydrogen bond from Lys-49 (2.89 A) rather than from Lys-84 (8.63 A), and N3 of thymine donates a hydrogen bond to the OH of Tyr-115 (3.16 A) which does not occur in the X-ray structure (5.28 A). These interactions have been tested by binding studies of 3',5'-pdTp, Ca2+, and Mn2+ to the K49A, K84A, and Y115A mutants of staphylococcal nuclease using water proton relaxation rate and EPR methods. Each mutant was fully active and structurally intact, as found by CD and two-dimensional NMR spectroscopy, but bound Ca2+ 9.1- to 9.9-fold more weakly than the wild-type enzyme. While the K84A mutation did not significantly weaken 3',5'-pdTp binding to the enzyme (1.5 +/- 0.7 fold), the K49A mutation weakened 3',5'-pdTp binding to the enzyme by the factor of 4.4 +/- 1.8-fold. Similarly, the Y115A mutation weakened 3',5'-pdTp binding to the enzyme 3.6 +/- 1.6-fold. Comparable weakening effects of these mutations were found on the binding of Ca(2+)-3',5'-pdTp. These results are more readily explained by the NMR-docked structure of staphylococcal nuclease-metal-3',5'-pdTp than by the X-ray structure.

Molecular dynamics study of the stability of staphylococcal nuclease mutants: component analysis of the free energy difference of denaturation.

The stability of two mutants G88V (Gly-88-->Val) and A69T (Ala-69-->Thr) of staphylococcal nuclease was analyzed by molecular dynamics simulations. The calculated free energy differences of denaturation for G88V and A69T were -1.1 and -2.8 kcal/mol, respectively. These values are in good agreement with the experimental values. The free energy differences divided into electrostatic and van der Waals components were analyzed. These two mutants are mainly destabilized due to van der Waals interactions. There is little difference between the electrostatic contribution to the free energy change in the native state and that in the denatured state. In each mutant structure, a small cavity appears in the vicinity of the perturbed residue. It is suggested that intramolecular van der Waals interactions of the mutants are weaker than those of the wild-type. Furthermore, analyses of the contributions of each residue near the perturbed residue and of water to the free energy difference of denaturation suggest that the interaction between water and the perturbed residue plays a very important role in the stability of staphylococcal nuclease, and that a small hydrophobic core consisting of the three aromatic rings (Tyr-27, Phe-34, Phe-76) and the side chain of Met-32 is also important for the stability.

Comparison of conformational features of staphylococcal nuclease in ternary complexes with pdTp, pdGp, and nitrophenyl-pdTp.

The conformations of wild-type staphylococcal nuclease (SNase) in the ternary complexes with thymidine 3',5'-bisphosphate (pdTp), 2'-deoxyguanine 3',5'-bisphosphate (pdGp), and thymidine 3'-phosphate 5'-(p-nitrophenylphosphate) (NpdTp) with Ca2+ were examined by two-dimensional NMR NOESY and ROESY experiments. The results of these experiments indicate that the conformational features of the SNase are quite similar in the three ternary complexes. This suggests that the conformational features of SNase, in these ternary complexes, are not strongly dependent on whether the 5'-phosphate is a mono- or diester. This is in contrast to our prior studies on substitutions of active site charged amino acids which indicated that the conformational features of SNase in the ternary complex are quite sensitive to substitutions for active site charged amino acids (Hibler et al., 1987; Wilde et al., 1988; Pourmotabbed et al., 1990). The similarity of the SNase conformational features in the ternary complexes with pdTp and pdGp indicates that the features of the nucleotide bound at the active site are not strong determinants of the enzyme conformation in the ternary complexes. These conclusions are in general agreement with the results on pdApdT ternary complexes with SNase which suggested that it is the conformational features of the bound nucleic acid which determine the differences in catalysis observed for SNase with different substrates (Weber et al., 1991), more so than the conformational features of the enzyme.

Fluorescence and conformational stability studies of Staphylococcus nuclease and its mutants, including the less stable nuclease-concanavalin A hybrids.

We report steady-state and time-resolved fluorescence studies with the single tryptophan protein, Staphylococcus aureus A, and several of its site-directed mutants. A couple of these mutants, nuclease-conA and nuclease-conA-S28G (which are hybrid proteins containing a six amino acid beta-turn substitute from concanavalin A), are found to have a much lower thermodynamic stability than the wild type. The thermal transition temperatures for nuclease-conA and S28G are 32.8 and 30.5 degrees C, which are about 20 degrees C lower than the Tm for wild-type nuclease A. These mutant proteins also are denatured by a much lower concentration of the denaturants urea and guanidine hydrochloride. We also show that an unfolding transition in the structure of the nuclease-conA hybrids can be induced by relatively low hydrostatic pressure (approximately 700 bar). The free energy for unfolding of nuclease-conA (and nuclease-conA-S28G) is found to be only 1.4 kcal/mol (and 1.2 kcal/mol) by thermal, urea, guanidine hydrochloride, and pressure unfolding. Time-resolved fluorescence intensity and anisotropy measurements with nuclease-conA-S28G show the temperature-, urea-, and pressure-perturbed states each to have a reduced average intensity decay time and to depolarize with a rotational correlation time of approximately 1.0 ns (as compared to a rotational correlation time of 11 ns for the native form of nuclease-conA-S28G at 20 degrees C).