Spectroscopic evidence for conformational relaxation in myoglobin

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Rebinding dynamics of NO to microperoxidase-8 probed by time-resolved vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding kinetics of NO to microperoxidase-8 (Mp), an ideal model system for the active site of ligand-binding heme proteins, including myoglobin and hemoglobin, after the photodeligation of MpNO in glycerol/water (G/W) solutions at 294 K. The geminate rebinding (GR) of NO to Mp in viscous solutions was highly efficient and ultrafast and negligibly dependent on the solution viscosity, which was adjusted by changing the glycerol content from 65% to 90% by volume in G/W mixtures. The kinetics of the GR of NO to Mp in viscous solutions was well represented by an exponential function with a time constant of ca. 11 ps. Although the kinetic traces of the GR of NO to Mp in solutions with three different viscosities (18, 81, and 252 cP) almost overlap, they show a slight difference early in the decay process. The kinetic traces were also described by the diffusion-controlled reaction theory with a Coulomb potential. Since the ligand is deligated in a neutral form, an ionic pair of NO(-) and Mp(+) may be produced before forming the Mp-NO bond by an electron transfer from Mp to NO as the deligated NO is sufficiently near to the Fe atom of Mp. The strong reactivity between NO and ferrous heme may arise from the Coulomb interaction between the reacting pair, which is consistent with the harpooning mechanism for NO binding to heme.

Probing the electronic and geometric structure of ferric and ferrous myoglobins in physiological solutions by Fe K-edge absorption spectroscopy

We present an iron K-edge X-ray absorption study of carboxymyoglobin (MbCO), nitrosylmyoglobin (MbNO), oxymyoglobin (MbO2), cyanomyoglobin (MbCN), aquomet myoglobin (metMb) and unligated myoglobin (deoxyMb) in physiological media. The analysis of the XANES region is performed using the full-multiple scattering formalism, implemented within the MXAN package. This reveals trends within the heme structure, absent from previous crystallographic and X-ray absorption analysis. In particular, the iron-nitrogen bond lengths in the porphyrin ring converge to a common value of about 2 Å, except for deoxyMb whose bigger value is due to the doming of the heme. The trends of the Fe-Nε (His93) bond length is found to be consistent with the effect of ligand binding to the iron, with the exception of MbNO, which is explained in terms of the repulsive trans effect. We derive a high resolution description of the relative geometry of the ligands with respect to the heme and quantify the magnitude of the heme doming in the deoxyMb form. Finally, time-dependent density functional theory is used to simulate the pre-edge spectra and is found to be in good agreement with the experiment. The XAS spectra typically exhibit one pre-edge feature which arises from transitions into the unoccupied dσ and dπ - πligand* orbitals. 1s → dπ transitions contribute weakly for MbO2, metMb and deoxyMb. However, despite this strong Fe d contribution these transitions are found to be dominated by the dipole (1s → 4p) moment due to the low symmetry of the heme environment.

Local structure and dynamics of hemeproteins by X-ray absorption near edge structure spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a synchrotron radiation technique sensitive to the local structure and dynamics around the metal site of a heme containing protein. Advances in detection techniques and theoretical/computational platforms in the last 15 years allowed the use of XANES as a quantitative probe of the key structural determinants driving functional changes, both in a concerted way with protein crystallography and EXAFS (extended X-ray absorption fine structure), or as a stand-alone method to apply in the crystal state as well as in solution. Moreover, the local dynamics of the heme site has been deeply investigated, on one hand, coupling XANES to classical photolysis experiments at cryogenic temperatures; on the other hand, the intrinsic property of the synchrotron radiation to induce radiolysis events, has been exploited to investigate specific cryotrapped intermediates, using X-rays both as a pump and a probe. Insights on the XANES method and some specific examples are presented to illustrate these topics.

Ligand migration and binding in nonsymbiotic hemoglobins of Arabidopsis thaliana

We have studied carbon monoxide (CO) migration and binding in the nonsymbiotic hemoglobins AHb1 and AHb2 of Arabidopsis thaliana using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures. Both proteins have similar amino acid sequences but display pronounced differences in ligand binding properties, at both physiological and cryogenic temperatures. Near neutral pH, the distal HisE7 side chain is close to the heme-bound ligand in the majority of AHb1-CO molecules, as indicated by a low CO stretching frequency at 1921 cm(-1). In this fraction, two CO docking sites can be populated, the primary site B and the secondary site C. When the pH is lowered, a high-frequency stretching band at approximately 1964 cm(-1) grows at the expense of the low-frequency band, indicating that HisE7 protonates and, concomitantly, moves away from the bound ligand. Geminate rebinding barriers are markedly different for the two conformations, and docking site C is not accessible in the low-pH conformation. Rebinding of NO ligands was observed only from site B of AHb1, regardless of conformation. In AHb2, the HisE7 side chain is removed from the bound ligand; rebinding barriers are low, and CO molecules can populate only primary docking site B. These results are interpreted in terms of differences in the active site structures and physiological functions.

Ligand dynamics in heme proteins observed by Fourier transform infrared-temperature derivative spectroscopy

