Viscosity dependence of O2 escape from respiratory proteins as a function of cosolvent molecular weight

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75 ± 0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Effect of viscogens on the kinetic response of a photoperturbed allosteric protein

By covalently binding a photoswitchable linker across the binding groove of the PDZ2 domain, a small conformational change can be photo-initiated that mimics the allosteric transition of the protein. The response of its binding groove is investigated with the help of ultrafast pump-probe IR spectroscopy from picoseconds to tens of microseconds. The temperature dependence of that response is compatible with diffusive dynamics on a rugged energy landscape without any prominent energy barrier. Furthermore, the dependence of the kinetics on the concentration of certain viscogens, sucrose, and glycerol, has been investigated. A pronounced viscosity dependence is observed that can be best fit by a power law, i.e., a fractional viscosity dependence. The change of kinetics when comparing sucrose with glycerol as viscogen, however, provides strong evidence that direct interactions of the viscogen molecule with the protein do play a role as well. This conclusion is supported by accompanying molecular dynamics simulations.

Viscosity Dependence of Some Protein and Enzyme Reaction Rates: Seventy-Five Years after Kramers

Kramers rate theory is a milestone in chemical reaction research, but concerns regarding the basic understanding of condensed phase reaction rates of large molecules in viscous milieu persist. Experimental studies of Kramers theory rely on scaling reaction rates with inverse solvent viscosity, which is often equated with the bulk friction coefficient based on simple hydrodynamic relations. Apart from the difficulty of abstraction of the prefactor details from experimental data, it is not clear why the linearity of rate versus inverse viscosity, k ∝ η(-1), deviates widely for many reactions studied. In most cases, the deviation simulates a power law k ∝ η(-n), where the exponent n assumes fractional values. In rate-viscosity studies presented here, results for two reactions, unfolding of cytochrome c and cysteine protease activity of human ribosomal protein S4, show an exceedingly overdamped rate over a wide viscosity range, registering n values up to 2.4. Although the origin of this extraordinary reaction friction is not known at present, the results indicate that the viscosity exponent need not be bound by the 0-1 limit as generally suggested. For the third reaction studied here, thermal dissociation of CO from nativelike cytochrome c, the rate-viscosity behavior can be explained using Grote-Hynes theory of time-dependent friction in conjunction with correlated motions intrinsic to the protein. Analysis of the glycerol viscosity-dependent rate for the CO dissociation reaction in the presence of urea as the second variable shows that the protein stabilizing effect of subdenaturing amounts of urea is not affected by the bulk viscosity. It appears that a myriad of factors as diverse as parameter uncertainty due to the difficulty of knowing the exact reaction friction and both mode and consequences of protein-solvent interaction work in a complex manner to convey as though Kramers rate equation is not absolute.

Viscosity dependence of passage through a fluctuating bottleneck

We generalize the model of a rate process involving the passage of an object through a fluctuating bottleneck. The rate of passage is considered to be proportional to a power function of the radius of the bottleneck with exponent α > 0. The fluctuations of the bottleneck are coupled to the motion of the surrounding medium and governed by Langevin dynamics. We show numerically and also explain analytically that for slow bottleneck fluctuations the long time decay rate of the process has a fractional power law dependence on the solvent viscosity with exponent α/(α + 2). The results are consistent with the experimental data on ligand binding to myoglobin, and might also be relevant to other reactions for which exponents between 0 and 1 were reported.

Internal friction in enzyme reactions

The empirical concept of internal friction was introduced 20 years ago. This review summarizes the results of experimental and theoretical studies that help to uncover the nature of internal friction. After the history of the concept, we describe the experimental challenges in measuring and interpreting internal friction based on the viscosity dependence of enzyme reactions. We also present speculations about the structural background of this viscosity dependence. Finally, some models about the relationship between the energy landscape and internal friction are outlined. Alternative concepts regarding the viscosity dependence of enzyme reactions are also discussed.

Temperature derivative spectroscopy to monitor the autoxidation decay of cytochromes P450

