Linkage of functional and structural heterogeneity in proteins: dynamic hole burning in carboxymyoglobin

Supramolecular Copolymer Constructed by Hierarchical Self-Assembly of Orthogonal Host-Guest, H-Bonding, and Coordination Interactions

Supramolecular copolymers with complex architectures and emergent functions constitute a class of challenging but enticing synthetic targets in polymer science. Individual building blocks can be tailored to endow a resulting supramolecular copolymer with increased structural and functional complexity. Herein, we describe the construction of a linear supramolecular copolymer comprising mechanically interlocked segments with hydrogen-bonding metallorhomboidal units. Specifically, a hierarchical supramolecular polymerization of a crown ether-based [2]rotaxane and a discrete organoplatinum(II) metallacycle driven by 2-ureido-4-pyrimidinone (UPy) quadruple hydrogen bonding provides the impetus for its formation. This system demonstrates enhanced structural complexity accessed by the unification of orthogonal noncovalent interactions: metal coordination, host-guest chemistry, and multiple hydrogen bonding interfaces.

Biopolymers and supramolecular polymers as biomaterials for biomedical applications

Protein- and peptide-based structural biopolymers are abundant building blocks of biological systems. Either in their natural forms, such as collagen, silk or fibronectin, or as related synthetic materials they can be used in various technologies. An emerging area is that of biomimetic materials inspired by protein-based biopolymers, which are made up of small molecules rather than macromolecules and can therefore be described as supramolecular polymers. These materials are very useful in biomedical applications because of their ability to imitate the extracellular matrix both in architecture and their capacity to signal cells. This article describes important features of the natural extracellular matrix and highlight how these features are being incorporated into biomaterials composed of biopolymers and supramolecular polymers. We particularly focus on the structures, properties, and functions of collagen, fibronectin, silk, and the supramolecular polymers inspired by them as biomaterials for regenerative medicine.

Supramolecular polymers with tunable topologies via hierarchical coordination-driven self-assembly and hydrogen bonding interfaces

A powerful strategy to obtain complex supramolecular materials is the bottom-up construction of noncovalently bound materials by hierarchical self-assembly. This assembly process involves stepwise, uniform increases to the architectural complexity of a substrate, starting from discrete precursors and growing in dimensionality through controlled reactivity to a final product. Herein, two orthogonal processes are exploited: coordination-driven self-assembly and hydrogen bonding. The former relies on the predictable formation of metal-ligand bonds wherein the directionalities of the rigid precursors used determines the structural outcome. The latter uses 2-ureido-4-pyrimidinone interfaces that are structurally robust by virtue of the quadruple hydrogen bonding that can occur between subunits. By combining these two processes into a single system, it is possible to generate hierarchical materials that preserve the attractive tunability associated with discrete supramolecular coordination complexes. For instance, the synthesis of a one-dimensional chain comprising linked metalla-rhomboids is readily adapted to a 2D cross-linked hexagonal network by simply selecting a different metal acceptor precursor as an assembly component. The specific interactions between subunits, in this case platinum(II)-pyridyl bonds and the quadruple H-bonding of ureidopyrimidinone, are unchanged, establishing a unique strategy to obtain supramolecular polymers with marked topological differences with minimal synthetic redesign. In addition, the structural rigidity imposed by the inclusion of the platinum metallacycles serves to minimize the formation of cyclic oligomers, increasing the efficacy of formation and improving the properties of the resultant materials. Furthermore, this study taps the potential of organoplatinum(II) metallacycles in materials science.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Ribonuclease S dynamics measured using a nitrile label with 2D IR vibrational echo spectroscopy

