Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Interdomain Interactions Modulate the Active Site Dynamics of Human Inducible Nitric Oxide Synthase

Nitric oxide synthase (NOS) is a homodimeric flavohemoprotein responsible for catalyzing the oxidation of l-arginine (l-Arg) to citrulline and nitric oxide. Electrons are supplied for the reaction via interdomain electron transfer between an N-terminal heme-containing oxygenase domain and a FMN-containing (sub)domain of a C-terminal reductase domain. Extensive attention has focused on elucidating how conformational dynamics regulate electron transfer between the domains. Here we investigate the impact of the interdomain FMN-heme interaction on the heme active site dynamics of inducible NOS (iNOS). Steady state linear and time-resolved two-dimensional infrared (2D IR) spectroscopy was applied to probe a CO ligand at the heme within the oxygenase domain for full-length and truncated or mutated constructs of human iNOS. Whereas the linear IR spectra of the CO ligand were identical among the constructs, 2D IR spectroscopy revealed variation in the frequency dynamics. The wild-type constructs that can properly form the FMN/oxygenase docked state due to the presence of both the FMN and oxygenase domains showed slower dynamics than the oxygenase domain alone. Introduction of the mutation (E546N) predicted to perturb electrostatic interactions between the domains resulted in measured dynamics intermediate between those for the full-length and individual oxygenase domain, consistent with perturbation to the docked/undocked equilibrium. These results indicate that docking of the FMN domain to the oxygenase domain not only brings the FMN cofactor within electron transfer distance of the heme domain but also modulates the dynamics sensed by the CO ligand within the active site in a way expected to promote efficient electron transfer.

Analysis of the effect of 1-Allyl-3-Methylimidazolium chloride on thermodynamic stability, folding kinetics, and motional dynamics of horse cytochrome c

1-allyl-3-methylimidazolium chloride (AMIMCl) acts as a potential green solvent for proteins. The present work provides a possible pathway by which the structural, kinetic, thermodynamic, and folding properties of horse cytochrome c (cyt c) are affected in green aqueous-AMIMCl systems. Analysis of the effect of AMIMCl on thermodynamic stability, refolding/unfolding kinetics, and motional dynamics of cyt c provided important information, (i) AMIMCl decreases the thermodynamic stability of reduced cyt c and also strengthens the guanidinium chloride (GdmCl)-mediated decrease in thermodynamic stability of protein, (ii) AMIMCl reduces the thermal-fluctuation of Met80-containing omega-loop of natively-folded compact state of carbonmonoxycytochrome c (MCO-state) due to polyfunctional interactions between the AMIM⁺ and different groups of protein, (iii) AMIMCl shifts the kinetic chevron plot, ln k_obs[GdmCl] to the lower concentration of GdmCl, (iv) AMIMCl shifts the refolding and unfolding limps to vertically downwards and upwards, respectively, and (v) AMIMCl reducing the unfolding free energy estimated by both thermodynamic and kinetic analysis.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

Ruthenium hydrides encapsulated in sol–gel glasses exhibit new ultrafast vibrational dynamics

Vibrational dynamics were measured by IR pump–probe spectroscopy and two-dimensional IR spectroscopy for triruthenium dodecacarbonyl and the undecacarbonyl hydride that forms when it is encapsulated in an alumina sol–gel glass. For comparison, a triruthenium undecacarbonyl hydride salt was also synthesized and studied in neat solution to identify the potential influence of the confined solvent environment on the dynamics experienced by carbon monoxide ligands. The vibrational lifetime was found to be significantly decreased for both hydride species relative to the dodecacarbonyl compound. Conversely, spectral diffusion of the CO vibrations was measured to be faster for the parent compound. The most significant dynamic changes occurred upon transformation from the starting compound to the hydride, while only minor differences were observed between the dynamics of the freely dissolved and sol–gel encapsulated hydrides. The results suggest that the structural change to the hydride has the largest impact on the dynamics and that its improved catalytic properties likely do not originate from confined solvent effects.

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

Proteins function as ensembles of interconverting structures. The motions span from picosecond bond rotations to millisecond and longer subunit displacements. Characterization of functional dynamics on all spatial and temporal scales remains challenging experimentally. Two-dimensional infrared spectroscopy (2D IR) is maturing as a powerful approach for investigating proteins and their dynamics. We outline the advantages of IR spectroscopy, describe 2D IR and the information it provides, and introduce vibrational groups for protein analysis. We highlight example studies that illustrate the power and versatility of 2D IR for characterizing protein dynamics and conclude with a brief discussion of the outlook for biomolecular 2D IR.

