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.     

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.     

Conformational substates and motions in myoglobin. External influences on structure and dynamics.

Myoglobin, a simppe dioxygen-storage protein, is a good laboratory for the investigation of the connection between protein structure, dynamics, and function. Fourier-transform infrared spectroscopy on carbon-monoxymyoglobin (MbCO) shows three major CO bands. These bands are excellent probes for the investigation of the structure-function relationship. They have different CO binding kinetics and their CO dipoles form different angles with respect to the heme normal, implying that MbCO exists in three major conformational substates, A0, A1, and A3. The entropies and enthalpies of these substates depend on temperature above approximately 180 K and are influenced by pH, solvent, and pressure. These results suggest that even a protein as simple as Mb can assume a small number of clearly different structures that perform the same function, but with different rates. Moreover, protein structure and dynamics depend strongly on the interaction of the protein with its environment.     

Kinetic evidence for three photolyzable taxonomic conformational substates in oxymyoglobin

The kinetics of oxygen geminate binding with the taxonomic substates of MbO2 are reported. The maximum entropy method was used to analyze the rebinding kinetics of MbCO and MbO2 monitored in the Soret. The resulting rate distributions were found to consist of a small number of overlapping bands. A global parametric fit of a series of rate distributions recorded at several temperatures was performed using a Gaussian basis set to resolve the individual enthalpy distributions P(H). This approach was first validated by showing that the well-documented taxonomic substates of MbCO could be recovered. The method was then applied to MbO2. Three taxonomic substates were identified at pH 4.8, whereas only two of them contribute to oxygen geminate rebinding at pH 7.0. These findings show that, similarly to MbCO, MbO2 also exists as three photolyzable and kinetically different taxonomic substates and suggest reconsidering the issue of the photolysis quantum yield of MbO2.     

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

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

Rebinding and relaxation in the myoglobin pocket

The infrared stretching bands of carboxymyoglobin (MbCO) and the rebinding of CO to Mb after photodissociation have been studied in the temperature range 10-300 K in a variety of solvents. Four stretching bands imply that MbCO can exist in four substates, A0-A3. The temperature dependences of the intensities of the four bands yield the relative binding enthalpies and and entropies. The integrated absorbances and pH dependences of the bands permit identification of the substates with the conformations observed in the X-ray data (Kuriyan et al., J. Mol. Biol. 192 (1986) 133). At low pH, A0 is hydrogen-bonded to His E7. The substates A0-A3 interconvert above about 180 K in a 75% glycerol/water solvent and above 270 K in buffered water. No major interconversion is seen at any temperature if MbCO is embedded in a solid polyvinyl alcohol matrix. The dependence of the transition on solvent characteristics is explained as a slaved glass transition. After photodissociation at low temperature the CO is in the heme pocket B. The resulting CO stretching bands which are identified as B substates are blue-shifted from those of the A substates. At 40 K, rebinding after flash photolysis has been studied in the Soret, the near-infrared, and the integrated A and B substates. All data lie on the same rebinding curve and demonstrate that rebinding is nonexponential in time from at least 100 ns to 100 ks. No evidence for discrete exponentials is found. Flash photolysis with monitoring in the infrared region shows four different pathways within the pocket B to the bound substates Ai. Rebinding in each of the four pathways B----A is nonexponential in time to at least 10 ks and the four pathways have different kinetics below 180 K. From the time and temperature dependence of the rebinding, activation enthalpy distributions g(HBA) and preexponentials ABA are extracted. No pumping from one A substate to another, or one B substate to another, is observed below the transition temperature of about 180 K. If MbCO is exposed to intense white light for 10-10(3) s before being fully photolyzed by a laser flash, the amplitude of the long-lived states increases. The effect is explained in terms of a hierarchy of substates and substate symmetry breaking. The characteristics of the CO stretching bands and of the rebinding processes in the heme pocket depend strongly on the external parameters of solvent, pH and pressure. This sensitivity suggests possible control mechanisms for protein reactions.     

Conformational transitions in protein-protein association: binding of fasciculin-2 to acetylcholinesterase

