A fluctuating quantum model of the CO vibration in carboxyhemoglobin

Quasi-Static Two-Dimensional Infrared Spectra of the Carboxyhemoglobin Subsystem under Electric Fields: A Theoretical Study

There is no doubt that electric fields of a specific frequency and intensity could excite certain vibrational modes of a macromolecule, which alters its mode coupling and conformation. Motivated by recent experiments and theories, we study the mode coupling between the Fe-CO mode and CO-stretch mode and vibration energy transfer among the active site and proteins in carboxyhemoglobin (HbCO) under different electric fields using the quasi-static two-dimensional infrared spectra. This study uses iron-porphyrin-imidazole-CO and two distal histidines in HbCO as the subsystem. The potential energy and dipole moment surfaces of the subsystem are calculated using an all-electron ab initio (B3LYP-D3(BJ)) method with the basis set Lanl2dz for the Fe atom and 6-31G(d,p) for C, H, O, and N atoms. Although the subsystem is reduced dimensionally, the anharmonic frequency and anharmonicity of the CO-stretch mode show excellent agreement with experimental values. We use the revealing noncovalent interaction method to confirm the hydrogen bond between the H_ε atom of the His63 and the CO molecule. Our study confirms that the mode coupling between the Fe-CO mode and CO-stretch mode does not exist when the subsystem is free of electric field perturbation, which is coupled when the electric field is -0.5142 V/nm. In addition, with the increases of distance between the active site and the His92, there is no vibrational energy transfer between them when the electric field is 1.028 V/nm. We believe that our work could provide new ideas for increasing the dissociation efficiency of the Fe-CO bond and theoretical references for experimental research.

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.

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.

Molecular mechanics force field-based general map for the solvation effect on amide I probe of peptide in different micro-environments

A general electrostatic potential map based on molecular mechanics force field for modeling the amide I frequency is presented. This map is applied to N-methylacetamide (NMA) and designed to be transferable in different micro-environments. The electrostatic potentials from solvent and peptide side chain are projected on the amide unit of NMA to induce the frequency shift of amide I mode. It is shown that the predicted amide I frequency reproduces the experimental data satisfactorily, especially when NMA in polar solvents. The amide I frequency shift is largely determined by the solvents in aqueous solution while it is dominated by the local structure of peptide in other solvent environments. The map parameters are further applied on NMA-MeOH system and the obtained IR spectra show doublet peak profile with negligible deviation from the experimental data, suggesting the usefulness of this general map for providing information about vibrational parameters of amide motions of peptide in different environments.