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.

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.

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.

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.

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.

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

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.

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.

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.

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.

Newly characterized interaction stabilizes DNA structure: oligoethylene glycols stabilize G-quadruplexes CH-π interactions.

Oligoethylene glycols are used as crowding agents in experiments that aim to understand the effects of intracellular environments on DNAs. Moreover, DNAs with covalently attached oligoethylene glycols are used as cargo carriers for drug delivery systems. To investigate how oligoethylene glycols interact with DNAs, we incorporated deoxythymidine modified with oligoethylene glycols of different lengths, such as tetraethylene glycol (TEG), into DNAs that form antiparallel G-quadruplex or hairpin structures such that the modified residues were incorporated into loop regions. Thermodynamic analysis showed that because of enthalpic differences, the modified G-quadruplexes were stable and the hairpin structures were slightly unstable relative to unmodified DNA. The stability of G-quadruplexes increased with increasing length of the ethylene oxides and the number of deoxythymidines modified with ethylene glycols in the G-quadruplex. Nuclear magnetic resonance analyses and molecular dynamics calculations suggest that TEG interacts with bases in the G-quartet and loop via CH-π and lone pair-π interactions, although it was previously assumed that oligoethylene glycols do not directly interact with DNAs. The results suggest that numerous cellular co-solutes likely affect DNA function through these CH-π and lone pair-π interactions.

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.

Elevated serum uric acid level was a notable adverse event during combination therapy with sofosbuvir and ribavirin.

AIM: Combination therapy with sofosbuvir and ribavirin (SOF/RBV) has been recently available for chronic hepatitis C patients with genotype 2 (CHG2) in Japan. The domestic phase III clinical trial showed a high antiviral effect with a relatively safe adverse event (AE) profile. Our aim was to report an important AE detected during treatment. METHODS: A prospective multi-institutional study of 12-week combination therapy with SOF/RBV for CHG2 was carried out to evaluate efficacy and safety. RESULTS: The eligible subjects included 142 patients. Out of 50 assessable patients, 16% of the patients were diagnosed with hyperuricemia. The proportions of subjects with grade 1, grade 3, and grade 4 hyperuricemia were 12, 2, and 2%, respectively. Serum uric acid (UA) levels at week 1 of the therapy (W1) were numerically the highest during therapy in patients with hyperuricemia, and the ratio of W1/baseline serum UA levels was significantly higher than that of post-treatment week 4 or 8/baseline serum UA levels in assessable patients. Serum UA levels at W1 were significantly correlated with body mass index. The difference between serum UA levels at W1 and baseline serum UA levels was significantly correlated with the difference between serum creatinine levels at W1 and baseline serum creatinine levels. CONCLUSIONS: Elevated serum UA level was a notable AE associated with SOF/RBV therapy for CHG2. However, because of the small number of subjects, the exact frequency of AEs should be re-evaluated with larger cohorts. We need to remember that elevated serum UA level might develop during the therapy, especially at W1.

Clinical and virological features of acute hepatitis E in Gunma prefecture, Japan between 2004 and 2015.

AIM: To evaluate the clinical and virological features of acute hepatitis E (AH-E) in Gunma prefecture and focus on the hepatitis E virus (HEV) infection in immunocompromised patients. METHODS: A total of 30 patients with AH-E diagnosed at our Gunma University Hospital, and located in 3-39-15 Showa-machi, Maebashi, Gunma 371-8511 Japan, and its affiliated hospitals from 2004 to 2015, were studied. We evaluated the detailed medical histories, laboratory examinations and virological features of these participants. RESULTS: Of the 30 patients, 21 patients were men, with a median age of 61 years. Three of these patients had a history of recent oversea travel. A total of 14 patients had eaten raw or undercooked meat/viscera from animals, and two patients had contracted transfusion-transmitted AH-E. Eight patients were immunocompromised, including those with hematological disease, cancer receiving systemic chemotherapy and kidney transplant or connective tissue disease undergoing immunosuppressive medications. The alanine aminotransferase and total bilirubin levels were more significantly reduced in these immunocompromised patients than in the non-immunocompromised patients. Severe thrombocytopenia, an extra-hepatic manifestation of AH-E, occurred in one case. Among the 22 HEV strains whose subgenotype was determined, two were imported strains (1a and 1f), and 11 strains formed four distinct phylogenetic clusters within subgenotype 3b. The remaining nine strains differed from each other by 9.8-22.4%, and were classified into four subgenotypes (3a, 3b, 3e and 3f). CONCLUSION: Markedly divergent HEV strains (3a, 3b, 3e and 3f) were found to circulate in Gunma. Although immunosuppression appears to play a crucial role in establishing chronic sequels, AH-E in eight immunocompromised patients, including transfusion-transmitted HEV infection in two patients, did not become chronic.

Transforming growth factor-α activates pancreatic stellate cells and may be involved in matrix metalloproteinase-1 upregulation.

The role that transforming growth factor-α (TGF-α) has in chronic pancreatitis and pancreatic cancer has not been fully elucidated. We evaluated the effects of TGF-α on the human pancreatic stellate cell (PSC) line RLT-PSC and primary human PSCs, and the expression levels of TGF-α and metalloproteinase-1 (MMP-1) in human chronic pancreatitis and pancreatic cancer tissues. TGF-α stimulated the proliferation and migration of PSCs. Although the mRNA expression levels of tissue inhibitor of metalloproteinase-1 and α1(I) collagen were unchanged, the mRNA expression levels of MMP-1 increased concomitant with increases in MMP-1 protein levels and collagenase activity. TGF-α-stimulated migration of RLT-PSC cells was partially blocked by tissue inhibitor of metalloproteinase-1 protein and MMP-1 small interfering RNA. MMP-1 was also observed to stimulate the migration of PSCs. TGF-α-induced MMP-1 expression was completely blocked by gefitinib in PSCs. The Ras-ERK and PI3/Akt pathways appear to be involved in the activation of MMP-1 in PSCs. Immunohistochemical analyses showed that MMP-1 expression was significantly increased in the pancreatic interstitial tissues in case of chronic pancreatitis or pancreatic cancer compared with those in case of normal pancreas. In conclusion, TGF-α increased proliferation and migration of PSCs. TGF-α-induced migration of cells may be partly due to upregulation of MMP-1. TGF-α and MMP-1 upregulation may contribute to the pathogenesis of chronic pancreatitis and pancreatic cancer.

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.