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

FURNA: A database for functional annotations of RNA structures

Despite the increasing number of 3D RNA structures in the Protein Data Bank, the majority of experimental RNA structures lack thorough functional annotations. As the significance of the functional roles played by noncoding RNAs becomes increasingly apparent, comprehensive annotation of RNA function is becoming a pressing concern. In response to this need, we have developed FURNA (Functions of RNAs), the first database for experimental RNA structures that aims to provide a comprehensive repository of high-quality functional annotations. These include Gene Ontology terms, Enzyme Commission numbers, ligand-binding sites, RNA families, protein-binding motifs, and cross-references to related databases. FURNA is available at https://seq2fun.dcmb.med.umich.edu/furna/ to enable quick discovery of RNA functions from their structures and sequences.

Structural characterization of two prototypical repressors of SorC family reveals tetrameric assemblies on DNA and mechanism of function

The SorC family of transcriptional regulators plays a crucial role in controlling the carbohydrate metabolism and quorum sensing. We employed an integrative approach combining X-ray crystallography and cryo-electron microscopy to investigate architecture and functional mechanism of two prototypical representatives of two sub-classes of the SorC family: DeoR and CggR from Bacillus subtilis. Despite possessing distinct DNA-binding domains, both proteins form similar tetrameric assemblies when bound to their respective DNA operators. Structural analysis elucidates the process by which the CggR-regulated gapA operon is derepressed through the action of two effectors: fructose-1,6-bisphosphate and newly confirmed dihydroxyacetone phosphate. Our findings provide the first comprehensive understanding of the DNA binding mechanism of the SorC-family proteins, shedding new light on their functional characteristics.

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.

Are kuravirus capsid diameters quantized? The first all-atom genome tracing method for double-stranded DNA viruses

The revolution in cryo-electron microscopy has resulted in unprecedented power to resolve large macromolecular complexes including viruses. Many methods exist to explain density corresponding to proteins and thus entire protein capsids have been solved at the all-atom level. However methods for nucleic acids lag behind, and no all-atom viral double-stranded DNA genomes have been published at all. We here present a method which exploits the spiral winding patterns of DNA in icosahedral capsids. The method quickly generates shells of DNA wound in user-specified, idealized spherical or cylindrical spirals. For transition regions, the method allows guided semiflexible fitting. For the kuravirus SU10, our method explains most of the density in a semiautomated fashion. The results suggest rules for DNA turns in the end caps under which two discrete parameters determine the capsid inner diameter. We suggest that other kuraviruses viruses may follow the same winding scheme, producing a discrete rather than continuous spectrum of capsid inner diameters. Our software may be used to explain the published density maps of other double-stranded DNA viruses and uncover their genome packaging principles.

Isotopic Depletion Increases the Spatial Resolution of FPOP Top-Down Mass Spectrometry Analysis

Protein radical labeling, like fast photochemical oxidation of proteins (FPOP), coupled to a top-down mass spectrometry (MS) analysis offers an alternative analytical method for probing protein structure or protein interaction with other biomolecules, for instance, proteins and DNA. However, with the increasing mass of studied analytes, the MS/MS spectra become complex and exhibit a low signal-to-noise ratio. Nevertheless, these difficulties may be overcome by protein isotope depletion. Thus, we aimed to use protein isotope depletion to analyze FPOP-oxidized samples by top-down MS analysis. For this purpose, we prepared isotopically natural (IN) and depleted (ID) forms of the FOXO4 DNA binding domain (FOXO4-DBD) and studied the protein-DNA interaction interface with double-stranded DNA, the insulin response element (IRE), after exposing the complex to hydroxyl radicals. As shown by comparing tandem mass spectra of natural and depleted proteins, the ID form increased the signal-to-noise ratio of useful fragment ions, thereby enhancing the sequence coverage by more than 19%. This improvement in the detection of fragment ions enabled us to detect 22 more oxidized residues in the ID samples than in the IN sample. Moreover, less common modifications were detected in the ID sample, including the formation of ketones and lysine carbonylation. Given the higher quality of ID top-down MSMS data set, these results provide more detailed information on the complex formation between transcription factors and DNA-response elements. Therefore, our study highlights the benefits of isotopic depletion for quantitative top-down proteomics. Data are available via ProteomeXchange with the identifier PXD044447.

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.

The Nucleic Acid Knowledgebase: a new portal for 3D structural information about nucleic acids

The Nucleic Acid Knowledgebase (nakb.org) is a new data resource, updated weekly, for experimentally determined 3D structures containing DNA and/or RNA nucleic acid polymers and their biological assemblies. NAKB indexes nucleic acid-containing structures derived from all major structure determination methods (X-ray, NMR and EM), including all held by the Protein Data Bank (PDB). As the planned successor to the Nucleic Acid Database (NDB), NAKB's design preserves all functionality of the NDB and provides novel nucleic acid-centric content, including structural and functional annotations, as well as annotations from and links to external resources. A variety of custom interactive tools have been developed to enable rapid exploration and drill-down of NAKB's content.