Structural Insights into the Heme Pocket and Oligomeric State of Non-Symbiotic Hemoglobins from Arabidopsis thaliana

Non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana are hexacoordinate heme-proteins that likely have different biological roles, in view of diverse tissue localization, expression pattern, and ligand binding properties. Herein, we expand upon previous biophysical studies on these isoforms, focusing on their oligomeric states and circular dichroism (CD) characteristics. We found that AHb1 exists in solution in a concentration-dependent monomer-dimer equilibrium, while AHb2 is present only as a monomer. The quaternary structure of AHb1 affects its degree of hexacoordination with the formation of the dimer that enhances pentacoordination. Accordingly, the mutant of a conserved residue within the dimeric interface, AHb1-T45A, which is mostly monomeric in solution, has an equilibrium that is shifted toward a hexacoordinate form compared to the wild-type protein. CD studies further support differences in the globin's structure and heme moiety. The Soret CD spectra for AHb2 are opposite in sense to those for AHb1, reflecting different patterns of heme-protein side chain contacts in the two proteins. Moreover, the smaller contribution of the heme to the near-UV CD in AHb2 compared to AHb1 suggests a weaker heme-protein association in AHb2. Our data corroborate the structural diversity of AHb1 and AHb2 and confirm the leghemoglobin-like structural properties of AHb2.

Two-dimensional infrared spectroscopic study of cytochrome c peroxidase activity in deep eutectic solvent

Deep eutectic solvents (DESs) prepared by mixing hydrogen-bond donor and acceptor molecules have been found to be of use in several applications. Recently, it was shown that DESs can enhance the peroxidation activity of cytochrome c. Here, to elucidate the effects of DESs on the peroxidase activity of cytochrome c, we carried out linear and nonlinear infrared spectroscopic studies of the CO stretch mode of carbon monoxide cytochrome c (COCytc) in ethylammonium chloride (EAC)/urea DES. The FTIR spectrum of COCytc shows a significant spectral shift upon addition of the DES. The broadening and red-shifting of the CO band are observed in both urea and DES solutions, which are induced by the change of the distal ligands around the heme. Although the FTIR study is sensitive to structural changes in the active site, it does not provide quantitative information about structural dynamics related to the catalytic activity itself. Thus, we carried out two-dimensional IR spectroscopy of the CO mode, which suggests that there is a different conformer that could be related to the enhanced catalytic activity in DES. In particular, the spectral diffusion dynamics of that conformer exhibits quite different behavior. The experimental results lead us to propose a hypothesis that the DES increases the population of the conformer with distal ligand lysines close to the reaction center through the combining effect of urea and EAC, which results in the enhancement of the peroxidase activity of cytochrome c. We anticipate that the present experimental work stimulates future investigations of the effects of DES on biocatalysis.

Enhanced vibrational solvatochromism and spectral diffusion by electron rich substituents on small molecule silanes

Fourier transform infrared and two-dimensional IR (2D-IR) spectroscopies were applied to two different silanes in three different solvents. The selected solutes exhibit different degrees of vibrational solvatochromism for the Si-H vibration. Density functional theory calculations confirm that this difference in sensitivity is the result of higher mode polarization with more electron withdrawing ligands. This mode sensitivity also affects the extent of spectral diffusion experienced by the silane vibration, offering a potential route to simultaneously optimize the sensitivity of vibrational probes in both steady-state and time-resolved measurements. Frequency-frequency correlation functions obtained by 2D-IR show that both solutes experience dynamics on similar time scales and are consistent with a picture in which weakly interacting solvents produce faster, more homogeneous fluctuations. Molecular dynamics simulations confirm that the frequency-frequency correlation function obtained by 2D-IR is sensitive to the presence of hydrogen bonding dynamics in the surrounding solvation shell.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

Unraveling the dynamics and structure of functionalized self-assembled monolayers on gold using 2D IR spectroscopy and MD simulations

