Proton nuclear magnetic resonance as a probe of differences in structure between the C102T and F82S,C102T variants of iso-1-cytochrome c from the yeast Saccharomyces cerevisiae

NMR Studies of Retroviral Genome Packaging

Nearly all retroviruses selectively package two copies of their unspliced RNA genomes from a cellular milieu that contains a substantial excess of non-viral and spliced viral RNAs. Over the past four decades, combinations of genetic experiments, phylogenetic analyses, nucleotide accessibility mapping, in silico RNA structure predictions, and biophysical experiments were employed to understand how retroviral genomes are selected for packaging. Genetic studies provided early clues regarding the protein and RNA elements required for packaging, and nucleotide accessibility mapping experiments provided insights into the secondary structures of functionally important elements in the genome. Three-dimensional structural determinants of packaging were primarily derived by nuclear magnetic resonance (NMR) spectroscopy. A key advantage of NMR, relative to other methods for determining biomolecular structure (such as X-ray crystallography), is that it is well suited for studies of conformationally dynamic and heterogeneous systems—a hallmark of the retrovirus packaging machinery. Here, we review advances in understanding of the structures, dynamics, and interactions of the proteins and RNA elements involved in retroviral genome selection and packaging that are facilitated by NMR.

Redox state dependence of axial ligand dynamics in Nitrosomonas europaea cytochrome c552

Analysis of NMR spectra reveals that the heme axial Met ligand orientation and dynamics in Nitrosomonas europaea cytochrome c552 (Ne cyt c) are dependent on the heme redox state. In the oxidized state, the heme axial Met is fluxional, interconverting between two conformers related to each other by inversion through the Met δS atom. In the reduced state, there is no evidence of fluxionality, with the Met occupying one conformation similar to that seen in the homologous Pseudomonas aeruginosa cytochrome c551. Comparison of the observed and calculated pseudocontact shifts for oxidized Ne cyt c using the reduced protein structure as a reference structure reveals a redox-dependent change in the structure of the loop bearing the axial Met (loop 3). Analysis of nuclear Overhauser effects (NOEs) and existing structural data provides further support for the redox state dependence of the loop 3 structure. Implications for electron transfer function are discussed.

Structural determinants and mechanism of HIV-1 genome packaging

Like all retroviruses, the human immunodeficiency virus selectively packages two copies of its unspliced RNA genome, both of which are utilized for strand-transfer-mediated recombination during reverse transcription-a process that enables rapid evolution under environmental and chemotherapeutic pressures. The viral RNA appears to be selected for packaging as a dimer, and there is evidence that dimerization and packaging are mechanistically coupled. Both processes are mediated by interactions between the nucleocapsid domains of a small number of assembling viral Gag polyproteins and RNA elements within the 5'-untranslated region of the genome. A number of secondary structures have been predicted for regions of the genome that are responsible for packaging, and high-resolution structures have been determined for a few small RNA fragments and protein-RNA complexes. However, major questions regarding the RNA structures (and potentially the structural changes) that are responsible for dimeric genome selection remain unanswered. Here, we review efforts that have been made to identify the molecular determinants and mechanism of human immunodeficiency virus type 1 genome packaging.

Second virial coefficients as a measure of protein--osmolyte interactions

The cytoplasm contains high concentrations of cosolutes. These cosolutes include macromolecules and small organic molecules called osmolytes. However, most biophysical studies of proteins are conducted in dilute solutions. Two broad classes of models have been used to describe the interaction between osmolytes and proteins. One class focuses on excluded volume effects, while the other focuses on binding between the protein and the osmolyte. To better understand protein--smolyte interactions, we have conducted sedimentation equilibrium analytical ultracentrifugation experiments using ferricytochrome c as a model protein. From these experiments, we determined the second virial coefficients for a series of osmolytes. We have interpreted the second virial coefficient as a measure of both excluded volume and protein--osmolyte binding. We conclude that simple models are not sufficient to understand the interactions between osmolytes and proteins.

Sequential main-chain proton and carbon resonance assignments for wild-type yeast iso-1 ferrocytochrome c and identification of secondary structural elements