Fourier transform infrared (FTIR) spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. Nitric oxide and carbon monoxide are attractive physiologically relevant ligands because their bond stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands display changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300K) using specific temperature protocols for sample photodissociation can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy (TDS) has proven to be a particularly powerful technique to study protein-ligand interactions. The FTIR-TDS technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. Here we describe infrared cryo-spectroscopy and present a variety of applications to the study of protein-ligand interactions in heme proteins. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Charge separation, stabilization, and protein relaxation in photosystem II core particles with closed reaction center

The fluorescence kinetics of cyanobacterial photosystem II (PSII) core particles with closed reaction centers (RCs) were studied with picosecond resolution. The data are modeled in terms of electron transfer (ET) and associated protein conformational relaxation processes, resolving four different radical pair (RP) states. The target analyses reveal the importance of protein relaxation steps in the ET chain for the functioning of PSII. We also tested previously published data on cyanobacterial PSII with open RCs using models that involved protein relaxation steps as suggested by our data on closed RCs. The rationale for this reanalysis is that at least one short-lived component could not be described in the previous simpler models. This new analysis supports the involvement of a protein relaxation step for open RCs as well. In this model the rate of ET from reduced pheophytin to the primary quinone Q(A) is determined to be 4.1 ns(-1). The rate of initial charge separation is slowed down substantially from approximately 170 ns(-1) in PSII with open RCs to 56 ns(-1) upon reduction of Q(A). However, the free-energy drop of the first RP is not changed substantially between the two RC redox states. The currently assumed mechanistic model, assuming the same early RP intermediates in both states of RC, is inconsistent with the presented energetics of the RPs. Additionally, a comparison between PSII with closed RCs in isolated cores and in intact cells reveals slightly different relaxation kinetics, with a approximately 3.7 ns component present only in isolated cores.

A physical picture of protein dynamics and conformational changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an alpha-helix. Molecules performing these types of motion are in the "flexible state". This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.

Characterization of ligand migration mechanisms inside hemoglobins from the analysis of geminate rebinding kinetics

The presence of internal hydrophobic cavities and packing defects has been demonstrated for several small globular proteins, including hemoglobins. The reduced thermodynamic stability appears to be compensated for by the capability of controlling ligand diffusion through the protein matrix to the active site, possibly by stocking more than one reactant molecule in selected sites. Photolysis of carbon monoxide complexes of hemoglobins encapsulated in silica gels leads to multiphasic geminate rebinding kinetics at room temperature, reflecting rebinding also from different temporary docking sites inside the protein matrix. A careful analysis of the ligand rebinding kinetics allows the determination of the microscopic rates for the underlying reactions, including those governing the migration to and from the docking sites. This chapter describes the experimental approach used to characterize the ligand rebinding kinetics for heme proteins in silica gels after nanosecond laser flash photolysis and the computational methods necessary to retrieve the kinetic parameters.

Identification of biomolecular conformations from incomplete torsion angle observations by hidden Markov models

We present a novel method for the identification of the most important conformations of a biomolecular system from molecular dynamics or Metropolis Monte Carlo time series by means of Hidden Markov Models (HMMs). We show that identification is possible based on the observation sequences of some essential torsion or backbone angles. In particular, the method still provides good results even if the conformations do have a strong overlap in these angles. To apply HMMs to angular data, we use von Mises output distributions. The performance of the resulting method is illustrated by numerical tests and by application to a hybrid Monte Carlo time series of trialanine and to MD simulation results of a DNA-oligomer.

Ligand dynamics in heme proteins observed by Fourier transform infrared spectroscopy at cryogenic temperatures

Fourier transform infrared spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. From the variety of ligands that bind to the heme iron, nitric oxide and carbon monoxide are particularly attractive, as their bond-stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, the ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand-stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands reveal changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites, and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300 K) using specific temperature protocols for sample photodissociation thus can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy has proven to be a particularly powerful technique to study protein-ligand interactions. This technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. This chapter describes infrared cryospectroscopy techniques and presents examples that demonstrate their applicability to nitric oxide binding to heme proteins.

Asymmetric band profile of the Soret band of deoxymyoglobin is caused by electronic and vibronic perturbations of the heme group rather than by a doming deformation

We measured the Soret band of deoxymyoglobin (deoxyMb), myoglobin cyanide (MbCN), and aquo-metmyoglobin (all from horse heart) with absorption and circular dichroism (CD) spectroscopies. A clear non-coincidence was observed between the absorption and CD profiles of deoxyMb and MbCN, with the CD profiles red- and blueshifted with respect to the absorption band position, respectively. On the contrary, the CD and absorption profiles of aquametMb were nearly identical. The observed noncoincidence indicates a splitting of the excited B state due to heme-protein interactions. CD and absorption profiles of deoxyMb and MbCN were self-consistently analyzed by employing a perturbation approach for weak vibronic coupling as well as the relative intensities and depolarization ratios of seven bands in the respective resonance Raman spectra measured with B-band excitation. The respective B(y) component was found to dominate the observed Cotton effect of both myoglobin derivatives. The different signs of the noncoincidences between CD and absorption bands observed for deoxyMb and MbCN are due to different signs of the respective matrix elements of A(1g) electronic interstate coupling, which reflects an imbalance of Gouterman's 50:50 states. The splitting of the B band reflects contributions from electronic and vibronic perturbations of B(1g) symmetry. The results of our analysis suggest that the broad and asymmetric absorption band of deoxyMb results from this band splitting rather than from its dependence on heme doming. Thus, we are able to explain recent findings that the temperature dependences of CO rebinding to myoglobin and the Soret band profile are uncorrelated[Ormos et al., Proc. Natl. Acad. Sci U.S.A. 95, 6762 (1998)].

