Electron trapping and reactions in rhamnose by ESR and ENDOR.

The main objective for a reinvestigation of rhamnose was to devise a mechanistic link between the trapped electron detected previously and the secondary radicals observed at 77 K and at room temperature. Single crystals of rhamnose were X-irradiated at temperatures between 15 and 300 K and examined using ESR, ENDOR, and field-swept ENDOR techniques. After low-temperature irradiation a C3 H-abstraction radical is formed following the visible light-induced decay of the trapped electron. This species was previously assigned erroneously to a C2 H-abstraction species. At temperatures above 120 K, this radical deprotonates at the C3 hydroxy group. Furthermore, a C2 H-abstraction radical is formed following the thermally induced decay of the trapped electron. The C2 and C3 H-abstraction radicals did not convert into each other. A third radical species formed at low temperatures is a C5 H-abstraction radical. It is unstable above 250 K and decays without any apparent successor. The C2 and C3 H-abstraction radicals are formed thermally and photochemically from the parent trapped electron. The conversions are mediated by hydrogen atoms formed intermediately or by elimination of hydride ions. The thermal decomposition pathway requires further studies, in particular with respect to the possible role of water. Recently, Box et al. analyzed the site of the trapped electron in rhamnose crystals. The present results support the results obtained by these authors (Radiat. Res. 121, 262 (1990)). In particular, trapped electron vs proton distances closely match the conversion mechanisms suggested.     

Hydrogenases in the "active" state: determination of g-matrix axes and electron spin distribution at the active site by 1H ENDOR spectroscopy

Hydrons and electrons are substrates for the enzyme hydrogenase, but cannot be observed in X-ray crystal structures. High-resolution 1H electron nuclear double resonance (ENDOR) spectroscopy offers a means to detect the distribution of protons and unpaired electrons. ENDOR spectra were recorded from frozen solutions of the nickel-iron hydrogenases of Desulfovibrio gigas and Desulfomicrobium baculatum, in the "active" state ("Ni-C" EPR signal) and analyzed by orientationally selective simulation methods. The experimental spectra were fitted using a structural model of the nickel-iron centre based on crystallographic results, allowing for differences in electron spin distribution as well as the spatial orientation of the g-matrix ( g-tensor), and anisotropic and isotropic hyperfine couplings of the protons nearest to the nickel ion. ENDOR signals, detected after complete deuterium exchange, were assigned to six protons of the cysteines bound to nickel. The assignment took advantage of the substitution of a selenium for a sulfur ligand, which occurs naturally between the [NiFeSe] and [NiFe] hydrogenases from Dm. baculatum and D. gigas, respectively, and was found to affect just two signals. The four signals with the largest hyperfine couplings, including isotropic contributions from 4.5 to 13.5 MHz, were assigned to the beta-methylene protons of the two terminal cysteine ligands, one of which is substituted by seleno-cysteine in [NiFeSe] hydrogenase. The electron spin is delocalized onto the nickel (50%) and its sulfur ligands, with a higher proportion on the terminal than the bridging ligands. The g-matrix was found to align with the active site in such a way that the g1- g2 plane is nearly coplanar (18.3 degrees) with the plane defined by nickel and three sulfur atoms, and the g2 axis deviates by 22.9 degrees from the vector between nickel and iron. Significantly for the reaction of the enzyme, direct evidence for the binding of hydrons at the active site was obtained by the detection of H/D-exchangeable ENDOR signals.     

Protein-cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: II. Geometry of the hydrogen bonds to the primary quinone formula by 1H and 2H ENDOR spectroscopy

