Molecular dynamics simulations of human glutathione transferase P1-1: conformational fluctuations of the apo-structure.

We have investigated by molecular dynamics simulations the conformational fluctuations of the monomer of human apo-glutathione transferase P1-1. After attainment of steady-state dynamics, the structural fluctuations involve mainly the protein segments that participate also in the holo-apo transition discussed in the accompanying article (Stella et al., 1999:37:1-9.). The most mobile region is the C-terminal segment of helix 2. In contrast, helices 1, 6, 7, and 8 constitute a relatively rigid protein core. An "essential dynamics" analysis of the simulation shows that the largest fluctuations involve specific regions of glutathione transferases. In such regions, atomic motions are correlated. Motions of helix 2 are accounted for by the second most prominent principal component, which reveals a fluctuation between two distinct conformations. The residues that constitute the H-site undergo a breathing motion, possibly relevant during the binding of hydrophobic cosubstrates. Based on our simulation, several experimental findings can be rationalized, including the viscosity-dependent reactivity of Cys 47 and Cys 101 as well as the selective proteolysis of the peptide bond between Lys 44 and Ala 45. We have also modeled the structural changes that lead to the formation of an intrachain disulfide bridge between cysteines 47 and 101 and to the inactivation of the enzyme. The resulting structure maintains essentially the native fold except for helix 2, which closes the G-site. Proteins 1999;37:10-19.

Molecular dynamics simulations of human glutathione transferase P1-1: analysis of the induced-fit mechanism by GSH binding.

We report here a 1-ns molecular dynamics simulation on the ligand-free monomer of human glutathione transferase P1-1 in bulk water. The average conformation obtained from the last 500 ps of simulation is taken as a model for the apo-structure of this protein and compared to the available crystallographic data. Remarkable changes in the tertiary structure take place during the simulation and are ascribed to the removal of the ligand. They support an induced fit mechanism occurring upon glutathione binding, whose major features can be described in detail. A portion of helix 2 (residues 42-50), which participates in the formation of the active site, undergoes the most prominent conformational changes. Other protein segments, such as the C-terminal loop and helix 4, also show relevant structural rearrangements. All these transitions cause a significant shielding from the solvent of the hydrophobic binding site of the co-substrate, whose exposed surface goes from 4.6 nm(2) in the holo-structure to about 3.1 nm(2) in the apo-conformation. The results of this simulation are consistent with numerous experimental observations previously obtained on GST P1-1 and provide new insights for their explanation at the molecular level. Proteins 1999;37:1-9.

Dynamic properties of monomeric insect erythrocruorin III from Chironomus thummi-thummi: relationships between structural flexibility and functional complexity.

We have investigated the kinetics of geminate carbon monoxide binding to the monomeric component III of Chironomus thummi-thummi erythrocruorin, a protein that undergoes pH-induced conformational changes linked to a pronounced Bohr effect. Measurements were performed from cryogenic temperatures to room temperature in 75% glycerol and either 0.1 M potassium phosphate (pH 7) or 0.1 potassium borate (pH 9) after nanosecond laser photolysis. The distributions of the low temperature activation enthalpy g(H) for geminate ligand binding derived from the kinetic traces are quite narrow and are influenced by temperature both below and above approximately 170 K, the glass transition temperature. The thermal evolution of the CO binding kinetics between approximately 50 K and approximately 170 K indicates the presence of some degree of structural relaxation, even in this temperature range. Above approximately 220 K the width of the g(H) progressively decreases, and at 280 K geminate CO binding becomes exponential in time. Based on a comparison with analogous investigations of the homodimeric hemoglobin from Scapharca inaequivalvis, we propose a link between dynamic properties and functional complexity.

Crystal structure of Asian elephant (Elephas maximus) cyano-metmyoglobin at 1.78-A resolution. Phe29(B10) accounts for its unusual ligand binding properties.

The crystal structure of Asian elephant cyano-metmyoglobin which has a glutamine instead of the usual distal site histidine has been determined to high resolution. In addition to this replacement, the substitution of a conserved leucine residue in position 29(B10) at the distal side by a phenylalanine was unambiguously identified based on the available electron density. The suspicion, that there were errors in the original sequence which has caused some confusion, is thus confirmed. Comparison with other myoglobin structures in various ligated forms reveals an essentially unchanged tertiary structure in elephant myoglobin despite the two amino acid substitutions in the heme pocket. Our current structural model shows that the N epsilon 2 atom of Gln64(E7) has moved with respect to the corresponding nitrogen position of His64(E7) in the CO complex of sperm whale myoglobin. The newly assigned residue Phe29(B10) penetrates into the distal side of the heme pocket approaching the ligand within van der Waals distance and causing a much more crowded heme pocket compared to other myoglobins. Kinetic properties of Asian elephant myoglobin, wild type, and recombinant sperm whale myoglobins are discussed in relation to the structural consequences of the two amino acid substitutions H64Q and L29F.

