Primary protein response after ligand photodissociation in carbonmonoxy myoglobin

Capturing ultrafast energy flow of a heme protein in crowded milieu

Energy flow in biomolecules is a dynamic process vital for understanding health, disease, and applications in biotechnology and medicine. In crowded environments, where biomolecular functions are modulated, comprehending energy flow becomes crucial for accurately understanding cellular processes like signaling and subsequent functions. This study employs ultrafast transient absorption spectroscopy to demonstrate energy funneling from the photoexcited heme of bovine heart cytochrome c to the protein exterior, in the presence of common synthetic (Dextran 40, Ficoll 70, PEG 8 and Dextran 70) and protein-based (BSA and β-LG) crowders. The through-space energy transfer mode for ferric and the methionine rebinding mode for ferrous cytochrome c show the strongest solvent coupling. The heterogeneous behaviour of crowders, influenced by crowder-protein interactions and caging effects at certain higher concentrations, reveal diverse trends. Notably, protein crowders perturb all transport routes of vibrational energy transfer, causing delays in energy transfer processes. These findings provide significant insights into the basic tenets of energy flow, one of the most fundamental processes, in crowded cellular environments.

Contact-Mediated Retinal-Opsin Coupling Enables Proton Pumping in Gloeobacter Rhodopsin

When a chromophore embedded in a photoreceptive protein undergoes a reaction upon photoexcitation, the photoreaction triggers structural changes in the protein moiety that are necessary for the function of the protein. It is thus essential to elucidate the coupling between the chromophore and protein moiety to understand the functional mechanism for photoreceptive proteins, but the mechanism by which this coupling occurs remains poorly understood. Here, we show that nonbonded atomic contacts play an essential role in driving functionally important structural changes following photoisomerization of the chromophore in Gloeobacter rhodopsin (GR). Time-resolved ultraviolet resonance Raman spectroscopy revealed that the substitution of Trp222, which contacts with methyl groups of the retinal chromophore, with a Phe residue reduced the extent of structural change. The proton-pumping activity of the GR mutant was as small as 9% of that of the wild type. Time-resolved visible absorption and resonance Raman spectra showed that the photocycle of the mutant proceeded to the L intermediate following the all-trans to 13-cis photoisomerization step but did not result in the deprotonation of the chromophore. The present results demonstrate that the atomic contacts between the chromophore and the Trp222 side chain induce the structural changes necessary for proton transfer. The requirement for dense atomic packing in a protein structure for the efficient propagation of structural changes through a coupling mechanism is discussed.

Origin of Protein Quake: Energy Waves Conducted by a Precise Mechanical Machine

A long-standing challenge in protein biophysics is to understand protein quake in myoglobin─the structural dynamics responsible for redistributing the excess heme energy after photolysis. Despite extensive efforts, the molecular mechanism of this process remains elusive. Using the energy flow theory, we uncovered a fundamental new phenomenon: the heme energy is redistributed by sinusoidal waves with a ubiquitous fundamental frequency and two overtones. The energy waves emanate from the heme into the myoglobin backbone via a conduit of five consecutive dihedrals of the proximal histidine and then travel quickly along the backbone to reach sidechains across the protein. This mechanism is far more effective than the diffusion-based mechanism from previous studies because waves are systematic while diffusion is random. To propagate energy waves, coordinates must cooperate, resulting in collective modes that are singular vectors of the generalized work functional. These modes show task partitioning: a handful of high-energy modes generate large-scale breathing motion, which loosens up the protein matrix to enable hundreds of low-energy vibrational modes for energy transduction.

Concerted Motions and Molecular Function: What Physical Chemistry We Can Learn from Light-Driven Ion-Pumping Rhodopsins

Transmembrane ion gradients are generated and maintained by ion-pumping proteins in cells. Light-driven ion-pumping rhodopsins are retinal-containing proteins found in archaea, bacteria, and eukarya. Photoisomerization of the retinal chromophore induces structural changes in the protein, allowing the transport of ions in a particular direction. Understanding unidirectional ion transport by ion-pumping rhodopsins is an exciting challenge for biophysical chemistry. Concerted changes in ion-binding affinities of the ion-binding sites in proteins are key to unidirectional ion transport, as is the coupling between the chromophore and the protein moiety to drive the concerted motions regulating ion-binding affinities. The commonality of ion-pumping rhodopsin protein structures and the diversity of their ion-pumping functions suggest universal principles governing ion transport, which would be widely applicable to molecular systems. In this Perspective, I review the insights obtained from previous studies on rhodopsins and discuss future perspectives.

