Ligand binding to heme proteins: II. Transitions in the heme pocket of myoglobin.

Phenomena occurring in the heme pocket after photolysis of carbonmonoxymyoglobin (MbCO) below about 100 K are investigated using temperature-derivative spectroscopy of the infrared absorption bands of CO. MbCO exists in three conformations (A substrates) that are distinguished by the stretch bands of the bound CO. We establish connections among the A substates and the substates of the photoproduct (B substates) using Fourier-transform infrared spectroscopy together with kinetic experiments on MbCO solution samples at different pH and on orthorhombic crystals. There is no one-to-one mapping between the A and B substates; in some cases, more than one B substate corresponds to a particular A substate. Rebinding is not simply a reversal of dissociation; transitions between B substates occur before rebinding. We measure the nonequilibrium populations of the B substates after photolysis below 25 K and determine the kinetics of B substate transitions leading to equilibrium. Transitions between B substates occur even at 4 K, whereas those between A substates have only been observed above about 160 K. The transitions between the B substates are nonexponential in time, providing evidence for a distribution of substates. The temperature dependence of the B substate transitions implies that they occur mainly by quantum-mechanical tunneling below 10 K. Taken together, the observations suggest that the transitions between the B substates within the same A substate reflect motions of the CO in the heme pocket and not conformational changes. Geminate rebinding of CO to Mb, monitored in the Soret band, depends on pH. Observation of geminate rebinding to the A substates in the infrared indicates that the pH dependence results from a population shift among the substates and not from a change of the rebinding to an individual A substate.     

Resonance raman investigations of site-directed mutants of myoglobin: effects of distal histidine replacement

The resonance Raman spectra of met-, deoxy-, and (carbonmonoxy)myoglobin (MbCO) are studied as a function of amino acid replacement at the distal histidine-E7 position. The synthetic wild type is found to be spectroscopically identical with the native material. The methionine and glycine replacements do not affect the met or deoxy spectra but do lead to distinct changes in the nu Fe-CO region of the MbCO spectrum. The native MbCO displays a pH-dependent population redistribution of the nu Fe-CO modes, while the analogous population in the mutant systems is found to be pH independent. This indicates that histidine-E7 is the titratable group in native MbCO. Moreover, the pH dependence of the population dynamics is found to be inconsistent with a simple two-state Henderson-Hasselbalch analysis. Instead, we suggest a four-state model involving the coupling of histidine protonation and conformational change. Within this model, the pK of the distal histidine is found to be 6.0 in the "open" configuration and 3.8 in the "closed" conformation. This corresponds to a 3 kcal/mol destabilization of the positively charged distal histidine within the hydrophobic pocket and suggests how protonation can lead to a larger population of the "open" conformation. At pH 7, the pocket is found to be "open" approximately 3% of the time. Further work, involving both IR and Raman measurements, allows the electron-nuclear coupling strengths of the various nu Fe-CO and nu C-O Raman modes to be determined. The slowly rebinding conformational state, corresponding to nu Fe-CO = 518 cm-1 (nu C-O = 1932 cm-1), displays unusually weak coupling of the Fe-CO mode to the Soret transition. Studies of the nu Fe-CO region as a function of temperature reveal that the equilibria between the conformational states are quenched in both the native and glycine mutant below the freezing point of the solvent. Unusual line narrowing of the nu Fe-CO modes at the phase transition is also observed in all samples studied. This line narrowing stands in marked contrast to the other heme Raman modes and suggests that Fe-CO librational motion and/or distal pocket vibrational (or conformational) excitations are involved in the line broadening at room temperature.     

Crystal structures of myoglobin-ligand complexes at near-atomic resolution.