Stereodynamic properties of the cooperative homodimeric Scapharca inaequivalvis hemoglobin studied through optical absorption spectroscopy and ligand rebinding kinetics.

The study of the thermal evolution of the Soret band in heme proteins has proved to be a useful tool to understand their stereodynamic properties; moreover, it enables one to relate protein matrix fluctuations and functional behavior when carried out in combination with kinetic experiments on carbon monoxide rebinding after flash photolysis. In this work, we report the thermal evolution of the Soret band of deoxy, carbonmonoxy, and nitric oxide derivatives of the cooperative homodimeric Scapharca inaequivalvis hemoglobin in the temperature range 10-300 K and the carbon monoxide rebinding kinetics after flash photolysis in the temperature range 60-200 K. The two sets of results indicate that Scapharca hemoglobin has a very rigid protein structure compared with other hemeproteins. This feature is brought out i) by the absence of nonharmonic contributions to the soft modes coupled to the Soret band in the liganded derivatives, and ii) by the almost "in plane" position of the iron atom in the photoproduct obtained approximately 10(-8) s after dissociating the bound carbon monoxide molecule at 15 K.

Structure-dynamics-function relationships in Asian elephant (Elephas maximus) myoglobin. An optical spectroscopy and flash photolysis study on functionally important motions.

In this work we report the thermal behavior (10-300 K) of the Soret band lineshape of deoxy and carbonmonoxy derivatives of Asian elephant (Elephas maximus) and horse myoglobins together with their carbon monoxide recombination kinetics after flash photolysis; the results are compared to analogous data relative to sperm whale myoglobin. The Soret band profile is modeled as a Voigt function that accounts for the coupling with high and low frequency vibrational modes, while inhomogeneous broadening is taken into account with suitable distributions of purely electronic transition frequencies. This analysis makes it possible to isolate the various contributions to the overall lineshape that; in turn, give information on structural and dynamic properties of the systems studied. The optical spectroscopy data point out sizable differences between elephant myoglobin on one hand and horse and sperm whale myoglobins on the other. These differences, more pronounced in deoxy derivatives, involve both the structure and dynamics of the heme pocket; in particular, elephant myoglobin appears to be characterized by larger anharmonic contributions to soft modes than the other two proteins. Flash photolysis data are analyzed as sums of kinetic processes with temperature-dependent fractional amplitudes, characterized by discrete pre-exponentials and either discrete or distributed activation enthalpies. In the whole temperature range investigated the behavior of elephant myoglobin appears to be more complex than that of horse and sperm whale myoglobins, which is in agreement with the increased anharmonic contributions to soft modes found in the former protein. Thus, to satisfactorily fit the time courses for CO recombination to elephant myoglobin five distinct processes are needed, only one of which is populated over the whole temperature range investigated. The remarkable convergence and complementarity between optical spectroscopy and flash photolysis data confirms the utility of combining these two experimental techniques in order to gain new and deeper insights into the functional relevance of protein fluctuations.

Protein dynamics in minimyoglobin: is the central core of myoglobin the conformational domain?

The kinetics of CO binding to the horse myoglobin fragment Mb-(32-139), the so-called "mini-Mb," were investigated by laser flash photolysis in 0.1 M phosphate buffer and in buffer with 75% (vol/vol) glycerol. The reaction displays complex time courses that can be approximated satisfactorily only with a sum of five exponentials. The features of the kinetic components and a comparison of the deoxy-minus-carbonyl difference spectra of mini-Mb and horse Mb obtained under equilibrium conditions, with the kinetic difference spectra resulting from the global analysis of the traces recorded between 400 and 450 nm, show that CO binding to mini-Mb is accompanied by large structural changes. In view of the fact that mini-Mb is an approximation of the Mb-(31-105) fragment encoded by the central exon of the Mb gene, this finding is particularly relevant. On the basis of our data and previous reports [De Sanctis, G., Falcioni, G., Giardina, B., Ascoli, F. & Brunori, M. (1988) J. Mol. Biol. 200, 725-733; De Sanctis, G., Falcioni, G., Grelloni, F., Desideri, A., Polizo, F., Giardina, B., Ascoli, F. & Brunori, M. (1992) J. Mol. Biol. 222, 637-643], we propose that the protein fragment encoded by the central exon of the Mb gene is the domain responsible for ligand-linked conformational transitions, while the two terminal fragments dampen the amplitude of the structural changes that accompany ligand binding, thus rendering the protein stable and kinetically more efficient in its physiological function.