Tracking Ultrafast Structural Dynamics by Time-Domain Raman Spectroscopy

In traditional Raman spectroscopy, narrow-band light is irradiated on a sample, and its inelastic scattering, i.e., Raman scattering, is detected. The energy difference between the Raman scattering and the incident light corresponds to the vibrational energy of the molecule, providing the Raman spectrum that contains rich information about the molecular-level properties of the materials. On the other hand, by using ultrashort optical pulses, it is possible to induce Raman-active coherent nuclear motion of the molecule and to observe the molecular vibration in real time. Moreover, this time-domain Raman measurement can be combined with femtosecond photoexcitation, triggering chemical changes, which enables tracking ultrafast structural dynamics in a form of “time-resolved” time-domain Raman spectroscopy, also known as time-resolved impulsive stimulated Raman spectroscopy. With the advent of stable, ultrashort laser pulse sources, time-resolved impulsive stimulated Raman spectroscopy now realizes high sensitivity and a wide detection frequency window from THz to 3000 cm^–1, and has seen success in unveiling the molecular mechanisms underlying the efficient functions of complex molecular systems. In this Perspective, we overview the present status of time-domain Raman spectroscopy, particularly focusing on its application to the study of femtosecond structural dynamics. We first explain the principle and a brief history of time-domain Raman spectroscopy and then describe the apparatus and recent applications to the femtosecond dynamics of complex molecular systems, including proteins, molecular assemblies, and functional materials. We also discuss future directions for time-domain Raman spectroscopy, which has reached a status allowing a wide range of applications.

Gas-Phase Resonance Raman Spectroscopy Combined with IR-Laser Ablation of a Droplet Beam: Local Structural Analysis of Myoglobin

Gas phase spectroscopy is a powerful tool for examining fundamental chemical structures and properties free from solvent molecules. We developed a gas-phase resonance Raman spectroscopy combined with IR-laser ablation of a droplet beam, which allowed us to elucidate local structures around chromophores in gas-phase proteins and DNAs. To demonstrate the potential of this approach, we applied this method to myoglobin, one of the heme proteins, and elucidated its structures in the gas phase and in aqueous solution. The experimental spectra are compared with calculated spectra of stable heme structures for the structural determination. These results show the oxidation/spin states of the Fe atom in myoglobin in the gas phase and were compared with the aqueous solution from the obtained resonant Raman spectra. The present method gives an important tool to investigate the gas-phase structure of large biomolecules.

Origin of complexity in haemoglobin evolution

Most proteins associate into multimeric complexes with specific architectures1,2, which often have functional properties such as cooperative ligand binding or allosteric regulation3. No detailed knowledge is available about how any multimer and its functions arose during evolution. Here we use ancestral protein reconstruction and biophysical assays to elucidate the origins of vertebrate haemoglobin, a heterotetramer of paralogous α- and β-subunits that mediates respiratory oxygen transport and exchange by cooperatively binding oxygen with moderate affinity. We show that modern haemoglobin evolved from an ancient monomer and characterize the historical 'missing link' through which the modern tetramer evolved-a noncooperative homodimer with high oxygen affinity that existed before the gene duplication that generated distinct α- and β-subunits. Reintroducing just two post-duplication historical substitutions into the ancestral protein is sufficient to cause strong tetramerization by creating favourable contacts with more ancient residues on the opposing subunit. These surface substitutions markedly reduce oxygen affinity and even confer cooperativity, because an ancient linkage between the oxygen binding site and the multimerization interface was already an intrinsic feature of the protein's structure. Our findings establish that evolution can produce new complex molecular structures and functions via simple genetic mechanisms that recruit existing biophysical features into higher-level architectures.

Molecular dynamics simulations and linear response theories jointly describe biphasic responses of myoglobin relaxation and reveal evolutionarily conserved frequent communicators

