In-Cell Stability Prediction of RNA/DNA Hybrid Duplexes for Designing Oligonucleotides Aimed at Therapeutics.

In cells, the formation of RNA/DNA hybrid duplexes regulates gene expression and modification. The environment inside cellular organelles is heterogeneously crowded with high concentrations of biomolecules that affect the structure and stability of RNA/DNA hybrid duplexes. However, the detailed environmental effects remain unclear. Therefore, the mechanistic details of the effect of such molecular crowding were investigated at the molecular level by using thermodynamic and nuclear magnetic resonance analyses, revealing structure-dependent destabilization of the duplexes under crowded conditions. The transition from B- to A-like hybrid duplexes due to a change in conformation of the DNA strand guided by purine-pyrimidine asymmetry significantly increased the hydration number, which resulted in greater destabilization by the addition of cosolutes. By quantifying the individual contributions of environmental factors and the bulk structure of the duplex, we developed a set of parameters that predict the stability of hybrid duplexes with conformational dissimilarities under diverse crowding conditions. A comparison of the effects of environmental conditions in living cells and in vitro crowded solutions on hybrid duplex formation using the Förster resonance energy transfer technique established the applicability of our parameters to living cells. Moreover, our derived parameters can be used to estimate the efficiency of transcriptional inhibition, genome editing, and silencing techniques in cells. This supports the usefulness of our parameters for the visualization of cellular mechanisms of gene expression and the development of nucleic acid-based therapeutics targeting different cells.

Theranostic approach to specifically targeting the interloop region of BCL2 i-motif DNA by crystal violet.

Ligands that recognise specific i-motif DNAs are helpful in cancer diagnostics and therapeutics, as i-motif formation can cause cancer. Although the loop regions of i-motifs are promising targets for ligands, the interaction between a ligand and the loop regions based on sequence information remains unexplored. Herein, we investigated the loop regions of various i-motif DNAs to determine whether these regions specifically interact with fluorescent ligands. Crystal violet (CV), a triphenylmethane dye, exhibited strong fluorescence with the i-motif derived from the promoter region of the human BCL2 gene in a sequence- and structure-specific manner. Our systematic sequence analysis indicated that CV was bound to the site formed by the first and third loops through inter-loop interactions between the guanine bases present in these loops. As the structural stability of the BCL2 i-motif was unaffected by CV, the local stabilisation of the loops by CV could inhibit the interaction of transcription factors with these loops, repressing the BCL2 expression of MCF-7 cells. Our finding suggests that the loops of the i-motif can act as a novel platform for the specific binding of small molecules; thus, they could be utilised for the theranostics of diseases associated with i-motif DNAs.

Quantum Chemical Interaction Analysis between SARS-CoV-2 Main Protease and Ensitrelvir Compared with Its Initial Screening Hit.

A non-covalent oral drug targeting SARS-CoV-2 main protease (Mpro), ensitrelvir (Xocova), has been developed using structure-based drug design (SBDD). To elucidate the factors responsible for enhanced inhibitory activities from an in silico screening hit compound to ensitrelvir, we analyzed the interaction energies of the inhibitors with each residue of Mpro using fragment molecular orbital (FMO) calculations. This analysis reveals that functional group conversion for P1' and P1 parts in the inhibitors increases the strength of existing interactions with Mpro and also provides novel interactions for ensitrelvir; the associated changes in the conformation of Mpro induce further interactions for ensitrelvir in other parts, including hydrogen bonds, a halogen bond, and π-orbital interactions. Thus, we illuminate the promising strategies of SBDD for leading ensitrelvir to get higher activity against Mpro by elucidating microscopic interactions through FMO-based analysis. These detailed mechanism findings, including water cross-linkings, will help to design novel inhibitors in SBDD.

Nearest-neighbor parameters for the prediction of RNA duplex stability in diverse in vitro and cellular-like crowding conditions.

RNA performs various spatiotemporal functions in living cells. As the solution environments significantly affect the stability of RNA duplexes, a stability prediction of the RNA duplexes in diverse crowded conditions is required to understand and modulate gene expression in heterogeneously crowded intracellular conditions. Herein, we determined the nearest-neighbor (NN) parameters for RNA duplex formation when subjected to crowding conditions with an ionic concentration relevant to that found in cells. Determination of the individual contributions of excluded volume effect and water activity to each of the NN parameters in crowded environments enabled prediction of the thermodynamic parameters and their melting temperatures for plenty of tested RNA duplex formation in vitro and in cell with significant accuracy. The parameters reported herein will help predicting RNA duplex stability in different crowded environments, which will lead to an improved understanding of the stability-function relationship for RNAs in various cellular organelles with different molecular environments.