A nitrile-labeled amino acid, p-cyanophenylalanine, is introduced near the active site of the semisynthetic enzyme ribonuclease S to serve as a probe of protein dynamics and fluctuations. Ribonuclease S is the limited proteolysis product of subtilisin acting on ribonuclease A, and consists of a small fragment including amino acids 1-20, the S-peptide, and a larger fragment including residues 21-124, the S-protein. A series of two-dimensional vibrational echo experiments performed on the nitrile-labeled S-peptide and the RNase S are described. The time-dependent changes in the two-dimensional infrared vibrational echo line shapes are analyzed using the center line slope method to obtain the frequency-frequency correlation function (FFCF). The observations show that the nitrile probe in the S-peptide has dynamics that are similar to, but faster than, those of the single amino acid p-cyanophenylalanine in water. In contrast, the dynamics of the nitrile label when the peptide is bound to form ribonuclease S are dominated by homogeneous dephasing (motionally narrowed) contributions with only a small contribution from very fast inhomogeneous structural dynamics. The results provide insights into the nature of the structural dynamics of the ribonuclease S complex. The equilibrium dynamics of the nitrile labeled S-peptide and the ribonuclease S complex are also investigated by molecular dynamics simulations. The experimentally determined FFCFs are compared to the FFCFs obtained from the molecular dynamics simulations, thereby testing the capacity of simulations to determine the amplitudes and time scales of protein structural fluctuations on fast time scales under thermal equilibrium conditions.

Differences between protein dynamics of hemoglobin upon dissociation of oxygen and carbon monoxide

Protein dynamics of human adult hemoglobin (HbA) upon ligand photolysis of oxygen (O(2)) and carbon monoxide (CO) was investigated using time-resolved resonance Raman (TR(3)) spectroscopy. The TR(3) spectra of the both photoproducts at 1-ns delay differed from that of the equilibrium deligated form (deoxy form) in the frequencies of the iron-histidine stretching [ν(Fe-His)] and methine wagging (γ(7)) modes, and the band intensity of pyrrole stretching and substituent bending (ν(8)) modes. Spectral changes of the O(2) photoproduct in the submicrosecond region were faster than those of the CO photoproduct, indicating that the structural dynamics following the photodissociation is ligand dependent for HbA. In contrast, no ligand dependence of the dynamics was observed for myoglobin, which has a structure similar to that of the subunit of HbA. The structural dynamics and relevance to the functionality of HbA also are discussed.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Imaging protein statistical substate occupancy in a spectrum-function phase space

Hemeprotein ligand rebinding studies reveal varying IR absorbance and rebinding functions across a cryogenic ensemble. Since IR-active vibrations and rebinding barriers couple to structural coordinates, spectral and functional heterogeneity arise from conformational heterogeneity. Modeling rebinding data as a spectrally resolved superposition of first-order rate processes and employing maximum entropy regularization, protein heterogeneity is imaged as an ensemble occupancy of a spectrum-function phase space. Results from myoglobin rebinding carbon monoxide are discussed.

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.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

Charge density-dependent modifications of hydration shell waters by Hofmeister ions

Gadolinium (Gd(3+)) vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe, as a function of added Hofmeister series salts, changes in the OH stretching frequency derived from first-shell waters of aqueous Gd(3+) and of Gd(3+) coordinated to three different types of molecules: (i) a chelate (EDTA), (ii) structured peptides (mSE3/SE2) of the lanthanide-binding tags (LBTs) family with a single high-affinity binding site, and (iii) a calcium-binding protein (calmodulin) with four binding sites. The vibronic sideband (VSB) corresponding to the OH stretching mode of waters coordinated to Gd(3+), whose frequency is inversely correlated with the strength of the hydrogen bonding to neighboring waters, exhibits an increase in frequency when Gd(3+) becomes coordinated to either EDTA, calmodulin, or mSE3 peptide. In all of these cases, the addition of cation chloride or acetate salts to the solution increases the frequency of the vibronic band originating from the OH stretching mode of the coordinated waters in a cation- and concentration-dependent fashion. The cation dependence of the frequency increase scales with charge density of the cations, giving rise to an ordering consistent with the Hofmeister ordering. On the other hand, water Raman spectroscopy shows no significant change upon addition of these salts. Additionally, it is shown that the cation effect is modulated by the specific anion used. The results indicate a mechanism of action for Hofmeister series ions in which hydrogen bonding among hydration shell waters is modulated by several factors. High charge density cations sequester waters in a configuration that precludes strong hydrogen bonding to neighboring waters. Under such conditions, anion effects emerge as anions compete for hydrogen-bonding sites with the remaining free waters on the surface of the hydration shell. The magnitude of the anion effect is both cation and Gd(3+)-binding site specific.

Dynamic fluctuation of proteins watched in real time

