Magnetic sensitivity of cryptochrome 4 from a migratory songbird.

Night-migratory songbirds are remarkably proficient navigators1. Flying alone and often over great distances, they use various directional cues including, crucially, a light-dependent magnetic compass2,3. The mechanism of this compass has been suggested to rely on the quantum spin dynamics of photoinduced radical pairs in cryptochrome flavoproteins located in the retinas of the birds4-7. Here we show that the photochemistry of cryptochrome 4 (CRY4) from the night-migratory European robin (Erithacus rubecula) is magnetically sensitive in vitro, and more so than CRY4 from two non-migratory bird species, chicken (Gallus gallus) and pigeon (Columba livia). Site-specific mutations of ErCRY4 reveal the roles of four successive flavin-tryptophan radical pairs in generating magnetic field effects and in stabilizing potential signalling states in a way that could enable sensing and signalling functions to be independently optimized in night-migratory birds.

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.

Upper bound for broadband radiofrequency field disruption of magnetic compass orientation in night-migratory songbirds.

Night-migratory songbirds have a light-dependent magnetic compass sense, the mechanism of which is thought to depend on the photochemical formation of radical pairs in cryptochrome (Cry) proteins located in the retina. The finding that weak radiofrequency (RF) electromagnetic fields can prevent birds from orienting in the Earth's magnetic field has been regarded as a diagnostic test for this mechanism and as a potential source of information on the identities of the radicals. The maximum frequency that could cause such disorientation has been predicted to lie between 120 and 220 MHz for a flavin-tryptophan radical pair in Cry. Here we show that the magnetic orientation capabilities of Eurasian blackcaps (Sylvia atricapilla) are not affected by RF noise in the frequency bands 140 to 150 MHz and 235 to 245 MHz. From a consideration of its internal magnetic interactions, we argue that RF field effects on a flavin-containing radical-pair sensor should be approximately independent of frequency up to 116 MHz and that birds' sensitivity to RF disorientation should fall by about two orders of magnitude when the frequency exceeds 116 MHz. Taken together with our earlier finding that 75 to 85 MHz RF fields disrupt the magnetic orientation of blackcaps, these results provide compelling evidence that the magnetic compass of migratory birds operates by a radical pair mechanism.

Adaptive evolution and loss of a putative magnetoreceptor in passerines.

Migratory birds possess remarkable accuracy in orientation and navigation, which involves various compass systems including the magnetic compass. Identifying the primary magnetosensor remains a fundamental open question. Cryptochromes (Cry) have been shown to be magnetically sensitive, and Cry4a from a migratory songbird seems to show enhanced magnetic sensitivity in vitro compared to Cry4a from resident species. We investigate Cry and their potential involvement in magnetoreception in a phylogenetic framework, integrating molecular evolutionary analyses with protein dynamics modelling. Our analysis is based on 363 bird genomes and identifies different selection regimes in passerines. We show that Cry4a is characterized by strong positive selection and high variability, typical characteristics of sensor proteins. We identify key sites that are likely to have facilitated the evolution of an optimized sensory protein for night-time orientation in songbirds. Additionally, we show that Cry4 was lost in hummingbirds, parrots and Tyranni (Suboscines), and thus identified a gene deletion, which might facilitate testing the function of Cry4a in birds. In contrast, the other avian Cry (Cry1 and Cry2) were highly conserved across all species, indicating basal, non-sensory functions. Our results support a specialization or functional differentiation of Cry4 in songbirds which could be magnetosensation.

Peculiar Differences between Two Copper Complexes Containing Similar Redox-Active Ligands: Density Functional and Multiconfigurational Calculations.

