Search

GK (glucokinase) is activated by glucose binding to its substrate site, is inhibited by GKRP (GK regulatory protein) and stimulated by GKAs (GK activator drugs). To explore further the mechanisms of these processes we studied pure recombinant human GK (normal enzyme and a selection of 31 mutants) using steady-state kinetics of the enzyme and TF (tryptophan fluorescence). TF studies of the normal binary GK–glucose complex corroborate recent crystallography studies showing that it exists in a closed conformation greatly different from the open conformation of the ligand-free structure, but indistinguishable from the ternary GK–glucose–GKA complex. GKAs did activate and GKRP did inhibit normal GK, whereas its TF was doubled by glucose saturation. However, the enzyme kinetics, GKRP inhibition, TF enhancement by glucose and responsiveness to GKA of the selected mutants varied greatly. Two predominant response patterns were identified accounting for nearly all mutants: (i) GK mutants with a normal or close to normal response to GKA, normally low basal TF (indicating an open conformation), some variability of kinetic parameters (kcat, glucose S0.5, h and ATP Km), but usually strong GKRP inhibition (13/31); and (ii) GK mutants that are refractory to GKAs, exhibit relatively high basal TF (indicating structural compaction and partial closure), usually show strongly enhanced catalytic activity primarily due to lowering of the glucose S0.5, but with reduced or no GKRP inhibition in most cases (14/31). These results and those of previous studies are best explained by envisioning a common allosteric regulator region with spatially non-overlapping GKRP- and GKA-binding sites.

We measured the low-wavenumber polarized resonance Raman spectra of horse heart (hhc), chicken (chc) and yeastC102T (yc) ferrocytochromes c with Soret excitation. We examined the out-of-plane (oop) deformations of the heme groups by virtue of relative intensities and depolarization ratios of a variety of oop and in-plane (ip) Raman active bands. Analysis of relative Raman intensities shows differences in deviation from planarity of the heme groups of yeast, horse heart and chicken cytochromes c. The heme groups in cytochrome c proteins have been shown by normal coordinate static deformation (NSD) analysis from crystal structures to exhibit a dominant ruffling (B1u) deformation. As a consequence the B1u modes, γ10 − γ12, become resonance Raman active. We used normalized Raman intensity ratios and depolarization ratios of oop Raman active modes, whose intensities are attributable to specific nonplanar deformations, to estimate and compare their Franck-Condon-type and Jahn-Teller-type coupling magnitudes for horse heart, chicken and yeast ferrocytochrome c at neutral pH. These coupling magnitudes allow for a quantitative comparison of oop deformations between individual heme groups. Chicken ferrocytochrome was found to have the largest ruffling deformation of the three investigated proteins, followed by horse heart and yeast cytochrome c. The heme group of the former is slightly more ruffled than the corresponding active site of the latter, while saddling in both proteins is substantially larger than in chicken ferrocytochrome c. The Raman data are sensitive enough to allow a comparison of lesser deformations. Doming, which is a kinetic coordinate in many heme proteins, is largest in chicken and smallest in yeast cytochrome c. Waving is largest in yeast, followed by horse heart and chicken cytochrome c. Propellering deformations could be compared for chicken and horse heart cytochrome c and were found to be substantially larger in the latter. A comparison with heme deformations obtained from X-ray structures (for horse heart and yeast cytochrome c) and from molecular dynamics simulations (MDS) (performed for all three proteins) yields some agreement with the main ruffling and saddling deformations derived from the crystal structures, whereas the heme conformations produced by MDS seem to account better for smaller deformations like doming and propellering. The present study demonstrates the usefulness of resonance Raman spectroscopy for the analysis of nonplanar deformations in heme proteins. Copyright © 2008 John Wiley & Sons, Ltd.

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

We report results of 12 ns, all-atom molecular dynamics simulation (MDS) and Poisson-Boltzmann free energy calculations (PBFE) on calmodulin (CaM) bound to two molecules of trifluoperazine (TFP) and of N-(3,3, diphenylpropyl)- N'-[1- R-(3,4-bis-butoxyphenyl)-ethyl]-propylenediamine (DPD). X-ray data show very similar structures for the two complexes, yet the antagonists significantly differ with respect to their CaM binding affinities, the neutral DPD is much more potent. The goal of the study was to unravel the reason why TFP is less potent although its positive charge should facilitate binding. The electrostatic energy terms in CHARMM and binding free energy terms of the PBFE approach showed TFP a better antagonist, while inspection of hydrophobic contacts supports DPD binding. Detailed inspection of the amino acid contributions of PBFE calculations unravel that steric reasons oppose the favorable binding of TFP. Structural conditions are given for a successful drug design strategy, which may benefit also from charge-charge interactions.

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations—in which interdomain interactions cannot occur—the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.

