A study of vibrational relaxation of B-state carbon monoxide in the heme pocket of photolyzed carboxymyoglobin.

The vibrational energy relaxation of dissociated carbon monoxide in the heme pocket of sperm whale myoglobin has been studied using equilibrium molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the delta- and epsilon-tautomers of the distal histidine, His64. Vibrational population relaxation times were estimated using the Landau-Teller model. For carbon monoxide (CO) in the myoglobin epsilon-tautomer, for a frequency of omega0 = 2131 cm-1 corresponding to the B1 state, T1epsilon(B1) = 640 +/- 185 ps, and for a frequency of omega0 = 2119 cm-1 corresponding to the B2 state, T1epsilon(B2) = 590 +/- 175 ps. Although the CO relaxation rates in both the epsilon- and delta-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the delta-tautomer. For CO in the myoglobin delta-tautomer, it was found that the relaxation times were identical within error for the two CO substate frequencies, T1delta(B1) = 335 +/- 115 ps and T1delta(B2) = 330 +/- 145 ps. These simulation results are in reasonable agreement with experimental results of Anfinrud and coworkers (unpublished results). Normal mode calculations were used to identify the dominant coupling between the protein and CO molecules. The calculations suggest that the residues of the myoglobin pocket, acting as a first solvation shell to the CO molecule, contribute the primary "doorway" modes in the vibrational relaxation of the oscillator.     

Raman spectra and vibrational assignments for deuterated membrane lipids. 1,2-Dipalmitoyl phosphatidylcholine-d9 and -d62

Vibrational Raman spectra of polycrystalline 1,2-dipalmitoyl phosphatidylcholine-d9 (fully deuterated choline methyl groups) and 1,2-dipalmitoyl phosphatidylcholine-d62 (fully deuterated acyl chains) were recorded in the 3050- 2800, 2250-2050 and 1800-700 cm-1 regions. The fundamental vibrational modes were assigned primarily on the basis of isotopic frequency shift ratios, group frequency correlations and comparisons with specific model compounds. Since deuterium-substituted lipids provide well-isolated spectral probes, particularly in the carbon-deuterium stretching region, the dependence of the 2250-2050 cm-1 region on lipid phase was examined for the dipalmitoyl phosphatidylcholine-d62 species. The methylene CD2 deformation and twisting modes at 984 and 919 cm-1, respectively, also exhibit intense, isolated vibrational transitions which should prove useful for monitoring molecular order in mixed dueterated and undeuterated lipid systems. Except for the relatively weak choline methyl C-D and C-H stretching modes, the spectrum of 1,2-dipalmitoyl phosphatidylcholine-d9 is not distinguishable from that of the undeuterated system. For both the d9 and undeuterated species, the vibrational modes associated with the lipid head group region are sensitive to slight hydration.     

Dynamics of vibrational frequency fluctuations in deuterated liquid ammonia: roles of fluctuating hydrogen bonds and free ND modes

We present an ab initio molecular dynamics study of the roles of fluctuating hydrogen bonds and free ND modes in the dynamics of ND stretch frequency fluctuations in deuterated liquid ammonia. We have also looked at some of the other dynamical quantities such as diffusion and orientational relaxation and also structural quantities such as pair correlations and hydrogen bonding properties which are relevant in the current context. The time correlation function of ND stretch frequencies is found to decay with primarily two time scales: A short-time decay with a time scale of less than 100 fs arising from intermolecular motion of intact hydrogen bonds and also from fast hydrogen bond breaking and a longer time scale of about 500 fs which can be assigned to the lifetime of free ND modes. Unlike water, in liquid ammonia an ND mode is found to remain free for a longer period than it stays hydrogen bonded and this longer lifetime of free ND modes determines the long-time behaviour of frequency fluctuations. Our hole dynamics calculations produced results of vibrational spectral diffusion that are similar to the decay of frequency time correlation. Inclusion of dispersion corrections is found to make the dynamics slightly faster.     