We have used x-ray crystallography to determine the structures of sperm whale myoglobin (Mb) in four different ligation states (unligated, ferric aquomet, oxygenated, and carbonmonoxygenated) to a resolution of better than 1.2 A. Data collection and analysis were performed in as much the same way as possible to reduce model bias in differences between structures. The structural differences among the ligation states are much smaller than previously estimated, with differences of <0.25 A root-mean-square deviation among all atoms. One structural parameter previously thought to vary among the ligation states, the proximal histidine (His-93) azimuthal angle, is nearly identical in all the ferrous complexes, although the tilt of the proximal histidine is different in the unligated form. There are significant differences, however, in the heme geometry, in the position of the heme in the pocket, and in the distal histidine (His-64) conformations. In the CO complex the majority conformation of ligand is at an angle of 18 +/- 3 degrees with respect to the heme plane, with a geometry similar to that seen in encumbered model compounds; this angle is significantly smaller than reported previously by crystallographic studies on monoclinic Mb crystals, but still significantly larger than observed by photoselection. The distal histidine in unligated Mb and in the dioxygenated complex is best described as having two conformations. Two similar conformations are observed in MbCO, in addition to another conformation that has been seen previously in low-pH structures where His-64 is doubly protonated. We suggest that these conformations of the distal histidine correspond to the different conformational substates of MbCO and MbO(2) seen in vibrational spectra. Full-matrix refinement provides uncertainty estimates of important structural parameters. Anisotropic refinement yields information about correlated disorder of atoms; we find that the proximal (F) helix and heme move approximately as rigid bodies, but that the distal (E) helix does not.     

Myoglobin-CO conformational substate dynamics: 2D vibrational echoes and MD simulations

Two-dimensional (2D) infrared vibrational echoes were performed on horse heart carbonmonoxymyoglobin (MbCO) in water over a range of temperatures. The A(1) and A(3) conformational substates of MbCO are found to have different dephasing rates with different temperature dependences. A frequency-frequency correlation function derived from molecular dynamics simulations on MbCO at 298 K is used to calculate the vibrational echo decay. The calculated decay shows substantial agreement with the experimentally measured decays. The 2D vibrational echo probes protein dynamics and provides an observable that can be used to test structural assignments for the MbCO conformational substates.     

Investigations of photolysis and rebinding kinetics in myoglobin using proximal ligand replacements

We use laser flash photolysis and time-resolved Raman spectroscopy of CO-bound H93G myoglobin (Mb) mutants to study the influence of the proximal ligand on the CO rebinding kinetics. In H93G mutants, where the proximal linkage with the protein is eliminated and the heme can bind exogenous ligands (e.g., imidazole, 4-bromoimidazole, pyridine, or dibromopyridine), we observe significant effects on the CO rebinding kinetics in the 10 ns to 10 ms time window. Resonance Raman spectra of the various H93G Mb complexes are also presented to aid in the interpretation of the kinetic results. For CO-bound H93G(dibromopyridine), we observe a rapid large-amplitude geminate phase with a fundamental CO rebinding rate that is approximately 45 times faster than for wild-type MbCO at 293 K. The absence of an iron proximal ligand vibrational mode in the 10 ns photoproduct Raman spectrum of CO-bound H93G(dibromopyridine) supports the hypothesis that proximal ligation has a significant influence on the kinetics of diatomic ligand binding to the heme.     

Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from M?ssbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.     

Alkali cation effect on carbonyl-hemoglobin's and -myoglobin's conformer populations when exposed to freeze-concentration of their phosphate-buffered aqueous solutions