The molecular details of the mechanism of action of allosteric effectors on hemoglobin oxygen affinity are not clearly understood. The global allostery model proposed by Yonetani et al. suggests that the binding of allosteric effectors can take place both in the R and T states and that they influence oxygen affinity through inducing global tertiary changes in the subunits. Recently published high pressure studies yielded dissociation constants at atmospheric pressure that showed a stabilizing effect of heterotropic allosteric effectors on the dimer interface in the R state, and a more pronounced destabilizing effect in a T state model. In the present work, we report on computational modeling used to interpret the high pressure experimental data. We show structural changes in the hemoglobin interdimeric interfaces, indicative of a global tertiary structural change induced by the binding of allosteric effectors. We also show that the number of water molecules bound at the interface is significantly influenced by binding effectors in the T state in accordance with the experimental data. Our results suggest that the binding of effectors at definite sites leads to tertiary changes that propagate to the interfaces and results in overall structural re-organizations.

Resonance Raman spectroscopy is used to probe the effect of calcium depletion on the heme group of horseradish peroxidase C at pH 8. Polarized Raman spectra are recorded with an argon ion laser at eight different wavelengths to provide a sound database for a reliable spectral decomposition. Upon calcium depletion, the spectrum is indicative of a predominantly pentacoordinated high spin state of the heme iron coexisting with small fractions of hexacoordinated high and low spin states. The dominant quantum mixed spin state of native ferric horseradish peroxidase, which is characteristic for class III peroxidases, is not detectable in the spectrum of the enzyme with partial distal Ca(2+) depletion. The quenching of the quantum mixed spin state and the predominance of the pentacoordinated high spin state are likely to arise from distortions induced by distal calcium depletion, which translates into a weaker Fe-N(epsilon)(His) bond and a more tilted imidazole. A correlation is proposed between the lower enzyme activity and the elimination of the pentacoordinated quantum mixed state upon Ca(2+) depletion.

The effect of trehalose on the interaction of human serum albumin (HSA) with neutral and negatively charged small unilamellar vesicles (SUVs) composed of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) or of mixtures of DMPC (19:1 w/w) with 1,2-dimyristoyl-sn-glycero-3-phosphatidylglycerol (DMPG) was studied by time-resolved fluorescence and dynamic light scattering measurements. The results are interpreted with supporting nonbond calculations describing the nonbond domains most likely to be involved in the protein-SUV interaction. In the absence of trehalose, lifetime measurements of the single Trp of HSA are indicative of two different SUV-HSA associative mechanisms depending on the [lipid]/[HSA] concentration ratios. At low ratios, depletion of phospholipid molecules from vesicles by HSA occurs independently of the lipid composition of the vesicles via favorable hydrophobic contacts. At higher ratios, vesicle-HSA assocation is favored by electrostatic interactions for the negatively charged SUNs. For neutral SUNs, hydrophobically driven penetration of HSA is proposed. All association mechanisms are damped in the presence of trehalose, due to its capacity to coat the interacting surfaces. The results of dynamic light scattering experiments clearly show that the aging of the liposomes is dependent on the lipid composition. The aging of DMPC vesicles is faster and not affected by the presence of either HSA or trehalose. The aging of DMPC/DMPG liposomes is more pronounced in the presence of HSA. These SUVs are stabilized by trehalose through different mechanisms depending on whether they are covered by HSA or not.

Horseradish peroxidase C binds a wide variety of small H-donor compounds such as benzohydroxamic acid (BHA) and 2-naphthohydroxamic acid (NHA). In this work, we use the Mg(II)-mesoporphyrin prosthetic group derivative as a spectroscopic probe of the active site and of the interaction with the substrates. We report on high-resolution fluorescence line-narrowed spectra which show that the effects of substrate binding on the electronic transitions are similar for both substrates and present data on the normal vibrational modes that are active in the vibronic spectra. Analysis of the vibrational frequencies shows that the Mg(II) ion is 5-coordinate in all cases, thus ruling out a solvent water as sixth ligand. The frequency shifts observed as a result of substrate binding are also indicative of a more rigid prosthetic group upon substrate binding. We present models for MgMP-HRP and its complexes with both substrates and compare the resulting structures on the basis of a modeling approach combining energy minimization to finite difference Poisson--Boltzmann calculations which partitions the various relative protein contributions to substrate binding. We show that the electrostatic potential of the prosthetic group is modified by the binding event. Analysis of the unbound and bound energy-minimized structures shows that the enzyme modulates substrate binding by subtle charge reorganization in the vicinity of the catalytic site and that this rearrangement is not attributable to significant secondary structure conformational changes but to side-chain reorganization.

Horseradish peroxidase C is an oxidoreductase which catalyzes in plant roots the oxidation of a remarkably wide variety of aromatic compounds by H2O2. The recently available X-ray structures of the enzyme bound to aromatic substrates are not indicative of significant structural rearrangements as a result of substrate binding when compared to the structure of the unbound enzyme. Our most recent spectral hole-burning studies on HRPC fluorescent derivatives provided direct experimental evidence indicative of different fields experienced at the heme as a result of substrate binding which can only originate in the protein matrix in which the heme is embedded. In this report, we present results on modeling the fluctuations of these protein matrix electrostatic rearrangements using a combination of molecular dynamics and electrostatic conformational sampling to compare the binding of BHA and NHA to HRPC and to the derivative used in our experimental studies. Our initial results suggest that HRPC does undergo conformational rearrangement on binding aromatic substrates, but of an electrostatic nature, as opposed to major conformational structural changes. This implies that the high versatility of substrate binding in horseradish peroxidase could be due to an electrostatic dynamic selectivity mediated by small-amplitude internal protein matrix motions. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 290–301, 2001