Temperature derivative spectroscopy (TDS), a type of relaxation spectroscopy, is a powerful tool to study protein dynamics (Berendzen, J.; Braunstein, D. Proc. Natl. Acad. Sci. U. S. A. 1990, 87, 1). We developed the version of temperature derivative spectroscopy to monitor kinetics of autoxidation of cytochromes P450 and applied it to study the properties of the oxy-ferrous complex of a human membrane bound P450, CYP19A1 (aromatase), and that of a bacterial soluble P450, CYP101 when bound with their most common substrates, androstenedione (AD) and camphor, respectively. TDS extends the panel of methods that can be used to monitor heme protein kinetics, providing a rapid measurement technique and enabling measurement of the autoxidation rate over a wide range of temperatures, yielding the activation energy as well as absolute reaction rate in a single experiment.

Viscosity scaling for the glassy phase of protein folding

Although commendable progress has been made in the understanding of the physics of protein folding, a key unresolved issue is whether Kramers' diffusion model of chemical reactions is generally applicable to activated barrier crossing events during folding. To examine the solvent viscosity effect on the folding transition of native-like trapped intermediates, laser flash photolysis has been used to measure the microsecond folding kinetics of a natively folded state of CO-liganded ferrocytochrome c (M-state) in the 1-250 cP range of glycerol viscosity at pH 7.0, 20 degrees C. The single rate coefficient for the folding of the M-state to the native state of the protein (i.e., the M --> N folding process) decreases initially when the solvent viscosity is low (<10 cP), but saturates at higher viscosity, indicating that Kramers model is not general enough for scaling the viscosity dependence of post-transition folding involving glassy dynamics. Analysis based on the Grote-Hynes idea of time dependent friction in conjunction with defect diffusion dynamics can account for the observed non-Kramers scaling.

Quaternary relaxations in sol-gel encapsulated hemoglobin studied via NIR and UV spectroscopy

In this work, we study the kinetics of the R --> T transition in hemoglobin using a combination of near-infrared and near-ultraviolet spectroscopy. We use a sol-gel encapsulation protocol to decelerate the conformational transitions and to avoid spectral perturbations arising from ligand migration and recombination. We monitor two spectroscopic markers: band III in the near-IR, which is a fine probe of the heme pocket conformation, and the tryptophan band in the near-UV, which probes the formation of the Trpbeta37-Aspalpha94 hydrogen bond, characteristic of the T structure, at the critical alpha1beta2 subunit interface. The time evolution of these two bands is monitored after deoxygenation of encapsulated oxyhemoglobin, obtained by diffusion of a reducing agent into the porous silica matrix. Characteristic spectral shifts are observed: comparison with myoglobin enables us to assign them to quaternary structure relaxations. Band III spectral relaxation is clearly nonexponential, and analysis with the Maximum Entropy Method enables us to identify three processes. On the other hand, near-UV spectral relaxation follows an exponential decay with a time constant closely corresponding to the second process observed in the near IR. Very interestingly, the rates of all processes markedly depend on the viscosity of the co-encapsulated solvent, following a power law. Our results reveal correlations between heme pocket relaxations, induced by the R --> T transition, and structural event(s) occurring at the alpha1beta2 interface and highlight their solvent dependence. The power law viscosity dependence of relaxation rates suggests that the observed protein relaxations are "slaved" to the co-encapsulated solvent. The stepwise character of the quaternary transition is also evidenced.

Protein folding is slaved to solvent motions

Proteins, the workhorses of living systems, are constructed from chains of amino acids, which are synthesized in the cell based on the instructions of the genetic code and then folded into working proteins. The time for folding varies from microseconds to hours. What controls the folding rate is hotly debated. We postulate here that folding has the same temperature dependence as the alpha-fluctuations in the bulk solvent but is much slower. We call this behavior slaving. Slaving has been observed in folded proteins: Large-scale protein motions follow the solvent fluctuations with rate coefficient k(alpha) but can be slower by a large factor. Slowing occurs because large-scale motions proceed in many small steps, each determined by k(alpha). If conformational motions of folded proteins are slaved, so a fortiori must be the motions during folding. The unfolded protein makes a Brownian walk in the conformational space to the folded structure, with each step controlled by k(alpha). Because the number of conformational substates in the unfolded protein is extremely large, the folding rate coefficient, k(f), is much smaller than k(alpha). The slaving model implies that the activation enthalpy of folding is dominated by the solvent, whereas the number of steps n(f) = k(alpha)/k(f) is controlled by the number of accessible substates in the unfolded protein and the solvent. Proteins, however, undergo not only alpha- but also beta-fluctuations. These additional fluctuations are local protein motions that are essentially independent of the bulk solvent fluctuations and may be relevant at late stages of folding.