The geometry of the hydrogen bonds to the two carbonyl oxygens of the semiquinone Q(A)(. -) in the reaction center (RC) from the photosynthetic purple bacterium Rhodobacter sphaeroides R-26 were determined by fitting a spin Hamiltonian to the data derived from (1)H and (2)H ENDOR spectroscopies at 35 GHz and 80 K. The experiments were performed on RCs in which the native Fe(2+) (high spin) was replaced by diamagnetic Zn(2+) to prevent spectral line broadening of the Q(A)(. -) due to magnetic coupling with the iron. The principal components of the hyperfine coupling and nuclear quadrupolar coupling tensors of the hydrogen-bonded protons (deuterons) and their principal directions with respect to the quinone axes were obtained by spectral simulations of ENDOR spectra at different magnetic fields on frozen solutions of deuterated Q(A)(. -) in H(2)O buffer and protonated Q(A)(. -) in D(2)O buffer. Hydrogen-bond lengths were obtained from the nuclear quadrupolar couplings. The two hydrogen bonds were found to be nonequivalent, having different directions and different bond lengths. The H-bond lengths r(OH) are 1.73 +/- 0.03 Angstrom and 1.60 +/- 0.04 Angstrom, from the carbonyl oxygens O(1) and O(4) to the NH group of Ala M260 and the imidazole nitrogen N(delta) of His M219, respectively. The asymmetric hydrogen bonds of Q(A)(. -) affect the spin density distribution in the quinone radical and its electronic structure. It is proposed that the H-bonds play an important role in defining the physical properties of the primary quinone, which affect the electron transfer processes in the RC.     

Electron paramagnetic resonance and electron nuclear double resonance spectroscopic identification and characterization of the tyrosyl radicals in prostaglandin H synthase 1

The tyrosyl radicals generated in reactions of ethyl hydrogen peroxide with both native and indomethacin-pretreated prostaglandin H synthase 1 (PGHS-1) were examined by low-temperature electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) spectroscopies. In the reaction of peroxide with the native enzyme at 0 degrees C, the tyrosyl radical EPR signal underwent a continuous reduction in line width and lost intensity as the incubation time increased, changing from an initial, 35-G wide doublet to a wide singlet of slightly smaller line width and finally to a 25-G narrow singlet. The 25-G narrow singlet produced by self-inactivation was distinctly broader than the 22-G narrow singlet obtained by indomethacin treatment. Analysis of the narrow singlet EPR spectra of self-inactivated and indomethacin-pretreated enzymes suggests that they reflect conformationally distinct tyrosyl radicals. ENDOR spectroscopy allowed more detailed characterization by providing hyperfine couplings for ring and methylene protons. These results establish that the wide doublet and the 22-G narrow singlet EPR signals arise from tyrosyl radicals with different side-chain conformations. The wide-singlet ENDOR spectrum, however, is best accounted for as a mixture of native wide-doublet and self-inactivated 25-G narrow-singlet species, consistent with an earlier EPR study [DeGray et al. (1992) J. Biol. Chem. 267, 23583-23588]. We conclude that a tyrosyl residue other than the catalytically essential Y385 species is most likely responsible for the indomethacin-inhibited, narrow-singlet spectrum. Thus, this inhibitor may function by redirecting radical formation to a catalytically inactive side chain. Either radical migration or conformational relaxation at Y385 produces the 25-G narrow singlet during self-inactivation. Our ENDOR data also indicate that the catalytically active, wide-doublet species is not hydrogen bonded, which may enhance its reactivity toward the fatty-acid substrate bound nearby.     

Role of the coordinating histidine in altering the mixed valency of Cu(A): an electron nuclear double resonance-electron paramagnetic resonance investigation

The binuclear Cu(A) site engineered into Pseudomonas aeruginosa azurin has provided a Cu(A)-azurin with a well-defined crystal structure and a CuSSCu core having two equatorial histidine ligands, His120 and His46. The mutations His120Asn and His120Gly were made at the equatorial His120 ligand to understand the histidine-related modulation to Cu(A), notably to the valence delocalization over the CuSSCu core. For these His120 mutants Q-band electron nuclear double resonance (ENDOR) and multifrequency electron paramagnetic resonance (EPR) (X, C, and S-band), all carried out under comparable cryogenic conditions, have provided markedly different electronic measures of the mutation-induced change. Q-band ENDOR of cysteine C(beta) protons, of weakly dipolar-coupled protons, and of the remaining His46 nitrogen ligand provided hyperfine couplings that were like those of other binuclear mixed-valence Cu(A) systems and were essentially unperturbed by the mutation at His120. The ENDOR findings imply that the Cu(A) core electronic structure remains unchanged by the His120 mutation. On the other hand, multifrequency EPR indicated that the H120N and H120G mutations had changed the EPR hyperfine signature from a 7-line to a 4-line pattern, consistent with trapped-valence, Type 1 mononuclear copper. The multifrequency EPR data imply that the electron spin had become localized on one copper by the His120 mutation. To reconcile the EPR and ENDOR findings for the His120 mutants requires that either: if valence localization to one copper has occurred, the spin density on the cysteine sulfurs and the remaining histidine (His46) must remain as it was for a delocalized binuclear Cu(A) center, or if valence delocalization persists, the hyperfine coupling for one copper must markedly diminish while the overall spin distribution on the CuSSCu core is preserved.     