The binding of acetyl­pepstatin to the Q7K/L33I/L63I mutant of HIV-1 protease was studied by fluorescence, phosphorescence, and 500-ps molecular dynamics. The protease is a homodimer with two tryptophans per monomer. Maximum entropy method (MEM) analysis and acrylamide quenching results show two tryptophyl, tryptophan (Trp) populations in the apoenzyme that merge into one in the complex. These results are in agreement with molecular dynamics simulations indicative of Trp asymmetry in the apoenzyme as revealed by the occurrence of nonequivalent Trp42 indole rotamer interconversions, not observed for the complex. Analysis of the local Trp42B environments of the apoenzyme with respect to possible quencher groups shows that the c2 interconversions do not influence the lifetime, while the c1 interconversions do. Upon binding the inhibitor, Trp42B acquires a single conformation with the same lifetime and orientation as that of Trp42, and also with less quenching accessibility. Thus, protein conformational dynamics become constrained with inhibitor binding. This conclusion is supported by red-edge effect experiments and phosphorescence lifetime measurements. The low temperature tp (~5.8 s) is quenched to ~200 ms as protein motions become activated around the glass transition temperature. In the case of the complex, the phosphorescence lifetime data show a more cooperative activation of the quenching mechanisms.

To clarify the role of metal ion coordination in horseradish peroxidase C (HRPC), the effect of pressure and of an externally applied electric field on spectral holes was compared for both metal-free and Mg-mesoporphyrin-substituted horseradish peroxidase C (MP-HRP and MgMP-HRP), as affected by the binding of 2-naphthohydroxamic acid (NHA). The data are compared to earlier studies performed on the same derivatives. Results obtained for MP-HRP show the presence of a predominant MP tautomer, as well as that of another small population with different pocket field and isothermal compressibility (0.12 vs 0.24 GPa(-1)). Binding NHA induces the formation of two new almost equal populations of MP-HRP tautomer complexes and the protein compressibility in both forms is increased to 0.50 and 0.36 GPa(-1). The protein structure becomes much softer than in the absence of NHA. Binding the same substrate to MgMP-HRP resulted in MgMP adopting a single conformation with no compressibility changes, while without NHA, two forms were possible. Stark effect results show charge rearrangement upon substrate binding in both cases. We propose that it is the presence of the metal that stabilizes the structure during the reorganization of the protein matrix induced by the substrate binding event. With the metal, only one conformation is adopted, without significant structural rearrangement but with charge redistribution. The dissociation constants determined for NHA binding to both derivatives and to native HRPC show that studies using mesoporphyrin and Mg-mesoporphyrin derivatives are relevant to investigating the specificity of the substrate-binding pocket in this enzyme.

The Q7K/L331/L631 HIV-1 protease mutant was expressed in Escherichia coli and the effect of binding a substrate-analog inhibitor, acetyl-pepstatin, was investigated by fluorescence spectroscopy and molecular dynamics. The dimeric enzyme has four intrinsic tryptophans, located at positions 6 and 42 in each monomer. Fluorescence spectra and acrylamide quenching experiments show two differently accessible Trp populations in the apoenzyme with k(q1) = 6.85 x 10(9) M(-1) s(-1) and k(q2) = 1.88 x 10(9) M(-1) s(-1), that merge into one in the complex with k(q) = 1.78 x 10(9) M(-1) s(-1). 500 ps trajectory analysis of Trp X1/X2 rotameric interconversions suggest a model to account for the observed Trp fluorescence. In the simulations, Trp6/Trp6B rotameric interconversions do not occur on this timescale for both HIV forms. In the apoenzyme simulations, however, both Trp42s and Trp42Bs are flipping between X1/X2 states; in the complexed form, no such interconverions occur. A detailed investigation of the local Trp environments sampled during the molecular dynamics simulation suggests that one of the apoenzyme Trp42B rotameric interconversions would allow indole-quencher contact, such as with nearby Tyr59. This could account for the short lifetime component. The model thus interprets the experimental data on the basis of the conformational fluctuations of Trp42s alone. It suggests that the rotameric interconversions of these Trps, located relatively far from the active site and at the very start of the flap region, becomes restrained when the apoenzyme binds the inhibitor. The model is thus consistent with associating components of the fluorescence decay in HIV-1 protease to ground state conformational heterogeneity.

We report on a comparative investigation of the heme pocket fields of two Zn-substituted c -type cytochromes—namely yeast and horse heart cytochromes c —using a combination of hole burning Stark spectroscopy and electrostatic calculations. The spectral hole burning experiments are consistent with different pocket fields experienced at the hemes of the respective cytochromes. In the case of horse heart Zn-cytochrome c , two distinguishable electronic origins with different electrostatic properties are observed. The yeast species, on the other hand, displays a single electronic origin. Electrostatic calculations and graphics modeling using the linearized finite-difference Poisson-Boltzmann equation performed at selected time intervals on nanosecond-molecular dynamics trajectories show that the hemes of the respective cytochromes sample different potentials as they explore conformational space. The electrostatic potentials generated by the protein matrix at the heme show different patterns in both cytochromes, and we suggest that the cytochromes differ by the number of "electrostatic substates" that they can sample, thus accounting for the different spectral populations observed in the two cytochromes.

