Dynamics of Hemoglobin Investigated by Mössbauer Spectroscopy

Abstract

Mössbauer Spectra of ⁵⁷Fe enriched horse hemoglobin and sperm whale myoglobin were measured in the temperature range from 80 K to 260 K. An analysis of the temperature dependence of the recoiless fraction (the Lamb-Mössbauer factor) shows it to be sensitive to conformational fluctuations which affect the mean square displacement of the iron. We have found that the protein conformation has a dramatic effect on these measurements. For hemoglobin greater conformational fluctuations at lower temperatures are observed for carbonmonoxyhemoglobin in the liganded conformation than for deoxyhemoglobin in the unliganded conformation. On the other hand, the Lamb-Mössbauer factor is insensitive to the binding of ligands to myoglobin and shows conformational fluctuations similar to deoxyhemoglobin even in the liganded state. It is also shown that a reversible complex with the distal histidine is formed in frozen deoxyhemoglobin solutions above 200 K where the Lamb-Mössbauer factor shows the excitation of new modes of conformational fluctuations. This complex is not formed with carbonmonoxyhemoglobin which already has a sixth ligand and with deoxymyoglobin which appears to undergo much more limited conformational fluctuations. A possible relationship between the formation of the distal histidine complex and the cooperative ligand binding reaction is suggested by results with partially liganded hemoglobin which indicate increased formation of the distal histidine complex.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Ligand-induced conformational change of a protein reproduced by a linear combination of displacement vectors obtained from normal mode analysis

The conformational change of a protein upon ligand binding was examined by normal mode analysis (NMA) based on an elastic-network model (ENM) for a full-atom system using dihedral angles as independent variables. Specifically, we investigated the extent to which conformational change vectors of atoms from an apo form to a holo form of a protein can be represented by a linear combination of the displacement vectors of atoms in the apo form calculated for the lowest-frequency m normal modes (m=1, 2,…, 20). In this analysis, the latter vectors were best fitted to the former ones by the least-squares method. Twenty-two paired proteins in the holo and apo forms, including three dimer pairs, were examined. The results showed that, in most cases, the conformational change vectors were reproduced well by a linear combination of the displacement vectors of a small number of low-frequency normal modes. The conformational change around an active site was reproduced as well as the entire conformational change, except for some proteins that only undergo significant conformational changes around active sites. The weighting factors for 20 normal modes optimized by the least-squares fitting characterize the conformational changes upon ligand binding for these proteins. The conformational changes sampled around the apo form of a protein by the linear combination of the displacement vectors obtained by ENM-based NMA may help solve the flexible-docking problem of a protein with another molecule because the results presented herein suggest that they have a relatively high probability of being involved in an actual conformational change.     

Quantifying allosteric effects in proteins

In allosteric regulation, protein activity is altered when ligand binding causes changes in the protein conformational distribution. Little is known about which aspects of protein design lead to effective allosteric regulation, however. To increase understanding of the relation between protein structure and allosteric effects, we have developed theoretical tools to quantify the influence of protein-ligand interactions on probability distributions of reaction rates and protein conformations. We define the rate divergence, Dk, and the allosteric potential, Dx, as the Kullback-Leibler divergence between either the reaction-rate distributions or protein conformational distributions with and without the ligand bound. We then define Dx as the change in the conformational distribution of the combined protein/ligand system, derive Dx in the harmonic approximation, and identify contributions from 3 separate terms: the first term, D[stackxomega], results from changes in the eigenvalue spectrum; the second term, D[stackxDeltax], results from changes in the mean conformation; and the third term, Dxv, corresponds to changes in the eigenvectors. Using normal modes analysis, we have calculated these terms for a natural interaction between lysozyme and the ligand tri-N-acetyl-D-glucosamine, and compared them with calculations for a large number of simulated random interactions. The comparison shows that interactions in the known binding-site are associated with large values of Dxv. The results motivate using allosteric potential calculations to predict functional binding sites on proteins, and suggest the possibility that, in Nature, effective ligand interactions occur at intrinsic control points at which binding induces a relatively large change in the protein conformational distribution.     

Spectroscopic methods for analysis of protein secondary structure

Several methods for determination of the secondary structure of proteins by spectroscopic measurements are reviewed. Circular dichroism (CD) spectroscopy provides rapid determinations of protein secondary structure with dilute solutions and a way to rapidly assess conformational changes resulting from addition of ligands. Both CD and Raman spectroscopies are particularly useful for measurements over a range of temperatures. Infrared (IR) and Raman spectroscopy require only small volumes of protein solution. The frequencies of amide bands are analyzed to determine the distribution of secondary structures in proteins. NMR chemical shifts may also be used to determine the positions of secondary structure within the primary sequence of a protein. However, the chemical shifts must first be assigned to particular residues, making the technique considerably slower than the optical methods. These data, together with sophisticated molecular modeling techniques, allow for refinement of protein structural models as well as rapid assessment of conformational changes resulting from ligand binding or macromolecular interactions. A selected number of examples are given to illustrate the power of the techniques in applications of biological interest.     

