Photoselection in polarized photolysis experiments on heme proteins.

Polarized photolysis experiments have been performed on the carbon monoxide complex of myoglobin to assess the effects of photoselection on the kinetics of ligand rebinding and to investigate the reorientational dynamics of the heme plane. The results are analyzed in terms of the optical theory developed in the preceding paper by Ansari and Szabo. Changes in optical density arising from rotational diffusion of the photoselected population produce large deviations from the true geminate ligand rebinding curves if measurements are made with only a single polarization. The apparent ligand rebinding curves are significantly distorted even at photolysis levels greater than 90%. These deviations are eliminated by obtaining isotropically-averaged optical densities from measurements using both parallel and perpendicular polarizations of the probe pulse. These experiments also yield the optical anisotropy, which gives a novel method for accurately determining the degree of photolysis, as well as important information on the reorientational dynamics of the heme plane. The correlation time for the overall rotational diffusion of the molecule is obtained from the decay of the anisotropy. The anisotropy prior to rotational diffusion is lower than that predicted for a rigidly attached, perfectly circular absorber, corresponding to an apparent order parameter of S = 0.95 +/- 0.02. Polarized absorption data on single crystals suggest that the decreased anisotropy results more from internal motions of the heme plane which take place on time scales shorter than the duration of the laser pulse (10 ns) than from out-of-plane polarized transitions.     

Long-timescale simulation methods

The outstanding challenges in computer simulations of biological macromolecules are related to their complexity. Part of the complexity of biological systems concerns their physical size. Enumerating atoms ranging from a few in small signal molecules to the millions of particles in biological complexes is an obvious example of biological hierarchy. Another aspect is the extremely broad range of timescales of life science processes (many orders of magnitude); this adds another dimension of complexity. This extended range of timescales may even be observed for a single biomolecular process. Consider, for example, the R to T transition in hemoglobin. The complete conformational change occurs in tens of microseconds. However, the system has more than one timescale. Considerable activity occurs on a range of timescales before the final event (heme relaxation, picoseconds; tertiary relaxation, nanoseconds; ligand escape from the protein matrix and rebinding, hundreds of nanoseconds and so on). Whereas the basic time-step of atomically detailed simulations is about a femtosecond, it is not difficult to find molecular processes in biology that span more than ten orders of magnitude of relevant times, making the straightforward simulation of these events very difficult. Several techniques have been developed in recent years to address these problems.     

Reorientational contact-weighted elastic network model for the prediction of protein dynamics: comparison with NMR relaxation

A new model for the prediction of protein backbone motions is presented. The model, termed reorientational contact-weighted elastic network model, is based on a multidimensional reorientational harmonic potential of the backbone amide bond vector orientations and it is applied to the interpretation of dynamics parameters obtained from NMR relaxation data. The individual energy terms are weighted as a function of the intervector distances and by the contact strengths of each bond vector with respect to its local environment. Correlated reorientational motional properties of the bond vectors are obtained by means of normal mode analysis. Application to a set of proteins with known three-dimensional structures yields good to excellent agreement between predicted and experimental NMR order parameters presenting an improvement over the local contact model. The reorientational eigenmodes of the reorientational contact-weighted elastic network model method provide direct information on the collective nature of protein backbone motions. The dominant eigenmodes have a notably low collectivity, which is consistent with the behavior found for reorientational eigenmodes from molecular dynamics simulations.     

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used ¹⁵N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the ¹⁵N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.     

Picosecond structural dynamics of myoglobin following photodissociation of carbon monoxide as revealed by ultraviolet time-resolved resonance Raman spectroscopy

Picosecond protein dynamics of myoglobin in response to structural changes in heme upon CO dissociation were observed in a site-specific fashion for the first time using time-resolved UV resonance Raman spectroscopy. Transient UV resonance Raman spectra showed several phases of intensity changes in both tryptophan and tyrosine Raman bands. Five picoseconds after dissociation, the W18, W16, and W3 bands of tryptophan residues and the Y8a band of tyrosine residues decreased in intensity, followed by recovery of the Y8a band intensity in hundreds of picoseconds and recovery of the tryptophan bands in nanoseconds. These spectral changes suggest that the change in heme structure impulsively drives concerted movement of the EF helical section and that rearrangements toward a deoxy structure occur in the heme vicinity and in the A helix within a time frame of sub-nanoseconds to nanoseconds.     

