Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

When will RNA get its AlphaFold moment?

The protein structure prediction problem has been solved for many types of proteins by AlphaFold. Recently, there has been considerable excitement to build off the success of AlphaFold and predict the 3D structures of RNAs. RNA prediction methods use a variety of techniques, from physics-based to machine learning approaches. We believe that there are challenges preventing the successful development of deep learning-based methods like AlphaFold for RNA in the short term. Broadly speaking, the challenges are the limited number of structures and alignments making data-hungry deep learning methods unlikely to succeed. Additionally, there are several issues with the existing structure and sequence data, as they are often of insufficient quality, highly biased and missing key information. Here, we discuss these challenges in detail and suggest some steps to remedy the situation. We believe that it is possible to create an accurate RNA structure prediction method, but it will require solving several data quality and volume issues, usage of data beyond simple sequence alignments, or the development of new less data-hungry machine learning methods.

Diffraction anisotropy and paired refinement: crystal structure of H33, a protein binder to interleukin 10.

Binder H33 is a small protein binder engineered by ribosome display to bind human interleukin 10. Crystals of binder H33 display severe diffraction anisotropy. A set of data files with correction for diffraction anisotropy based on different local signal-to-noise ratios was prepared. Paired refinement was used to find the optimal anisotropic high-resolution diffraction limit of the data: 3.13-2.47 Å. The structure of binder H33 belongs to the 2% of crystal structures with the highest solvent content in the Protein Data Bank.

Regulation of IL-24/IL-20R2 complex formation using photocaged tyrosines and UV light.

Human interleukin 24 (IL-24) is a multifunctional cytokine that represents an important target for autoimmune diseases and cancer. Since the biological functions of IL-24 depend on interactions with membrane receptors, on-demand regulation of the affinity between IL-24 and its cognate partners offers exciting possibilities in basic research and may have applications in therapy. As a proof-of-concept, we developed a strategy based on recombinant soluble protein variants and genetic code expansion technology to photocontrol the binding between IL-24 and one of its receptors, IL-20R2. Screening of non-canonical ortho-nitrobenzyl-tyrosine (NBY) residues introduced at several positions in both partners was done by a combination of biophysical and cell signaling assays. We identified one position for installing NBY, tyrosine70 of IL-20R2, which results in clear impairment of heterocomplex assembly in the dark. Irradiation with 365-nm light leads to decaging and reconstitutes the native tyrosine of the receptor that can then associate with IL-24. Photocaged IL-20R2 may be useful for the spatiotemporal control of the JAK/STAT phosphorylation cascade.

Conformation-based refinement of 18-mer DNA structures.

Nine new crystal structures of CG-rich DNA 18-mers with the sequence 5'-GGTGGGGGC-XZ-GCCCCACC-3', which are related to the bacterial repetitive extragenic palindromes, are reported. 18-mer oligonucleotides with the central XZ dinucleotide systematically mutated to all 16 sequences show complex behavior in solution, but all ten so far successfully crystallized 18-mers crystallized as A-form duplexes. The refinement protocol benefited from the recurrent use of geometries of the dinucleotide conformer (NtC) classes as refinement restraints in regions of poor electron density. The restraints are automatically generated at the dnatco.datmos.org web service and are available for download. This NtC-driven protocol significantly helped to stabilize the structure refinement. The NtC-driven refinement protocol can be adapted to other low-resolution data such as cryo-EM maps. To test the quality of the final structural models, a novel validation method based on comparison of the electron density and conformational similarity to the NtC classes was employed.

Genetically encoded non-canonical amino acids reveal asynchronous dark reversion of chromophore, backbone, and side-chains in EL222.