Resonance Raman spectra with Q-band excitation are reported for microperoxidase-11, the cytochrome c analog. Spectra were acquired in the mid-frequency range for the oxidized, and reduced forms of the undecapeptide, as well as for the imidazole and carbonyl complexes. Oxidation and spin state marker bands of the undecapeptides are consistent with a six-coordinate, low spin iron in both oxidation states. Porphyrin core size correlations yield a porphyrin-centre to pyrrole-nitrogen distance of 2.00 A for MP11, suggestive of a six-coordinate species in a distorted heme environment. Molecular dynamics results show that the non-planarity of the heme of the parent cytochrome is conserved in the microperoxidase and its carbonmonoxy analog.

We performed hole burning Stark effect experiments on a cytochrome c type protein. From the behavior of the holes in an external electric field conclusions on local electric fields in the heme pocket can be drawn. We varied the pH value in order to study the influence of different groups on the pocket field: (a) titratable amino acid residues, (b) the porphyrin side chains (especially the propionates) and (c) the two axial ligands His18 and Met80. Surprisingly, a variation of the pH between 4.1 and 12.7 has no influence on the Stark effect. However, at a pH-level

The effect of the orientation of CO bound to the heme of two c-type cytochromes on the CO-stretch frequency νCO has been investigated using molecular mechanics and finite difference Poisson−Boltzmann calculations. Our approach treats the charge distribution of the protein as an external electric field capable of inducing Stark frequency shifts. Results show that modifying the Fe−C−O bending angle (β) does not change the CO stretch frequency within a range of 100−175°, equivalent to a bending motion away from the propionic acid chains. The calculated Stark shifts range from 5.4 to 8.6 cm-1 and are in good agreement with the experimentally observed shift (6 cm-1). However, motion of the CO toward the propionics exerts significant influence on the calculated shifts (−4.5 to −1.8 cm-1), which are then in total disagreement with experiment, not only in magnitude but also in predicting the wrong direction for the shift. The Stark shifts calculated for the tilt angle (τ) show that it has no significant effect on...

We report electrostatic potential and electric field calculations at the CO ligand of the heme in both horse heart and yeast cytochrome c, obtained from a finite difference solution to the Poisson−Boltzmann equation. This method takes into account the protein shape and charge distribution, as well as the solvent and generalized ionic strength effects. The calculations support recent experimental and theoretical evidence suggesting that polar interactions can significantly affect the vibrational frequency of the ligand. Our work shows that the ionizable amino acid residues and polar contributions of the protein matrix can induce a Stark effect on the CO stretch frequency of the carbonmonoxycytochromes. The observed CO stretches, at 1965.9 and 1960.1 cm-1 for horse and yeast cytochromes c, respectively, show that the experimental shift is of the order of 6 cm-1. This is in good agreement with a calculated value of 8 cm-1 (±1 cm-1) for a vibrational Stark shift due to the different charge distributions in bo...

Fluorescence-line-narrowing (FLN) spectra are presented for Sn cytochrome c (Snc) obtained under Q0-0 and Q0-1 excitation and for Zn cytochrome c (Znc) under Q0-1 excitation. Vibrational frequencies of the ground and first excited states are reported for Snc. Additional excited-state Znc frequencies are also presented, thus contributing to a previous FLN Znc investigation. The spectra of Snc shows more phonon broadening than for Znc, but vibrational frequencies of the S0 and S1 states could be obtained and the distribution function of the 0, 0 transition was determined. Snc and Znc showed good agreement in their excited-state vibrational frequencies. Some S1 vibrational frequencies of Znc are seen to downshift with respect to ground-state vibrations, suggesting that the porphyrin S0 to S1 expansion observed in other porphyrin systems also occurs in cytochrome c derivatives. © 1995 John Wiley & Sons, Inc.

Molecular dynamics simulation of explicitly solvated horseradish peroxidase and of its Ca-depleted form have been carried out, and the trajectories have been analyzed by the essential dynamics method. The results indicate that the motion of the native species is defined by a few preferred directions identified by the first four eigenvectors. The eigenvectors are significantly sampled and reveal that collective motions are perturbed in the absence of calcium. The destabilization of HRP and the corresponding decrease in the catalytic activity of the enzyme is due to perturbed collective motions primarily in the region located around the proximal calcium site.

We measured the low‐frequency polarized resonance Raman spectra of horse heart, chicken, and yeast(C102T) ferrocytochromes c with Soret excitation. We examined the out‐of‐plane deformations of the heme groups by determining the relative intensities and depolarization ratios of a variety of out‐of‐plane and in‐plane Raman active bands. Analysis of relative Raman intensities shows differences in non‐planarity of the heme groups of yeast(C102T), horse heart and chicken cytochrome c. Cytochrome c has been shown to have a dominant ruffling (B1u) deformation by means of normal coordinate structural decomposition (NSD) analysis of the heme group in crystal structures. The presence and intensity of B1u modes, γ10‐γ12, support the indication of ruffling being the major contribution to the non‐planar deformations in cytochrome c. Other types of non‐planar deformations like doming (A2U) and waving (Eg) can be deduced from the Raman activity of γ5 (A2u), γ21 and γ22 (Eg). The depolarization ratios of γ5, γ10, γ11 and...