Engineering out motion: a surface disulfide bond alters the mobility of tryptophan 22 in cytochrome b5 as probed by time-resolved fluorescence and 1H NMR experiments

In the accompanying paper [Storch et al. (1999) Biochemistry 38, 5054-5064] equilibrium denaturation studies and molecular dynamics (MD) simulations were used to investigate localized dynamics on the surface of cytochrome b5 (cyt b5) that result in the formation of a cleft. In those studies, an S18C:R47C disulfide mutant was engineered to inhibit cleft mobility. Temperature- and urea-induced denaturation studies revealed significant differences in Trp 22 fluorescence between the wild-type and mutant proteins. On the basis of the results, it was proposed that wild type populates a conformational ensemble that is unavailable to the disulfide mutant and is mediated by cleft mobility. As a result, the solvent accessibility of Trp 22 is decreased in S18C:R47C, suggesting that the local environment of this residue is less mobile due to the constraining effects of the disulfide on cleft dynamics. To further probe the structural effects on the local environment of Trp 22 caused by inhibition of cleft formation, we report here the results of steady-state and time-resolved fluorescence quenching, differential phase/modulation fluorescence anisotropy, and 1H NMR studies. In Trp fluorescence experiments, the Stern-Volmer quenching constant increases in wild type versus the oxidized disulfide mutant with increasing temperature. At 50 degrees C, KSV is nearly 1.5-fold greater in wild type compared to the oxidized disulfide mutant. In the reduced disulfide mutant, KSV was the same as wild type. The bimolecular collisional quenching constant, kq, for acrylamide quenching of Trp 22 increases 2.7-fold for wild type and only 1.8-fold for S18C:R47C, upon increasing the temperature from 25 to 50 degrees C. The time-resolved anisotropy decay at 25 degrees C was fit to a double-exponential decay for both the wild type and S18C:R47C. Both proteins exhibited a minor contribution from a low-amplitude fast decay, consistent with local motion of Trp 22. This component was more prevalent in the wild type, and the fractional contribution increased significantly upon raising the temperature. The fast rotational component of the S18C:R47C mutant was less sensitive to increasing temperature. A comparison of the 1H NMR monitored temperature titration of the delta-methyl protons of Ile 76 for wild type and oxidized disulfide mutant, S18C:R47C, showed a significantly smaller downfield shift for the mutant protein, suggesting that Trp 22 in the mutant protein experiences comparatively decreased cleft dynamics in core 2 at higher temperatures. Furthermore, comparison of the delta'-methyl protons of Leu 25 in the two proteins revealed a difference in the ratio of the equilibrium heme conformers of 1.2:1 for S18C:R47C versus 1.5:1 for wild type at 40 degrees C. The difference in equilibrium heme orientations between wild type and S18C:R47C suggests that the disulfide bond affects heme binding within core 1, possibly through damped cleft fluctuations. Taken together, the NMR and fluorescence studies support the proposal that an engineered disulfide bond inhibits the formation of a dynamic cleft on the surface of cyt b5.     

Molecular dynamics simulations of the orientation properties of cytochrome <Emphasis Type="Italic">c</Emphasis> on the surface of single-walled carbon nanotubes

Electron transfer between cytochrome c (Cytc) and electrodes can be influenced greatly by the orientation of protein on the surface of the electrodes. In the present study, different initial orientations of Cytc on the surface of five types of single-walled carbon nanotubes (SWNTs), with different diameters and chirality, were constructed. Properties of the orientations of proteins on the surface of these tubes were first investigated through molecular dynamics simulations. It was shown that variations in SWNT diameter do not significantly affect the orientation; however, the chirality of the SWNTs is crucial to the orientation of the heme embedded in Cytc, and the orientation of the protein can consequently be influenced by the heme orientation. A new electron pathway between Cytc and SWNT, which hopefully benefits electron transfer efficiency, has also been proposed. This study promises to provide theoretical guidance for the rational design of bio-sensors or bio-fuel cells by using Cytc-decorated carbon nanotube electrodes.     

Role of proximal His93 in nitric oxide binding to metmyoglobin. Application of continuum solvation in Monte Carlo protein simulations.

