Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome

Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Importance of Polarizable Embedding for Absorption Spectrum Calculations of Arabidopsis thaliana Cryptochrome 1

Cryptochromes are essential flavoproteins for circadian rhythms and avian magnetoreception. Flavin adenine dinucleotide (FAD), a chromophore within cryptochromes, absorbs blue light, initiating electron transfer processes that lead to a biological signaling cascade. A key step in this cascade is the formation of the FAD semiquinone radical (FADH•), characterized through a specific red-light absorption. The absorption spectra of FADH• in cryptochromes are, however, significantly different from those recorded for the cofactor in solution, primarily due to protein-induced shifts in the absorption peaks. This study employs a multiscale approach, combining molecular dynamics (MD) simulations with quantum mechanical/molecular mechanical (QM/MM) methodologies, to investigate the influence of protein dynamics on embedded FADH• absorption. We emphasize the role of the protein's polarizable environment in the shaping of the absorption spectrum, crucial for accurate spectral predictions in cryptochromes. Our findings provide valuable insights into the absorption process, advancing our understanding of cryptochrome functioning.

Activation of Cryptochrome 4 from Atlantic Herring

Marine fish migrate long distances up to hundreds or even thousands of kilometers for various reasons that include seasonal dependencies, feeding, or reproduction. The ability to perceive the geomagnetic field, called magnetoreception, is one of the many mechanisms allowing some fish to navigate reliably in the aquatic realm. While it is believed that the photoreceptor protein cryptochrome 4 (Cry4) is the key component for the radical pair-based magnetoreception mechanism in night migratory songbirds, the Cry4 mechanism in fish is still largely unexplored. The present study aims to investigate properties of the fish Cry4 protein in order to understand the potential involvement in a radical pair-based magnetoreception. Specifically, a computationally reconstructed atomistic model of Cry4 from the Atlantic herring (Clupea harengus) was studied employing classical molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) methods to investigate internal electron transfers and the radical pair formation. The QM/MM simulations reveal that electron transfers occur similarly to those found experimentally and computationally in Cry4 from European robin (Erithacus rubecula). It is therefore plausible that the investigated Atlantic herring Cry4 has the physical and chemical properties to form radical pairs that in turn could provide fish with a radical pair-based magnetic field compass sensor.

Long-Time Oxygen and Superoxide Localization in Arabidopsis thaliana Cryptochrome

Cryptochromes are proteins that are highly conserved across species and in many instances bind the flavin adenine dinucleotide (FAD) cofactor within their photolyase-homology region (PHR) domain. The FAD cofactor has multiple redox states that help catalyze reactions, and absorbs photons at about 450 nm, a feature linked to the light-related functions of cryptochrome proteins. Reactive oxygen species (ROS) are produced from redox reactions involving molecular oxygen and are involved in a myriad of biological processes. Superoxide O2•- is an exemplary ROS that may be formed through electron transfer from FAD to O2, generating an electron radical pair. Although the formation of a superoxide-FAD radical pair has been speculated, it is still unclear if the required process steps could be realized in cryptochrome. Here, we present results from molecular dynamics (MD) simulations of oxygen interacting with the PHR domain of Arabidopsis thaliana cryptochrome 1 (AtCRY1). Using MD simulation trajectories, oxygen binding locations are characterized through both the O2-FAD intermolecular distance and the local protein environment. Oxygen unbinding times are characterized through replica simulations of the bound oxygen. Simulations reveal that oxygen molecules can localize at certain sites within the cryptochrome protein for tens of nanoseconds, and superoxide molecules can localize for significantly longer. This relatively long-duration molecule binding suggests the possibility of an electron-transfer reaction leading to superoxide formation. Estimates of electron-transfer rates using the Marcus theory are performed for the identified potential binding sites. Molecular oxygen binding results are compared with recent results demonstrating long-time oxygen binding within the electron-transfer flavoprotein (ETF), another FAD binding protein.

Mutational Study of the Tryptophan Tetrad Important for Electron Transfer in European Robin Cryptochrome 4a

The ability of migratory birds to sense magnetic fields has been known for decades, although the understanding of the underlying mechanism is still elusive. Currently, the strongest magnetoreceptor candidate in birds is a protein called cryptochrome 4a. The cryptochrome 4a protein has changed through evolution, apparently endowing some birds with a more pronounced magnetic sensitivity than others. Using phylogenetic tools, we show that a specific tryptophan tetrad and a tyrosine residue predicted to be essential for cryptochrome activation are highly conserved in the avian clade. Through state-of-the-art molecular dynamics simulations and associated analyses, we also studied the role of these specific residues and the associated mutants on the overall dynamics of the protein. The analyses of the single residue mutations were used to judge how far a local change in the protein structure can impact specific dynamics of European robin cryptochrome 4a. We conclude that the replacements of each of the tryptophans one by one with a phenylalanine do not compromise the overall stability of the protein.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Exploring Post-activation Conformational Changes in Pigeon Cryptochrome 4

