Analyzing and building nucleic acid structures with 3DNA

Computational Identification and Illustrative Standard for Representation of Unimolecular G-Quadruplex Secondary Structures (CIIS-GQ)

G-quadruplexes refer to a large group of nucleic acid-based structures. In recent years, they have been attracting attention due to their biological roles in the telomeres and promoter regions. These structures show wide diversity in topology, however, development of methods for structural classification of G-quadruplexes has been evaded for a long time. There has been a limited number of studies aiming to bring forth a secondary structure classification method. The situation was even more complex than imagined, since the discovery of bulged and mismatched G-quadruplexes while most of the available tools fail to distinguish these non-canonical G-quadruplex motifs. Moreover, the interpretation of their analysis output still requires expert knowledge. In this study, we propose a new method for identification of unimolecular G-Quadruplexes and classification by secondary structures based on three-dimensional structural data. Briefly, coordinates of guanines are processed to identify tetrads, loops and bulges. Then, we present the secondary structure in the form of a depiction which shows the loop types, bulges, and guanines that participate in each tetrad. Moreover, CIIS-GQ identifies non-guanine nucleotides that joins the G-tetrads and forms multiplets. Finally, the results of our study are compared with DSSR and ElTetrado classification methods, and the advantages of the proposed depiction method for representing secondary structures were discussed. The source code of the method can be accessed via https://github.com/TugayDirek/CIIS-GQ .

Investigating the Effect of Chemical Modifications on the Ribose Sugar Conformation, Watson-Crick Base Pairing, and Intrastrand Stacking Interactions: A Theoretical Approach

Over the last few decades, chemically modified sugars have been incorporated into nucleic acid-based therapeutics to improve their pharmacological potential. Chemical modification can influence the sugar conformation, Watson-Crick hydrogen (W-C) bonding, and nucleobase stacking interactions, which play major roles in the structural integrity and dynamic properties of nucleic acid duplexes. In this study, we categorized 33 uridine (U*) and cytidine (C*) sugar modifications and calculated their sugar conformational parameters. We also calculated the Watson-Crick hydrogen bond energies of the modified RNA-type base pairs (U*:A and C*:G) using DFT and sSAPT0 methods. The W-C base pairing energy calculations suggested that the South-type modified sugar strengthens the C*:G base pair and weakens the U*:A base pair compared to the unmodified one. In contrast, the North-type sugar modifications form weaker C*:G base pair and marginally stronger U*:A base pair compared to the South-type modified sugars. Moreover, intrastrand base stacking energies were calculated for 15 modifications incorporated at the fourth position in 7-mer non-self-complementary RNA duplexes [(GCAU*GAC)2 and (GCAC*GAC)2], utilizing molecular dynamics simulation and quantum mechanical (DFT and sSAPT0) methods. The sugar modifications were found to have minimal effect on the intrastrand base-stacking interactions. However, the glycol nucleic acid modification disturbs the intrastrand base-stacking significantly, corroborating the experimental data.

3dDNA: A Computational Method of Building DNA 3D Structures

Considerable progress has been made in the prediction methods of 3D structures of RNAs. In contrast, no such methods are available for DNAs. The determination of 3D structures of the latter is also increasingly needed for understanding their functions and designing new DNA molecules. Since the number of experimental structures of DNA is limited at present, here, we propose a computational and template-based method, 3dDNA, which combines DNA and RNA template libraries to predict DNA 3D structures. It was benchmarked on three test sets with different numbers of chains, and the results show that 3dDNA can predict DNA 3D structures with a mean RMSD of about 2.36 Å for those with one or two chains and fewer than 4 Å with three or more chains.

Computational normal mode analysis accurately replicates the activity and specificity profiles of CRISPR-Cas9 and high-fidelity variants

The CRISPR-Cas system has transformed the field of gene-editing and created opportunities for novel genome engineering therapeutics. The field has significantly progressed, and recently, CRISPR-Cas9 was utilized in clinical trials to target disease-causing mutations. Existing tools aim to predict the on-target efficacy and potential genome-wide off-targets by scoring a particular gRNA according to an array of gRNA design principles or machine learning algorithms based on empirical results of large numbers of gRNAs. However, such tools are unable to predict the editing outcome by variant Cas enzymes and can only assess potential off-targets related to reference genomes. Here, we employ normal mode analysis (NMA) to investigate the structure of the Cas9 protein complexed with its gRNA and target DNA and explore the function of the protein. Our results demonstrate the feasibility and validity of NMA to predict the activity and specificity of SpyCas9 in the presence of mismatches by comparison to empirical data. Furthermore, despite the absence of their exact structures, this method accurately predicts the enzymatic activity of known high-fidelity engineered Cas9 variants.

Deciphering the mechanism of p73 recognition of p53 response elements using the crystal structure of p73-DNA complexes and computational studies