Monte Carlo protein simulations with continuum solvation were used to explore the conformational mobility of NO within the active site of metmyoglobin. To the best of our knowledge this is the first application of a continuum solvation model for exploring protein binding sites. The usefulness of the Monte Carlo conformational analysis was demonstrated in comparative molecular dynamics simulations. Analysis of conformer populations revealed that Monte Carlo conformational analysis is more effective in mapping the relevant region of the potential surface. Distribution of low-energy conformers obtained for the metmyoglobin-NO complex was found to depend on the orientation of proximal His93. Different charge distributions corresponding to the two experimentally verified possible torsions of this proximal residue result in strong binding of NO or its release to a nearby hydrophobic trap. Conformer populations obtained by Monte Carlo conformational analysis were grouped into three main families: one, with the NO directly bound to the iron, appears when the CA-CB-CG-CD2 torsion of His93 is at its ligand binding value (-113 degrees); and two conformers exist where NO is trapped in a nearby hydrophobic pocket, the same cavity that was determined to be the geminate trap of CO in ferrous Mb as a result of the torsional flip of His93 to its ligand releasing state (-125 degrees). Based on this analysis, we suggest that the electrostatic rearrangement coupled to the conformational fluctuation of the proximal His leads to the geminate trapping of the ligand. Conformational rearrangement of the proximal side would provide the possibility of rebinding of the ligand to Fe.     

On the convergence of the conformational coordinates basis set obtained by the essential dynamics analysis of proteins' molecular dynamics simulations.

In this article we present a quantitative evaluation of the convergence of the conformational coordinates of proteins, obtained by the Essential Dynamics method. Using a detailed analysis of long molecular dynamics trajectories in combination with a statistical assessment of the significance of the measured convergence, we obtained that simulations of a few hundreds of picoseconds are in general sufficient to provide a stable and statistically reliable definition of the essential and near constraints subspaces, at least within the nanoseconds time range. Proteins 1999;36:419-424.     

Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

Globins are respiratory proteins that reversibly bind dioxygen and other small ligands at the iron of a heme prosthetic group. Hemoglobin and myoglobin are the most prominent members of this protein family. Unexpectedly a few years ago a new member was discovered and called neuroglobin (Ngb), being predominantly expressed in the brain. Ngb is a single polypeptide of 151 amino acids and despite the small sequence similarity with other globins, it displays the typical globin fold. Oxygen, nitric oxide, or carbon monoxide can displace the distal histidine which, in ferrous Ngb as well as in ferric Ngb, is bound to the iron, yielding a reversible adduct. Recent crystallographic data on carboxy Ngb show that binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities; in particular, the huge internal tunnel that connects the bulk with the active site, peculiar to Ngb, is heavily reorganized. We report the results of extended (90 ns) molecular dynamics simulations in water of ferrous deoxy and carboxy murine neuroglobin, which are both coordinated on the distal site, in the latter case by CO and in the former one by the distal His(64)(E7). The long timescale of the simulations allowed us to characterize the equilibrated protein dynamics and to compare protein structure and dynamical behavior coupled to the binding of an exogenous ligand. We have characterized the heme sliding motion, the topological reorganization of the internal cavities, the dynamics of the distal histidine, and particularly the conformational change of the CD loop, whose flexibility depends ligand binding.     

Ligand preference and orientation in b- and c-type heme-binding proteins

Hemes are often incorporated into designed proteins. The importance of the heme ligand type and its orientation is still a matter of debate. Here, heme ligands and ligand orientation were investigated using a nonredundant (87 structures) and a redundant (1503 structures) set of structures to compare and contrast design features of natural b- and c-type heme-binding proteins. Histidine is the most common ligand. Marked differences in ligation motifs between b- and c-type hemes are higher occurrence of His-Met in c-type heme binding motifs (16.4% vs. 1.4%) and higher occurrence of exchangeable, small molecules in b-type heme binding motifs (67.6% vs. 9.9%). Histidine ligands that are part of the c-type CXXCH heme-binding motif show a distinct asymmetric distribution of orientation. They tend to point between either the heme propionates or between the NA and NB heme nitrogens. Molecular mechanics calculations show that this asymmetry is due to the bonded constraints of the covalent attachment between the heme and the protein. In contrast, the orientations of b-type hemes histidine ligands are found evenly distributed with no preference. Observed histidine heme ligand orientations show no dominating influence of electrostatic interactions between the heme propionates and the ligands. Furthermore, ligands in bis-His hemes are found more frequently perpendicular rather than parallel to each other. These correlations support energetic constraints on ligands that can be used in designing proteins.     