Conformational dependence of hemoglobin reactivity under high viscosity conditions: the role of solvent slaved dynamics

The concept of protein dynamic states is introduced. This concept is based on (i) protein dynamics being organized hierarchically with respect to solvent slaving and (ii) which tier of dynamics is operative over the time window of a given measurement. The protein dynamic state concept is used to analyze the kinetic phases derived from the recombination of carbon monoxide to sol-gel-encapsulated human adult hemoglobin (HbA) and select recombinant mutants. The temperature-dependent measurements are made under very high viscosity conditions obtained by bathing the samples in an excess of glycerol. The results are consistent with a given tier of solvent slaved dynamics becoming operative at a time delay (with respect to the onset of the measurement) that is primarily solvent- and temperature-dependent. However, the functional consequences of the dynamics are protein- and conformation-specific. The kinetic traces from both equilibrium populations and trapped allosteric intermediates show a consistent progression that exposes the role of both conformation and hydration in the control of reactivity. Iron-zinc symmetric hybrid forms of HbA are used to show the dramatic difference between the kinetic patterns for T state alpha and beta subunits. The overall results support a model for allostery in HbA in which the ligand-binding-induced transition from the deoxy T state to the high -affinity R state proceeds through a progression of T state intermediates.

Temperature-dependent heme kinetics with nonexponential binding and barrier relaxation in the absence of protein conformational substates

We present temperature-dependent kinetic measurements of ultrafast diatomic ligand binding to the "bare" protoheme (L(1)-FePPIX-L(2), where L(1) = H(2)O or 2-methyl imidazole and L(2) = CO or NO). We found that the binding of CO is temperature-dependent and nonexponential over many decades in time, whereas the binding of NO is exponential and temperature-independent. The nonexponential nature of CO binding to protoheme, as well as its relaxation above the solvent glass transition, mimics the kinetics of CO binding to myoglobin (Mb) but on faster time scales. This demonstrates that the nonexponential kinetic response observed for Mb is not necessarily due to the presence of protein conformational substates but rather is an inherent property of the solvated heme. The nonexponential kinetic data were analyzed by using a linear coupling model with a distribution of enthalpic barriers that fluctuate on slower time scales than the heme-CO recombination time. Below the solvent glass transition (T(g) approximately 180 K), the average enthalpic rebinding barrier for H(2)O-PPIX-CO was found to be approximately 1 kJ/mol. Above T(g), the barrier relaxes and is approximately 6 kJ/mol at 290 K. Values for the first two moments of the heme doming coordinate distribution extracted from the kinetic data suggest significant anharmonicity above T(g). In contrast to Mb, the protoheme shows no indication of the presence of "distal" enthalpic barriers. Moreover, the wide range of Arrhenius prefactors (10(9) to 10(11) s(-1)) observed for CO binding to heme under differing conditions suggests that entropic barriers may be an important source of control in this class of biochemical reactions.

Ligand recombination and a hierarchy of solvent slaved dynamics: the origin of kinetic phases in hemeproteins

Ligand recombination studies play a central role both for characterizing different hemeproteins and their conformational states but also for probing fundamental biophysical processes. Consequently, there is great importance to providing a foundation from which one can understand the physical processes that give rise to and modulate the large range of kinetic patterns associated with ligand recombination in myoglobins and hemoglobins. In this work, an overview of cryogenic and solution phase recombination phenomena for COMb is first reviewed and then a new paradigm is presented for analyzing the temperature and viscosity dependent features of kinetic traces in terms of multiple phases that reflect which tier(s) of solvent slaved protein dynamics is (are) operative on the photoproduct population during the time course of the measurement. This approach allows for facile inclusion of both ligand diffusion among accessible cavities and conformational relaxation effects. The concepts are illustrated using kinetic traces and MEM populations derived from the CO recombination process for wild type and mutant myoglobins either in sol-gel matrices bathed in glycerol or in trehalose-derived glassy matrices.

Spectroscopic study of substrate binding to the carbonmonoxy form of dehaloperoxidase from Amphitrite ornata

Dehaloperoxidase (DHP) is a globular heme enzyme found in the marine worm Amphitrite ornata that can catalyze the dehalogenation of halophenols to the corresponding quinones by using hydrogen peroxide as a cosubstrate. Its three-dimensional fold is surprisingly similar to that of the oxygen storage protein myoglobin (Mb). A key structural feature common to both DHP and Mb is the existence of multiple conformations of the distal histidine. In DHP, the conformational flexibility may be involved in promotion of substrate and cosubstrate entry and exit. Here we have explored the dynamics of substrate binding in DHP using Fourier transform infrared spectroscopy and flash photolysis. A number of discrete conformations at the active site were identified from the appearance of multiple CO absorbance bands in the infrared region of the spectrum. Upon photolysis at cryogenic temperatures, the CO molecules are trapped at docking sites within the protein matrix, as inferred from the appearance of several photoproduct bands characteristic of each site. Substrate binding stabilizes the protein by approximately 20 kJ/mol. The low yield of substrate-bound DHP at ambient temperature points toward a steric inhibition of substrate binding by carbon monoxide.

