Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins

HIPPO: HIstogram-based Pseudo-POtential for scoring protein-ssRNA fragment-based docking poses

BACKGROUND

The RNA-Recognition motif (RRM) is a protein domain that binds single-stranded RNA (ssRNA) and is present in as much as 2% of the human genome. Despite this important role in biology, RRM-ssRNA interactions are very challenging to study on the structural level because of the remarkable flexibility of ssRNA. In the absence of atomic-level experimental data, the only method able to predict the 3D structure of protein-ssRNA complexes with any degree of accuracy is ssRNA'TTRACT, an ssRNA fragment-based docking approach using ATTRACT. However, since ATTRACT parameters are not ssRNA-specific and were determined in 2010, there is substantial opportunity for enhancement.

RESULTS

Here we present HIPPO, a composite RRM-ssRNA scoring potential derived analytically from contact frequencies in near-native versus non-native docking models. HIPPO consists of a consensus of four distinct potentials, each extracted from a distinct reference pool of protein-trinucleotide docking decoys. To score a docking pose with one potential, for each pair of RNA-protein coarse-grained bead types, each contact is awarded or penalised according to the relative frequencies of this contact distance range among the correct and incorrect poses of the reference pool. Validated on a fragment-based docking benchmark of 57 experimentally solved RRM-ssRNA complexes, HIPPO achieved a threefold or higher enrichment for half of the fragments, versus only a quarter with the ATTRACT scoring function. In particular, HIPPO drastically improved the chance of very high enrichment (12-fold or higher), a scenario where the incremental modelling of entire ssRNA chains from fragments becomes viable. However, for the latter result, more research is needed to make it directly practically applicable. Regardless, our approach already improves upon the state of the art in RRM-ssRNA modelling and is in principle extendable to other types of protein-nucleic acid interactions.

Rapid Long-distance Migration of RPA on Single Stranded DNA Occurs Through Intersegmental Transfer Utilizing Multivalent Interactions

Replication Protein A (RPA) is asingle strandedDNA(ssDNA)binding protein that coordinates diverse DNA metabolic processes including DNA replication, repair, and recombination. RPA is a heterotrimeric protein with six functional oligosaccharide/oligonucleotide (OB) domains and flexible linkers. Flexibility enables RPA to adopt multiple configurations andis thought to modulate its function. Here, usingsingle moleculeconfocal fluorescencemicroscopy combinedwith optical tweezers and coarse-grained molecular dynamics simulations, we investigated the diffusional migration of single RPA molecules on ssDNA undertension.The diffusioncoefficientDis the highest (20,000nucleotides2/s) at 3pNtension and in 100 mMKCl and markedly decreases whentensionor salt concentrationincreases. We attribute the tension effect to intersegmental transfer which is hindered by DNA stretching and the salt effect to an increase in binding site size and interaction energy of RPA-ssDNA. Our integrative study allowed us to estimate the size and frequency of intersegmental transfer events that occur through transient bridging of distant sites on DNA by multiple binding sites on RPA. Interestingly, deletion of RPA trimeric core still allowed significant ssDNA binding although the reduced contact area made RPA 15-fold more mobile. Finally, we characterized the effect of RPA crowding on RPA migration. These findings reveal how the high affinity RPA-ssDNA interactions are remodeled to yield access, a key step in several DNA metabolic processes.

DNA flexibility can shape the preferential hypermutation of antibody genes

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Single molecule technique unveils the role of electrostatic interactions in ssDNA-gp32 molecular complex stability

The exploration of single-strand DNA-binding protein (SSB)-ssDNA interactions and their crucial roles in essential biological processes lagged behind other types of protein-nucleic acid interactions, such as protein-dsDNA and protein-RNA interactions. The ssDNA binding protein gene product 32 (gp32) of the T4 bacteriophage is a central integrating component of the replication complex that must continuously bind to and unbind from transiently exposed template strands during the DNA synthesis. To gain deeper insights into the electrostatic conditions influencing the stability of the ssDNA-gp32 molecular complex, like the salt concentration or some metal ions proven to specifically bind to gp32, we employed a method that performs rapid measurements of the DNA-protein stability using an α-Hemolysin (α-HL) protein nanopore. We indirectly probed the stability of a protein-nucleic acid complex by monitoring the dissociation process between the gp32 protein and the ssDNA molecular complex in single-molecular electrophysiology experiments, but also through fluorescence spectroscopy techniques. We have shown that the complex is more stable in 0.5 M KCl solution than in 2 M KCl solution and that the presence of Zn2+ ions further increases this stability for any salt used in the present study. This method can be applied to other nucleic acid-protein molecular complexes, as well as for an accurate determination of the drug-protein carrier stability.

Heterogeneity in Neutrophil Extracellular Traps from Healthy Human Subjects

Neutrophil extracellular traps (NETs), a key component of early defense against microbial infection, are also associated with tissue injury. NET composition has been reported to vary with some disease states, but the composition and variability of NETs across many healthy subjects provide a critical comparison that has not been well investigated. We evaluated NETs from twelve healthy subjects of varying ages isolated from multiple blood draws over a three-and-one-half-year period to delineate the variability in extracellular DNA, protein, enzymatic activities, and susceptibility to protease inhibitors. We calculated correlations for NET constituents and loss of human bronchial epithelial barrier integrity, measured by transepithelial electrical resistance, after NET exposure. We found that although there was some variability within the same subject over time, the mean NET total DNA, dsDNA, protein, LDH, neutrophil elastase (NE), and proteinase 3 (PR3) in isolated NETs were consistent across subjects. NET serine protease activity varied considerably within the same donor from day to day. The mean NET cathepsin G and MPO were significantly different across donors. IL-8 > IL-1RA > G-CSF were the most abundant cytokines in NETs. There was no significant difference in the mean concentration or variability of IL-8, IL-1RA, G-CSF, IL-1α, IL-1β, or TNF-α in different subjects' NETs. NET DNA concentration was correlated with increased NET neutrophil elastase activity and higher NET IL-1RA concentrations. The mean reduction in protease activity by protease inhibitors was significantly different across donors. NET DNA concentration correlated best with reductions in the barrier integrity of human bronchial epithelia. Defining NET concentration by DNA content correlates with other NET components and reductions in NET-driven epithelial barrier dysfunction, suggesting DNA is a reasonable surrogate measurement for these complex structures in healthy subjects.

