Diffusion of small molecules through the structure of myoglobin. Environmental effects

The effect of the ambient solvent viscosity on the mobility of small molecules within myoglobin was studied by substituting Zn-protoporphyrin (ZnPP) for the native Fe-protoporphyrin and using it as an optical probe in the protein (ZnPPMb). The quenching of the ZnPPMb triplet state by oxygen, by anthraquinonesulfonate, and by methyl viologen was followed by exciting it with a laser flash and measuring its decay rate as a function of quencher concentration. The quenching rate constants were taken to measure the diffusion rate of the quencher within the protein. At room temperature, these constants were determined in aqueous and in 37% and 55% (by weight) glycerol-water solutions by measuring the ZnPPMb-delayed fluorescence at 606 nm. It was found that although the quenching rate constants varied the activation energies in the protein were very similar for the different quenchers. In aqueous solution, Ea = 6.0-7.4 kcal/mol; in 37% glycerol, Ea = 6.8-7.5 kcal/mol; and in 55% glycerol, Ea = 8.5-9.2 kcal/mol. The quenching rate of ZnPPMb by oxygen was also measured between 190K and 293K in 80% glycerol, and its triplet decay in the absence of oxygen was determined down to 120K in 88% glycerol. In all experiments, the quenching rates in the protein were compared to those of Zn-hematoporphyrin in the same solvent. The results are discussed in terms of Northrup and McCammon's gated reaction theory.(ABSTRACT TRUNCATED AT 250 WORDS)     

The effect of tryptophanyl substitution on folding and structure of myoglobin.

Mammalian myoglobins contain two tryptophanyl residues at the invariant positions 7 (A-5) and 14 (A-12) in the N-terminal region (A helix) of the protein molecule. The crucial role of tryptophanyl residues has been investigated by site-directed mutagenesis and molecular dynamics simulation. The apomyoglobin mutants with a double W-->F substitution were found to be not correctly folded and therefore not expressed as holoprotein. The introduction of a tyrosyl residue at position 7, that is, W7YW14F, resulted in the expression of a correctly folded myoglobin. Not correctly folded apomyoglobins were found with the following mutants: W7FW14Y, W7EW14F, W7FW14E, W7KW14F, W7FW14K. Moreover, in all these cases, very low levels of expression were observed. The acid-induced denaturation curves of wild-type and folded mutant W7YW14F, obtained following the fluorescence variation of the extrinsic fluorophore 1-anilino-8-naphthalenesulfonate, revealed that the stability of the native state of mutant apoprotein is decreased, thus indicating that the replacement W-->Y in position 7 is able to restore a correct folding but not the same stability. Molecular dynamics simulation indicated that both tryptophans are involved in forming favorable, specific tertiary interactions in the native apomyoglobin structure. The lack of some of these interactions caused by tryptophanyl replacement affects the overall protein structure and may provide an explanation for the observed stability decrease. In the case of the double W-->F substitution, the simulated structure shows conclusively the domain formed by helices A, G and H to be not correctly folded. This effect is attenuated if at least one of the two residues is conserved or a tyrosyl residue replaces W7.     

The effect of solvent viscosity and temperature on DNA viscoelastic behavior

The effect of solvent viscosity (η_s) and temperature (T) on the shape of the concentration dependence of the principal and total recoils in creep-recovery viscoelastometry experiments has been studied for T4 DNA solutions. The range of DNA concentration (c) was 2 – 40 μg/ml; glycerol, 70–80% v/v, sucrose, 60% v/v; NaCl, 5 mM – 1M; and T, 275 – 323 K. A linear proportionality between recoil and c was obtained at high η_s/T. At low η_s/T, the c-dependence was nonlinear, approaching saturation at higher c. At low c, the slope of both curves was the same. Transition between “linear” and “nonlinear” values occurred over a narrow range of η_s/T (a width of 1–5 K if η_s/T was changed by varying T). (η_s/T)_tr, the midpoint of the transition, was independent of solvent properties other than viscosity. Also, (η_s/T)_tr increased with c. For a given c, η_s/T values above this transitional value yield linear behavior; below this, nonlinear behavior. The ratio of linear to nonlinear recoil values is a linear function of c with K_c, the slope of this dependence, independent of η_s and T. A kinetic model for the observed nonlinearity of recoil with c is presented. It explains the independence of K_c on η_s and T. An attempt has been made to explain the linear–nonlinear transitions by comparison of τ₁ and T_R, the lifetime of the contact points of the polymer network in the de Gennes theory. The nonlinear values are consistent with a pseudogel that exists when τ₁ < T_R. At τ₁ > T_R, the DNA behavior is similar to that in dilute solutions (linear values). Thus, the condition for transition is τ₁ = T_R. However, some unsolved problems remain.     

