Extreme phenotypic heterogeneity in non-expansion spinocerebellar ataxias.

Although the best-known spinocerebellar ataxias (SCAs) are triplet repeat diseases, many SCAs are not caused by repeat expansions. The rarity of individual non-expansion SCAs, however, has made it difficult to discern genotype-phenotype correlations. We therefore screened individuals who had been found to bear variants in a non-expansion SCA-associated gene through genetic testing, and after we eliminated genetic groups that had fewer than 30 subjects, there were 756 subjects bearing single-nucleotide variants or deletions in one of seven genes: CACNA1A (239 subjects), PRKCG (175), AFG3L2 (101), ITPR1 (91), STUB1 (77), SPTBN2 (39), or KCNC3 (34). We compared age at onset, disease features, and progression by gene and variant. There were no features that reliably distinguished one of these SCAs from another, and several genes-CACNA1A, ITPR1, SPTBN2, and KCNC3-were associated with both adult-onset and infantile-onset forms of disease, which also differed in presentation. Nevertheless, progression was overall very slow, and STUB1-associated disease was the fastest. Several variants in CACNA1A showed particularly wide ranges in age at onset: one variant produced anything from infantile developmental delay to ataxia onset at 64 years of age within the same family. For CACNA1A, ITPR1, and SPTBN2, the type of variant and charge change on the protein greatly affected the phenotype, defying pathogenicity prediction algorithms. Even with next-generation sequencing, accurate diagnosis requires dialogue between the clinician and the geneticist.

A novel in-frame deletion in MYOT causes an early adult onset distal myopathy.

Missense mutations in MYOT encoding the sarcomeric Z-disk protein myotilin cause three main myopathic phenotypes including proximal limb-girdle muscular dystrophy, spheroid body myopathy, and late-onset distal myopathy. We describe a family carrying a heterozygous MYOT deletion (Tyr4_His9del) that clinically was characterized by an early-adult onset distal muscle weakness and pathologically by a myofibrillar myopathy (MFM). Molecular modeling of the full-length myotilin protein revealed that the 4-YERPKH-9 amino acids are involved in local interactions within the N-terminal portion of myotilin. Injection of in vitro synthetized mutated human MYOT RNA or of plasmid carrying its cDNA sequence in zebrafish embryos led to muscle defects characterized by sarcomeric disorganization of muscle fibers and widening of the I-band, and severe motor impairments. We identify MYOT novel Tyr4_His9 deletion as the cause of an early-onset MFM with a distal myopathy phenotype and provide data supporting the importance of the amino acid sequence for the structural role of myotilin in the sarcomeric organization of myofibers.

Structure-function relationships of the disease-linked A218T oxytocin receptor variant.

Various single nucleotide polymorphisms (SNPs) in the oxytocin receptor (OXTR) gene have been associated with behavioral traits, autism spectrum disorder (ASD) and other diseases. The non-synonymous SNP rs4686302 results in the OXTR variant A218T and has been linked to core characteristics of ASD, trait empathy and preterm birth. However, the molecular and intracellular mechanisms underlying those associations are still elusive. Here, we uncovered the molecular and intracellular consequences of this mutation that may affect the psychological or behavioral outcome of oxytocin (OXT)-treatment regimens in clinical studies, and provide a mechanistic explanation for an altered receptor function. We created two monoclonal HEK293 cell lines, stably expressing either the wild-type or A218T OXTR. We detected an increased OXTR protein stability, accompanied by a shift in Ca2+ dynamics and reduced MAPK pathway activation in the A218T cells. Combined whole-genome and RNA sequencing analyses in OXT-treated cells revealed 7823 differentially regulated genes in A218T compared to wild-type cells, including 429 genes being associated with ASD. Furthermore, computational modeling provided a molecular basis for the observed change in OXTR stability suggesting that the OXTR mutation affects downstream events by altering receptor activation and signaling, in agreement with our in vitro results. In summary, our study provides the cellular mechanism that links the OXTR rs4686302 SNP with genetic dysregulations associated with aspects of ASD.

Low Molecular Weight Inhibitors Targeting the RNA-Binding Protein HuR.