Functionalized self-assembled monolayers (SAMs) are the focus of ongoing investigations because they can be chemically tuned to control their structure and dynamics for a wide variety of applications, including electrochemistry, catalysis, and as models of biological interfaces. Here we combine reflection 2D infrared vibrational echo spectroscopy (R-2D IR) and molecular dynamics simulations to determine the relationship between the structures of functionalized alkanethiol SAMs on gold surfaces and their underlying molecular motions on timescales of tens to hundreds of picoseconds. We find that at higher head group density, the monolayers have more disorder in the alkyl chain packing and faster dynamics. The dynamics of alkanethiol SAMs on gold are much slower than the dynamics of alkylsiloxane SAMs on silica. Using the simulations, we assess how the different molecular motions of the alkyl chain monolayers give rise to the dynamics observed in the experiments.

Histidine Orientation Modulates the Structure and Dynamics of a de Novo Metalloenzyme Active Site

The ultrafast dynamics of a de novo metalloenzyme active site is monitored using two-dimensional infrared spectroscopy. The homotrimer of parallel, coiled coil α-helices contains a His3-Cu(I) metal site where CO is bound and serves as a vibrational probe of the hydrophobic interior of the self-assembled complex. The ultrafast spectral dynamics of Cu-CO reveals unprecedented ultrafast (2 ps) nonequilibrium structural rearrangements launched by vibrational excitation of CO. This initial rapid phase is followed by much slower ∼40 ps vibrational relaxation typical of metal-CO vibrations in natural proteins. To identify the hidden coupled coordinate, small molecule analogues and the full peptide were studied by QM and QM/MM calculations, respectively. The calculations show that variation of the histidines' dihedral angles in coordinating Cu controls the coupling between the CO stretch and the Cu-C-O bending coordinates. Analysis of different optimized structures with significantly different electrostatic field magnitudes at the CO ligand site indicates that the origin of the stretch-bend coupling is not directly due to through-space electrostatics. Instead, the large, ∼3.6 D dipole moments of the histidine side chains effectively transduce the electrostatic environment to the local metal coordination orientation. The sensitivity of the first coordination sphere to the protein electrostatics and its role in altering the potential energy surface of the bound ligands suggests that long-range electrostatics can be leveraged to fine-tune function through enzyme design.

Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions representing radiation-matter interactions. Here we apply this method to an anharmonic chromophore coupled to a harmonic bath. A forward-backward trajectory implementation of the OMT method is described that addresses the numerical challenges of applying the OMT to large systems with disparate frequency scales. The OMT is shown to well reproduce line shapes and waiting time dynamics in the pure dephasing limit of weak coupling to an off-resonant bath. The OMT is also shown to describe a case where energy transfer is the predominant source of line broadening.

Ultrafast Dynamics of Carboxy-Hemoglobin: Two-Dimensional Infrared Spectroscopy Experiments and Simulations

This Letter presents a comparison between experimental and simulated 2D mid-infrared spectra of carboxy-hemoglobin in the spectral region of the carbon monoxide stretching mode. The simulations rely on a fluctuating potential energy surface that includes both the effect of heme and the protein surroundings computed from molecular dynamics simulations. A very good agreement between theory and experiment is obtained with no adjustable parameters. The simulations show that the effect of the distal histidine through the hydrogen bond is strong and is directly responsible for the slow decay of the frequency-frequency correlation function on a 10 ps time scale. This study confirms that fluctuations in carboxy-hemoglobin are more inhomogeneous than those in the more frequently studied carboxy-myoglobin. The comparison between simulations and experiments brings valuable information on the complex relation between protein structure and spectral diffusion.

Vibrational ladder climbing in carboxy-hemoglobin: effects of the protein environment

We present simulations on vibrational ladder climbing in carboxy-hemoglobin. Motivated by recent experiments, we study the influence of different realistic pump probe parameters. To allow for a direct comparison with experimental results, transient absorption spectra obtained by a weak probe pulse following the strong, shaped pump pulse are calculated. The influence of the protein fluctuations is taken into account using a recently developed microscopic model. This model consists of a quantum Hamiltonian describing the CO vibration in carboxy-hemoglobin, together with a fluctuating potential, which is obtained by electronic structure calculation based on a large number of protein configurations. Using realistic pulse parameters, vibrational excitations to very high-lying states are possible, in qualitative agreement with experimental observations.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

