Photochemistry and transient intermediates in a bacteriophytochrome photocycle revealed by multiscale simulations

09 March 2022, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Phytochromes are photoreceptors responsible for sensing light in plants, fungi and bacteria. Their photoactivation is initiated by the photoisomerization of an embedded chromophore, which triggers a large conformational change in the structure of the entire protein. Although phytochromes have been subject of numerous studies, the photoisomerization mechanism and the following reaction path leading to the final active state remain elusive. Here, we use an integrated computational approach that combines non-adiabatic surface hopping and adiabatic ground-state molecular dynamics simulations to gain atomistic details on the photoactivation mechanism of Deinococcus radiodurans bacteriophytochrome. Our simulations show that the ps-scale photoisomerization of the chromophore proceeds through a hula-twist mechanism that forces a counterclockwise rotation of the D-ring. The initial photoproduct rapidly evolves in an early intermediate which we characterize through IR spectroscopy simulation. The early intermediate then evolves on the nanosecond-to-microsecond scale to a late intermediate, characterized by a more disordered binding pocket and a clear weakening of the aspartate-to-arginine salt bridge interaction, whose cleavage is essential to interconvert to the final active state.

Keywords

Molecular Dynamics
phytochromes
Deinococcus radiodurans
QM/MM MD simulations
IR spectroscopy simulations
Surface Hopping
photosensory proteins
photoswitch families
photoreceptors
biliverdin
hula-twist

Supplementary materials

Title
Description
Actions
Title
Supplementary Information
Description
Selection of the QM semiempirical method; Computational details on the set-up of the non-adiabatic/adiabatic molecular dynamics simulations; Clustering details; Computational details on the IR spectroscopy simulation; Supplementary Figures.
Actions
Title
Supplementary Movie 2
Description
Each dot represents a trajectory in the conformational space defined by the dihedrals D5 and D6. In this video only some representative trajectories are represented. The left panel represents the first electronic excited state of the biliverdin chromophore, while the right panel the ground state. Blue dots are the reactive trajectories, while the orange dots are the non-reactive trajectories.
Actions
Title
Supplementary Movie 1
Description
Hula-Twist mechanism. We highlighted the interactions with Tyr263, His290 and a water molecule.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.