In this study, we provide a time-dependent mechanical model, taking advantage of molecular dynamics simulations, quasiharmonic analysis of molecular dynamics trajectories, and time-dependent linear response theories to describe vibrational energy redistribution within the protein matrix. The theoretical description explained the observed biphasic responses of specific residues in myoglobin to CO-photolysis and photoexcitation on heme. The fast responses were found to be triggered by impulsive forces and propagated mainly by principal modes <40 cm⁻¹. The predicted fast responses for individual atoms were then used to study signal propagation within the protein matrix and signals were found to propagate ~8 times faster across helices (4076 m/s) than within the helices, suggesting the importance of tertiary packing in the sensitivity of proteins to external perturbations. We further developed a method to integrate multiple intramolecular signal pathways and discover frequent “communicators”. These communicators were found to be evolutionarily conserved including those distant from the heme.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Ultrafast anisotropic protein quake propagation after CO photodissociation in myoglobin

"Protein quake" denotes the dissipation of excess energy across a protein, in response to a local perturbation such as the breaking of a chemical bond or the absorption of a photon. Femtosecond time-resolved small- and wide-angle X-ray scattering (TR-SWAXS) is capable of tracking such ultrafast protein dynamics. However, because the structural interpretation of the experiments is complicated, a molecular picture of protein quakes has remained elusive. In addition, new questions arose from recent TR-SWAXS data that were interpreted as underdamped oscillations of an entire protein, thus challenging the long-standing concept of overdamped global protein dynamics. Based on molecular-dynamics simulations, we present a detailed molecular movie of the protein quake after carbon monoxide (CO) photodissociation in myoglobin. The simulations suggest that the protein quake is characterized by a single pressure peak that propagates anisotropically within 500 fs across the protein and further into the solvent. By computing TR-SWAXS patterns from the simulations, we could interpret features in the reciprocal-space SWAXS signals as specific real-space dynamics, such as CO displacement and pressure wave propagation. Remarkably, we found that the small-angle data primarily detect modulations of the solvent density but not oscillations of the bare protein, thereby reconciling recent TR-SWAXS experiments with the notion of overdamped global protein dynamics.

NO binding kinetics in myoglobin investigated by picosecond Fe K-edge absorption spectroscopy

Diatomic ligands in hemoproteins and the way they bind to the active center are central to the protein's function. Using picosecond Fe K-edge X-ray absorption spectroscopy, we probe the NO-heme recombination kinetics with direct sensitivity to the Fe-NO binding after 532-nm photoexcitation of nitrosylmyoglobin (MbNO) in physiological solutions. The transients at 70 and 300 ps are identical, but they deviate from the difference between the static spectra of deoxymyoglobin and MbNO, showing the formation of an intermediate species. We propose the latter to be a six-coordinated domed species that is populated on a timescale of ∼ 200 ps by recombination with NO ligands. This work shows the feasibility of ultrafast pump-probe X-ray spectroscopic studies of proteins in physiological media, delivering insight into the electronic and geometric structure of the active center.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks

The ligand migration network for O2-diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k-means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k-means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.

Insights into Protein Structure and Dynamics by Ultraviolet and Visible Resonance Raman Spectroscopy

Raman spectroscopy is a form of vibrational spectroscopy based on inelastic scattering of light. In resonance Raman spectroscopy, the wavelength of the incident light falls within an absorption band of a chromophore, and this overlap of excitation and absorption energy greatly enhances the Raman scattering efficiency of the absorbing species. The ability to probe vibrational spectra of select chromophores within a complex mixture of molecules makes resonance Raman spectroscopy an excellent tool for studies of biomolecules. In this Current Topic, we discuss the type of molecular insights obtained from steady-state and time-resolved resonance Raman studies of a prototypical photoactive protein, rhodopsin. We also review recent efforts in ultraviolet resonance Raman investigations of soluble and membrane-associated biomolecules, including integral membrane proteins and antimicrobial peptides. These examples illustrate that resonance Raman is a sensitive, selective, and practical method for studying the structures of biological molecules, and the molecular bonding, geometry, and environments of protein cofactors, the backbone, and side chains.