Two-dimensional spectroscopy of coupled vibrations with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical representation of the nonlinear vibrational response function used to compute multidimensional infrared spectra. In this method, response functions are calculated from a sequence of classical trajectories linked by discontinuities representing the effects of radiation-matter interactions, thus providing an approximation to quantum dynamics using classical inputs. This approach was previously formulated and assessed numerically for a single anharmonic degree of freedom. Our previous work is generalized here in two respects. First, the derivation of the OMT is extended to any number of coupled anharmonic vibrations by determining semiclassical approximations for pairs of double-sided Feynman diagrams. Second, an efficient numerical procedure is developed for calculating two-dimensional infrared spectra of coupled anharmonic vibrations in the OMT approximation. The OMT approximation is shown to reproduce the fundamental features of the quantum response function including both coherence and population dynamics.

pH-dependent picosecond structural dynamics in the distal pocket of nitrophorin 4 investigated by 2D IR spectroscopy

Nitrophorin 4 (NP4) belongs to a family of pH-sensitive, nitric oxide (NO) transporter proteins that undergo a large structural change from a closed to an open conformation at high pH to allow for NO delivery. Measuring the pH-dependent structural dynamics in NP4-NO around the ligand binding site is crucial for developing a mechanistic understanding of NO binding and release. In this study, we use coherent two-dimensional infrared (2D IR) spectroscopy to measure picosecond structural dynamics sampled by the nitrosyl stretch in NP4-NO as a function of pH at room temperature. Our results show that both the closed and open conformers of the protein are present at low (pD 5.1) and high (pD 7.9) pH conditions. The closed and open conformers are characterized by two frequencies of the nitrosyl stretching vibration labeled A0 and A1, respectively. Analysis of the 2D IR line shapes reveals that at pD 5.1, the closed conformer experiences structural fluctuations arising from solvation dynamics on a ∼3 ps time scale. At pD 7.9, both the open and closed conformers exhibit fluctuations on a ∼1 ps time scale. At both pD conditions, the closed conformers maintain a static distribution of structures within the experimental time window of 100 ps. This is in contrast to the open conformer, which is able to interconvert among its substates on a ∼100 ps time scale. Our results directly measure the time scales of solvation dynamics in the distal pocket, the flexibility of the open conformation at high pH, and the rigidity of the closed conformers at both pH conditions. We discuss how the pH-dependent equilibrium structural fluctuations of the nitrosyl ligand measured in this study are related to the uptake and delivery of nitric oxide in NP4.

Computational infrared and two-dimensional infrared photon echo spectroscopy of both wild-type and double mutant myoglobin-CO proteins

The CO stretching mode of both wild-type and double mutant ( T67R / S92D ) MbCO (carbonmonoxymyoglobin) proteins is an ideal infrared (IR) probe for studying the local electrostatic environment inside the myoglobin heme pocket. Recently, to elucidate the conformational switching dynamics between two distinguishable states, extensive IR absorption, IR pump-probe, and two-dimensional (2D) IR spectroscopic studies for various mutant MbCO's have been performed by the Fayer group. They showed that the 2D IR spectroscopy of the double mutant, which has a peroxidase enzyme activity, reveals a rapid chemical exchange between two distinct states, whereas that of the wild-type does not. Despite the fact that a few simulation studies on these systems were already performed and reported, such complicated experimental results have not been fully reproduced nor described in terms of conformational state-to-state transition processes. Here, we first develop a distributed vibrational solvatochromic charge model for describing the CO stretch frequency shift reflecting local electric potential changes. Then, by carrying out molecular dynamic simulations of the two MbCO's and examining their CO frequency trajectories, it becomes possible to identify a proper reaction coordinate consisting of His64 imidazole ring rotation and its distance to the CO ligand. From the 2D surfaces of the resulting potential of mean forces, the spectroscopically distinguished A1 and A3 states of the wild-type as well as two more substates of the double mutant are identified and their vibrational frequencies and distributions are separately examined. Our simulated IR absorption and 2D IR spectra of the two MbCO's are directly compared with the previous experimental results reported by the Fayer group. The chemical exchange rate constants extracted from the two-state kinetic analyses of the simulated 2D IR spectra are in excellent agreement with the experimental values. On the basis of the quantitative agreement between the simulated spectra and experimental ones, we further examine the conformational differences in the heme pockets of the two proteins and show that the double mutation, T67R / S92D , suppresses the A1 population, restricts the imidazole ring rotation, and increases hydrogen-bond strength between the imidazole Nε-H and the oxygen atom of the CO ligand. It is believed that such delicate change of distal His64 imidazole ring dynamics induced by the double mutation may be responsible for its enhanced peroxidase catalytic activity as compared to the wild-type myoglobin.