The effect of solvent viscosity on the rate-determining step of fatty acid synthetase

The overall activity of an animal fatty acid synthetase at the saturation level of substrate concentration decreased when the solvent viscosity, eta, of the reaction mixture was increased with viscogens such as glycerol, sucrose, and polyethylene glycol. The activity of the enzyme changed roughly proportional to eta-P, where p = 1.0 for glycerol, p = 0.66 for sucrose, and p less than 0.6 for polyethylene glycol with different molecular sizes. The thioesterase activity, which catalyzes the final partial reaction in the multifunctional enzyme, was not affected by 5-fold increase of solvent viscosity with sucrose. These results suggested that the rate-determining step of the enzyme other than the thioesterase reaction involves a microscopic transport step, the rate of which is influenced by the solvent viscosity. The microscopic transport step may be related to the transfer of the reaction intermediate from one active site to another or to the motion of a larger part of the enzyme requisite for the catalytic reaction. In the solution containing glycerol, the rate-determining motion was primarily diffusion limited since the inverse of the initial rate was proportional to eta, i.e., p = 1. Since the substrate concentration was at a saturation level in this experiment, the viscosity-dependent step cannot be the encounter between the enzyme and substrates, but must be intramolecular in origin, most probably the reaction catalyzed by beta-ketoacyl synthetase. In solutions containing other viscogens, however, p was less than 1.0, indicating a significant involvement of chemical steps in the rate-determining step as well. Bovine serum albumin, when used as a proteinic viscogen, also decreased the initial rate.(ABSTRACT TRUNCATED AT 250 WORDS)     

A novel and simple interpretation of the three-dimensional structure of globular proteins based on quantum mechanical computations on small model molecules. II. The clusters of myoglobin

We show that the bridges and clusters of the kind found earlier in rubredoxin also occur in myoglobin. The familiar acid–base(salt) bridges are at the heart of many of the clusters. The data reported by Takano [(1977) J. Mol. Biol., 110 , 537–508] on sperm whale metmyoglobin reveals ten clusters together with the heme structure and some unidentified fragments. The data reported by Hanson and Schoenborn [(1981) J. Mol. Biol. 153 , 117–146] on carbonmonoxymyoglobin shows numerous bridges that generally fit into the cluster picture derived from Takano's data. The data reported by S. E. V. Phillips (personal communication) on oxymyoglobin also fits into the same pattern of cluster. All but one of the myoglobin clusters are formed from two pieces of local structure each of which covers up to fifteen consecutive aminoacid residues. A “long-range” union joins the two pieces of local structure. The clusters are generally involved with the nonhelical structure and with the junctions of the two helices. Two applications of the theory are suggested. First, the clusters offer an explanation for the formation of the α-helices and second the unfolding of the myoglobin molecule at pH 3.8 to 5.3 is readily interpreted as the protonation at this pH of the carboxylate ions that generate the clusters that in turn support the α-helices.     

The covalent structure of dog myoglobin.

The primary structure of the myoglobin of the domestic dog (German shepherd) was studied. Tryptic and thermolytic peptides were compared with the sequence of other known myoglobins; the stepwise automatic Edman's degradation of the whole globin and also the chymotryptic digestion of the median fragment obtained by CNBr cleavage completed this sequence. Comparison of the established dog myoglobin structure with those from other carnivora shows 16 differences versus badger, 20 versus harbour seal and 15 versus California sea lion.     

The viscosity of neutrophils and their transit times through small pores