Optical and resonance Raman spectroscopy of the heme groups of the quinol-oxidizing cytochrome aa3 of Bacillus subtilis.

The cytochrome aa3-type terminal quinol oxidase of Bacillus subtilis catalyzes the four-electron reduction of dioxygen to water. It resembles the aa3-type cytochrome-c oxidase in using heme A as its active-site chromophores but lacks the CuA center and the cytochrome-c oxidizing activity of the mitochondrial enzyme. We have used optical and resonance Raman spectroscopies to study the B. subtilis oxidase in detail. The alpha-band absorption maximum of the reduced minus oxidized enzyme is shifted by 5-7 nm to the blue relative to most other aa3-type oxidases, and accordingly, we designate the Bacillus enzyme as cytochrome aa3-600. The shifted optical spectrum cannot be ascribed to an alteration in the strength of the hydrogen bond between the formyl group of the low-spin heme and its environment, as the Raman line assigned to this mode in aa3-600 has the same frequency and degree of resonance enhancement as the low-spin heme a formyl mode in most other aa3-type oxidases. Raman modes arise at 194 and 214 cm-1 in aa3-600, whereas a single band at about 214 cm-1 is assigned to the iron-histidine stretch for the other aa3-type oxidases. Possible explanations for the occurrence of these two modes are discussed. Comparison of formyl and vinyl modes and heme skeletal vibrational modes in different oxidation states of aa3-600 and of beef heart cytochrome-c oxidase shows a strong similarity, which suggests conservation of essential features of the heme environments in these oxidases.     

Cytochrome c oxidase: evidence for interaction of water molecules with cytochrome a

The resonance Raman spectra of cytochrome c oxidase in protonated buffer compared to that in deuterated buffer indicate that water molecules are near the heme of cytochrome a. Differences in widths of the heme line at 1610 cm-1, after short exposure to D2O, and, additionally, of the heme line at 1625 cm-1, after long exposure, can be accounted for by changes in resonance vibrational energy transfer between modes of cytochrome a2+ and the bending mode of water molecules in the heme pocket. On the basis of the assignment of these modes, we place one water molecule near the vinyl group and one water molecule near the formyl group of the cytochrome a heme. These water molecules may play several possible functional roles.     

FTIR spectroscopic evidence for the involvement of an acidic residue in quinone binding in cytochrome bd from Escherichia coli

In this work, FTIR difference spectroscopy is used to search for possible binding partners and protonable groups involved in the binding of the quinol to cytochrome bd from Escherichia coli. In addition, the electrochemically induced FTIR difference spectra are compared for preparations of the enzyme isolated from cells grown at different oxygen levels in which the quinone content of the membrane is altered. On this basis, difference signals can be tentatively attributed to the vibrational modes of the different quinones types that are associated with the enzyme depending on growth conditions. Furthermore, vibrational modes due to the redox-dependent reorganization of the protein vary depending on the quinone associated with the isolated enzyme. Of particular interest are the observations that a mode at 1738 cm(-1) is decreased and a mode at 1595 cm(-1) is increased as observed in direct comparison to the data obtained from samples grown anaerobically. These signals indicate a change in the protonation state of an aspartic or glutamic acid. Since these changes are observed when the ubiquinone ratio in the preparation increases, the data provide evidence for the modulation of the binding site by the interacting quinone and the involvement of an acidic group in the binding site. The tentative assignments of the vibrational modes are supported by electrochemically induced FTIR difference spectra of cytochrome bd in the presence of the specific quinone binding site inhibitors heptylhydroxyquinoline-N-oxide (HQNO) or 2-methyl-3-undecylquinolone-4. Whereas HQNO leads to strong shifts in the FTIR redox difference spectrum, 2-methyl-3-undecylquinolone-4 induces a specific shift of a mode at 1635 cm(-1), which likely originates from the displacement of the C=O group of the bound quinone.     