Numerical Simulation of the Effect of Solvent Viscosity on the Motions of a β-Peptide Heptamer

This report examines the effect of a decrease in solvent viscosity on the simulated folding behaviour of a beta-peptide heptamer in methanol. Simulations of the molecular dynamics of the heptamer H-beta3-HVal-beta3-HAla-beta3-HLeu-(S,S)-beta3-HAla(alphaMe)-beta3-HVal-beta3-HAla-beta3-HLeu-OH in methanol, with an explicit representation of the methanol molecules, were performed for 80 ns at various solvent viscosities. The simulations indicate that at a solvent viscosity of one third of that of methanol, only the dynamic aspects of the folding process are altered, and that the rate of folding is increased. At a viscosity of one tenth of that of methanol, insufficient statistics are obtained within the 80 ns period. We suggest that 80 ns is an insufficient time to reach conformational equilibrium at very low viscosity because the dependence of the folding rate of a beta-peptide on solvent viscosity has two regimes; a result that was observed in another computational study for alpha-peptides.

Dominant features of protein reaction dynamics: conformational relaxation and ligand migration

Here, we review the dominant aspects of protein dynamics as revealed by studying hemoproteins using the combination of laser flash photolysis, kinetic spectroscopy and low temperature. The first breakthrough was the finding that geminate ligand rebinding with myoglobin is highly non-exponential at temperature T<200 K, providing evidence for the trapping of a large number of protein statistical substates. Another major advance was the introduction of a "model free" approach to analyze polychromatic kinetics in terms of their rate spectrum rather than to fit the data to some arbitrarily predefined kinetic scheme. Kinetic processes are identified and quantified directly from the rate spectrum without a priori assumptions. In recent years, further progresses were achieved by using xenon gas as a soft external perturbing agent that competes with ligand rebinding pathways by occupying hydrophobic protein cavities. The first part of this paper introduces several basic principles that are spread throughout a vast literature. The second part describes the main conclusions regarding conformational relaxation and ligand migration in hemoproteins obtained by combining these approaches.

Fast dynamics and stabilization of proteins: binary glasses of trehalose and glycerol

We present elastic and inelastic incoherent neutron scattering data from a series of trehalose glasses diluted with glycerol. A strong correlation with recently published protein stability data in the same series of glasses illustrates that the dynamics at Q >or= 0.71 A(-1) and omega > 200 MHz are important to stabilization of horseradish peroxidase and yeast alcohol dehydrogenase in these glasses. To the best of our knowledge, this is the first direct evidence that enzyme stability in a room temperature glass depends upon suppressing these short-length scale, high-frequency dynamics within the glass. We briefly discuss the coupling of protein motions to the local dynamics of the glass. Also, we show that T(g) alone is not a good indicator for the protein stability in this series of glasses; the glass that confers the maximum room-temperature stability does not have the highest T(g).

Internal friction controls the speed of protein folding from a compact configuration

Several studies have found millisecond protein folding reactions to be controlled by the viscosity of the solvent: Reducing the viscosity allows folding to accelerate. In the limit of very low solvent viscosity, however, one expects a different behavior. Internal interactions, occurring within the solvent-excluded interior of a compact molecule, should impose a solvent-independent upper limit to folding speed once the bulk diffusional motions become sufficiently rapid. Why has this not been observed? We have studied the effect of solvent viscosity on the folding of cytochrome c from a highly compact, late-stage intermediate configuration. Although the folding rate accelerates as the viscosity declines, it tends toward a finite limiting value approximately 10(5) s(-1) as the viscosity tends toward zero. This limiting rate is independent of the cosolutes used to adjust solvent friction. Therefore, interactions within the interior of a compact denatured polypeptide can limit the folding rate, but the limiting time scale is very fast. It is only observable when the solvent-controlled stages of folding are exceedingly rapid or else absent. Interestingly, we find a very strong temperature dependence in these "internal friction"-controlled dynamics, indicating a large energy scale for the interactions that govern reconfiguration within compact, near-native states of a protein.