The RNA-binding protein human antigen R (HuR) regulates stability, translation, and nucleus-to-cytoplasm shuttling of its target mRNAs. This protein has been progressively recognized as a relevant therapeutic target for several pathologies, like cancer, neurodegeneration, as well as inflammation. Inhibitors of mRNA binding to HuR might thus be beneficial against a variety of diseases. Here, we present the rational identification of structurally novel HuR inhibitors. In particular, by combining chemoinformatic approaches, high-throughput virtual screening, and RNA-protein pulldown assays, we demonstrate that the 4-(2-(2,4,6-trioxotetrahydropyrimidin-5(2H)-ylidene)hydrazineyl)benzoate ligand exhibits a dose-dependent HuR inhibition effect in binding experiments. Importantly, the chemical scaffold is new with respect to the currently known HuR inhibitors, opening up a new avenue for the design of pharmaceutical agents targeting this important protein.

Robust principal component analysis-based prediction of protein-protein interaction hot spots.

Proteins often exert their function by binding to other cellular partners. The hot spots are key residues for protein-protein binding. Their identification may shed light on the impact of disease associated mutations on protein complexes and help design protein-protein interaction inhibitors for therapy. Unfortunately, current machine learning methods to predict hot spots, suffer from limitations caused by gross errors in the data matrices. Here, we present a novel data pre-processing pipeline that overcomes this problem by recovering a low rank matrix with reduced noise using Robust Principal Component Analysis. Application to existing databases shows the predictive power of the method.

MERMAID: dedicated web server to prepare and run coarse-grained membrane protein dynamics.

Atomistic molecular dynamics simulations of membrane proteins have been shown to be extremely useful for characterizing the molecular features underlying their function, but require high computational power, limiting the understanding of complex events in membrane proteins, e.g. ion channels gating, GPCRs activation. To overcome this issue, it has been shown that coarse-grained approaches, although requiring less computational power, are still capable of correctly describing molecular events underlying big conformational changes in biological systems. Here, we present the Martini coarse-grained membrane protein dynamics (MERMAID), a publicly available web interface that allows the user to prepare and run coarse-grained molecular dynamics (CGMD) simulations and to analyse the trajectories.

A Mechanistic Model of NMDA and AMPA Receptor-Mediated Synaptic Transmission in Individual Hippocampal CA3-CA1 Synapses: A Computational Multiscale Approach.

Inside hippocampal circuits, neuroplasticity events that individual cells may undergo during synaptic transmissions occur in the form of Long-Term Potentiation (LTP) and Long-Term Depression (LTD). The high density of NMDA receptors expressed on the surface of the dendritic CA1 spines confers to hippocampal CA3-CA1 synapses the ability to easily undergo NMDA-mediated LTP and LTD, which is essential for some forms of explicit learning in mammals. Providing a comprehensive kinetic model that can be used for running computer simulations of the synaptic transmission process is currently a major challenge. Here, we propose a compartmentalized kinetic model for CA3-CA1 synaptic transmission. Our major goal was to tune our model in order to predict the functional impact caused by disease associated variants of NMDA receptors related to severe cognitive impairment. Indeed, for variants Glu413Gly and Cys461Phe, our model predicts negative shifts in the glutamate affinity and changes in the kinetic behavior, consistent with experimental data. These results point to the predictive power of this multiscale viewpoint, which aims to integrate the quantitative kinetic description of large interaction networks typical of system biology approaches with a focus on the quality of a few, key, molecular interactions typical of structural biology ones.

The Interplay of Cholesterol and Ligand Binding in hTSPO from Classical Molecular Dynamics Simulations.

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol-binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N-terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C-terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150-X-Y152-X(3)-R156 and R135-X(2)-Y138-X(2)-L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC-resembling motif that we observed in TM I (residues L17-X(2)-F20-X(3)-R24) and to CRAC in TM V. We expect that the CRAC-like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.

Effect of the Heteroaromatic Antenna on the Binding of Chiral Eu(III) Complexes to Bovine Serum Albumin.

