Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.     

Primary processes in heme-based sensor proteins

A wide and still rapidly increasing range of heme-based sensor proteins has been discovered over the last two decades. At the molecular level, these proteins function as bistable switches in which the catalytic activity of an enzymatic domain is altered mostly by binding or dissociation of small gaseous ligands (O₂, NO or CO) to the heme in a sensor domain. The initial “signal” at the heme level is subsequently transmitted within the protein to the catalytic site, ultimately leading to adapted expression levels of specific proteins. Making use of the photolability of the heme-ligand bond that mimics thermal dissociation, early processes in this intra-protein signaling pathway can be followed using ultrafast optical spectroscopic techniques; they also occur on timescales accessible to molecular dynamics simulations. Experimental studies performed over the last decade on proteins including the sensors FixL (O₂), CooA (CO) and soluble guanylate cyclase (NO) are reviewed with an emphasis on emerging general mechanisms. After heme-ligand bond breaking, the ligand can escape from the heme pocket and eventually from the protein, or rebind directly to the heme. Remarkably, in all sensor proteins the rebinding, specifically of the sensed ligand, is highly efficient. This ”ligand trap” property possibly provides means to smoothen the effects of fast environmental fluctuations on the switching frequency. For 6-coordinate proteins, where exchange between an internal heme-bound residue and external gaseous ligands occurs, the study of early processes starting from the unliganded form indicates that mobility of the internal ligand may facilitate signal transfer. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

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.     

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Ligand preference and orientation in b- and c-type heme-binding proteins

Hemes are often incorporated into designed proteins. The importance of the heme ligand type and its orientation is still a matter of debate. Here, heme ligands and ligand orientation were investigated using a nonredundant (87 structures) and a redundant (1503 structures) set of structures to compare and contrast design features of natural b- and c-type heme-binding proteins. Histidine is the most common ligand. Marked differences in ligation motifs between b- and c-type hemes are higher occurrence of His-Met in c-type heme binding motifs (16.4% vs. 1.4%) and higher occurrence of exchangeable, small molecules in b-type heme binding motifs (67.6% vs. 9.9%). Histidine ligands that are part of the c-type CXXCH heme-binding motif show a distinct asymmetric distribution of orientation. They tend to point between either the heme propionates or between the NA and NB heme nitrogens. Molecular mechanics calculations show that this asymmetry is due to the bonded constraints of the covalent attachment between the heme and the protein. In contrast, the orientations of b-type hemes histidine ligands are found evenly distributed with no preference. Observed histidine heme ligand orientations show no dominating influence of electrostatic interactions between the heme propionates and the ligands. Furthermore, ligands in bis-His hemes are found more frequently perpendicular rather than parallel to each other. These correlations support energetic constraints on ligands that can be used in designing proteins.     

Ultrafast heme-residue bond formation in six-coordinate heme proteins: implications for functional ligand exchange

A survey is presented of picosecond kinetics of heme-residue bond formation after photolysis of histidine, methionine, or cysteine, in a broad range of ferrous six-coordinate heme proteins. These include human neuroglobin, a bacterial heme-binding superoxide dismutase (SOD), plant cytochrome b 559, the insect nuclear receptor E75, horse heart cytochrome c and the heme domain of the bacterial sensor protein Dos. We demonstrate that the fastest and dominant phase of binding of amino acid residues to domed heme invariably takes place with a time constant in the narrow range of 5-7 ps. Remarkably, this is also the case in the heme-binding SOD, where the heme is solvent-exposed. We reason that this fast phase corresponds to barrierless formation of the heme-residue bond from a configuration close to the bound state. Only in proteins where functional ligand exchange occurs, additional slower rebinding takes place on the time scale of tens of picoseconds after residue dissociation. We propose that the presence of these slower phases reflects flexibility in the heme environment that allows external ligands (O2, CO, NO, . . .) to functionally replace the internal residue after thermal dissociation of the heme-residue bond.     

Mechanisms of ligand discrimination by heme proteins

Several classes of heme proteins have been identified whose primary role is to sense and transport one of the three gaseous ligands: O(2), NO, or CO. A common feature shared by all of these proteins is the need to not only recognize its target ligand, but also discriminate against other heme ligands of similar size and shape. This review describes the mechanisms that each class of heme protein sensor and transporter utilizes to promote this discrimination. The common factors utilized in the recognition of each ligand are discussed.     