The contribution of heme propionate groups to the conformational dynamics associated with CO photodissociation from horse heart myoglobin

Photoacoustic calorimetry and transient absorption spectroscopy were used to study conformational dynamics associated with CO photodissociation from horse heart myoglobin (Mb) reconstituted with either Fe protoporphyrin IX dimethylester (FePPDME), Fe octaethylporphyrin (FeOEP), or with native Fe protoporphyrin IX (FePPIX). The volume and enthalpy changes associated with the Fe-CO bond dissociation and formation of a transient deoxyMb intermediate for the reconstituted Mbs were found to be similar to those determined for native Mb (DeltaV1 = -2.5+/-0.6 ml mol(-1) and DeltaH1 = 8.1+/-3.0 kcal mol(-1)). The replacement of FePPIX by FeOEP significantly alters the conformational dynamics associated with CO release from protein. Ligand escape from FeOEP reconstituted Mb was determined to be roughly a factor of two faster (tau=330 ns) relative to native protein (tau=700 ns) and accompanying reaction volume and enthalpy changes were also found to be smaller (DeltaV2 = 5.4+/-2.5 ml mol(-1) and DeltaH2 = 0.7+/-2.2 kcal mol(-1)) than those for native Mb (DeltaV2 = 14.3+/-0.8 ml mol(-1) and DeltaH2 = 7.8+/-3.5 kcal mol(-1)). On the other hand, volume and enthalpy changes for CO release from FePPIX or FePPDME reconstituted Mb were nearly identical to those of the native protein. These results suggest that the hydrogen bonding network between heme propionate groups and nearby amino acid residues likely play an important role in regulating ligand diffusion through protein matrix. Disruption of this network leads to a partially open conformation of protein with less restricted ligand access to the heme binding pocket.

Bacteriorhodopsin photocycle at cryogenic temperatures reveals distributed barriers of conformational substates

The time course of thermal reactions after illumination of 100% humidified bacteriorhodopsin films was followed with FTIR spectroscopy between 125 and 195 K. We monitored the conversion of the initial photoproduct, K, to the next, L intermediate, and a shunt reaction of the L state directly back to the initial BR state. Both reactions can be described by either multiexponential kinetics, which would lead to apparent end-state mixtures that contain increasing amounts of the product, i.e., L or BR, with increasing temperature, or distributed kinetics. Conventional kinetic schemes that could account for the partial conversion require reversible reactions, branching, or parallel cycles. These possibilities were tested by producing K or L and monitoring their interconversion at a single temperature and by shifting the temperature upward or downward after an initial incubation and after their redistribution. The results are inconsistent with any conventional scheme. Instead, we attribute the partial conversions to the other alternative, distributed kinetics, observed previously in myoglobin, which arise from an ensemble of frozen conformational substates at the cryogenic temperatures. In this case, the time course of the reactions reflects the progressive depletion of distinct microscopic substates in the order of their increasing activation barriers, with a distribution width for K to L reaction of approximately 7 kJ/mol.

Myoglobin with chlorophyllous chromophores: influence on protein stability

The stabilities of myoglobin, apo-myoglobin, and of two myoglobins with chlorophyllous chromophores (Zn-pheophorbide a and Zn-bacteriopheophorbide a), have been studied by thermal and chemical denaturation. With guanidinium chloride, the stability order is myoglobin>Zn-pheophorbide-myoglobin>Zn-bacteriopheophorbide-myoglobin approximately apo-myoglobin. The thermal behavior is more complex. The transition temperature of thermal unfolding of the apoprotein (62.4 degrees C) is increased by Zn-pheophorbide a (83.9 degrees C) and Zn-bacteriopheophorbide a (82.6 degrees C) to a similar degree as by the native chromophore, heme (83.5 degrees C). The recovery with Zn-pheophorbide (92-98%) is even higher than with heme (74-76%), while with Zn-bacteriopheophorbide (40%) it is as low as with the apoprotein (42%). Recovery also depends on the rates of heating, and in particular the time spent at high temperatures. It is concluded that irreversibility of unfolding is related to loss of the chromophores, which are required for proper re-folding.

Ligand migration pathway and protein dynamics in myoglobin: a time-resolved crystallographic study on L29W MbCO