Electron nuclear double resonance and electron paramagnetic resonance study on the structure of the NO-ligated heme alpha 3 in cytochrome c oxidase

The techniques of EPR and electron nuclear double resonance (ENDOR) were used to probe structure and electronic distribution at the nitric oxide (NO)-ligated heme alpha 3 in the nitrosylferrocytochrome alpha 3 moiety of fully reduced cytochrome c oxidase. Hyperfine and quadrupole couplings to NO (in both 15NO and 14NO forms), to histidine nitrogens, and to protons near the heme site were obtained. Parallel studies were also performed on NO-ligated myoglobin and model NO-heme-imidazole systems. The major findings and interpretations on nitrosylferrocytochrome alpha 3 were: 1) compared to other NO-heme-imidazole systems, the nitrosylferrocytochrome alpha3 gave better resolution of EPR and ENDOR signals; 2) at the maximal g value (gx = 2.09), particularly well resolved NO nitrogen hyperfine and quadrupole couplings and mesoproton hyperfine couplings were seen. These hyperfine and quadrupole couplings gave information on the electronic distribution on the NO, on the orientation of the g tensor with respect to the heme, and possibly on the orientation of the FeNO plane; 3) a combination of experimental EPR-ENDOR results and EPR spectral simulations evidenced a rotation of the NO hyperfine tensor with respect to the electronic g tensor; this implied a bent Fe-NO bond; 4) ENDOR showed a unique proton not seen in the other NO heme systems studied. The magnitude of this proton's hyperfine coupling was consistent with this proton being part of a nearby protein side chain that perturbs an axial ligand like NO or O2.     

The NMR structure of cyclosporin A bound to cyclophilin in aqueous solution

Cyclosporin A bound to the presumed receptor protein cyclophilin was studied in aqueous solution at pH 6.0 by nuclear magnetic resonance spectroscopy using uniform 15N- or 13C-labeling of cyclosporin A and heteronuclear spectral editing techniques. Sequence-specific assignments were obtained for all but one of the cyclosporin A proton resonances. With an input of 108 intramolecular NOEs and four vicinal 3JHN alpha coupling constants, the three-dimensional structure of cyclosporin A bound to cyclophilin was calculated with the distance geometry program DISMAN, and the structures resulting from 181 converged calculations were energy refined with the program FANTOM. A group of 120 conformers was selected on the basis of the residual constraint violations and energy criteria to represent the solution structure. The average of the pairwise root-mean-square distances calculated for the backbone atoms of the 120 structures was 0.58 A. The structure represents a novel conformation of cyclosporin A, for which the backbone conformation is significantly different from the previously reported structures in single crystals and in chloroform solution. The structure has all peptide bonds in the trans form, contains no elements of regular secondary structure and no intramolecular hydrogen bonds, and exposes nearly all polar groups to its environment. The root-mean-square distance between the backbone atoms of the crystal structure of cyclosporin A and the mean of the 120 conformers representing the NMR structure of cyclosporin A bound to cyclophilin is 2.5 A.     

EPR and ENDOR studies of the water oxidizing complex of Photosystem II