Protein dynamics. An overview on flash-photolysis over broad temperature ranges.

Ligand binding kinetics to heme-proteins between 40 and 300 K point to a regulatory role of protein dynamics. A protein-specific susceptibility of the heme-iron reactivity to dynamic fluctuations emerges from the distribution of reaction enthalpies derived from flash-photolysis measurements below ca. 180 K; we quantify it in terms of 'intramolecular viscosity', postulating that narrow low-temperature enthalpy distributions correspond to low internal viscosity and vice versa. The thermal evolution of ligand binding kinetics suggests, with other results, an interplay between high-frequency transitions of the amino acid side chains and low-frequency collective motions as a possible regulatory mechanism of protein dynamics.

Isoelectric focusing of cytochrome P450: isolation of six phenobarbital-inducible rat liver microsomal isoenzymes.

A procedure for the isolation of native proteins from membranes by isoelectric focusing is described. It was used to resolve into six components the major fraction of cytochrome P450, obtained from liver microsomes of phenobarbital-treated rats, after chromatography on DE-52 cellulose. When eluted from the gel, these proteins are in a native form as shown by (a) the light absorption spectra of the Soret region of their reduced carbonyl derivatives, all characterized by maxima around 450 nm, and (b) their enzymatic activities toward three different substrates. Characterization by a monoclonal antibody and partial sequence analysis of tryptic peptides reveal that three of the IEF-purified proteins have P450IIB1 character, whereas the other three are related to P450IIB2.

Protein dynamics. Comparative investigation on heme-proteins with different physiological roles.

We report the low temperature carbon monoxide recombination kinetics after photolysis and the temperature dependence of the visible absorption spectra of the isolated alpha SH-CO and beta SH-CO subunits from human hemoglobin A in ethylene glycol/water and in glycerol/water mixtures. Kinetic measurements on sperm whale (Physeter catodon) myoglobin and previously published optical spectroscopy data on the latter protein and on human hemoglobin A, in both solvents, (Cordone, L., A. Cupane, M. Leone, E. Vitrano, and D. Bulone. 1988. J. Mol. Biol. 199:312-218) are taken as reference. Low temperature flash photolysis data are analyzed within the multiple substates model proposed by Frauenfelder and co-workers (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355-5373). Within this model a distribution of activation enthalpies for ligand binding accounts for the structural heterogeneity of the protein, while the preexponential factor, containing also the entropic contribution to the free energy of the process, is considered to be constant for all conformational substates. Optical spectra are deconvoluted in gaussian components and the temperature dependence of the moments of the resulting bands is analyzed, within the harmonic Frank-Condon approximation, to obtain information on the stereodynamic properties of the heme pocket. The kinetic and spectral parameters thus obtained are found to be protein dependent also with respect to their sensitivity to changes in the composition of the external medium. A close correlation between the kinetic and spectral features is observed for the proteins examined under all experimental conditions studied. The results reported are discussed in terms of differences in the heme pocket structure and in the conformational heterogeneity among the various proteins, as related to their different capability to accommodate constraints imposed by the external medium.

Isocyanide binding kinetics to monomeric hemoproteins. A study on the ligand partition between solvent and heme pocket.

The kinetics of methyl-, ethyl-, iso-propyl-, and ter-butyl-isocyanide binding to Aplysia limacina myoglobin (distal His----Lys) and the isolated beta chains from hemoglobin Zurich (distal His----Arg) have been investigated by flash photolysis at various temperatures above 0 degrees C. Sperm whale (Physter catodon) myoglobin and the isolated beta chains from normal adult hemoglobin have been used as references. In most reaction systems investigated the apparent extent of photolysis increases with temperature. For sperm whale myoglobin and the normal beta chains the increase is of the same magnitude and not correlated to the type of ligand used. On the contrary, for the two proteins lacking the distal histidine, the phenomenon is dependent on the size of the alkyl side chain of the ligand. The results, analyzed on the basis of the multibarrier model (Austin, R.H., K.W. Beeson, L. Eisenstein, H. Frauenfelder, and I.C. Gunsalus, 1975, Biochemistry, 16:5355-5373), suggest that the partition of the ligand molecules between the solvent and the heme pocket, occurring during the photolysis process, is primarily determined by interactions between the ligand and residues in the heme cavity rather than by diffusion through the protein matrix.