The recently discovered new members of the globin family, neurogobin and cytoglobin, are the object of sustained structural and functional studies aimed at understanding their physiological role and elucidating the impact of their bis-his heme hexacoordination. However, no studies have yet considered the dynamics of this protein family, an essential link between structure and function. In this communication, we present normal mode analysis results for neuroglobin, cytoglobin, hemoglobin and myoglobin to provide exploratory insights into globin characteristic motions. Our results show a clear correlation in the protein dynamics of this family. All four globins exhibit a high degree of correlated displacements involving residues in the C, E and F helices and link regions. They suggest that these motions play an important role in the reversible oxygen binding function of these proteins. Further, our results may help rationalize some functional features of the 6c-globins in that they alone exhibit correlated displacements of the G-helix region.

We have measured the low temperature (T = 20 K) absorption spectra of the N52A, N52V, N52I, Y67F, and N52AY67F mutants of ferrous Saccharomyces cerevisiae (baker's yeast) cytochrome c. All the bands in the Q0- and Q(v)-band region are split, and the intensity distributions among the split bands are highly asymmetric. The spectra were analyzed by a decomposition into Voigtian profiles. The spectral parameters thus obtained were further analyzed in terms of the vibronic coupling model of Schweitzer-Stenner and Bigman (Schweitzer-Stenner, R.; Bigman, D. J. Phys. Chem. B 2001, 7064-7073) to identify parameters related to electronic and vibronic perturbations of the heme macrocycle. We report that the electronic perturbation is of B(1g) symmetry and reflects the heterogeneity of the electric field at the heme, that is, the difference between the gradients along the perpendicular N-Fe-N axis of the heme core. We found that all the investigated mutations substantially increase this electronic perturbation, so that the spectral properties become similar to those of horse heart cytochrome c. Moreover, the electronic perturbation was found to correlate nonlinearly with the enthalpy changes associated with the reduction of the heme iron. Group theoretical arguments are invoked to propose a simple model which explains how a perturbation of the obtained symmetry can stabilize the reduced state of the heme iron. Finally, vibronic coupling parameters obtained from the analysis of the Q(v)-band region suggest that the investigated mutations decrease the nonplanar deformations of the heme group. This finding was reproduced by a normal mode structural decomposition (NSD) analysis of the N52V and N52VY67F heme conformations obtained from a 1 ns molecular dynamics simulation. We argue that the reduced nonplanarity contributes to the stabilization of the reduced state.

We have measured and analyzed the low-temperature (T=10 K) absorption spectrum of reduced horse heart and yeast cytochrome c. Both spectra show split and asymmetric Q(0) and Q(upsilon) bands. The spectra were first decomposed into the individual split vibronic sidebands assignable to B(1g) (nu15) and A(2g) (nu19, nu21, and nu22) Herzberg-Teller active modes due to their strong intensity in resonance Raman spectra acquired with Q(0) and Q(upsilon) excitations. The measured band splittings and asymmetries cannot be rationalized solely in terms of electronic perturbations of the heme macrocycle. On the contrary, they clearly point to the importance of considering not only electronic perturbations but vibronic perturbations as well. The former are most likely due to the heterogeneity of the electric field produced by charged side chains in the protein environment, whereas the latter reflect a perturbation potential due to multiple heme-protein interactions, which deform the heme structure in the ground and excited states. Additional information about vibronic perturbations and the associated ground-state deformations are inferred from the depolarization ratios of resonance Raman bands. The results of our analysis indicate that the heme group in yeast cytochrome c is more nonplanar and more distorted along a B(2g) coordinate than in horse heart cytochrome c. This conclusion is supported by normal structural decomposition calculations performed on the heme extracted from molecular-dynamic simulations of the two investigated proteins. Interestingly, the latter are somewhat different from the respective deformations obtained from the x-ray structures.

The lipid bilayer is the basic component of biological membranes. It functions as a relatively impermeable barrier for most water-soluble molecules. A large amount of proteins are intrinsic to the membrane or partially bound and mediate transport of specific molecules or act as receptors. The molecular modeling of lipid bilayers and proteins has now truly come of age with the widespread availability of parallel processing and advances in modeling methodologies. This chapter will discuss recent advances in specifically modeling protein–lipid interactions focusing on two computational approaches: Molecular Dynamics (MD) and Monte Carlo (MC) simulation. MD is a choice method because the protein and lipid atoms are represented explicitly as well as solvation water and the whole system is followed in time, thus providing insights on both spatial organization and temporal dynamics. Even though MC cannot provide time-dependent quantities, it represents a powerful method for sampling conformations and locations of a single peptide and also for simulating lateral diffusion of membrane lipids and proteins as well as the coupling of these distributions.

We performed a docking study followed by a 500-ps molecular dynamics simulation of R-state human adult hemoglobin (HbA) complexed to different heterotropic effectors [2,3-diphosphoglycerate (DPG), inositol hexaphosphate (IHP), and 2-[4-[(3,5-dichlorophenylcarbamoyl)-]methyl]-phenoxy]-2-methylpropionic acid (RSR13)) to propose a molecular basis for recently reported interactions of effectors with oxygenated hemoglobin. The simulations were carried out with counterions and explicit solvation. As reported for T-state HbA, the effector binding sites are also located in the central cavity of the R-state and differ depending on effector anionic character. DPG and IHP bind between the alpha-subunits and the RSR13 site spans the alpha1-, alpha2- and beta2-subunits. The generated models provide the first report of the molecular details of R-state HbA bound to heterotropic effectors.