Yeast iso-1 cytochrome c is a naturally occurring protein that possesses an unusually reactive Cys102 that imbues iso-1 with a complicated solution chemistry which includes spontaneous dimerization and poorly characterized redox reactions. For this reason previous studies of this typical member of the c-type cytochromes have been relegated to variant proteins in which the 102 position has been mutated, with most common changes involving serine and threonine. However, we have determined sequential 1H resonance assignments for the wild-type native protein because it is the actual participant in yeast mitochondrial electron transfer processes and because the wild-type native protein should be the fundamental assignment basis. In addition to 1H resonance assignments for 97 of 106 amino acids, we have also provided an extensive database of long-range NOEs. Comparison of these NOEs and a chemical shift index based upon alpha-H resonances has lead to identification of solution secondary structural elements that are consistent with the solid-state crystal structure. Although there is currently no efficient expression system that would facilitate isotope labeling of iso-1 cytochrome c, we tried to assess the usefulness of future heteronuclear experiments by using natural-abundance 1H/13C HMQC experiments to unambiguously assign 35 alpha-C resonances.

Pseudocontact shifts as constraints for energy minimization and molecular dynamics calculations on solution structures of paramagnetic metalloproteins

The pseudocontact shifts of NMR signals, which arise from the magnetic susceptibility anisotropy of paramagnetic molecules, have been used as structural constraints under the form of a pseudopotential in the SANDER module of the AMBER 4.1 molecular dynamics software package. With this procedure, restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations can be performed on structural models by using pseudocontact shifts. The structure of the cyanide adduct of the Met80Ala mutant of the yeast iso-1-cytochrome c has been used for successfully testing the calculations. For this protein, a family of structures is available, which was obtained by using NOE and pseudocontact shifts as constraints in a distance geometry program. The structures obtained by REM and RMD calculations with the inclusion of pseudocontact shifts are analyzed.

Use of paramagnetic NMR probes for structural analysis in cytochrome c3 from Desulfovibrio vulgaris