P73 belongs to p53 family transcription factor activating more than 50% of cell fate p53 target genes involved in cell cycle, apoptosis, DNA damage response alongside neuronal system development and differentiation by binding to 20-bp response elements (REs) having sequence motif (PPPCA/T-T/AGYYY) where P-purines and Y-pyrimidines with each 10-bp separated by minimum 0 to 13-bp spacer. The promiscuous nature of recognizing both cell fate and development genes and the underlying RE selectivity mechanism by p73 is not well understood. Here, we report the molecular details of p73 recognizing the REs using the crystal structure of p73 DNA binding domain (DBD) in complex with 12 base pair DNA sequence 5'-cAGGCATGCCTg-3' and molecular dynamics simulations with six different p53 natural promoter sequences. Each 20-base pair natural promoter forms a different major/minor groove due to the presence of nucleotides A/T, A/C, G/G, T/T and G/T at positions 3, 8, 13, 18 uniquely recognized by p73 key residues Lys138 and Arg268. The loops L1 and L3 bearing these residues influence inter-and intra-dimer interfaces interactions and hence p73 forms a unique tetramer with each natural promoter sequence. Structural features of the DNA and the spacing between half-sites influence p73 tetramerization and its transactivation function.

A reproducibility analysis-based statistical framework for residue-residue evolutionary coupling detection

Direct coupling analysis (DCA) has been widely used to infer evolutionary coupled residue pairs from the multiple sequence alignment (MSA) of homologous sequences. However, effectively selecting residue pairs with significant evolutionary couplings according to the result of DCA is a non-trivial task. In this study, we developed a general statistical framework for significant evolutionary coupling detection, referred to as irreproducible discovery rate (IDR)-DCA, which is based on reproducibility analysis of the coupling scores obtained from DCA on manually created MSA replicates. IDR-DCA was applied to select residue pairs for contact prediction for monomeric proteins, protein-protein interactions and monomeric RNAs, in which three different versions of DCA were applied. We demonstrated that with the application of IDR-DCA, the residue pairs selected using a universal threshold always yielded stable performance for contact prediction. Comparing with the application of carefully tuned coupling score cutoffs, IDR-DCA always showed better performance. The robustness of IDR-DCA was also supported through the MSA downsampling analysis. We further demonstrated the effectiveness of applying constraints obtained from residue pairs selected by IDR-DCA to assist RNA secondary structure prediction.

Coarse-grained molecular dynamics simulations of base-pair mismatch recognition protein MutS sliding along DNA

DNA mismatches are frequently generated by various intrinsic and extrinsic factors including DNA replication errors, oxygen species, ultraviolet, and ionizing radiation. These mismatches should be corrected by the mismatches repair (MMR) pathway to maintain genome integrity. In the Escherichia coli (E. coli) MMR pathway, MutS searches and recognizes a base-pair mismatch from millions of base-pairs. Once recognized, ADP bound to MutS is exchanged with ATP, which induces a conformational change in MutS. Previous single-molecule fluorescence microscopy studies have suggested that ADP-bound MutS temporarily slides along double-stranded DNA in a rotation-coupled manner to search a base-pair mismatch and so does ATP-bound MutS in a rotation-uncoupled manner. However, the detailed structural dynamics of the sliding remains unclear. In this study, we performed coarse-grained molecular dynamics simulations of the E. coli MutS bound on DNA in three different conformations: ADP-bound (MutS^ADP), ATP-bound open clamp (MutSOpenATP), and ATP-bound closed clamp (MutSClosedATP) conformations. In the simulations, we observed conformation-dependent diffusion of MutS along DNA. MutS^ADP and MutSClosedATP diffused along DNA in a rotation-coupled manner with rare and frequent groove-crossing events, respectively. In the groove-crossing events, MutS overcame an edge of a groove and temporarily diffused in a rotation-uncoupled manner. It was also indicated that mismatch searches by MutSOpenATP is inefficient in terms of mismatch checking even though it diffuses along DNA and reaches unchecked regions more rapidly than MutS^ADP.

Vfold2D-MC: A Physics-Based Hybrid Model for Predicting RNA Secondary Structure Folding

Accurate prediction of RNA structure and folding stability has a far-reaching impact on our understanding of RNA functions. Here we develop Vfold2D-MC, a new physics-based model, to predict RNA structure and folding thermodynamics from the sequence. The model employs virtual bond-based coarse-graining of RNA backbone conformation and generates RNA conformations through Monte Carlo sampling of the bond angles and torsional angles of the virtual bonds. Using a coarse-grained statistical potential derived from the known structures, we assign each conformation with a statistical weight. The weighted average over the conformational ensemble gives the entropy and free energy parameters for the hairpin, bulge, and internal loops, and multiway junctions. From the thermodynamic parameters, we predict RNA structures, melting curves, and structural changes from the sequence. Theory-experiment comparisons indicate that Vfold2D-MC not only gives improved structure predictions but also enables the interpretation of thermodynamic results for different RNA structures, including multibranched junctions. This new model sets a promising framework to treat more complicated RNA structures, such as pseudoknotted and intramolecular kissing loops, for which experimental thermodynamic parameters are often unavailable.

Effective Drug Concentration and Selectivity Depends on Fraction of Primitive Cells

