Toward the quantum chemical calculation of nuclear magnetic resonance chemical shifts of proteins

Despite the many protein structures solved successfully by nuclear magnetic resonance (NMR) spectroscopy, quality control of NMR structures is still by far not as well established and standardized as in crystallography. Therefore, there is still the need for new, independent, and unbiased evaluation tools to identify problematic parts and in the best case also to give guidelines that how to fix them. We present here, quantum chemical calculations of NMR chemical shifts for many proteins based on our fragment-based quantum chemical method: the adjustable density matrix assembler (ADMA). These results show that (13)C chemical shifts of reasonable accuracy can be obtained that can already provide a powerful measure for the structure validation. (1)H and even more (15)N chemical shifts deviate more strongly from experiment due to the insufficient treatment of solvent effects and conformational averaging.     

1H, 13C, and 15N resonance assignment of photoactive yellow protein

Photoactive yellow protein (PYP) is involved in the negative phototactic response towards blue light of the bacterium Halorhodospira halophila. Here, we report nearly complete backbone and side chain ¹H, ¹³C and ¹⁵N resonance assignments at pH 5.8 and 20 °C of PYP in its electronic ground state.     

Relationship between nuclear magnetic resonance chemical shift and protein secondary structure

An analysis of the 1H nuclear magnetic resonance chemical shift assignments and secondary structure designations for over 70 proteins has revealed some very strong and unexpected relationships. Similar studies, performed on smaller databases, for 13C and 15N chemical shifts reveal equally strong correlations to protein secondary structure. Among the more interesting results to emerge from this work is the finding that all 20 naturally occurring amino acids experience a mean alpha-1H upfield shift of 0.39 parts per million (from the random coil value) when placed in a helical configuration. In a like manner, the alpha-1H chemical shift is found to move downfield by an average of 0.37 parts per million when the residue is placed in a beta-strand or extended configuration. Similar changes are also found for amide 1H, carbonyl 13C, alpha-13C and amide 15N chemical shifts. Other relationships between chemical shift and protein conformation are also uncovered; in particular, a correlation between helix dipole effects and amide proton chemical shifts as well as a relationship between alpha-proton chemical shifts and main-chain flexibility. Additionally, useful relationships between alpha-proton chemical shifts and backbone dihedral angles as well as correlations between amide proton chemical shifts and hydrogen bond effects are demonstrated.     

Light induces destabilization of photoactive yellow protein

To understand the effect of visible light on the stability of photoactive yellow protein (PYP), urea denaturation experiments were performed with PYP in the dark and with PYP(M) under continuous illumination. The urea concentrations at the midpoint of denaturation were 5.26 +/- 0.29 and 3.77 +/- 0.19 M for PYP and PYP(M), respectively, in 100 mM acetate buffer, and 5.26 +/- 0.24 and 4.11 +/- 0.12 M for PYP and PYP(M), respectively, in 100 mM citrate buffer. The free energy change upon denaturation (DeltaG(D)(H2O)), obtained from the denaturation curve, was 11.0 +/- 0.4 and 7.6 +/- 0.2 kcal/mol for PYP and PYP(M), respectively, in acetate buffer, and 11.5 +/- 0.3 and 7.8 +/- 0.1 kcal/mol for PYP and PYP(M), respectively, in citrate buffer. Even though the DeltaG(D)(H2O) value for PYP(M) is almost identical in the two buffer systems, the urea concentration at the midpoint of denaturation is lower in acetate buffer than in citrate buffer. Although their CD spectra indicate that the protein conformations of the denatured states of PYP and PYP(M) are indistinguishable, the configurations of the chromophores in their denatured structures are not necessarily identical. Both denatured states are interconvertible through PYP and PYP(M). Therefore, the free energy difference between PYP and PYP(M) is 3.4-3.7 kcal/mol for the protein moiety, plus the additional contribution from the difference in configuration of the chromophore.     

Static and transient hydrogen-bonding interactions in recombinant desulfatohirudin studied by 1H nuclear magnetic resonance measurements of amide proton exchange rates and pH-dependent chemical shifts