The cationic enantiopure (R,R) and luminescent Eu(III) complex [Eu(bisoQcd)(H2O)2] OTf (with bisoQcd = N,N'-bis(2-isoquinolinmethyl)-trans-1,2-diaminocyclohexane N,N'-diacetate and OTf = triflate) was synthesized and characterized. At physiological pH, the 1:1 [Eu(bisoQcd)(H2O)2]+ species, possessing two water molecules in the inner coordination sphere, is largely dominant. The interaction with bovine serum albumin (BSA) was studied by means of several experimental techniques, such as luminescence spectroscopy, isothermal titration calorimetry (ITC), molecular docking (MD), and molecular dynamics simulations (MDS). In this direction, a ligand competition study was also performed by using three clinically established drugs (i.e., ibuprofen, warfarin, and digitoxin). The nature of this interaction is strongly affected by the type of the involved heteroaromatic antenna in the Eu(III) complexes. In fact, the presence of isoquinoline rings drives the corresponding complex toward the protein superficial area containing the tryptophan residue 134 (Trp134). As the main consequence, the metal center undergoes the loss of one water molecule upon interaction with the side chain of a glutamic acid residue. On the other hand, the similar complex containing pyridine rings ([Eu(bpcd)(H2O)2]Cl with bpcd = N,N'-bis(2-pyridylmethyl)-trans-1,2-diaminocyclohexane N,N'-diacetate) interacts more weakly with the protein in a different superficial cavity, without losing the coordinated water molecules.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach.

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Genetic and Clinical Heterogeneity in Thirteen New Cases with Aceruloplasminemia. Atypical Anemia as a Clue for an Early Diagnosis.

Aceruloplasminemia is a rare autosomal recessive genetic disease characterized by mild microcytic anemia, diabetes, retinopathy, liver disease, and progressive neurological symptoms due to iron accumulation in pancreas, retina, liver, and brain. The disease is caused by mutations in the Ceruloplasmin (CP) gene that produce a strong reduction or absence of ceruloplasmin ferroxidase activity, leading to an impairment of iron metabolism. Most patients described so far are from Japan. Prompt diagnosis and therapy are crucial to prevent neurological complications since, once established, they are usually irreversible. Here, we describe the largest series of non-Japanese patients with aceruloplasminemia published so far, including 13 individuals from 11 families carrying 13 mutations in the CP gene (7 missense, 3 frameshifts, and 3 splicing mutations), 10 of which are novel. All missense mutations were studied by computational modeling. Clinical manifestations were heterogeneous, but anemia, often but not necessarily microcytic, was frequently the earliest one. This study confirms the clinical and genetic heterogeneity of aceruloplasminemia, a disease expected to be increasingly diagnosed in the Next-Generation Sequencing (NGS) era. Unexplained anemia with low transferrin saturation and high ferritin levels without inflammation should prompt the suspicion of aceruloplasminemia, which can be easily confirmed by low serum ceruloplasmin levels. Collaborative joint efforts are needed to better understand the pathophysiology of this potentially disabling disease.

Impact of Cholesterol on the Stability of Monomeric and Dimeric Forms of the Translocator Protein TSPO: A Molecular Simulation Study.

The translocator protein (TSPO) is a transmembrane protein present across the three domains of life. Its functional quaternary structure consists of one or more subunits. In mice, the dimer-to-monomer equilibrium is shifted in vitro towards the monomer by adding cholesterol, a natural component of mammalian membranes. Here, we present a coarse-grained molecular dynamics study on the mouse protein in the presence of a physiological content and of an excess of cholesterol. The latter turns out to weaken the interfaces of the dimer by clusterizing mostly at the inter-monomeric space and pushing the contact residues apart. It also increases the compactness and the rigidity of the monomer. These two factors might play a role for the experimentally observed incremented stability of the monomeric form with increased content of cholesterol. Comparison with simulations on bacterial proteins suggests that the effect of cholesterol is much less pronounced for the latter than for the mouse protein.

CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations.

Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations.

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Unifying view of mechanical and functional hotspots across class A GPCRs.