Drug-excipient interactions in ketoprofen: a vibrational spectroscopy study

Ketoprofen (3-benzoyl-alpha-methylbenzeneacetic acid) is a widely used nonsteroidal anti-inflammatory drug (NSAID), always administered in the form of drug-excipient physical mixtures (PMs). The occurrence of possible interactions between ketoprofen and two commonly used excipients-lactose (LAC) and polyvinylpyrrolidone (PVP)-was evaluated, through vibrational spectroscopy techniques [both Raman and Inelastic Neutron Scattering (INS)]. Spectral evidence of drug:excipient close contacts, which were enhanced by aging, was verified for the (1:1) (w:w) (ketoprofen:PVP) and (ketoprofen:LAC) PMs, both by Raman and INS. These interactions were found to involve mainly the central carbonyl and the terminal methyl-carboxylic moieties of the ketoprofen molecule, this being reflected in particular vibrational modes, such as the methyl torsion, the out-of-plane C-OH bending, and the inter-ring C=O stretching.     

Picosecond dynamics of the glutamate receptor in response to agonist-induced vibrational excitation

Conformational changes of proteins are dominated by the excitation and relaxation processes of their vibrational states. To elucidate the mechanism of receptor activation, the conformation dynamics of receptors must be analyzed in response to agonist-induced vibrational excitation. In this study, we chose the bending vibrational mode of the guanidinium group of Arg485 of the glutamate receptor subunit GluR2 based on our previous studies, and we investigated picosecond dynamics of the glutamate receptor caused by the vibrational excitation of Arg485 via molecular dynamics simulations. The vibrational excitation energy in Arg485 in the ligand-binding site initially flowed into Lys730, and then into the J-helix at the subunit interface of the ligand-binding domain. Consequently, the atomic displacement in the subunit interface around an intersubunit hydrogen bond was evoked in about 3 ps. This atomic displacement may perturb the subunit packing of the receptor, triggering receptor activation.     

Deep-UV resonance Raman analysis of the Rhodobacter capsulatus cytochrome bc₁complex reveals a potential marker for the transmembrane peptide backbone

Classical strategies for structure analysis of proteins interacting with a lipid phase typically correlate ensemble secondary structure content measurements with changes in the spectroscopic responses of localized aromatic residues or reporter molecules to map regional solvent environments. Deep-UV resonance Raman (DUVRR) spectroscopy probes the vibrational modes of the peptide backbone itself, is very sensitive to the ensemble secondary structures of a protein, and has been shown to be sensitive to the extent of solvent interaction with the peptide backbone [ Wang , Y. , Purrello , R. , Georgiou , S. , and Spiro , T. G. ( 1991 ) J. Am. Chem. Soc. 113 , 6368 - 6377 ]. Here we show that a large detergent solubilized membrane protein, the Rhodobacter capsulatus cytochrome bc(1) complex, has a distinct DUVRR spectrum versus that of an aqueous soluble protein with similar overall secondary structure content. Cross-section calculations of the amide vibrational modes indicate that the peptide backbone carbonyl stretching modes differ dramatically between these two proteins. Deuterium exchange experiments probing solvent accessibility confirm that the contribution of the backbone vibrational mode differences are derived from the lipid solubilized or transmembrane α-helical portion of the protein complex. These findings indicate that DUVRR is sensitive to both the hydration status of a protein's peptide backbone, regardless of primary sequence, and its secondary structure content. Therefore, DUVRR may be capable of simultaneously measuring protein dynamics and relative water/lipid solvation of the protein.     

NMR-derived models of amidopyrine and its metabolites in complexes with rabbit cytochrome P450 2B4 reveal a structural mechanism of sequential N-dealkylation