Investigations of the myoglobin cavity mutant H93G with unnatural imidazole proximal ligands as a modular peroxide O-O bond cleavage model system

A general inability to elucidate extensive variations in the electronic characteristics of proximal heme iron ligands in heme proteins has hampered efforts to obtain a clear understanding of the role of the proximal heme iron ligand in the activation of oxygen and peroxide. The disadvantage of the frequently applied site-directed mutagenesis technique is that it is limited by the range of natural ligands available within the genetic code. The myoglobin cavity mutant H93G [Barrick, D. (1994) Biochemistry 33, 6546-6554] has its proximal histidine ligand replaced with glycine, a mutation which leaves an open cavity capable of accommodating a variety of unnatural potential proximal ligands. We have carried out investigations of the effect of changing the electron donor characteristics of a variety of substituted imidazole proximal ligands on the rate of formation of myoglobin compound II and identified a correlation between the substituted imidazole N-3 pK(a) (which provides a measure of the electron donor ability of N-3) and the apparent rate of formation of compound II. A similar rate dependence correlation is not observed upon binding of azide. This finding indicates that O-O bond cleavage and not the preceding peroxide binding step is being influenced by the electron donor characteristics of the substituted imidazole ligands. The proximal ligand effects are clearly visible, but their overall magnitude is quite low (1.7-fold increase in the O-O bond cleavage rate per pK(a) unit). This appears to provide support for recent commentaries which concluded that the partial ionization of the proximal histidine ligand in typical heme peroxidases may not be enough of an influence to provide a mechanistically critical push effect [Poulos, T. L. (1996) JBIC, J. Biol. Inorg. Chem. 1, 356-359]. Further attempts were made to define the mechanism of the influence of N-3 pK(a) on O-O bond cleavage by using peracetic acid and cumene hydroperoxide as mechanistic probes. The observation of heme destruction in these reactions indicates that displacement of the proximal imidazole ligands by peracetic acid or cumene hydroperoxide has occurred. A combination mutation (H64D/H93G) was prepared with the objective of observing compound I of H64D/H93G with substituted imidazoles as proximal ligands upon reaction with H(2)O(2). This double mutant was found to simultaneously bind imidazole to both axial positions, an arrangement which prevents a reaction with H(2)O(2).     

Covalent heme attachment in Synechocystis hemoglobin is required to prevent ferrous heme dissociation

Synechocystis hemoglobin contains an unprecedented covalent bond between a nonaxial histidine side chain (H117) and the heme 2-vinyl. This bond has been previously shown to stabilize the ferric protein against denaturation, and also to affect the kinetics of cyanide association. However, it is unclear why Synechocystis hemoglobin would require the additional degree of stabilization accompanying the His117-heme 2-vinyl bond because it also displays endogenous bis-histidyl axial heme coordination, which should greatly assist heme retention. Furthermore, the mechanism by which the His117-heme 2-vinyl bond affects ligand binding has not been reported, nor has any investigation of the role of this bond on the structure and function of the protein in the ferrous oxidation state. Here we report an investigation of the role of the Synechocystis hemoglobin His117-heme 2-vinyl bond on structure, heme coordination, exogenous ligand binding, and stability in both the ferrous and ferric oxidation states. Our results reveal that hexacoordinate Synechocystis hemoglobin lacking this bond is less stable in the ferrous oxidation state than the ferric, which is surprising in light of our understanding of pentacoordinate Hb stability, in which the ferric protein is always less stable. It is also demonstrated that removal of the His117-heme 2-vinyl bond increases the affinity constant for intramolecular histidine coordination in the ferric oxidation state, thus presenting greater competition for the ligand binding site and lowering the observed rate and affinity constants for exogenous ligands.     

Solution 1H NMR investigation of the seating and rotational "hopping" of centrosymmetric etioheme-I in myoglobin: effect of globin origin and its oxidation/spin state on heme dynamics