Observing heme doming in myoglobin with femtosecond X-ray absorption spectroscopy

We report time-resolved X-ray absorption measurements after photolysis of carbonmonoxy myoglobin performed at the LCLS X-ray free electron laser with nearly 100 fs (FWHM) time resolution. Data at the Fe K-edge reveal that the photoinduced structural changes at the heme occur in two steps, with a faster (∼70 fs) relaxation preceding a slower (∼400 fs) one. We tentatively attribute the first relaxation to a structural rearrangement induced by photolysis involving essentially only the heme chromophore and the second relaxation to a residual Fe motion out of the heme plane that is coupled to the displacement of myoglobin F-helix.

Ultrafast myoglobin structural dynamics observed with an X-ray free-electron laser

Light absorption can trigger biologically relevant protein conformational changes. The light-induced structural rearrangement at the level of a photoexcited chromophore is known to occur in the femtosecond timescale and is expected to propagate through the protein as a quake-like intramolecular motion. Here we report direct experimental evidence of such 'proteinquake' observed in myoglobin through femtosecond X-ray solution scattering measurements performed at the Linac Coherent Light Source X-ray free-electron laser. An ultrafast increase of myoglobin radius of gyration occurs within 1 picosecond and is followed by a delayed protein expansion. As the system approaches equilibrium it undergoes damped oscillations with a ~3.6-picosecond time period. Our results unambiguously show how initially localized chemical changes can propagate at the level of the global protein conformation in the picosecond timescale.

Ligand-induced protein responses and mechanical signal propagation described by linear response theories

In this study, a general linear response theory (LRT) is formulated to describe time-dependent and -independent protein conformational changes upon CO binding with myoglobin. Using the theory, we are able to monitor protein relaxation in two stages. The slower relaxation is found to occur from 4.4 to 81.2 picoseconds and the time constants characterized for a couple of aromatic residues agree with those observed by UV Resonance Raman (UVRR) spectrometry and time resolved x-ray crystallography. The faster "early responses", triggered as early as 400 femtoseconds, can be best described by the theory when impulse forces are used. The newly formulated theory describes the mechanical propagation following ligand-binding as a function of time, space and types of the perturbation forces. The "disseminators", defined as the residues that propagate signals throughout the molecule the fastest among all the residues in protein when perturbed, are found evolutionarily conserved and the mutations of which have been shown to largely change the CO rebinding kinetics in myoglobin.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

Using THz time-scale infrared spectroscopy to examine the role of collective, thermal fluctuations in the formation of myoglobin allosteric communication pathways and ligand specificity

In this investigation we use THz time-scale spectroscopy to conduct an initial set of studies on myoglobin with the aim of providing further insight into the global, collective thermal fluctuations in the protein that have been hypothesized to play a prominent role in the dynamic formation of transient ligand channels as well as in shaping the molecular level basis for ligand discrimination. Using the two ligands O2 and CO, we have determined that the perturbation from the heme-ligand complex has a strong influence on the characteristics of the myoglobin collective dynamics that are excited upon binding. Further, the differences detected in the collective protein motions in Mb-O2 compared with those in Mb-CO appear to be intimately tied with the pathways of long-range allosteric communication in the protein, which ultimately determine the trajectories selected by the respective ligands on the path to and from the heme-binding cavity.

Identification of functionally key residues in AMPA receptor with a thermodynamic method

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

Two-Dimensional Stimulated Ultraviolet Resonance Raman Spectra of Tyrosine and Tryptophan; A Simulation Study

We report an ab-initio simulation study of the ultrafast broad bandwidth ultraviolet (UV) stimulated resonance Raman spectra (SRRS) of L-tyrosine, L-tryptophan and trans-L-tryptophan-L-tyrosine (WY) dipeptide. Two-pulse one-dimensional (1D) SRRS and three-pulse 2D SRRS that reveal inter- and intra-residue vibrational coorelations are simulated using electronically resonant or preresonant pulse configurations that select the Raman signal and discriminate against excited state pathways. Multimode effects are incorporated via the cumulant expansion. The 2D SRRS technique is more sensitive to residue couplings than spontaneous Raman.