Solution 1H NMR study of the active site structure for the double mutant H64Q/V68F cyanide complex from mouse neuroglobin

Solution (1)H NMR spectroscopy has been carried out to investigate the molecular and electronic structures of the active site in H64Q/V68F double mutant mouse neuroglobin in the cyanomet form. Two heme orientations resulting from a 180 degrees rotation about the alpha-gamma-meso axis were observed with a population ratio about 1:1, and the clearly distinguished B isomer was used to perform the study. Based on the analysis of the dipolar shifts and paramagnetic relaxation constants, the distal Gln(64)(E7) side chain is obtained to adopt an orientation that may produce hydrogen bond between the N(epsilon)H(1) and the Fe-bound cyanide. The side chain of Phe(68)(E11) is oriented out of the heme pocket just like that in triple mutant of cyanide complex of sperm whale myoglobin. A 15 degrees rotation of the imidazole ring in axial His(96) is observed, which is close to the varphi angle determined from the crystal structure of NgbCO. The quantitative determinations of the orientation and anisotropies of the paramagnetic susceptibility tensor reveal that cyanide is tilted by 8 degrees from the heme normal which allows for contact to the Gln(64)(E7) N(epsilon)H(1). The E7 and E11 residues appear to control the direction and the extent of tilt of the bound ligand. Furthermore, the tilt of the ligand has no obvious influence on the heme heterogeneity of cyanide ligation for isomer A/B of the wild type and mutant protein, indicating that factors other than steric effects, such as polarity of heme pocket, impacts on ligand binding affinity.     

Modelling of lysozyme binding to a cation exchange surface at atomic detail: the role of flexibility

Different approaches were made to predict the adsorbed orientation based on rigid, flexible, or a mixture of both models. To determine the role of flexibility during adsorption, the orientation of lysozyme adsorbed to a negatively charged ligand surface was predicted by a rigid and a flexible model based on two differing protein structures at atomic resolution. For the rigid model, the protein structures were placed at different distances from the ligand surface and the electrostatic interaction energy was calculated for all possible orientations. The results were compared to a flexible model where the binding to the ligand surface was modeled by multiple molecular dynamics simulations starting with 14 initial orientations. Different aspects of the adsorption process were not covered by the rigid model and only detectable by the flexible model. Whereas the results of the rigid model depended sensitively on the protein-surface distance and the protein structure, the preferred orientation obtained by the flexible model was closer to a previous experimental determined orientation, robust toward the initial orientation and independent of the initial protein structure. Additionally, it was possible to obtain insights into the preferred binding process of lysozyme on a negatively charged surface by the flexible model.     

Electrostatic orientation of the electron-transfer complex between plastocyanin and cytochrome c.

To understand the specificity and efficiency of protein-protein interactions promoting electron transfer, we evaluated the role of electrostatic forces in precollision orientation by the development of two new methods, computer graphics alignment of protein electrostatic fields and a systematic orientational search of intermolecular electrostatic energies for two proteins at present separation distances. We applied these methods to the plastocyanin/cytochrome c interaction, which is faster than random collision, but too slow for study by molecular dynamics techniques. Significant electrostatic potentials were concentrated on one-fourth (969 A2) of the plastocyanin surface, with the greatest negative potential centered on the Tyr-83 hydroxyl within the acidic patch, and on one-eighth (632 A2) of the cytochrome c surface, with the greatest positive potential centered near the exposed heme edge. Coherent electrostatic fields occurred only over these regions, suggesting that local, rather than global, charge complementarity controls productive recognition. The three energetically favored families of pre-collision orientations all directed the positive region surrounding the heme edge of cytochrome c toward the acidic patch of plastocyanin but differed in heme plane orientation. Analysis of electrostatic fields, electrostatic energies of precollision orientations with 12 and 6 A separation distances, and surface topographies suggested that the favored orientations should converge to productive complexes promoting a single electron-transfer pathway from the cytochrome c heme edge to Tyr-83 of plastocyanin. Direct interactions of the exposed Cu ligand in plastocyanin with the cytochrome c heme edge are not unfavorable sterically or electrostatically but should occur no faster than randomly, indicating that this is not the primary pathway for electron transfer.     

Conformational relaxation in hemoproteins: the cytochrome P-450cam case