To understand the molecular basis of sequential N-dealkylation by cytochrome P450 2B enzymes, we studied the binding of amidopyrine (AP) as well as the metabolites of this reaction, desmethylamidopyrine (DMAP) and aminoantipyrine (AAP), using the X-ray crystal structure of rabbit P450 2B4 and two nuclear magnetic resonance (NMR) techniques: saturation transfer difference (STD) spectroscopy and longitudinal (T(1)) relaxation NMR. Results of STD NMR of AP and its metabolites bound to P450 2B4 were similar, suggesting that they occupy similar niches within the enzyme's active site. The model-dependent relaxation rates (R(M)) determined from T(1) relaxation NMR of AP and DMAP suggest that the N-linked methyl is closest to the heme. To determine the orientation(s) of AP and its metabolites within the P450 2B4 active site, we used distances calculated from the relaxation rates to constrain the metabolites to the X-ray crystal structure of P450 2B4. Simulated annealing of the complex revealed that the metabolites do indeed occupy similar hydrophobic pockets within the active site, while the N-linked methyls are free to rotate between two binding modes. From these bound structures, a model of N-demethylation in which the N-linked methyl functional groups rotate between catalytic and noncatalytic positions was developed. This study is the first to provide a structural model of a drug and its metabolites complexed to a cytochrome P450 based on NMR and to provide a structural mechanism for how a drug can undergo sequential oxidations without unbinding. The rotation of the amide functional group might represent a common structural mechanism for N-dealkylation reactions for other drugs such as the local anesthetic lidocaine.     

Nuclear resonance vibrational spectroscopy--NRVS

The recent, synchrotron-based vibrational technique nuclear resonance vibrational spectroscopy (NRVS) is introduced. The method can be used for a number of M?ssbauer active isotopes including 57Fe, which has yielded most of the results to date. The NRVS experiment can be thought of as M?ssbauer spectroscopy with vibrational sidebands. Importantly, the NRVS experiment provides the complete set of bands corresponding to modes that involve motion of the iron atom. The method has a selectivity reminiscent of that of resonance Raman spectroscopy, but with the significant advantage that NRVS is not subject to the optical selection rules of Raman or infrared spectroscopy. Indeed, NRVS provides the ultimate limit in selectivity because only the vibrational dynamics of the probe nucleus contribute to the observed signal. All iron-ligand modes will be observed, including many that had not been previously observed. For hemes, these include in-plane iron vibrations that have not yet been reported by resonance Raman studies and the iron-imidazole stretch that has not been identified in six-coordinate porphyrins. Other modes that can be investigated include that of heme doming that is expected to be a low-frequency mode. The experimental setup at a beam line and sample requirements for iron-based derivatives are presented. Both powder and polarized single-crystal measurements can be made. The general features of data extraction and analysis are given. Data for heme and heme proteins are given. Examples of assignment of spectra for nitrosyl and carbonyl derivatives are given. These data demonstrate the importance of peripheral substituents on the vibrational spectrum of heme derivatives. Delocalization of modes appears to be common. Although this technique has only been available for a relatively short time, this early progress report indicates that NRVS has significant potential for probing the dynamics of Fe-containing molecules of biological interest.     

A vibrational study of the CD2 stretching bands of selectively deuterated palmitic and stearic acids

The C-D stretching regions of the infrared and Raman spectra of 14 different palmitic and stearic acids containing isolated CD₂ groups are reported. Anomalous behaviour is observed when substitution occurs near the terminal methyl group, which behaviour can not be explained in terms of crystal field effects.     

Relaxation dynamics of a photoinduced di-cobalt-tetraoxolene valence tautomer

The temperature dependence of the relaxation dynamics of a dinuclear cobalt-tetraoxolene complex undergoing photoinduced valence tautomerism is reported and discussed. Two different regimes were detected, one at low (5-20 K) and a second at higher temperatures (20-35 K). The activation energy for the low temperature regime suggests a tunnelling process assisted by lattice phonons. In the thermally activated region the activation energy is far from that expected for the total-symmetric Co-O breathing mode, thus suggesting an important role for other internal vibrational modes.     