Solution 1H NMR spectroscopy was used to investigate the heme active-site structure and dynamics of rotation about the Fe-His bond of centrosymmetric etioheme-I reconstituted into sperm whale and horse myoglobin (Mb). Comparison of the NOESY cross-peak pattern and paramagnetic relaxation properties of the cyanomet complexes confirm a heme pocket that is essentially the same as Mb with either native protoheme or etioheme-I. Dipolar contacts between etioheme and the conserved heme pocket residues establish a unique seating of etioheme that conserves the orientation of the N-Fe-N vector relative to the axial His plane, with ethyl groups occupying the vinyl positions of protoheme. Saturation transfer between methyls on adjacent pyrroles in etioheme-reconstituted horse Mb in all accessible oxidation/spin states reveals rotational hopping rates that decrease dramatically with either loss of ligands or reduction of the heme, and correlate qualitatively with expectations based on the Fe-His bond strength and the rate of heme dissociation from Mb. The rate of hopping for etioheme in metMbCN, in contrast to hemes with propionates, is the same in the sperm whale and horse proteins.     

Electrochemical and ligand binding studies of a de novo heme protein

Heme proteins can perform a variety of electrochemical functions. While natural heme proteins carry out particular functions selected by biological evolution, artificial heme proteins, in principle, can be tailored to suit specified technological applications. Here we describe initial characterization of the electrochemical properties of a de novo heme protein, S824C. Protein S824C is a four-helix bundle derived from a library of sequences that was designed by binary patterning of polar and nonpolar amino acids. Protein S824C was immobilized on a gold electrode and the formal potential of heme-protein complex was studied as a function of pH and ionic strength. The binding of exogenous N-donor ligands to heme/S824C was monitored by measuring shifts in the potential that occurred upon addition of various concentrations of imidazole or pyridine derivatives. The response of heme/S824C to these ligands was then compared to the response of isolated heme (without protein) to the same ligands. The observed shifts in potential depended on both the concentration and the structure of the added ligand. Small changes in structure of the ligand (e.g. pyridine versus 2-amino pyridine) produced significant shifts in the potential of the heme-protein. The observed shifts correlate to the differential binding of the N-donor molecules to the oxidized and reduced states of the heme. Further, it was observed that the electrochemical response of the buried heme in heme/S824C differed significantly from that of isolated heme. These studies demonstrate that the structure of the de novo protein modulates the binding of N-donor ligands to heme.     

Picosecond kinetics of cytochromes b5 and c

Ligand photolysis and subsequent recombination in cytochromes b5 and c have been studied with picosecond resolution. In both proteins, an iron-histidine bond is broken after excitation with 314-nm light, and recombination occurs with a rate constant of about 1.4 x 10(11) s-1. Photolysis and reformation of the iron-histidine bond may be surprising as these hemoproteins do not reversibly bind ligands in nature. The findings are explained using results both from experiments on model hemes and from computer investigations with atomic resolution on the three-dimensional structure of the protein. After photolysis, the formation and recombination of the geminate contact pair are attributed to simple low amplitude ligand bond rotations, a result that can be applied to geminate processes in other hemoproteins and model heme compounds as well.     

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.     

Time-resolved methods in Biophysics. 2. Monitoring haem proteins at work with nanosecond laser flash photolysis.

Haem proteins have long been the most studied proteins in biophysics, and have become paradigms for the characterization of fundamental biomolecular processes as ligand binding and regulatory conformational transitions. The presence of the haem prosthetic group, the absorbance spectrum of which has a ligation sensitive region conveniently located in the UV-visible range, has offered a powerful and sensitive tool for the investigation of molecular functions. The central Fe atom is capable of reversibly binding diatomic ligands, including O(2), CO, and NO. The Fe-ligand bond is photolabile, and a reactive unligated state can be transiently generated with a pulsed laser. The photodissociated ligands quickly rebind to the haem and the process can be monitored by transient absorbance methods. The ligand rebinding kinetics reflects protein dynamics and ligand migration within the protein inner cavities. The characterization of these processes was done in the past mainly by low temperature experiments. The use of silica gels to trap proteins allows the characterization of internal ligand dynamics at room temperature. In order to show the potential of the laser flash photolysis techniques, combined with modern numerical analysis methods, we report experiments conducted on two non-symbiotic haemoglobins from Arabidopsis thaliana. The comparison between time courses recorded on haemoglobins in solution and encapsulated in silica gels allows for the highlighting of different interplays between protein dynamics and ligand migration.     

