From RNA sequence to its three-dimensional structure: geometrical structure, stability and dynamics of selected fragments of SARS-CoV-2 RNA.

In this computational study, we explore the folding of a particular sequence using various computational tools to produce two-dimensional structures, which are then transformed into three-dimensional structures. We then study the geometry, energetics and dynamics of these structures using full electron quantum-chemical and classical molecular dynamics calculations. Our study focuses on the SARS-CoV-2 RNA fragment GGaGGaGGuguugcaGG and its various structures, including a G-quadruplex and five different hairpins. We examine the impact of two types of counterions (K+ and Na+) and flanking nucleotides on their geometrical characteristics, relative stability and dynamic properties. Our results show that the G-quadruplex structure is the most stable among the constructed hairpins. We confirm its topological stability through molecular dynamics simulations. Furthermore, we observe that the nucleotide loop consisting of seven nucleotides is the most flexible part of the RNA fragment. Additionally, we find that RNA networks of intermolecular hydrogen bonds are highly sensitive to the surrounding environment. Our findings reveal the loss of 79 old hydrogen bonds and the formation of 91 new ones in the case when the G-quadruplex containing flanking nucleotides is additionally stabilized by Na+ counterions.

4-(Azolyl)-Benzamidines as a Novel Chemotype for ASIC1a Inhibitors.

Acid-sensing ion channels (ASICs) play a key role in the perception and response to extracellular acidification changes. These proton-gated cation channels are critical for neuronal functions, like learning and memory, fear, mechanosensation and internal adjustments like synaptic plasticity. Moreover, they play a key role in neuronal degeneration, ischemic neuronal injury, seizure termination, pain-sensing, etc. Functional ASICs are homo or heterotrimers formed with (ASIC1-ASIC3) homologous subunits. ASIC1a, a major ASIC isoform in the central nervous system (CNS), possesses an acidic pocket in the extracellular region, which is a key regulator of channel gating. Growing data suggest that ASIC1a channels are a potential therapeutic target for treating a variety of neurological disorders, including stroke, epilepsy and pain. Many studies were aimed at identifying allosteric modulators of ASIC channels. However, the regulation of ASICs remains poorly understood. Using all available crystal structures, which correspond to different functional states of ASIC1, and a molecular dynamics simulation (MD) protocol, we analyzed the process of channel inactivation. Then we applied a molecular docking procedure to predict the protein conformation suitable for the amiloride binding. To confirm the effect of its sole active blocker against the ASIC1 state transition route we studied the complex with another MD simulation run. Further experiments evaluated various compounds in the Enamine library that emerge with a detectable ASIC inhibitory activity. We performed a detailed analysis of the structural basis of ASIC1a inhibition by amiloride, using a combination of in silico approaches to visualize its interaction with the ion pore in the open state. An artificial activation (otherwise, expansion of the central pore) causes a complex modification of the channel structure, namely its transmembrane domain. The output protein conformations were used as a set of docking models, suitable for a high-throughput virtual screening of the Enamine chemical library. The outcome of the virtual screening was confirmed by electrophysiological assays with the best results shown for three hit compounds.

Integrated workflow for the identification of new GABAA R positive allosteric modulators based on the in†silico screening with further in†vitro validation. Case study using Enamine's stock chemical space.

Numerous studies reported an association between GABAA R subunit genes and epilepsy, eating disorders, autism spectrum disorders, neurodevelopmental disorders, and bipolar disorders. This study was aimed to find some potential positive allosteric modulators and was performed by combining the in silico approach with further in vitro evaluation of its real activity. We started from the GABAA R-diazepam complexes and assembled a lipid embedded protein ensemble to refine it via molecular dynamics (MD) simulation. Then we focused on the interaction of α1ß2γ2 with some Z-drugs (non-benzodiazepine compounds) using an Induced Fit Docking (IFD) into the relaxed binding site to generate a pharmacophore model. The pharmacophore model was validated with a reference set and applied to decrease the pre-filtered Enamine database before the main docking procedure. Finally, we succeeded in identifying a set of compounds, which met all features of the docking model. The aqueous solubility and stability of these compounds in mouse plasma were assessed. Then they were tested for the biological activity using the rat Purkinje neurons and CHO cells with heterologously expressed human α1ß2γ2 GABAA receptors. Whole-cell patch clamp recordings were used to reveal the GABA induced currents. Our study represents a convenient and tunable model for the discovery of novel positive allosteric modulators of GABAA receptors. A High-throughput virtual screening of the largest available database of chemical compounds resulted in the selection of 23 compounds. Further electrophysiological tests allowed us to determine a set of 3 the most outstanding active compounds. Considering the structural features of leader compounds, the study can develop into the MedChem project soon.