Peptides designed from a bacteriophage capsid protein function as synthetic transcription repressors

The bacteriophage capsid protein, Psu (polarity suppression), inhibits the bacterial transcription terminator, Rho. In an effort to find nontraditional antibacterial agents, we previously designed peptides from the Psu C terminus that function as inhibitors of Rho. Here, we demonstrated that these peptides have positive surface-charge densities, and they downregulate many genes in Escherichia coli. We hypothesized that these peptides could bind to nucleic acids and repress gene expression. One of these peptides, peptide 33, represses in vitro transcription from the T7A1 and Plac promoters efficiently by blocking the access of RNA polymerase to the promoter, a mode of transcription repression akin to many bacterial repressors. In vivo, expressions of the peptides reduce the total RNA level as well as transcription from Plac and Posm promoters significantly. However, they are less efficient in repressing transcription from the rRNA promoters with a very high turnover of RNA polymerase. The peptide 33 binds to both single and dsDNA as well as to RNA with dissociation constants ranging from 1 to 5 μM exhibiting preferences for the single-stranded DNA and RNAs. These interactions are salt-resistant and not sequence-specific. Interactions with dsDNA are entropy-driven, while it is enthalpy-driven for the ssDNA. This mode of interaction with nucleic acids is similar to many nonspecific ssDNA-binding proteins. Expression of peptide 33 induces cell elongation and impaired cell division, possibly due to the dislodging of the DNA-binding proteins. Overall, we surmised that these synthetic transcription repressors would function like bacterial nucleoid-associated proteins.

Heterogeneity in Neutrophil Extracellular Traps from Healthy Human Subjects

Neutrophil Extracellular Traps (NETs), a key component of early defense against microbial infection, are also associated with tissue injury. NET composition has been reported to vary with some disease states, but the composition and variability of NETs across many healthy subjects provides a critical comparison that has not been well investigated. We evaluated NETs from twelve healthy subjects of varying ages isolated from multiple blood draws over a three and one half-year period to delineate the variability in extracellular DNA, protein, enzymatic activities, and susceptibility to protease inhibitors. We calculated correlations for NET constituents and loss of human bronchial epithelial barrier integrity, measured by transepithelial electrical resistance, after NET exposure. We found that although there was some variability within the same subject over time, the mean numbers of neutrophils, protein, LDH, serine protease activities, and cytokines IL-8, IL-1RA, and G-CSF in isolated NETs were consistent across subjects. Total DNA and double stranded DNA content in NETs were different across donors. NETs had little or no TNFα, IL-17A, or GM-CSF. NET DNA concentration correlated with increased NET neutrophil elastase activity and higher NET IL-1RA concentrations. NET serine protease activity varied considerably within the same donor from day-to-day. Mean response to protease inhibitors was significantly different across donors. NET DNA concentration correlated best with reductions in barrier integrity of human bronchial epithelia. Defining NET concentration by DNA content correlates with other NET components and reductions in NET-driven epithelial barrier dysfunction, suggesting DNA is a reasonable surrogate measurement for these complex structures in healthy subjects.

Analysis and Prediction of Pathogen Nucleic Acid Specificity for Toll-like Receptors in Vertebrates

Identification of key sequence, expression and function related features of nucleic acid-sensing host proteins is of fundamental importance to understand the dynamics of pathogen-specific host responses. To meet this objective, we considered toll-like receptors (TLRs), a representative class of membrane-bound sensor proteins, from 17 vertebrate species covering mammals, birds, reptiles, amphibians, and fishes in this comparative study. We identified the molecular signatures of host TLRs that are responsible for sensing pathogen nucleic acids or other pathogen-associated molecular patterns (PAMPs), and potentially play important roles in host defence mechanism. Interestingly, our findings reveal that such host-specific features are directly related to the strand (single or double) specificity of nucleic acid from pathogens. However, during host-pathogen interactions, such features were unable to explain the pathogenic PAMP (i.e., DNA, RNA or other) selectivity, suggesting a more complex mechanism. Using these features, we developed a number of machine learning models, of which Random Forest achieved a high performance (94.57% accuracy) to predict strand specificity of TLRs from protein-derived features. We applied the trained model to propose strand specificity of some previously uncharacterized distinct fish-specific novel TLRs (TLR18, TLR23, TLR24, TLR25, TLR27).

β-Sitosterol could serve as a dual inhibitor of Trypanosoma congolense sialidase and phospholipase A2: in vitro kinetic analyses and molecular dynamic simulations