Horseradish peroxidase C (HRPC) is a class III peroxidase whose structure is stabilized by the presence of two endogenous calcium atoms. Calcium removal has been shown to decrease the enzymatic activity of the enzyme. The spin state of the iron, a mixture of high spin (HS) and mixed quantum spin state (QS) consisting of intermediate spin (IS) 3/2 + (HS) 5/2, is also significantly affected by calcium removal, going from a predominant QS component to a predominant HS component upon removal of one calcium. Removal of both calcium ions, however, results in the appearance of a significant LS contribution, easily monitored in the charge transfer (CT) band region by low-T absorption. Normal structural decomposition (NSD) calculations of the in-plane (ip) modes of the heme extracted from HRPC native and Ca2+-depleted models show that removal of the proximal calcium is associated with perturbed Eu and increased A1g ip distortions of the heme. The effect of complete or distal calcium removal on the heme also results in increased A1g ip distortions, but in significantly decreased Eu distortions. The overall effect is to decrease the nonplanarity of the heme: the total ip distortion of the native HRPC heme is 0.200 and 0.134 A for the Ca2+-depleted species. Our NSD results corroborate the role proposed for the protein matrix, namely to fine-tune the active site by inducing subtle changes in heme planarity and spin state of the iron. © 2004 Wiley Periodicals, Inc. Biopolymers, 2004

The distortions of the alpha1, alpha2, beta1, and beta2 hemes of human hemoglobin (HbA) in various quaternary states and as affected by the presence of allosteric effectors was investigated by subjecting CHARMM energy-minimized models to normal coordinate structural decomposition (NSD) analysis. NSD was applied to the individual hemes extracted from the R, T, and R2-state models of HbA and to HbA bound to DPG and to IHP. Overall, NSD results are indicative of characteristic distortions, not only for the hemes of the different HbA quaternary states, but also for the hemes of the HbA models bound to allosteric effectors. Comparing the distortions of the inequivalent alpha and beta hemes in T-state HbA, we show good correlation between NSD and the experimentally observed low-frequency nu52 (Eg) and gamma7 (A2u) modes reported in the literature for alpha and beta HbA hemes while noting substantial differences between these types for B2u and B1u distortions. For the R2 hemes, NSD yields heme distortions that are more comparable to those of the R-state, especially in magnitude. However, the R2 hemes do not exhibit inequivalence of alpha and beta heme distortions, a result that may contribute to an understanding of the functional importance of this state. Relative to T-state heme distortions, NSD results on the effector-bound hemes show that tertiary changes induced in T-state HbA as a result of binding DPG and IHP drastically affect heme distortions. In the alpha hemes extracted from the HbA-DPG model, most noteworthy are the increased wav(x) and wav(y) distortions and enhancement of ruf and dom deformations. In the beta hemes, the wav(y) is the most affected distortion with increase in sad. The NSD results are also different for the hemes of the HbA-IHP model, in that the beta sad and ruf deformations are more enhanced with increase of doming in the alpha hemes. Our results describe the effect of the subtle protein-induced changes on the nonplanarity of the HbA hemes that may play a role in the regulation of their oxygen affinity.

Horseradish peroxidase C (HRPC) binds 2 mol calcium per mol of enzyme with binding sites located distal and proximal to the heme group. The effect of calcium depletion on the conformation of the heme was investigated by combining polarized resonance Raman dispersion spectroscopy with normal coordinate structural decomposition analysis of the hemes extracted from models of Ca(2+)-bound and Ca(2+)-depleted HRPC generated and equilibrated using molecular dynamics simulations. Results show that calcium removal causes reorientation of heme pocket residues. We propose that these rearrangements significantly affect both the in-plane and out-of-plane deformations of the heme. Analysis of the experimental depolarization ratios are clearly consistent with increased B(1g)- and B(2g)-type distortions in the Ca(2+)-depleted species while the normal coordinate structural decomposition results are indicative of increased planarity for the heme of Ca(2+)-depleted HRPC and of significant changes in the relative contributions of three of the six lowest frequency deformations. Most noteworthy is the decrease of the strong saddling deformation that is typical of all peroxidases, and an increase in ruffling. Our results confirm previous work proposing that calcium is required to maintain the structural integrity of the heme in that we show that the preferred geometry for catalysis is lost upon calcium depletion.

Knowledge of the interactions involving charged, polar and polarizable groups in proteins is fundamental, not only because they are important determinants for gaining insight into biophysical molecular recognition and assembly processes, but also for understanding how the matrix of a protein can be viewed as an electric field capable of inducing Stark perturbations on the spectral properties of biological optical centers. This review describes the essential features of non-covalent interactions in protein systems and discusses the concept of the dielectric constant of a protein in the context of different microscopic and macroscopic modeling approaches. It also provides an account of a specific type of high resolution vibrational and optical Stark spectroscopy attempting to correlate the observed spectral properties of biological optical centers to the intrinsic protein fields induced by the matrix in which they reside.

