Structure and Dynamics of RNA Guanine Quadruplexes in SARS-CoV-2 Genome. Original Strategies against Emerging Viruses

Guanine quadruplex (G4) structures in the viral genome have a key role in modulating viruses' biological activity. While several DNA G4 structures have been experimentally resolved, RNA G4s are definitely less explored. We report the first calculated G4 structure of the RG-1 RNA sequence of SARS-CoV-2 genome, obtained by using a multiscale approach combining quantum and classical molecular modeling and corroborated by the excellent agreement between the corresponding calculated and experimental circular dichroism spectra. We prove the stability of the RG-1 G4 arrangement as well as its interaction with G4 ligands potentially inhibiting viral protein translation

Resolving a guanine-quadruplex structure in the SARS-CoV-2 genome through circular dichroism and multiscale molecular modeling

The genome of SARS-CoV-2 coronavirus is made up of a single-stranded RNA fragment that can assume a specific secondary structure, whose stability can influence the virus's ability to reproduce. Recent studies have identified putative guanine quadruplex sequences in SARS-CoV-2 genome fragments that are involved in coding for both structural and non-structural proteins. In this contribution, we focus on a specific G-rich sequence referred to as RG-2, which codes for the non-structural protein 10 (Nsp10) and assumes a guanine-quadruplex (G4) arrangement. We provide the secondary structure of RG-2 G4 at atomistic resolution by molecular modeling and simulation, validated by the superposition of experimental and calculated electronic circular dichroism spectra. Through both experimental and simulation approaches, we have demonstrated that pyridostatin (PDS), a widely recognized G4 binder, can bind to and stabilize RG-2 G4 more strongly than RG-1, another G4 forming sequence that was previously proposed as a potential target for antiviral drug candidates. Overall, this study highlights RG-2 as a valuable target to inhibit the translation and replication of SARS-CoV-2, paving the way towards original therapeutic approaches against emerging RNA viruses.Parallel or hybrid? A combination of multiscale molecular modeling and circular dichroism is used to predict a G-quadruplex structure at atomistic resolution in the SARS-CoV-2 genome

Forever Young: Structural Stability of Telomeric Guanine Quadruplexes in the Presence of Oxidative DNA Lesions**

Human telomeric DNA, in G-quadruplex (G4) conformation, is characterized by a remarkable structural stability that confers it the capacity to resist to oxidative stress producing one or even clustered 8-oxoguanine (8oxoG) lesions. We present a combined experimental/computational investigation, by using circular dichroism in aqueous solutions, cellular immunofluorescence assays and molecular dynamics simulations, that identifies the crucial role of the stability of G4s to oxidative lesions, related also to their biological role as inhibitors of telomerase, an enzyme overexpressed in most cancers associated to oxidative stress

Molecular Basis of SARS-CoV-2 Infection and Rational Design of Potential Antiviral Agents: Modeling and Simulation Approaches

The emergence in late 2019 of the coronavirus SARS-CoV-2 has resulted in the breakthrough of the COVID-19 pandemic that is presently affecting a growing number of countries. The development of the pandemic has also prompted an unprecedented effort of the scientific community to understand the molecular bases of the virus infection and to propose rational drug design strategies able to alleviate the serious COVID-19 morbidity. In this context, a strong synergy between the structural biophysics and molecular modeling and simulation communities has emerged, resolving at the atomistic level the crucial protein apparatus of the virus and revealing the dynamic aspects of key viral processes. In this Review, we focus on how in silico studies have contributed to the understanding of the SARS-CoV-2 infection mechanism and the proposal of novel and original agents to inhibit the viral key functioning. This Review deals with the SARS-CoV-2 spike protein, including the mode of action that this structural protein uses to entry human cells, as well as with nonstructural viral proteins, focusing the attention on the most studied proteases and also proposing alternative mechanisms involving some of its domains, such as the SARS unique domain. We demonstrate that molecular modeling and simulation represent an effective approach to gather information on key biological processes and thus guide rational molecular design strategies

Genetic analysis of the flowering date and number of petals in rose

Rose is the ornamental species with the highest financial impact. Floral traits such as the number of petals and the date of flowering are major characteristics of ornamental plants. Our objective is to study the genetic determinism of floral traits: date of flowering and number of petals, which are a major issue for rose breeders. The study was conducted on two interspecific populations interconnected by the male parent: H190 x hybrid of Rosa wichurana (referred to as HW) and “The Fairy” x hybrid of Rosa wichurana (referred to as FW). The number of petals and the date of flowering were scored over 2 and 8 years, respectively. A new HW genetic map covering 468 cM and the already available genetic map of the FW population (Kawamura et al. TAG Theor Appl Genet 122:661–675, 2011) were used for the genetic determinism studies. In each population, half of the hybrids exhibited single flowers (less than 10 petals), whereas the other half revealed double flowers. The number of petals is controlled by the NP gene located on LG3. Additionally, we detected two new major quantitative trait loci (QTLs) on LG2 and LG5, close to RoAP1b and RoRAG, respectively, two genes involved in the control of floral identity. For the date of flowering, three QTLs with a major effect and high stability between years were found on linkage groups 3, 4, and 6, indicating a high stability of QTLs to the changing environment. Candidate genes underlying these QTLs were investigated and key genes were identified. These major QTLs were linked to candidate genes, i.e., the identified QTL on LG4 was linked to RoFT, the one on LG3 to genes involved in gibberellin pathways, and the one on LG6 to RoFD. These QTLs, which are very stable over time, are good candidates to develop markers applicable in marker-assisted selection (MAS)

