Determination of microscopic rate constants for CO binding and migration in myoglobin encapsulated in silica gels

CO rebinding kinetics after nanosecond photolysis of myoglobin encapsulated in wet silica gels exhibits an enhanced geminate phase that allows the determination of the microscopic rate constants and the activation barriers for distinct ligand docking sites inside the protein matrix. Using a maximum entropy method, we demonstrate that the geminate phase can be well-described by a biphasic lifetime distribution, reflecting rebinding from the distal and proximal sites. Microscopic rates and activation barriers were estimated using a four-state model.     

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

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

Evidence for two geminate rebinding states following laser photolysis of R state hemoglobin encapsulated in wet silica gels

In this letter we report the first experimental evidence for CO rebinding to human hemoglobin from multiple geminate states. The analysis of the rebinding kinetics using a maximum entropy method allowed the identification of two distinct rebinding states within the protein matrix, which become populated under conditions of increased viscosity in a silica gel at high glycerol concentration. Our findings suggest the presence of at least two distinct docking sites for the photolyzed ligand. Assuming a minimal four-state model, we estimate the microscopic rates and the activation energies for the elementary processes.     

Free-energy barriers in MbCO rebinding

The rebinding of CO to myoglobin (Mb) from locations around the active site is studied using a combination of molecular dynamics and stochastic simulations for native and L29F mutant Mb. The interaction between the dissociated ligand and the protein environment is described by the recently developed fluctuating three-point charge model for the CO molecule. Umbrella sampling along trajectories, previously found to sample the binding site (B) and the Xe4 pocket, is used to construct free-energy profiles for the ligand escape. On the basis of the Smoluchowski equation, the relaxation of different initial population distributions is followed in space and time. For native Mb at room temperature, the calculated rebinding times are in good agreement with experimental values and give an inner barrier of 4.3 kcal/mol between the docking site B (Mb...CO) and the A state (bound MbCO), compared to an effective barrier, Heff, of 4.5 kcal/mol and barriers into the majority conformation A1 and the minority conformation A3 of 2.4 and 4.3 kcal/mol, respectively. In the case of the L29F mutant, the free-energy surface is flatter and the dynamics is much more rapid. As was found in experiment, escape to the Xe4 pocket is facile for L29F whereas, for native Mb, the barriers to this site are larger. At lower temperatures, the rebinding dynamics is delayed by orders of magnitude also due to increased barriers between the docking sites.     

Ligand migration in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

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

The reactivity with CO of AHb1 and AHb2 from Arabidopsis thaliana is controlled by the distal HisE7 and internal hydrophobic cavities

The nonsymbiotic hemoglobins, AHb1 and AHb2, have recently been isolated from Arabidopsis thaliana. Using steady-state and time-resolved spectroscopic methods, we show that Fe2+ AHb1 contains a mixture of penta- and hexacoordinated heme, while Fe2+ AHb2 is fully hexacoordinated. In the CO complexes, polar interactions and H-bonds with the ligand are stronger for AHb1 than for AHb2. The ligand binding kinetics are substantially different, reflecting the distribution between the penta- and hexacoordinated species, and indicate that protein dynamics and ligand migration pathways are very specific for each of the two proteins. In particular, a very small, non-exponential geminate rebinding observed in AHb1 suggests that the distal heme cavity is connected with the exterior by a relatively open channel. The large, temperature-dependent geminate rebinding observed for AHb2 implies a major role of protein dynamics in the ligand migration from the distal cavity to the solvent. The structures of AHb1 and AHb2, modeled on the basis of the homologous rice hemoglobin, exhibit a different cavity system that is fully compatible with the observed ligand binding kinetics. Overall, these kinetic and structural data are consistent with the putative NO-dioxygenase activity previously attributed to AHb1, whereas the role of AHb2 remains elusive.     

Temperature-dependent studies of NO recombination to heme and heme proteins