Two-Dimensional Stimulated Ultraviolet Resonance Raman Spectra of Tyrosine and Tryptophan; A Simulation Study

We report an ab-initio simulation study of the ultrafast broad bandwidth ultraviolet (UV) stimulated resonance Raman spectra (SRRS) of L-tyrosine, L-tryptophan and trans-L-tryptophan-L-tyrosine (WY) dipeptide. Two-pulse one-dimensional (1D) SRRS and three-pulse 2D SRRS that reveal inter- and intra-residue vibrational coorelations are simulated using electronically resonant or preresonant pulse configurations that select the Raman signal and discriminate against excited state pathways. Multimode effects are incorporated via the cumulant expansion. The 2D SRRS technique is more sensitive to residue couplings than spontaneous Raman.

Oxygen migration pathways in NO-bound truncated hemoglobin

Atomistic simulations of dioxygen (O(2)) dynamics and migration in nitric oxide-bound truncated Hemoglobin N (trHbN) of Mycobacterium tuberculosis are reported. From more than 100 ns of simulations the connectivity network involving the metastable states for localization of the O(2) ligand is built and analyzed. It is found that channel I is the primary entrance point for O(2) whereas channel II is predominantly an exit path although access to the protein active site is also possible. For O(2) a new site compared to nitric oxide, from which reaction with the heme group can occur, was found. As this site is close to the heme iron, it could play an important role in the dioxygenation mechanism as O(2) can remain there for hundreds of picoseconds after which it can eventually leave the protein, while NO is localized in Xe2. The present study supports recent experimental work which proposed that O(2) docks in alternative pockets than Xe close to the reactive site. Similar to other proteins, a phenylalanine residue (Phe62) plays the role of a gate along the access route in channel I. The most highly connected site is the Xe3 pocket which is a "hub" and free energy barriers between the different metastable states are ≈1.5 kcal mol(-1) which allows facile O(2) migration within the protein.

Correlation of TrpGly and GlyTrp Rotamer Structure with W7 and W10 UV Resonance Raman Modes and Fluorescence Emission Shifts

Tryptophyl glycine (TrpGly) and glycyl tryptophan (GlyTrp) dipeptides at pH 5.5 and pH 9.3 show a pattern of fluorescence emission shifts with the TrpGly zwitterion emission solely blue shifted. This pattern is matched by shifts in the UV resonance Raman (UVRR) W10 band position and the W7 Fermi doublet band ratio. Ab initio calculations show that the 1340 cm⁻¹ band of the W7 doublet is composed of three modes, two of which determine the W7 band ratios for the dipeptides. Molecular dynamics simulations show that the dipeptides take on two conformations: one with the peptide backbone extended; one with the backbone curled over the indole. The dihedral angle critical to these conformations is χ¹ and takes on three discrete values. Only the TrpGly zwitterion spends an appreciable amount of time in the extended backbone conformation as this is stabilized by two hydrogen bonds with the terminal amine cation. According to a Stark effect model, a positive charge near the pyrrole keeps the ¹L_a transition at high energy, limiting fluorescence emission red shift, as observed for the TrpGly zwitterion. The hydrogen bond stabilized backbone provides a rationale for the C_methylene-C_α-C_carbonyl W10 symmetric stretch that is unique to the TrpGly zwitterion.

Site-specific protein dynamics in communication pathway from sensor to signaling domain of oxygen sensor protein, HemAT-Bs: Time-resolved Ultraviolet Resonance Raman Study