Poor efficiency of chemotherapeutics in the eradication of Cancer Stem Cells (CSCs) has been driving the search for more active and specific compounds. In this work, we show how cell density-dependent stage culture profiles can be used in drug development workflows to achieve more robust drug activity (IC₅₀ and EC₅₀) results. Using flow cytometry and light microscopy, we characterized the cytological stage profiles of the HL-60-, A-549-, and HEK-293-derived sublines with a focus on their primitive cell content. We then used a range of cytotoxic substances—C-123, bortezomib, idarubicin, C-1305, doxorubicin, DMSO, and ethanol—to highlight typical density-related issues accompanying drug activity determination. We also showed that drug EC₅₀ and selectivity indices normalized to primitive cell content are more accurate activity measurements. We tested our approach by calculating the corrected selectivity index of a novel chemotherapeutic candidate, C-123. Overall, our study highlights the usefulness of accounting for primitive cell fractions in the assessment of drug efficiency.

Distinct mechanisms control genome recognition by p53 at its target genes linked to different cell fates

The tumor suppressor p53 integrates stress response pathways by selectively engaging one of several potential transcriptomes, thereby triggering cell fate decisions (e.g., cell cycle arrest, apoptosis). Foundational to this process is the binding of tetrameric p53 to 20-bp response elements (REs) in the genome (RRRCWWGYYYN_0-13RRRCWWGYYY). In general, REs at cell cycle arrest targets (e.g. p21) are of higher affinity than those at apoptosis targets (e.g., BAX). However, the RE sequence code underlying selectivity remains undeciphered. Here, we identify molecular mechanisms mediating p53 binding to high- and low-affinity REs by showing that key determinants of the code are embedded in the DNA shape. We further demonstrate that differences in minor/major groove widths, encoded by G/C or A/T bp content at positions 3, 8, 13, and 18 in the RE, determine distinct p53 DNA-binding modes by inducing different Arg248 and Lys120 conformations and interactions. The predictive capacity of this code was confirmed in vivo using genome editing at the BAX RE to interconvert the DNA-binding modes, transcription pattern, and cell fate outcome.

Structural Studies of HNA Substrate Specificity in Mutants of an Archaeal DNA Polymerase Obtained by Directed Evolution

Archaeal DNA polymerases from the B-family (polB) have found essential applications in biotechnology. In addition, some of their variants can accept a wide range of modified nucleotides or xenobiotic nucleotides, such as 1,5-anhydrohexitol nucleic acid (HNA), which has the unique ability to selectively cross-pair with DNA and RNA. This capacity is essential to allow the transmission of information between different chemistries of nucleic acid molecules. Variants of the archaeal polymerase from Thermococcus gorgonarius, TgoT, that can either generate HNA from DNA (TgoT_6G12) or DNA from HNA (TgoT_RT521) have been previously identified. To understand how DNA and HNA are recognized and selected by these two laboratory-evolved polymerases, we report six X-ray structures of these variants, as well as an in silico model of a ternary complex with HNA. Structural comparisons of the apo form of TgoT_6G12 together with its binary and ternary complexes with a DNA duplex highlight an ensemble of interactions and conformational changes required to promote DNA or HNA synthesis. MD simulations of the ternary complex suggest that the HNA-DNA hybrid duplex remains stable in the A-DNA helical form and help explain the presence of mutations in regions that would normally not be in contact with the DNA if it were not in the A-helical form. One complex with two incorporated HNA nucleotides is surprisingly found in a one nucleotide-backtracked form, which is new for a DNA polymerase. This information can be used for engineering a new generation of more efficient HNA polymerase variants.

Atypical structures of GAA/TTC trinucleotide repeats underlying Friedreich's ataxia: DNA triplexes and RNA/DNA hybrids

Expansion of the GAA/TTC repeats in the first intron of the FXN gene causes Friedreich's ataxia. Non-canonical structures are linked to this expansion. DNA triplexes and R-loops are believed to arrest transcription, which results in frataxin deficiency and eventual neurodegeneration. We present a systematic in silico characterization of the possible DNA triplexes that could be assembled with GAA and TTC strands; the two hybrid duplexes [r(GAA):d(TTC) and d(GAA):r(UUC)] in an R-loop; and three hybrid triplexes that could form during bidirectional transcription when the non-template DNA strand bonds with the hybrid duplex (collapsed R-loops, where the two DNA strands remain antiparallel). For both Y·R:Y and R·R:Y DNA triplexes, the parallel third strand orientation is more stable; both parallel and antiparallel protonated d(GA+A)·d(GAA):d(TTC) triplexes are stable. Apparent contradictions in the literature about the R·R:Y triplex stability is probably due to lack of molecular resolution, since shifting the third strand by a single nucleotide alters the stability ranking. In the collapsed R-loops, antiparallel d(TTC+)·d(GAA):r(UUC) is unstable, while parallel d(GAA)·r(GAA):d(TTC) and d(GA+A)·r(GAA):d(TTC) are stable. In addition to providing new structural perspectives for specific therapeutic aims, our results contribute to a systematic structural basis for the emerging field of quantitative R-loop biology.

Topology-based classification of tetrads and quadruplex structures

Motivation