The dynamic nature of protein function is a fundamental concept in the physics of proteins. Although the basic general ideas are well accepted most experimental evidence has an indirect nature. The detailed characterization of the dynamics is necessary for the understanding in detail. The dynamic fluctuations thought crucial for the function span an extremely broad time, starting from the picosecond regime. Recently, a few new experimental techniques emerged that permit the observation of dynamical phenomena directly. Notably, pulsed infrared (IR) spectroscopy has been applied with great success to observe structural changes with picosecond time resolution. Using two-dimensional-IR vibrational echo chemical exchange spectroscopy Ishikawa and co-workers [Ishikawa et al. (2008), Proc. Natl. Acad. Sci. U.S.A. 101, 14402-14407] managed to observe the transition between well defined conformational substrates of carbonmonoxy myoglobin directly. This is an important step in improving our insight into the details of protein function.

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.

Influence of distal residue B10 on CO dynamics in myoglobin and neuroglobin

For many years, myoglobin has served as a paradigm for structure-function studies in proteins. Ligand binding and migration within myoglobin has been studied in great detail by crystallography and spectroscopy, showing that gaseous ligands such as O(2), CO, and NO not only bind to the heme iron but may also reside transiently in three internal ligand docking sites, the primary docking site B and secondary sites C and D. These sites affect ligand association and dissociation in specific ways. Neuroglobin is another vertebrate heme protein that also binds small ligands. Ligand migration pathways in neuroglobin have not yet been elucidated. Here, we have used Fourier transform infrared temperature derivative spectroscopy at cryogenic temperatures to compare the influence of the side chain volume of amino acid residue B10 on ligand migration to and rebinding from docking sites in myoglobin and neuroglobin.

Ligand migration in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

AHb1 is a hexacoordinated type 1 nonsymbiotic hemoglobin recently discovered in Arabidopsis thaliana. To gain insight into the ligand migration inside the protein, we studied the CO rebinding kinetics of AHb1 encapsulated in silica gels, in the presence of glycerol. The CO rebinding kinetics after nanosecond laser flash photolysis exhibits complex ligand migration patterns, consistent with the existence of discrete docking sites in which ligands can temporarily be stored before rebinding to the heme at different times. This finding may be of relevance to the physiological NO dioxygenase activity of this protein, which requires sequential binding of two substrates, NO and O2, to the heme.

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.

Conformational substates of horse heart cytochrome c exhibit different thermal unfolding of the heme cavity

The charge transfer (CT) band at 695 nm in the spectrum of ferri-cytochrome c is highly asymmetric, indicating conformational heterogeneity due to the coexistence of different conformational substates. We have measured the respective band profile of horse heart ferri-cytochrome c as a function of temperature between 283 K (10 degrees C) and 333 K (60 degrees C) and found that the well-known decrease of the absorptivity is wavenumber-dependent and exhibits a biphasic behavior. This indicates that the underlying conformational substates differ in their thermodynamic stability with respect to the structural changes associated with the disappearance of the 695 nm band, which eventually (at high temperatures) involves the replacement of M80 by a nearby lysine residue. Our data further indicates that the thermal unfolding process involves two structurally different intermediate states.

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.

Optical band splitting and electronic perturbations of the heme chromophore in cytochrome C at room temperature probed by visible electronic circular dichroism spectroscopy

We have measured the electronic circular dichroism (ECD) of the ferri- and ferro-states of several natural cytochrome c derivatives (horse heart, chicken, bovine, and yeast) and the Y67F mutant of yeast in the region between 300 and 750 nm. Thus, we recorded the ECD of the B- and Q-band region as well as the charge-transfer band at approximately 695 nm. The B-band region of the ferri-state displays a nearly symmetric couplet at the B0-position that overlaps with a couplet 790 cm-1 higher in energy, which we assigned to a vibronic side-band transition. For the ferro-state, the couplet is greatly reduced, but still detectable. The B-band region is dominated by a positive Cotton effect at energies lower than B0 that is attributed to a magnetically allowed iron-->heme charge-transfer transition as earlier observed for nitrosyl myoglobin and hemoglobin. The Q-band region of the ferri-state is poorly resolved, but displays a pronounced positive signal at higher wavenumbers. This must result from a magnetically allowed transition, possibly from the methionine ligand to the dxy-hole of Fe3+. For the ferro-state, the spectra resolve the vibronic structure of the Qv-band. A more detailed spectral analysis reveals that the positively biased spectrum can be understood as a superposition of asymmetric couplets of split Q0 and Qv-states. Substantial qualitative and quantitative differences between the respective B-state and Q-state ECD spectra of yeast and horse heart cytochrome c can clearly be attributed to the reduced band splitting in the former, which results from a less heterogeneous internal electric field. Finally, we investigated the charge-transfer band at 695 nm in the ferri-state spectrum and found that it is composed of at least three bands, which are assignable to different taxonomic substates. The respective subbands differ somewhat with respect to their Kuhn anisotropy ratio and their intensity ratios are different for horse and yeast cytochrome c. Our data therefore suggests different substate populations for these proteins, which is most likely assignable to a structural heterogeneity of the distal Fe-M80 coordination of the heme chromophore.