The dipolar field generated by each of the four haems in the tetrahaem ferricytochrome c3 from Desulfovibrio vulgaris (Hildenborough) (c3DvH) is determined by means of a novel procedure. In this method the 13C chemical shifts of the nuclei directly bound to the haems are used to determine the in-plane orientations of the rhombic perturbation in each of the four haems with respect to a model of molecular orbitals of e(g) symmetry which are subject to a rhombic perturbation [Turner, D. L., Salgueiro, C. A., Schenkels, P., LeGall, J. & Xavier, A. V. (1995) Biochim. Biophys. Acta 1246, 24-281. These orientations, together with the components of the magnetic susceptibility tensors obtained from the EPR g values and the crystal structure of c3DvH, can be used to calculate the dipolar shifts induced by each haem throughout the protein. Thus the observed 13C paramagnetic shifts of the c3DvH haem substituents were fitted considering both the pseudocontact and contact shifts of each haem simultaneously. The dipolar shifts calculated by this method were tested against the observed dipolar shifts for some amino acid residues strategically placed in the protein and also for the haem propionate groups. The effect of considering the calculated dipolar extrinsic shifts on the behaviour of the chemical shifts of the haem methyl groups in the intermediate stages of oxidation at different pH values was also analysed. The several tests applied to the calculated dipolar shifts have shown that the method is extremely useful for predicting chemical shifts as an aid to complete proton assignment, and to add further constraints in the refinement of solution structures of paramagnetic proteins and hence to probe subtle structural rearrangements around the haem pocket.

Analysis of the 1H-NMR chemical shifts of Cu(I)-, Cu(II)- and Cd-substituted pea plastocyanin. Metal-dependent differences in the hydrogen-bond network around the copper site

To compare cadmium-substituted plastocyanin with copper plastocyanin, the 1H-NMR spectra of CuI-, CuII- and Cd-plastocyanin from pea have been analyzed. Full assignments of the spectra of CuI- and Cd-plastocyanin indicate chemical shift differences up to 1 ppm. The affected protons are located in the four loops that surround the Cu site. The largest differences were found for protons in the hydrogen bond network which stabilizes this part of the protein. This suggests that the chemical shift differences are caused by very small but extensive structural changes in the network upon replacement of CuI by Cd. For CuII-plastocyanin the resonances of 72% of the protons observed in the CuI form have been identified. Protons within approximately 0.9 nm of the CuII were not observed due to fast paramagnetic relaxation. The protons between 0.9-1.7 nm from the CuII showed chemical shift differences up to 0.4 ppm compared to both CuI- and Cd-plastocyanin. These differences can be predicted assuming that they represent pseudocontact shifts. When corrected for the pseudocontact shift contribution, the CuII-plastocyanin chemical shifts were nearly all identical within error to those of the Cd form, but not of the CuI-plastocyanin, indicating that the CuII-plastocyanin structure, in as far as it can be observed, resembles Cd-rather than CuI-plastocyanin. In a single stretch of residues (64-69) chemical shift differences remained between all three forms after correction. The fact that pseudocontact shifts were observed for protons which were not broadened may be attributable to the weaker distance dependence of the pseudocontact shift effect compared to paramagnetic relaxation. This results in two shells around the Cu atom, an inner paramagnetic shell (0-0.9 nm), in which protons are not observed due to broadening, and an outer paramagnetic shell (0.9-1.7 nm), in which protons can be observed and show pseudocontact shifts. It is concluded that Cd-plastocyanin is a suitable redox-inactive substitute for Cu-plastocyanin.

Assignment of 1H and 13C hyperfine-shifted resonances for tuna ferricytochrome c

Tuna ferricytochrome c has been used to demonstrate the potential for completely assigning 1H and 13C strongly hyperfine-shifted resonances in metalloprotein paramagnetic centers. This was done by implementation of standard two-dimensional NMR experiments adapted to take advantage of the enhanced relaxation rates of strongly hyperfine-shifted nuclei. The results show that complete proton assignments of the heme and axial ligands can be achieved, and that assignments of several strongly shifted protons from amino acids located close to the heme can also be made. Virtually all proton-bearing heme 13C resonances have been located, and additional 13C resonances from heme vicinity amino acids are also identified. These results represent an improvement over previous proton resonance assignment efforts that were predicated on the knowledge of specific assignments in the diamagnetic protein and relied on magnetization transfer experiments in heterogeneous solutions composed of mixtures of diamagnetic ferrocytochrome c and paramagnetic ferricytochrome c. Even with that more complicated procedure, complete heme proton assignments for ferricytochrome c have never been demonstrated by a single laboratory. The results presented here were achieved using a more generally applicable strategy with a solution of the uniformly oxidized protein, thereby eliminating the requirement of fast electron self-exchange, which is a condition that is frequently not met.

Hydrogen isotope effects on the proton nuclear magnetic resonance spectrum of bovine ferricytochrome b5: axial hydrogen bonding involving the axial His-39 imidazole ligand

The potential role of hydrogen bonding interactions in modulating the molecular and electronic structure of the active site of solubilized bovine ferricytochrome b5 has been investigated by monitoring solvent isotope effects on proton-NMR spectral parameters. It is observed that the hyperfine shifts of both the heme prosthetic group and one coordinated His are sensitive, while those for the other axial His and non-coordinated residues are insensitive, to 2H for 1H exchange. Two types of isotope influences are characterized; one whose chemical shift influence is time-resolved on the NMR time scale, and involves a single proton on one axial ligand, and a second effect which involves multiple protons, is not time resolved, and influences primarily the heme. A large isotope effect on the hyperfine shift is identified for the C beta H signals of His-39 but not His-63. The exchangeable ring NH of His-39 is assigned, and the pH influence on the exchange properties of heme pocket labile protons, when compared to the rate of base catalyzed averaging of the His-39 C beta H isotope effect, lead to the conclusion that the axial hydrogen bond which is responsible for this isotope effect is that between His-39 ring NH and Gly-42 carbonyl. The more rapid exchange of labile protons with solvent for His-63 than His-39 confirms a less solvent accessible and stronger hydrogen bonded His-39 than His-63. The stronger His-39-Gly-42 than His-63-Phe-58 hydrogen bond involving the ring NH leads to more extensive His-39 imidazolate character and hence a stronger iron-His-39 than iron-His-63 bond. The much larger hyperfine shifts for His-39 than His-63 imidazole ring non-labile protons support the stronger bonding of the former ligand, and account for the orientation of the rhombic magnetic axes by His-39 rather than His-63. The solvent isotope effect on the heme leads to rotation of the prosthetic group about the His-Fe-His bond by approximately 0.5 degrees so as to shorten the 7-propionate link to Ser-64. This suggests that the hydrogen bonds between the 7-propionate group and Ser-64 are responsible for the effect.

Probing weakly polar interactions in cytochrome c

Theoretical, statistical, and model studies suggest that proteins are stabilized by weakly polar attractions between sulfur atoms and properly oriented aromatic rings. The two sulfur-containing amino acids, methionine and cysteine, occur frequently among functional alleles in random mutant libraries of Saccharomyces cerevisiae iso-1-cytochrome c genes at positions that form a weakly polar aromatic-aromatic interaction, the wild-type protein. To determine if a weakly polar sulfur-aromatic interaction replaced the aromatic-aromatic interaction, the structure and stability of two variants were examined. Phenylalanine 10, which interacts with tyrosine 97, was replaced by methionine and cysteine. The cysteine was modified to form the methionine and cysteine analog, S-methyl cysteine (CysSMe). Proton NMR studies indicate that changing Phe 10 to Met or CysSMe affects only local structure and that the structures of sulfur-containing variants are nearly identical. Analysis of chemical shifts and nuclear Overhauser effect data indicates that both sulfur-containing side chains are in position to form a weakly polar interaction with Tyr 97. The F10M and F10CSMe variants are 2-3 kcal mol-1 less stable than iso-1-cytochrome c at 300 K. Comparison of the stabilities of the F10M and F10CSMe variants allows evaluation of the potential weakly polar interaction between the additional sulfur atom of F10CSMe and the aromatic moiety of Tyr 97. The F10CSMe;C102T variant is 0.7 +/- 0.3 kcal mol-1 more stable than the F10M;C102T protein. The increased stability is explained by the difference in hydrophobicity of the sulfur-containing side chains. We conclude that any weakly polar interaction between the additional sulfur and the aromatic ring is too weak to detect or is masked by destabilizing contributions to the free energy of denaturation.

Changes in global stability and local structure of cytochrome c upon substituting phenylalanine-82 with tyrosine

We have examined the F82Y;C102T variant of Saccharomyces cerevisiae iso-1-cytochrome c using high-resolution proton nuclear magnetic resonance spectroscopy, chemical denaturation, and differential scanning calorimetry. Comparison of proton chemical shifts, paramagnetic shifts, and nuclear Overhauser effects indicates structural changes are localized to the vicinity of position 82. One alteration involves the rearrangement of the side chain of leucine-85. Using many more proton assignments than were available in the initial report [G. J. Pielak, R. A. Atkinson, J. Boyd, and R. J. P. Williams, Eur. J. Biochem. 177, 179-185 (1988)], a second alteration involving an interaction between arginine-13 and tyrosine-82 is observed. The interaction appears to involve a hydrogen bond with the eta-protons of arginine's guanido group acting as donor and tyrosine's phenolic eta-oxygen as acceptor. In spite of this potentially-stabilizing interaction, the free energy of denaturation decreases by approximately 2.4 kcal mol-1. Results are discussed with respect to alterations in the native and denatured states.

1H- and 13C-NMR investigation of redox-state-dependent and temperature-dependent conformation changes in horse cytochrome c

The redox-state dependent changes in chemical shift, which have been measured for almost 100 CHn groups in the 13C-NMR spectra of horse cytochrome c [Santos, H., and Turner, D. L. (1992) Eur. J. Biochem. 206, 721-728], have been used to investigate the nature of the redox-related change in conformation. Apart from the haem and its axial ligands, the shifts are found to be dominated by the electron-nuclear dipolar coupling in the oxidised form, as was the case in 1H-NMR studies. These pseudocontact shifts are well described by using an empirically determined magnetic susceptibility tensor in conjunction with atomic coordinates for the horse cytochrome c. The groups which fit least well are located in the vicinity of Trp59. Comparison between 1H and 13C shifts and their temperature dependence shows that the differences from expectation based on a single structure for both oxidation states are caused largely by changes in the diamagnetic contribution to the chemical shifts. Since these are different for 1H and 13C resonances they indicate, independently from crystal structure data, some redox-related movement of the protein under the haem. The significance of these results for understanding electron transfer pathways is discussed. Finally, the temperature dependence of the pseudocontact shifts in the range 30-50 degrees C is shown to be anomalous. Approximately half of the anomalous effect may be attributed to Zeeman mixing of the electronic wavefunctions with a spin-orbit coupling constant lambda = 241 cm-1, while the other half is attributed to thermal expansion of the protein.

Laser flash photolysis studies of electron transfer mechanisms in cytochromes: an aromatic residue at position 82 is not required for cytochrome c reduction by flavin semiquinones or electron transfer from cytochrome c to cytochrome oxidase

The influence of an aromatic side chain at position 82 of yeast iso-1-cytochrome c on the kinetics of its electron transfer reactions has been investigated using laser flash photolysis methods to compare a series of site-specific mutant cytochromes in their reduction by free flavin semiquinone and in electron transfer from reduced cytochrome to bovine cytochrome c oxidase. Although small (approximately 10%) but significant differences are observed between some of the mutants (S82, Y82, I82) and wild-type (F82) or G82 cytochrome in the second-order rate constant for reduction by lumiflavin semiquinone, these do not correlate with side-chain aromaticity. In the reaction between the ferrocytochromes and cytochrome c oxidase, significantly larger deviations from exponentiality are found for those mutants having aliphatic residues at position 82 than for wild type or Y82. We interpret the nonexponential behavior in terms of multiple orientations of the cytochromes within the oxidase binding site; the extent to which this occurs is apparently influenced by the character of the residue at position 82. However, a comparison of the average rate constants for electron transfer to cytochrome oxidase for the various mutants reveals that all are closely comparable to WT, except for I82 which is significantly slower (approximately threefold). These results, combined with those obtained previously from steady-state kinetic and thermodynamic measurements, suggest that the observed differences among the mutants are due to alterations in the mode of binding of the cytochrome to the oxidase, rather than to a specific requirement for the presence of an aromatic group at position 82.

Electron - proton coupling in cytochrome c studied using protein variants

An NMR study of the cytochrome c variant Asn52Ile is used to show how the redox state change in native cytochrome c is coupled to a rearrangement of a proton network which runs through the cytochrome c molecule. The substitution breaks the H-bond network and removes the coupling. The uncovering of this putative proton channel and the connection of changes to it with redox state changes of the iron centre of the protein allows a possible description of the way in which redox energy state changes can be coupled to energization and gating of protons in membranes.

Temperature-sensitive variants of Saccharomyces cerevisiae iso-1-cytochrome c produced by random mutagenesis of codons 43 to 54

In vitro random mutagenesis within the CYC1 gene from the yeast Saccharomyces cerevisiae was used to produce a library of mutants encompassing codons 43 to 54 of iso-1-cytochrome c. This region consists of an evolutionarily conserved structure within an evolutionarily diverse sequence. The library, on a low-copy-number yeast shuttle phagemid, was introduced into a yeast strain lacking cytochrome c. The ability of transformants harboring a functional cytochrome c to grow on the non-fermentable carbon source glycerol at 30 degrees C and 37 degrees C was used to determine the phenotype of nearly 1000 transformants. Approximately 90% of the missense mutants present in the library give rise to the wild-type phenotype, 7% result in the temperature-sensitive (Cycts) phenotype, and 3% give rise to the non-functional (Cyc-) phenotype. Phagemids from 20 Cycts and 30 Cyc- transformants were subjected to DNA sequence analysis. All the mutations occur within the targeted region. One-third of the mutants from Cyc- transformants and all the mutants from Cycts transformants are missense mutants. The remaining mutants from Cyc- transformants are nonsense or frame-shift mutants. Missense mutations within the codons for Gly45, Tyr46, Thr49, Asn52 or Ile53 alone are sufficient to produce temperature-sensitive behavior both in vivo and in the variant proteins. The deduced amino acid substitutions correlate remarkably well with side-chain dynamics, secondary structure and tertiary structure of the wild-type protein.

Constraints on amino acid substitutions in the N-terminal helix of cytochrome c explored by random mutagenesis

The interaction of the N- and C-terminal helices is a hallmark of the cytochrome c family. Oligodeoxyribonucleotide-directed random mutagenesis within the gene encoding the C102T protein variant of Saccharomyces cerevisiae iso-1-cytochrome c was used to generate a library of mutations at the evolutionary invariant residues Gly-6 and Phe-10 in the N-terminal helix. Transformation of this library (contained on a low-copy-number yeast shuttle phagemid) into a yeast strain lacking a functional cytochrome c, followed by selection for cytochrome c function, reveals that 4-10% of the 400 possible amino acid substitutions are compatible with function. DNA sequence analysis of phagemids isolated from transformants exhibiting the functional phenotype elucidates the requirements for a stable helical interface. Basic residues are not tolerated at position 6 or 10. There is a broad volume constraint for amino acids at position 6. The amino acid substitutions observed to be compatible with function at Phe-10 show that the hydrophobic effect alone is sufficient to promote helical association. There are severe constraints that limit the combinations consistent with function, but the number of functionally consistent combinations observed exemplifies the plasticity of proteins.