Quadruplexes attract the attention of researchers from many fields of bio-science. Due to a specific structure, these tertiary motifs are involved in various biological processes. They are also promising therapeutic targets in many strategies of drug development, including anticancer and neurological disease treatment. The uniqueness and diversity of their forms cause that quadruplexes show great potential in novel biological applications. The existing approaches for quadruplex analysis are based on sequence or 3D structure features and address canonical motifs only.

Results

In our study, we analyzed tetrads and quadruplexes contained in nucleic acid molecules deposited in Protein Data Bank. Focusing on their secondary structure topology, we adjusted its graphical diagram and proposed new dot-bracket and arc representations. We defined the novel classification of these motifs. It can handle both canonical and non-canonical cases. Based on this new taxonomy, we implemented a method that automatically recognizes the types of tetrads and quadruplexes occurring as unimolecular structures. Finally, we conducted a statistical analysis of these motifs found in experimentally determined nucleic acid structures in relation to the new classification.

Availability and implementation

https://github.com/tzok/eltetrado/

Supplementary information

Supplementary data are available at Bioinformatics online.

Different views of the dynamic landscape covered by the 5'-hairpin of the 7SK small nuclear RNA

The 7SK small nuclear RNA (7SKsnRNA) plays a key role in the regulation of RNA polymerase II by sequestrating and inhibiting the positive transcription elongation factor b (P-TEFb) in the 7SK ribonucleoprotein complex (7SKsnRNP), a process mediated by interaction with the protein HEXIM. P-TEFb is also an essential cellular factor recruited by the viral protein Tat to ensure the replication of the viral RNA in the infection cycle of the human immunodeficiency virus (HIV-1). Tat promotes the release of P-TEFb from the 7SKsnRNP and subsequent activation of transcription, by displacing HEXIM from the 5'-hairpin of the 7SKsnRNA. This hairpin (HP1), comprising the signature sequence of the 7SKsnRNA, has been the subject of three independent structural studies aimed at identifying the structural features that could drive the recognition by the two proteins, both depending on arginine-rich motifs (ARM). Interestingly, four distinct structures were determined. In an attempt to provide a comprehensive view of the structure-function relationship of this versatile RNA, we present here a structural analysis of the models, highlighting how HP1 is able to adopt distinct conformations with significant impact on the compactness of the molecule. Since these models are solved under different conditions by nuclear magnetic resonance (NMR) and crystallography, the impact of the buffer composition on the conformational variation was investigated by complementary biophysical approaches. Finally, using isothermal titration calorimetry, we determined the thermodynamic signatures of the Tat-ARM and HEXIM-ARM peptide interactions with the RNA, showing that they are associated with distinct binding mechanisms.

3dRNA: Building RNA 3D structure with improved template library

Most of computational methods of building RNA tertiary structure are template-based. The template-based methods usually can give more accurate 3D structures due to the use of native 3D templates, but they cannot work if the 3D templates are not available. So, a more complete library of the native 3D templates is very important for this type of methods. 3dRNA is a template-based method for building RNA tertiary structure previously proposed by us. In this paper we report improved 3D template libraries of 3dRNA by using two different schemes that give two libraries 3dRNA_Lib1 and 3dRNA_Lib2. These libraries expand the original one by nearly ten times. Benchmark shows that they can significantly increase the accuracy of 3dRNA, especially in building complex and large RNA 3D structures.

Structural mechanism for replication origin binding and remodeling by a metazoan origin recognition complex and its co-loader Cdc6

Eukaryotic DNA replication initiation relies on the origin recognition complex (ORC), a DNA-binding ATPase that loads the Mcm2–7 replicative helicase onto replication origins. Here, we report cryo-electron microscopy (cryo-EM) structures of DNA-bound Drosophila ORC with and without the co-loader Cdc6. These structures reveal that Orc1 and Orc4 constitute the primary DNA binding site in the ORC ring and cooperate with the winged-helix domains to stabilize DNA bending. A loop region near the catalytic Walker B motif of Orc1 directly contacts DNA, allosterically coupling DNA binding to ORC’s ATPase site. Correlating structural and biochemical data show that DNA sequence modulates DNA binding and remodeling by ORC, and that DNA bending promotes Mcm2–7 loading in vitro. Together, these findings explain the distinct DNA sequence-dependencies of metazoan and S. cerevisiae initiators in origin recognition and support a model in which DNA geometry and bendability contribute to Mcm2–7 loading site selection in metazoans.

The origin recognition complex (ORC) is essential for loading the Mcm2–7 replicative helicase onto DNA during DNA replication initiation. Here, the authors describe several cryo-electron microscopy structures of Drosophila ORC bound to DNA and its cofactor Cdc6 and also report an in vitro reconstitution system for Drosophila Mcm2–7 loading, revealing unexpected features of ORC’s DNA binding and remodeling mechanism during Mcm2–7 loading.

Canonical DNA minor groove insertion of bisbenzamidine-Ru(ii) complexes with chiral selectivity