Environment-Controlled Interchromophore Charge Transfer Transitions in Dipeptides Probed by UV Absorption and Electronic Circular Dichroism Spectroscopy

Charge transfer (CT) transitions between the C-terminal carboxylate and peptide group have been investigated for alanyl-X and X-alanine dipeptides by far-UV absorption and electronic circular dichroism (ECD) spectroscopy (where X represents different amino acid residues). The spectra used in the present study were obtained by subtracting the spectrum of the cationic species from that of the corresponding zwitterionic peptide spectrum. These spectra displayed three bands, e.g., band I between 44 and 50 kK (kK = 10(3) cm(-1)), band II at 53 kK, and band III above 55 kK, which were, respectively, assigned to a n(COO-) --> pi* CT transition, a pi(COO-) --> pi* CT transition, and a carboxylate pi --> pi* (NV1) transition, respectively By comparison of the intensity, bandwidth, and wavenumber position of band I of some of the investigated dipeptides, we found that positive charges on the N-terminal side chain (for X = K), and to a minor extent also the N-terminal proton, reduce its intensity. This can be understood in terms of attractive Coulomb interactions that stabilize the ground state over the charge transfer state. For alanylphenylalanine, we assigned band I to a n(COO-) --> pi* CT transition into the aromatic side chain, indicating that aromatic side chains interact electronically with the backbone. We also performed ECD measurements at different pH values (pH 1-6) for a selected subset of XA and AX peptides. By subtraction of the pH 1 spectrum from that observed at pH 6, the ECD spectrum of the CT transition was obtained. A titration curve of their spectra reveals a substantial dependence on the protonation state of the aspartic acid side chain of AD, which is absent in DA and AE. This most likely reflects a conformational transition of the C-terminus into a less extended state, though the involvement of a side chain --> peptide CT transition cannot be completely ruled out.

The Position 68(E11) Side Chain in Myoglobin Regulates Ligand Capture, Bond Formation with Heme Iron, and Internal Movement into the Xenon Cavities

After photodissociation, ligand rebinding to myoglobin exhibits complex kinetic patterns associated with multiple first-order geminate recombination processes occurring within the protein and a simpler bimolecular phase representing second-order ligand rebinding from the solvent. A smooth transition from cryogenic-like to solution phase properties can be obtained by using a combination of sol-gel encapsulation, addition of glycerol as a bathing medium, and temperature tuning (-15 --> 65 degrees C). This approach was applied to a series of double mutants, myoglobin CO (H64L/V68X, where X = Ala, Val, Leu, Asn, and Phe), which were designed to examine the contributions of the position 68(E11) side chain to the appearance and disappearance of internal rebinding phases in the absence of steric and polar interactions with the distal histidine. Based on the effects of viscosity, temperature, and the stereochemistry of the E11 side chain, the three major phases, B --> A, C --> A, and D --> A, can be assigned, respectively, to ligand rebinding from the following: (i) the distal heme pocket, (ii) the xenon cavities prior to large amplitude side chain conformational relaxation, and (iii) the xenon cavities after significant conformational relaxation of the position 68(E11) side chain. The relative amplitudes of the B --> A and C --> A phases depend markedly on the size and shape of the E11 side chain, which regulates sterically both ligand return to the heme iron atom and ligand migration to the xenon cavities. The internal xenon cavities provide a transient docking site that allows side chain relaxations and the entry of water into the vacated distal pocket, which in turn slows ligand recombination markedly.