With proton nuclear magnetic resonance spectroscopy at 22 degrees C and pD 4.5, individual exchange rates in the range from 2 X 10(-5) to 1 X 10(-1) min-1 were observed for 23 amide protons in recombinant desulfatohirudin. The remaining 38 backbone amide protons exchange more rapidly than 1 X 10(-1) min-1. All 23 slowly exchanging protons are located in the polypeptide segment from residue 4 to residue 42, which forms a well-defined globular domain. Three different breathing modes of this molecular region are manifested in the exchange data, which appear to be correlated with the location of the three disulfide bonds. Chemical shift changes larger than 0.15 ppm between pH 2.5 and pH 5.0 arising from through-space interactions with carboxyl groups were observed for seven backbone amide protons. Two of these shifts can be explained by hydrogen bonds in the core of the protein, Gly 25 NH-Glu 43 O epsilon and Ser 32 NH-Asp 33 O delta, and two others by intraresidual NH-O epsilon interactions in Glu 61 and Glu 62. The remaining three pH shifts for Glu 35, Cys 39, and Ile 59 imply the existence of transient interactions between the molecular core and the flexible C-terminal segment 49-65, which have so far not been characterized by nuclear Overhauser effects or other conformational constraints.     

Redox-dependent structure change and hyperfine nuclear magnetic resonance shifts in cytochrome c

Proton nuclear magnetic resonance assignments for reduced and oxidized equine cytochrome c show that many individual protons exhibit different chemical shifts in the two protein forms, reflecting diamagnetic shift effects due to structure change, and in addition contact and pseudocontact shifts that occur only in the paramagnetic oxidized form. To evaluate the chemical shift differences (delta delta) for structure change, we removed the pseudocontact shift contribution by a calculation based on knowledge of the electron spin g tensor. The g-tensor parameters were determined from the delta delta values of a large set (64) of C alpha H protons at well-defined spatial positions in the oxidized horse protein. The g-tensor calculation, when repeated using only 12 available C alpha H proton resonances for cytochrome c from tuna, proved to be remarkably stable. The largest principal value of the g tensor (gz) falls precisely along the ligand bond between the heme iron and methionine-80 sulfur, while gx and gy closely match the natural heme axes defined by the pyrrole nitrogens. The derived g tensor was then used together with spatial coordinates for the oxidized form to calculate the pseudocontact shift contribution (delta pc) to proton resonances at 400 identifiable sites throughout the protein, so that the redox-dependent chemical shift discrepancy, delta delta-delta pc, could be evaluated. Large residual changes in chemical shift define the Fermi contact shifts, which are found as expected to be limited to the immediate covalent structure of the heme and its ligands and to be asymmetrically distributed over the heme. Smaller chemical shift discrepancies point to a concerted change, involving residues 39-43 and 50-60 (bottom of the protein), and to other changes in the immediate vicinity of the heme ligands. Also, the three internal water molecules are implicated in redox sensitivity. The residues found to change are in good but not perfect agreement with prior X-ray diffraction observations of subangstrom redox-related displacements in the tuna protein. The chemical shift discrepancies observed appear in the main to reflect structure-dependent diamagnetic shifts rather than hyperfine effects due to displacements in the pseudocontact shift field. Although 51 protons in 29 different residues exhibit significant chemical shift changes, the general impression is one of small structural adjustments to redox-dependent strain rather than sizeable structural displacements or rearrangements.     

Sequential resonance assignments in protein 1H nuclear magnetic resonance spectra. Basic pancreatic trypsin inhibitor

The assignment of the ¹H nuclear magnetic resonance spectrum of the basic pancreatic trypsin inhibitor with the use of two-dimensional ¹H nuclear magnetic resonance techniques at 500 MHz is described. The assignments are based entirely on the known amino acid sequence and the nuclear magnetic resonance data. Individual resonance assignments were obtained for all backbone and C^β protons, with the exception of those of Arg1, Pro2, Pro13 and the amide proton of Gly37. The side-chain resonance assignments are complete, with the exception of Pro2 and Pro13, the N^δ protons of Asn44 and the peripheral protons of the lysine residues and all but two of the arginine residues.     

Predicting the signaling state of photoactive yellow protein