Transition metal complexes featuring redox-active ligands often exhibit multiple redox states, influenced by the interplay between the metal center and the ligand. This study delves into the electronic structures of two mononuclear complexes of copper with two similar redox-active urea azine ligands. The ligands differ by the replacement of an NCH3 moiety by an S atom in the ligand backbone. Experimental analysis yields pronounced electronic structural disparities between these complexes, observable in both the solution and solid phases. Conventional quantum chemical methods, such as density functional theory using different functionals (B3LYP, TPSSh, and CAM-B3LYP), remain inadequate to rationalize the observed spectroscopic anomalies. However, a multiconfigurational approach elucidates the disparate behaviors of these complexes. Multireference perturbation theory, based on complete active space self-consistent field computations, identifies Cu(I) in the case of the complex with the NCH3 containing ligands and a state with substantial Cu(II) contributions in the case of the complex with the S atom containing ligands. In contrast, DFT indicates Cu(I) in both scenarios.

Different receptor models show differences in ligand binding strength and location: a computational drug screening for the tick-borne encephalitis virus.

The tick-borne encephalitis virus (TBE) is a neurotrophic disease that has spread more rapidly throughout Europe and Asia in the past few years. At the same time, no cure or specific therapy is known to battle the illness apart from vaccination. To find a pharmacologically relevant drug, a computer-aided drug screening was initiated. Such a procedure probes a possible binding of a drug to the RNA Polymerase of TBE. The crystal structure of the receptor, however, includes missing and partially modeled regions, which rendered the structure incomplete and of questionable use for a thorough drug screening procedure. The quality of the receptor model was addressed by studying three putative structures created. We show that the choice of receptor models greatly influences the binding affinity of potential drug molecules and that the binding location could also be significantly impacted. We demonstrate that some drug candidates are unsuitable for one model but show decent results for another. Without any prejudice on the three employed receptor models, the study reveals the imperative need to investigate the receptor structure before drug binding is probed whether experimentally or computationally.

Tracking the Electron Transfer Cascade in European Robin Cryptochrome 4 Mutants.

The primary step in the mechanism by which migratory birds sense the Earth's magnetic field is thought to be the light-induced formation of long-lived magnetically sensitive radical pairs within cryptochrome flavoproteins located in the birds' retinas. Blue-light absorption by the non-covalently bound flavin chromophore triggers sequential electron transfers along a chain of four tryptophan residues toward the photoexcited flavin. The recently demonstrated ability to express cryptochrome 4a from the night-migratory European robin (Erithacus rubecula), ErCry4a, and to replace each of the tryptophan residues by a redox-inactive phenylalanine offers the prospect of exploring the roles of the four tryptophans. Here, we use ultrafast transient absorption spectroscopy to compare wild type ErCry4a and four mutants having a phenylalanine at different positions in the chain. We find that each of the three tryptophan residues closest to the flavin adds a distinct relaxation component (time constants: 0.5, 30, and 150 ps) in the transient absorption data. The dynamics of the mutant containing a phenylalanine at the fourth position, furthest from the flavin, are very similar to those of wild type ErCry4a, except for a reduced concentration of long-lived radical pairs. The experimental results are evaluated and discussed in the framework of real-time quantum mechanical/molecular mechanical electron transfer simulations based on the density functional-based tight binding approach. This comparison between simulation results and experimental measurements provides a detailed microscopic insight into the sequential electron transfers along the tryptophan chain. Our results offer a route to the study of spin transport and dynamical spin correlations in flavoprotein radical pairs.

Conformationally Restricted Glycopeptide Backbone Inhibits Gas-Phase H/D Scrambling between Glycan and Peptide Moieties.