CO rebinding to protoheme: investigations of the proximal and distal contributions to the geminate rebinding barrier

The rebinding kinetics of CO to protoheme (FePPIX) in the presence and absence of a proximal imidazole ligand reveals the magnitude of the rebinding barrier associated with proximal histidine ligation. The ligation states of the heme under different solvent conditions are also investigated using both equilibrium and transient spectroscopy. In the absence of imidazole, a weak ligand (probably water) is bound on the proximal side of the FePPIX-CO adduct. When the heme is encapsulated in micelles of cetyltrimethylammonium bromide (CTAB), photolysis of FePPIX-CO induces a complicated set of proximal ligation changes. In contrast, the use of glycerol-water solutions leads to a simple two-state geminate kinetic response with rapid (10-100 ps) CO recombination and a geminate amplitude that can be controlled by adjusting the solvent viscosity. By comparing the rate of CO rebinding to protoheme in glycerol solution with and without a bound proximal imidazole ligand, we find the enthalpic contribution to the proximal rebinding barrier, H(p), to be 11 +/- 2 kJ/mol. Further comparison of the CO rebinding rate of the imidazole bound protoheme with the analogous rate in myoglobin (Mb) leads to a determination of the difference in their distal free energy barriers: DeltaG(D) approximately 12 +/- 1 kJ/mol. Estimates of the entropic contributions, due to the ligand accessible volumes in the distal pocket and the xenon-4 cavity of myoglobin ( approximately 3 kJ/mol), then lead to a distal pocket enthalpic barrier of H(D) approximately 9 +/- 2 kJ/mol. These results agree well with the predictions of a simple model and with previous independent room-temperature measurements of the enthalpic MbCO rebinding barrier (18 +/- 2 kJ/mol).

Dynamics of hemoglobin in human erythrocytes and in solution: influence of viscosity studied by ultrafast vibrational echo experiments

Ultrafast spectrally resolved stimulated vibrational echo experiments are used to measure the vibrational dephasing of the CO stretching mode of hemoglobin-CO (HbCO) inside living human erythrocytes (red blood cells), in liquid solutions, and in a glassy matrix. A method is presented to overcome the adverse impact on the vibrational echo signal from the strong light scattering caused by the cells. The results from the cytoplasmic HbCO are compared to experiments on aqueous HbCO samples prepared in different buffers, solutions containing low and high concentrations of glycerol, and in a solid trehalose matrix. Measurements are also presented that provide an accurate determination of the viscosity at the very high Hb concentration that is found inside the cells. It is demonstrated that the dynamics of the protein, as sensed by the CO ligand, are the same inside the erythrocytes and in aqueous solution and are independent of the viscosity. In solutions that are predominantly glycerol, the dynamics are modified somewhat but are still independent of viscosity. The experiments in trehalose give the dynamics at infinite viscosity and are used to separate the viscosity-dependent dynamics from the viscosity-independent dynamics. Although the HbCO dynamics are the same in the red blood cell and in the equivalent aqueous solutions, differences in the absorption spectra show that the distribution of a protein's equilibrium substates is sensitive to small pH differences.

Coupling of protein relaxation to ligand binding and migration in myoglobin

Protein relaxation, ligand binding, and ligand migration into a hydrophobic cavity in myoglobin are unified by a bounded diffusion model which produces an accurate fit to complex ligand rebinding data over eight decades in time and a 160 K temperature range, in qualitative agreement with time-resolved x-ray crystallography. Protein relaxation operates in a cyclic manner to move the ligand away from the binding site.

Different relaxations in myoglobin after photolysis