Deformability of the hemoglobin molecule as the basis of its functional behavior.

It has been proposed that the process of ligand binding to hemoproteins is governed by sequential potential barriers as depicted in the following scheme [Austin et al., Biochemistry 14, pp. 5355-5373, 1975; Ansari et al., Proc. natn. Acad. Sci. USA. pp. 5000-5004, 1985]: (formula; see text) At sufficiently high temperatures (approximately 300 K), the binding process is exponential in time and its velocity proportional to the concentration of the ligand in the solvent. At low temperatures (typically less than 210 K), the combination velocity is not correlated with the amount of free ligand in the solvent and the time course is nonexponential (approximately power law). This implies that (a) proteins are flexible structures and fluctuate between many substates, each having a different ligand binding velocity, and (b) that the rate of interconversion between the various conformational substates is high compared to that of ligand binding at and above room temperature and extremely slow below approximately 200 K. We have investigated the temperature dependence of CO binding after a photolysing pulse to a variety of myoglobins and hemoglobin subunits. The results show that for all myoglobins investigated, the most probable activation energy for the inner process (B----A) is in the order of 10 kJ/mol, while for the hemoglobin subunits it varies between 2.5 and 5 kJ/mol. The distribution of barrier heights is very broad in myoglobins, somewhat narrower in the non-alpha-chains and narrowest in the alpha-hemoglobin subunits. This finding implies that different degrees of flexibility exist between storage and transport proteins and that hemoglobin cooperativity may be related to the presence of subunits with different flexibility within the tetramer.

Kinetics of carbon monoxide binding to phenobarbital-induced cytochrome P-450 from rat liver microsomes: a simple bimolecular process.

The kinetics of carbon monoxide binding to phenobarbital-induced cytochrome P-450 (P-450PB) and to its enzymatically inactive form P-420PB have been investigated by both stopped-flow and flash-photolysis spectrophotometry. When the simultaneous presence of both forms of the enzyme is taken into account, the binding of CO to these two proteins can be described in terms of two bimolecular processes with rate constants of 4.5 X 10(6) M-1.S-1 and 4.7 X 10(5) M-1.S-1 for P-450PB and 1.7 X 10(7) M-1.S-1 and 1.5 X 10(6) M-1.S-1 for P-420PB. From kinetic studies of the binding of CO to P-450PB under different experimental conditions, investigations of the homogeneity of our P-450PB preparations, and comparative kinetic investigations of P-450s from different sources, we conclude that CO binding to reduced P-450PB is a simple bimolecular process and that the observed biphasic traces are due to heterogeneity of the proteins. This conclusion is in contrast with previous reports of complex reaction mechanisms for the binding of CO to P-450PB. Optical spectroscopy studies indicate the existence of an equilibrium between P-450PB and P-420PB, at least for the reduced carbonyl derivatives of the enzymes. The interconversion is strongly influenced by the aggregation state of the protein. Large differences between the CO binding properties of P-450PB and those of P-420PB are found. These are discussed in terms of possible effects of the proximal ligation state of the iron on heme reactivity.

Encapsulation of hemoglobin in phospholipid liposomes: characterization and stability.

Hemoglobin is encapsulated in liposomes of different lipid composition. The resulting dispersion consists primarily of multilamellar liposomes (hemosomes) of a wide particle size distribution (diameter ranging mainly between 0.1 and 1 micron). The encapsulation efficiency is significantly larger with liposomes containing negatively charged lipids as compared to liposomes made of phosphatidylcholine. The integrity of the phospholipid bilayer is maintained in the presence of hemoglobin. The reaction rate of CO binding to encapsulated hemoglobin is reduced compared to that of free hemoglobin, but it is still greater than that observed in red blood cells. Hemoglobin encapsulated in liposomes made from negatively charged phospholipids is less stable than hemoglobin entrapped in isoelectric phosphatidylcholine. The instability of hemoglobin is due to the protein interacting with the negatively charged lipid bilayer. This interaction leads in turn to hemoglobin denaturation, possibly involving the dissociation of the heme group from the heme-globin complex. The nature of the negatively charged phospholipid is important in promoting the interaction with hemoglobin, the effect being in the order phosphatidic acid greater than phosphatidylinositol congruent to phosphatidylglycerol greater than phosphatidylserine. The presence of equimolar amounts of cholesterol in the phospholipid bilayer has a stabilizing effect on hemoglobin. This effect is pronounced with saturated phospholipids, but it is also observed, though to a lesser extent, with unsaturated ones, indicating that the bilayer fluidity has a modulating effect. The presence of cholesterol possibly interferes with secondary interactions following the binding of hemoglobin to the negatively charged lipid bilayer.