G protein-coupled receptors (GPCRs) are the largest superfamily of signaling proteins. Their activation process is accompanied by conformational changes that have not yet been fully uncovered. Here, we carry out a novel comparative analysis of internal structural fluctuations across a variety of receptors from class A GPCRs, which currently has the richest structural coverage. We infer the local mechanical couplings underpinning the receptors' functional dynamics and finally identify those amino acids whose virtual deletion causes a significant softening of the mechanical network. The relevance of these amino acids is demonstrated by their overlap with those known to be crucial for GPCR function, based on static structural criteria. The differences with the latter set allow us to identify those sites whose functional role is more clearly detected by considering dynamical and mechanical properties. Of these sites with a genuine mechanical/dynamical character, the top ranking is amino acid 7x52, a previously unexplored, and experimentally verifiable key site for GPCR conformational response to ligand binding.

Hypomorphic mutations in POLR3A are a frequent cause of sporadic and recessive spastic ataxia.

Despite extensive efforts, half of patients with rare movement disorders such as hereditary spastic paraplegias and cerebellar ataxias remain genetically unexplained, implicating novel genes and unrecognized mutations in known genes. Non-coding DNA variants are suspected to account for a substantial part of undiscovered causes of rare diseases. Here we identified mutations located deep in introns of POLR3A to be a frequent cause of hereditary spastic paraplegia and cerebellar ataxia. First, whole-exome sequencing findings in a recessive spastic ataxia family turned our attention to intronic variants in POLR3A, a gene previously associated with hypomyelinating leukodystrophy type 7. Next, we screened a cohort of hereditary spastic paraplegia and cerebellar ataxia cases (n = 618) for mutations in POLR3A and identified compound heterozygous POLR3A mutations in ∼3.1% of index cases. Interestingly, >80% of POLR3A mutation carriers presented the same deep-intronic mutation (c.1909+22G>A), which activates a cryptic splice site in a tissue and stage of development-specific manner and leads to a novel distinct and uniform phenotype. The phenotype is characterized by adolescent-onset progressive spastic ataxia with frequent occurrence of tremor, involvement of the central sensory tracts and dental problems (hypodontia, early onset of severe and aggressive periodontal disease). Instead of the typical hypomyelination magnetic resonance imaging pattern associated with classical POLR3A mutations, cases carrying c.1909+22G>A demonstrated hyperintensities along the superior cerebellar peduncles. These hyperintensities may represent the structural correlate to the cerebellar symptoms observed in these patients. The associated c.1909+22G>A variant was significantly enriched in 1139 cases with spastic ataxia-related phenotypes as compared to unrelated neurological and non-neurological phenotypes and healthy controls (P = 1.3 × 10-4). In this study we demonstrate that (i) autosomal-recessive mutations in POLR3A are a frequent cause of hereditary spastic ataxias, accounting for about 3% of hitherto genetically unclassified autosomal recessive and sporadic cases; and (ii) hypomyelination is frequently absent in POLR3A-related syndromes, especially when intronic mutations are present, and thus can no longer be considered as the unifying feature of POLR3A disease. Furthermore, our results demonstrate that substantial progress in revealing the causes of Mendelian diseases can be made by exploring the non-coding sequences of the human genome.

Structural Prediction of the Dimeric Form of the Mammalian Translocator Membrane Protein TSPO: A Key Target for Brain Diagnostics.

Positron emission tomography (PET) radioligands targeting the human translocator membrane protein (TSPO) are broadly used for the investigations of neuroinflammatory conditions associated with neurological disorders. Structural information on the mammalian protein homodimers-the suggested functional state of the protein-is limited to a solid-state nuclear magnetic resonance (NMR) study and to a model based on the previously-deposited solution NMR structure of the monomeric mouse protein. Computational studies performed here suggest that the NMR-solved structure in the presence of detergents is not prone to dimer formation and is furthermore unstable in its native membrane environment. We, therefore, propose a new model of the functionally-relevant dimeric form of the mouse protein, based on a prokaryotic homologue. The model, fully consistent with solid-state NMR data, is very different from the previous predictions. Hence, it provides, for the first time, structural insights into this pharmaceutically-important target which are fully consistent with experimental data.

Role of Extracellular Loops and Membrane Lipids for Ligand Recognition in the Neuronal Adenosine Receptor Type 2A: An Enhanced Sampling Simulation Study.