Measuring electrostatic fields in both hydrogen-bonding and non-hydrogen-bonding environments using carbonyl vibrational probes

Vibrational probes can provide a direct readout of the local electrostatic field in complex molecular environments, such as protein binding sites and enzyme active sites. This information provides an experimental method to explore the underlying physical causes of important biomolecular processes such as binding and catalysis. However, specific chemical interactions such as hydrogen bonds can have complicated effects on vibrational probes and confound simple electrostatic interpretations of their frequency shifts. We employ vibrational Stark spectroscopy along with infrared spectroscopy of carbonyl probes in different solvent environments and in ribonuclease S to understand the sensitivity of carbonyl frequencies to electrostatic fields, including those due to hydrogen bonds. Additionally, we carried out molecular dynamics simulations to calculate ensemble-averaged electric fields in solvents and in ribonuclease S and found excellent correlation between calculated fields and vibrational frequencies. These data enabled the construction of a robust field-frequency calibration curve for the C═O vibration. The present results suggest that carbonyl probes are capable of quantitatively assessing the electrostatics of hydrogen bonding, making them promising for future study of protein function.

An optimized semiclassical approximation for vibrational response functions

The observables of multidimensional infrared spectroscopy may be calculated from nonlinear vibrational response functions. Fully quantum dynamical calculations of vibrational response functions are generally impractical, while completely classical calculations are qualitatively incorrect at long times. These challenges motivate the development of semiclassical approximations to quantum mechanics, which use classical mechanical information to reconstruct quantum effects. The mean-trajectory (MT) approximation is a semiclassical approach to quantum vibrational response functions employing classical trajectories linked by deterministic transitions representing the effects of the radiation-matter interaction. Previous application of the MT approximation to the third-order response function R(3)(t3, t2, t1) demonstrated that the method quantitatively describes the coherence dynamics of the t3 and t1 evolution times, but is qualitatively incorrect for the waiting-time t2 period. Here we develop an optimized version of the MT approximation by elucidating the connection between this semiclassical approach and the double-sided Feynman diagrams (2FD) that represent the quantum response. Establishing the direct connection between 2FD and semiclassical paths motivates a systematic derivation of an optimized MT approximation (OMT). The OMT uses classical mechanical inputs to accurately reproduce quantum dynamics associated with all three propagation times of the third-order vibrational response function.

Effect of solvent polarity on the vibrational dephasing dynamics of the nitrosyl stretch in an Fe(II) complex revealed by 2D IR spectroscopy

The vibrational dephasing dynamics of the nitrosyl stretching vibration (ν(NO)) in sodium nitroprusside (SNP, Na2[Fe(CN)5NO]·2H2O) are investigated using two-dimensional infrared (2D IR) spectroscopy. The ν(NO) in SNP acts as a model system for the nitrosyl ligand found in metalloproteins which play an important role in the transportation and detection of nitric oxide (NO) in biological systems. We perform a 2D IR line shape study of the ν(NO) in the following solvents: water, deuterium oxide, methanol, ethanol, ethylene glycol, formamide, and dimethyl sulfoxide. The frequency of the ν(NO) exhibits a large vibrational solvatochromic shift of 52 cm(-1), ranging from 1884 cm(-1) in dimethyl sulfoxide to 1936 cm(-1) in water. The vibrational anharmonicity of the ν(NO) varies from 21 to 28 cm(-1) in the solvents used in this study. The frequency-frequency correlation functions (FFCFs) of the ν(NO) in SNP in each of the seven solvents are obtained by fitting the experimentally obtained 2D IR spectra using nonlinear response theory. The fits to the 2D IR line shape reveal that the spectral diffusion time scale of the ν(NO) in SNP varies from 0.8 to 4 ps and is negatively correlated with the empirical solvent polarity scales. We compare our results with the experimentally determined FFCFs of other charged vibrational probes in polar solvents and in the active sites of heme proteins. Our results suggest that the vibrational dephasing dynamics of the ν(NO) in SNP reflect the fluctuations of the nonhomogeneous electric field created by the polar solvents around the nitrosyl and cyanide ligands. The solute solvent interactions occurring at the trans-CN ligand are sensed through the π-back-bonding network along the Fe-NO bond in SNP.