Photoreceptors containing the light-oxygen-voltage (LOV) domain elicit biological responses upon excitation of their flavin mononucleotide (FMN) chromophore by blue light. The mechanism and kinetics of dark-state recovery are not well understood. Here we incorporated the non-canonical amino acid p-cyanophenylalanine (CNF) by genetic code expansion technology at 45 positions of the bacterial transcription factor EL222. Screening of light-induced changes in infrared (IR) absorption frequency, electric field and hydration of the nitrile groups identified residues CNF31 and CNF35 as reporters of monomer/oligomer and caged/decaged equilibria, respectively. Time-resolved multi-probe UV/visible and IR spectroscopy experiments of the lit-to-dark transition revealed four dynamical events. Predominantly, rearrangements around the A'α helix interface (CNF31 and CNF35) precede FMN-cysteinyl adduct scission, folding of α-helices (amide bands), and relaxation of residue CNF151. This study illustrates the importance of characterizing all parts of a protein and suggests a key role for the N-terminal A'α extension of the LOV domain in controlling EL222 photocycle length.

Combined in vitro and cell-based selection display method producing specific binders against IL-9 receptor in high yields.

We combined cell-free ribosome display and cell-based yeast display selection to build specific protein binders to the extracellular domain of the human interleukin 9 receptor alpha (IL-9Rα). The target, IL-9Rα, is the receptor involved in the signalling pathway of IL-9, a pro-inflammatory cytokine medically important for its involvement in respiratory diseases. The successive use of modified protocols of ribosome and yeast displays allowed us to combine their strengths-the virtually infinite selection power of ribosome display and the production of (mostly) properly folded and soluble proteins in yeast display. The described experimental protocol is optimized to produce binders highly specific to the target, including selectivity to common proteins such as BSA, and proteins potentially competing for the binder such as receptors of other cytokines. The binders were trained from DNA libraries of two protein scaffolds called 57aBi and 57bBi developed in our laboratory. We show that the described unconventional combination of ribosome and yeast displays is effective in developing selective small protein binders to the medically relevant molecular target.

Sub-Millisecond Photoinduced Dynamics of Free and EL222-Bound FMN by Stimulated Raman and Visible Absorption Spectroscopies.

