Human myoglobin recognition of oxygen: dynamics of the energy landscape

Ultrafast dynamics of heme distortion in the O2-sensor of a thermophilic anaerobe bacterium

Heme-Nitric oxide and Oxygen binding protein domains (H-NOX) are found in signaling pathways of both prokaryotes and eukaryotes and share sequence homology with soluble guanylate cyclase, the mammalian NO receptor. In bacteria, H-NOX is associated with kinase or methyl accepting chemotaxis domains. In the O₂-sensor of the strict anaerobe Caldanaerobacter tengcongensis (Ct H-NOX) the heme appears highly distorted after O₂ binding, but the role of heme distortion in allosteric transitions was not yet evidenced. Here, we measure the dynamics of the heme distortion triggered by the dissociation of diatomics from Ct H-NOX using transient electronic absorption spectroscopy in the picosecond to millisecond time range. We obtained a spectroscopic signature of the heme flattening upon O₂ dissociation. The heme distortion is immediately (<1 ps) released after O₂ dissociation to produce a relaxed state. This heme conformational change occurs with different proportions depending on diatomics as follows: CO < NO < O₂. Our time-resolved data demonstrate that the primary structural event of allostery is the heme distortion in the Ct H-NOX sensor, contrastingly with hemoglobin and the human NO receptor, in which the primary structural events are respectively the motion of the proximal histidine and the rupture of the iron-histidine bond.

Effect of Single-Point Mutations on Nitric Oxide Rebinding and the Thermodynamic Stability of Myoglobin

The effect of single amino acid mutations on the rebinding dynamics of nitrogen monoxide (NO) to myoglobin is investigated using reactive molecular dynamics simulations. In particular, mutations of residues surrounding the heme-active site (Leu29, His64, Val68) were considered. Consistent with experiments, all mutations studied here have a significant effect on the kinetics of the NO-rebinding process, which consists of a rapid (several 10 ps) and a slow (100s of ps) time scale. For all modifications considered, the time scales and rebinding fractions agree to within a few percents with results from experiments by adjusting one single, physically meaningful, conformationally averaged quantity: the asymptotic energy separation between the NO-bound (²A) and photodissociated (⁴A) states. It is furthermore shown that the thermodynamic stability of wild-type versus mutant Mb for the ligand-free and ligand-bound variants of the protein can be described by the same computational model. Therefore, ligand kinetics and thermodynamics are related in a direct fashion akin to Φ-value analysis, which establishes a relationship between protein folding rates and thermal stability of proteins.

Observation time scale, free-energy landscapes, and molecular symmetry

When structures that interconvert on a given time scale are lumped together, the corresponding free-energy surface becomes a function of the observation time. This view is equivalent to grouping structures that are connected by free-energy barriers below a certain threshold. We illustrate this time dependence for some benchmark systems, namely atomic clusters and alanine dipeptide, highlighting the connections to broken ergodicity, local equilibrium, and "feasible" symmetry operations of the molecular Hamiltonian.

Heterogeneity of single molecule FRET signals reveals multiple active ribosome subpopulations

Single molecule methods have revealed that heterogeneity is common in biological systems. However, interpretations of the complex signals are challenging. By tracking the fluorescence resonance energy transfer (FRET) signals between the A-site tRNA and L27 protein in single ribosomes, we attempt to develop a qualitative method to subtract the inherent patterns of the heterogeneous single molecule FRET data. Seven ribosome subpopulations are identified using this method and spontaneous exchanges among these subpopulations are observed. All of the pretranslocation subpopulations are competent in real-time translocation, but via distinguished pathways. These observations suggest that the ribosome may function through multiple reaction pathways.

Single-molecule study of ribosome hierarchic dynamics at the peptidyl transferase center

During protein biosynthesis the ribosome moves along mRNA in steps of precisely three nucleotides. The mechanism for this ribosome motion remains elusive. Using a classification algorithm to sort single-molecule fluorescence resonance energy transfer data into subpopulations, we found that the ribosome dynamics detected at the peptidyl transferase center are highly inhomogeneous. The pretranslocation complex has at least four subpopulations that sample two hybrid states, whereas the posttranslocation complex is mainly static. We observed transitions among the ribosome subpopulations under various conditions, including 1), in the presence of EF-G; 2), spontaneously; 3), in different buffers, and 4), bound to antibiotics. Therefore, these subpopulations represent biologically active ribosomes. One key observation indicates that the Hy2 hybrid state only exists in a fluctuating ribosome subpopulation, which prompts us to propose that ribosome dynamics are hierarchically arranged. This proposal may have important implications for the regulation of cellular translation rates.

