Introducing VIKING: A Novel Online Platform for Multiscale Modeling

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.

Magnetic orientation in juvenile Atlantic herring (Clupea harengus) could involve cryptochrome 4 as a potential magnetoreceptor

The Earth's magnetic field can provide reliable directional information, allowing migrating animals to orient themselves using a magnetic compass or estimate their position relative to a target using map-based orientation. Here we show for the first time that young, inexperienced herring (Clupea harengus, Ch) have a magnetic compass when they migrate hundreds of kilometres to their feeding grounds. In birds, such as the European robin (Erithacus rubecula), radical pair-based magnetoreception involving cryptochrome 4 (ErCRY4) was demonstrated; the molecular basis of magnetoreception in fish is still elusive. We show that cry4 expression in the eye of herring is upregulated during the migratory season, but not before, indicating a possible use for migration. The amino acid structure of herring ChCRY4 shows four tryptophans and a flavin adenine dinucleotide-binding site, a prerequisite for a magnetic receptor. Using homology modelling, we successfully reconstructed ChCRY4 of herring, DrCRY4 of zebrafish (Danio rerio) and StCRY4 of brown trout (Salmo trutta) and showed that ChCRY4, DrCRY4 and ErCRY4a, but not StCRY4, exhibit very comparable dynamic behaviour. The electron transfer could take place in ChCRY4 in a similar way to ErCRY4a. The combined behavioural, transcriptomic and simulation experiments provide evidence that CRY4 could act as a magnetoreceptor in Atlantic herring.

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.

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.

How do Antimicrobial Peptides Interact with the Outer Membrane of Gram-Negative Bacteria? Role of Lipopolysaccharides in Peptide Binding, Anchoring, and Penetration

Gram-negative bacteria possess a complex structural cell envelope that constitutes a barrier for antimicrobial peptides that neutralize the microbes by disrupting their cell membranes. Computational and experimental approaches were used to study a model outer membrane interaction with an antimicrobial peptide, melittin. The investigated membrane included di[3-deoxy-d-manno-octulosonyl]-lipid A (KLA) in the outer leaflet and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) in the inner leaflet. Molecular dynamics simulations revealed that the positively charged helical C-terminus of melittin anchors rapidly into the hydrophilic headgroup region of KLA, while the flexible N-terminus makes contacts with the phosphate groups of KLA, supporting melittin penetration into the boundary between the hydrophilic and hydrophobic regions of the lipids. Electrochemical techniques confirmed the binding of melittin to the model membrane. To probe the peptide conformation and orientation during interaction with the membrane, polarization modulation infrared reflection absorption spectroscopy was used. The measurements revealed conformational changes in the peptide, accompanied by reorientation and translocation of the peptide at the membrane surface. The study suggests that melittin insertion into the outer membrane affects its permeability and capacitance but does not disturb the membrane's bilayer structure, indicating a distinct mechanism of the peptide action on the outer membrane of Gram-negative bacteria.

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.

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.

A marine cryptochrome with an inverse photo-oligomerization mechanism

Cryptochromes (CRYs) are a structurally conserved but functionally diverse family of proteins that can confer unique sensory properties to organisms. In the marine bristle worm Platynereis dumerilii, its light receptive cryptochrome L-CRY (PdLCry) allows the animal to discriminate between sunlight and moonlight, an important requirement for synchronizing its lunar cycle-dependent mass spawning. Using cryo-electron microscopy, we show that in the dark, PdLCry adopts a dimer arrangement observed neither in plant nor insect CRYs. Intense illumination disassembles the dimer into monomers. Structural and functional data suggest a mechanistic coupling between the light-sensing flavin adenine dinucleotide chromophore, the dimer interface, and the C-terminal tail helix, with a likely involvement of the phosphate binding loop. Taken together, our work establishes PdLCry as a CRY protein with inverse photo-oligomerization with respect to plant CRYs, and provides molecular insights into how this protein might help discriminating the different light intensities associated with sunlight and moonlight.

DataPype: A Fully Automated Unified Software Platform for Computer-Aided Drug Design