Passive neutrophils from five different individuals are rapidly aspirated at constant suction pressure and at room temperature into a pipet with a diameter of 4 microns. The excess suction pressures (i.e., the pressures in excess of the small threshold pressure required to produce continuous flow into the pipet) are 5000, 10,000 and 20,000 dyn/cm2 (0.5, 1 and 2 kPa) and are comparable to those encountered in the microcirculation. The rate of entry into the pipet is modeled with a linearized version of a theory by Yeung and Evans for the newtonian flow of a neutrophil into a pipet or pore. From this theory and measurements of the cell size and its rate of entry into the pipet, we can calculate a value for the cytoplasmic viscosity. A linear (newtonian) fit of the theory to the experimental data gives a value for the viscosity of 1050 poise. A non-linear fit predicts a decrease in the "apparent viscosity" from about 1500 poise at zero excess pressure to 1000 poise at an excess aspiration pressure of 20,000 dyn/cm2. Our experiments and analysis also allow us to calculate a value for the transit time through short pores over a wide range of excess aspiration pressures and pore diameters. For example, for a pore diameter of 3 microns and an aspiration pressure of 1250 dyn/cm2, we predict a transit time of about 70 s. At 6 microns and 20,000 dyn/cm2, the predicted transit time is only about 0.04 s.     

The effect of regulating molecules on the structure of the PPAR-RXR complex

The PPAR-RXR complex is one of the most significant and prevalent regulatory systems, controlling lipid metabolism by gene expression. Both proteins are members of the nuclear hormone receptor family, consisting of a ligand-binding domain (LBD), a hinge and a DNA binding domain (DBD). The two proteins form a heterodimer in the nucleus. The ligand-free complex interacts with corepressor proteins and blocks the expression of the genes. With the activating ligands and coactivator segments of regulating proteins, the heterodimer becomes active and allows translation of the genes under its control. We implemented model-independent all-atom molecular dynamics simulations for clarifying the structure changes that the activating ligand and the regulatory peptides impose on the PPAR-RXR system, starting with an LBD up to the PPAR-RXR-DNA complex. The simulations were carried out first with an active state of the protein. Once the relaxed state was attained, it was transformed into the inactive-state, the resulting structure was simulated. As the complex alternates between the active-inactive conformations, most of the changes are noticed at the junction area between the two subunits, located on the surface of a long fused helical structure made of H10–H11 of the proteins. The significant differences between the states included enhanced rigidity of the inactive complex, enhancement of tight contacts. The main drive for the transformation is the relocation of the tip of H12 of the PPAR that drives the carboxylate of the C-terminal towards the junction between H10–H11 of the RXR, leading to rearrangement of the main contact zone of the proteins.     

The influence of solvent viscosity on the threshold values of primary tastes

To study the effect of viscosity on the stimulus and recognitionthresholds of the four primary tastes, four thickeners whichare relatively free of any specific taste and odour (methylcellulose, guar seed flour, carob seed flour and tara seed flour)were selected from a greater number of such substances availableand the stimulus and recognition thresholds at the four viscositylevels, 1, 10, 100 and 1 000 cP determined for sucrose, sodiumchloride, caffeine and citric acid. It was found that with increasingviscosity also the stimulus and recognition thresholds increasedwhich means that higher concentrations of the taste substancesare necessary to produce the same taste intensity at higherviscosity. Significant differences were found particularly atthe highest viscosity level of 1 000 cP as compared to 1 cP.By relating the individual values to the corresponding valueat 1 cP it is possible to determine, bymeans of the so-calledamplification factors, the concentration increase that is necessarytoproduce the same taste intensity as at 1 cP. On the basisof regression straight lines, this amplification factor canbe determined for all substances and taste qualities. It becameobvious in this way that not only viscosity is important butalso the kind of thickener used.     

The covalent structure of dog myoglobin

The primary structure of the myoglobin of the domestic dog (German shepherd) was studied. Tryptic and thermolytic peptides were compared with the sequence of other known myoglobins; the stepwise automatic Edman's degradation of the whole globin and also the chymotryptic digestion of the median fragment obtained by CNBr cleavage completed this sequence. Comparison of the established dog myoglobin structure with those from other carnivora shows 16 differences versus badger, 20 versus harbour seal and 15 versus California sea lion.     

Structural dynamics of myoglobin: effect of internal cavities on ligand migration and binding

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation. Photoproduct intermediates, characterized by CO in different locations, were selectively enhanced by laser illumination at specific temperatures. Measurements were performed on the wild-type protein and a series of mutants (L104W, I107W, I28F, and I28W) in which bulky amino acid side chains were introduced to block passageways between cavities or to fill these sites. Binding of xenon was also employed as an alternative means of filling cavities. In all samples, photolyzed CO ligands were observed to initially bind at primary docking site B in the vicinity of the heme iron, from where they migrate to the secondary docking sites, the Xe4 and/or Xe1 cavities. To examine the relevance of these internal docking sites for physiological ligand binding, we have performed room-temperature flash photolysis on the entire set of proteins in the CO- and O(2)-bound form. Together with the cryospectroscopic results, these data provide a clear picture of the role of the internal sites for ligand escape from and binding to myoglobin.     