The rebinding kinetics of NO to the heme iron of myoglobin (Mb) is investigated as a function of temperature. Below 200 K, the transition-state enthalpy barrier associated with the fastest (approximately 10 ps) recombination phase is found to be zero and a slower geminate phase (approximately 200 ps) reveals a small enthalpic barrier (approximately 3 +/- 1 kJ/mol). Both of the kinetic rates slow slightly in the myoglobin (Mb) samples above 200 K, suggesting that a small amount of protein relaxation takes place above the solvent glass transition. When the temperature dependence of the NO recombination in Mb is studied under conditions where the distal pocket is mutated (e.g., V68W), the rebinding kinetics lack the slow phase. This is consistent with a mechanism where the slower (approximately 200 ps) kinetic phase involves transitions of the NO ligand into the distal heme pocket from a more distant site (e.g., in or near the Xe4 cavity). Comparison of the temperature-dependent NO rebinding kinetics of native Mb with that of the bare heme (PPIX) in glycerol reveals that the fast (enthalpically barrierless) NO rebinding process observed below 200 K is independent of the presence or absence of the proximal histidine ligand. In contrast, the slowing of the kinetic rates above 200 K in MbNO disappears in the absence of the protein. Generally, the data indicate that, in contrast to CO, the NO ligand binds to the heme iron through a "harpoon" mechanism where the heme iron out-of-plane conformation presents a negligible enthalpic barrier to NO rebinding. These observations strongly support a previous analysis (Srajer et al. J. Am. Chem. Soc. 1988, 110, 6656-6670) that primarily attributes the low-temperature stretched exponential rebinding of MbCO to a quenched distribution of heme geometries. A simple model, consistent with this prior analysis, is presented that explains a variety of MbNO rebinding experiments, including the dependence of the kinetic amplitudes on the pump photon energy.     

Dynamics of nitric oxide rebinding and escape in horseradish peroxidase

Ultrafast kinetic measurements of NO rebinding to horseradish peroxidase (HRP) are reported for the first time. The geminate kinetics are found to be exponential for all HRP samples studied. The ferric forms of HRP have NO geminate recombination time constants in the range of 15-30 ps, while the ferrous form has a time constant of approximately 7 ps. The simple exponential NO geminate kinetics found for HRP demonstrate that heme relaxation is not the underlying source of the nonexponential NO rebinding in myoglobin (Mb). The NO ligand escape rates from HRP are also determined, and they are found to depend dramatically on the presence or absence of the competitive inhibitor benzohydroxamic acid (BHA). The kinetic results indicate that, in contrast to Mb, there is direct solvent access to the distal heme pocket of HRP.     

The escape process of carbon monoxide from myoglobin to solution at physiological temperature

The carbon monoxide (CO) docking sites involved in the ligand escape process from the iron atom in hem of myoglobin (Mb) to solution at physiological temperature were studied on the basis of the effect of xenon (Xe) on the ligand escape rate by the transient grating (TG) technique. The TG method provides a direct measurement of the changes in molecular volume. The apparent CO escaping rate and the volume contraction increase with increasing Xe pressure. The pressure dependence of the rate is consistent with that of the Xe population at the Xe(1) site. This result clearly shows that CO is trapped at the Xe(1) site before escaping to solvent in a Xe-free solution at room temperature. It is shown that only CO but not the trapped Xe is released by the photoexcitation of the Xe-trapped MbCO. A dissociation scheme is proposed to explain the enhancement of the escaping rate by the presence of Xe(1). There are two branches for the CO escaping pathway. The dominant part of the dissociated CO escapes to the solvent through the Xe(1) trapping site under the Xe-free condition, and there are at least three intermediate states along this pathway. When a Xe atom blocks the Xe(1) site, the CO escapes through another route.     

Hydrophobic distal pocket affects NO-heme geminate recombination dynamics in dehaloperoxidase and H64V myoglobin

The recombination dynamics of NO with dehaloperoxidase (DHP) from Amphitrite ornata following photolysis were measured by femtosecond time-resolved absorption spectroscopy. Singular value decomposition (SVD) analysis reveals two important basis spectra. The first SVD basis spectrum reports on the population of photolyzed NO molecules and has the appearance of the equilibrium difference spectrum between the deoxy and NO forms of DHP. The first basis time course has two kinetic components with time constants of tau(11) approximately 9 ps and tau(12) approximately 50 ps that correspond to geminate recombination. The fast geminate process tau(11) arises from a contact pair with the heme iron in a bound state with S = 3/2 spin. The slow geminate process tau(12) corresponds to the recombination from a more remote docking site >3 A from the heme iron with the greater barrier corresponding to a S = 5/2 spin state. The second SVD basis spectrum represents a time-dependent Soret band shift indicative of heme photophysical processes and protein relaxation with time constants of tau(21) approximately 3 ps and tau(22) approximately 17 ps, respectively. A comparison between the more rapid rate constant of the slow geminate phase in DHP-NO and horse heart myoglobin (HHMbNO) or sperm whale myoglobin (SWMbNO) suggests that protein interactions with photolyzed NO are weaker in DHP than in the wild-type MbNOs, consistent with the hydrophobic distal pocket of DHP. The slower protein relaxation rate tau(22) in DHP-NO relative to HHMbNO implies less effective trapping in the docking site of the distal pocket and is consistent with a greater yield for the fast geminate process. The trends observed for DHP-NO also hold for the H64V mutant of SWMb (H64V MbNO), consistent with a more hydrophobic distal pocket for that protein as well. We examine the influence of solution viscosity on NO recombination by varying the glycerol content in the range from 0% to 90% (v/v). The dominant effect of increasing viscosity is the increase of the rate of the slow geminate process, tau(12), coupled with a population decrease of the slow geminate component. Both phenomena are similar to the effect of viscosity on wild-type Mb due to slowing of protein relaxation resulting from an increased solution viscosity and protein surface dehydration.     

