Molecular dynamics recipes for genome research

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain–domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains

Kv Channel Gating Requires a Compatible S4-S5 Linker and Bottom Part of S6, Constrained by Non-interacting Residues

Voltage-dependent K+ channels transfer the voltage sensor movement into gate opening or closure through an electromechanical coupling. To test functionally whether an interaction between the S4-S5 linker (L45) and the cytoplasmic end of S6 (S6T) constitutes this coupling, the L45 in hKv1.5 was replaced by corresponding hKv2.1 sequence. This exchange was not tolerated but could be rescued by also swapping S6T. Exchanging both L45 and S6T transferred hKv2.1 kinetics to an hKv1.5 background while preserving the voltage dependence. A one-by-one residue substitution scan of L45 and S6T in hKv1.5 further shows that S6T needs to be α-helical and forms a “crevice” in which residues I422 and T426 of L45 reside. These residues transfer the mechanical energy onto the S6T crevice, whereas other residues in S6T and L45 that are not involved in the interaction maintain the correct structure of the coupling

Novel phenotype associated with a mutation in the KCNA1(Kv1.1) gene

Episodic ataxia type 1 (EA1) is an autosomal dominant K+ channelopathy which manifests with short attacks of cerebellar ataxia and dysarthria, and may also show interictal myokymia. Episodes can be triggered by emotional or physical stress, startle response, sudden postural change or fever. Here we describe a 31-year-old man displaying markedly atypical symptoms, including long-lasting attacks of jerking muscle contractions associated with hyperthermia, severe migraine, and a relatively short-sleep phenotype. A single nucleotide change in KCNA1 (c.555C>G) was identified that changes a highly conserved residue (p.C185W) in the first transmembrane segment of the voltage-gated K+ channel Kv1.1. The patient is heterozygous and the mutation was inherited from his asymptomatic mother. Next generation sequencing revealed no variations in the CACNA1A, CACNB4, KCNC3, KCNJ10, PRRT2 or SCN8A genes of either the patient or mother, except for a benign variant in SLC1A3. Functional analysis of the p.C185W mutation in KCNA1 demonstrated a deleterious dominant-negative phenotype where the remaining current displayed slower activation kinetics, subtle changes in voltage-dependence and faster recovery from slow inactivation. Structural modeling also predicts the C185W mutation to be functionally deleterious. This description of novel clinical features, associated with a Kv1.1 mutation highlights a possibly unrecognized relationship between K+ channel dysfunction, hyperthermia and migraine in EA1, and suggests that thorough assessments for these symptoms should be carefully considered for all patients affected by EA1.peer-reviewe

Gain-of-function defects of astrocytic Kir4.1 channels in children with autism spectrum disorders and epilepsy

Dysfunction of the inwardly-rectifying potassium channels Kir4.1 (KCNJ10) represents a pathogenic mechanism contributing to Autism-Epilepsy comorbidity. To define the role of Kir4.1 variants in the disorder, we sequenced KCNJ10 in a sample of affected individuals, and performed genotype-phenotype correlations. The effects of mutations on channel activity, protein trafficking, and astrocyte function were investigated in Xenopus laevis oocytes, and in human astrocytoma cell lines. An in vivo model of the disorder was also explored through generation of kcnj10a morphant zebrafish overexpressing the mutated human KCNJ10. We detected germline heterozygous KCNJ10 variants in 19/175 affected children. Epileptic spasms with dysregulated sensory processing represented the main disease phenotype. When investigated on astrocyte-like cells, the p.R18Q mutation exerted a gain-of-function effect by enhancing Kir4.1 membrane expression and current density. Similarly, the p.R348H variant led to gain of channel function through hindrance of pH-dependent current inhibition. The frequent polymorphism p.R271C seemed, instead, to have no obvious functional effects. Our results confirm that variants in KCNJ10 deserve attention in autism-epilepsy, and provide insight into the molecular mechanisms of autism and seizures. Similar to neurons, astrocyte dysfunction may result in abnormal synaptic transmission and electrical discharge, and should be regarded as a possible pharmacological target in autism-epilepsy. Supplementary information accompanies this paper in the files section.peer-reviewe

A Calsequestrin-1 mutation associated with a skeletal muscle disease alters sarcoplasmic Ca2+ release