As a bacterial blue light sensor the photoactive yellow protein (PYP) undergoes conformational changes upon signal transduction. The absorption of a photon triggers a series of events that are initially localized around the protein chromophore, extends to encompass the whole protein within microseconds, and leads to the formation of the transient pB signaling state. We study the formation of this signaling state pB by molecular simulation and predict its solution structure. Conventional straightforward molecular dynamics is not able to address this formation process due to the long (microsecond) timescales involved, which are (partially) caused by the presence of free energy barriers between the metastable states. To overcome these barriers, we employed the parallel tempering (or replica exchange) method, thus enabling us to predict qualitatively the formation of the PYP signaling state pB. In contrast to the receptor state pG of PYP, the characteristics of this predicted pB structure include a wide open chromophore-binding pocket, with the chromophore and Glu(46) fully solvent-exposed. In addition, loss of alpha-helical structure occurs, caused by the opening motion of the chromophore-binding pocket and the disruptive interaction of the negatively charged Glu(46) with the backbone atoms in the hydrophobic core of the N-terminal cap. Recent NMR experiments agree very well with these predictions.     

Absorption spectra of photoactive yellow protein chromophores in vacuum

The absorption spectra of two photoactive yellow protein model chromophores have been measured in vacuum using an electrostatic ion storage ring. The absorption spectrum of the isolated chromophore is an important reference for deducing the influence of the protein environment on the electronic energy levels of the chromophore and separating the intrinsic properties of the chromophore from properties induced by the protein environment. In vacuum the deprotonated trans-thiophenyl-p-coumarate model chromophore has an absorption maximum at 460 nm, whereas the photoactive yellow protein absorbs maximally at 446 nm. The protein environment thus only slightly blue-shifts the absorption. In contrast, the absorption of the model chromophore in aqueous solution is significantly blue-shifted (lambda(max) = 395 nm). A deprotonated trans-p-coumaric acid has also been studied to elucidate the effect of thioester formation and phenol deprotonation. The sum of these two changes on the chromophore induces a red shift both in vacuum and in aqueous solution.     

Time-resolved resonance raman structural studies of the pB' intermediate in the photocycle of photoactive yellow protein

Time-resolved resonance Raman spectroscopy is used to obtain chromophore vibrational spectra of the pR, pB', and pB intermediates during the photocycle of photoactive yellow protein. In the pR spectrum, the C8-C9 stretching mode at 998 cm(-1) is approximately 60 cm(-1) lower than in the dark state, and the combination of C-O stretching and C7H=C8H bending at 1283 cm(-1) is insensitive to D2O substitution. These results indicate that pR has a deprotonated, cis chromophore structure and that the hydrogen bonding to the chromophore phenolate oxygen is preserved and strengthened in the early photoproduct. However, the intense C7H=C8H hydrogen out-of-plane (HOOP) mode at 979 cm(-1) suggests that the chromophore in pR is distorted at the vinyl and adjacent C8-C9 bonds. The formation of pB' involves chromophore protonation based on the protonation state marker at 1174 cm(-1) and on the sensitivity of the COH bending at 1148 cm(-1) as well as the combined C-OH stretching and C7H=C8H bending mode at 1252 cm(-1) to D2O substitution. The hydrogen out-of-plane Raman intensity at 985 cm(-1) significantly decreases in pB', suggesting that the pR-to-pB' transition is the stage where the stored photon energy is transferred from the distorted chromophore to the protein, producing a more relaxed pB' chromophore structure. The C=O stretching mode downshifts from 1660 to 1651 cm(-1) in the pB'-to-pB transition, indicating the reformation of a hydrogen bond to the carbonyl oxygen. Based on reported x-ray data, this suggests that the chromophore ring flips during the transition from pB' to pB. These results confirm the existence and importance of the pB' intermediate in photoactive yellow protein receptor activation.     

Initial events in the photocycle of photoactive yellow protein

The light-induced isomerization of a double bond is the key event that allows the conversion of light energy into a structural change in photoactive proteins for many light-mediated biological processes, such as vision, photosynthesis, photomorphogenesis, and photo movement. Cofactors such as retinals, linear tetrapyrroles, and 4-hydroxy-cinnamic acid have been selected by nature that provide the essential double bond to transduce the light signal into a conformational change and eventually, a physiological response. Here we report the first events after light excitation of the latter chromophore, containing a single ethylene double bond, in a low temperature crystallographic study of the photoactive yellow protein. We measured experimental phases to overcome possible model bias, corrected for minimized radiation damage, and measured absorption spectra of crystals to analyze the photoproducts formed. The data show a mechanism for the light activation of photoactive yellow protein, where the energy to drive the remainder of the conformational changes is stored in a slightly strained but fully cis-chromophore configuration. In addition, our data indicate a role for backbone rearrangements during the very early structural events.     