A comparative study of X-band EPR and ENDOR of the S₂ state of photosystem II membrane fragments and core complexes in the frozen state is presented. The S₂ state was generated either by continuous illumination at T=200 K or by a single turn-over light flash at T=273 K yielding entirely the same S₂ state EPR signals at 10 K. In membrane fragments and core complex preparations both the multiline and the g=4.1 signals were detected with comparable relative intensity. The absence of the 17 and 23 kDa proteins in the core complex preparation has no effect on the appearance of the EPR signals. ¹H-ENDOR experiments performed at two different field positions of the S₂ state multiline signal of core complexes permitted the resolution of four hyperfine (hf) splittings. The hf coupling constants obtained are 4.0, 2.3, 1.1 and 0.6 MHz, in good agreement with results that were previously reported (Tang et al. (1993) J Am Chem Soc 115: 2382–2389). The intensities of all four line pairs belonging to these hf couplings are diminished in D₂O. A novel model is presented and on the basis of the two largest hfc's distances between the manganese ions and the exchangeable protons are deduced. The interpretation of the ENDOR data indicates that these hf couplings might arise from water which is directly ligated to the manganese of the water oxidizing complex in redox state S₂.Abbreviations cw continuous wave - ENDOR electron nuclear double resonance - EPR electron paramagnetic resonance - hf hyperfine - hfc hyperfine coupling - MLS multiline signal - PS II Photosystem II - rf radio frequency - WOC water oxidizing complex     

Q-band ENDOR (electron nuclear double resonance) of the high-affinity ubisemiquinone center in cytochrome bo3 from Escherichia coli

Electron nuclear double resonance (ENDOR) was performed on the protein-bound, stabilized, high-affinity ubisemiquinone radical, QH*-, of bo3 quinol oxidase to determine its electronic spin distribution and to probe its interaction with its surroundings. Until this present work, such ENDOR studies of protein-stabilized ubisemiquinone centers have only been done on photosynthetic reaction centers whose function is to reduce a ubiquinol pool. In contrast, QH*- serves to oxidize a ubiquinol pool in the course of electron transfer from the ubiquinol pool to the oxygen-consuming center of terminal bo3 oxidase. As documented by large hyperfine couplings (>10 MHz) to nonexchangeable protons on the QH*- ubisemiquinone ring, we provide evidence for an electronic distribution on QH*- that is different from that of the semiquinones of reaction centers. Since the ubisemiquinone itself is physically nearly identical in both QH*- and the bacterial photosynthetic reaction centers, this electronic difference is evidently a function of the local protein environment. Interaction of QH*- with this local protein environment was explicitly shown by exchangeable deuteron ENDOR that implied hydrogen bonding to the quinone and by weak proton hyperfine couplings to the local protein matrix.     

Probing magnetic properties of the reduced [2Fe-2S] cluster of the ferredoxin from Arthrospira platensis by 1H ENDOR spectroscopy

The 1H electron nuclear double resonance (ENDOR) spectra in frozen solutions of the reduced [2Fe-2S] cluster in ferredoxin from Arthrospira (Spirulina) platensis have been measured at low temperatures (5-20 K) and simulated using orientational selection methods. The analysis confirmed the existence of a single paramagnetic species with iron valence states II and III connected uniquely to the cluster irons. The experimental ENDOR spectra were fitted to a model including the spin distribution on the centre, the orientation of the g-matrix, and the isotropic and anisotropic hyperfine couplings of the nearest protons in the crystallographically determined structure. In order to partially simulate ENDOR line shapes, a statistical distribution of the corresponding torsion angles between the Fe(III) centre and one of the beta-CH2 protons was introduced. From the analysis, four of the larger hyperfine couplings found were assigned to the cysteine beta-protons near the Fe(III) ion of the cluster, with isotropic hyperfine couplings ranging from 1.6 to 4.1 MHz. The spin distribution on the two iron ions was estimated to be +1.85 for the Fe(III) ion and -0.9 for the Fe(II) ion. The Fe(III) ion was identified as being coordinated to the cysteine ligands Cys49 and Cys79, confirming previous NMR results. The direction of the g-tensor with respect to the cluster was deduced. The g1-g2 plane is parallel to the planes through each iron and its adjacent cysteine sulfurs; the g2-g3 plane is nearly perpendicular to the latter planes and deviates by 25 degrees from the FeSSFe plane. The g1 direction is dominated by the bonding geometry of Fe(II) and does not align with the Fe(II)-Fe(III) vector.     

Crystallographic structure refinement of Chromatium high potential iron protein at two Angstroms resolution.