To clarify the interplay of kinetic hole-burning (KHB), structural relaxation, and ligand migration in myoglobin (Mb), we measured time-resolved absorption spectra in the Soret region after photolysis of carbon monoxide Mb (MbCO) in the temperature interval 120-260 K and in the time window 350 ns to 200 ms. The spectral contributions of both photolyzed (Mb*) and liganded Mb (MbCO) have been analyzed by taking into account homogeneous bandwidth, coupling to vibrational modes, and static conformational heterogeneity. We succeeded in separating the "time-dependent" spectral changes, and this work provides possibilities to identify the events in the process of ligand rebinding. KHB is dominant at T <190 K in both the Mb* and the MbCO components. For MbCO, conformational substates interconversion at higher temperatures tends to average out the KHB effect. At 230-260 K, whereas almost no shift is observed in the MbCO spectrum, a shift of the order of approximately 80 cm(-1) is observed in Mb*. We attribute this shift to protein relaxation coupled to ligand migration. The time dependence of the Mb* spectral shift is interpreted with a model that enables us to calculate the highly nonexponential relaxation kinetics. Fits of stretched exponentials to this kinetics yield Kohlrausch parameter values of 0.25, confirming the analogy between proteins and glasses.

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.

Multiple monovalent ion-dependent pathways for the folding of the L-21 Tetrahymena thermophila ribozyme

Synchrotron hydroxyl radical (*OH) footprinting is a technique that monitors the local changes in solvent accessibility of the RNA backbone on milliseconds to minutes time-scales. The Mg(2+)-dependent folding of the L-21 Sca 1 Tetrahymena thermophila ribozyme has been followed using this technique at an elevated concentration of monovalent ion (200 mM NaCl) and as a function of the initial annealing conditions and substrate. Previous studies conducted at low concentrations of monovalent ion displayed sequential folding of the P4-P6 domain, the peripheral helices and the catalytic core, with each protection displaying monophasic kinetics. For ribozyme annealed in buffer containing 200 mM NaCl and folded by the addition of 10 mM MgCl(2), multiple kinetic phases are observed for *OH protections throughout the ribozyme. The independently folding P4-P6 domain is the first to fold with its protections displaying 50-90% burst phase amplitudes. That the folding of P4-P6 within the ribozyme does not display the 100% burst phase of isolated P4-P6 at 200 mM NaCl shows that interactions with the remainder of the ribozyme impede this domain's folding. In addition, *OH protections constituting each side of a tertiary contact are not coincident in some cases, consistent with the formation of transient non-native interactions. While the peripheral contacts and triple helical scaffold exhibit substantial burst phases, the slowest protection to appear is J8/7 in the catalytic core, which displays a minimal burst amplitude and whose formation is coincident with the recovery of catalytic activity. The number of kinetic phases as well as their amplitudes and rates are different when the ribozyme is annealed in low-salt buffer and folded by the concomitant addition of monovalent and divalent cations. Annealed substrate changes the partitioning of the ribozyme among the multiple folding populations. These results provide a map of the early steps in the ribozyme's folding landscape and the degree to which the preferred pathways are dependent upon the initial reaction conditions.

Ligand binding and protein dynamics in neuroglobin

Neuroglobin (Ngb) is a recently discovered protein in vertebrate brain tissue that belongs to the globin family of proteins. It has been implicated in the neuronal response to hypoxia or ischemia, although its physiological role has been hitherto unknown. Ngb is hexacoordinate in the ferrous deoxy form under physiological conditions. To bind exogenous ligands like O(2) and CO, the His E7 endogenous ligand is displaced from the sixth coordination. By using infrared spectroscopy and nanosecond time-resolved visible spectroscopy, we have investigated the ligand-binding reaction over a wide temperature range (3-353 K). Multiple, intrinsically heterogeneous distal heme pocket conformations exist in NgbCO. Photolysis at cryogenic temperatures creates a five-coordinate deoxy species with very low geminate-rebinding barriers. The photodissociated CO is observed to migrate within the distal heme pocket even at 20 K. Flash photolysis near physiological temperature (275-353 K) exhibits four sequential kinetic features: (i) geminate rebinding (t < 1 micros); (ii) extremely fast bimolecular exogenous ligand binding (10 micros < t < 1 ms) with a nontrivial temperature dependence; (iii) endogenous ligand binding (100 micros < t < 10 ms), which can be studied by using flash photolysis on deoxy Ngb; and (iv) displacement of the endogenous by the exogenous ligand (10 ms < t < 10 ks). All four processes are markedly nonexponential, suggesting that Ngb fluctuates among different conformations on surprisingly long time scales.

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.

Geminate rebinding in trehalose-glass embedded myoglobins reveals residue-specific control of intramolecular trajectories