Towards understanding non-equivalence of α and β subunits within human hemoglobin in conformational relaxation and molecular oxygen rebinding

Picosecond to millisecond laser time-resolved transient absorption spectroscopy was used to study molecular oxygen (O₂) rebinding and conformational relaxation following O₂ photodissociation in the α and β subunits within human hemoglobin in the quaternary R-like structure. Oxy-cyanomet valency hybrids, α₂(Fe²⁺–O₂)β₂(Fe³⁺–CN) and α₂(Fe³⁺–CN)β₂(Fe²⁺–O₂), were used as models for oxygenated R-state hemoglobin. An extended kinetic model for geminate O₂ rebinding in the ferrous hemoglobin subunits, ligand migration between the primary and secondary docking site(s), and nonexponential tertiary relaxation within the R quaternary structure, was introduced and discussed. Significant functional non-equivalence of the α and β subunits in both the geminate O₂ rebinding and concomitant structural relaxation was revealed. For the β subunits, the rate constant for the geminate O₂ rebinding to the unrelaxed tertiary structure and the tertiary transition rate were found to be greater than the corresponding values for the α subunits. The conformational relaxation following the O₂ photodissociation in the α and β subunits was found to decrease the rate constant for the geminate O₂ rebinding, this effect being more than one order of magnitude greater for the β subunits than for the α subunits. Evidence was provided for the modulation of the O₂ rebinding to the individual α and β subunits within human hemoglobin in the R-state structure by the intrinsic heme reactivity through a change in proximal constraints upon the relaxation of the tertiary structure on a picosecond to microsecond time scale. Our results demonstrate that, for native R-state oxyhemoglobin, O₂ rebinding properties and spectral changes following the O₂ photodissociation can be adequately described as the sum of those for the α and β subunits within the valency hybrids. The isolated β chains (hemoglobin H) show similar behavior to the β subunits within the valency hybrids and can be used as a model for the β subunits within the R-state oxyhemoglobin. At the same time, the isolated α chains behave differently to the α subunits within the valency hybrids.

O₂ rebinding and conformational relaxation following O₂ photodissociation were studied on picosecond to millisecond time scale in the α and β subunits within human hemoglobin in the quaternary R-like structure.     

Effects of solvent viscosity on ligand interconversion dynamics in the primary docking site of heme proteins

Interconversion dynamics of the ligand in the primary docking site of myoglobin (Mb) and hemoglobin (Hb) in trehalose and glycerol/D2O mixtures at 283 K was investigated by probing time-resolved vibrational spectra of CO photolyzed from these proteins. The interconversion dynamics in viscous media are similar to those in aqueous solution, indicating that it is minimally coupled to the solvent-coupled large-scale protein motion. Interconversion rates in the heme pocket of Hb in water solution are slower than those of Mb in trehalose glass, suggesting that the interconversion barrier in Hb is intrinsically higher than that in Mb and is not modified by the solvent viscosity.     

Diffusive dynamics on multidimensional rough free energy surfaces