The structure of Chromatium high potential iron protein (HiPIP) has been refined by semiautomatic Fo-Fc (observed minus calculated structure amplitude Fourier methods to a convential R index, R=sum of the absolute value of Fo-Fc divided by the sum of Fo, of 24.7% for a model in which bond distances and angles are constrained to standard values. Bond length and angle constraints were applied only intermittenly during the computations. At a late stage of the refinement, atomic parameters for only the Fe4S4 cluster plus the 4 associated cystein S-gamma atoms were adjusted by least squares methods and kept fixed during the rest of the refinement. The refined model consists of 625 of the 632 nonhydrogen atoms in the protein plus 75 water molecules. Seven side chain atoms could not be located in the final electron density map. A computer program rather than visual inspection was used wherever possible in the refinement: for locating water molecules, for removing water molecules that too closely approach other atoms, for deleting atoms that lay in regions of low electron density, and for evaluating the progress of refinement. Fo-Fc Fourier refinement is sufficiently economical to be applied routinely in protein crystal structure determinations. The complete HiPIP refinement required approximately 12 hours of CDC 3600 computer time and cost less than $3000 starting from a "trial structure," based upon multipe isomorphoous replacement phases, which gave an R of 43%...     

Electron spin echo envelope modulation spectroscopy in photosystem I.

The applications of electron spin echo envelope modulation (ESEEM) spectroscopy to study paramagnetic centers in photosystem I (PSI) are reviewed with special attention to the novel spectroscopic techniques applied and the structural information obtained. We briefly summarize the physical principles and experimental techniques of ESEEM, the spectral shapes and the methods for their analysis. In PSI, ESEEM spectroscopy has been used to the study of the cation radical form of the primary electron donor chlorophyll species, P(700)(+), and the phyllosemiquinone anion radical, A(1)(-), that acts as a low-potential electron carrier. For P(700)(+), ESEEM has contributed to a debate concerning whether the cation is localized on a one or two chlorophyll molecules. This debate is treated in detail and relevant data from other methods, particularly electron nuclear double resonance (ENDOR), are also discussed. It is concluded that the ESEEM and ENDOR data can be explained in terms of five distinct nitrogen couplings, four from the tetrapyrrole ring and a fifth from an axial ligand. Thus the ENDOR and ESEEM data can be fully accounted for based on the spin density being localized on a single chlorophyll molecule. This does not eliminate the possibility that some of the unpaired spin is shared with the other chlorophyll of P(700)(+); so far, however, no unambiguous evidence has been obtained from these electron paramagnetic resonance methods. The ESEEM of the phyllosemiquinone radical A(1)(-) provided the first evidence for a tryptophan molecule pi-stacked over the semiquinone and for a weaker interaction from an additional nitrogen nucleus. Recent site-directed mutagenesis studies verified the presence of the tryptophan close to A(1), while the recent crystal structure showed that the tryptophan was indeed pi-stacked and that a weak potential H-bond from an amide backbone to one of the (semi)quinone carbonyls is probably the origin of the to the second nitrogen coupling seen in the ESEEM. ESEEM has already played an important role in the structural characterization on PSI and since it specifically probes the radical forms of the chromophores and their protein environment, the information obtained is complimentary to the crystallography. ESEEM then will continue to provide structural information that is often unavailable using other methods.     

Compound I radical in site-directed mutants of cytochrome c peroxidase as probed by electron paramagnetic resonance and electron-nuclear double resonance.

The reaction of ferric cytochrome c peroxidase (CcP) from Saccharomyces cerevisiae with peroxide produces compound I, characterized by both an oxyferryl iron center and a protein-based free radical. The electron paramagnetic resonance (EPR) signal of the CcP compound I radical can be resolved into a broad majority component which accounts for approximately 90% of the spin intensity and a narrow minority component which accounts for approximately 10% of the integrated spin intensity [Hori, H., & Yonetani, T. (1985) J. Biol. Chem. 260, 3549-3555]. It was shown previously that the broad component of the compound I radical signal is eliminated by mutation of Trp-191 to Phe [Scholes, C. P., Liu, Y., Fishel, L. F., Farnum, M. F., Mauro, J. M., & Kraut, J. (1989) Isr. J. Chem. 29, 85-92]. The present work probed the effect of mutations in the vicinity of this residue by EPR and electron-nuclear double resonance (ENDOR). These mutations were obtained from a plasmid-encoded form of S. cerevisiae expressed in Escherichia coli [Fishel, L. A., Villafranca, J. E., Mauro, J. M., & Kraut, J. (1987) Biochemistry 26, 351-360]. The EPR line shape and ENDOR signals of the compound I radical were perturbed only by mutations that alter Trp-191 or residues in its immediate vicinity: namely, Met-230 and Met-231, which have sulfur atoms within 4 A of the indole ring, and Asp-235, which forms a hydrogen bond with the indole nitrogen of Trp-191. Mutations of other potential oxidizable sites (tryptophan, tyrosine, methionine, and cysteine) did not alter the EPR line shapes of the compound I radical, although the integrated spin intensities were weaker in some of these mutants. Mutations at Met-230 and/or -231 perturbed the EPR line shapes of the compound I radical signal but did not eliminate it. ENDOR of these two methionine mutants showed alteration to the hyperfine couplings of several strongly coupled protons, which are characteristic of the majority compound I radical electronic structure, and a change in weaker hyperfine couplings, which suggests a different orientation of the radical with respect to its surroundings in the presence of these methionine mutations. Besides the Trp-191----Phe mutation, only the Asp-235----Asn mutation eliminated the broad component of the compound I signal. Loss of the broad compound I EPR signal coincides with both the loss of the Asp----Trp-191 hydrogen-bonding interaction and alteration of the position of the indole ring of Trp-191.(ABSTRACT TRUNCATED AT 400 WORDS)     

Density functional calculated spin densities and hyperfine couplings for hydrogen bonded 1,4-naphthosemiquinone and phyllosemiquinone anion radicals: a model for the A1 free radical formed in photosystem I

The B3LYP hybrid density functional method is used to calculate spin densities and hyperfine couplings for the 1,4-naphthosemiquinone anion radical and a model of the phyllosemiquinone anion radical. The effect of hydrogen bonding on the spin density distribution is shown to lead to a redistribution of pi spin density from the semiquinone carbonyl oxygens to the carbonyl carbon atoms. The effect of in plane and out of plane hydrogen bonding is examined. Out of plane hydrogen bonding is shown to give rise to a significant delocalisation of spin density on to the hydrogen bond donor heavy atom. Excellent agreement is observed between calculated and experimental hyperfine couplings. Comparison of calculated hyperfine couplings with experimental determinations for the A1 phyllosemiquinone anion radical present in Photosystem I (PS I) of higher plant photosynthesis indicates that the in vivo radical may have a hydrogen bond to the O4 atom only as opposed to hydrogen bonds to each oxygen atom in alcohol solvents. The hydrogen bonding situation appears to be the reverse of that observed for QA in the bacterial type II reaction centres where the strong hydrogen bond occurs to the quinone O1 oxygen atom. For different types of reaction centre the presence or absence of the non-heme Fe(II) atom may well determine which type of hydrogen bonding situation prevails at the primary quinone site which in turn may influence the direction of subsequent electron transfer.     

The metal-binding site of imidazole glycerol phosphate dehydratase; EPR and ENDOR studies of the oxo-vanadyl enzyme

 The apo protein of imidazole glycerol phosphate dehydratase (IGPD) from Saccharomyces cerevisiae combines stoichiometrically with certain specific divalent metal cations to assemble the catalytically active form comprising 24 protein subunits and tightly bound metal. VO²⁺ ions react similarly but, uniquely, result in a metallo-protein (VO-IGPD) with neither catalytic activity nor the ability to bind to the reaction intermediate analogue, 2-hydroxy-3-(1,2,4-triazol-1-yl) propylphosphonate. Since VO²⁺ apparently assembles the quaternary structure correctly, it is used in the present study as a spin probe to investigate the metal centre coordination environment by electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) spectroscopy. At neutral pH, the EPR spectrum of VO-IGPD reveals at least three distinct VO²⁺ sub-spectra with one predominant at low pH. The spin Hamiltonian parameters for some of the sub-spectra are consistent with ⁵¹V having nitrogen in the inner-sphere equatorial coordination environment from, most probably, multiple coordinating histidines. Further evidence for inner-sphere nitrogen ligands is obtained from ENDOR spectroscopy. The spectra of the low rf region show signals from interactions with ¹⁴N which are consistent with couplings to the imino nitrogen of coordinated histidine residues. In addition a number of proton ENDOR line pairs are resolved. Of the few that disappear upon exchange of the protein into D₂O, one most likely originates from the exchangeable proton of the N-H group of a coordinated histidine imidazole. ¹H-ENDOR line pairs from non-exchangeable protons with splittings of approximately 3 MHz can be attributed to imidazole carbon protons. Thus, most of the couplings observed by ENDOR are consistent with being from the imidazole heterocycle of one or more histidine ligands. Received: 27 June 1996 / Accepted: 14 March 1997     