Kinetics of oxygen and carbon monoxide binding to liver fluke (Dicrocoelium dendriticum) hemoglobin. An extreme case?

The kinetics of oxygen and carbon monoxide binding to the monomeric liver fluke (Dicrocoelium dendriticum) hemoglobin have been studied. The ligand association rates are approximately 1 X 10(8) and approximately 3 X 10(8) M-1 s-1, respectively, for CO and O2 and show no pH dependence. On the contrary the ligand dissociation rates decrease by lowering the pH below 7, the pK of the transition being around 5.5. These findings, together with spectroscopic properties of the protein, are discussed in relation to the fact that, in this hemoglobin, the distal histidine is replaced by a glycine.

Studies on the oxidation of hemoglobin Zurich (beta 63 E7 Arg).

Autoxidation and chemically-induced oxidation of hemoglobin Zurich (beta 63 E7 Arg) have been investigated by electron paramagnetic resonance and optical absorption spectroscopy. The results show that the replacement of the distal histidine of the hemoglobin beta chains by an arginine greatly enhances the susceptibility of the heme-iron to oxidative challenge. Both the kinetics and the products of the oxidation are pH dependent. Thus, at acidic and neutral pH, treatment of the protein with ferricyanide leads to a fast conversion of the oxy-protein to aquo-methemoglobin, which, eventually, is slowly converted to hemichromes. In contrast, the hydroxy-met derivative, formed upon chemical oxidation at high pH, is rapidly converted to hemichromes. The electron paramagnetic resonance features of the ferric derivatives of hemoglobin Zurich are somewhat singular, reflecting the modifications of the heme environment in the distal region of the abnormal chains. However, they can be related to heme complexes having their structural counterparts in oxidation products of hemoglobin A.

Electron paramagnetic resonance properties of liver fluke (Dicrocoelium dendriticum) nitrosyl hemoglobin.

The electron paramagnetic resonance properties of the nitric oxide derivative of liver fluke (Dicrocoelium dendriticum) hemoglobin (DD-Hb) have been investigated in the pH range from 4.8 to 7.8. In the neutral and alkaline regions the spectra have a rhombic shape, with gx = 2.09, gy = 1.99 and gz = 2.009, and a triplet hyperfine structure of 2.2 mT, due to the nitrogen of the bound NO molecule, in the center resonance. No superhyperfine lines in the gz region, related to the interaction of the iron with the proximal histidine, are detected, suggesting a large distance between the metal and the N epsilon of the imidazole. By lowering the pH the EPR spectrum undergoes a reversible change showing a 3-line pattern in the high-field region. Such a spectrum is fully formed at pH 4.8 and is interpreted in terms of a dissociation of the proximal histidine from the heme iron.

Kinetic properties of cobalt--iron hybrid hemoglobins.

The replacement of O2 with CO was studied on cobalt-iron hemoglobin hybrids. Both proto- and mesocobalt hemes were used for the reconstitution. In the oxy quaternary conformation no difference is observed between alpha- and beta-subunits when only proto hemes are present in the hybrid (k4 = 30 s-1, k'4/l'4 = 2.5). If Co-meso heme is present on the beta-chains the binding properties of the partner subunit are modified (k'4/l'4 = 4).

Erythrocytic proteases: preferential degradation of alpha hemoglobin chains.

Proteolytic activity against native hemoglobin polypeptide chains is demonstrated, under strictly physiological conditions, in human reticulocytes of both normal subjects and individuals suffering from a variety of pathologic conditions involving erythrocytes, including beta-thalassemia. Two thirds of the activity are found in the cytoplasm and the remainder of it is associated with the reticulocyte membrane. That this proteolytic activity is due to contamination by WBC is excluded. The activity preferentially degrades the alpha-hemoglobin chains. An increase in this substrate within the erythroid cells, as observed in beta-thalassemia, does not enhance proteolysis. Protease inhibitors produce a variable decrease in proteolysis. None inhibit completely, thus showing that several enzymes, with different specificities, are involved.