The dynamics of processes relevant to chemistry and biophysics on rough free energy landscapes is investigated using a recently developed algorithm to solve the Smoluchowski equation. Two different processes are considered: ligand rebinding in MbCO and protein folding. For the rebinding dynamics of carbon monoxide (CO) to native myoglobin (Mb) from locations around the active site, the two-dimensional free energy surface (FES) is constructed using extensive molecular dynamics simulations. The surface describes the minima in the A state (bound MbCO), CO in the distal pocket and in the Xe4 pocket, and the transitions between these states and allows to study the diffusion of CO in detail. For the folding dynamics of protein G, a previously determined two-dimensional FES was available. To follow the diffusive dynamics on these rough free energy surfaces, the Smoluchowski equation is solved using the recently developed hierarchical discrete approximation method. From the relaxation of the initial nonequilibrium distribution, experimentally accessible quantities such as the rebinding time for CO or the folding time for protein G can be calculated. It is found that the free energy barrier for CO in the Xe4 pocket and in the distal pocket (B state) closer to the heme iron is approximately 6 kcal/mol which is considerably larger than the inner barrier which separates the bound state and the B state. For the folding of protein G, a barrier of approximately 10 kcal/mol between the unfolded and the folded state is consistent with folding times of the order of milliseconds.     

Structural Identification of Spectroscopic Substates in Neuroglobin

The structural origins of infrared absorptions of photodissociated CO in murine neuroglobin (Ngb) are determined by combining Fourier transform infrared (FTIR) spectroscopy and molecular dynamics (MD) simulations. Such an approach allows to identify and characterize both the different conformations of the Ngb active site and the transient ligand docking sites. To capture the influence of the protein environment on the spectroscopy and dynamics, experiments and simulations are carried out for the wild type protein and its F28L and F28W mutants. It is found that a voluminous side chain at position 28 divides site B into two subsites, B’ and B”. At low temperatures, CO in wt Ngb only migrates to site B’ from where it can rebind, and B” is not populated. The spectra of CO in site B’ for wt Ngb from simulations and experiments are very similar in spectral shift and shape. They both show doublets, red‐shifted with respect to gas‐phase CO and split by≈8 cm^?1. The FTIR spectra of the F28L mutant show additional bands which are also found in the simulations and can be attributed to CO located in substate B”. The different bands are mainly related to different orientations of the His64 side chain with respect to the CO ligand. Large red‐shifts arise from strong interactions between the Histidine? NH and the CO oxygen. After dissociation from the heme iron, the CO ligand visits multiple docking sites. The locations of the primary docking site B and a secondary site C, which corresponds to the Mb Xe4 cavity, could be identified unambiguously. Finally, by comparing experiment and simulations it is also possible to identify protonation of its ε position (His_ε64 NgbCO) as the preferred heme‐bound conformation in the wild type protein with a signal at 1935 cm^?1.     

The origin of stark splitting in the initial photoproduct state of MbCO

Ligand migration and binding in heme proteins have been measured by X-ray diffraction and time-resolved spectroscopy of photoproduct intermediates. In myoglobin (Mb), internal cavities serve as docking sites for carbon monoxide (CO) ligands. In these sites, the CO ligands display characteristic infrared (IR) stretching bands due to interactions with the local electrical field. In the primary docking site, a CO can reside in two opposite orientations, characterized by a doublet of infrared bands, B1 at approximately 2130 and B2 at approximately 2120 cm-1. To assign these bands to the specific orientations, we have reexamined the effects of mutating His64 and Val68 on the infrared stretching bands associated with the B1 and B2 photoproduct states. Wild-type, H64L, V68F, and H64L-V68F MbCO were selected for experimental and theoretical analyses. Fourier transform infrared (FTIR) spectroscopy and density functional theory (DFT) calculations were used to interpret the effects of the electrostatic environment on the B state bands. The imidazole side chain of His64 appears to be the primary cause of the observed Stark splitting. The high-frequency B1 band is assigned to the CO orientation in which the carbon (white atom) is directed toward the heme iron and the Nepsilon-H proton of His64. At low temperatures, CO molecules in the opposite orientational conformer, B2 with the O atom (red) toward His64, first rotate by 180 degrees into the more stable B1 state and then rebind.     

Orientational distribution of CO before and after photolysis of MbCO and HbCO: a determination using time-resolved polarized Mid-IR spectroscopy

The technique of time-resolved polarized mid-IR spectroscopy was used to probe the orientational distribution of carbon monoxide (CO) bound to and docked within horse myoglobin, sperm whale myoglobin, and human hemoglobin A in neutral pH solution at 283 K. An accurate determination of the orientation required that the experimentally measured polarization anisotropy be corrected for the effects of fractional photolysis in an optically thick sample. The experimental method measures the direction of the transition dipole, which is parallel to the CO bond axis when docked and nearly parallel when bound to the heme. The polarization anisotropy of bound CO is virtually the same for all protein systems investigated and is unchanging across its inhomogeneously broadened mid-IR absorption spectrum. From these results, it was concluded that the transition dipole moment of bound CO is oriented 相似文献   