The electrostatic potential for the phosphodiester group determined from X-ray diffraction.

The charge density distribution in the crystal structure of ammonium dimethylphosphate at 123 K has been determined from x-ray diffraction data (MoK alpha) using 8437 reflections with sin theta/lambda less than 1.33 A-1 [NH4+.(CH3)2PO4-, M(r) = 143.08, monoclinic, P2(1)/c, a = 10.007(1), b = 6.926(1), c = 9.599(2) A, beta = 105.40(1) degrees, V = 641.4(3) A3, Z = 4, F000 = 304, Dx = 1.4815 g.cm-3, mu = 3.726 cm-1]. Least-squares structure refinement assuming Stewart's rigid pseudoatom model (variables including Slater-type radial exponents and electron populations for multipole terms extending to octapoles for C, N, O, and P, and dipoles for H) gave R(F2) = 0.039 for all reflections. The dimethylphosphate anion is in the gauche-gauche conformation and has approximate twofold symmetry. One phosphoryl O atom forms three hydrogen bonds and the other forms one. Neither of the ester O atoms is hydrogen bonded. For the dimethylphosphate anion isolated from the crystal structure, a map of the electrostatic potential obtained using the pseudoatom charge parameters shows that the phosphoryl O atoms are considerably more electronegative than the ester O atoms. The electrostatic potential distribution obtained in this way has been fitted by least squares to a system of atom-centered point charges. The potential calculated from these point charges agrees with the experimental result. It also agrees reasonably well with potentials obtained from three other systems of point charges that are widely used as part of the semiempirical force field for molecular mechanics and molecular dynamics calculations involving nucleic acids.     