Human G-protein coupled receptors (GPCRs) are important targets for pharmaceutical intervention against neurological diseases. Here, we use molecular simulation to investigate the key step in ligand recognition governed by the extracellular domains in the neuronal adenosine receptor type 2A (hA2AR), a target for neuroprotective compounds. The ligand is the high-affinity antagonist (4-(2-(7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-ylamino)ethyl)phenol), embedded in a neuronal membrane mimic environment. Free energy calculations, based on well-tempered metadynamics, reproduce the experimentally measured binding affinity. The results are consistent with the available mutagenesis studies. The calculations identify a vestibular binding site, where lipids molecules can actively participate to stabilize ligand binding. Bioinformatic analyses suggest that such vestibular binding site and, in particular, the second extracellular loop, might drive the ligand toward the orthosteric binding pocket, possibly by allosteric modulation. Taken together, these findings point to a fundamental role of the interaction between extracellular loops and membrane lipids for ligands' molecular recognition and ligand design in hA2AR.

Identification of new BMP6 pro-peptide mutations in patients with iron overload.

Hereditary Hemochromatosis (HH) is a genetically heterogeneous disorder caused by mutations in at least five different genes (HFE, HJV, TFR2, SLC40A1, HAMP) involved in the production or activity of the liver hormone hepcidin, a key regulator of systemic iron homeostasis. Nevertheless, patients with an HH-like phenotype that remains completely/partially unexplained despite extensive sequencing of known genes are not infrequently seen at referral centers, suggesting a role of still unknown genetic factors. A compelling candidate is Bone Morphogenetic Protein 6 (BMP6), which acts as a major activator of the BMP-SMAD signaling pathway, ultimately leading to the upregulation of hepcidin gene transcription. A recent seminal study by French authors has described three heterozygous missense mutations in BMP6 associated with mild to moderate late-onset iron overload (IO). Using an updated next-generation sequencing (NGS)-based genetic test in IO patients negative for the classical HFE p.Cys282Tyr mutation, we found three BMP6 heterozygous missense mutations in four patients from three different families. One mutation (p.Leu96Pro) has already been described and proven to be functional. The other two (p.Glu112Gln, p.Arg257His) were novel, and both were located in the pro-peptide domain known to be crucial for appropriate BMP6 processing and secretion. In silico modeling also showed results consistent with their pathogenetic role. The patients' clinical phenotypes were similar to that of other patients with BMP6-related IO recently described. Our results independently add further evidence to the role of BMP6 mutations as likely contributing factors to late-onset moderate IO unrelated to mutations in the established five HH genes.

A comprehensive picture of the mutations associated with aromatic amino acid decarboxylase deficiency: from molecular mechanisms to therapy implications.

Dopa decarboxylase (DDC), or aromatic amino acid decarboxylase (AADC), is a pyridoxal 5'-phosphate enzyme responsible for the production of the neurotransmitters dopamine and serotonin. Deficit of this enzyme causes AADC deficiency, an inherited neurometabolic disorder. To date, 18 missense homozygous mutations have been identified through genetic screening in ∼80 patients. However, little is known about the mechanism(s) by which mutations cause disease. Here we investigated the impact of these pathogenic mutations and of an artificial one on the conformation and the activity of wild-type DDC by a combined approach of bioinformatic, spectroscopic and kinetic analyses. All mutations reduce the kcat value, and, except the mutation R347Q, alter the tertiary structure, as revealed by an increased hydrophobic surface and a decreased near-UV circular dichroism signal. The integrated analysis of the structural and functional consequences of each mutation strongly suggests that the reason underlying the pathogenicity of the majority of disease-causing mutations is the incorrect apo-holo conversion. In fact, the most remarkable effects are seen upon mutation of residues His70, His72, Tyr79, Phe80, Pro81, Arg462 and Arg447 mapping to or directly interacting with loop1, a structural key element involved in the apo-holo switch. Instead, different mechanisms are responsible for the pathogenicity of R347Q, a mere catalytic mutation, and of L38P and A110Q mutations causing structural-functional defects. These are due to local perturbation transmitted to the active site, as predicted by molecular dynamic analyses. Overall, the results not only give comprehensive molecular insights into AADC deficiency, but also provide an experimental framework to suggest appropriate therapeutic treatments.