Competition with xenon elicits ligand migration and escape pathways in myoglobin

Evidence for ligand migration toward the xenon-binding cavities in myoglobin comes from a number of laser photolysis studies of MbO2 including mutants and from cryo- and time-resolved crystallography of MbCO. To explore ligand migration in greater detail, we investigated the rebinding kinetics of both MbO2 and MbCO under a xenon partial pressure ranging from 1 to 16 atm over the temperature range (293-77 K). Below 180 K xenon affects to a significant, but minor, extent the thermodynamic parameters for rebinding from the primary docking site in each Mb taxonomic substate. Above 200 K the ligand migrates to the proximal Xe1 site but when the latter is occupied by xenon a new kinetic process appears. It is attributed to rebinding from transient docking sites located on the path between the primary and the secondary docking site of both ligands. Ligand escape exhibits a more complicated pattern than expected. At room temperature O2 and CO escape appears to take place exclusively from the primary site. In contrast, at T approximately 250 K, roughly 50% of the CO molecules that have escaped from the protein originate from the Xe1 secondary site.

Effect of solvent diffusion on the apomyoglobin-water interface

Few techniques can identify interactions between proteins and individual water molecules when the protein is in solution. The present work has sought to bridge the gap between the molecular level studies and the search for a physical property of the solution (bathing the proteins) that would regulate the protein hydration level. The properties of the solution were varied by adding nondenaturing solutes and solvents to the protein solutions and then studying their effect on the intrinsic fluorescence of apomyoglobin. The resolution of the tryptophan emission into the two component spectra corresponding to tryptophans W7 (accessible to the solvent) and W14 (buried in the protein matrix) has allowed us to probe two specific parts of the protein. Whereas W14 is not affected when the medium is altered, the analysis of W7 fluorescence has shown that cosolvent diffusion plays a dominant role in the mobility of water molecules near the protein surface.

Hydration, slaving and protein function

Protein dynamics is crucial for protein function. Proteins in living systems are not isolated, but operate in networks and in a carefully regulated environment. Understanding the external control of protein dynamics is consequently important. Hydration and solvent viscosity are among the salient properties of the environment. Dehydrated proteins and proteins in a rigid environment do not function properly. It is consequently important to understand the effect of hydration and solvent viscosity in detail. We discuss experiments that separate the two effects. These experiments have predominantly been performed with wild-type horse and sperm whale myoglobin, using the binding of carbon monoxide over a broad range of temperatures as a tool. The experiments demonstrate that data taken only in the physiological temperature range are not sufficient to understand the effect of hydration and solvent on protein relaxation and function. While the actual data come from myoglobin, it is expected that the results apply to most or all globular proteins.

The effect of water on the rate of conformational change in protein allostery

The influence of solvation on the rate of quaternary structural change is investigated in human hemoglobin, an allosteric protein in which reduced water activity destabilizes the R state relative to T. Nanosecond absorption spectroscopy of the heme Soret band was used to monitor protein relaxation after photodissociation of aqueous HbCO complex under osmotic stress induced by the nonbinding cosolute poly(ethylene glycol) (PEG). Photolysis data were analyzed globally for six exponential time constants and amplitudes as a function of osmotic stress and viscosity. Increases in time constants associated with geminate rebinding, tertiary relaxation, and quaternary relaxation were observed in the presence of PEG, along with a decrease in the fraction of hemes rebinding CO with the slow rate constant characteristic of the T state. An analysis of these results along with those obtained by others for small cosolutes showed that both osmotic stress and solvent viscosity are important determinants of the microscopic R --> T rate constant. The size and direction of the osmotic stress effect suggests that at least nine additional water molecules are required to solvate the allosteric transition state relative to the R-state hydration, implying that the transition state has a greater solvent-exposed area than either end state.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

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

The O(2) binding pocket of myohemerythrin: role of a conserved leucine