Dynamics in the interior of AOT lamellae investigated with two-dimensional infrared spectroscopy

The dynamics inside the organic regions of aerosol-OT (AOT)/water mixtures in the lamellar mesophase, bicontinuous cubic (BC) phase, and in an analogous molecule without the charged sulfonate headgroup are investigated by observing spectral diffusion, orientational relaxation and population relaxation using ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy and IR pump-probe experiments on the asymmetric CO stretch of a vibrational probe, tungsten hexacarbonyl (W(CO)6). The water layer thickness between the bilayer planes in the lamellar phase was varied. For comparison, the dynamics of W(CO)6 in the normal liquid bis(2-ethylhexyl) succinate (EHS), which is analogous to AOT but has no charged sulfonate headgroup, were also studied. The 2D IR experiments measure spectral diffusion, which results from the structural evolution of the system. Spectral diffusion is quantified by the frequency-frequency correlation function (FFCF). In addition to a homogeneous component, the FFCFs are biexponential decays with fast and slow time components of ∼12.5 and ∼150 ps in the lamellar phase. Both components of the FFCF are independent of the number of water molecules per headgroup for the lamellae, but they slow somewhat in the BC phase. The dynamics in the ordered phases are in sharp contrast to the dynamics in EHS, which displays fast and slow components of the FFCF of 5 and 80 ps, respectively. As the hydration level of AOT increases, vibrational lifetime decreases, suggesting some change in the local environment of W(CO)6 with water content.

Protein dynamics studied with ultrafast two-dimensional infrared vibrational echo spectroscopy

Proteins, enzymes, and other biological molecules undergo structural dynamics as an intrinsic part of their biological functions. While many biological processes occur on the millisecond, second, and even longer time scales, the fundamental structural dynamics that eventually give rise to such processes occur on much faster time scales. Many decades ago, chemical kineticists focused on the inverse of the reaction rate constant as the important time scale for a chemical reaction. However, through transition state theory and a vast amount of experimental evidence, we now know that the key events in a chemical reaction can involve structural fluctuations that take a system of reactants to its transition state, the crossing of a barrier, and the eventual relaxation to product states. Such dynamics occur on very fast time scales. Today researchers would like to investigate the fast structural fluctuations of biological molecules to gain an understanding of how biological processes proceed from simple structural changes in biomolecules to the final, complex biological function. The study of the fast structural dynamics of biological molecules requires experiments that operate on the appropriate time scales, and in this Account, we discuss the application of ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy to the study of protein dynamics. The 2D IR vibrational echo experiment is akin to 2D NMR, but it operates on time scales many orders of magnitude faster. In the experiments, a particular vibrational oscillator serves as a vibrational dynamics probe. As the structure of the protein evolves in time, the structural changes are manifested as time-dependent changes in the frequency of the vibrational dynamics probe. The 2D IR vibrational echo experiments can track the vibrational frequency evolution, which we then relate to the time evolution of the protein structure. In particular, we measured protein substate interconversion for mutants of myoglobin using 2D IR chemical exchange spectroscopy and observed well-defined substate interconversion on a sub-100 ps time scale. In another study, we investigated the influence of binding five different substrates to the enzyme cytochrome P450(cam). The various substrates affect the enzyme dynamics differently, and the observed dynamics are correlated with the enzyme's selectivity of hydroxylation of the substrates and with the substrate binding affinity.

Fast dynamics of HP35 for folded and urea-unfolded conditions

The changes in fast dynamics of HP35 with a double CN vibrational dynamics label (HP35-P(2)) as a function of the extent of denaturation by urea were investigated with two-dimensional infrared (2D IR) vibrational echo spectroscopy. Cyanophenylalanine (PheCN) replaces the native phenylalanine at two residues in the hydrophobic core of HP35, providing vibrational probes. NMR data show that HP35-P(2) maintains the native folded structure similar to wild type and that both PheCN residues share essentially the same environment within the peptide. A series of time-dependent 2D IR vibrational echo spectra were obtained for the folded peptide and the increasingly unfolded peptide. Analysis of the time dependence of the 2D spectra yields the system's spectral diffusion, which is caused by the sampling of accessible structures of the peptide under thermal equilibrium conditions. The structural dynamics become faster as the degree of unfolding is increased.

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

A solvatochromic model calibrates nitriles' vibrational frequencies to electrostatic fields