Photodissociation of (CO)P-450(cam)(substrate) complexes was found to trigger a conformational relaxation process that interferes with ligand rebinding at temperatures as low as 140 K even though the protein conformational substates (CS(1)) remain frozen. To analyze the rebinding and relaxation kinetics, we developed a model that takes the distribution of relaxation rates explicitly into account and in which rebinding and relaxation rates are connected by a linear free energy relation. In all complexes heme relaxation occurs first and is probably faster than 100 ns even at 77 K. This is the only process found in substrate-free P-450(cam). Above 140 K and in the presence of a substrate, this initial, fast rebinding state (P) progressively relaxes to another state (P degrees ) in which rebinding is slower. The relaxation rate is independent of solvent rigidity and is governed by the protein's internal dynamics. Rebinding enthalpies in P and P degrees as well as the enthalpy shift brought about by relaxation correlate with the substrate propensity to block access to the iron site. In P degrees the barrier is higher because the substrate is closer to the heme normal and exerts more steric repulsion for CO binding. The relaxation process implies the return of substrate and heme to their ligand-free positions in which access to the heme is reduced.     

Time-resolved resonance Raman study on ultrafast structural relaxation and vibrational cooling of photodissociated carbonmonoxy myoglobin

A localized small structural change is converted to a higher order conformational change of protein and extends to a mesoscopic scale to induce a physiological function. To understand such features of protein, ultrafast dynamics of myoglobin (Mb) following photolysis of carbon monoxide were investigated. Recent results are summarized here with a stress on structural and vibrational energy relaxation. The core expansion of heme takes place within 2 ps but the out of plane displacement of the heme iron and the accompanying protein conformational change occur in 10 and 100 s of the picosecond regimes, respectively. Unexpectedly, it was found from UV resonance Raman spectra that Trp7 in the N-terminal region and Tyr151 in the C-terminal region undergo appreciable structural changes upon ligand binding-dissociation while Tyr104, Tyr146, and Trp14 do not. Because of the communication between the movements of these surface residues and the heme iron, the rate of spectral change of the iron-histidine (Fe- His) stretching band after CO photodissociation is influenced by the viscosity of solvent. Temporal changes of the anti-Stokes Raman intensity demonstrated immediate generation of vibrationally excited heme upon photodissociation and its decay with a time constant of 1-2 ps.     

Increasing the conformational stability by replacement of heme axial ligand in c-type cytochrome

To investigate the role of the heme axial ligand in the conformational stability of c-type cytochrome, we constructed M58C and M58H mutants of the red alga Porphyra yezoensis cytochrome c(6) in which the sixth heme iron ligand (Met58) was replaced with Cys and His residues, respectively. The Gibbs free energy change for unfolding of the M58H mutant in water (DeltaG degrees (unf)=1.48 kcal/mol) was lower than that of the wild-type (2.43 kcal/mol), possibly due to the steric effects of the mutation on the apoprotein structure. On the other hand, the M58C mutant exhibited a DeltaG degrees (unf) of 5.45 kcal/mol, a significant increase by 3.02 kcal/mol compared with that of wild-type. This increase was possibly responsible for the sixth heme axial bond of M58C mutant being more stable than that of wild-type according to the heme-bound denaturation curve. Based on these observations, we propose that the sixth heme axial ligand is an important key to determine the conformational stability of c-type cytochromes, and the sixth Cys heme ligand will give stabilizing effects.     

Structural analysis of cytochromes P450 shows differences in flexibility of heme 2- and 4-vinyls

Structural analysis of the orientations of heme vinyl side chains was carried out using the published crystallographic data for different cytochromes P450. Torsional angles (tau, C(alpha)C(beta)-C(a)C(b)) show different distributions for the vinyls in positions 2 and 4. Whereas the orientation of 2-vinyls is rather restricted (tau between -120 degrees and -180 degrees ), the 4-vinyls have a much higher mobility over almost the entire conformational space. On the basis of the empirical correlation recently reported for peroxidases (M.P. Marzocchi, G. Smulevich, Relationship between heme vinyl conformation and the protein matrix in peroxidases, J. Raman Spectrosc. 34 (2003), 725-736), an attempt has been made to compare the observed vinyl orientations with the experimental frequencies of the vinyl stretching vibrational modes. The data for P450 proteins do not exactly match the peroxidase-derived function, although a qualitatively similar relationship is likely to exist. Differences between P450 forms suggest a variability in heme-region flexibility and in communication with the rest of enzyme.