Mechanism of unfolding of a model helical peptide

Synthetic model helical peptides, Acetyl-W(EAAAR)(5)A-amide with (13)C=O specifically labeled alanine segments in repeats n = 1,2 or 4,5 were studied in aqueous D(2)O solution as a function of temperature using Fourier transform infrared spectroscopy and two-dimensional correlation analysis. The (13)C==O provided a probe which was sensitive to the carbonyl stretch in the peptide bonds of the alanine residues at the amino terminal end in one peptide as compared to the probe in the carboxy terminal end of the other peptide during thermal perturbation. The relative stability of each terminal end was examined; the more stable terminal was determined to be the amino terminal end. Also studied were the glutamate and arginine side-chain modes involved in the salt bridging interaction. Two-dimensional correlation analysis enabled enhanced resolution in the spectral region of 1520--1700 cm(-1), and thus, the order in which these vibrational modes were perturbed as a function of increasing temperature were established.     

Evidence of protein collective motions on the picosecond timescale

We investigate the presence of structural collective motions on a picosecond timescale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations, and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ∼25% (g H₂O/g protein). Quasiharmonic vibrational modes and dipole-dipole correlation functions were calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence, providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of buried internal water molecules. This source for the observed oxidation dependence is consistent with the lack of an oxidation dependence in nuclear resonant vibrational spectroscopy measurements.     

Quantum chemical study of agonist-receptor vibrational interactions for activation of the glutamate receptor

To understand the mechanism of activation of a receptor by its agonist, the excitation and relaxation processes of the vibrational states of the receptor should be examined. As a first approach to this problem, we calculated the normal vibrational modes of agonists (glutamate and kainate) and an antagonist (6-cyano-7-nitroquinoxaline-2,3-dione: CNQX) of the glutamate receptor, and then investigated the vibrational interactions between kainate and the binding site of glutamate receptor subunit GluR2 by use of a semiempirical molecular orbital method (MOPAC2000-PM3). We found that two local vibrational modes of kainate, which were also observed in glutamate but not in CNQX, interacted through hydrogen bonds with the vibrational modes of GluR2: (i) the bending vibration of the amine group of kainate, interacting with the stretching vibration of the carboxyl group of Glu705 of GluR2, and (ii) the symmetric stretching vibration of the carboxyl group of kainate, interacting with the bending vibration of the guanidinium group of Arg485. We also found collective modes with low frequency at the binding site of GluR2 in the kainate-bound state. The vibrational energy supplied by an agonist may flow from the high-frequency local modes to the low-frequency collective modes in a receptor, resulting in receptor activation.     

Fast structural dynamics in reduced and oxidized cytochrome c

The sub‐nanosecond structural dynamics of reduced and oxidized cytochrome c were characterized. Dynamic properties of the protein backbone measured by amide ¹⁵N relaxation and side chains measured by the deuterium relaxation of methyl groups change little upon change in the redox state. These results imply that the solvent reorganization energy associated with electron transfer is small, consistent with previous theoretical analyses. The relative rigidity of both redox states also implies that dynamic relief of destructive electron transfer pathway interference is not operational in free cytochrome c.     

Nonlinear electron paramagnetic resonance studies of the interaction of cytochrome c oxidase with spin-labeled lipids in gel-phase membranes

The interaction of lipids, spin-labeled at different positions in the sn-2 chain, with cytochrome c oxidase reconstituted in gel-phase membranes of dimyristoylphosphatidylglycerol has been studied by electron paramagnetic resonance (EPR) spectroscopy. Nonlinear EPR methods, both saturation transfer EPR and progressive saturation EPR, were used. Interaction with the protein largely removes the flexibility gradient of the lipid chains in gel-phase membranes. The rotational mobility of the chain segments is reduced, relative to that for gel-phase lipids, by the intramembranous interaction with cytochrome c oxidase. This holds for all positions of chain labeling, but the relative effect is greater for chain segments closer to the terminal methyl ends. Modification of the paramagnetic metal-ion centers in the protein by binding azide has a pronounced effect on the spin-lattice relaxation of the lipid spin labels. This demonstrates that the centers modified are sufficiently close to the first-shell lipids to give appreciable dipolar interactions and that their vertical location in the membrane is closer to the 5-position than to the 14-position of the lipid chains.     