We report the first Ru(ii) coordination compounds that interact with DNA through a canonical minor groove insertion mode and with selectivity for A/T rich sites. This was made possible by integrating a bis-benzamidine minor groove DNA-binding agent with a ruthenium(ii) complex. Importantly, one of the enantiomers (Δ-[Ru(bpy)2 b4bpy]2+, Δ-4Ru) shows a considerably higher DNA affinity than the parent organic ligand and the other enantiomer, particularly for the AATT sequence, while the other enantiomer preferentially targets long AAATTT sites with overall lower affinity. Finally, we demonstrate that the photophysical properties of these new binders can be exploited for DNA cleavage using visible light.

Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures

Web 3DNA (w3DNA) 2.0 is a significantly enhanced version of the widely used w3DNA server for the analysis, visualization, and modeling of 3D nucleic-acid-containing structures. Since its initial release in 2009, the w3DNA server has continuously served the community by making commonly-used features of the 3DNA suite of command-line programs readily accessible. However, due to the lack of updates, w3DNA has clearly shown its age in terms of modern web technologies and it has long lagged behind further developments of 3DNA per se. The w3DNA 2.0 server presented here overcomes all known shortcomings of w3DNA while maintaining its battle-tested characteristics. Technically, w3DNA 2.0 implements a simple and intuitive interface (with sensible defaults) for increased usability, and it complies with HTML5 web standards for broad accessibility. Featurewise, w3DNA 2.0 employs the most recent version of 3DNA, enhanced with many new functionalities, including: the automatic handling of modified nucleotides; a set of 'simple' base-pair and step parameters for qualitative characterization of non-Watson-Crick double-helical structures; new structural parameters that integrate the rigid base plane and the backbone phosphate group, the two nucleic acid components most reliably determined with X-ray crystallography; in silico base mutations that preserve the backbone geometry; and a notably improved module for building models of single-stranded RNA, double-helical DNA, Pauling triplex, G-quadruplex, or DNA structures 'decorated' with proteins. The w3DNA 2.0 server is freely available, without registration, at http://web.x3dna.org.

Physical binding of the tobacco smoke carcinogen NNK diazonium ion to the human tumor suppressor gene TP53 Exon 5

The tobacco smoke N-nitrosamine, NNK, is an important carcinogen. It has been shown to induce lung, liver, and pancreatic cancer in animals. Its metabolites are associated with lung cancer in tobacco smokers. Our work focuses upon the physical interaction of NNK diazonium ion with DNA. This species is implicated in the formation of pyridyloxobutyl adducts, reacting with DNA bases and phosphate groups. Past research has investigated the metabolic activation of NNK by various enzymes, subsequent adduct formation with DNA, and the role of these adducts in mutagenesis. We present the first study of the physical interaction of NNK diazonium ion with TP53 (exon 5), a frequently mutated human tumor suppressor gene. We identify physical binding sites found via free energy minimization in computational docking simulations. These structures represent local potential energy minima in this system and suggest plausible sites for adduct formation.

The yeast C/D box snoRNA U14 adopts a "weak" K-turn like conformation recognized by the Snu13 core protein in solution

Non-coding RNAs associate with proteins to form ribonucleoproteins (RNPs), such as ribosome, box C/D snoRNPs, H/ACA snoRNPs, ribonuclease P, telomerase and spliceosome to ensure cell viability. The assembly of these RNA-protein complexes relies on the ability of the RNA to adopt the correct bound conformation. K-turn motifs represent ubiquitous binding platform for proteins found in several cellular environment. This structural motif has an internal three-nucleotide bulge flanked on its 3' side by a G•A/A•G tandem pairs followed by one or two non-Watson-Crick pairs, and on its 5' side by a classical RNA helix. This peculiar arrangement induces a strong curvature of the phosphodiester backbone, which makes it conducive to multiple tertiary interactions. SNU13/Snu13p (Human/Yeast) binds specifically the U14 C/D box snoRNA K-turn sequence motif. This event is the prerequisite to promote the assembly of the RNP, which contains NOP58/Nop58 and NOP56/Nop56 core proteins and the 2'-O-methyl-transferase, Fibrillarin/Nop1p. The U14 small nucleolar RNA is a conserved non-coding RNA found in yeast and vertebrates required for the pre-rRNA maturation and ribose methylation. Here, we report the solution structure of the free U14 snoRNA K-turn motif (kt-U14) as determined by Nuclear Magnetic Resonance. We demonstrate that a major fraction of free kt-U14 adopts a pre-folded conformation similar to protein bound K-turn, even in the absence of divalent ions. In contrast to the kt-U4 or tyrS RNA, kt-U14 displays a sharp bent in the phosphodiester backbone. The U•U and G•A tandem base pairs are formed through weak hydrogen bonds. Finally, we show that the structure of kt-U14 is stabilized upon Snu13p binding. The structure of the free U14 RNA is the first reference example for the canonical motifs of the C/D box snoRNA family.

NanoModeler: A Webserver for Molecular Simulations and Engineering of Nanoparticles