Probing protein dynamics by nuclear magnetic resonance

Proteins are dynamic molecules that often undergo conformational changes while performing their specific functions, such as target recognition, ligand binding and catalysis. NMR spectroscopy is uniquely suited to study protein dynamics, because site-specific information can be obtained for motions that span a broad range of time scales. The information obtained from NMR dynamics experiments has provided insights into specific structural changes or conformational energetics associated with molecular function. In the last decade, a number of new advancements in NMR methodologies have further extended our ability to characterize protein dynamics. Here, we present an overview of current NMR technology that is used to monitor the dynamic properties of proteins.     

Light-induced global conformational change of photoactive yellow protein in solution

The light-induced global conformational change of photoactive yellow protein was directly observed by small-angle X-ray scattering (SAXS). The N-terminal 6, 15, or 23 amino acid residues were enzymatically truncated (T6, T15, or T23, respectively), and their near-UV intermediates were accumulated under continuous illumination for SAXS measurements. The Kratky plot demonstrated that illumination induced partial loss of globularity. The change in globularity was marked in T6 but very small in T15 and T23, suggesting that structural change in positions 7-15 mainly reduces the globularity. The radius of gyration (R(g)) estimated by Guinier plot was increased by 1.1 A for T6 and 0.7 A for T15 and T23 upon illumination. As T23 lacks most of the N-terminal loop, structural change in the main part composed of the PAS core, helical connector, and beta-scaffold caused an increase of R(g) by 0.7 A. The structural change of positions 7-15 caused an additional increase by 0.4 A. The decrease of R(g) upon truncation of positions 7-15 for dark state was 0.3 A, while that for the intermediate was 0.7 A, suggesting that this region moves outward on formation of the intermediate. These results indicate that a light-induced structural change of PYP takes place in the main part and N-terminal 15 amino acid residues. The former induces only dimensional increase, but the latter results in additional change in shape.     

Resonance Raman spectroscopy and quantum chemical calculations reveal structural changes in the active site of photoactive yellow protein

Photoactive yellow protein (PYP) is a bacterial photoreceptor containing a 4-hydroxycinnamyl chromophore. Photoexcitation of PYP triggers a photocycle that involves at least two intermediate states: an early red-shifted PYP(L) intermediate and a long-lived blue-shifted PYP(M) intermediate. In this study, we have explored the active site structures of these intermediates by resonance Raman spectroscopy. Quantum chemical calculations based on a density functional theory are also performed to simulate the observed spectra. The obtained structure of the chromophore in PYP(L) has cis configuration and no hydrogen bond at the carbonyl oxygen. In PYP(M), the cis chromophore is protonated at the phenolic oxygen and forms the hydrogen bond at the carbonyl group. These results allow us to propose structural changes of the chromophore during the photocycle of PYP. The chromophore photoisomerizes from trans to cis configuration by flipping the carbonyl group to form PYP(L) with minimal perturbation of the tightly packed protein interior. Subsequent conversion to PYP(M) involves protonation on the phenolic oxygen, followed by rotation of the chromophore as a whole. This large motion of the chromophore is potentially correlated with the succeeding global conformational changes in the protein, which ultimately leads to transduction of a biological signal.     

Mimic of photocycle by a protein folding reaction in photoactive yellow protein

The blue light receptor photoactive yellow protein (PYP) displays rhodopsin-like photochemistry based on the trans to cis photoisomerization of its p-coumaric acid chromophore. Here, we report that protein refolding from the acid-denatured state of PYP mimics the last photocycle transition in PYP. This implies a direct link between transient protein unfolding and photosensory signal transduction. We utilize this link to study general issues in protein folding. Chromophore trans to cis photoisomerization in the acid-denatured state strongly decelerates refolding, and converts the pH dependence of the barrier for refolding from linear to nonlinear. We propose transition state movement to explain this phenomenon. The cis chromophore significantly stabilizes the acid-denatured state, but acidification of PYP results in the accumulation of the acid-denatured state containing a trans chromophore. This provides a clear example of kinetic control in a protein unfolding reaction. These results demonstrate the power of PYP as a light-triggered model system to study protein folding.