Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Abstract

Significance In this multidisciplinary study, we investigated local and global structural responses upon carbon monoxide (CO) dissociation from the heme iron(II) in myoglobin in solution using ultrafast X-ray absorption and scattering. CO dissociation was optically triggered and pulsed X-rays probed structural and electronic changes as a function of time, revealing an iron(II) spin state transition and nuclear motions propagating from the active site throughout the entire protein. Quantum mechanical computational results were used to assign the iron K-edge absorption features, and to map out electronic and nuclear trajectories during Fe–CO bond elongation and heme macrocycle doming. The study has important implications for heme–protein function and the interaction of small molecule substrates with the metal active centers in enzymes. Ultrafast structural dynamics with different spatial and temporal scales were investigated during photodissociation of carbon monoxide (CO) from iron(II)–heme in bovine myoglobin during the first 3 ps following laser excitation. We used simultaneous X-ray transient absorption (XTA) spectroscopy and X-ray transient solution scattering (XSS) at an X-ray free electron laser source with a time resolution of 80 fs. Kinetic traces at different characteristic X-ray energies were collected to give a global picture of the multistep pathway in the photodissociation of CO from heme. In order to extract the reaction coordinates along different directions of the CO departure, XTA data were collected with parallel and perpendicular relative polarizations of the laser pump and X-ray probe pulse to isolate the contributions of electronic spin state transition, bond breaking, and heme macrocycle nuclear relaxation. The time evolution of the iron K-edge X-ray absorption near edge structure (XANES) features along the two major photochemical reaction coordinates, i.e., the iron(II)–CO bond elongation and the heme macrocycle doming relaxation were modeled by time-dependent density functional theory calculations. Combined results from the experiments and computations reveal insight into interplays between the nuclear and electronic structural dynamics along the CO photodissociation trajectory. Time-resolved small-angle X-ray scattering data during the same process are also simultaneously collected, which show that the local CO dissociation causes a protein quake propagating on different spatial and temporal scales. These studies are important for understanding gas transport and protein deligation processes and shed light on the interplay of active site conformational changes and large-scale protein reorganization.