Fluorescence line narrowing is a high resolution spectroscopic technique that uses low temperature and laser excitation to optically select specific subpopulations from the inhomogeneously broadened absorption band of the sample. When applied to the study of fluorescent groups in proteins one can obtain vibronically resolved spectra, which can be analyzed to give information on spectral line shapes, vibrational energies of both the ground and excited state molecule, and the inhomogeneous distribution function of the electronic transitions. These parameters reveal information about the chromophoric prosthetic group and the protein matrix and are functions of geometric strains and local electric fields imposed by the protein. Examples of the use of fluorescence line narrowing are discussed in investigations of heme proteins, photosynthetic systems and tryptophan-containing proteins.

The absorption and emission properties of hybrid Zn-substituted human hemoglobin (Hb) were used to monitor differences in interaction between the porphyrin and the polypeptide chain for the two subunits. Although alpha-substituted (alpha-ZnHb), beta-substituted (beta-ZnHb), or totally substituted Hb all show optical properties characteristic of Zn porphyrins, the spectra are also indicative of specific interactions between the polypeptide chain and the porphyrin. The Q0,0 absorption band of alpha-ZnHb at 5 K shows a splitting of approximately 300 cm-1, comparable to the largest split ever reported for a heme protein. This value is approximately 140 cm-1 for beta-ZnHb. The possible origin of the split is discussed in terms of the local electric field imposed by the amino acids of the respective heme pockets, different configurations of the porphyrin, and/or influences of the liganding histidine. The Zn derivatives show quasiline spectra under fluorescence line narrowing conditions, and the resolved excitation spectrum reveals differences in the vibrational levels of the Zn porphyrin in the two subunits. Broad underlying emission in the fluorescence line-narrowed emission spectrum can be accounted for, in part, by the existence of the two closely spaced electronic origins and also by the extent of phonon coupling between the porphyrin and the protein matrix.

The FTIR spectra of the asymmetric carboxylate absorption region of three c -type cytochromes—namely horse heart, yeast and bonito cytochromes c —as well as continuum electrostatic calculations performed on their respective protein matrices, show that these combined methods can target specific protein regions and yield pertinent protein charge information that correlates with the observed spectral data. Deconvolution of the IR carboxylate stretch frequency region (1525–1675 cm −1 ) in the three cytochromes yield different ν oco a distributions. In the case of the bonito cytochrome c carboxylates, two ν oco a populations are clearly distinguishable in the deconvoluted spectra—which is not the case for the more complex ν oco a deconvolutions of the other two cytochromes. The frequency distributions of the calculated potentials are consistent with the experimental observations and we conclude that the IR carboxylate absorption in proteins can be modified by the electrostatic environment.

The influence of the protein on the nonplanarity of the macrocycle for nickel(II)-reconstituted cytochrome c (NiCyt-c) has been investigated with pH-dependent resonance Raman and UV-visible absorption spectroscopy and molecular mechanics calculations. The spectra reveal that NiCyt-c near neutral pH has axially coordinated Ni, but below pH 3 and above pH 12, four-coordinate species predominate. The shape of the structure-sensitive Raman line nu10 of NiCyt-c is asymmetric and broad and it changes with pH. This broad line can be decomposed well into at least two sublines, a low-frequency line that results from a nonplanar conformer and a high-frequency line that arises from a nearly planar conformer. Upon lowering the pH from 3.0 to 1.0, the amount of the nonplanar conformer decreases relative to that of the planar conformer. The decreased nonplanarity can be accounted for in terms of the disruption of a hydrogen-bonding network in the peptide backbone upon lowering the pH. Molecular mechanics (MM) calculations on iron(III) and nickel(II) microperoxidase 5 (MP-5) as well as some model heme derivatives have been carried out in order to locate the part of the protein that causes the heme distortion observed in the X-ray crystal structures of cytochromes c. The energy-optimized structures of MP-5 and the model compounds were analyzed using the normal-coordinate structural decomposition method to specify and quantify the out-of-plane macrocyclic distortions. MM calculations for MP-5 show that two hydrogen bonds formed between the amide groups in the peptide backbone are important in maintaining the ruffled deformation of the macrocycle. All evidence presented supports the hypothesis that the nonplanar distortion of the porphyrin of cytochromes c is largely maintained by a relatively small protein segment including the cysteines, the amino acids between the cysteines, and the adjacent histidine ligand. Hydrogen bonding within the backbone of this segment is important in maintaining the conformation of the peptide that induces the porphyrin distortion.