Ultrafast proteinquake dynamics in cytochrome c

We report here our systematic studies of the heme dynamics and induced protein conformational relaxations in two redox states of ferric and ferrous cytochrome c upon femtosecond excitation. With a wide range of probing wavelengths from the visible to the UV and a site-directed mutation we unambiguously determined that the protein dynamics in the two states are drastically different. For the ferrous state the heme transforms from 6-fold to 5-fold coordination with ultrafast ligand dissociation in less than 100 fs, followed by vibrational cooling within several picoseconds, but then recombining back to its original 6-fold coordination in 7 ps. Such impulsive bond breaking and late rebinding generate proteinquakes and strongly perturb the local heme site and shake global protein conformation, which were found to completely recover in 13 and 42 ps, respectively. For the ferric state the heme however maintains its 6-fold coordination. The dynamics mainly occur at the local site, including ultrafast internal conversion in hundreds of femtoseconds, vibrational cooling on the similar picosecond time scale, and complete ground-state recovery in 10 ps, and no global conformation relaxation was observed.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

Dynamics of ligand substitution in labile cobalt complexes resolved by ultrafast T-jump

Ligand exchange of hydrated metal complexes is common in chemical and biological systems. Using the ultrafast T-jump, we examined this process, specifically the transformation of aqua cobalt (II) complexes to their fully halogenated species. The results reveal a stepwise mechanism with time scales varying from hundreds of picoseconds to nanoseconds. The dynamics are significantly faster when the structure is retained but becomes rate-limited when the octahedral-to-tetrahedral structural change bottlenecks the transformation. Evidence is presented, from bimolecular kinetics and energetics (enthalpic and entropic), for a reaction in which the ligand assists the displacement of water molecules, with the retention of the entering ligand in the activated state. The reaction time scale deviates by one to two orders of magnitude from that of ionic diffusion, suggesting the involvement of a collisional barrier between the ion and the much larger complex.

Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin

Myoglobin is a globular protein involved in oxygen storage and transport. No consensus yet exists on the atomic level mechanism by which oxygen and other small nonpolar ligands move between the myoglobin's buried heme, which is the ligand binding site, and surrounding solvent. This study uses room temperature molecular dynamics simulations to provide a complete atomic level picture of ligand migration in myoglobin. Multiple trajectories--providing a cumulative total of 7 micros of simulation--are analyzed. Our simulation results are consistent with and tie together previous experimental findings. Specifically, we characterize: (i) Explicit full trajectories in which the CO ligand shuttles between the internal binding site and the solvent and (ii) pattern and structural origins of transient voids available for ligand migration. The computations are performed both in sperm whale myoglobin wild-type and in sperm whale V68F myoglobin mutant, which is experimentally known to slow ligand-binding kinetics. On the basis of these independent, but mutually consistent ligand migration and transient void computations, we find that there are two discrete dynamical pathways for ligand migration in myoglobin. Trajectory hops between these pathways are limited to two bottleneck regions. Ligand enters and exits the protein matrix in common identifiable portals on the protein surface. The pathways are located in the "softer" regions of the protein matrix and go between its helices and in its loop regions. Localized structural fluctuations are the primary physical origin of the simulated CO migration pathways inside the protein.

Toward quantitative simulations of carbon monoxide escape pathways in myoglobin