HemAT-Bs is a heme-based signal transducer protein responsible for aerotaxis. Time-resolved ultraviolet resonance Raman (UVRR) studies of wild-type and Y70F mutant of the full-length HemAT-Bs and the truncated sensor domain were performed to determine the site-specific protein dynamics following carbon monoxide (CO) photodissociation. The UVRR spectra indicated two phases of intensity changes for Trp, Tyr, and Phe bands of both full-length and sensor domain proteins. The W16 and W3 Raman bands of Trp, the F8a band of Phe, and the Y8a band of Tyr increased in intensity at hundreds of nanoseconds after CO photodissociation, and this was followed by recovery in ∼50 μs. These changes were assigned to Trp-132 (G-helix), Tyr-70 (B-helix), and Phe-69 (B-helix) and/or Phe-137 (G-helix), suggesting that the change in the heme structure drives the displacement of B- and G-helices. The UVRR difference spectra of the sensor domain displayed a positive peak for amide I in hundreds of nanoseconds after photolysis, which was followed by recovery in ∼50 μs. This difference band was absent in the spectra of the full-length protein, suggesting that the isolated sensor domain undergoes conformational changes of the protein backbone upon CO photolysis and that the changes are restrained by the signaling domain. The time-resolved difference spectrum at 200 μs exhibited a pattern similar to that of the static (reduced - CO) difference spectrum, although the peak intensities were much weaker. Thus, the rearrangements of the protein moiety toward the equilibrium ligand-free structure occur in a time range of hundreds of microseconds.

Quaternary structure controls ligand dynamics in soluble guanylate cyclase

Soluble guanylate cyclase (sGC) is the mammalian endogenous nitric oxide (NO) receptor. The mechanisms of activation and deactivation of this heterodimeric enzyme are unknown. For deciphering them, functional domains can be overexpressed. We have probed the dynamics of the diatomic ligands NO and CO within the isolated heme domain β(1)(190) of human sGC by piconanosecond absorption spectroscopy. After photo-excitation of nitrosylated sGC, only NO geminate rebinding occurs in 7.5 ps. In β(1)(190), both photo-dissociation of 5c-NO and photo-oxidation occur, contrary to sGC, followed by NO rebinding (7 ps) and back-reduction (230 ps and 2 ns). In full-length sGC, CO geminate rebinding to the heme does not occur. In contrast, CO geminately rebinds to β(1)(190) with fast multiphasic process (35, 171, and 18 ns). We measured the bimolecular association rates k(on) = 0.075 ± 0.01 × 10(6) M(-1) · S(-1) for sGC and 0.83 ± 0.1 × 10(6) M(-1) · S(-1) for β(1)(190). These different dynamics reflect conformational changes and less proximal constraints in the isolated heme domain with respect to the dimeric native sGC. We concluded that the α-subunit and the β(1)(191-619) domain exert structural strains on the heme domain. These strains are likely involved in the transmission of the energy and relaxation toward the activated state after Fe(2+)-His bond breaking. This also reveals the heme domain plasticity modulated by the associated domains and subunit.

Heme reactivity is uncoupled from quaternary structure in gel-encapsulated hemoglobin: a resonance Raman spectroscopic study

Encapsulation of hemoglobin (Hb) in silica gel preserves structure and function but greatly slows protein motion, thereby providing access to intermediates along the allosteric pathway that are inaccessible in solution. Resonance Raman (RR) spectroscopy with visible and ultraviolet laser excitation provides probes of heme reactivity and of key tertiary and quaternary contacts. These probes were monitored in gels after deoxygenation of oxyHb and after CO binding to deoxyHb, which initiate conformational change in the R-T and T-R directions, respectively. The spectra establish that quaternary structure change in the gel takes a week or more but that the evolution of heme reactivity, as monitored by the Fe-histidine stretching vibration, ν(FeHis), is completed within two days, and is therefore uncoupled from the quaternary structure. Within each quaternary structure, the evolving ν(FeHis) frequencies span the full range of values between those previously associated with the high- and low-affinity end states, R and T. This result supports the tertiary two-state (TTS) model, in which the Hb subunits can adopt high- and low-affinity tertiary structures, r and t, within each quaternary state. The spectra also reveal different tertiary pathways, involving the breaking and reformation of E and F interhelical contacts in the R-T direction but not the T-R direction. In the latter, tertiary motions are restricted by the T quaternary contacts.

Protein conformational changes of the oxidative stress sensor, SoxR, upon redox changes of the [2Fe-2S] cluster probed with ultraviolet resonance Raman spectroscopy