With the advent of computer-aided drug design (CADD), traditional physical testing of thousands of molecules has now been replaced by target-focused drug discovery, where potentially bioactive molecules are predicted by computer software before their physical synthesis. However, despite being a significant breakthrough, CADD still faces various limitations and challenges. The increasing availability of data on small molecules has created a need to streamline the sourcing of data from different databases and automate the processing and cleaning of data into a form that can be used by multiple CADD software applications. Several standalone software packages are available to aid the drug designer, each with its own specific application, requiring specialized knowledge and expertise for optimal use. These applications require their own input and output files, making it a challenge for nonexpert users or multidisciplinary discovery teams. Here, we have developed a new software platform called DataPype, which wraps around these different software packages. It provides a unified automated workflow to search for hit compounds using specialist software. Additionally, multiple virtual screening packages can be used in the one workflow, and if different ways of looking at potential hit compounds all predict the same set of molecules, we have higher confidence that we should make or purchase and test the molecules. Importantly, DataPype can run on computer servers, speeding up the virtual screening for new compounds. Combining access to multiple CADD tools within one interface will enhance the early stage of drug discovery, increase usability, and enable the use of parallel computing.

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.

Dimerization of European Robin Cryptochrome 4a

Homo-dimer formation is important for the function of many proteins. Although dimeric forms of cryptochromes (Cry) have been found by crystallography and were recently observed in vitro for European robin Cry4a, little is known about the dimerization of avian Crys and the role it could play in the mechanism of magnetic sensing in migratory birds. Here, we present a combined experimental and computational investigation of the dimerization of robin Cry4a resulting from covalent and non-covalent interactions. Experimental studies using native mass spectrometry, mass spectrometric analysis of disulfide bonds, chemical cross-linking, and photometric measurements show that disulfide-linked dimers are routinely formed, that their formation is promoted by exposure to blue light, and that the most likely cysteines are C317 and C412. Computational modeling and molecular dynamics simulations were used to generate and assess a number of possible dimer structures. The relevance of these findings to the proposed role of Cry4a in avian magnetoreception is discussed.

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.

Computational Reconstruction and Analysis of Structural Models of Avian Cryptochrome 4

A recent study by Xu et al. (Nature, 2021, 594, 535-540) provided strong evidence that cryptochrome 4 (Cry4) is a key protein to endow migratory birds with the magnetic compass sense. The investigation compared the magnetic field response of Cry4 from migratory and nonmigratory bird species and suggested that a difference in magnetic sensitivity could exist. This finding prompted an in-depth investigation into Cry4 protein differences on the structural and dynamic levels. In the present study, the pigeon Cry4 (ClCry4) crystal structure was used to reconstruct the missing avian Cry4 protein structures via homology modeling for carefully selected bird species. The reconstructed Cry4 structure from European robin, Eurasian blackcap, zebra finch, chicken, and pigeon were subsequently simulated dynamically and analyzed. The studied avian Cry4 structures show flexibility in analogous regions pointing to similar activation mechanisms and/or signaling interaction partners. It can be concluded that the experimentally recorded difference in the magnetic field sensitivity of Cry4 from different birds is unlikely to be due to solely intrinsic dynamics of the proteins but requires additional factors that have not yet been identified.

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 Ca²⁺-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.

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 Energy Landscape of Side Reactions in the Cytochrome bc1 Complex

Much of the metabolic molecular machinery responsible for energy transduction processes in living organisms revolves around a series of electron and proton transfer processes. The highly redox active enzymes can, however, also pose a risk of unwanted side reactions leading to reactive oxygen species, which are harmful to cells and are a factor in aging and age-related diseases. Using extensive quantum and classical computational modeling, we here show evidence of a particular superoxide production mechanism through stray reactions between molecular oxygen and a semiquinone reaction intermediate bound in the mitochondrial complex III of the electron transport chain, also known as the cytochrome b c 1 complex. Free energy calculations indicate a favorable electron transfer from semiquinone occurring at low rates under normal circumstances. Furthermore, simulations of the product state reveal that superoxide formed at the Q _o -site exclusively leaves the b c 1 complex at the positive side of the membrane and escapes into the intermembrane space of mitochondria, providing a critical clue in further studies of the harmful effects of mitochondrial superoxide production.

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.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

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.