Straightforward molecular dynamics trajectories have been computed to explore the diffusion of carbon monoxide through myoglobin. The classical equations of motion were integrated for 2 ns and the resulting pathways analyzed. Two types of runs were examined. Type i: Myoglobin and a ligand embedded in a periodic box with 9996 water molecules; the water molecules are rigid but the bonds of the protein are flexible. Type ii: Myoglobin with a solvation shell (153 water molecules) in which all bond lengths are fixed. In trajectories of type i, the diffusing ligand visits a significant part of the protein matrix and was not constrained to the proximity of the heme pocket before escaping. The maximum time of the trajectories was 2 ns. It was shorter if the ligand escaped earlier. Two ligands (from a total of 88) escape to the solvent from nonclassical gates (non-E-helix gates). In trajectories of type ii, the overall fluctuations of the protein are smaller and the ligand explores significantly smaller internal space. The escape rate from type ii trajectories (11 of 400) is comparable to type i and is not dramatically different from experiment (1 of 100). Interestingly, the two simulations with comparable rates sampled different pathways. In trajectories of type ii, we observe escapes from the classical gate (His 64) and from the Xe4 cavity. Further studies (that are underway) are required to define the escape pathways and the overall rate.

Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.

4D ultrafast electron diffraction, crystallography, and microscopy

In this review, we highlight the progress made in the development of 4D ultrafast electron diffraction (UED), crystallography (UEC), and microscopy (UEM) with a focus on concepts, methodologies, and prototypical applications. The joint atomic-scale resolutions in space and time, and sensitivity reached, make it possible to determine complex transient structures and assemblies in different phases. These applications include studies of isolated chemical reactions (molecular beams), interfaces, surfaces and nanocrystals, self-assembly, and 2D crystalline fatty-acid bilayers. In 4D UEM, we are now able, using timed, single-electron packets, to image nano-to-micro scale structures of materials and biological cells. Future applications of these methods are foreseen across areas of physics, chemistry, and biology.

Protein surface hydration mapped by site-specific mutations

Water motion at protein surfaces is fundamental to protein structure, stability, dynamics, and function. By using intrinsic tryptophans as local optical probes, and with femtosecond resolution, it is possible to probe surface-water motions in the hydration layer. Here, we report our studies of local hydration dynamics at the surface of the enzyme Staphylococcus nuclease using site-specific mutations. From these studies of the WT and four related mutants, which change local charge distribution and structure, we are able to ascertain the contribution to solvation by protein side chains as relatively insignificant. We determined the time scales of hydration to be 3-5 ps and 100-150 ps. The former is the result of local librational/rotational motions of water near the surface; the latter is a direct measure of surface hydration assisted by fluctuations of the protein. Experimentally, these hydration dynamics of the WT and the four mutants are also consistent with results of the total dynamic Stokes shifts and fluorescence emission maxima and are correlated with their local charge distribution and structure. We discuss the role of protein fluctuation on the time scale of labile hydration and suggest reexamination of recent molecular dynamics simulations.

Dynamics of clusters: from elementary to biological structures

Between isolated atoms or molecules and bulk materials there lies a class of unique structures, known as clusters, that consist of a few to hundreds of atoms or molecules. Within this range of "nanophase," many physical and chemical properties of the materials evolve as a function of cluster size, and materials may exhibit novel properties due to quantum confinement effects. Understanding these phenomena is in its own rights fundamental, but clusters have the additional advantage of being controllable model systems for unraveling the complexity of condensed-phase and biological structures, not to mention their vanguard role in defining nanoscience and nanotechnology. Over the last two decades, much progress has been made, and this short overview highlights our own involvement in developing cluster dynamics, from the first experiments on elementary systems to model systems in the condensed phase, and on to biological structures.

Dynamics of nitric oxide rebinding and escape in horseradish peroxidase

Ultrafast kinetic measurements of NO rebinding to horseradish peroxidase (HRP) are reported for the first time. The geminate kinetics are found to be exponential for all HRP samples studied. The ferric forms of HRP have NO geminate recombination time constants in the range of 15-30 ps, while the ferrous form has a time constant of approximately 7 ps. The simple exponential NO geminate kinetics found for HRP demonstrate that heme relaxation is not the underlying source of the nonexponential NO rebinding in myoglobin (Mb). The NO ligand escape rates from HRP are also determined, and they are found to depend dramatically on the presence or absence of the competitive inhibitor benzohydroxamic acid (BHA). The kinetic results indicate that, in contrast to Mb, there is direct solvent access to the distal heme pocket of HRP.

Studying reactive processes with classical dynamics: rebinding dynamics in MbNO

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mus of total simulation time) for the rebinding of NO to myoglobin. The unbound surface (Fe...NO) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.