By using time-resolved x-ray crystallography at room temperature, structural relaxations and ligand migration were examined in myoglobin (Mb) mutant L29W from nanoseconds to seconds after photodissociation of carbon monoxide (CO) from the heme iron by nanosecond laser pulses. The data were analyzed in terms of transient kinetics by fitting trial functions to integrated difference electron density values obtained from select structural moieties, thus allowing a quantitative description of the processes involved. The observed relaxations are linked to other investigations on protein dynamics. At the earliest times, the heme has already completely relaxed into its domed deoxy structure, and there is no photo-dissociated CO visible at the primary docking site. Initial relaxations of larger globin moieties are completed within several hundred nanoseconds. They influence the concomitant migration of photo-dissociated CO to the Xe1 site, where it appears at approximately 300 ns and leaves again at approximately 1.5 ms. The extremely long residence time in Xe1 as compared with wild-type MbCO implies that, in the latter protein, the CO exits the protein from Xe1 predominantly via the distal pocket. A well-defined deligated state is populated between approximately 2 micros and approximately 1 ms; its structure is very similar to the equilibrium deoxy structure. Between 1.5 and 20 ms, no CO is visible in the protein interior; it is either distributed among many sites within the protein or has escaped to the solvent. Finally, recombination at the heme iron occurs after >20 ms.

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

Picosecond thermometer in the amide I band of myoglobin

The amide I and II bands in myoglobin show a heterogeneous temperature dependence, with bands at 6.17 and 6.43 microm which are more intense at low temperatures. The amide I band temperature dependence is on the long wavelength edge of the band, while the short wavelength side has almost no temperature dependence. We compare concepts of anharmonic solid-state crystal physics and chemical physics for the origins of these bands. We suggest that the long wavelength side is composed of those amino acids which hydrogen bond to the hydration shell of the protein, and that temperature dependent bands can be used to determine the time it takes vibrational energy to flow into the hydration shell. We determine that vibrational energy flow to the hydration shell from the amide I takes approximately 20 ps to occur.

The structural dynamics of myoglobin

Conformational fluctuations in proteins were initially invoked to explain the observation that diffusion of small ligands through the matrix is a global phenomenon. Small globular proteins contain internal cavities that play a role not only in matrix dynamics but also in controlling function, tracing a pathway for the diffusion of the ligand to and from the active site. This is the main point addressed in this Review, which presents pertinent information obtained on myoglobin (Mb). Mb, a simple globular heme protein which binds reversibly oxygen and other ligands. The bond between the heme Fe(II) and gaseous ligands can be photodissociated by a laser pulse, generating a non-equilibrium population of protein structures that relaxes on a picosecond to millisecond time range. This process is associated with migration of the ligand to internal cavities of the protein, which are known to bind xenon. Some of the results obtained by laser photolysis, molecular dynamics simulations, and X-ray diffraction of intermediate states of wild-type and mutant myoglobins are summarized. The extended relaxation of the globin moiety directly observed by Laue crystallography reflects re-equilibration among conformational substates known to play an essential role in controlling protein function.

Structural Dynamics Controls Nitric Oxide Affinity in Nitrophorin 4

Nitrophorin 4 (NP4) is one of seven nitric oxide (NO) transporting proteins in the blood-sucking insect Rhodnius prolixus. In its physiological function, NO binds to a ferric iron centered in a highly ruffled heme plane. Carbon monoxide (CO) also binds after reduction of the heme iron. Here we have used Fourier transform infrared spectroscopy at cryogenic temperatures to study CO and NO binding and migration in NP4, complemented by x-ray cryo-crystallography on xenon-containing NP4 crystals to identify cavities that may serve as ligand docking sites. Multiple infrared stretching bands of the heme-bound ligands indicate different active site conformations with varying degrees of hydrophobicity. Narrow infrared stretching bands are observed for photodissociated CO and NO; temperature-derivative spectroscopy shows that these bands are associated with ligand docking sites close to the extremely reactive heme iron. No rebinding from distinct secondary sites was detected, although two xenon binding cavities were observed in the x-ray structure. Photolysis studies at approximately 200 K show efficient NO photoproduct formation in the more hydrophilic, open NP4 conformation. This result suggests that ligand escape is facilitated in this conformation, and blockage of the active site by water hinders immediate reassociation of NO to the ferric iron. In the closed, low-pH conformation, ligand escape from the active site of NP4 is prevented by an extremely reactive heme iron and the absence of secondary ligand docking sites.

Structural Dynamics in the Active Site of Murine Neuroglobin and Its Effects on Ligand Binding

We have examined the effects of active site residues on ligand binding to the heme iron of mouse neuroglobin using steady-state and time-resolved visible spectroscopy. Absorption spectra of the native protein, mutants H64L and K67L and double mutant H64L/K67L were recorded for the ferric and ferrous states over a wide pH range (pH 4-11), which allowed us to identify a number of different species with different ligands at the sixth coordination, to characterize their spectroscopic properties, and to determine the pK values of active site residues. In flash photolysis experiments on CO-ligated samples, reaction intermediates and the competition of ligands for the sixth coordination were studied. These data provide insights into structural changes in the active site and the role of the key residues His64 and Lys67. His64 interferes with exogenous ligand access to the heme iron. Lys67 sequesters the distal pocket from the solvent. The heme iron is very reactive, as inferred from the fast ligand binding kinetics and the ability to bind water or hydroxyl ligands to the ferrous heme. Fast bond formation favors geminate rebinding; yet the large fraction of bimolecular rebinding observed in the kinetics implies that ligand escape from the distal pocket is highly efficient. Even slight pH variations cause pronounced changes in the association rate of exogenous ligands near physiological pH, which may be important in functional processes.

Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants

We have combined Fourier transform infrared/temperature derivative (FTIR-TDS) spectroscopy at cryogenic temperatures and flash photolysis at ambient temperature to examine the effects of polar and bulky amino acid replacements of the highly conserved distal valine 68 in sperm whale myoglobin. In FTIR-TDS experiments, the CO ligand can serve as an internal voltmeter that monitors the local electrostatic field not only at the active site but also at intermediate ligand docking sites. Mutations of residue 68 alter size, shape, and electric field of the distal pocket, especially in the vicinity of the primary docking site (state B). As a consequence, the infrared bands associated with the ligand at site B are shifted. The effect is most pronounced in mutants with large aromatic side chains. Polar side chains (threonine or serine) have only little effect on the peak frequencies. Ligands that migrate toward more remote sites C and D give rise to IR bands with altered frequencies. TDS experiments separate the photoproducts according to their recombination temperatures. The rates and extent of ligand migration among internal cavities at cryogenic temperatures can be used to interpret geminate and bimolecular O2 and CO recombination at room temperature. The kinetics of geminate recombination can be explained by steric arguments alone, whereas both the polarity and size of the position 68 side chain play major roles in regulating bimolecular ligand binding from the solvent.

Spectroscopically and kinetically distinct conformational populations of sol-gel-encapsulated carbonmonoxy myoglobin. A comparison with hemoglobin

We have used sol-gel encapsulation protocols to trap kinetically and spectroscopically distinct conformational populations of native horse carbonmonoxy myoglobin. The method allows for direct comparison of functional and spectroscopic properties of equilibrium and non-equilibrium populations under the same temperature and viscosity conditions. The results implicate tertiary structure changes that include the proximal heme environment in the mechanism for population-specific differences in the observed rebinding kinetics. Differences in the resonance Raman frequency of nu(Fe-His), the iron-proximal histidine stretching mode, are attributed to differences in the positioning of the F helix. For myoglobin, the degree of separation between the F helix and the heme is assigned as the conformational coordinate that modulates both this frequency and the innermost barrier controlling CO rebinding. A comparison with the behavior of encapsulated derivatives of human adult hemoglobin indicates that these CO binding-induced conformational changes are qualitatively similar to the tertiary changes that occur within both the R and T quaternary states. Protein-specific differences in the time scale for the proposed F helix relaxation are attributed to variations in the intra-helical hydrogen bonding patterns that help stabilize the position of the F helix.

Heme charge-transfer band III is vibronically coupled to the Soret band

A complete resonance Raman excitation profile of the heme charge-transfer band known as band III is presented. The data obtained throughout the near-infrared region show preresonance with the Q-band, but the data also clearly show the enhancement of a number of modes in the spectral region of band III. Only nontotally symmetric modes are observed to have resonance enhancement in the band III region. The observed resonance enhancements in modes of B(1g) symmetry are compared with the enhancements of those same modes in the excitation profiles of the Q-band of deoxy myoglobin, also presented here for this first time. The Q-band data agree well with the theory of vibronic coupling in metalloporphyrins (Shelnutt, J. A. J. Chem. Phys. 1981, 74, 6644-6657). The strong vibronic coupling of the Q-band of the deoxy form of hemes is discussed in terms of the enhancement of modes with both B(1g) and A(2g) symmetry. The comparison between the Q-band and band III reveals that, consistent with the theory, only modes of B(1g) symmetry are enhanced in the vicinity of band III. These results show that band III is vibronically coupled to the Soret band. The coupling of band III to modes with strong rhombic distortion of the heme macrocycle calls into question the hypothesis that the axial iron out-of-plane displacement is primarily responsible for the structure-dynamics correlations observed in myoglobin.

The effect of ligand dynamics on heme electronic transition band III in myoglobin

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.

Influence of static and dynamic disorder on the visible and infrared absorption spectra of carbonmonoxy horseradish peroxidase

Spectroscopy of horseradish peroxidase with and without the substrate analog, benzohydroxamic acid, was monitored in a glycerol/water solvent as a function of temperature. It was determined from the water infrared (IR) absorption that the solvent has a glass transition at 170-180 K. In the absence of substrate, both the heme optical Q(0,0) absorption band and the IR absorption band of CO bound to heme broaden markedly upon heating from 10-300 K. The Q(0,0) band broadens smoothly in the whole temperature interval, whereas the IR bandwidth is constant in the glassy matrix and increases from 7 to 16 cm(-1) upon heating above the glass transition. Binding of substrate strongly diminishes temperature broadening of both the bands. The results are consistent with the view that the substrate strongly reduces the amplitude of motions of amino acids forming the heme pocket. The main contribution to the Q(0,0) bandwidth arises from the heme vibrations that are not affected by the phase transition. The CO band thermal broadening stems from the anharmonic coupling with motions of the heme environment, which, in the glassy state, are frozen in. Unusually strong temperature broadening of the CO band is interpreted to be caused by thermal population of a very flexible excited conformational substrate. Analysis of literature data on the thermal broadening of the A(0) band of Mb(CO) (Ansari et al., 1987. Biophys. Chem. 26:337-355) shows that such a state presents itself also in myoglobin.