The neurotoxin fasciculin-2 (FAS2) is a picomolar inhibitor of synaptic acetylcholinesterase (AChE). The dynamics of binding between FAS2 and AChE is influenced by conformational fluctuations both before and after protein encounter. Submicrosecond molecular dynamics trajectories of apo forms of fasciculin, corresponding to different conformational substates, are reported here with reference to the conformational changes of loop I of this three-fingered toxin. This highly flexible loop exhibits an ensemble of conformations within each substate corresponding to its functions. The high energy barrier found between the two major substates leads to transitions that are slow on the timescale of the diffusional encounter of noninteracting FAS2 and AChE. The more stable of the two apo substates may not be the one observed in the complex with AChE. It seems likely that the more stable apo form binds rapidly to AChE and conformational readjustments then occur in the resulting encounter complex.     

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.     

Low temperature optical absorption spectroscopy: an approach to the study of stereodynamic properties of hemeproteins

In this short review we show how suitable analysis of the temperature dependence of the optical absorption spectra of metalloproteins can give insight into their stereodynamic properties in the region of the chromophore. To this end, the theory of coupling between an intense allowed electronic transition of a chromophore and Franck-Condon active vibrations of the nearby atoms is applied to the Soret band of hemeproteins to obtain an analytical expression suitable for fitting the spectral profile at various temperatures. The reported approach enables one to separate the various contributions to the overall bandwidth together with the parameters that characterize the vibrational coupling. The thermal behavior of these quantities gives information on the dynamic properties of the active site and on their dependence upon protein structure and ligation state. The Soret band of hemeproteins appears to be coupled to high frequency vibrational modes of the heme group (as already shown by resonance Raman spectroscopy) and to a bath of low frequency modes most likely deriving from the bulk of the protein. For the deoxy derivatives inhomogeneous broadening arising from conformational heterogeneity appears to contribute substantially to the linewidth. The data indicate the onset; at temperatures near 180 K, of large scale anharmonic motions that can be attributed to jumping among different conformational substates of the protein.Abbreviations MbCO Carbonmonoxy-myoglobin - Mb Deoxymyoglobin - Mb³⁺ Aquomet-myoglobin - SWMbCO Spermwhale carbonmonoxy-myoglobin - SWMb Spermwhale deoxy-myoglobin Correspondence to: A. Cupane     

Starting‐structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD‐derived order parameters

Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin‐binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295‐4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR‐derived order parameters. Thus, MD‐derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Statistical properties of the dichotomous noise generated in biochemical processes

Dichotomous noise detected with the help of various single-molecule techniques convincingly reveals the actual occurrence of a multitude of conformational substates composing the native state of proteins. The nature of the stochastic dynamics of transitions between these substates is determined by the particular statistical properties of the noise observed. These involve nonexponential and possibly oscillatory time decay of the second order autocorrelation function, its relation to the third order autocorrelation function, and a relationship to dwell-time distribution densities and their correlations. Processes gated by specific conformational substates are distinguished from those with fluctuating barriers. This study throws light on the intriguing matter of the possibility of multiple stepping of the myosin motor along the actin filament per ATP molecule hydrolyzed. Paper authored by participants of the international conference: International Workshop on Ionic Channels, Szczyrk, Poland, May 27 – June 01, 2007. Publication cost was covered by the organisers of this meeting.     

The sodium channel activator Brevetoxin-3 uncovers a multiplicity of different open states of the cardiac sodium channel.

The interaction of Brevetoxin 3 (Pbtx-3), a sodium channel activator, with the cardiac sodium channel was studied at the single channel level. It was found that Pbtx-3 (20 microM) shifted steady-state activation to negative potentials, without major effects on the time course of macroscopic activation or macroscopic currents decay, as calculated from averaged single-channel records. Single-channel open times were found to be prolonged. Under the influence of the toxin, sodium channel openings could be observed frequently even at maintained depolarisation. These openings occurred to at least nine different subconductance levels of the open state with smaller conductivities than the maximal one and differed in their open times. Current amplitudes of these open substates were found to cluster around certain amplitude values. Appearance of substates at maintained depolarisation was dependent on the transmembrane potential (Em): Substates with smaller conductivity appeared more frequently at lower Em values whereas at higher Em values substates with higher conductivity values dominated. Furthermore, it was demonstrated that appearance of substates did not result from incomplete recovery from inactivation. From these observations it was concluded that the open substates observed correspond to different conformational states of the channel's activation gates. Under physiological conditions, when the sodium channel opens directly from its closed state these 'incomplete'-open states of the cardiac sodium channel are obscured by fast gating transitions between the corresponding, electrically silent, preopen states. Thus, Pbtx-3 acts mainly via stabilisation of the channel's preopen and different open states. A classification of sodium channel modifiers, based on their interaction with different conformational states of the channel is suggested.     

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.     