Molecular docking using the molecular lipophilicity potential as hydrophobic descriptor: impact on GOLD docking performance

GOLD is a molecular docking software widely used in drug design. In the initial steps of docking, it creates a list of hydrophobic fitting points inside protein cavities that steer the positioning of ligand hydrophobic moieties. These points are generated based on the Lennard-Jones potential between a carbon probe and each atom of the residues delimitating the binding site. To thoroughly describe hydrophobic regions in protein pockets and properly guide ligand hydrophobic moieties toward favorable areas, an in-house tool, the MLP filter, was developed and herein applied. This strategy only retains GOLD hydrophobic fitting points that match the rigorous definition of hydrophobicity given by the molecular lipophilicity potential (MLP), a molecular interaction field that relies on an atomic fragmental system based on 1-octanol/water experimental partition coefficients (log P(oct)). MLP computations in the binding sites of crystallographic protein structures revealed that a significant number of points considered hydrophobic by GOLD were actually polar according to the MLP definition of hydrophobicity. To examine the impact of this new tool, ligand-protein complexes from the Astex Diverse Set and the PDB bind core database were redocked with and without the use of the MLP filter. Reliable docking results were obtained by using the MLP filter that increased the quality of docking in nonpolar cavities and outperformed the standard GOLD docking approach.     

Geminate carbon monoxide rebinding to a c-type haem

A chemically modified form of cytochrome c(cyt. c), termed carboxymethyl cytochrome c(cm cyt. c), possesses a vacant sixth coordination site to the haem iron that is available to bind external ligands. We present data on the rapid flash photolysis of CO from the ferrous haem iron of cm cyt. c and describe the kinetics and spectral transitions that accompany the recombination. This was achieved using 30-femtosecond laser pulses and a white light continuum to monitor spectral transitions. Whereas the photo-dissociation quantum yield is close to 1, the yield of CO escape from the protein (the apparent quantum yield, varphi) relative to myoglobin (varphi=1) is small due to rapid geminate recombination of CO. On ligand photo-dissociation the haem undergoes a spin-state transition from low-spin ferrous CO bound to penta-coordinate high-spin. Subsequently the system reverts to the CO bound form. The data were fitted with a minimum number of exponentials using global analysis. Recombination of CO with the haem iron of cm cyt. c is multiphasic (tau=16 ps, 120 ps and 1 ns), involving three spectrally distinct components. The fraction of haem (0.11) not recombining with CO within 4 ns is similar to the value of varphi(0.12) measured on the same preparation by the "pulse method" (M. Brunori, G. Giacometti, E. Antonini and J. Wyman, Proc. Natl. Acad. Sci. USA, 1973, 70, 3141-3144, ). This implies that no further geminate recombination occurs at t>4 ns. This unusually efficient CO-haem geminate recombination indicates the sterically hindered ("caged") nature of the distal haem pocket in cm cyt. c from which it is difficult for CO to escape. The large geminate phase may be contrasted with the behaviour of myoglobin in which geminate recombination is small. This is in general agreement with the well-documented extensive structural dynamics in myoglobin that allow ligand passage, and a higher structural rigidity in cyt. c imposed by the restraints of minimising reorganisation energy for electron transfer (M. Brunori, D. Bourgeois and D. Vallone, J. Struct. Biol., 2004, 147, 223-234, ). The high pH ferrous form of cm cyt. c is a low-spin species having a lysine bound to the central iron atom of the haem (M. Brunori, M. Wilson and E. Antonini, J. Biol. Chem., 1972, 247, 6076-6081; G. Silkstone, G. Stanway, P. Brzezinski and M. Wilson, Biophys. Chem., 2002, 98, 65-77, ). This high pH (pH approximately 8) form of deoxy cm cyt. c undergoes photo-dissociation of lysine (although the proximal histidine is possible) after photo-excitation. Recombination occurs with a time constant (tau) of approximately 7 ps. This is similar to that observed for the geminate rebinding of the Met80 residue in native ferrous cyt. c(tau approximately 6 ps) following its photo-dissociation (S. Cianetti, M. Negrerie, M. Vos, J.-L. Martin and S. Kruglik, J. Am. Chem. Soc., 2004, 126, 13 932-13 933; W. Wang, X. Ye, A. Demidov, F. Rosca, T. Sjodin, W. Cao, M. Sheeran and P. Champion, J. Phys. Chem., 2000, 104, 10 789-10 801, ).     

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO.

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.