Slow relaxation process in DNA

A dynamic transition at temperatures ∼200-230K is observed in manyhydrated bio-polymers. It shows up as a sharp increase of the mean-squaredatomic displacements above this temperature range. We present neutronscattering data of DNA at different levels of hydration. The analysis showsthat the dynamic transition in DNA is related to a slow relaxation processin the MHz-GHz frequency range. This slow relaxation process iscompletely suppressed in the dry DNA sample where no dynamic transitionwas observed. The nature of the slow process is discussed. We ascribe it toa global relaxation of DNA molecule that involves cooperative motion ofmany base-pairs and backbone.

Diffractive optics-based heterodyne-detected four-wave mixing signals of protein motion: from "protein quakes" to ligand escape for myoglobin

Ligand transport through myoglobin (Mb) has been observed by using optically heterodyne-detected transient grating spectroscopy. Experimental implementation using diffractive optics has provided unprecedented sensitivity for the study of protein motions by enabling the passive phase locking of the four beams that constitute the experiment, and an unambiguous separation of the Real and Imaginary parts of the signal. Ligand photodissociation of carboxymyoglobin (MbCO) induces a sequence of events involving the relaxation of the protein structure to accommodate ligand escape. These motions show up in the Real part of the signal. The ligand (CO) transport process involves an initial, small amplitude, change in volume, reflecting the transit time of the ligand through the protein, followed by a significantly larger volume change with ligand escape to the surrounding water. The latter process is well described by a single exponential process of 725 +/- 15 ns. at room temperature. The overall dynamics provide a distinctive signature that can be understood in the context of segmental protein fluctuations that aid ligand escape via a few specific cavities, and they suggest the existence of discrete escape pathways.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Slow relaxation process in DNA at different levels of hydration

There are many speculations about the dynamic transition observed in hydrated bio-polymers at temperatures T ∼ 200 - 230 K being an important factor for enabling of their functions. The transition shows up as a sharp increase of atomic mean-squared displacements above this temperature. The nature of the dynamic transition is not yet clear. Using inelastic neutron scattering we show in this Note that the transition in DNA is related to the appearance of a slow relaxation process. Decrease in the hydration level suppresses the process and the dynamic transition. It is found that, in terms of dynamics, the decrease in water content is similar in effect to a decrease in temperature. The obtained results support the idea that the dynamic transition is mediated by the water of hydration since bulk water has a dynamic transition around the same temperature.

Protein dynamics in an intermediate state of myoglobin: optical absorption, resonance Raman spectroscopy, and x-ray structure analysis

A metastable state of myoglobin is produced by reduction of metmyoglobin at low temperatures. This is done either by irradiation with x-rays at 80 K or by electron transfer from photoexcited tris(2, 2'-bipyridine)-ruthenium(II) at 20 K. At temperatures above 150 K, the conformational transition toward the equilibrium deoxymyoglobin is observed. X-ray crystallography, Raman spectroscopy, and temperature-dependent optical absorption spectroscopy show that the metastable state has a six-ligated iron low-spin center. The x-ray structure at 115K proves the similarity of the metastable state with metmyoglobin. The Raman spectra yield the high-frequency vibronic modes and give additional information about the distortion of the heme. Analysis of the temperature dependence of the line shape of the Soret band reveals that a relaxation within the metastable state starts at approximately 120 K. Parameters representative of static properties of the intermediate state are close to those of CO-ligated myoglobin, while parameters representative of dynamics are close to deoxymyoglobin. Thus within the metastable state the relaxation to the equilibrium is initiated by changes in the dynamic properties of the active site.

Ligand binding and conformational motions in myoglobin

Myoglobin, a small globular haem protein that binds gaseous ligands such as O2, CO and NO reversibly at the haem iron, serves as a model for studying structural and dynamic aspects of protein reactions. Time-resolved spectroscopic measurements after photodissociation of the ligand revealed a complex ligand-binding reaction with multiple kinetic intermediates, resulting from protein relaxation and movements of the ligand within the protein. To observe the structural changes induced by ligand dissociation, we have carried out X-ray crystallographic investigations of carbon monoxy-myoglobin (MbCO mutant L29W) crystals illuminated below and above 180 K, complemented by time-resolved infrared spectroscopy of CO rebinding. Here we show that below 180 K photodissociated ligands migrate to specific sites within an internal cavity--the distal haem pocket--of an essentially immobilized, frozen protein, from where they subsequently rebind by thermally activated barrier crossing. Upon photodissociation above 180 K, ligands escape from the distal pocket, aided by protein fluctuations that transiently open exit channels. We recover most of the ligands in a cavity on the opposite side of the haem group.

Iron-histidine resonance Raman band of deoxyheme proteins: effects of anharmonic coupling and glass-liquid phase transition

Weak anharmonic coupling of two soft molecular vibrations is shown to cause pronounced temperature dependence of the corresponding resonance Raman bands. The developed theory is used to interpret the temperature dependence of the iron-histidine band of deoxyheme proteins and model compounds. It is shown that anharmonic coupling of the iron-histidine and heme doming vibrations must cause pronounced broadening of the band, its asymmetry, and shift of its maximum to the red upon heating. It also can lead to a structured shape of this band at room temperature. Proper consideration of the anharmonic coupling allows simulation of the temperature dependence of the iron-histidine band shape of horse heart myoglobin in the temperature interval of 10-300 K, using the minimum number of necessary parameters. Analysis of this temperature dependence clearly shows that the iron-histidine band of deoxyheme proteins is sensitive to the glass-liquid phase transition in the protein hydration shell, which takes place at 160-190 K.

Dehydration and crystallization of trehalose and sucrose glasses containing carbonmonoxy-myoglobin

We report a study wherein we contemporarily measured 1) the dehydration process of trehalose or sucrose glasses embedding carbonmonoxy-myoglobin (MbCO) and 2) the evolution of the A substates in saccharide-coated MbCO. Our results indicate that microcrystallization processes, sizeably different in the two saccharides, take place during dehydration; moreover, the microcrystalline structure is maintained unless the dry samples are equilibrated with a humidity >/=75% (>/=60%) at 25 degrees C for the trehalose (sucrose) sample. The evolution of the parameters that characterize the A substates of MbCO indicates that 1) the effects of water withdrawal are analogous in samples dried in the presence or in the absence of sugars, although much larger effects are observed in the samples without sugar; 2) the distribution of A substates is determined by the overall matrix structure and not only by the sample water content; and 3) the population of A0 substate (i. e., the substate currently put in relation with MbCO molecules having the distal histidine out of the heme pocket) is largely enhanced during the dehydration process. However, after rehumidification its population is largely decreased with respect to the values obtained, at similar water content, during the first dehydration run.