Freeze-concentrated aqueous phosphate-buffered (pH 6.8) solutions of carbonyl-hemoglobin (HbCO) and -myoglobin (MbCO) were investigated by Fourier-transform infrared spectroscopy for the effect of alkali cation on the population of conformers. When using sodium phosphates as buffer components, HbCO was transformed from conformer III (at approximately 1951 cm-1) which is the dominant form at ambient temperatures, into conformer IV (at buffer concentration at a given temperature. The conformational changes started slightly below the temperature where ice began to crystallize and the remaining solution became freeze-concentrated, and they were reversible for HbCO. For MbCO in 0.5 M sodium phosphate buffer solution, however, they were irreversible and MbCO denatured completely. When potassium phosphate salts were used for preparing the buffer at the same pH of 6.8, little or no transformation of conformer III into conformer IV was observed. The conformational changes induced by sodium salts are attributed to a decrease in pH and it is shown by infrared spectroscopy that during freeze concentration drastic changes in composition of the two buffer components H2PO4-/HPO(2)4- occur, the acid component increasing strongly relative to the base component. Supersaturation is also important because change from conformer III to IV requires a minimum concentration of sodium salts: whereas 0.1 M sodium phosphate buffer concentration shows a strong effect, 0.03 M concentration does not and therefore behaves like a potassium phosphate buffer.     

Low pH myoglobin photoproducts.

Recently, there has been interest in determining the conditions under which the iron-histidine bond ruptures in myoglobin at low pH, so that the effect of proximal heme ligation can be studied. A 220-cm-1 Raman mode, assigned to iron-histidine stretching, is clearly visible after photolysis of aqueous MbCO samples below pH4 at room temperature (Sage et al. Biochemistry. 30:1237-1247). In contrast, Iben et al. (Biophys. J. 59:908-919) do not observe this mode upon photolysis of a pH3 MbCO sample in a glycerol/water glass at low temperature. In order to account for both the low temperature and the room temperature experiments, Iben et al. suggest a scheme involving an unusual protonation state of the proximal histidine. Here, we discuss some inconsistencies in their explanation of the room temperature results and offer instead a simple modification of an earlier model. In addition, circular dichroism data are presented that indicate partial unfolding of MbCO in aqueous solution below pH4, and raise questions about the claim of Iben et al. that MbCO remains folded in 75% glycerol at pH3.     

Co-rebinding in myoglobin as seen by time-resolved X-ray absorption spectroscopy

The dynamics of selected conformational coordinates, key roles in the understanding of the CO-rebinding process, are investigated in horse heart carbonmonoxy myoglobin (MbCO) through time-resolved X-ray absorption spectroscopy. We present here the results obtained at 90 K in the second time scale. The approach of the CO molecule towards the Fe atom in the active site pocket is speculated to act as a natural precursor to the Fe displacement with the consequent undoming of the protein porphyrin plane. The arrangement of the Fe-C-O bonding angle geometry follows and the final MbCO active site configuration is completely reached within 1 min.     

Pressure effects on the proximal heme pocket in myoglobin probed by Raman and near-infrared absorption spectroscopy.

The influence of high pressure on the heme protein conformation of myoglobin in different ligation states is studied using Raman spectroscopy over the temperature range from 30 to 295 K. Photostationary experiments monitoring the oxidation state marker bands demonstrate the change of rebinding rate with pressure. While frequency changes of vibrational modes associated with rigid bonds of the porphyrin ring are <1 cm(-1), we investigate a significant shift of the iron-histidine mode to higher frequency with increasing pressure (approximately 3 cm(-1) for deltaP = 190 MPa in Mb). The observed frequency shift is interpreted structurally as a conformational change affecting the tilt angle between the heme plane and the proximal histidine and the out-of-plane iron position. Independent evidence for iron motion comes from measurements of the redshift of band III in the near-infrared with pressure. This suggests that at high pressure the proximal heme pocket and the protein are altered toward the bound state conformation, which contributes to the rate increase for CO binding. Raman spectra of Mb and photodissociated MbCO measured at low temperature and variable pressure further support changes in protein conformation and are consistent with glasslike properties of myoglobin below 160 K.     

FTIR spectroscopic study of the dynamics of conformational substates in hydrated carbonyl-myoglobin films via temperature dependence of the CO stretching band parameters.

Two hydrated carbonyl myoglobin (MbCO) films, one containing (0.30 g water)/(g MbCO) from MbCO solution in water at pH 5.5 and the other (0.32 g water)/(gMbCO) from 0.1 M potassium phosphate buffer solution at pH 6.8, were studied by FTIR spectroscopy from 293 K to 78 K at selected temperatures on cooling and reheating. Above approximately 180 K the general trend in temperature dependence of half-bandwidths, peak maxima, and band area ratios of the A1 and A3 conformer bands is similar to those reported by Ansari et al. (1987. Biophys. J. 26:337) for MbCO in 75% glycerol/water solution, but abrupt changes in slopes at approximately 180-200 K and freezing-in of conformer populations, which could be taken as indicator for glass transition of the solvent or the protein, are absent for the hydrated MbCO films. This is interpreted in terms of an exceptionally broad distribution of relaxation times, and is in accord with conclusions from recent calorimetric annealing studies of hydrated protein powders (Sartor et al. 1994. Biophys. J. 66:249). Exchange between the three A conformers does not stop at approximately 180-200 K but occurs over the whole temperature region studied. These results are then discussed with respect to MbCO's behavior in the glass-->liquid transition region of glass-forming solvents, and it is concluded that, in analogy to the behavior of low-molecular-weight compounds with a distribution of rapidly interconverting conformers, freezing-in of MbCO's A conformer populations by the solvent should not be mistaken for a glass transition of MbCO.     

Temperature dependence of the structure and dynamics of myoglobin. A simulation approach

The results of simulations of the structure and internal motions of carbonomonoxymyoglobin (MbCO) at two different temperatures (325 and 80 K) are presented and compared with experimental data. Properties calculated from the 120 ps trajectory at 325 K are used as a reference in the analysis of the motion of the protein at 80 K. Three separate 80 K molecular dynamics trajectories were calculated; they were started with different coordinate sets from the 325 K simulation and the lower temperature was achieved by scaling the velocities. The simulations yield results for the structural changes between 325 and 80 K that are in general accord with those from X-ray data. Both the experimental and calculated radii of gyration, distances from the center of mass and main-chain difference distance matrices show that there is a significant but inhomogeneous shrinkage with decreasing temperature. For the atomic fluctuations, by contrast, the calculated temperature dependence is very different from the X-ray results; i.e. the calculated root-mean-square backbone fluctuations decrease to 0.11 A at 80 K from 0.51 A at 325 K, while the fluctuations obtained from the X-ray B factors go from 0.56 A at 260 K to 0.47 A at 80 K. The smaller temperature dependence of the B factors suggests that there is significant conformational disorder in MbCO crystals at lower temperatures. This is in accord with the simulation results, which show that the protein is trapped in restricted regions of conformational space at 80 K, while at 325 K a much larger region is accessible to the protein. Analysis of the fluctuations at 325 K and 80 K shows that the room temperature flexibility of the protein is determined by the mobility of the loop regions and by side-chain torsional motions (in accord with earlier simulation results), while the low temperature fluctuations involve motion within a single well. Examination of the calculated iron atom fluctuations and comparison with Mossbauer data show good agreement. It is found that the dominant contribution to the iron motion arises from heme sliding; motion of the iron relative to the heme are much smaller.     

Single-site mutations in a hyperthermophilic variant of the B1 domain of protein G result in self-assembled oligomers

We have characterized two homologous, single-point core mutants of a 57-residue, hyperthermophilic variant of the B1 domain of protein G (HTB1). These single-point mutations in HTB1 replace a Phe residue in the hydrophobic core with either a Glu or Asp residue. Both of these homologous core-variant mutants undergo significant structural rearrangement from the native, monomeric fold and exist as stable soluble oligomeric species of 5 and 30 nm in diameter. Gel-filtration, dynamic light scattering, circular dichroism spectroscopy, fluorescence spectroscopy, along with Congo Red and Thioflavin T binding clearly demonstrated that these core-variants undergo significant structural rearrangement from the native, monomeric ubiquitin fold. The two oligomeric species did not equilibrate over extended periods of time and displayed distinct secondary structures. The larger of the two species was found to possess structural features that are reminiscent of an emerging class of protein assemblies prone to beta-sheet-mediated aggregation. These results are significant as there are very few examples of extensive conformational or oligomerization switching brought about by single-point mutations in a stable protein-fold.     

Carboxy Mb at pH 3. Time-resolved resonance Raman study at cryogenic temperatures.

Cryogenic samples of MbCO at pH3 are studied using nanosecond and picosecond time-resolved resonance Raman spectroscopy. It is observed that under excitation conditions sufficient to completely photodissociate MbCO at pH7, the pH3 sample at 10 ns remains substantially unphotolyzed even at 15 K. The similarity in the optical and resonance Raman spectra of MbCO at pH3 with that of pH7 indicates that at pH3 the iron remains six-coordinate and low-spin. The Fe-CO stretch frequency is consistent with a more upright CO orientation. The absence of the v(Fe-His) band in the 30 ps photoproduct Raman spectrum suggests that the Fe-His(F8) bond is broken within 30 ps of photodissociation. Other Raman bands, though, are not consistent with a normal four-coordinate heme for the photoproduct, Mb*. Suggested possible interpretations include a four-coordinate heme highly perturbed by the close lying protonated proximal histidine or a five-coordinate heme with the Fe-His bond significantly weakened. The partial photolysis monitored at 30 ps and 100 K indicates either a significant amount of geminate recombination within 30 ps or low quantum yield or photolysis. The time course for CO recombination is monitored via the Raman spectra from 30 ps to 3 ns at 100 K and 160 K. Of the fraction of protein-ligand pairs that remain photodissociated at 30 ps, 50% recombine by approximately 250 ps at 100 K and 160 K, supporting the flash photolysis rebinding data of Cowen et al. (Cowen, B. R. 1990. Ph. D. thesis. University of Illinois at Urbana-Champaign; Cowen, B. R., D. Braunstein, H. Frauenfelder, P. J. Steinbach, and R. D. Young. 1989. Biophys. J. 55:55a. [Abstr.].) The conclusions from these resonance Raman studies are extended to solution phase studies at ambient temperatures.     

Mechanism of chaperone function in small heat shock proteins. Two-mode binding of the excited states of T4 lysozyme mutants by alphaA-crystallin

To elucidate the mechanism of alphaA-crystallin chaperone function, a detailed thermodynamic analysis of its binding to destabilized, site-directed mutants of T4 lysozyme was carried out. The selected mutants form a ladder of stabilities spanning the 5-10 kcal/mol range of free energy of unfolding. The crystal structures of the majority of the mutants have been previously determined and found to be similar to that of the wild type with no evidence of static local unfolding. Complex formation between alphaA-crystallin and T4 lysozyme was observed directly via the changes in the electron paramagnetic resonance lineshape of a nitroxide introduced at a non-destabilizing, solvent exposed site in T4 lysozyme. AlphaA-crystallin differentially interacts with the mutants, binding the more destabilized ones to a larger extent despite the similar structure of their native states. Our results suggest that the states recognized by alphaA-crystallin are non-native excited states distinct from the unfolded state. Stable complexes are formed when the free energy of binding to alphaA-crystallin is on the order of the free energy associated with the transition from the excited state to the native state. Biphasic binding isotherms reveal two modes of interactions with distinct affinities and stoichiometries. Highly destabilized mutants preferentially bind to the high capacity mode, suggesting conformational preference in the use of each mode. Furthermore, binding can be enhanced by increased temperature and pH, which may be reflecting conformational changes in alphaA-crystallin oligomeric structure.     

Structural dynamics of myoglobin: FTIR-TDS study of NO migration and binding

By using Fourier transform infrared photolysis difference spectroscopy combined with temperature derivative spectroscopy at cryogenic temperatures, we have measured infrared spectra of the stretching absorption on nitric oxide (NO) in the heme-bound and photodissociated states of ferrous and ferric nitrosyl myoglobin (MbNO) and a few site-specific Mb mutants. The NO absorption was utilized as a sensitive local probe of ligand interactions with active-site residues and movements within the protein. By comparison with results obtained in previous spectroscopic and structural studies of carbonmonoxy myoglobin (MbCO), the MbNO data were interpreted in structural terms. In the NO-bound state, conformational heterogeneity was inferred from the appearance of multiple bands, arising from different electrostatic interactions with active site residues, most importantly, His-64. In ferrous MbNO, a primary photoproduct site similar to site B of MbCO was found, as indicated by a characteristic NO stretching spectrum. In ferric MbNO, the His-64 side chain appears to interfere with trapping of NO in this site; only a very weak photoproduct spectrum was observed in Mb variants in which His-64 was present. Upon extended illumination, the photoproduct spectrum changed in a characteristic way, indicating that NO readily migrates to a secondary docking site C, the Xe4 cavity, in which the ligand performs librational motions on the picosecond time scale. This docking site may play a role in the physiological NO scavenging reaction. Surprisingly, NO cannot be trapped at all in secondary docking site D, the Xe1 cavity.     

Effective light-induced hydroxylamine reactions occur with C13 = C14 nonisomerizable bacteriorhodopsin pigments.

The light-driven proton pump bacteriorhodopsin (bR) undergoes a bleaching reaction with hydroxylamine in the dark, which is markedly catalyzed by light. The reaction involves cleavage of the (protonated) Schiff base bond, which links the retinyl chromophore to the protein. The catalytic light effect is currently attributed to the conformational changes associated with the photocycle of all-trans bR, which is responsible for its proton pump mechanism and is initiated by the all-trans --> 13-cis isomerization. This hypothesis is now being tested in a series of experiments, at various temperatures, using three artificial bR molecules in which the essential C13==C14 bond is locked by a rigid ring structure into an all-trans or 13-cis configuration. In all three cases we observe an enhancement of the reaction by light despite the fact that, because of locking of the C13==C14 bond, these molecules do not exhibit a photocycle, or any proton-pump activity. An analysis of the rate parameters excludes the possibility that the light-catalyzed reaction takes place during the approximately 20-ps excited state lifetimes of the locked pigments. It is concluded that the reaction is associated with a relatively long-lived (micros-ms) light-induced conformational change that is not reflected by changes in the optical spectrum of the retinyl chromophore. It is plausible that analogous changes (coupled to those of the photocycle) are also operative in the cases of native bR and visual pigments. These conclusions are discussed in view of the light-induced conformational changes recently detected in native and artificial bR with an atomic force sensor.     

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.     

Picosecond resonance Raman evidence for unrelaxed heme in the (carbonmonoxy)myoglobin photoproduct

An actively and passively mode-locked Nd:YAG laser, producing 30-ps pulses of 1-mJ energy at 532 nm, has been used to photolyze (carbonmonoxy)myoglobin (MbCO) and generate its resonance Raman spectrum, which was recorded with a vidicon multichannel analyzer. The photoproduct spectrum was obtained by subtraction of the MbCO spectrum, obtained at lower incident power levels. Comparison with the spectrum of deoxyMb, obtained with the same apparatus, revealed frequency downshifts of approximately 4 cm-1, for bands at 1604, 1554, and 1542 cm-1, which are identified with porphyrin skeletal modes v10, v19, and v11. These frequencies are known to correlate inversely with the core size of the porphyrin ring, and the shifts imply a larger core size for the photoproduct than for deoxyMb. Similar shifts have been observed for the (carbonmonoxy)hemoglobin (HbCO) photoproduct; in that case, the shifts persist for longer than 20 ns, whereas they are absent in the MbCO photoproduct spectrum within 7 ns of photolysis. The unrelaxed state of the heme group region is therefore suggested to be maintained by protein forces, which relax more rapidly for Mb than Hb. This may reflect a tighter coupling in Hb of the out-of-plane movement of the Fe atom with the proximal histidine-containing F helix.     

Protein conformational changes in the bacteriorhodopsin photocycle

We report a comprehensive electron crystallographic analysis of conformational changes in the photocycle of wild-type bacteriorhodopsin and in a variety of mutant proteins with kinetic defects in the photocycle. Specific intermediates that accumulate in the late stages of the photocycle of wild-type bacteriorhodopsin, the single mutants D38R, D96N, D96G, T46V, L93A and F219L, and the triple mutant D96G/F171C/F219L were trapped by freezing two-dimensional crystals in liquid ethane at varying times after illumination with a light flash. Electron diffraction patterns recorded from these crystals were used to construct projection difference Fourier maps at 3.5 A resolution to define light-driven changes in protein conformation.Our experiments demonstrate that in wild-type bacteriorhodopsin, a large protein conformational change occurs within approximately 1 ms after illumination. Analysis of structural changes in wild-type and mutant bacteriorhodopsins under conditions when either the M or the N intermediate is preferentially accumulated reveals that there are only small differences in structure between M and N intermediates trapped in the same protein. However, a considerably larger variation is observed when the same optical intermediate is trapped in different mutants. In some of the mutants, a partial conformational change is present even prior to illumination, with additional changes occurring upon illumination. Selected mutations, such as those in the D96G/F171C/F219L triple mutant, can sufficiently destabilize the wild-type structure to generate almost the full extent of the conformational change in the dark, with minimal additional light-induced changes. We conclude that the differences in structural changes observed in mutants that display long-lived M, N or O intermediates are best described as variations of one fundamental type of conformational change, rather than representing structural changes that are unique to the optical intermediate that is accumulated. Our observations thus support a simplified view of the photocycle of wild-type bacteriorhodopsin in which the structures of the initial state and the early intermediates (K, L and M1) are well approximated by one protein conformation, while the structures of the later intermediates (M2, N and O) are well approximated by the other protein conformation. We propose that in wild-type bacteriorhodopsin and in most mutants, this conformational change between the M1 and M2 states is likely to make an important contribution towards efficiently switching proton accessibility of the Schiff base from the extracellular side to the cytoplasmic side of the membrane.