Volumetric Strategy for Quantitatively Elucidating a Local Hydration Network around a G-Quadruplex.

Hydration around nucleic acids, such as DNA and RNA, is an important factor not only for the stability of nucleic acids but also for their interaction with binding molecules. Thus, it is necessary to quantitatively elucidate the hydration properties of nucleic acids around a certain structure. In this study, volumetric changes in G-quadruplex (G4) RNA formation were investigated by systematically changing the number of G-quartet stacks under high pressure. The volumetric contribution at the level of each G4 structural unit revealed that the core G4 helix was significantly more dehydrated than the other parts, including the edges of G-quartets and loops. These findings will help in predicting the binding of G4 ligands on the surface of G4, depending on the chemical structure of the ligand and solution environment. Therefore, the preset volumetric parameter provides information that can predict molecular interactions in G4 formations during molecular crowding in cells.

Intramolecular G-quadruplex-hairpin loop structure competition of a GC-rich exon region in the TMPRSS2 gene.

We identified cytosine-rich regions adjacent to guanine-rich regions in protease genes. A typical GC-rich sequence derived from the TMPRSS2 gene showed structural competition between a G-quadruplex and a hairpin loop, and this competition significantly affected transcription efficiency. These results suggest an impact of neighboring sequences on the gene expression of guanine-rich sequences.

Special Features of COVID-19 in the FMODB: Fragment Molecular Orbital Calculations and Interaction Energy Analysis of SARS-CoV-2-Related Proteins.

SARS-CoV-2 is the causative agent of coronavirus (known as COVID-19), the virus causing the current pandemic. There are ongoing research studies to develop effective therapeutics and vaccines against COVID-19 using various methods and many results have been published. The structure-based drug design of SARS-CoV-2-related proteins is promising, however, reliable information regarding the structural and intra- and intermolecular interactions is required. We have conducted studies based on the fragment molecular orbital (FMO) method for calculating the electronic structures of protein complexes and analyzing their quantitative molecular interactions. This enables us to extensively analyze the molecular interactions in residues or functional group units acting inside the protein complexes. Such precise interaction data are available in the FMO database (FMODB) (https://drugdesign.riken.jp/FMODB/). Since April 2020, we have performed several FMO calculations on the structures of SARS-CoV-2-related proteins registered in the Protein Data Bank. We have published the results of 681 structures, including three structural proteins and 11 nonstructural proteins, on the COVID-19 special page (as of June 8, 2021). In this paper, we describe the entire COVID-19 special page of the FMODB and discuss the calculation results for various proteins. These data not only aid the interpretation of experimentally determined structures but also the understanding of protein functions, which is useful for rational drug design for COVID-19.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence.

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K+ and Na+, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

Ipragliflozin-induced improvement of liver steatosis in obese mice may involve sirtuin signaling.

BACKGROUND: Sodium glucose cotransporter 2 (SGLT2) inhibitors are newly developed oral antidiabetic drugs. SGLT2 is primarily expressed in the kidneys and reabsorbs approximately 90% of the glucose filtered by the renal glomeruli. SGLT2 inhibitors lower glucose levels independently of insulin action by facilitating urinary glucose excretion. The SGLT2 inhibitor ipragliflozin has reportedly improved liver steatosis in animal models and clinical studies. However, the mechanisms by which SGLT2 inhibitors improve liver steatosis are not fully understood. AIM: To investigate the ameliorative effects of ipragliflozin on liver steatosis and the mechanisms of these effects in obese mice. METHODS: We analyzed 8-wk-old male obese (ob/ob) mice that were randomly divided into a group receiving a normal chow diet and a group receiving a normal chow diet supplemented with ipragliflozin (3 mg/kg or 10 mg/kg) for 4 wk. We also analyzed their lean sex-matched littermates receiving a normal chow diet as another control group. Body weight and liver weight were evaluated, and liver histology, immunoblotting, and reverse transcription-polymerase chain reaction analyses were performed. RESULTS: Hepatic lipid accumulation was significantly ameliorated in ob/ob mice treated with 10 mg/kg ipragliflozin compared to untreated ob/ob mice irrespective of body weight changes. Ipragliflozin had no appreciable effects on hepatic oxidative stress-related gene expression levels or macrophage infiltration, but significantly reduced hepatic interleukin-1ß (IL-1ß) mRNA expression levels. Ipragliflozin increased both the mRNA and protein expression levels of sirtuin 1 (SIRT1) in the liver. The hepatic mRNA levels of peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α), peroxisome proliferator-activated receptor α (PPARα), and fibroblast growth factor-21 (FGF21) were also significantly higher in ipragliflozin-treated ob/ob mice than in untreated ob/ob mice. CONCLUSION: Our study suggests that the liver steatosis-ameliorating effects of ipragliflozin in ob/ob mice may be mediated partly by hepatic SIRT1 signaling, possibly through the PGC-1α/PPARα-FGF21 pathway.