The involvement of Trypanosoma congolense sialidase alongside phospholipase A2 has been widely accepted as the major contributing factor to anemia during African animal trypanosomiasis. The enzymes aid the parasite in scavenging sialic acid and fatty acids necessary for survival in the infected host, but there are no specific drug candidates against the two enzymes. This study investigated the inhibitory effects of β-sitosterol on the partially purified T. congolense sialidase and phospholipase A2. Purification of the enzymes using DEAE cellulose column led to fractions with highest specific activities of 8016.41 and 39.26 µmol/min/mg for sialidase and phospholipase A2, respectively. Inhibition kinetics studies showed that β-sitosterol is non-competitive and an uncompetitive inhibitor of sialidase and phospholipase A2 with inhibition binding constants of 0.368 and 0.549 µM, respectively. Molecular docking of the compound revealed binding energies of - 8.0 and - 8.6 kcal/mol against the sialidase and phospholipase A2, respectively. Furthermore, 100 ns molecular dynamics simulation using GROMACS revealed stable interaction of β-sitosterol with both enzymes. Hydrogen bond interactions between the ligand and Glu284 and Leu102 residues of the sialidase and phospholipase A2, respectively, were found to be the major stabilizing forces. In conclusion, β-sitosterol could serve as a dual inhibitor of T. congolense sialidase and phospholipase A2; hence, the compound could be exploited further in the search for newer trypanocides.

Targeting RNA Structure to Inhibit Editing in Trypanosomes

Mitochondrial RNA editing in trypanosomes represents an attractive target for developing safer and more efficient drugs for treating infections with trypanosomes because this RNA editing pathway is not found in humans. Other workers have targeted several enzymes in this editing system, but not the RNA. Here, we target a universal domain of the RNA editing substrate, which is the U-helix formed between the oligo-U tail of the guide RNA and the target mRNA. We selected a part of the U-helix that is rich in G-U wobble base pairs as the target site for the virtual screening of 262,000 compounds. After chemoinformatic filtering of the top 5000 leads, we subjected 50 representative complexes to 50 nanoseconds of molecular dynamics simulations. We identified 15 compounds that retained stable interactions in the deep groove of the U-helix. The microscale thermophoresis binding experiments on these five compounds show low-micromolar to nanomolar binding affinities. The UV melting studies show an increase in the melting temperatures of the U-helix upon binding by each compound. These five compounds can serve as leads for drug development and as research tools to probe the role of the RNA structure in trypanosomal RNA editing.

Transcriptome-wide identification of single-stranded RNA binding proteins

RNA-protein interactions are precisely regulated by RNA secondary structures in various biological processes. Large-scale identification of proteins that interact with particular RNA structure is important to the RBPome. Herein, a kethoxal assisted single-stranded RNA interactome capture (KASRIC) strategy was developed to globally identify single-stranded RNA binding proteins (ssRBPs). This approach combines RNA secondary structure probing technology with the conventional method of RNA-binding proteins profiling, realizing the transcriptome-wide identification of ssRBPs. Applying KASRIC, we identified 3180 candidate RBPs and 244 candidate ssRBPs in HeLa cells. Importantly, the 244 candidate ssRBPs contained 55 previously reported ssRBPs and 189 novel ssRBPs. Function analysis of the candidate ssRBPs exhibited enrichment in cellular processes related to RNA splicing and RNA degradation. The KASRIC strategy will facilitate the investigation of RNA-protein interactions.

Thymine DNA Glycosylase is an RNA-Binding Protein with High Selectivity for G-Rich Sequences

Thymine DNA glycosylase (TDG) is a multifaceted enzyme involved in several critical biological pathways, including transcriptional activation, DNA demethylation, and DNA repair. Recent studies have established regulatory relationships between TDG and RNA, but the molecular interactions underlying these relationships is poorly understood. Herein, we now demonstrate that TDG binds directly to RNA with nanomolar affinity. Using synthetic oligonucleotides of defined length and sequence, we show that TDG has a strong preference for binding G-rich sequences in single-stranded RNA but binds weakly to single-stranded DNA and duplex RNA. TDG also binds tightly to endogenous RNA sequences. Studies with truncated proteins indicate that TDG binds RNA primarily through its structured catalytic domain and that its disordered C-terminal domain plays a key role in regulating TDG's affinity and selectivity for RNA. Finally, we show that RNA competes with DNA for binding to TDG, resulting in inhibition of TDG-mediated excision in the presence of RNA. Together, this work provides support for and insights into a mechanism wherein TDG-mediated processes (e.g., DNA demethylation) are regulated through the direct interactions of TDG with RNA.

Influence of diverse storage conditions of double-stranded RNA in vitro on the RNA interference efficiency in vivo insect Tribolium castaneum

BACKGROUND

A significant variation in RNA interference (RNAi) efficiency hinders further functional gene studies and pest control application in many insects. The available double-stranded RNA (dsRNA) molecules introduced into the target cells are regarded as the crucial factor for efficient RNAi response. However, numerous studies have only focused on dsRNA stability in vivo; it is uncertain whether different dsRNA storage conditions in vitro play a role in variable RNAi efficiency among insects.

RESULTS

A marker gene cardinal, which leads to white eyes when knocked-down in the red flour beetle Tribolium castaneum, was used to evaluate the effects of RNAi efficiency under different dsRNA storage conditions. We demonstrated that the dsRNA molecule is very stable under typical cryopreservation temperatures (-80 and -20 °C) within 180 days, and RNAi efficiency shows no significant differences under either low temperature. Unexpectedly, while dsRNA molecules were treated with multiple freeze-thaw cycles up to 50 times between -80/-20 °C and room temperature, we discovered that dsRNA integrity and RNAi efficiency were comparable with fresh dsRNA. Finally, when the stability of dsRNA was further measured under refrigerated storage conditions (4 °C), we surprisingly found that dsRNA is still stable within 180 days and can induce an efficient RNAi response as that of initial dsRNA.

CONCLUSION

Our results indicate that dsRNA is extraordinarily stable under various temperature storage conditions that did not significantly impact RNAi efficiency in vivo insects. © 2022 Society of Chemical Industry.

Interdomain dynamics in human Replication Protein A regulates kinetics and thermodynamics of its binding to ssDNA