Dynamics of hemoglobin investigated by M?ssbauer spectroscopy

M?ssbauer Spectra of Fe enriched horse hemoglobin and sperm whale myoglobin were measured in the temperature range from 80 K to 260 K. An analysis of the temperature dependence of the recoiless fraction (the Lamb-M?ssbauer factor) shows it to be sensitive to conformational fluctuations which affect the mean square displacement of the iron. We have found that the protein conformation has a dramatic effect on these measurements. For hemoglobin greater conformational fluctuations at lower temperatures are observed for carbonmonoxyhemoglobin in the liganded conformation than for deoxyhemoglobin in the unliganded conformation. On the other hand, the Lamb-M?ssbauer factor is insensitive to the binding of ligands to myoglobin and shows conformational fluctuations similar to deoxyhemoglobin even in the liganded state. It is also shown that a reversible complex with the distal histidine is formed in frozen deoxyhemoglobin solution above 200 K where the Lamb-M?ssbauer factor shows the excitation of new modes of conformational fluctuations. This complex is not formed with carbonmonoxyhemoglobin which already has a sixth ligand and with deoxymyoglobin which appears to undergo much more limited conformational fluctuations. A possible relationship between the formation of the distal histidine complex and the cooperative ligand binding reaction is suggested by results with partially liganded hemoglobin which indicate increased formation of the distal histidine complex.     

Importance of hyaluronan length in a hyaladherin-based assay for hyaluronan

Specific hyaladherin-based assays have been set up to measure the concentration of hyaluronan in biological fluids. Hyaluronectin (HN; a hyaladherin extracted from ovine brain) binds to hyaluronan (HA) that must be 10 units (HA10) or more long. It was therefore of interest to determine whether HN would continue to bind to HA10 in full-length HA since conformational changes might mask potential binding sites. We used the enzyme-linked sorbent assay (ELSA) to assay HA and hyaluronan-derived oligosaccharides, with different standard HAs, and the results were compared to results obtained with the carbazole technique. Oligosaccharide length was calculated from the ratio glucuronic acid/reducing N-acetylglucosamine in fractions of hyaluronidase-digested macromolecular hyaluronan prepared by chromatography; the size of the HA12 oligosaccharide was confirmed by matrix-assisted laser desorption ionization mass spectrometry. During the digestion of macromolecular HA with hyaluronidase, the binding of HN to HA first increased and then decreased as shown using the ELSA. The concentration of HA fragments of HA60 and below was overestimated when intact macromolecular HA was used as the reference for the ELSA, while the concentration of HA100 and above was underestimated when HA10 was used as the reference. The binding of HN to HA20, HA40, and HA60 saccharides was consistent with binding to multiples of HA10 sites. In conclusion, the level of HN binding is determined by the conformation of HA, which may mask binding sites. Hence, calibration HA used in the ELSA must be adapted to the size of HA to assay.     

Why are large conformational changes well described by harmonic normal modes?

Low-frequency normal modes generated by elastic network models tend to correlate strongly with large conformational changes of proteins, despite their reliance on the harmonic approximation, which is only valid in close proximity of the native structure. We consider 12 variants of the torsional network model (TNM), an elastic network model in torsion angle space, that adopt different sets of torsion angles as degrees of freedom and reproduce with similar quality the thermal fluctuations of proteins but present drastic differences in their agreement with conformational changes. We show that these differences are related to the extent of the deviations from the harmonic approximation, assessed through an anharmonic energy function whose harmonic approximation coincides with the TNM. Our results indicate that mode anharmonicity is more strongly related to its collectivity, i.e., the number of atoms displaced by the mode, than to its amplitude; low-frequency modes can remain harmonic even at large amplitudes, provided they are sufficiently collective. Finally, we assess the potential benefits of different strategies to minimize the impact of anharmonicity. The reduction of the number of degrees of freedom or their regularization by a torsional harmonic potential significantly improves the collectivity and harmonicity of normal modes and the agreement with conformational changes. In contrast, the correction of normal mode frequencies to partially account for anharmonicity does not yield substantial benefits. The TNM program is freely available at https://github.com/ugobas/tnm.     

Lipid fluidity and membrane protein dynamics