Functionalized nanoparticles (NPs) are at the frontier of nanoscience. They hold the promise of innovative applications for human health and technology. In this context, molecular dynamics (MD) simulations of NPs are increasingly employed to understand the fundamental structural and dynamical features of NPs. While informative, such simulations demand a laborious two-step process for their setup. In-house scripts are required to (i) construct complex 3D models of the inner metal core and outer layer of organic ligands, and (ii) correctly assign force-field parameters to these composite systems. Here, we present NanoModeler ( www.nanomodeler.it ), the first Webserver designed to automatically generate and parametrize model systems of monolayer-protected gold NPs and gold nanoclusters. The only required input is a structure file of one or two ligand(s) to be grafted onto the gold core, with the option of specifying homogeneous or heterogeneous NP morphologies. NanoModeler then generates 3D models of the nanosystem and the associated topology files. These files are ready for use with the Gromacs MD engine, and they are compatible with the AMBER family of force fields. We illustrate NanoModeler's capabilities with MD simulations of selected representative NP model systems. NanoModeler is the first platform to automate and standardize the construction and parametrization of realistic models for atomistic simulations of gold NPs and gold nanoclusters.

Structure and Biophysics for a Six Letter DNA Alphabet that Includes Imidazo[1,2-a]-1,3,5-triazine-2(8H)-4(3H)-dione (X) and 2,4-Diaminopyrimidine (K)

A goal of synthetic biology is to develop new nucleobases that retain the desirable properties of natural nucleobases at the same time as expanding the genetic alphabet. The nonstandard Watson-Crick pair between imidazo[1,2-a]-1,3,5-triazine-2(8H)-4(3H)-dione (X) and 2,4-diaminopyrimidine (K) does exactly this, pairing via complementary arrangements of hydrogen bonding in these two nucleobases, which do not complement any natural nucleobase. Here, we report the crystal structure of a duplex DNA oligonucleotide in B-form including two consecutive X:K pairs in GATCXK DNA determined as a host-guest complex at 1.75 Å resolution. X:K pairs have significant propeller twist angles, similar to those observed for A:T pairs, and a calculated hydrogen bonding pairing energy that is weaker than that of A:T. Thus, although inclusion of X:K pairs results in a duplex DNA structure that is globally similar to that of an analogous G:C structure, the X:K pairs locally and energetically more closely resemble A:T pairs.

SupeRNAlign: a new tool for flexible superposition of homologous RNA structures and inference of accurate structure-based sequence alignments

RNA has been found to play an ever-increasing role in a variety of biological processes. The function of most non-coding RNA molecules depends on their structure. Comparing and classifying macromolecular 3D structures is of crucial importance for structure-based function inference and it is used in the characterization of functional motifs and in structure prediction by comparative modeling. However, compared to the numerous methods for protein structure superposition, there are few tools dedicated to the superimposing of RNA 3D structures. Here, we present SupeRNAlign (v1.3.1), a new method for flexible superposition of RNA 3D structures, and SupeRNAlign-Coffee—a workflow that combines SupeRNAlign with T-Coffee for inferring structure-based sequence alignments. The methods have been benchmarked with eight other methods for RNA structural superposition and alignment. The benchmark included 151 structures from 32 RNA families (with a total of 1734 pairwise superpositions). The accuracy of superpositions was assessed by comparing structure-based sequence alignments to the reference alignments from the Rfam database. SupeRNAlign and SupeRNAlign-Coffee achieved significantly higher scores than most of the benchmarked methods: SupeRNAlign generated the most accurate sequence alignments among the structure superposition methods, and SupeRNAlign-Coffee performed best among the sequence alignment methods.

How proteins bind to DNA: target discrimination and dynamic sequence search by the telomeric protein TRF1

Target search as performed by DNA-binding proteins is a complex process, in which multiple factors contribute to both thermodynamic discrimination of the target sequence from overwhelmingly abundant off-target sites and kinetic acceleration of dynamic sequence interrogation. TRF1, the protein that binds to telomeric tandem repeats, faces an intriguing variant of the search problem where target sites are clustered within short fragments of chromosomal DNA. In this study, we use extensive (>0.5 ms in total) MD simulations to study the dynamical aspects of sequence-specific binding of TRF1 at both telomeric and non-cognate DNA. For the first time, we describe the spontaneous formation of a sequence-specific native protein-DNA complex in atomistic detail, and study the mechanism by which proteins avoid off-target binding while retaining high affinity for target sites. Our calculated free energy landscapes reproduce the thermodynamics of sequence-specific binding, while statistical approaches allow for a comprehensive description of intermediate stages of complex formation.

Crystal structure of an engineered, HIV-specific recombinase for removal of integrated proviral DNA

As part of the HIV infection cycle, viral DNA inserts into the genome of host cells such that the integrated DNA encoding the viral proteins is flanked by long terminal repeat (LTR) regions from the retrovirus. In an effort to develop novel genome editing techniques that safely excise HIV provirus from cells, Tre, an engineered version of Cre recombinase, was designed to target a 34-bp sequence within the HIV-1 LTR (loxLTR). The sequence targeted by Tre lacks the symmetry present in loxP, the natural DNA substrate for Cre. We report here the crystal structure of a catalytically inactive (Y324F) mutant of this engineered Tre recombinase in complex with the loxLTR DNA substrate. We also report that 17 of the 19 amino acid changes relative to Cre contribute to the altered specificity, even though many of these residues do not contact the DNA directly. We hypothesize that some mutations increase the flexibility of the Cre tetramer and that this, along with flexibility in the DNA, enable the engineered enzyme and DNA substrate to adopt complementary conformations.