Structure of scorpion toxin variant-3 at 1.2 A resolution.

The crystal structure of the variant-3 protein neurotoxin from the scorpion Centruroides sculpturatus Ewing has been refined at 1.2 A resolution using restrained least-squares. The final model includes 492 non-hydrogen protein atoms, 453 protein hydrogen atoms, eight 2-methyl-2,4-pentanediol (MPD) solvent atoms, and 125 water oxygen atoms. The variant-3 protein model geometry deviates from ideal bond lengths by 0.024 A and from ideal angles by 3.6 degrees. The crystallographic R-factor for structure factors calculated from the final model is 0.192 for 17,706 unique reflections between 10.0 to 1.2 A. A comparison between the models of the initial 1.8 A and the 1.2 A refinement shows a new arrangement of the previously poorly defined residues 31 to 34. Multiple conformations are observed for four cysteine residues and an MPD oxygen atom. The electron density indicates that disulfide bonds between Cys12 and Cys65 and between Cys29 and Cys48 have two distinct side-chain conformations. A molecule of MPD bridges neighboring protein molecules in the crystal lattice, and both MPD enantiomers are present in the crystal. A total of 125 water molecules per molecule of protein are included in the final model with B-values ranging from 11 to 52 A2 and occupancies from unity down to 0.4. Comparisons between the 1.2 A and 1.8 A models, including the bound water structure and crystal packing contacts, are emphasized.     

Electron nuclear double resonance of cytochrome oxidase: nitrogen and proton ENDOR from the 'copper' EPR signal.

Proton ENDOR resonances have been found from at least two different protons with fairly large and isotropic couplings of about 12 and 19 MHz. It is possible that such protons are attached to carbons that are one bond removed from the point of ligation to copper. A number of weakly coupled protons with anisotropic couplings have also been seen. None of the protons, either weakly or strongly coupled, appears to exchange with 2H2O. We have obtained nitrogen ENDOR from at least one nitrogen with a hyperfine coupling large enough for the nitrogen to be a ligand of copper. We have not yet demonstrated experimentally ENDOR characteristic of the copper nucleus itself.     

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.