Protein glycosylation is a common post-translational modification on extracellular proteins. The conformational dynamics of several glycoproteins have been characterized by hydrogen/deuterium exchange mass spectrometry (HDX-MS). However, it is, in most cases, not possible to extract information about glycan conformation and dynamics due to the general difficulty of separating the deuterium content of the glycan from that of the peptide (in particular, for O-linked glycans). Here, we investigate whether the fragmentation of protonated glycopeptides by collision-induced dissociation (CID) can be used to determine the solution-specific deuterium content of the glycan. Central to this concept is that glycopeptides can undergo a facile loss of glycans upon CID, thereby allowing for the determination of their masses. However, an essential prerequisite is that hydrogen and deuterium (H/D) scrambling can be kept in check. Therefore, we have measured the degree of scrambling upon glycosidic bond cleavage in glycopeptides that differ in the conformational flexibility of their backbone and glycosylation pattern. Our results show that complete scrambling precedes the glycosidic bond cleavage in normal glycopeptides derived from a glycoprotein; i.e., all labile hydrogens have undergone positional randomization prior to loss of the glycan. In contrast, the glycosidic bond cleavage occurs without any scrambling in the glycopeptide antibiotic vancomycin, reflecting that the glycan cannot interact with the peptide moiety due to a conformationally restricted backbone as revealed by molecular dynamics simulations. Scrambling is also inhibited, albeit to a lesser degree, in the conformationally restricted glycopeptides ristocetin and its pseudoaglycone, demonstrating that scrambling depends on an intricate interplay between the flexibility and proximity of the glycan and the peptide backbone.

Modeling spin relaxation in complex radical systems using MolSpin.

Spin relaxation is an important aspect of the spin dynamics of free radicals and can have a significant impact on the outcome of their spin-selective reactions. Examples range from the use of radicals as spin qubits in quantum information processing to the radical pair reactions in proteins that may allow migratory birds to sense the direction of the Earth's magnetic field. Accurate modeling of spin relaxation, however, is non-trivial. Bloch-Redfield-Wangsness theory derives a quantum mechanical master equation from system-bath interactions in the Markovian limit that provides a comprehensive framework for describing spin relaxation. Unfortunately, the construction of the master equation is system-specific and often resource-heavy. To address this challenge, we introduce a generalized and efficient implementation of BRW theory as a new feature of the spin dynamics toolkit MolSpin which offers an easy-to-use approach for studying systems of reacting radicals of varying complexity.

Introducing the Automated Ligand Searcher.

The Automated Ligand Searcher (ALISE) is designed as an automated computational drug discovery tool. To approximate the binding free energy of ligands to a receptor, ALISE includes a three-stage workflow, with each stage involving an increasingly sophisticated computational method: molecular docking, molecular dynamics, and free energy perturbation, respectively. To narrow the number of potential ligands, poorly performing ligands are gradually segregated out. The performance and usability of ALISE are benchmarked for a case study containing known active ligands and decoys for the HIV protease. The example illustrates that ALISE filters the decoys successfully and demonstrates that the automation, comprehensiveness, and user-friendliness of the software make it a valuable tool for improved and faster drug development workflows.

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.

Spin Dynamics of Flavoproteins.

This short review reports the surprising phenomenon of nuclear hyperpolarization occurring in chemical reactions, which is called CIDNP (chemically induced dynamic nuclear polarization) or photo-CIDNP if the chemical reaction is light-driven. The phenomenon occurs in both liquid and solid-state, and electron transfer systems, often carrying flavins as electron acceptors, are involved. Here, we explain the physical and chemical properties of flavins, their occurrence in spin-correlated radical pairs (SCRP) and the possible involvement of flavin-carrying SCRPs in animal magneto-reception at earth's magnetic field.

Effects of Dynamical Degrees of Freedom on Magnetic Compass Sensitivity: A Comparison of Plant and Avian Cryptochromes.