A widespread hypothesis ascribes the ability of migratory birds to navigate over large distances to an inclination compass realized by the protein cryptochrome in the birds' retinae. Cryptochromes are activated by blue light, which induces a radical pair state, the spin dynamics of which may become sensitive to earth's weak magnetic fields. The magnetic information is encoded and passed on to downstream processes by structural rearrangements of the protein, the details of which remain vague. We utilize extensive all-atom molecular dynamics simulations to probe the conformational changes of pigeon cryptochrome 4 upon light activation. The structural dynamics are analyzed based on principal component analysis and with the help of distance matrices, which reveal significant changes in selected inter-residue distances. The results are evaluated and discussed with reference to the protein structure and its putative function as a magnetoreceptor. It is suggested that the phosphate-binding loop could act as a gate controlling the access to the flavin adenine dinucleotide cofactor depending on the redox state of the protein.

Modeling the effect of ion-induced shock waves and DNA breakage with the reactive CHARMM force field

Ion-induced DNA damage is an important effect underlying ion beam cancer therapy. This article introduces the methodology of modeling DNA damage induced by a shock wave caused by a projectile ion. Specifically it is demonstrated how single- and double strand breaks in a DNA molecule could be described by the reactive CHARMM (rCHARMM) force field implemented in the program MBN Explorer. The entire workflow of performing the shock wave simulations, including obtaining the crucial simulation parameters, is described in seven steps. Two exemplary analyses are provided for a case study simulation serving to: (a) quantify the shock wave propagation and (b) describe the dynamics of formation of DNA breaks. The article concludes by discussing the computational cost of the simulations and revealing the possible maximal computational time for different simulation set-ups.

Structural Explanations of Flavin Adenine Dinucleotide Binding in Drosophila melanogaster Cryptochrome

Cryptochrome proteins are thought to be involved in light-sensitive magnetoreception in migratory birds triggered by flavin adenine dinucleotide (FAD) light absorption. A recent study, however, calls into question the ability of vertebrate cryptochrome proteins to bind FAD, rendering them unlikely to function as magnetoreceptive proteins. In this Letter, we investigate the structural changes occurring in Drosophila melanogaster cryptochrome, upon key amino acid mutations, which reduce FAD binding. Through computational analysis we have now suggested why some mutations do not preclude FAD binding in all vertebrate cryptochrome proteins.

Biological relevance of charge transfer branching pathways in photolyases

The repair of sun-induced DNA lesions by photolyases is driven by a photoinduced electron transfer from a fully reduced FAD to the damaged DNA. A chain of several aromatic residues connecting FAD to solvent ensures the prior photoreduction of the FAD cofactor. In PhrA, a class III CPD photolyase, two branching tryptophan charge transfer pathways have been characterized. According to previous experiments, both pathways play a role in the FAD photoreduction. To provide a molecular insight to the charge transfer abilities of both pathways, we perform multiscales simulations where the protein motion and the positive charge are simultaneously propagated. Our computational approach reveals that one pathway drives a very fast charge transfer whereas the other pathway provides a very good thermodynamic stabilization of the positive charge. During the simulations, the positive charge firstly moves on the fast triad, while a reorganization of the close FAD˙^- environment occurs. Then, backward transfers can lead to the propagation of the positive charge on the second pathway. After one nanosecond, we observe a nearly equal probability to find the charge at ending tryptophan of either pathway; eventually the charge distribution will likely evolve towards a charge stabilization on the last tryptophan of the slowest pathway. Our results highlight the role the protein environment, which manages the association of a kinetic and a thermodynamic pathways to trigger a fast and efficient FAD photoreduction.

Modeling MesoBioNano systems with MBN Studio made easy

This paper introduces MesoBioNano (MBN) Studio - a graphical user interface for a popular multiscale simulation package MBN Explorer. MBN Studio has been developed to facilitate setting up and starting MBN Explorer calculations, monitoring their progress and examining the calculation results. It is tailored for any calculations that are supported by MBN Explorer, such as for example the single-point energy calculations, structure optimization, molecular dynamics, and kinetic Monte Carlo simulations. Apart from that MBN Studio has built-in tools allowing the calculation and analysis of specific characteristics that are determined by the output of the simulations, such as the diffusion coefficients of molecular species, melting temperatures and associated heat capacities, radial distribution function; a dedicated modeling plug-in allows constructing molecular systems in a quick and efficient manner. Employing this plug-in, one can easily construct molecular systems of different geometries (e.g., spherical or ellipsoidal nanoparticles, cubic crystalline samples) with various atomic composition. The paper presents the first public release of MBN Studio and provides an overview of its significant capabilities, as well as the reference point for further extensions.

Absorption Spectra of FAD Embedded in Cryptochromes

The magnetic compass sense utilized by migratory birds for long-distance navigation functions only once light of a certain wavelength is present. This piece of evidence fits partially with the popular hypothesis of chemical magnetoreception in cryptochrome proteins, located in the bird retina. According to this hypothesis a magnetosensitive radical pair is produced after photoexcitation of an FAD cofactor inside cryptochrome, and as such the absorption properties of FAD are of crucial importance for cryptochrome activation. However, we reveal that absorption spectra of FAD show very little variation between six different cryptochromes, suggesting that the electronic transitions are barely affected by the chemical differences in the proteins. This conclusion hints on the presence of a secondary photoreceptor or cofactor that could be necessary to explain green-light-activated magnetoreception in birds.