A conserved O(2) binding pocket residue in Phascolopsis gouldii myohemerythrin (myoHr), namely, L104, was mutated to several other residues, and the effects on O(2) association and dissociation rates, O(2) affinity, and autoxidation were examined. The L104V, -F, and -Y myoHrs formed stable O(2) adducts whose UV-vis and resonance Raman spectra closely matched those of wild-type oxymyoHr. The L104V mutation produced only minimal effects on either O(2) association or dissociation, whereas the L104F and -Y mutations resulted in 100-300-fold decreases in both O(2) association and dissociation rates. These decreases are attributed to introduction of steric restrictions into the O(2) binding pocket, which are not present in either wild-type or L104V myoHrs. The failure to observe increased O(2) association or dissociation rates for L104V indicates that the side chain of leucine at position 104 does not sterically "gate" O(2) entry into or exit from the binding pocket in the rate-determining step(s). L104V myoHr autoxidized approximately 3 times faster than did wild type, whereas L104T autoxidized >10(6) times faster than did wild type. The latter large increase is attributed to increased side chain polarity, thereby increasing water occupancy in the oxymyoHr binding pocket. These results indicate that L104 contributes a hydrophobic barrier that restricts water entry into the oxymyoHr binding pocket. Thus, a leucine at position 104 in myoHr appears to have the optimal combination of size and hydrophobicity to facilitate O(2) binding while simultaneously inhibiting autoxidation.

Microscopic viscosity and rotational diffusion of proteins in a macromolecular environment

The Stokes-Einstein-Debye equation is currently used to obtain information on protein size or on local viscosity from the measurement of the rotational correlation time. However, the implicit assumptions of a continuous and homogeneous solvent do not hold either in vivo, because of the high density of macromolecules, or in vitro, where viscosity is adjusted by adding viscous cosolvents of various size. To quantify the consequence of nonhomogeneity, we have measured the rotational Brownian motion of three globular proteins with molecular mass from 66 to 4000 kD in presence of 1.5 to 2000 kD dextrans as viscous cosolvents. Our results indicate that the linear viscosity dependence of the Stokes-Einstein relation must be replaced by a power law to describe the rotational Brownian motion of proteins in a macromolecular environment. The exponent of the power law expresses the fact that the protein experiences only a fraction of the hydrodynamic interactions of macromolecular cosolvents. An explicit expression of the exponent in terms of protein size and cosolvent's mass is obtained, permitting definition of a microscopic viscosity. Experimental data suggest that a similar effective microviscosity should be introduced in Kramers' equation describing protein reaction rates.

Flavin fluorescence dynamics and photoinduced electron transfer in Escherichia coli glutathione reductase

Time-resolved polarized flavin fluorescence was used to study the active site dynamics of Escherichia coli glutathione reductase (GR). Special consideration was given to the role of Tyr177, which blocks the access to the NADPH binding-site in the crystal structure of the enzyme. By comparing wild-type GR with the mutant enzymes Y177F and Y177G, a fluorescence lifetime of 7 ps that accounts for approximately 90% of the fluorescence decay could be attributed to quenching by Y177. Based on the temperature invariance for this lifetime, and the very high quenching rate, electron transfer from Y177 to the light-excited isoalloxazine part of flavin adenine dinucleotide (FAD) is proposed as the mechanism of flavin fluorescence quenching. Contrary to the mutant enzymes, wild-type GR shows a rapid fluorescence depolarization. This depolarization process is likely to originate from a transient charge transfer interaction between Y177 and the light-excited FAD, and not from internal mobility of the flavin, as has previously been proposed. Based on the fluorescence lifetime distributions, the mutants Y177F and Y177G have a more flexible protein structure than wild-type GR: in the range of 223 K to 277 K in 80% glycerol, both tyrosine mutants mimic the closely related enzyme dihydrolipoyl dehydrogenase. The fluorescence intensity decays of the GR enzymes can only be explained by the existence of multiple quenching sites in the protein. Although structural fluctuations are likely to contribute to the nonexponential decay and the probability of quenching by a specific site, the concept of conformational substates need not be invoked to explain the heterogeneous fluorescence dynamics.

Reduction of protein volume and compressibility by macromolecular cosolvents: dependence on the cosolvent molecular weight

The partial specific volume (V) and adiabatic compressibility (beta) of myoglobin have been shown to be reduced by small cosolvents such as glycerol (A. Priev, A. Almagor, S. Yedgar, B. Gavish, Biochemistry 35 (1996) 2061-2066). To elucidate the effect of the cosolvent size on these protein properties, in the present study we determined V and beta of myoglobin in solutions containing a homologous cosolvent series from sucrose to dextran--500 (M.W. 500,000). It was found that in addition to the expected effect of the cosolvent concentration, V and beta decrease with increasing cosolvent M.W. This suggests that structural properties of the cosolvent contribute to its effect on the protein interior.