The magnetic compass of migratory birds is thought to rely on a radical pair reaction inside the blue-light photoreceptor protein cryptochrome. The sensitivity of such a sensor to weak external magnetic fields is determined by a variety of magnetic interactions, including electron-nuclear hyperfine interactions. Here, we investigate the implications of thermal motion, focusing on fluctuations in the dihedral and librational angles of flavin adenine dinucleotide (FAD) and tryptophan (Trp) radicals in cryptochrome 4a from European robin (Erithacus rubecula, ErCry4a) and pigeon (Columba livia, ClCry4a) and cryptochrome 1 from the plant Arabidopsis thaliana (AtCry1). Molecular dynamics simulations and density functional theory-derived hyperfine interactions are used to calculate the quantum yield of radical pair recombination dependent on the direction of the geomagnetic field. This quantity and various dynamical parameters are compared for [FAD•- Trp•+] in ErCry4a, ClCry4a, and AtCry1, with TrpC or TrpD being the third and fourth components of the tryptophan triad/tetrad in the respective proteins. We find that (i) differences in the average dihedral angles in the radical pairs are small, (ii) the librational motions of TrpC•+ in the avian cryptochromes are appreciably smaller than in AtCry1, (iii) the rapid vibrational motions of the radicals leading to strong fluctuations in the hyperfine couplings affect the spin dynamics depending on the usage of instantaneous or time-averaged interactions. Future investigations of radical pair compass sensitivity should therefore not be based on single snapshots of the protein structure but should include the ensemble properties of the hyperfine interactions.

Multiscale modeling of stochastic dynamics processes with MBN Explorer.

Stochastic dynamics describes processes in complex systems having the probabilistic nature. They can involve very different dynamical systems and occur on very different temporal and spatial scale. This paper discusses the concept of stochastic dynamics and its implementation in the popular program MBN Explorer. Stochastic dynamics in MBN Explorer relies on the Monte Carlo approach and permits simulations of physical, chemical, and biological processes. The paper presents the basic theoretical concepts underlying stochastic dynamics implementation and provides several examples highlighting its applicability to different systems, such as diffusing proteins seeking an anchor point on a cell membrane, deposition of nanoparticles on a surface leading to structures with fractal morphologies, and oscillations of compounds in an autocatalytic reaction. The chosen examples illustrate the diversity of applications that can be modeled by means of stochastic dynamics with MBN Explorer.

Long-Time Oxygen Localization in Electron Transfer Flavoprotein.

Reactive oxygen species (ROS) exert a wide range of biological effects from beneficial regulatory function to deleterious oxidative stress. The electron transfer flavoprotein (ETF) is ubiquitous to life and is associated with aerobic metabolism and ROS production due to its location in the mitochondria. Quantifying oxygen localization within the ETF complex is critical for understanding the potential for electron transfer and radical pair formation between flavin adenine dinucleotide (FAD) cofactor and superoxide during ROS formation. Our study employed all-atom molecular dynamics simulations and identified several novel, long-lived oxygen binding sites within the ETF complex that appear near the FAD cofactor. Site locations, the local electrostatic environment, and characteristic oxygen binding times for each site were evaluated to establish factors that may lead to possible charge transfer reactions and superoxide formation within the ETF complex. The study revealed that some oxygen binding sites are naturally linked to protein domain features, suggesting opportunities to engineer and control ROS production and subsequent dynamics.

The Transition of Photoreceptor Guanylate Cyclase Type 1 to the Active State.

Membrane-bound guanylate cyclases (GCs), which synthesize the second messenger guanosine-3', 5'-cyclic monophosphate, differ in their activation modes to reach the active state. Hormone peptides bind to the extracellular domain in hormone-receptor-type GCs and trigger a conformational change in the intracellular, cytoplasmic part of the enzyme. Sensory GCs that are present in rod and cone photoreceptor cells have intracellular binding sites for regulatory Ca2+-sensor proteins, named guanylate-cyclase-activating proteins. A rotation model of activation involving an α-helix rotation was described as a common activation motif among hormone-receptor GCs. We tested whether the photoreceptor GC-E underwent an α-helix rotation when reaching the active state. We experimentally simulated such a transitory switch by integrating alanine residues close to the transmembrane region, and compared the effects of alanine integration with the point mutation V902L in GC-E. The V902L mutation is found in patients suffering from retinal cone-rod dystrophies, and leads to a constitutively active state of GC-E. We analyzed the enzymatic catalytic parameters of wild-type and mutant GC-E. Our data showed no involvement of an α-helix rotation when reaching the active state, indicating a difference in hormone receptor GCs. To characterize the protein conformations that represent the transition to the active state, we investigated the protein dynamics by using a computational approach based on all-atom molecular dynamics simulations. We detected a swinging movement of the dimerization domain in the V902L mutant as the critical conformational switch in the cyclase going from the low to high activity state.