Improved nearest-neighbor parameters for the stability of RNA/DNA hybrids under a physiological condition.

The stability of Watson-Crick paired RNA/DNA hybrids is important for designing optimal oligonucleotides for ASO (Antisense Oligonucleotide) and CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas9 techniques. Previous nearest-neighbour (NN) parameters for predicting hybrid stability in a 1 M NaCl solution, however, may not be applicable for predicting stability at salt concentrations closer to physiological condition (e.g. ∼100 mM Na+ or K+ in the presence or absence of Mg2+). Herein, we report measured thermodynamic parameters of 38 RNA/DNA hybrids at 100 mM NaCl and derive new NN parameters to predict duplex stability. Predicted ΔG°37 and Tm values based on the established NN parameters agreed well with the measured values with 2.9% and 1.1°C deviations, respectively. The new results can also be used to make precise predictions for duplexes formed in 100 mM KCl or 100 mM NaCl in the presence of 1 mM Mg2+, which can mimic an intracellular and extracellular salt condition, respectively. Comparisons of the predicted thermodynamic parameters with published data using ASO and CRISPR-Cas9 may allow designing shorter oligonucleotides for these techniques that will diminish the probability of non-specific binding and also improve the efficiency of target gene regulation.

Nearest-neighbor parameters for predicting DNA duplex stability in diverse molecular crowding conditions.

The intracellular environment is crowded and heterogeneous. Although the thermodynamic stability of nucleic acid duplexes is predictable in dilute solutions, methods of predicting such stability under specific intracellular conditions are not yet available. We recently showed that the nearest-neighbor model for self-complementary DNA is valid under molecular crowding condition of 40% polyethylene glycol with an average molecular weight of 200 (PEG 200) in 100 mM NaCl. Here, we determined nearest-neighbor parameters for DNA duplex formation under the same crowding condition to predict the thermodynamics of DNA duplexes in the intracellular environment. Preferential hydration of the nucleotides was found to be the key factor for nearest-neighbor parameters in the crowding condition. The determined parameters were shown to predict the thermodynamic parameters (∆H°, ∆S°, and ∆G°37) and melting temperatures (Tm) of the DNA duplexes in the crowding condition with significant accuracy. Moreover, we proposed a general method for predicting the stability of short DNA duplexes in different cosolutes based on the relationship between duplex stability and the water activity of the cosolute solution. The method described herein would be valuable for investigating biological processes that occur under specific intracellular crowded conditions and for the application of DNA-based biotechnologies in crowded environments.

Effect of Molecular Crowding on the Stability of RNA G-Quadruplexes with Various Numbers of Quartets and Lengths of Loops.

G-Quadruplexes are noncanonical structures formed by guanine-rich regions of not only DNA but also RNA. RNA G-quadruplexes are widely present in the transcriptome as mRNAs and noncoding RNAs and take part in various essential functions in cells. Furthermore, stable RNA G-quadruplexes control the extent of biological functions, such as mRNA translation and antigen presentation. To understand and regulate the functions controlled by RNA G-quadruplexes in cellular environments, which are molecularly crowded, we would be required to investigate the stability of G-quadruplexes in molecular crowding. Here, we systematically investigated the thermodynamic stability of RNA G-quadruplexes with different numbers of G-quartets and lengths of loops. The molecular crowding conditions of polyethylene glycol with an average molecular weight of 200 (PEG200) were found to stabilize RNA G-quadruplexes with three and four G-quartets, while G-quadruplexes with two G-quartets did not exhibit any stabilization upon addition of PEG200. On the other hand, no difference in stabilization by PEG200 was observed among the G-quadruplexes with different loop lengths. Thermodynamic analysis of the RNA G-quadruplexes revealed more appropriate motifs for identifying G-quadruplex-forming sequences. The informatics analysis with new motifs demonstrated that the distributions of G-quadruplexes in human noncoding RNAs differed depending on the number of G-quartets. Therefore, RNA G-quadruplexes with different numbers of G-quartets may play different roles in response to environmental changes in cells.