The carboxylate groups of organic acids give strong absorption in the infrared between approximately 1550 and 1650 cm-1. For acetate and chloroacetate derivatives, the infrared (IR) frequency of the carboxylate antisymmetric stretching mode (v(a)OCO) is related to the square root of the pK of the acid, with a shift of approximately 20 cm-1 to higher frequency for a pK drop in the range 5-3. It follows that v(a)OCO may respond to conditions on the protein surface. In this paper, the IR amide I' and carboxylate absorptions of cytochrome c from horse, yeast, and tuna are compared with model compounds such as Val-Glu and microperoxidase-11, the 11 amino acid fragment of horse cytochrome c containing the covalently bound heme. For microperoxidase-11, the contribution from all four carboxylates can be accounted for and the 1567 cm-1 absorption is assigned to the heme propionates. For the proteins, the carboxylate absorption band is inhomogeneous, i.e., there is a distribution of frequencies. Both the amide I' and carboxylate bands are sensitive to protein conformation as shown by their different pH, salt, and redox dependence.

Perhaps the most important contribution of FLN is that it provides an experimental approach to relate physical changes in the protein to predicted dynamical behavior. It is clear that the sample is inhomogeneously broadened in a continuous manner, consistent with the damped motion of proteins. At the same time configurational substates can be selected, suggesting that there is indeed a hierarchy of protein motion and structure. As yet, identification of the structure, and relating it to the spectra, has not been achieved. It is clear that the electric field exerted by neighboring atoms shifts the electronic transition, and the inhomogeneity is greater when the surrounding disorder is greater. The inhomogeneity for the chromophore in the protein is dependent on the protein conformation and is intermediate between that of a crystal and a glass. The phonon coupling also depends on the chromophore and the protein. Fluorescence line narrowing provides in addition ground- and excited-state vibrational frequencies, thereby allowing for structural differences between the excited-state and the ground-state molecule to be detected.

The alpha or Q0,0 absorption band of horse iron(II) cytochrome c splits and shifts to the blue as temperature decreases over the temperature range of 290-10 K. At room temperature, its maximum is at 18 150 cm-1 and the spectral width is 273 cm-1, whereas at 10 K, the two bands of the Q0,0 transition occur at 18 364 and 18 253 cm-1 and the width of the lowest-energy band is 96 cm-1. Temperature dependent splitting also occurs for zinc cytochrome c, a derivative in which Fe has been replaced by Zn; at 10 K, the peaks in the Q0,0 band region occur at 17 106 and 16 996 cm-1. The peak positions are independent of the cryosolvent (aqueous ethylene glycol or glycerol mixtures). The splitting of the Q0,0 band seen in the protein (approximately 110 cm-1 for iron and zinc cytochrome c) is comparable to the crystal field splitting observed for metalloporphyrins in mixed crystals. In contrast, the Q0,0 band of zinc coproporphyrin III in a glassy solvent (dimethylformamide/ethylene glycol) or in poly(vinyl chloride) shows a blue shift with temperature decrease but no evidence of Q0,0 splitting. Available spectral data show that the Q0,0 band is composed of two nearly degenerate electronic transitions and the split is due to the asymmetry in the heme pocket of the protein that arises from the surrounding polypeptide chain. This asymmetry results in the stabilization of one form of the excited state over the other, according to a Jahn-Teller mechanism.

A simple, rapid method is described for the preparation of cytochrome c peroxidase from baker's yeast. The procedure involves lysis of the yeast in the presence of ethyl acetate, extraction of the peroxidase in 0.05 M sodium acetate buffer, pH 5.0, and the concentration of the crude extract on a DEAE-agarose column. The DEAE eluate is further concentrated by ultrafiltration, and gel filtration of the concentrate results in a highly purified form of the enzyme. Consistent yields with 80% recovery are easily obtained. Protein isolated by this method in the presence or absence of the protease inhibitor, phenylmethylsulfonyl fluoride (PMSF), contains purely high-spin Fe(III) heme as monitored by its resonance Raman spectrum.

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

1. Porphyrin cytochrome c, the iron-free derivative of cytochrome c, has been used extensively as a fluorescent analog of cytochrome c. It appears as though its fluorescence intensity but not its relative quantum yield is affected by pH in the physiological range; an apparent pK of about 6.2 is found suggesting a histidine close to the porphyrin. 2. The fluorescence intensity of the porphyrin cytochrome c in the presence of cytochrome c oxidase is independent of pH; this suggests that the oxidase has the capacity to control the pK of whichever group is responsible for the pH sensitivity of the free porphyrin cytochrome c. The most likely candidate for this pH-sensitive group is histidine-18. The N-3 nitrogen of this residue forms one of the axial ligands to the iron in the intact cytochrome c but it is uncoordinated in the iron-free derivative.

Horseradish peroxidase C is a class III peroxidase whose structure is stabilised by the presence of two endogenous calcium atoms. The effect of calcium depletion on protein conformational dynamics was investigated by performing 20-ns molecular dynamics simulations on explicitely solvated hrpc and Ca-depleted hrpc models. Cross-correlation analysis identified correlated motions that are perturbed in the absence of calcium.The trajectories were also analyzed using the essential dynamics method to describe the conformational space sampled by the respective models. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. The results indicate that the motion of the native species is defined by a few preferred directions identified by the first four eigenvectors. The eigenvectors are significantly sampled, suggesting that, on average, large motions involving different subdomains of the protein occur. The analysis reveals that the calcium depletion affects the most important components of the hrpc motions and modify the overall dynamics in regions where functionally significant residues are located -notably in the heme pocket.View Large Image | View Hi-Res Image | Download PowerPoint Slide