Reconciling Mediating and Slaving Roles of Water in Protein Conformational Dynamics

Proteins accomplish their physiological functions with remarkably organized dynamic transitions among a hierarchical network of conformational substates. Despite the essential contribution of water molecules in shaping functionally important protein dynamics, their exact role is still controversial. Water molecules were reported either as mediators that facilitate or as masters that slave protein dynamics. Since dynamic behaviour of a given protein is ultimately determined by the underlying energy landscape, we systematically analysed protein self energies and protein-water interaction energies obtained from extensive molecular dynamics simulation trajectories of barstar. We found that protein-water interaction energy plays the dominant role when compared with protein self energy, and these two energy terms on average have negative correlation that increases with increasingly longer time scales ranging from 10 femtoseconds to 100 nanoseconds. Water molecules effectively roughen potential energy surface of proteins in the majority part of observed conformational space and smooth in the remaining part. These findings support a scenario wherein water on average slave protein conformational dynamics but facilitate a fraction of transitions among different conformational substates, and reconcile the controversy on the facilitating and slaving roles of water molecules in protein conformational dynamics.     

Backbone dynamics and energetics of a calmodulin domain mutant exchanging between closed and open conformations.

Previous studies have suggested that the Ca2+-saturated E140Q mutant of the C-terminal domain of calmodulin exhibits equilibrium exchange between "open" and "closed" conformations similar to those of the Ca2+-free and Ca2+-saturated states of wild-type calmodulin. The backbone dynamics of this mutant were studied using15N spin relaxation experiments at three different temperatures. Measurements at each temperature of the15N rate constants for longitudinal and transverse auto-relaxation, longitudinal and transverse cross-correlation relaxation, and the1H-15N cross-relaxation afforded unequivocal identification of conformational exchange processes on microsecond to millisecond time-scales, and characterization of fast fluctuations on picosecond to nanosecond time-scales using model-free approaches. The results show that essentially all residues of the protein are involved in conformational exchange. Generalized order parameters of the fast internal motions indicate that the conformational substates are well folded, and exclude the possibility that the exchange involves a significant population of unfolded or disordered species. The temperature dependence of the order parameters offers qualitative estimates of the contribution to the heat capacity from fast fluctuations of the protein backbone, revealing significant variation between the well-ordered secondary structure elements and the more flexible regions. The temperature dependence of the conformational exchange contributions to the transverse auto-relaxation rate constants directly demonstrates that the microscopic exchange rate constants are greater than 2.7x10(3)s-1at 291 K. The conformational exchange contributions correlate with the chemical shift differences between the Ca2+-free and Ca2+-saturated states of the wild-type protein, thereby substantiating that the conformational substates are similar to the open and closed states of wild-type calmodulin. Taking the wild-type chemical shifts to represent the conformational substates of the mutant and populations estimated previously, the microscopic exchange rate constants could be estimated as 2x10(4)to 3x10(4)s-1at 291 K for a subset of residues. The temperature depen dence of the exchange allows the characterization of apparent energy barriers of the conformational transition, with results suggesting a complex process that does not correspond to a single global transition between substates.     

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).     

Solvent isotope effect on macromolecular dynamics in <Emphasis Type="Italic">E. coli</Emphasis>