Human Replication Protein A (hRPA) is a multidomain protein that interacts with ssDNA intermediates to provide the latter much-needed stability during DNA metabolism and maintain genomic integrity. Although the ssDNA organization with hRPA was studied recently through experimental means, characterizing the underlying mechanism at the atomic level remains challenging because of the dynamic domain architecture of hRPA and poorly understood heterogeneity of ssDNA-protein interactions. Here, we used a computational framework, precisely tailored to capture protein-ssDNA interactions, and investigated the binding of hRPA with a 60 nt ssDNA. Two distinct binding mechanisms are realized based on the hRPA domain flexibility. For a rigid domain architecture of hRPA, ssDNA binds sequentially with hRPA domains, resulting in slow association kinetics. The binding pathway involves the formation of stable and distinct intermediate states. On contrary, for a flexible domain architecture of hRPA, ssDNA binds synergistically to the A and B domains followed by the rest of hRPA. The domain dynamics in hRPA alleviates the free energy cost of domain orientation necessary for specific binding with ssDNA, leading to fast association kinetics along a downhill binding free energy landscape. An ensemble of free energetically degenerate intermediate states is encountered that makes it arduous to characterize them structurally. An excellent match between our results with the available experimental observations provides new insights into the rich dynamics of hRPA binding to ssDNA and in general paves the way to investigate intricate details of ssDNA-protein interactions, crucial for cellular functioning.

Structure-Guided Design of a Potent and Specific Inhibitor against the Genomic Mutator APOBEC3A

Nucleic acid structure plays a critical role in governing the selectivity of DNA- and RNA-modifying enzymes. In the case of the APOBEC3 family of cytidine deaminases, these enzymes catalyze the conversion of cytosine (C) to uracil (U) in single-stranded DNA, primarily in the context of innate immunity. DNA deamination can also have pathological consequences, accelerating the evolution of viral genomes or, when the host genome is targeted by either APOBEC3A (A3A) or APOBEC3B (A3B), promoting tumor evolution leading to worse patient prognosis and chemotherapeutic resistance. For A3A, nucleic acid secondary structure has emerged as a critical determinant of substrate targeting, with a predilection for DNA that can form stem loop hairpins. Here, we report the development of a specific nanomolar-level, nucleic acid-based inhibitor of A3A. Our strategy relies on embedding the nucleobase 5-methylzebularine, a mechanism-based inhibitor, into a DNA dumbbell structure, which mimics the ideal substrate secondary structure for A3A. Structure-activity relationship studies using a panel of diverse inhibitors reveal a critical role for the stem and position of the inhibitor moiety in achieving potent inhibition. Moreover, we demonstrate that DNA dumbbell inhibitors, but not nonstructured inhibitors, show specificity against A3A relative to the closely related catalytic domain of A3B. Overall, our work demonstrates the feasibility of leveraging secondary structural preferences in inhibitor design, offering a blueprint for further development of modulators of DNA-modifying enzymes and potential therapeutics to circumvent APOBEC-driven viral and tumor evolution.

Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

Harnessing nucleic acid technologies for human health on earth and in space

Nucleic acid therapeutics are a versatile class of sequence-programmable drugs that offer a robust and clinically viable strategy to modulate expression or correct genetic defects contributing to disease. The majority of drugs currently on the market target proteins; however, proteins only represent a subset of possible disease targets. Nucleic acid therapeutics allow intuitive engagement with genome sequences providing a more direct way to target many diseases at their genetic root cause. Their clinical success depends on platform technologies which can support durable and well tolerated pharmacological activity in a given tissue. Nucleic acid drugs possess a potent combination of target specificity and adaptability required to advance drug development for many diseases. As these therapeutic technologies mature, their clinical applications can also expand access to personalized therapies for patients with rare or solo genetic diseases. Spaceflight crew members exposed to the unique hazards of spaceflight, especially those related to galactic cosmic radiation (GCR) exposure, represent another patient subset who may also benefit from nucleic acid drugs as countermeasures. In this review, we will discuss the various classes of RNA- and DNA-targeted nucleic acid drugs, provide an overview of their present-day clinical applications, and describe major strategies to improve their delivery, safety, and overall efficacy.

Chloride Enhances DNA Reactivity with Chlorine under Conditions Relevant to Water Treatment

Free available chlorine (FAC) is widely used to inactivate viruses by oxidizing viral components, including genomes. It is commonly assumed that hypochlorous acid (HOCl) is the chlorinating agent responsible for virus inactivation; however, recent studies have underscored that minor constituents of FAC existing in equilibrium with HOCl, such as molecular chlorine (Cl₂), can influence FAC reactivity toward select organic compounds. This study measures the FAC reaction kinetics with dsDNA and ssDNA extracted from representative bacteriophages (T3 and ϕX174) in samples augmented with chloride. Herein, chloride enhances FAC reactivity toward dsDNA and, to a lesser extent, toward ssDNA, especially at pH < 7.5. The enhanced reactivity can be attributed to the formation of Cl₂. Second-order rate constants were determined for reactions of ssDNA and dsDNA with HOCl and Cl₂. DNA chlorination kinetics followed the reactivity-selectivity principle, where the more-reactive nucleophilic species (ssDNA, ∼100× more reactive than dsDNA) reacted less selectively with electrophilic FAC species. The addition of chloride was also shown to enhance the inactivation of bacteriophage T3 (dsDNA genome) by FAC but did not enhance the inactivation of bacteriophage ϕX174 (ssDNA genome). Overall, the results suggest that Cl₂ is an important chlorinating agent of nucleic acids and viruses.

Mechanistic insight into the RNA stimulated-ATPase activity of tick-borne encephalitis virus helicase