Structure of the 5′ untranslated region in SARS-CoV-2 genome and its specific recognition by innate immune system via the human oligoadenylate synthase 1

2′-5′-Oligoadenylate synthetase 1 (OAS1) is one of the key enzymes driving the innate immune system response to SARS-CoV-2 infection whose activity has been related to COVID-19 severity. OAS1 is a sensor of endogenous RNA that triggers the 2′-5′-oligoadenylate/RNase L pathway. Upon SARS-CoV-2 infection, OAS1 is responsible for the recognition of viral RNA and has been shown to possess a particularly high sensitivity for the 5′-untranslated (5′-UTR) RNA region, which is organized in a double-strand stem loop motif (SL1). Here we report the structure of the SL1/OAS1 complex also rationalizing the high affinity for OAS1

Specific Recognition of the 5′-Untranslated Region of West Nile Virus Genome by Human Innate Immune System

In the last few years, the sudden outbreak of COVID-19 caused by SARS-CoV-2 proved the crucial importance of understanding how emerging viruses work and proliferate, in order to avoid the repetition of such a dramatic sanitary situation with unprecedented social and economic costs. West Nile Virus is a mosquito-borne pathogen that can spread to humans and induce severe neurological problems. This RNA virus caused recent remarkable outbreaks, notably in Europe, highlighting the need to investigate the molecular mechanisms of its infection process in order to design and propose efficient antivirals. Here, we resort to all-atom Molecular Dynamics simulations to characterize the structure of the 5′-untranslated region of the West Nile Virus genome and its specific recognition by the human innate immune system via oligoadenylate synthetase. Our simulations allowed us to map the interaction network between the viral RNA and the host protein, which drives its specific recognition and triggers the host immune response. These results may provide fundamental knowledge that can assist further antivirals’ design, including therapeutic RNA strategies

Role of RNA Guanine Quadruplexes in Favoring the Dimerization of SARS Unique Domain in Coronaviruses

Coronaviruses may produce severe acute respiratory syndrome (SARS). As a matter of fact, a new SARS-type virus, SARS-CoV-2, is responsible for the global pandemic in 2020 with unprecedented sanitary and economic consequences for most countries. In the present contribution we study, by all-atom equilibrium and enhanced sampling molecular dynamics simulations, the interaction between the SARS Unique Domain and RNA guanine quadruplexes, a process involved in eluding the defensive response of the host thus favoring viral infection of human cells. Our results evidence two stable binding modes involving an interaction site spanning either the protein dimer interface or only one monomer. The free energy profile unequivocally points to the dimer mode as the thermodynamically favored one. The effect of these binding modes in stabilizing the protein dimer was also assessed, being related to its biological role in assisting the SARS viruses to bypass the host protective response. This work also constitutes a first step in the possible rational design of efficient therapeutic agents aiming at perturbing the interaction between SARS Unique Domain and guanine quadruplexes, hence enhancing the host defenses against the virus

How Fragile We Are: Influence of Stimulator of Interferon Genes (STING) Variants on Pathogen Recognition and Immune Response Efficiency

The stimulator of interferon genes (STING) protein is a cornerstone of the human immune response. Its activation by cGAMP in the presence of cytosolic DNA stimulates the production of type I interferons and inflammatory cytokines. In the human population, several STING variants exist and exhibit dramatic differences in their activity, impacting the efficiency of the host defense against infections. Understanding the molecular mechanisms of these variants opens perspectives for personalized medicine treatments against diseases such as viral infections, cancers, or autoinflammatory diseases. Through microsecond-scale molecular modeling simulations, contact analyses, and machine learning techniques, we reveal the dynamic behavior of four STING variants (wild type, G230A, R293Q, and G230A/R293Q) and rationalize the variability of efficiency observed experimentally. Our results show that the decrease in STING activity is linked to a stiffening of key structural elements of the binding cavity together with changes in the interaction patterns within the protein

G-Quadruplex Recognition by DARPIns through Epitope/Paratope Analogy**

We investigated the mechanisms leading to the specific recognition of Guanine Guadruplex (G4) by DARPins peptides, which can lead to the design of G4s specific sensors. To this end we carried out all-atom molecular dynamics simulations to unravel the interactions between specific nucleic acids, including human-telomeric (h-telo), Bcl-2, and c-Myc, with different peptides, forming a DARPin/G4 complex. By comparing the sequences of DARPin with that of a peptide known for its high affinity for c-Myc, we show that the recognition cannot be ascribed to sequence similarity but, instead, depends on the complementarity between the three-dimensional arrangement of the molecular fragments involved: the α-helix/loops domain of DARPin and the G4 backbone. Our results reveal that DARPins tertiary structure presents a charged hollow region in which G4 can be hosted, thus the more complementary the structural shapes, the more stable the interaction