Fully Quantum Chemical Treatment of Chromophore–Protein Interactions in Phytochromes

Hydrogen Bonding and Noncovalent Electric Field Effects in the Photoconversion of a Phytochrome

A profound understanding of protein structure and mechanism requires dedicated experimental and theoretical tools to elucidate electrostatic and hydrogen bonding interactions in proteins. In this work, we employed an approach to disentangle noncovalent and hydrogen-bonding electric field changes during the reaction cascade of a multidomain protein, i.e., the phytochrome Agp2. The approach exploits the spectroscopic properties of nitrile probes commonly used as reporter groups of the vibrational Stark effect. These probes were introduced into the protein through site-specific incorporation of noncanonical amino acids resulting in four variants with different positions and orientations of the nitrile groups. All substitutions left structures and the reaction mechanism unchanged. Structural models of the dark states (Pfr) were used to evaluate the total electric field at the nitrile label and its transition dipole moment. These quantities served as an internal standard to calculate the respective properties of the photoinduced products (Lumi-F, Meta-F, and Pr) based on the relative intensities of the nitrile stretching bands. In most cases, the spectral analysis revealed two substates with a nitrile in a hydrogen-bonded or hydrophobic environment. Using frequencies and intensities, we managed to extract the noncovalent contribution of the electric field from the individual substates. This analysis resulted in profiles of the noncovalent and hydrogen-bond-related electric fields during the photoinduced reaction cascade of Agp2. These profiles, which vary significantly among the four variants due to the different positions and orientations of the nitrile probes, were discussed in the context of the molecular events along the Pfr → Pr reaction cascade.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method

Fast parameterized methods such as density-functional tight-binding (DFTB) facilitate realistic calculations of large molecular systems, which can be accelerated by the fragment molecular orbital (FMO) method. Fragmentation facilitates interaction analyses between functional parts of molecular systems. In addition to DFTB, other parameterized methods combined with FMO are also described. Applications of FMO methods to biochemical and inorganic systems are reviewed.

The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores

The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.

Fragment molecular orbital calculations for biomolecules

Exploring biomolecule behavior, such as proteins and nucleic acids, using quantum mechanical theory can identify many life science phenomena from first principles. Fragment molecular orbital (FMO) calculations of whole single particles of biomolecules can determine the electronic state of the interior and surface of molecules and explore molecular recognition mechanisms based on intermolecular and intramolecular interactions. In this review, we summarized the current state of FMO calculations in drug discovery, virology, and structural biology, as well as recent developments from data science.

Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2

Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that bind a linear tetrapyrrole as a chromophore. They show photochromicity by having two stable states that can be interconverted by the photoisomerization of the chromophore. These photochemical properties make them an attractive target for biotechnological applications. However, their application is impeded by structural heterogeneity that reduces the yield of the photoconversion. The heterogeneity can originate either from the chromophore structure or the protein environment. Here, we study the origin of the heterogeneity in AnPixJg2, a representative member of the red/green cyanobacteriochrome family, that has a red absorbing parental state and a green absorbing photoproduct state. Using molecular dynamics simulations and umbrella sampling we have identified the protonation state of a conserved histidine residue as a trigger for structural heterogeneity. When the histidine is in a neutral form, the chromophore structure is homogenous, while in a positively charged form, the chromophore is heterogeneous with two different conformations. We have identified a correlation between the protonation of the histidine and the structural heterogeneity of the chromophore by detailed characterization of the interactions in the protein binding site. Our findings reconcile seemingly contradicting spectroscopic studies that attribute the heterogeneity to different sources. Furthermore, we predict that circular dichroism can be used as a diagnostic tool to distinguish different substates.

Influence of the Environment on Shaping the Absorption of Monomeric Infrared Fluorescent Proteins

Infrared fluorescent proteins (iRFPs) are potential candidates for deep-tissue in vivo imaging. Here, we provide molecular-level insights into the role of the protein environment in the structural stability of the chromophore within the protein binding pocket through the flexible hydrogen-bonding network using molecular dynamics simulation. Furthermore, we present systematic excited-state analysis to characterize the nature of the first two excited states and the role of the environment in shaping the nature of the chromophore's excited states within the hybrid quantum mechanics/molecular mechanics framework. Our results reveal that the environment red-shifts the absorption of the chromophore by about 0.32 eV compared to the isolated counterpart, and besides the structural stability, the protein environment does not alter the nature of the excited state of the chromophore significantly. Our study contributes to the fundamental understanding of the excited-state processes of iRFPs in a complex environment and provides a design principle for developing iRFPs with desired spectral properties.

The structural changes in the signaling mechanism of bacteriophytochromes in solution revealed by a multiscale computational investigation

Phytochromes are red-light sensing proteins, with important light-regulatory roles in different organisms, which are capturing an increasing interest in bioimaging and optogenetics. Upon absorption of light by the embedded bilin chromophore, they undergo structural changes that extend from the chromophore to the protein and finally drive the biological function. Up to now, the underlying mechanism still has to be characterized fully. Here we investigate the Pfr activated form of a bacterial phytochrome, by combining extensive molecular dynamics simulations with a polarizable QM/MM description of the spectroscopic properties, revealing a large structure relaxation in solution, compared to the crystal structure, both in the chromophore-binding pocket and in the overall structure of the phytochrome. Our results indicate that the final opening of the dimeric structure is preceded by an important internal reorganization of the phytochrome specific (PHY) domain involving a bend of the helical spine connecting the PHY domain with the chromophore-binding domain, opening the way to a new understanding of the activation pathway.

Modeling photophysical properties of the bacteriophytochrome-based fluorescent protein IFP1.4

An enhanced interest in the phytochrome-based fluorescent proteins is explained by their ability to absorb and emit light in the far-red and infra-red regions particularly suitable for bioimaging. The fluorescent protein IFP1.4 was engineered from the chromophore-binding domain of a bacteriophytochrome in attempts to increase the fluorescence quantum yield. We report the results of simulations of structures in the ground S₀ and excited S₁ electronic states of IFP1.4 using the methods of quantum chemistry and quantum mechanics/molecular mechanics. We construct different protonation states of the biliverdin (BV) chromophore in the red-absorbing form of the protein by moving protons from the BV pyrrole rings to a suitable acceptor within the system and show that these structures are close in energy but differ by absorption bands. For the first time, we report structures of the minimum energy conical intersection points S₁/S₀ on the energy surfaces of BV in the protein environment and describe their connection to the local minima in the excited S₁ state. These simulations allow us to characterize the deactivation routes in IFP1.4.

Transient IR spectroscopy identifies key interactions and unravels new intermediates in the photocycle of a bacterial phytochrome

Infra-red spectroscopy advances our understanding of how photosensory proteins carry their function.