Elastic incoherent neutron scattering was used to explore solvent isotope effects on average macromolecular dynamics in vivo. Measurements were performed on living E. coli bacteria containing H(2)O and D(2)O, respectively, close to physiological conditions of temperature. Global macromolecular flexibility, expressed as mean square fluctuation (MSF) values, and structural resilience in a free energy potential, expressed as a mean effective force constant, [Symbol: see text]k'[Symbol: see text], were extracted in the two solvent conditions. They referred to the average contribution of all macromolecules inside the cell, mostly dominated by the internal motions of the protein fraction. Flexibility and resilience were both found to be smaller in D(2)O than in H(2)O. A difference was expected because the driving forces behind macromolecular stabilization and dynamics are different in H(2)O and D(2)O. In D(2)O, the hydrophobic effect is known to be stronger than in H(2)O: it favours the burial of non-polar surfaces as well as their van der Waals' packing in the macromolecule cores. This may lead to the observed smaller MSF values. In contrast, in H(2)O, macromolecules would present more water-exposed surfaces, which would give rise to larger MSF values, in particular at the macromolecular surface. The smaller [Symbol: see text]k'[Symbol: see text] value suggested a larger entropy content in the D(2)O case due to increased sampling of macromolecular conformational substates.     

Conformational substates in enzyme mechanism: the 120 K structure of alpha-lytic protease at 1.5 A resolution.

Insight into the dynamic properties of alpha-lytic protease (alpha LP) has been obtained through the use of low-temperature X-ray crystallography and multiple-conformation refinement. Previous studies of alpha LP have shown that the residues around the active site are able to move significantly to accommodate substrates of different sizes. Here we show a link between the ability to accommodate ligands and the dynamics of the binding pocket. Although the structure of alpha LP at 120 K has B-factors with a uniformly low value of 4.8 A2 for the main chain, four regions stand out as having significantly higher B-factors. Because thermal motion should be suppressed at cryogenic temperatures, the high B-factors are interpreted as the result of trapped conformational substates. The active site residues that are perturbed during accommodation of different substrates are precisely those showing conformational substates, implying that substrate binding selects a subset of conformations from the ensemble of accessible states. To better characterize the precise nature of these substates, a protein model consisting of 16 structures has been refined and evaluated. The model reveals a number of features that could not be well-described by conventional B-factors: for example, 40% of the main-chain residue conformations are distributed asymmetrically or in discrete clusters. Furthermore, these data demonstrate an unexpected correlation between motions on either side of the binding pocket that we suggest is a consequence of "dynamic close packing." These results provide strong evidence for the role of protein dynamics in substrate binding and are consistent with the results of dynamic studies of ligand binding in myoglobin and ribonuclease A.     

Active site conformation in myoglobin as determined by X-ray absorption spectroscopy

X-ray absorption fine structure experiments were performed to study structural and dynamic aspects of the active site of various forms of myoglobin. The structures determined for deoxyMb, MbCO, and MbO2 are consistent with the structure established by X-ray absorption fine structure experiment and X-ray crystallography. The first shell of ferrous MbNO determined contains 5 nitrogens located at 2.02 A and a short NO bond length of 1.76 A. This study focuses on the change of the XAFS Debye-Waller factor with temperature, which is a measure of thermal and static disorder. It was found that the changes of Debye-Waller factor with temperature for the Mb proteins, except deoxyMb, are consistent with a simple Einstein model, in which a single frequency was assumed for the bond stretching modes. In contrast, the temperature dependence of deoxyMb cannot be fitted to the Einstein model and a large disorder was found at low temperatures, which indicates the existence of conformational substates of the active site.     

Ionic channels with conformational substates.

Recent studies of protein dynamics suggest that ionic channels can assume many conformational substates. Long-lived substates have been directly observed in single-channel current records. In many cases, however, the lifetimes of conformational states will be far below the theoretical limit of time resolution of single-channel experiments. The existence of such hidden substates may strongly influence the observable (time-averaged) properties of a channel, such as the concentration dependence of conductance. A channel exhibiting fast, voltage-dependent transitions between different conductance states may behave as an intrinsic rectifier. In the presence of more than one permeable ion species, coupling between ionic fluxes may occur, even when the channel has only a single ion-binding site. In special situations the rate of ion translocation becomes limited by the rate of conformational transitions, meaning that the channel approaches the kinetic behavior of a carrier. As a result of the strong coulombic interaction between an ion in a binding site and polar groups of the protein, rate constants of conformational transitions may depend on the occupancy of the binding site. Under this condition a nonequilibrium distribution of conformational states is created when ions are driven through the channel by an external force. This may lead to an apparent violation of microscopic reversibility, i.e., to a situation in which the frequency of transitions from state A to state B is no longer equal to the transition frequency from state B to state A.