Role of distal arginine in early sensing intermediates in the heme domain of the oxygen sensor FixL

FixL is a bacterial heme-based oxygen sensor, in which release of oxygen from the sensing PAS domain leads to activation of an associated kinase domain. Static structural studies have suggested an important role of the conserved residue arginine 220 in signal transmission at the level of the heme domain. To assess the role of this residue in the dynamics and properties of the initial intermediates in ligand release, we have investigated the effects of R220X (X = I, Q, E, H, or A) mutations in the FixLH heme domain on the dynamics and spectral properties of the heme upon photolysis of O(2), NO, and CO using femtosecond transient absorption spectroscopy. Comparison of transient spectra for CO and NO dissociation with steady-state spectra indicated less strain on the heme in the ligand dissociation species for all mutants compared to the wild type (WT). For CO and NO, the kinetics were similar to those of the wild type, with the exception of (1) a relatively low yield of picosecond NO rebinding to R220A, presumably related to the increase in the free volume of the heme pocket, and (2) substantial pH-dependent picosecond to nanosecond rebinding of CO to R220H, related to formation of a hydrogen bond between CO and histidine 220. Upon excitation of the complex bound with the physiological sensor ligand O(2), a 5-8 ps decay phase and a nondecaying (>4 ns) phase were observed for WT and all mutants. The strong distortion of the spectrum associated with the decay phase in WT is substantially diminished in all mutant proteins, indicating an R220-induced role of the heme in the primary intermediate in signal transmission. Furthermore, the yield of dissociated oxygen after this phase ( approximately 10% in WT) is increased in all mutants, up to almost unity in R220A, indicating a key role of R220 in caging the oxygen near the heme through hydrogen bonding. Molecular dynamics simulations corroborate these findings and suggest motions of O(2) and arginine 220 away from the heme pocket as a second step in the signal pathway on the 50 ps time scale.     

Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

Globins are respiratory proteins that reversibly bind dioxygen and other small ligands at the iron of a heme prosthetic group. Hemoglobin and myoglobin are the most prominent members of this protein family. Unexpectedly a few years ago a new member was discovered and called neuroglobin (Ngb), being predominantly expressed in the brain. Ngb is a single polypeptide of 151 amino acids and despite the small sequence similarity with other globins, it displays the typical globin fold. Oxygen, nitric oxide, or carbon monoxide can displace the distal histidine which, in ferrous Ngb as well as in ferric Ngb, is bound to the iron, yielding a reversible adduct. Recent crystallographic data on carboxy Ngb show that binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities; in particular, the huge internal tunnel that connects the bulk with the active site, peculiar to Ngb, is heavily reorganized. We report the results of extended (90 ns) molecular dynamics simulations in water of ferrous deoxy and carboxy murine neuroglobin, which are both coordinated on the distal site, in the latter case by CO and in the former one by the distal His(64)(E7). The long timescale of the simulations allowed us to characterize the equilibrated protein dynamics and to compare protein structure and dynamical behavior coupled to the binding of an exogenous ligand. We have characterized the heme sliding motion, the topological reorganization of the internal cavities, the dynamics of the distal histidine, and particularly the conformational change of the CD loop, whose flexibility depends ligand binding.     

Can the axial ligand strength be monitored through spectroscopic measurements?

 The influence of the properties of the axial ligands in heme proteins on the electronic structure of the iron(III) ion has been studied through angular overlap calculations. In particular, a correlation between (a) ligand field parameters and spectroscopic data and (b) the reduction potentials of the iron ion is proposed. The results of this approach are discussed with respect to their relevance for the biological function of different heme proteins. Received and accepted: 7 May 1996     

The reactivity of cytochrome c with soft ligands.

The spectral changes caused by binding soft ligands to the cytochrome c iron and their correlation to ligand affinities support the hypothesis that the iron—methionine sulfur bond of this heme protein is enhanced by delocalization of the metal l₂, electrons into the empty 3d orbitals of the ligand atom. These findings also explain the unique spectrum of cytochrome c in the far red.