Hydroxyl groups in cosolutes regulate the G-quadruplex topology of telomeric DNA.

Telomeric G-quadruplex topology has the ability to regulate telomerase activity, which counteracts the shortening of telomere with successive cell divisions, thereby causing genomic longevity. However, the detailed mechanism of G-quadruplexes topologies formed by telomeric sequences requires further investigation. In this study, we quantitatively investigated the effect of cosolutes, particularly the varying number of hydroxyl groups, on the structural transition between hybrid type and parallel G-quadruplexes formed by telomeric DNA sequences. Cosolutes with one or no hydroxyl groups in the vicinal position more efficiently induced the transition to parallel G-quadruplex from hybrid G-quadruplex than those with more hydroxyl groups. We also examined the effect of cosolute structures on the hydration of G-quadruplex formation; the results indicated that cosolutes with fewer hydroxyl groups lead to the release of greater amount of water during G-quadruplex formation. Molecular dynamics results showed that the parallel G-quadruplex was more dehydrated than the hybrid type G-quadruplex. Generally, a dehydrated structure is favored under crowding condition. Thus, depending on the surrounding cosolutes, the G-quadruplex topology can be controlled by the G-quadruplex hydration state.

New Modified Deoxythymine with Dibranched Tetraethylene Glycol Stabilizes G-Quadruplex Structures.

Methods for stabilizing G-quadruplex formation is a promising therapeutic approach for cancer treatment and other biomedical applications because stable G-quadruplexes efficiently inhibit biological reactions. Oligo and polyethylene glycols are promising biocompatible compounds, and we have shown that linear oligoethylene glycols can stabilize G-quadruplexes. Here, we developed a new modified deoxythymine with dibranched or tribranched tetraethylene glycol (TEG) and incorporated these TEG-modified deoxythymines into a loop region that forms an antiparallel G-quadruplex. We analyzed the stability of the modified G-quadruplexes, and the results showed that the tribranched TEG destabilized G-quadruplexes through entropic contributions, likely through steric hindrance. Interestingly, the dibranched TEG modification increased G-quadruplex stability relative to the unmodified DNA structures due to favorable enthalpic contributions. Molecular dynamics calculations suggested that dibranched TEG interacts with the G-quadruplex through hydrogen bonding and CH-π interactions. Moreover, these branched TEG-modified deoxythymine protected the DNA oligonucleotides from degradation by various nucleases in human serum. By taking advantage of the unique interactions between DNA and branched TEG, advanced DNA materials can be developed that affect the regulation of DNA structure.

Effect of Potassium Concentration on Triplex Stability under Molecular Crowding Conditions.

The properties of non-canonical DNA structures, like G-quadruplexes and triplexes, change under cell-mimicking molecular crowding conditions relative to dilute aqueous solutions. The analysis of environmental effects on their stability is crucial since they play important roles in gene expression and regulation. In this study, three intramolecular and intermolecular triplex-forming sequences of different C+*G-C triplet content (*: Hoogsteen base pair; - : Watson-Crick base pair) were designed and their stability measured in the absence and presence of a crowding agent with different K+ concentrations. In dilute solution, the stability of the triplexes was reduced by decreasing the concentration of KCl. This reduction became smaller as the number of C+*G-C triplets increased. Under molecular crowding conditions, Watson-Crick base pairs and Hoogsteen base pairs were destabilized and stabilized, respectively. Interestingly, with lower KCl concentrations (≤1 M), the destabilization of the triplexes due to reduction of KCl concentration was significantly smaller than in dilute solutions. In addition, the C+*G-C content had greater influence on triplex stability under molecular crowding conditions. Our work provides quantitative information about the effects of K+ concentration on triplex stability under molecular crowding conditions and should further our understanding of the function and regulation of triplexes in bioprocesses.

Change in binding states between catabolite activating protein and DNA induced by ligand-binding: molecular dynamics and ab initio fragment molecular orbital calculations.