Electrostatic interactions provide a primary connection between a protein's three-dimensional structure and its function. Infrared probes are useful because vibrational frequencies of certain chemical groups, such as nitriles, are linearly sensitive to local electrostatic field and can serve as a molecular electric field meter. IR spectroscopy has been used to study electrostatic changes or fluctuations in proteins, but measured peak frequencies have not been previously mapped to total electric fields, because of the absence of a field-frequency calibration and the complication of local chemical effects such as H-bonds. We report a solvatochromic model that provides a means to assess the H-bonding status of aromatic nitrile vibrational probes and calibrates their vibrational frequencies to electrostatic field. The analysis involves correlations between the nitrile's IR frequency and its (13)C chemical shift, whose observation is facilitated by a robust method for introducing isotopes into aromatic nitriles. The method is tested on the model protein ribonuclease S (RNase S) containing a labeled p-CN-Phe near the active site. Comparison of the measurements in RNase S against solvatochromic data gives an estimate of the average total electrostatic field at this location. The value determined agrees quantitatively with molecular dynamics simulations, suggesting broader potential for the use of IR probes in the study of protein electrostatics.

Ribonuclease S dynamics measured using a nitrile label with 2D IR vibrational echo spectroscopy

A nitrile-labeled amino acid, p-cyanophenylalanine, is introduced near the active site of the semisynthetic enzyme ribonuclease S to serve as a probe of protein dynamics and fluctuations. Ribonuclease S is the limited proteolysis product of subtilisin acting on ribonuclease A, and consists of a small fragment including amino acids 1-20, the S-peptide, and a larger fragment including residues 21-124, the S-protein. A series of two-dimensional vibrational echo experiments performed on the nitrile-labeled S-peptide and the RNase S are described. The time-dependent changes in the two-dimensional infrared vibrational echo line shapes are analyzed using the center line slope method to obtain the frequency-frequency correlation function (FFCF). The observations show that the nitrile probe in the S-peptide has dynamics that are similar to, but faster than, those of the single amino acid p-cyanophenylalanine in water. In contrast, the dynamics of the nitrile label when the peptide is bound to form ribonuclease S are dominated by homogeneous dephasing (motionally narrowed) contributions with only a small contribution from very fast inhomogeneous structural dynamics. The results provide insights into the nature of the structural dynamics of the ribonuclease S complex. The equilibrium dynamics of the nitrile labeled S-peptide and the ribonuclease S complex are also investigated by molecular dynamics simulations. The experimentally determined FFCFs are compared to the FFCFs obtained from the molecular dynamics simulations, thereby testing the capacity of simulations to determine the amplitudes and time scales of protein structural fluctuations on fast time scales under thermal equilibrium conditions.

Two-dimensional IR spectroscopy of protein dynamics using two vibrational labels: a site-specific genetically encoded unnatural amino acid and an active site ligand

Protein dynamics and interactions in myoglobin (Mb) were characterized via two vibrational dynamics labels (VDLs): a genetically incorporated site-specific azide (Az) bearing unnatural amino acid (AzPhe43) and an active site CO ligand. The Az-labeled protein was studied using ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy. CO bound at the active site of the heme serves as a second VDL located nearby. Therefore, it was possible to use Fourier transform infrared (FT-IR) and 2D IR spectroscopic experiments on the Az in unligated Mb and in Mb bound to CO (MbAzCO) and on the CO in MbCO and MbAzCO to investigate the environment and motions of different states of one protein from the perspective of two spectrally resolved VDLs. A very broad bandwidth 2D IR spectrum, encompassing both the Az and CO spectral regions, found no evidence of direct coupling between the two VDLs. In MbAzCO, both VDLs reported similar time scale motions: very fast homogeneous dynamics, fast, ∼1 ps dynamics, and dynamics on a much slower time scale. Therefore, each VDL reports independently on the protein dynamics and interactions, and the measured dynamics are reflective of the protein motions rather than intrinsic to the chemical nature of the VDL. The AzPhe VDL also permitted study of oxidized Mb dynamics, which could not be accessed previously with 2D IR spectroscopy. The experiments demonstrate that the combined application of 2D IR spectroscopy and site-specific incorporation of VDLs can provide information on dynamics, structure, and interactions at virtually any site throughout any protein.

Oxygen binding to Arabidopsis thaliana AHb2 nonsymbiotic hemoglobin: evidence for a role in oxygen transport