Membrane fluidity plays an important role in cellular functions. Membrane proteins are mobile in the lipid fluid environment; lateral diffusion of membrane proteins is slower than expected by theory, due to both the effect of protein crowding in the membrane and to constraints from the aqueous matrix. A major aspect of diffusion is in macromolecular associations: reduction of dimensionality for membrane diffusion facilitates collisional encounters, as those concerned with receptor-mediated signal transduction and with electron transfer chains. In mitochondrial electron transfer, diffusional control is prevented by the excess of collisional encounters between fast-diffusing ubiquinone and the respiratory complexes. Another aspect of dynamics of membrane proteins is their conformational flexibility. Lipids may induce the optimal conformation for catalytic activity. Breaks in Arrhenius plots of membrane-bound enzymes may be related to lipid fluidity: the break could occur when a limiting viscosity is reached for catalytic activity. Viscosity can affect protein conformational changes by inhibiting thermal fluctuations to the inner core of the protein molecule.     

Comparative molecular dynamics—Similar folds and similar motions?

Proteins possessing the same fold may undergo similar motions, particularly if these motions involve large conformational transitions. The increasing amounts of structural data provide a useful starting point with which to test this hypothesis. We have performed a total of 0.29 micros of molecular dynamics across a series of proteins within the same fold family (periplasmic binding proteinlike) in order to address to what extent similarity of motion exists. Analysis of the local conformational space on these timescales (10-20 ns) revealed that the behavior of the proteins could be readily distinguished between an apo-state and a ligand-bound state. Moreover, analysis of the root-mean-square fluctuations reveals that the presence of the ligand exerts a stabilizing effect on the protein, with similar motions occurring, but with reduced magnitude. Furthermore, the conformational space in the presence of the ligand appears to be dictated by sequence but not by the type of ligand present. In contrast, apo-simulations showed considerable overlap of conformational space across the fold as a result of their ability to undergo larger fluctuations. Indeed, we observed several transitions from different simulations between states corresponding to the closed-cleft and open-cleft forms of the fold, with the predominant motions being conserved across the different proteins. Thus, large-scale conformational changes do indeed appear to be conserved across this fold architecture, but smaller conformational motions appear to reflect the differences in sequence and local fold.     

Conformational changes and catalysis by alcohol dehydrogenase

As shown by X-ray crystallography, horse liver alcohol dehydrogenase undergoes a global conformational change upon binding of NAD⁺ or NADH, involving a rotation of the catalytic domain relative to the coenzyme binding domain and the closing up of the active site to produce a catalytically efficient enzyme. The conformational change requires a complete coenzyme and is affected by various chemical or mutational substitutions that can increase the catalytic turnover by altering the kinetics of the isomerization and rate of dissociation of coenzymes. The binding of NAD⁺ is kinetically limited by a unimolecular isomerization (corresponding to the conformational change) that is controlled by deprotonation of the catalytic zinc-water to produce a negatively-charged zinc-hydroxide, which can attract the positively-charged nicotinamide ring. The deprotonation is facilitated by His-51 acting through a hydrogen-bonded network to relay the proton to solvent. Binding of NADH also involves a conformational change, but the rate is very fast. After the enzyme binds NAD⁺ and closes up, the substrate displaces the hydroxide bound to the catalytic zinc; this exchange may involve a double displacement reaction where the carboxylate group of a glutamate residue first displaces the hydroxide (inverting the tetrahedral coordination of the zinc), and then the exogenous ligand displaces the glutamate. The resulting enzyme-NAD⁺-alcoholate complex is poised for hydrogen transfer, and small conformational fluctuations may bring the reactants together so that the hydride ion is transferred by quantum mechanical tunneling. In the process, the nicotinamide ring may become puckered, as seen in structures of complexes of the enzyme with NADH. The conformational changes of alcohol dehydrogenase demonstrate the importance of protein dynamics in catalysis.     

Conformational dynamics and thermodynamics of protein-ligand binding studied by NMR relaxation

Protein conformational dynamics can be critical for ligand binding in two ways that relate to kinetics and thermodynamics respectively. First, conformational transitions between different substates can control access to the binding site (kinetics). Secondly, differences between free and ligand-bound states in their conformational fluctuations contribute to the entropy of ligand binding (thermodynamics). In the present paper, I focus on the second topic, summarizing our recent results on the role of conformational entropy in ligand binding to Gal3C (the carbohydrate-recognition domain of galectin-3). NMR relaxation experiments provide a unique probe of conformational entropy by characterizing bond-vector fluctuations at atomic resolution. By monitoring differences between the free and ligand-bound states in their backbone and side chain order parameters, we have estimated the contributions from conformational entropy to the free energy of binding. Overall, the conformational entropy of Gal3C increases upon ligand binding, thereby contributing favourably to the binding affinity. Comparisons with the results from isothermal titration calorimetry indicate that the conformational entropy is comparable in magnitude to the enthalpy of binding. Furthermore, there are significant differences in the dynamic response to binding of different ligands, despite the fact that the protein structure is virtually identical in the different protein-ligand complexes. Thus both affinity and specificity of ligand binding to Gal3C appear to depend in part on subtle differences in the conformational fluctuations that reflect the complex interplay between structure, dynamics and ligand interactions.     