The transcription mechanism of genetic information from DNA to RNA is efficiently controlled by regulatory proteins, such as catabolite activator protein (CAP), and their ligands. When cyclic AMP (cAMP) binds to CAP, the complex forms a dimer and binds specifically to DNA to activate the transcription mechanism. On the other hand, when cyclic GMP (cGMP) binds to CAP, the complex has no marked effect on the mechanism. In our previous study, based on molecular dynamics (MD) and ab initio fragment molecular orbital (FMO) methods, we elucidated which residues of CAP are important for the specific interactions between CAP and DNA in the CAP-monomer+DNA + cAMP complex. However, this monomer model for CAP cannot describe real interactions between the CAP-dimer and DNA because CAPs form a dimer before binding to DNA. Accordingly, here, we investigated stable structures and their electronic states for the CAP-dimer+DNA complex with cAMP or cGMP ligand, to clarify the influence of ligand-binding on the interactions between CAP-dimer and DNA. The MD simulations elucidated that the DNA-binding domains of CAP-dimer behave differently depending on the ligand bound to the CAP-dimer. In addition, FMO calculations revealed that the binding energy between CAP-dimer and DNA for the CAP-dimer+DNA + cAMP complex is larger than that for the CAP-dimer+DNA + cGMP complex, being consistent with experiments. It was also highlighted that the Arg185 and Lys188 residues of CAP-dimer are important for the binding between CAP-dimer and DNA. These results provide useful information for proposing new compounds that efficiently control the transcription mechanism.

RNA-Capturing Microsphere Particles (R-CAMPs) for Optimization of Functional Aptamers.

RNA aptamers are useful building blocks for constructing functional nucleic acid-based nanoarchitectures. The abilities of aptamers to recognize specific ligands have also been utilized for various biotechnological applications. Solution conditions, which can differ depending on the application, impact the affinity of the aptamers, and thus it is important to optimize the aptamers for the solution conditions to be employed. To simplify the aptamer optimization process, an efficient method that enables re-selection of an aptamer from a partially randomized library is developed. The process relies on RNA-capturing microsphere particles (R-CAMPs): each particle displays different clones of identical DNA and RNA sequences. Using a fluorescence-activated cell sorter, the R-CAMPs that are linked to functional aptamers are sorted. It is demonstrated that after a single round of reselection, several functional aptamers, including the wild-type, are selected from a library of 16 384 sequences. The selection using R-CAMPs is further performed under the solution containing high concentration of ethylene glycol, suggesting applicability in various conditions to optimize an aptamer for a particular application. As any type of RNA clone can be displayed on the microspheres, the technology demonstrated here will be useful for the selection of RNAs based on diverse functions.

Pursuing origins of (poly)ethylene glycol-induced G-quadruplex structural modulations.

Molecular crowding conditions provided by high concentration of cosolutes are utilized for characterization of biomolecules in cell-mimicking environment and development of drug-delivery systems. In this context, (poly)ethylene glycols are often used for studying non-canonical DNA structures termed G-quadruplexes, which came into focus by emerging structural biology findings and new therapeutic drug design approaches. Recently, several reports were made arguing against using (poly)ethylene glycols in role of molecular crowding agents due to their direct impact on DNA G-quadruplex stability and topology. However, the available data on structural details underlying DNA interaction is very scarce and thus limits in-depth comprehension. Herein, structural and thermodynamic analyses were strategically combined to assess G-quadruplex-cosolute interactions and address previously reported variances regarding the driving forces of G-rich DNA structural transformations under molecular crowding conditions. With the use of complementary (CD, NMR and UV) spectroscopic methods and model approach we characterized DNA G-quadruplex in the presence of the smallest and one of the largest typically used (poly)ethylene glycols. Dehydration effect is the key contributor to ethylene-glycol-induced increased stability of the G-quadruplex, which is in the case of the large cosolute mainly guided by the subtle direct interactions between PEG 8000 and the outer G-quartet regions.

Drastic stability change of X-X mismatch in d(CXG) trinucleotide repeat disorders under molecular crowding condition.

The trinucleotide repeat d(CXG) (X = A, C, G or T) is the most common sequence causing repeat expansion disorders. The formation of non-canonical structures, such as hairpin structures with X-X mismatches, has been proposed to affect gene expression and regulation, which are important in pathological studies of these devastating neurological diseases. However, little information is available regarding the thermodynamics of the repeat sequence under crowded cellular conditions where many non-canonical structures such as G-quadruplexes are highly stabilized, while duplexes are destabilised. In this study, we investigated the different stabilities of X-X mismatches in the context of internal d(CXG) self-complementary sequences in an environment with a high concentration of cosolutes to mimic the crowding conditions in cells. The stabilities of full-matched duplexes and duplexes with A-A, G-G, and T-T mismatched base pairs under molecular crowding conditions were notably decreased compared to under dilute conditions. However, the stability of the DNA duplex with a C-C mismatch base pair was only slightly destabilised. Investigating different stabilities of X-X mismatches in d(CXG) sequences is important for improving our understanding of the formation and transition of multiple non-canonical structures in trinucleotide repeat diseases, and may provide insights for pathological studies and drug development.