FURNA: a database for function annotations of RNA structures

Despite the increasing number of 3D RNA structures in the Protein Data Bank, the majority of experimental RNA structures lack thorough functional annotations. As the significance of the functional roles played by non-coding RNAs becomes increasingly apparent, comprehensive annotation of RNA function is becoming a pressing concern. In response to this need, we have developed FURNA (Functions of RNAs), the first database for experimental RNA structures that aims to provide a comprehensive repository of high-quality functional annotations. These include Gene Ontology terms, Enzyme Commission numbers, ligand binding sites, RNA families, protein binding motifs, and cross-references to related databases. FURNA is available at https://seq2fun.dcmb.med.umich.edu/furna/ to enable quick discovery of RNA functions from their structures and sequences.

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

Database-derived water probability densities around structurally and sequentially distinct DNA dinucleotide fragments reproduce the known hydration motifs, which thus can be used as building blocks to predict DNA hydration.

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/.

Structural insight into DNA recognition by bacterial transcriptional regulators of the SorC/DeoR family

The SorC/DeoR family is a large family of bacterial transcription regulators that are involved in the control of carbohydrate metabolism and quorum sensing. To understand the structural basis of DNA recognition, structural studies of two functionally characterized SorC/DeoR family members from Bacillus subtilis were performed: the deoxyribonucleoside regulator bsDeoR and the central glycolytic genes regulator bsCggR. Each selected protein represents one of the subgroups that are recognized within the family. Crystal structures were determined of the N-terminal DNA-binding domains of bsDeoR and bsCggR in complex with DNA duplexes representing the minimal operator sequence at resolutions of 2.3 and 2.1 Å, respectively. While bsDeoR_DBD contains a homeodomain-like HTH-type domain, bsCggR_DBD contains a winged helix-turn-helix-type motif. Both proteins form C2-symmetric dimers that recognize two consecutive major grooves, and the protein-DNA interactions have been analyzed in detail. The crystal structures were used to model the interactions of the proteins with the full DNA operators, and a common mode of DNA recognition is proposed that is most likely to be shared by other members of the SorC/DeoR family.

Understanding the Origin of Structural Diversity of DNA Double Helix

Deciphering the contribution of DNA subunits to the variability of its 3D structure represents an important step toward the elucidation of DNA functions at the atomic level. In the pursuit of that goal, our previous studies revealed that the essential conformational characteristics of the most populated “canonic” BI and AI conformational families of Watson–Crick duplexes, including the sequence dependence of their 3D structure, preexist in the local energy minima of the elemental single-chain fragments, deoxydinucleoside monophosphates (dDMPs). Those computations have uncovered important sequence-dependent regularity in the superposition of neighbor bases. The present work expands our studies to new minimal fragments of DNA with Watson–Crick nucleoside pairs that differ from canonic families in the torsion angles of the sugar-phosphate backbone (SPB). To address this objective, computations have been performed on dDMPs, cdDMPs (complementary dDMPs), and minimal fragments of SPBs of respective systems by using methods of molecular and quantum mechanics. These computations reveal that the conformations of dDMPs and cdDMPs having torsion angles of SPB corresponding to the local energy minima of separate minimal units of SPB exhibit sequence-dependent characteristics representative of canonic families. In contrast, conformations of dDMP and cdDMP with SPB torsions being far from the local minima of separate SPB units exhibit more complex sequence dependence.

New restraints and validation approaches for nucleic acid structures in PDB-REDO

The quality of macromolecular structure models crucially depends on refinement and validation targets, which optimally describe the expected chemistry. Commonly used software for these two procedures has been designed and developed in a protein-centric manner, resulting in relatively few established features for the refinement and validation of nucleic acid-containing structure models. Here, new nucleic acid-specific approaches implemented in PDB-REDO are described, including a new restraint model using noncovalent geometries (base-pair hydrogen bonding and base-pair stacking) as refinement targets. New validation routines are also presented, including a metric for Watson-Crick base-pair geometry normality (ZbpG). Applying the PDB-REDO pipeline with the new restraint model to the whole Protein Data Bank (PDB) demonstrates an overall positive effect on the quality of nucleic acid-containing structure models. Finally, we discuss examples of improvements in the geometry of specific nucleic acid structures in the PDB. The new PDB-REDO models and pipeline are available at https://pdb-redo.eu/.

Beyond the double helix: DNA structural diversity and the PDB

The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.

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.