The binding affinity of human serum albumin and paclitaxel through MMPBSA based on docked complex

The distribution, free concentration and metabolism of drugs can be significantly altered as a result of binding to albumin. At the same time, the conformational of serum albumin was also changed by interaction with low molecular weight drugs. In present work, we first equilibrated HSA in aqueous solution to obtain the solvated-HSA model. Further solvated-HSA was performed molecular docking with paclitaxel to find the binding sites. The two docking HSA-paclitaxel complexes were obtained and further equilibrated by a 12 ns MD simulation. Then, MMPBSA method was used to investigate the binding free energy of them. Finally, we correlated the fluctuations of residues with corresponding changes in the secondary structure by dssp method. Two binding sites of paclitaxel were found on HSA having considered the solvation effect. More hydrogen bonds were formed at site I respected to site II. A larger binding energy for primary binding also indicated that paclitaxel showed higher binding affinity mainly due to the stronger hydrogen bonding interactions. There was a significant difference between the two complexes on structure according to the dssp results. Moreover, structure of the binding sites exhibited more fluctuations after binding paclitaxel compared with other regions. Paclitaxel binding also induced distinct conformational changes in drug binding site even when it was empty and have contributed to a reduced binding capacity of HSA towards adriamycin.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

On the link between conformational changes,ligand binding and heat capacity

BackgroundConformational changes coupled to ligand binding constitute the structural and energetics basis underlying cooperativity, allostery and, in general, protein regulation. These conformational rearrangements are associated with heat capacity changes. ITC is a unique technique for studying binding interactions because of the simultaneous determination of the binding affinity and enthalpy, and for providing the best estimates of binding heat capacity changes.Scope of reviewStill controversial issues in ligand binding are the discrimination between the “conformational selection model” and the “induced fit model”, and whether or not conformational changes lead to temperature dependent apparent binding heat capacities. The assessment of conformational changes associated with ligand binding by ITC is discussed. In addition, the “conformational selection” and “induced fit” models are reconciled, and discussed within the context of intrinsically (partially) unstructured proteins.Major conclusionsConformational equilibrium is a major contribution to binding heat capacity changes. A simple model may explain both conformational selection and induced fit scenarios. A temperature-independent binding heat capacity does not necessarily indicate absence of conformational changes upon ligand binding. ITC provides information on the energetics of conformational changes associated with ligand binding (and other possible additional coupled equilibria).General significancePreferential ligand binding to certain protein states leads to an equilibrium shift that is reflected in the coupling between ligand binding and additional equilibria. This represents the structural/energetic basis of the widespread dependence of ligand binding parameters on temperature, as well as pH, ionic strength and the concentration of other chemical species. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Direct product-matrix method treatment of macromolecular binding

A major limitation of the well-known matrix method treatment of linear polymer systems is the large size of the transfer matrix required to account simultaneously for such factors as polymer conformation, polymer sequence, ligand binding, etc. This paper shows how to combine such factors by using a direct product of small matrices, one for each factor, to generate the large transfer matrix. The properties of the direct product operation can be used to simplify considerably the calculation of eigenvalues and eigenvectors of the transfer matrix. This permits one to develop fairly simply the statistical thermodynamics of relatively complicated macromolecular binding systems involving a large number of states for the individual polymer sites. Applications include treatments of sequence-specific binding of ligands, sequence-specific conformational changes, and the combination of the two processes.     

Conformational transitions in adenylate kinase. Allosteric communication reduces misligation

Large conformational changes in the LID and NMP domains of adenylate kinase (AKE) are known to be key to ligand binding and catalysis, yet the order of binding events and domain motion is not well understood. Combining the multiple available structures for AKE with the energy landscape theory for protein folding, a theoretical model was developed for allostery, order of binding events, and efficient catalysis. Coarse-grained models and nonlinear normal mode analysis were used to infer that intrinsic structural fluctuations dominate LID motion, whereas ligand-protein interactions and cracking (local unfolding) are more important during NMP motion. In addition, LID-NMP domain interactions are indispensable for efficient catalysis. LID domain motion precedes NMP domain motion, during both opening and closing. These findings provide a mechanistic explanation for the observed 1:1:1 correspondence between LID domain closure, NMP domain closure, and substrate turnover. This catalytic cycle has likely evolved to reduce misligation, and thus inhibition, of AKE. The separation of allosteric motion into intrinsic structural fluctuations and ligand-induced contributions can be generalized to further our understanding of allosteric transitions in other proteins.