A synthetic picture of intramolecular dynamics of proteins. Towards a contemporary statistical theory of biochemical processes

An increasing body of experimental evidence indicates the slow character of internal dynamics of native proteins. The important consequence of this is that theories of chemical reactions, used hitherto, appear inadequate for description of most biochemical reactions. Construction of a contemporary, truly advanced statistical theory of biochemical processes will need simple but realistic models of microscopic dynamics of biomolecules. In this review, intended to be a contribution towards this direction, three topics are considered. First, an intentionally simplified picture of dynamics of native proteins which emerges from recent investigations is presented. Fast vibrational modes of motion, of periods varying from 10(-14) to 10(-11) s, are contrasted with purely stochastic conformational transitions. Significant evidence is adduced that the relaxation time spectrum of the latter spreads in the whole range from 10(-11) to 10(5) s or longer, and up to 10(-7) s it is practically quasi-continuous. Next, the essential ideas of the theory of reaction rates based on stochastic models of intramolecular dynamics are outlined. Special attention is paid to reactions involving molecules in the initial conformational substrates confirmed to the transition state, which is realized in actual experimental situations. And finally, the two best experimentally justified classes of models of conformational transition dynamics, symbolically referred to as "protein glass" and "protein machine", are described and applied to the interpretation of a few simple biochemical processes, perhaps the most important result reported is the demonstration of the possibility of predominance of the short initial condition-dependent stage of protein involved reactions over the main stage described by the standard kinetics. This initial stage, and not the latter, is expected to be responsible for the coupling of component reactions in the complete enzymatic cycles as well as more complex processes of biological free energy transduction.

Structural factors controlling ligand binding to myoglobin: a kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as "kinetic hole burning." The study of these correlations gives direct insights into the structure-function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Electron transfer and protein dynamics in the photosynthetic reaction center

We have measured the kinetics of electron transfer (ET) from the primary quinone (Q(A)) to the special pair (P) of the reaction center (RC) complex from Rhodobacter sphaeroides as a function of temperature (5-300 K), illumination protocol (cooled in the dark and under illumination from 110, 160, 180, and 280 K), and warming rate (1.3 and 13 mK/s). The nonexponential kinetics are interpreted with a quantum-mechanical ET model (Fermi's golden rule and the spin-boson model), in which heterogeneity of the protein ensemble, relaxations, and fluctuations are cast into a single coordinate that relaxes monotonically and is sensitive to all types of relaxations caused by ET. Our analysis shows that the structural changes that occur in response to ET decrease the free energy gap between donor and acceptor states by 120 meV and decrease the electronic coupling between donor and acceptor states from 2.7 x 10(-4) cm(-1) to 1.8 x 10(-4) cm(-1). At cryogenic temperatures, conformational changes can be slowed or completely arrested, allowing us to monitor relaxations on the annealing time scale (approximately 10(3)-10(4) s) as well as the time scale of ET (approximately 100 ms). The relaxations occur within four broad tiers of conformational substates with average apparent Arrhenius activation enthalpies of 17, 50, 78, and 110 kJ/mol and preexponential factors of 10(13), 10(15), 10(21), and 10(25) s(-1), respectively. The parameterization provides a prediction of the time course of relaxations at all temperatures. At 300 K, relaxations are expected to occur from 1 ps to 1 ms, whereas at lower temperatures, even broader distributions of relaxation times are expected. The weak dependence of the ET rate on both temperature and protein conformation, together with the possibility of modeling heterogeneity and dynamics with a single conformational coordinate, make RC a useful model system for probing the dynamics of conformational changes in proteins.