The [2Fe-2S] transcription factor, SoxR, a member of the MerR family, functions as a bacterial sensor of oxidative stress in Escherichia coli. SoxR is activated by reversible one-electron oxidation of the [2Fe-2S] cluster and enhances the production of various antioxidant proteins through the SoxRS regulon. Ultraviolet resonance Raman (UVRR) spectroscopic analysis of SoxR revealed conformational changes upon reduction of the [2Fe-2S] cluster in the absence and presence of promoter oligonucleotide. UVRR spectra reflected the environmental or structural changes of Trp following reduction. Notably, the environment around Trp91 contacting the [2Fe-2S] cluster was altered to become more hydrophilic, whereas that around Trp98 exhibited a small change to become more hydrophobic. In addition, changes in cation-π interactions between the [2Fe-2S] cluster and Trp91 were suggested. On the other hand, the environment around Tyr was barely affected by the [2Fe-2S] reduction. Binding of the promoter oligonucleotide triggered changes in Tyr located in the DNA-binding domain, but not Trp. Furthermore, conformational changes induced upon reduction of DNA-bound SoxR were not significantly different from those of DNA-free SoxR.

Design of an infrared laser pulse to control the multiphoton dissociation of the Fe-CO bond in CO-heme compounds

Optimal control theory is used to design a laser pulse for the multiphoton dissociation of the Fe-CO bond in the CO-heme compounds. The study uses a hexacoordinated iron-porphyrin-imidazole-CO complex in its ground electronic state as a model for CO liganded to the heme group. The potential energy and dipole moment surfaces for the interaction of the CO ligand with the heme group are calculated using density functional theory. Optimal control theory, combined with a time-dependent quantum dynamical treatment of the laser-molecule interaction, is then used to design a laser pulse capable of efficiently dissociating the CO-heme complex model. The genetic algorithm method is used within the mathematical framework of optimal control theory to perform the optimization process. This method provides good control over the parameters of the laser pulse, allowing optimized pulses with simple time and frequency structures to be designed. The dependence of photodissociation yield on the choice of initial vibrational state and of initial laser field parameters is also investigated. The current work uses a reduced dimensionality model in which only the Fe-C and C-O stretching coordinates are explicitly taken into account in the time-dependent quantum dynamical calculations. The limitations arising from this are discussed in Sec. IV.

Picosecond primary structural transition of the heme is retarded after nitric oxide binding to heme proteins

We investigated the ultrafast structural transitions of the heme induced by nitric oxide (NO) binding for several heme proteins by subpicosecond time-resolved resonance Raman and femtosecond transient absorption spectroscopy. We probed the heme iron motion by the evolution of the iron-histidine Raman band intensity after NO photolysis. Unexpectedly, we found that the heme response and iron motion do not follow the kinetics of NO rebinding. Whereas NO dissociation induces quasi-instantaneous iron motion and heme doming (<0.6 ps), the reverse process results in a much slower picosecond movement of the iron toward the planar heme configuration after NO binding. The time constant for this primary domed-to-planar heme transition varies among proteins (approximately 30 ps for myoglobin and its H64V mutant, approximately 15 ps for hemoglobin, approximately 7 ps for dehaloperoxidase, and approximately 6 ps for cytochrome c) and depends upon constraints exerted by the protein structure on the heme cofactor. This observed phenomenon constitutes the primary structural transition in heme proteins induced by NO binding.

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.

Protein collective motions coupled to ligand migration in myoglobin

Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.

Why does binding of proteins to DNA or proteins to proteins not necessarily spell function?

Studies of binding are often question: first, is the observed binding functional, and second, if it is, which function? Is it activation or repression? The first question relates to binding at different sites; the second relates to binding at similar sites. These questions apply to transcription factors binding to genomic DNA and to protein interaction domains binding to their partners. Here, we explain that both can be understood in terms of allostery and the cellular (or in vitro) environment. The idea is simple yet powerful; it emphasizes the role of allostery in defining whether binding between transcription factors and (cognate or noncognate) DNA sequences will lead to function and to the type of function. Allosteric effects are the outcome of dynamically shifting populations; thus binding to even slightly different DNA sequences will lead to different transcription factor conformations that can be reflected in the binding sites to their co-regulators. Currently, allostery is not considered when trying to understand how binding phenomena determine the functional outcome. Allosteric effects can enhance the binding specificity in a function-oriented manner. Here we provide a biological rationale that considers cellular crowding effects.