Modulation of DNA structure formation using small molecules

Genome integrity is essential for proper cell function such that genetic instability can result in cellular dysfunction and disease. Mutations in the human genome are not random, and occur more frequently at “hotspot” regions that often co-localize with sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures. Non-B DNA-forming sequences are mutagenic, can stimulate the formation of DNA double-strand breaks, and are highly enriched at mutation hotspots in human cancer genomes. Thus, small molecules that can modulate the conformations of these structure-forming sequences may prove beneficial in the prevention and/or treatment of genetic diseases. Further, the development of molecular probes to interrogate the roles of non-B DNA structures in modulating DNA function, such as genetic instability in cancer etiology are warranted. Here, we discuss reported non-B DNA stabilizers, destabilizers, and probes, recent assays to identify ligands, and the potential biological applications of these DNA structure-modulating molecules.     

Migration of small molecules through the structure of hemoglobin: evidence for gating in a protein electron-transfer reaction

It has previously been shown that the rates and activation energies for migration molecules of different sizes through myoglobin are very similar. The results were interpreted in terms of conformational changes in the protein structure that facilitate the passage of the different molecules to a similar extent. Here we ask whether the quaternary structural changes that accompany the binding of ligands (O2 or CO) to hemoglobin might influence the migration rate from the solution into the protein's binding site. As a model for the R state of hemoglobin, we used the protein in which the Fe protoporphyrin (FePP) in the alpha subunit was substituted by Zn protoporphyrin (ZnPP) and the oxidized heme was ligated by CN-. The T state of hemoglobin was represented by the protein in which all four FePP groups were substituted by ZnPP. The quenching rate of the excited ZnPP triplet state within the hemoglobin by oxygen, methyl viologen, and anthraquinonesulfonate served as a measure of the migration rate through the protein into the binding site. It was found that the activation energies for all three quenchers were very similar and closely resembled those in myoglobin, suggesting that the migration rates are determined by the subunit structure only and that the quaternary configurational changes do not influence the quenching rates. The implications of the results for electron transfer in proteins are briefly discussed.     

The effect of viscosity on the perception of flavour

A trained sensory panel assessed flavour and sweetness intensity in solutions containing varying concentrations of hydroxy propyl methylcellulose (HPMC), sugar and flavour volatile. The flavour and sweetness of the viscous solutions were rated using magnitude estimation with a controlled modulus. In addition, the concentration of volatile released on the breath was measured using MS Nose. For low concentrations of HPMC (<0.5 g/100 g), perceived flavour intensity remained the same; however, a steady decrease was noted at higher concentrations (>0.6 g/100 g). The change in perceived intensity occurred at the point of random coil overlap (c(*)) for this hydrocolloid. The perceived sweetness of the solution showed a similar pattern with increasing HPMC concentration, although the inflection at c(*) was not so obvious. Despite the change in perceived flavour intensity, the actual concentration of volatile measured on the breath was not affected by the change in HPMC concentration. Low-order polynomial models were produced to describe perceived flavour intensity and sweetness in viscous solutions containing HPMC and potential explanations for the changes in perception are discussed.     

Quantifying the importance of protein conformation on ligand migration in myoglobin

Myoglobin (Mb) is a model system for ligand binding and migration. The energy barriers (ΔG) for ligand migration in Mb have been studied in the past by experiment and theory and significant differences between different approaches were found. From experiment, it is known that Mb can assume a large number of conformational substates. In this work, these substates are investigated as a possible source of the differences in migration barriers. We show that the initial structure significantly affects the calculated ΔG for a particular transition and that fluctuations in barrier heights δΔG are of similar magnitude as the free energy barriers themselves. The sensitivity of ΔG to the initial structure is compared to other sources of errors. Different protein structures can affect the calculated ΔG by up to 4 kcal/mol, whereas differences between simple point charge models and more elaborate multipolar charge models for the CO-ligand are smaller by a factor of two. Analysis of the structural changes underlying the large effect of the conformational substate reveals the importance of coupling between protein and ligand motion for migration.