Nonsymbiotic hemoglobins AHb1 and AHb2 discovered in Arabidopsis thaliana are likely to carry out distinct physiological roles, in consideration of their differences in sequence, structure, expression pattern, and tissue localization. Despite a relatively fast autoxidation in the presence of O(2) , we were able to collect O(2) -binding curves for AHb2 in the presence of a reduction enzymatic system. AHb2 binds O(2) noncooperatively with a p50 of 0.021 ± 0.003 Torr, a value consistent with a recently proposed role in O(2) transport. The analysis of the internal cavities derived from the structures sampled in molecular dynamics simulations confirms strong differences with AHb1, proposed to work as a NO deoxygenase in vivo. Overall, our results are consistent with a role for AHb2 as an oxygen carrier, as recently proposed on the basis of experiments on AHb2-overexpressing mutants of A. thaliana.

Influence of histidine tag attachment on picosecond protein dynamics

Polyhistidine affinity tags are routinely employed as a convenient means of purifying recombinantly expressed proteins. A tacit assumption is commonly made that His tags have little influence on protein structure and function. Attachment of a His tag to the N-terminus of the robust globular protein myoglobin leads to only minor changes to the electrostatic environment of the heme pocket, as evinced by the nearly unchanged Fourier transform infrared spectrum of CO bound to the heme of His-tagged myoglobin. Experiments employing two-dimensional infrared vibrational echo spectroscopy of the heme-bound CO, however, find that significant changes occur to the short time scale (picoseconds) dynamics of myoglobin as a result of His tag incorporation. The His tag mainly reduces the dynamics on the 1.4 ps time scale and also alters protein motions of myoglobin on the slower, >100 ps time scale, as demonstrated by the His tag's influence on the fluctuations of the CO vibrational frequency, which reports on protein structural dynamics. The results suggest that affinity tags may have effects on protein function and indicate that investigators of affinity-tagged proteins should take this into consideration when investigating the dynamics and other properties of such proteins.

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.

Dynamics of the folded and unfolded villin headpiece (HP35) measured with ultrafast 2D IR vibrational echo spectroscopy

A series of two-dimensional infrared vibrational echo experiments performed on nitrile-labeled villin headpiece [HP35-(CN)(2)] is described. HP35 is a small peptide composed of three alpha helices in the folded configuration. The dynamics of the folded HP35-(CN)(2) are compared to that of the guanidine-induced unfolded peptide, as well as the nitrile-functionalized phenylalanine (PheCN), which is used to differentiate the peptide dynamic contributions to the observables from those of the water solvent. Because the viscosity of solvent has a significant effect on fast dynamics, the viscosity of the solvent is held constant by adding glycerol. For the folded peptide, the addition of glycerol to the water solvent causes observable slowing of the peptide's dynamics. Holding the viscosity constant as GuHCl is added, the dynamics of unfolded peptide are much faster than those of the folded peptide, and they are very similar to that of PheCN. These observations indicate that the local environment of the nitrile in the unfolded peptide resembles that of PheCN, and the dynamics probed by the CN are dominated by the fluctuations of the solvent molecules, in contrast to the observations on the folded peptide.

Relationship between protein structural fluctuations and rebinding dynamics in ferric haem nitrosyls

The interaction of nitric oxide (NO) with haem proteins is widespread in biology. In the current paper, we present the first ultrafast 2D-IR (two-dimensional infrared) spectroscopic analysis of haem nitrosylation, which has been combined with time-resolved IR pump–probe studies to investigate the relationship between equilibrium vibrational dynamics of the haem environment and ligand rebinding behaviour following photolysis of NO from the Fe(III)–NO site. Studies of two haem proteins, Mb (myoglobin) and Cc (cytochrome c), which play different physiological roles, reveal marked contrasts in the ultrafast fluctuations of the protein pockets containing the haem, showing that the Mb pocket is somewhat more flexible than that of Cc. This correlates strongly with slower observed photolysis rebinding kinetics of Mb–NO compared with Cc–NO, and indicates a direct link between ultrafast fluctuations and biological functionality. Furthermore, this indicates the validity of linear response theories in relation to protein ligand binding. Finally, 2D-IR shows that Cc–NO displays two distinct structural sub-sites at room temperature that do not exchange on the timescales accessible via the NO vibrational lifetime.