Design of SARS-CoV-2 Main Protease Inhibitors Using Artificial Intelligence and Molecular Dynamic Simulations.

Drug design is a time-consuming and cumbersome process due to the vast search space of drug-like molecules and the difficulty of investigating atomic and electronic interactions. The present paper proposes a computational drug design workflow that combines artificial intelligence (AI) methods, i.e., an evolutionary algorithm and artificial neural network model, and molecular dynamics (MD) simulations to design and evaluate potential drug candidates. For the purpose of illustration, the proposed workflow was applied to design drug candidates against the main protease of severe acute respiratory syndrome coronavirus 2. From the ∼140,000 molecules designed using AI methods, MD analysis identified two molecules as potential drug candidates.

Theoretical Description of the Primary Proton-Coupled Electron Transfer Reaction in the Cytochrome bc1 Complex.

The cytochrome bc1 complex is a transmembrane enzymatic protein complex that plays a central role in cellular energy production and is present in both photosynthetic and respiratory chain organelles. Its reaction mechanism is initiated by the binding of a quinol molecule to an active site, followed by a series of charge transfer reactions between the quinol and protein subunits. Previous work hypothesized that the primary reaction was a concerted proton-coupled electron transfer (PCET) reaction because of the apparent absence of intermediate states associated with single proton or electron transfer reactions. In the present study, the kinetics of the primary bc1 complex PCET reaction is investigated with a vibronically nonadiabatic PCET theory in conjunction with all-atom molecular dynamics simulations and electronic structure calculations. The computed rate constants and relatively high kinetic isotope effects are consistent with experimental measurements on related biomimetic systems. The analysis implicates a concerted PCET mechanism with significant hydrogen tunneling and nonadiabatic effects in the bc1 complex. Moreover, the employed theoretical framework is shown to serve as a general strategy for describing PCET reactions in bioenergetic systems.

Introducing Pep McConst-A user-friendly peptide modeler for biophysical applications.

We are introducing Pep McConst-a software that employs a Monte-Carlo algorithm to construct 3D structures of polypeptide chains which could subsequently be studied as stand-alone macromolecules or complement the structure of known proteins. Using an approach to avoid steric clashes, Pep McConst allows to create multiple structures for a predefined primary sequence of amino acids. These structures could then effectively be used for further structural analysis and investigations. The article introduces the algorithm and describes its user-friendly approach that was made possible through the VIKING online platform. Finally, the manuscript provides several highlight examples where Pep McConst was used to predict the structure of the C-terminal of a known protein, generate a missing bit of already crystallized protein structures and simply generate short polypeptide chains.

Inhibition Mechanism of Antimalarial Drugs Targeting the Cytochrome bc1 Complex.

Plasmodium falciparum (P. falciparum) is the main parasite known to cause malaria in humans. The antimalarial drug atovaquone is known to inhibit the Qo-site of the cytochrome bc1 complex of P. falciparum, which ultimately blocks ATP synthesis, leading to cell death. Through the years, mutations of the P. falciparum cytochrome bc1 complex, causing resistance to atovaquone, have emerged. The present investigation applies molecular dynamics (MD) simulations to study how the specific mutations Y279S and L282V, known to cause atovaquone resistance in malarial parasites, affect the inhibition mechanism of two known inhibitors. Binding free energy estimates were obtained through free energy perturbation calculations but were unable to confidently resolve the effects of mutations due to the great complexity of the binding environment. Meanwhile, basic mechanistic considerations from the MD simulations provide a detailed characterization of inhibitor binding modes and indicate that the Y279S mutation weakens the natural binding of the inhibitors, while no conclusive effect of the L282V mutation could be observed.