Cytochrome c reconstituted from two peptide fragments displays native-like redox properties.

Recombination of two fragments of horse cytochrome c (the heme-containing N-fragment, residues 1-56, and the C-fragment, residues 57-104), which are substantially unstructured at neutral pH, gives rise to a 1:1 fragment complex with a compact conformation, in which the alpha helical structure and the native Met80-Fe(III) axial bond are recovered. With respect to the native protein, the ferric complex shows a less rigid atomic packing and a decreased stability [Delta(DeltaG(o))D = 14.7 kJ.mol(-1)], ascribed to perturbations involving the Trp59 microenvironment and, to a lower extent, the heme pocket region. The redox potential, E1/2 = 234 +/- 5 mV vs. normal hydrogen electrode at 25 degrees C, is close to that of the intact protein, consistent with recovery of the native Met80-heme Fe(III) axial bond. Furthermore, the fragment complex shows reactivity similar to intact cytochrome c, in the reaction with cytochrome c oxidase. We conclude that the absence in the complex of some native cross-links and interlocked packing important for protein rigidity and stability is not as relevant for maintaining the native redox properties of the protein, provided that some structural requirements (i.e. recovering of the native-like alpha helical structure) are fulfilled and coordination of Met80 to the heme-iron is restored.     

N-Isotope effects on the Raman spectra of Fe2S2 ferredoxin and Rieske ferredoxin: evidence for structural rigidity of metal sites

The diiron ferredoxins have a common diamond-core structure with two bridging sulfides, but differ in the nature of their terminal ligands: either four cysteine thiolates in the Fe(2)S(2) ferredoxins or two cysteine thiolates and two histidine imidazoles in the Rieske ferredoxins. Contributions of the bridging (b) and terminal (t) ligands to the resonance Raman spectra of the Fe(2)S(2) ferredoxins have been distinguished previously by isotopic substitution of the bridging sulfides. We now find that uniform (15)N-labeling of Anabaena Fe(2)S(2) ferredoxin results in shifts of -1 cm(-1) in the Fe-S(t) stretching modes at 282, 340, and 357 cm(-1). The (15)N dependence is ascribed to kinematic coupling of the Fe-S(Cys) stretch with deformations of the cysteine backbone, including the amide nitrogen. No (15)N dependence occurs for the nu(Fe-S(b)) modes at 395 and 426 cm(-1). Similar effects are observed for the Rieske center in T4MOC ferredoxin from the toluene-4-monooxygenase system of Pseudomonas mendocina. Upon selective (15)N-labeling of the alpha-amino group of cysteine, the vibrational modes at 321, 332, 350, and 362 cm(-1) all undergo shifts of -1 to -2 cm(-1), thereby identifying them as combinations of nu(Fe-S(t)) and delta(Cys). These same four modes undergo similar isotope shifts when T4MOC ferredoxin is selectively labeled with (15)N-histidine ((15)N in either the alpha1,delta1 or delta1,epsilon2 positions). Thus, the Fe-S(Cys) stretch must also be undergoing kinematic coupling with vibrations of the Fe-His moiety. The extensive kinematic coupling of iron ligand vibrations observed in both the Fe(2)S(2) and Rieske ferredoxins presumably arises from the rigidity of the protein framework and is reminiscent of the behavior of cupredoxins. In both cases, the structural rigidity is likely to play a role in minimizing the reorganization energy for electron transfer.