The helicase domain of nonstructural protein 3 (NS3H) unwinds the double-stranded RNA replication intermediate in an ATP-dependent manner during the flavivirus life cycle. While the ATP hydrolysis mechanism of Dengue and Zika viruses NS3H has been extensively studied, little is known in the case of the tick-borne encephalitis virus NS3H. We demonstrate that ssRNA binds with nanomolar affinity to NS3H and strongly stimulates the ATP hydrolysis cycle, whereas ssDNA binds only weakly and inhibits ATPase activity in a noncompetitive manner. Thus, NS3H is an RNA-specific helicase, whereas DNA might act as an allosteric inhibitor. Using modeling, we explored plausible allosteric mechanisms by which ssDNA inhibits the ATPase via nonspecific binding in the vicinity of the active site and ATP repositioning. We captured several structural snapshots of key ATP hydrolysis stages using X-ray crystallography. One intermediate, in which the inorganic phosphate and ADP remained trapped inside the ATPase site after hydrolysis, suggests that inorganic phosphate release is the rate-limiting step. Using structure-guided modeling and molecular dynamics simulation, we identified putative RNA-binding residues and observed that the opening and closing of the ATP-binding site modulates RNA affinity. Site-directed mutagenesis of the conserved RNA-binding residues revealed that the allosteric activation of ATPase activity is primarily communicated via an arginine residue in domain 1. In summary, we characterized conformational changes associated with modulating RNA affinity and mapped allosteric communication between RNA-binding groove and ATPase site of tick-borne encephalitis virus helicase.

Long-range electrostatic interactions significantly modulate the affinity of dynein for microtubules

The dynein family of microtubule minus-end-directed motor proteins drives diverse functions in eukaryotic cells, including cell division, intracellular transport, and flagellar beating. Motor protein processivity, which characterizes how far a motor walks before detaching from its filament, depends on the interaction between its microtubule-binding domain (MTBD) and the microtubule. Dynein's MTBD switches between high- and low-binding affinity states as it steps. Significant structural and functional data show that specific salt bridges within the MTBD and between the MTBD and the microtubule govern these affinity state shifts. However, recent computational work suggests that nonspecific, long-range electrostatic interactions between the MTBD and the microtubule may also play an important role in the processivity of dynein. To investigate this hypothesis, we mutated negatively charged amino acids remote from the dynein MTBD-microtubule-binding interface to neutral residues and measured the binding affinity using microscale thermophoresis and optical tweezers. We found a significant increase in the binding affinity of the mutated MTBDs for microtubules. Furthermore, we found that charge screening by free ions in solution differentially affected the binding and unbinding rates of MTBDs to microtubules. Together, these results demonstrate a significant role for long-range electrostatic interactions in regulating dynein-microtubule affinity. Moreover, these results provide insight into the principles that potentially underlie the biophysical differences between molecular motors with various processivities and protein-protein interactions more generally.

Single-molecule biophysics experiments in silico: Toward a physical model of a replisome

The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments outright with computational mimics. Here, we offer a blueprint for performing single-molecule studies in silico using a DNA-binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic single-molecule experiments. We apply the model to recapitulate, in silico, force-extension characterization of protein binding to single-stranded DNA and protein and DNA replacement assays, providing a detailed portrait of the underlying mechanics. Finally, we use the model to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.

•

Coarse-grained model derived from all-atom simulation recapitulates experiments

•

Model reproduces the elastic response to force and exchange dynamics

•

Model reveals structure of intermediate states usually inaccessible to experiment

•

Model applied to viral replisome with trombone loop containing tens of SSB proteins

Bicyclo-DNA mimics with enhanced protein binding affinities: insights from molecular dynamics simulations