It is becoming increasingly apparent that hydrophobic cavities (also referred to as xenon cavities) within proteins have significant functional implications. The potential functional role of these cavities in modulating the internal dynamics of carbon monoxide in myoglobin (Mb) is explored in the present study by using glassy matrices derived from trehalose to limit protein dynamics and to eliminate ligand exchange between the solvent and the protein. By varying the temperature (-15 to 65 degrees C) and humidity for samples of carbonmonoxy myoglobin embedded in trehalose-glass, it is possible to observe a hierarchy of distinct geminate recombination phases that extend from nanosecond to almost seconds that can be directly associated with rebinding from specific hydrophobic cavities. The use of mutant forms of Mb reveals the role of key residues in modulating ligand access between these cavities and the distal hemepocket.

Trehalose effect on low temperature protein dynamics: fluctuation and relaxation phenomena

We performed spectral diffusion experiments in trehalose-enriched glycerol/buffer-glass on horseradish peroxidase where the heme was replaced by metal-free mesoporphyrin IX, and compared them with the respective behavior in a pure glycerol/buffer-glass (Schlichter et al., J. Chem. Phys. 2000, 112:3045-3050). Trehalose has a significant influence: spectral diffusion broadening speeds up compared to the trehalose-free glass. This speeding up is attributed to a shortening of the correlation time of the frequency fluctuations most probably by preventing water molecules from leaving the protein interior. Superimposed to the frequency fluctuation dynamics is a relaxation dynamics that manifests itself as an aging process in the spectral diffusion broadening. Although the trehalose environment speeds up the fluctuations, it does not have any influence on the relaxation. Both relaxation and fluctuations are governed by power laws in time. The respective exponents do not seem to change with the protein environment. From the spectral dynamics, the mean square displacement in conformation space can be determined. It is governed by anomalous diffusion. The associated frequency correlation time is incredibly long, demonstrating that proteins at low temperatures are truly nonergodic systems.

Sampling field heterogeneity at the heme of c-type cytochromes by spectral hole burning spectroscopy and electrostatic calculations

We report on a comparative investigation of the heme pocket fields of two Zn-substituted c-type cytochromes-namely yeast and horse heart cytochromes c-using a combination of hole burning Stark spectroscopy and electrostatic calculations. The spectral hole burning experiments are consistent with different pocket fields experienced at the hemes of the respective cytochromes. In the case of horse heart Zn-cytochrome c, two distinguishable electronic origins with different electrostatic properties are observed. The yeast species, on the other hand, displays a single electronic origin. Electrostatic calculations and graphics modeling using the linearized finite-difference Poisson-Boltzmann equation performed at selected time intervals on nanosecond-molecular dynamics trajectories show that the hemes of the respective cytochromes sample different potentials as they explore conformational space. The electrostatic potentials generated by the protein matrix at the heme show different patterns in both cytochromes, and we suggest that the cytochromes differ by the number of "electrostatic substates" that they can sample, thus accounting for the different spectral populations observed in the two cytochromes.

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.

Spectroscopic effects of polarity and hydration in the distal heme pocket of deoxymyoglobin

Distal pocket mutations at the E7 position (His64) of sperm whale deoxymyoglobin (deoxyMb) are used as a probe of distal pocket polarity and hydration. Changes of two key spectroscopic markers, the Fe-His(F8) stretch in the resonance Raman spectrum and the position of band III in the absorption spectrum, are monitored as the His64Tyr, His64Phe, His64Leu, and His64Gly mutations alter the distal heme pocket environment. The Fe-His vibration for the Phe, Leu, and Gly mutants is shifted to a lower frequency by 1-2 cm-1 relative to the Tyr mutant, wild type (WT), and native deoxyMb. Band III shifts to the red by approximately 4 nm ( approximately 70 cm-1) relative to WT and native deoxyMb for all the His64 mutants examined in this work. We correlate the small shift in the Fe-His frequency to the local electrostatic environment directly above the heme iron, affected by the presence of a localized water molecule in the heme pocket that is hydrogen-bonded to the E7 residue. The position of band III is roughly correlated to the displacement of the iron from the heme plane; however, the relatively large scatter in this correlation, as well as its dependence on distal pocket mutations, suggests that the heme pocket environment, particularly the E7 residue, also affects the energy of this transition.