Solvent composition and viscosity effects on the kinetics of CO binding to horse myoglobin

Ligand binding to myoglobin in aqueous solution involves two kinetic components, one extramolecular and one intramolecular, which have been interpreted in terms of two sequential kinetic barriers. In mixed solvents and sub-zero temperatures, the outer barrier increases and the inner barrier splits into several components, giving rise to fast intramolecular recombination. The nature of these barriers and their relation to structural relaxation are examined using the effect of solvent composition and viscosity on the kinetics of CO binding to horse myoglobin in 60% ethylene glycol/water, 75% and 90% glycerol/water, 80% and 92% sucrose/water solutions. Measurements of the corresponding solvent structural relaxation rates by frequency resolved calorimetry allow us to discriminate between solvent composition and viscosity-related effects. The outer kinetic barrier controlling ligand entry and release depends on the viscosity consistent with Kramers-Stokes law of activated escape in the presence of friction. At high cosolvent concentration, we observe deviations from Stokes law, implying a smaller microviscosity at the protein-solvent interface as compared to the bulk. The inner barrier and its coupling to structural relaxation appears to be independent of viscosity but changes with solvent composition. As a possible explanation, we discuss the role of distal water molecules in the formation of the effective inner barrier. At low temperatures, this barrier has a distributed height, depending only slightly on the nature of the cosolvent and temperature at low cosolvent concentrations. In contrast, myoglobin embedded in a sucrose glass (92% sucrose/water) exhibits a temperature-dependent and bimodal enthalpy distribution. This result demonstrates that the exchange between protonation states of His64, A0 left and right arrow A1, can take place in the glass and at temperatures as low as 80 K.

Structures of wild-type chloromet and L103N hydroxomet Themiste zostericola myohemerythrins at 1.8 A resolution

Myohemerythrin (Mhr) is a nonheme iron oxygen carrier found in the retractor muscles of marine "peanut" worms. The X-ray crystal structures of two recombinant Themiste zostericola Mhrs are reported to a resolution of 1.8 A. Surprisingly, the met wild-type structure (R = 17.8%) was found to contain chloride bound to Fe2, while coordinated hydroxide was found in the met L103N structure (R = 18.3%). An internal water molecule was also found distal to the Fe-O-Fe center of the mutant protein, forming hydrogen bonds with the coordinated hydroxide and the OD1 atom of Asn-103. This finding is consistent with the kinetic and spectroscopic results reported for the L103N mutant Mhr [Raner, G. M., Martins, L. J., & Ellis, W. R., Jr. (1997) Biochemistry 36, 7037-7043]. Possible roles for the side chain of residue 103 (Leu in wild-type Mhr) in gating ligand binding are also discussed.

Cell membrane fluctuations are regulated by medium macroviscosity: evidence for a metabolic driving force

Extracellular fluid macroviscosity (EFM), modified by macromolecular cosolvents as occurs in body fluids, has been shown to affect cell membrane protein activities but not isolated proteins. In search for the mechanism of this phenomenon, we examined the effect of EFM on mechanical fluctuations of the cell membrane of human erythrocytes. The macroviscosity of the external medium was varied by adding to it various macromolecules [dextrans (70, 500, and 2,000 kDa), polyethylene glycol (20 kDa), and carboxymethyl-cellulose (100 kDa)], which differ in size, chemical nature, and in their capacity to increase fluid viscosity. The parameters of cell membrane fluctuations (maximal amplitude and half-width of amplitude distribution) were diminished with the elevation of solvent macroviscosity, regardless of the cosolvent used to increase EFM. Because thermally driven membrane fluctuations cannot be damped by elevation of EFM, the existence of a metabolic driving force is suggested. This is supported by the finding that in ATP-depleted red blood cells elevation of EMF did not affect cell membrane fluctuations. This study demonstrates that (i) EFM is a regulator of membrane dynamics, providing a possible mechanism by which EFM affects cell membrane activities; and (ii) cell membrane fluctuations are driven by a metabolic driving force in addition to the thermal one.