DNA-protein interactions occur at all levels of DNA expression and replication and are crucial determinants for the survival of a cell. Several modified nucleotides have been utilized to manipulate these interactions and have implications in drug discovery. In the present article, we evaluated the binding of bicyclo-nucleotides (generated by forming a methylene bridge between C1' and C5' in sugar, leading to a bicyclo system with C2' axis of symmetry at the nucleotide level) to proteins. We utilized four ssDNA-protein complexes with experimentally known binding free energies and investigated the binding of modified nucleotides to proteins via all-atom explicit solvent molecular dynamics (MD) simulations (200 ns), and compared the binding with control ssDNA-protein systems. The modified ssDNA displayed enhanced binding to proteins as compared to the control ssDNA, as seen by means of MD simulations followed by MM-PBSA calculations. Further, the Delphi-based electrostatic estimation revealed that the high binding of modified ssDNA to protein might be related to the enhanced electrostatic complementarity displayed by the modified ssDNA molecules in all the four systems considered for the study. The improved binding achieved with modified nucleotides can be utilized to design and develop anticancer/antisense molecules capable of targeting proteins or ssRNAs.Communicated by Ramaswamy H. Sarma.

Evaluation of FRET X for single-molecule protein fingerprinting

Single-molecule protein identification is an unrealized concept with potentially ground-breaking applications in biological research. We propose a method called FRET X (Förster Resonance Energy Transfer via DNA eXchange) fingerprinting, in which the FRET efficiency is read out between exchangeable dyes on protein-bound DNA docking strands and accumulated FRET efficiencies constitute the fingerprint for a protein. To evaluate the feasibility of this approach, we simulated fingerprints for hundreds of proteins using a coarse-grained lattice model and experimentally demonstrated FRET X fingerprinting on model peptides. Measured fingerprints are in agreement with our simulations, corroborating the validity of our modeling approach. In a simulated complex mixture of >300 human proteins of which only cysteines, lysines, and arginines were labeled, a support vector machine was able to identify constituents with 95% accuracy. We anticipate that our FRET X fingerprinting approach will form the basis of an analysis tool for targeted proteomics.

•

We propose a FRET-based single-molecule protein identification method

•

Peptides are experimentally distinguishable by their fingerprints

•

Our approach can classify the constituents of complex samples with 95% accuracy

Conformational variability in proteins bound to single-stranded DNA: A new benchmark for new docking perspectives

We explored the Protein Data Bank (PDB) to collect protein-ssDNA structures and create a multi-conformational docking benchmark including both bound and unbound protein structures. Due to ssDNA high flexibility when not bound, no ssDNA unbound structure is included in the benchmark. For the 91 sequence-identity groups identified as bound-unbound structures of the same protein, we studied the conformational changes in the protein induced by the ssDNA binding. Moreover, based on several bound or unbound protein structures in some groups, we also assessed the intrinsic conformational variability in either bound or unbound conditions and compared it to the supposedly binding-induced modifications. To illustrate a use case of this benchmark, we performed docking experiments using ATTRACT docking software. This benchmark is, to our knowledge, the first one made to peruse available structures of ssDNA-protein interactions to such an extent, aiming to improve computational docking tools dedicated to this kind of molecular interactions.

Destruction of DNA-Binding Proteins by Programmable Oligonucleotide PROTAC (O'PROTAC): Effective Targeting of LEF1 and ERG

DNA-binding proteins, including transcription factors (TFs), play essential roles in various cellular processes and pathogenesis of diseases, deeming to be potential therapeutic targets. However, these proteins are generally considered undruggable as they lack an enzymatic catalytic site or a ligand-binding pocket. Proteolysis-targeting chimera (PROTAC) technology has been developed by engineering a bifunctional molecule chimera to bring a protein of interest (POI) to the proximity of an E3 ubiquitin ligase, thus inducing the ubiquitination of POI and further degradation through the proteasome pathway. Here, the development of oligonucleotide-based PROTAC (O'PROTACs), a class of noncanonical PROTACs in which a TF-recognizing double-stranded oligonucleotide is incorporated as a binding moiety of POI is reported. It is demonstrated that O'PROTACs of lymphoid enhancer-binding factor 1 (LEF1) and ETS-related gene (ERG), two highly cancer-related transcription factors, successfully promote degradation of these proteins, impede their transcriptional activity, and inhibit cancer cell growth in vitro and in vivo. The programmable nature of O'PROTACs indicates that this approach is also applicable to destruct other TFs. O'PROTACs not only can serve as a research tool but also can be harnessed as a therapeutic arsenal to target DNA binding proteins for effective treatment of diseases such as cancer.

Qki activates Srebp2-mediated cholesterol biosynthesis for maintenance of eye lens transparency

Defective cholesterol biosynthesis in eye lens cells is often associated with cataracts; however, how genes involved in cholesterol biosynthesis are regulated in lens cells remains unclear. Here, we show that Quaking (Qki) is required for the transcriptional activation of genes involved in cholesterol biosynthesis in the eye lens. At the transcriptome level, lens-specific Qki-deficient mice present downregulation of genes associated with the cholesterol biosynthesis pathway, resulting in a significant reduction of total cholesterol level in the eye lens. Mice with Qki depletion in lens epithelium display progressive accumulation of protein aggregates, eventually leading to cataracts. Notably, these defects are attenuated by topical sterol administration. Mechanistically, we demonstrate that Qki enhances cholesterol biosynthesis by recruiting Srebp2 and Pol II in the promoter regions of cholesterol biosynthesis genes. Supporting its function as a transcription co-activator, we show that Qki directly interacts with single-stranded DNA. In conclusion, we propose that Qki-Srebp2-mediated cholesterol biosynthesis is essential for maintaining the cholesterol level that protects lens from cataract development.

A comparative study of protein-ssDNA interactions

Single-stranded DNA-binding proteins (SSBs) play crucial roles in DNA replication, recombination and repair, and serve as key players in the maintenance of genomic stability. While a number of SSBs bind single-stranded DNA (ssDNA) non-specifically, the others recognize and bind specific ssDNA sequences. The mechanisms underlying this binding discrepancy, however, are largely unknown. Here, we present a comparative study of protein-ssDNA interactions by annotating specific and non-specific SSBs and comparing structural features such as DNA-binding propensities and secondary structure types of residues in SSB-ssDNA interactions, protein-ssDNA hydrogen bonding and π-π interactions between specific and non-specific SSBs. Our results suggest that protein side chain-DNA base hydrogen bonds are the major contributors to protein-ssDNA binding specificity, while π-π interactions may mainly contribute to binding affinity. We also found the enrichment of aspartate in the specific SSBs, a key feature in specific protein-double-stranded DNA (dsDNA) interactions as reported in our previous study. In addition, no significant differences between specific and non-specific groups with respect of conformational changes upon ssDNA binding were found, suggesting that the flexibility of SSBs plays a lesser role than that of dsDNA-binding proteins in conferring binding specificity.

Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer

Long noncoding RNAs (lncRNAs) influence cellular function through binding events that often depend on the lncRNA secondary structure. One such lncRNA, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), is upregulated in many cancer types and has a myriad of protein- and miRNA-binding sites. Recently, a secondary structural model of MALAT1 in noncancerous cells was proposed to form 194 hairpins and 13 pseudoknots. That study postulated that, in cancer cells, the MALAT1 structure likely varies, thereby influencing cancer progression. This work analyzes how that structural model is expected to change in K562 cells, which originated from a patient with chronic myeloid leukemia (CML), and in HeLa cells, which originated from a patient with cervical cancer. Dimethyl sulfate-sequencing (DMS-Seq) data from K562 cells and psoralen analysis of RNA interactions and structure (PARIS) data from HeLa cells were compared to the working structural model of MALAT1 in noncancerous cells to identify sites that likely undergo structural alterations. MALAT1 in K562 cells is predicted to become more unstructured, with almost 60% of examined hairpins in noncancerous cells losing at least half of their base pairings. Conversely, MALAT1 in HeLa cells is predicted to largely maintain its structure, undergoing 18 novel structural rearrangements. Moreover, 50 validated miRNA-binding sites are affected by putative secondary structural changes in both cancer types, such as miR-217 in K562 cells and miR-20a in HeLa cells. Structural changes unique to K562 cells and HeLa cells provide new mechanistic leads into how the structure of MALAT1 may mediate cancer in a cell-type specific manner.

Genome Editing by CRISPR-Cas: A Game Change in the Genetic Manipulation of Chlamydomonas

Microalgae are promising photosynthetic unicellular eukaryotes among the most abundant on the planet and are considered as alternative sustainable resources for various industrial applications. Chlamydomonas is an emerging model for microalgae to be manipulated by multiple biotechnological tools in order to produce high-value bioproducts such as biofuels, bioactive peptides, pigments, nutraceuticals, and medicines. Specifically, Chlamydomonas reinhardtii has become a subject of different genetic-editing techniques adapted to modulate the production of microalgal metabolites. The main nuclear genome-editing tools available today include zinc finger nucleases (ZFNs), transcriptional activator-like effector nucleases (TALENs), and more recently discovered the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR associated protein (Cas) nuclease system. The latter, shown to have an interesting editing capacity, has become an essential tool for genome editing. In this review, we highlight the available literature on the methods and the applications of CRISPR-Cas for C. reinhardtii genetic engineering, including recent transformation methods, most used bioinformatic tools, best strategies for the expression of Cas protein and sgRNA, the CRISPR-Cas mediated gene knock-in/knock-out strategies, and finally the literature related to CRISPR expression and modification approaches.

Mechanism of Dynamic Binding of Replication Protein A to ssDNA

Replication protein A (RPA) serves as a hub protein inside eukaryotic cells, where it coordinates crucial DNA metabolic processes and activates the DNA-damage response system. A characteristic feature of its action is to associate with single-stranded DNA (ssDNA) intermediates before handing them over to downstream proteins. The length of ssDNA intermediates differs for different pathways. This means that RPA must have mechanisms for selective processing of ssDNA intermediates based on their length, the knowledge of which is fundamental to elucidate when and how DNA repair and replication processes are symphonized. By employing extensive molecular dynamics simulations, we investigated the mechanism of binding of RPA to ssDNA of different lengths. We show that the binding involves dynamic equilibrium with a stable intermediate, the population of which increases with the length of ssDNA. The vital underlying factors are decoded through collective variable principal component analysis. It suggests a differently orchestrated set of interactions that define the action of RPA based on the length of ssDNA intermediates. We further estimated the association kinetics that matches excellently well with previous experimental studies and probed the diffusion mechanism of RPA to ssDNA. RPA diffuses on short ssDNA through progressive "bulge" formation. With long ssDNA, we observed a conformational change in ssDNA coupled with its binding to RPA in a cooperative fashion. This unanticipated binding mechanism successfully explains how the "short-lived", long ssDNA intermediates are processed quickly in vivo. This study thus reveals the molecular basis of several recent experimental observations related to RPA binding to ssDNA and provides novel insights into the RPA functioning in DNA repair and replication.

Open-tubular nanoelectrochromatography (OT-NEC): gel-free separation of single stranded DNAs (ssDNAs) in thermoplastic nanochannels

Electrophoresis or electrochromatography carried out in nanometer columns (width and depth) offers some attractive benefits compared to microscale columns. These advantages include unique separation mechanisms that are scale dependent, fast separation times, and simpler workflow due to the lack of a need for column packing and/or wall coatings to create a stationary phase. We report the use of thermoplastics, in this case PMMA, as the substrate for separating single-stranded DNAs (ssDNAs). Electrophoresis nanochannels were created in PMMA using nanoimprint lithography (NIL), which can produce devices at lower cost and in a higher production mode compared to the fabrication techniques required for glass devices. The nanochannel column in PMMA was successful in separating ssDNAs in free solution that was not possible using microchip electrophoresis in PMMA. The separation could be performed in <1 s with resolution >1.5 when carried out using at an electric field strength of 280 V/cm and an effective column length of 60 μm (100 nm × 100 nm, depth and width). The ssDNAs transport through the PMMA column was driven electrokinetically under the influence of an EOF. The results indicated that the separation was dominated by chromatographic effects using an open tubular nano-electrochromatography (OT-NEC) mode of separation. Interesting to these separations was that no column packing was required nor a wall coating to create the stationary phase; the separation was affected using the native polymer that was UV/O₃ activated and an aqueous buffer mobile phase.

The evolution of universal adaptations of life is driven by universal properties of matter: energy, entropy, and interaction

The evolution of multicellular eukaryotes expresses two sorts of adaptations: local adaptations like fur or feathers, which characterize species in particular environments, and universal adaptations like microbiomes or sexual reproduction, which characterize most multicellulars in any environment. We reason that the mechanisms driving the universal adaptations of multicellulars should themselves be universal, and propose a mechanism based on properties of matter and systems: energy, entropy, and interaction. Energy from the sun, earth and beyond creates new arrangements and interactions. Metabolic networks channel some of this energy to form cooperating, interactive arrangements. Entropy, used here as a term for all forces that dismantle ordered structures (rather than as a physical quantity), acts as a selective force. Entropy selects for arrangements that resist it long enough to replicate, and dismantles those that do not. Interactions, energy-charged and dynamic, restrain entropy and enable survival and propagation of integrated living systems. This fosters survival-of-the-fitted - those entities that resist entropic destruction - and not only of the fittest - the entities with the greatest reproductive success. The "unit" of evolution is not a discrete entity, such as a gene, individual, or species; what evolves are collections of related interactions at multiple scales. Survival-of-the-fitted explains universal adaptations, including resident microbiomes, sexual reproduction, continuous diversification, programmed turnover, seemingly wasteful phenotypes, altruism, co-evolving environmental niches, and advancing complexity. Indeed survival-of-the-fittest may be a particular case of the survival-of-the-fitted mechanism, promoting local adaptations that express reproductive advantages in addition to resisting entropy. Survival-of-the-fitted accounts for phenomena that have been attributed to neutral evolution: in the face of entropy, there is no neutrality; all variations are challenged by ubiquitous energy and entropy, retaining those that are "fit enough". We propose experiments to test predictions of the survival-of-the-fitted theory, and discuss implications for the wellbeing of humans and the biosphere.

The evolution of universal adaptations of life is driven by universal properties of matter: energy, entropy, and interaction

The evolution of multicellular eukaryotes expresses two sorts of adaptations: local adaptations like fur or feathers, which characterize species in particular environments, and universal adaptations like microbiomes or sexual reproduction, which characterize most multicellulars in any environment. We reason that the mechanisms driving the universal adaptations of multicellulars should themselves be universal, and propose a mechanism based on properties of matter and systems: energy, entropy, and interaction. Energy from the sun, earth and beyond creates new arrangements and interactions. Metabolic networks channel some of this energy to form cooperating, interactive arrangements. Entropy, used here as a term for all forces that dismantle ordered structures (rather than as a physical quantity), acts as a selective force. Entropy selects for arrangements that resist it long enough to replicate, and dismantles those that do not. Interactions, energy-charged and dynamic, restrain entropy and enable survival and propagation of integrated living systems. This fosters survival-of-the-fitted - those entities that resist entropic destruction - and not only of the fittest - the entities with the greatest reproductive success. The "unit" of evolution is not a discrete entity, such as a gene, individual, or species; what evolves are collections of related interactions at multiple scales. Survival-of-the-fitted explains universal adaptations, including resident microbiomes, sexual reproduction, continuous diversification, programmed turnover, seemingly wasteful phenotypes, altruism, co-evolving environmental niches, and advancing complexity. Indeed survival-of-the-fittest may be a particular case of the survival-of-the-fitted mechanism, promoting local adaptations that express reproductive advantages in addition to resisting entropy. Survival-of-the-fitted accounts for phenomena that have been attributed to neutral evolution: in the face of entropy, there is no neutrality; all variations are challenged by ubiquitous energy and entropy, retaining those that are "fit enough". We propose experiments to test predictions of the survival-of-the-fitted theory, and discuss implications for the wellbeing of humans and the biosphere.

The evolution of universal adaptations of life is driven by universal properties of matter: energy, entropy, and interaction

The evolution of multicellular eukaryotes expresses two sorts of adaptations: local adaptations like fur or feathers, which characterize species in particular environments, and universal adaptations like microbiomes or sexual reproduction, which characterize most multicellulars in any environment. We reason that the mechanisms driving the universal adaptations of multicellulars should themselves be universal, and propose a mechanism based on properties of matter and systems: energy, entropy, and interaction. Energy from the sun, earth and beyond creates new arrangements and interactions. Metabolic networks channel some of this energy to form cooperating, interactive arrangements. Entropy, used here as a term for all forces that dismantle ordered structures (rather than as a physical quantity), acts as a selective force. Entropy selects for arrangements that resist it long enough to replicate, and dismantles those that do not. Interactions, energy-charged and dynamic, restrain entropy and enable survival and propagation of integrated living systems. This fosters survival-of-the-fitted - those entities that resist entropic destruction - and not only of the fittest - the entities with the greatest reproductive success. The "unit" of evolution is not a discrete entity, such as a gene, individual, or species; what evolves are collections of related interactions at multiple scales. Survival-of-the-fitted explains universal adaptations, including resident microbiomes, sexual reproduction, continuous diversification, programmed turnover, seemingly wasteful phenotypes, altruism, co-evolving environmental niches, and advancing complexity. Indeed survival-of-the-fittest may be a particular case of the survival-of-the-fitted mechanism, promoting local adaptations that express reproductive advantages in addition to resisting entropy. Survival-of-the-fitted accounts for phenomena that have been attributed to neutral evolution: in the face of entropy, there is no neutrality; all variations are challenged by ubiquitous energy and entropy, retaining those that are "fit enough". We propose experiments to test predictions of the survival-of-the-fitted theory, and discuss implications for the wellbeing of humans and the biosphere.

Comprehensive Analysis of the Safety Profile of a Single-Stranded RNA Nano-Structure Adjuvant

Adjuvants enhance the efficacy of vaccines by stimulating immune response-related gene expression and pathways. Although some adjuvants have been approved for commercial use in human vaccines (e.g., Alum, MF59, and AS03), they might elicit adverse side effects, such as autoimmune diseases. Recently, we developed a novel single-stranded RNA (ssRNA) nano-structure adjuvant, which can stimulate both Th1 and Th2 responses. In this study, we evaluated the safety and toxicological profiles of this ssRNA nano-structure adjuvant in vitro and in vivo. Mice were intramuscularly immunized with the ssRNA nano-structure adjuvant three times, once every 2 weeks. The results indicate no significant differences in hematological and serum biochemistry parameters between the ssRNA-treated groups and the control group. From a histopathological perspective, no evidence of tissue damage was found in any group. The levels of IgE and anti-nuclear antibodies, which are markers of autoimmune disease, were not different between the ssRNA-treated groups and the control group. The findings of this study suggest that the ssRNA nano-structure can be used as a safe adjuvant to increase vaccine efficacies.