An autosomal dominant protein aggregate myopathy, characterized by high plasma creatine kinase and calsequestrin-1 (CASQ1) accumulation in skeletal muscle, has been recently associated with a missense mutation in CASQ1 gene. The mutation replaces an evolutionarily-conserved aspartic acid with glycine at position 244 (p.D244G) of CASQ1, the main sarcoplasmic reticulum (SR) Ca2+ binding and storage protein localized at the terminal cisternae of skeletal muscle cells. Here, immunocytochemical analysis of myotubes, differentiated from muscle-derived primary myoblasts, shows that sarcoplasmic vacuolar aggregations positive for CASQ1 are significantly larger in CASQ1-mutated cells than control cells. A strong co-immuno staining of both RyR1 and CASQ1 was also noted in the vacuoles of myotubes and muscle biopsies derived from patients. Electrophysiological recordings and sarcoplasmic Ca2+ measurements provide evidence for less Ca2+ release from the SR of mutated myotubes when compared to that of controls. These findings further clarify the pathogenic nature of the p.D244G variant and point out defects in sarcoplasmic Ca2+ homeostasis as a mechanism underlying this human disease, which could be distinctly classified as "CASQ1-couplonopathy".peer-reviewe

Модельные характеристики кардиореспираторной системы высококвалифицированных гребцов на байдарках и каноэ

The most frequently used approach for protein structure prediction is currently homology modeling. The 3D model building phase of this methodology is critical for obtaining an accurate and biologically useful prediction. The most widely employed tool to perform this task is MODELLER. This program implements the "modeling by satisfaction of spatial restraints" strategy and its core algorithm has not been altered significantly since the early 1990s. In this work, we have explored the idea of modifying MODELLER with two effective, yet computationally light strategies to improve its 3D modeling performance. Firstly, we have investigated how the level of accuracy in the estimation of structural variability between a target protein and its templates in the form of σ values profoundly influences 3D modeling. We show that the σ values produced by MODELLER are on average weakly correlated to the true level of structural divergence between target-template pairs and that increasing this correlation greatly improves the program's predictions, especially in multiple-template modeling. Secondly, we have inquired into how the incorporation of statistical potential terms (such as the DOPE potential) in the MODELLER's objective function impacts positively 3D modeling quality by providing a small but consistent improvement in metrics such as GDT-HA and lDDT and a large increase in stereochemical quality. Python modules to harness this second strategy are freely available at https://github.com/pymodproject/altmod. In summary, we show that there is a large room for improving MODELLER in terms of 3D modeling quality and we propose strategies that could be pursued in order to further increase its performance

La biologia computazionale come strumento di indagine delle malattie genetiche

I canali ionici sono proteine di membrana che regolano il passaggio di ioni specifici verso l’interno e l’esterno delle cellule. Tra questi, i canali del potassio rivestono un’importanza fondamentale nel contesto di molte patologie di natura genetica: le malattie metaboliche sono spesso causate, infatti, da un malfunzionamento di specifici canali ionici, la cui funzione biologica risulta alterata in seguito ad una mutazione del relativo gene che la codifica. La biologia computazionale consente di indagare e studiare, a livello molecolare, le cause che alterano la funzione dei canali ionici; fornisce a medici e biologi uno strumento di interpretazione dei loro esperimenti e permette la progettazione di nuovi farmaci volti a curare e migliorare la qualità della vita di molti pazienti. Le moderne risorse di calcolo parallelo, basate sui nuovi processori GPU ad architettura multicore, consentono di simulare il comportamento dinamico dei canali ionici su scale di tempo molto maggiori

STRUMENTI COMPUTAZIONALI PER LA MEDICINA E LA BIOLOGIA: DINAMICA MOLECOLARE DI PROTEINE E DNA

Le simulazioni di dinamica molecolare consentono di analizzare la struttura e la flessibilità di macromolecole importanti dal punto di vista biomedico, come proteine e DNA. Questo permette di studiare le cause di molte malattie genetiche, spesso provocate da una mutazione di un gene che codifica per una proteina, ma anche di progettare nuovi farmaci più efficaci. In particolare, lo studio delle mutazioni che insorgono su proteine di membrana come i canali ionici voltaggio-dipendenti permettono di delucidare i meccanismi molecolari alla base delle più importanti patologie neurologiche e muscolari. Le moderne risorse di calcolo parallelo, basate sui nuovi processori GPU ad architettura multi-core, consentono di estendere la simulazione di questi sistemi biologici su scale di tempo molto maggiori e quindi di modellare funzioni biologiche più complesse. Ma lo sviluppo non è solo nelle tecnologie. Nuovi algoritmi e campi di forza che regolano le interazioni inter-atomiche vengono continuamente proposti e implementati nell’ambiente di calcolo del CASPUR