Time-resolved femtosecond-stimulated Raman spectroscopy (FSRS) provides valuable information on the structural dynamics of biomolecules. However, FSRS has been applied mainly up to the nanoseconds regime and above 700 cm-1, which covers only part of the spectrum of biologically relevant time scales and Raman shifts. Here we report on a broadband (~200-2200 cm-1) dual transient visible absorption (visTA)/FSRS set-up that can accommodate time delays from a few femtoseconds to several hundreds of microseconds after illumination with an actinic pump. The extended time scale and wavenumber range allowed us to monitor the complete excited-state dynamics of the biological chromophore flavin mononucleotide (FMN), both free in solution and embedded in two variants of the bacterial light-oxygen-voltage (LOV) photoreceptor EL222. The observed lifetimes and intermediate states (singlet, triplet, and adduct) are in agreement with previous time-resolved infrared spectroscopy experiments. Importantly, we found evidence for additional dynamical events, particularly upon analysis of the low-frequency Raman region below 1000 cm-1. We show that fs-to-sub-ms visTA/FSRS with a broad wavenumber range is a useful tool to characterize short-lived conformationally excited states in flavoproteins and potentially other light-responsive proteins.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks.

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the `spine of hydration' in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Developing Community Resources for Nucleic Acid Structures.

In this review, we describe the creation of the Nucleic Acid Database (NDB) at Rutgers University and how it became a testbed for the current infrastructure of the RCSB Protein Data Bank. We describe some of the special features of the NDB and how it has been used to enable research. Plans for the next phase as the Nucleic Acid Knowledgebase (NAKB) are summarized.

De novo developed protein binders mimicking Interferon lambda signaling.

We hereby describe the process of design and selection of nonantibody protein binders mimicking cytokine signaling. We chose to mimic signaling of IFN-λ1, type 3 interferon (also known as IL-29) for its novelty and the importance of its biological functions. All four known interferons λ signal through binding to the extracellular domains of IL-28 receptor 1 (IL-28R1) and IL-10 receptor 2 (IL-10R2). Our binders were therefore trained to bind both receptors simultaneously. The bifunctional binder molecules were developed by yeast display, a method of directed evolution. The signaling capacity of the bivalent binders was tested by measuring phosphorylation of the JAK/STAT signaling pathway and production of mRNA of six selected genes naturally induced by IFN- λ1 in human cell lines. The newly developed bivalent binders offer opportunities to study cytokine-related biological functions and modulation of the cell behavior by receptor activation on the cell surfaces alternative to the use of natural IFN-λ.

QM calculations predict the energetics and infrared spectra of transient glutamine isomers in LOV photoreceptors.

Photosensory receptors containing the flavin-binding light-oxygen-voltage (LOV) domain are modular proteins that fulfil a variety of biological functions ranging from gene expression to phototropism. The LOV photocycle is initiated by blue-light and involves a cascade of intermediate species, including an electronically excited triplet state, that leads to covalent bond formation between the flavin mononucleotide (FMN) chromophore and a nearby cysteine residue. Subsequent conformational changes in the polypeptide chain arise due to the remodelling of the hydrogen bond network in the cofactor binding pocket, whereby a conserved glutamine residue plays a key role in coupling FMN photochemistry with LOV photobiology. Although the dark-to-light transition of LOV photosensors has been previously addressed by spectroscopy and computational approaches, the mechanistic basis of the underlying reactions is still not well understood. Here we present a detailed computational study of three distinct LOV domains: EL222 from Erythrobacter litoralis, AsLOV2 from the second LOV domain of Avena sativa phototropin 1, and RsLOV from Rhodobacter sphaeroides LOV protein. Extended protein-chromophore models containing all known crucial residues involved in the initial steps (femtosecond-to-microsecond) of the photocycle were employed. Energies and rotational barriers were calculated for possible rotamers and tautomers of the critical glutamine side chain, which allowed us to postulate the most energetically favoured glutamine orientation for each LOV domain along the assumed reaction path. In turn, for each evolving species, infrared difference spectra were constructed and compared to experimental EL222 and AsLOV2 transient infrared spectra, the former from original work presented here and the latter from the literature. The good agreement between theory and experiment permitted the assignment of the majority of observed bands, notably the ∼1635 cm-1 transient of the adduct state to the carbonyl of the glutamine side chain after rotation. Moreover, both the energetic and spectroscopic approaches converge in suggesting a facile glutamine flip at the adduct intermediate for EL222 and more so for AsLOV2, while for RsLOV the glutamine keeps its initial configuration. Additionally, the computed infrared shifts of the glutamine and interacting residues could guide experimental research addressing early events of signal transduction in LOV proteins.

Protein Binder (ProBi) as a New Class of Structurally Robust Non-Antibody Protein Scaffold for Directed Evolution.

Engineered small non-antibody protein scaffolds are a promising alternative to antibodies and are especially attractive for use in protein therapeutics and diagnostics. The advantages include smaller size and a more robust, single-domain structural framework with a defined binding surface amenable to mutation. This calls for a more systematic approach in designing new scaffolds suitable for use in one or more methods of directed evolution. We hereby describe a process based on an analysis of protein structures from the Protein Data Bank and their experimental examination. The candidate protein scaffolds were subjected to a thorough screening including computational evaluation of the mutability, and experimental determination of their expression yield in E. coli, solubility, and thermostability. In the next step, we examined several variants of the candidate scaffolds including their wild types and alanine mutants. We proved the applicability of this systematic procedure by selecting a monomeric single-domain human protein with a fold different from previously known scaffolds. The newly developed scaffold, called ProBi (Protein Binder), contains two independently mutable surface patches. We demonstrated its functionality by training it as a binder against human interleukin-10, a medically important cytokine. The procedure yielded scaffold-related variants with nanomolar affinity.

Demonstration of electron diffraction from membrane protein crystals grown in a lipidic mesophase after lamella preparation by focused ion beam milling at cryogenic temperatures.

Electron crystallography of sub-micrometre-sized 3D protein crystals has emerged recently as a valuable field of structural biology. In meso crystallization methods, utilizing lipidic mesophases, particularly lipidic cubic phases (LCPs), can produce high-quality 3D crystals of membrane proteins (MPs). A major step towards realizing 3D electron crystallography of MP crystals, grown in meso, is to demonstrate electron diffraction from such crystals. The first task is to remove the viscous and sticky lipidic matrix that surrounds the crystals without damaging the crystals. Additionally, the crystals have to be thin enough to let electrons traverse them without significant multiple scattering. In the present work, the concept that focused ion beam milling at cryogenic temperatures (cryo-FIB milling) can be used to remove excess host lipidic mesophase matrix is experimentally verified, and then the crystals are thinned to a thickness suitable for electron diffraction. In this study, bacteriorhodopsin (BR) crystals grown in a lipidic cubic mesophase of monoolein were used as a model system. LCP from a part of a hexagon-shaped plate-like BR crystal (∼10 µm in thickness and ∼70 µm in the longest dimension), which was flash-frozen in liquid nitro-gen, was milled away with a gallium FIB under cryogenic conditions, and a part of the crystal itself was thinned into a ∼210 nm-thick lamella with the ion beam. The frozen sample was then transferred into an electron cryo-microscope, and a nanovolume of ∼1400 × 1400 × 210 nm of the BR lamella was exposed to 200 kV electrons at a fluence of ∼0.06 e Å-2. The resulting electron diffraction peaks were detected beyond 2.7 Šresolution (with an average peak height to background ratio of >2) by a CMOS-based Ceta 16M camera. The results demonstrate that cryo-FIB milling produces high-quality lamellae from crystals grown in lipidic mesophases and pave the way for 3D electron crystallography on crystals grown or embedded in highly viscous media.

Structural variability of CG-rich DNA 18-mers accommodating double T-T mismatches.

Solution and crystal data are reported for DNA 18-mers with sequences related to those of bacterial noncoding single-stranded DNA segments called repetitive extragenic palindromes (REPs). Solution CD and melting data showed that the CG-rich, near-palindromic REPs from various bacterial species exhibit dynamic temperature-dependent and concentration-dependent equilibria, including architectures compatible with not only hairpins, which are expected to be biologically relevant, but also antiparallel duplexes and bimolecular tetraplexes. Three 18-mer oligonucleotides named Hpar-18 (PDB entry 6rou), Chom-18 (PDB entry 6ros) and its brominated variant Chom-18Br (PDB entry 6ror) crystallized as isomorphic right-handed A-like duplexes. The low-resolution crystal structures were solved with the help of experimental phases for Chom-18Br. The center of the duplexes is formed by two successive T-T noncanonical base pairs (mismatches). They do not deform the double-helical geometry. The presence of T-T mismatches prompted an analysis of the geometries of these and other noncanonical pairs in other DNA crystals in terms of their fit to the experimental electron densities (RSCC) and their geometric fit to the NtC (dinucleotide conformational) classes (https://dnatco.datmos.org/). Throughout this work, knowledge of the NtC classes was used to refine and validate the crystal structures, and to analyze the mismatches.

Structural alphabets for conformational analysis of nucleic acids available at dnatco.datmos.org.

A detailed description of the dnatco.datmos.org web server implementing the universal structural alphabet of nucleic acids is presented. It is capable of processing any mmCIF- or PDB-formatted files containing DNA or RNA molecules; these can either be uploaded by the user or supplied as the wwPDB or PDB-REDO structural database access code. The web server performs an assignment of the nucleic acid conformations and presents the results for the intuitive annotation, validation, modeling and refinement of nucleic acids.

A community proposal to integrate structural bioinformatics activities in ELIXIR (3D-Bioinfo Community).

Structural bioinformatics provides the scientific methods and tools to analyse, archive, validate, and present the biomolecular structure data generated by the structural biology community. It also provides an important link with the genomics community, as structural bioinformaticians also use the extensive sequence data to predict protein structures and their functional sites. A very broad and active community of structural bioinformaticians exists across Europe, and 3D-Bioinfo will establish formal platforms to address their needs and better integrate their activities and initiatives. Our mission will be to strengthen the ties with the structural biology research communities in Europe covering life sciences, as well as chemistry and physics and to bridge the gap between these researchers in order to fully realize the potential of structural bioinformatics. Our Community will also undertake dedicated educational, training and outreach efforts to facilitate this, bringing new insights and thus facilitating the development of much needed innovative applications e.g. for human health, drug and protein design. Our combined efforts will be of critical importance to keep the European research efforts competitive in this respect. Here we highlight the major European contributions to the field of structural bioinformatics, the most pressing challenges remaining and how Europe-wide interactions, enabled by ELIXIR and its platforms, will help in addressing these challenges and in coordinating structural bioinformatics resources across Europe. In particular, we present recent activities and future plans to consolidate an ELIXIR 3D-Bioinfo Community in structural bioinformatics and propose means to develop better links across the community. These include building new consortia, organising workshops to establish data standards and seeking community agreement on benchmark data sets and strategies. We also highlight existing and planned collaborations with other ELIXIR Communities and other European infrastructures, such as the structural biology community supported by Instruct-ERIC, with whom we have synergies and overlapping common interests.

A unified dinucleotide alphabet describing both RNA and DNA structures.

By analyzing almost 120 000 dinucleotides in over 2000 nonredundant nucleic acid crystal structures, we define 96+1 diNucleotide Conformers, NtCs, which describe the geometry of RNA and DNA dinucleotides. NtC classes are grouped into 15 codes of the structural alphabet CANA (Conformational Alphabet of Nucleic Acids) to simplify symbolic annotation of the prominent structural features of NAs and their intuitive graphical display. The search for nontrivial patterns of NtCs resulted in the identification of several types of RNA loops, some of them observed for the first time. Over 30% of the nearly six million dinucleotides in the PDB cannot be assigned to any NtC class but we demonstrate that up to a half of them can be re-refined with the help of proper refinement targets. A statistical analysis of the preferences of NtCs and CANA codes for the 16 dinucleotide sequences showed that neither the NtC class AA00, which forms the scaffold of RNA structures, nor BB00, the DNA most populated class, are sequence neutral but their distributions are significantly biased. The reported automated assignment of the NtC classes and CANA codes available at dnatco.org provides a powerful tool for unbiased analysis of nucleic acid structures by structural and molecular biologists.

Femtosecond-to-nanosecond dynamics of flavin mononucleotide monitored by stimulated Raman spectroscopy and simulations.

Flavin mononucleotide (FMN) belongs to the large family of flavins, ubiquitous yellow-coloured biological chromophores that contain an isoalloxazine ring system. As a cofactor in flavoproteins, it is found in various enzymes and photosensory receptors, like those featuring the light-oxygen-voltage (LOV) domain. The photocycle of FMN is triggered by blue light and proceeds via a cascade of intermediate states. In this work, we have studied isolated FMN in an aqueous solution in order to elucidate the intrinsic electronic and vibrational changes of the chromophore upon excitation. The ultrafast transitions of excited FMN were monitored through the joint use of femtosecond stimulated Raman spectroscopy (FSRS) and transient absorption spectroscopy encompassing a time window between 0 ps and 6 ns with 50 fs time resolution. Global analysis of the obtained transient visible absorption and transient Raman spectra in combination with extensive quantum chemistry calculations identified unambiguously the singlet and triplet FMN populations and addressed solvent dynamics effects. The good agreement between the experimental and theoretical spectra facilitated the assignment of electronic transitions and vibrations. Our results represent the first steps towards more complex experiments aimed at tracking structural changes of FMN embedded in light-inducible proteins upon photoexcitation.