Predicting Ion Effects in an RNA Conformational Equilibrium

We develop a partial charge-based tightly bound ion (PCTBI) model for the ion effects in RNA folding. On the basis of the Monte Carlo tightly bound ion (MCTBI) approach, the model can account for ion fluctuation and correlation effects, and can predict the ion distribution around the RNA. Furthermore, unlike the previous coarse-grained RNA charge models, where negative charges are placed on the phosphates only, the current new model considers the detailed all-atom partial charge distribution on the RNA. Thus, the model not only keeps the advantage of the MCTBI model, but also has the potential to provide important detailed information unattainable by the previous MCTBI models. For example, the model predicts the reduction in ion binding upon protein binding and ion-induced conformational switches. For hepatitis C virus genomic RNA, the model predicts a Mg²⁺-induced stabilization of a kissing motif for a cis-acting regulatory element in the genomic RNA. Extensive theory-experiment comparisons support the reliability of the theoretical predictions. Therefore, the model may serve as a robust starting point for further development of an accurate method for ion effects in an RNA conformational equilibrium and RNA-cofactor interactions.

Solution structure of the 5'-terminal hairpin of the 7SK small nuclear RNA

The small nuclear 7SK RNA regulates RNA polymerase II (RNA Pol II) transcription, by sequestering and inhibiting the positive transcription elongation factor b (P-TEFb). P-TEFb is stored in the 7SK ribonucleoprotein (RNP) that contains the three nuclear proteins Hexim1, LaRP7, and MePCE. P-TEFb interacts with the protein Hexim1 and the 7SK RNA. Once P-TEFb is released from the 7SK RNP, it activates transcription by phosphorylating the C-terminal domain of RNA Pol II. P-TEFb also plays a crucial role in the replication of the human immunodeficiency virus HIV-1, through its recruitment by the viral transactivator Tat. Previous work demonstrated that the protein Tat promotes the release of P-TEFb from the 7SK RNP through direct binding to the 7SK RNA. Hexim1 and Tat proteins both comprise conserved and similar arginine-rich motifs that were identified to bind the 7SK RNA at a repeated GAUC site located at the top of the 5'-terminal hairpin (HPI). Here, we report the solution structure of this region as determined by nuclear magnetic resonance, to identify HPI structural features recognized by Hexim1 and Tat. The HPI solution structure displays an elongated shape featuring four helical segments interrupted by one internal loop and three bulges with distinct folds. In particular, the repeated GAUC motif adopts a pre-organized geometry. Our results suggest that the binding of Hexim1 and Tat to the 7SK RNA could originate from a conformational selection of this motif, highlighting how RNA local structure could lead to an adaptive recognition of their partners.

Crystal structure of a poly(rA) staggered zipper at acidic pH: evidence that adenine N1 protonation mediates parallel double helix formation

We have solved at 1.07 Å resolution the X-ray crystal structure of a polyriboadenylic acid (poly(rA)) parallel and continuous double helix. Fifty-nine years ago, double helices of poly(rA) were first proposed to form at acidic pH. Here, we show that 7-mer oligo(rA), i.e. rA₇, hybridizes and overlaps in all registers at pH 3.5 to form stacked double helices that span the crystal. Under these conditions, rA₇ forms well-ordered crystals, whereas rA₆ forms fragile crystalline-like structures, and rA₅, rA₈ and rA₁₁ fail to crystallize. Our findings support studies from ∼50 years ago: one showed using spectroscopic methods that duplex formation at pH 4.5 largely starts with rA₇ and begins to plateau with rA₈; another proposed a so-called ‘staggered zipper’ model in which oligo(rA) strands overlap in multiple registers to extend the helical duplex. While never shown, protonation of adenines at position N1 has been hypothesized to be critical for helix formation. Bond angles in our structure suggest that N1 is protonated on the adenines of every other rAMP−rAMP helix base pair. Our data offer new insights into poly(rA) duplex formation that may be useful in developing a pH sensor.

Numeric promoter description - A comparative view on concepts and general application

Nucleic acid molecules play a key role in a variety of biological processes. Starting from storage and transfer tasks, this also comprises the triggering of biological processes, regulatory effects and the active influence gained by target binding. Based on the experimental output (in this case promoter sequences), further in silico analyses aid in gaining new insights into these processes and interactions. The numerical description of nucleic acids thereby constitutes a bridge between the concrete biological issues and the analytical methods. Hence, this study compares 26 descriptor sets obtained by applying well-known numerical description concepts to an established dataset of 38 DNA promoter sequences. The suitability of the description sets was evaluated by computing partial least squares regression models and assessing the model accuracy. We conclude that the major importance regarding the descriptive power is attached to positional information rather than to explicitly incorporated physico-chemical information, since a sufficient amount of implicit physico-chemical information is already encoded in the nucleobase classification. The regression models especially benefited from employing the information that is encoded in the sequential and structural neighborhood of the nucleobases. Thus, the analyses of n-grams (short fragments of length n) suggested that they are valuable descriptors for DNA target interactions. A mixed n-gram descriptor set thereby yielded the best description of the promoter sequences. The corresponding regression model was checked and found to be plausible as it was able to reproduce the characteristic binding motifs of promoter sequences in a reasonable degree. As most functional nucleic acids are based on the principle of molecular recognition, the findings are not restricted to promoter sequences, but can rather be transferred to other kinds of functional nucleic acids. Thus, the concepts presented in this study could provide advantages for future nucleic acid-based technologies, like biosensoring, therapeutics and molecular imaging.

The systematic approach to describing conformational rearrangements in G-quadruplexes

Conformational changes in DNA G-quadruplex (GQ)-forming regions affect genome function and, thus, compose an interesting research topic. Computer modelling may yield insight into quadruplex folding and rearrangement, particularly molecular dynamics simulations. Here, we show that specific parameters, which are distinct from those commonly used in DNA conformational analyses, must be introduced for adequate interpretation and, most importantly, convenient visual representation of the quadruplex modelling results. We report a set of parameters that comprehensively and systematically describe GQ geometry in dynamics. The parameters include those related to quartet planarity, quadruplex twist, and quartet stacking; they are used to quantitatively characterise various types of quadruplexes and rearrangements, such as quartet distortion/disruption or deviation/bulging of a single nucleotide from the quartet plane. Our approach to describing conformational changes in quadruplexes using the new parameters is exemplified by telomeric quadruplex rearrangement, and the benefits of applying this approach to analyse other structures are discussed.

What controls DNA looping?

The looping of DNA provides a means of communication between sequentially distant genomic sites that operate in tandem to express, copy, and repair the information encoded in the DNA base sequence. The short loops implicated in the expression of bacterial genes suggest that molecular factors other than the naturally stiff double helix are involved in bringing the interacting sites into close spatial proximity. New computational techniques that take direct account of the three-dimensional structures and fluctuations of protein and DNA allow us to examine the likely means of enhancing such communication. Here, we describe the application of these approaches to the looping of a 92 base-pair DNA segment between the headpieces of the tetrameric Escherichia coli Lac repressor protein. The distortions of the double helix induced by a second protein--the nonspecific nucleoid protein HU--increase the computed likelihood of looping by several orders of magnitude over that of DNA alone. Large-scale deformations of the repressor, sequence-dependent features in the DNA loop, and deformability of the DNA operators also enhance looping, although to lesser degrees. The correspondence between the predicted looping propensities and the ease of looping derived from gene-expression and single-molecule measurements lends credence to the derived structural picture.

The structure of an authentic spore photoproduct lesion in DNA suggests a basis for recognition

The spore photoproduct lesion (SP; 5-thymine-5,6-dihydrothymine) is the dominant photoproduct found in UV-irradiated spores of some bacteria such as Bacillus subtilis. Upon spore germination, this lesion is repaired in a light-independent manner by a specific repair enzyme: the spore photoproduct lyase (SP lyase). In this work, a host-guest approach in which the N-terminal fragment of Moloney murine leukemia virus reverse transcriptase (MMLV RT) serves as the host and DNA as the guest was used to determine the crystal structures of complexes including 16 bp oligonucleotides with and without the SP lesion at 2.14 and 1.72 Å resolution, respectively. In contrast to other types of thymine-thymine lesions, the SP lesion retains normal Watson-Crick hydrogen bonding to the adenine bases of the complementary strand, with shorter hydrogen bonds than found in the structure of the undamaged DNA. However, the lesion induces structural changes in the local conformation of what is otherwise B-form DNA. The region surrounding the lesion differs significantly in helical form from B-DNA, and the minor groove is widened by almost 3 Å compared with that of the undamaged DNA. Thus, these unusual structural features associated with SP lesions may provide a basis for recognition by the SP lyase.

Weak operator binding enhances simulated Lac repressor-mediated DNA looping

The 50th anniversary of Biopolymers coincides closely with the like celebration of the discovery of the Escherichia coli (lac) lactose operon, a classic genetic system long used to illustrate the influence of biomolecular structure on function. The looping of DNA induced by the binding of the Lac repressor protein to sequentially distant operator sites on DNA continues to serve as a paradigm for understanding long-range genomic communication. Advances in analyses of DNA structures and in incorporation of proteins in computer simulations of DNA looping allow us to address long-standing questions about the role of protein-mediated DNA loop formation in transcriptional control. Here we report insights gained from studies of the sequence-dependent contributions of the natural lac operators to Lac repressor-mediated DNA looping. Novel superposition of the ensembles of protein-bound operator structures derived from NMR measurements reveals variations in DNA folding missed in conventional structural alignments. The changes in folding affect the predicted ease with which the repressor induces loop formation and the ways that DNA closes between the protein headpieces. The peeling of the auxiliary operators away from the repressor enhances the formation of loops with the 92-bp wildtype spacing and hints of a structural reason behind their weak binding.