Current Challenges and Limitations in the Studies of Intrinsically Disordered Proteins in Neurodegenerative Diseases by Computer Simulations

The impacts of the mitochondrial myopathy-associated G58R mutation on the dynamic structural properties of CHCHD10

The mitochondria are responsible for producing energy within the cell, and in mitochondrial myopathy, there is a defect in the energy production process. The CHCHD10 gene codes for a protein called coiled-coil-helix-coiled-coil-helix domain-containing protein 10 (CHCHD10), which is found in the mitochondria and is involved in the regulation of mitochondrial function. G58R mutation has been shown to disrupt the normal function of CHCHD10, leading to mitochondrial dysfunction and ultimately to the development of mitochondrial myopathy. The structures of G58R mutant CHCHD10 and how G58R mutation impacts the wild-type CHCHD10 protein at the monomeric level are unknown. To address this problem, we conducted homology modeling, multiple run molecular dynamics simulations and bioinformatics calculations. We represent herein the structural ensemble properties of the G58R mutant CHCHD10 (CHCHD10G58R) in aqueous solution. Moreover, we describe the impacts of G58R mutation on the structural ensembles of wild-type CHCHD10 (CHCHD10WT) in aqueous solution. The dynamics properties as well as structural properties of CHCHD10WT are impacted by the mitochondrial myopathy-related G58R mutation. Specifically, the secondary and tertiary structure properties, root mean square fluctuations, Ramachandran diagrams and results from principal component analysis demonstrate that the CHCHD10WT and CHCHD10G58R proteins possess different structural ensemble characteristics and describe the impacts of G58R mutation on CHCHD10WT. These findings may be helpful for designing new treatments for mitochondrial myopathy.Communicated by Ramaswamy H. Sarma.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Balanced Three-Point Water Model OPC3-B for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) play a critical role in many biological processes. Due to the inherent structural flexibility of IDPs, experimental methods present significant challenges for sampling their conformational information at the atomic level. Therefore, molecular dynamics (MD) simulations have emerged as the primary tools for modeling IDPs whose accuracy depend on force field and water model. To enhance the accuracy of physical modeling of IDPs, several force fields have been developed. However, current water models lack precision and underestimate the interaction between water molecules and proteins. Here, we used Monte-Carlo re-weighting method to re-parameterize a three-point water model based on OPC3 for IDPs (named OPC3-B). We benchmarked the performance of OPC3-B compared with nine different water models for 10 IDPs and three ordered proteins. The results indicate that the performance of OPC3-B is better than other water models for both IDPs and ordered proteins. At the same time, OPC3-B possess the power of transferability with other force field to simulate IDPs. This newly developed water model can be used to insight into the research of sequence-disordered-function paradigm for IDPs.

Conformational dynamics of A30G α-synuclein that causes familial Parkinson disease

The first gene shown to be responsible for autosomal-dominant Parkinson's disease (PD) is the SNCA gene, which encodes for alpha synuclein (α-Syn). Recently, a novel heterozygous A30G mutation of the SNCA gene associated with familial PD has been reported. However, little research has been done on how the A30G mutation affects the structure of α-Syn. So, using atomistic molecular dynamics (MD) simulation, we demonstrate here the key structural characteristics of A30G α-Syn in the free monomer form and membrane associated state. From the MD trajectory analysis, the structure of A30G α-Syn was noticed to exhibit rapid conformational change, increase in backbone flexibility near the site of mutation and decrease in α-helical propensity. The typical torsion angles in residues (Val26 and Glu28) near the mutation site were observed to deviate significantly in A30G α-Syn. In the case of membrane bound A30G α-Syn, the regions that were submerged in the lipid bilayer (N-helix (3-37) and turn region (38-44)) found to contain higher helical content than the elevated region above the lipid surface. The bending angle in the helix-N and helix-C regions were noticed to be relatively higher in the free form of A30G α-Syn (38.50) than in the membrane bound form (370). The A30G mutation in α-Syn was predicted to have an impact on the stability and function of the protein based on ΔΔG values obtained from the online servers. Our results demonstrate that the A30G mutation in α-Syn altered the protein's α-helical structure and slightly altered the membrane binding.Communicated by Ramaswamy H. Sarma.

Phosphorylation of Tau R2 Repeat Destabilizes Its Binding to Microtubules: A Molecular Dynamics Simulation Study

Phosphorylation, the most popular post-translational modification of tau protein, plays an important role in regulating tau physiological functions. However, aberrant phosphorylation attenuates the binding affinity of tau to a microtubule (MT), resulting in MT destabilization followed by accumulation of neurofibrillary tangles in the brain. There are in total 85 potential phosphorylation sites in a full-length tau protein, and about half of them are abnormal as they occur in tau of Alzheimer's disease (AD) brain only. In this work, we investigated the impact of abnormal Ser289, Ser293, and Ser289/Ser293 phosphorylation on tau R2-MT binding and the conformation of tau R2 using molecular dynamics simulation. We found that the phosphorylation significantly affected R2-MT interaction and reduced the binding affinity of tau R2 peptides to MTs. Free energy decomposition analysis suggested that the post-translational modified residues themselves made a significant contribution to destabilize tau repeat R2-MT binding. Therefore, the phosphorylation may attenuate the binding affinity of tau to MTs. Additionally, the phosphorylation also enhanced helix-coil transition of monomeric R2 peptides, which may result in the acceleration of tau aggregation. Since these phosphorylated sites have not been examined in previous experimental studies, our finding through all-atom molecular dynamics simulations and free energy analysis can inspire experimental scientists to investigate the impact of the phosphorylation on MT binding and aggregation of full-length tau and the pathological roles of the phosphorylation at those sites in AD development through in vitro/in vivo assays.

Effects of the Jokela type of spinal muscular atrophy-related G66V mutation on the structural ensemble characteristics of CHCHD10

The G66V pathological variant of the coiled-coil-helix-coiled-coil-helix domain-containing protein 10 (CHCHD10), mitochondrial, plays a role in Jokela type spinal muscular atrophy. The wild-type and G66V mutant-type CHCHD10 proteins contain intrinsically disordered regions, and therefore, their structural ensemble studies have been experiencing difficulties using conventional tools. Here, we show our results regarding the first characterization of the structural ensemble characteristics of the G66V mutant form of CHCHD10 and the first comparison of these characteristics with the structural ensemble properties of wild-type CHCHD10. We find that the structural properties, potential of mean force surfaces, and principal component analysis show stark differences between these two proteins. These results are important for a better pathology, biochemistry and structural biology understanding of CHCHD10 and its G66V genetic variant and it is likely that these reported structural properties are important for designing more efficient treatments for the Jokela type of spinal muscular atrophy disease.

Personal Precise Force Field for Intrinsically Disordered and Ordered Proteins Based on Deep Learning

Intrinsically disordered proteins (IDPs) are proteins without a fixed three-dimensional (3D) structure under physiological conditions and are associated with Parkinson's disease, Alzheimer's disease, cancer, cardiovascular disease, amyloidosis, diabetes, and other diseases. Experimental methods can hardly capture the ensemble of diverse conformations for IDPs. Molecular dynamics (MD) simulations can sample continuous conformations that might provide a valuable complement to experimental data. However, the accuracy of MD simulations depends on the quality of force field. In particular, the evolutionary conservation and coevolution of IDPs introduce that current force fields could not precisely reproduce the conformation of IDPs. In order to improve the performance of force field, deep learning and reweighting methods were used to automatically generate personal force field parameters for intrinsically disordered and ordered proteins. At first, the deep learning method predicted more accuracy φ/ψ dihedral of residue than the previous method. Then, reweighting optimized the personal force field parameters for each residue. Finally, typical representative systems such as IDPs, structure protein, and fast-folding protein were used to evaluate this force field. The results indicate that two personal force field parameters (named PPFF1 and PPFF1_af2) could better reproduce the experimental observables than ff03CMAP force field. In summary, this strategy will provide feasibility for the development of precise personal force fields.

Combination of structure-based virtual screening, molecular docking and molecular dynamics approaches for the discovery of anti-prion fibril flavonoids

Prion diseases are a group of rare neurodegenerative diseases caused by the structural conversion of cellular prion into Scrapie prion resulting aggregated fibrils. Therapy of prion diseases has been developed for several decades, especially drug designs based on the structure of prion monomers. Unfortunately, none of the designed anti-prion drugs function well clinically. To fight against prion fibrils, a drug design based on the precise structure of mammalian prion fibrils is highly required. Fortunately, based on the advantage of newly advanced cryo-electron microscopy (cryo-EM) in the deconvolution of large complexes, three prion fibril structures were resolved in the last 2 years. Based on the cryo-EM solved prion fibril structures, we are able to find some molecules fighting against prion fibrils. Quercetin, one flavonoid molecule in the polyphenol group, has been found to disaggregate the prion fibrils in vitro. In this study, we performed the molecular docking and molecular dynamics simulation on quercetin-like molecules possessing pharmacological properties to evaluate the anti-prion ability of tested molecules. As a result, four quercetin-like molecules interact with prion fibril and decrease the β-strand content by converting some β-strands into loop and helical structures to disintegrate the existing fibril structure. The results of this study are significant in the treatment of prion diseases, and the approaches used in this study are applicable to other amyloid diseases.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Structures of the Wild-Type and S59L Mutant CHCHD10 Proteins Important in Amyotrophic Lateral Sclerosis-Frontotemporal Dementia

The S59L genetic mutation of the mitochondrial coiled-coil-helix-coiled-coil-helix domain-containing protein 10 (CHCHD10) is involved in the pathogenesis of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The wild-type and mutant forms of this protein contain intrinsically disordered regions, and their structural characterization has been facing challenges. Here, for the first time in the literature, we present the structural ensemble properties of the wild-type and S59L mutant form of CHCHD10 in an aqueous solution environment at the atomic level with dynamics. Even though available experiments suggested that the S59L mutation may not change the structure of the CHCHD10 protein, our structural analysis clearly shows that the structure of this protein is significantly affected by the S59L mutation. We present here the secondary structure components with their abundances per residue, the tertiary structure properties, the free energy surfaces based on the radius of gyration and end-to-end distance values, the Ramachandran plots, the quantity of intramolecular hydrogen bonds, and the principal component analysis results. These results may be crucial in designing more efficient treatment for ALS and FTD diseases.

Protein folding, misfolding, and un/non-folding: overview of the SP16 Session at the 20th IUPAB congress, 45th Annual Meeting of SBBf, and 50th Annual Meeting of SBBq

The classical and idealistic consideration of the protein universe as a set of well-ordered proteinaceous machines with unique spatial organizations conducting unique biological functions is changed because of the recognition that proteins can fold, misfold, or be disordered to different degree. Functional repertoires of intrinsically disordered proteins complement functions of ordered proteins, whereas protein misfolding and concomitant oligomerization, aggregation, and fibrillation are related to the pathogenesis of numerous human diseases. Flexible proteins can undergo liquid-liquid phase separation, which is at the heart of the biogenesis of numerous membrane-less organelles. Many of these aspects were highlighted in four talks delivered at the SP16 Session at the 20th IUPAB congress, 45th Annual Meeting of SBBf, and 50th Annual Meeting of SBBq.

An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure

Amyloid beta (Aβ) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aβ peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aβ experiments, as it allows to explore the binding mechanism between metal ions and Aβ oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aβ peptides. Aβ oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aβ binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aβ. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aβ dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.

Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective

Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.

Secondary structure dependence of amyloid-β(1-40) on simulation techniques and force field parameters

Our recent studies revealed that none of the selected widely used force field parameters and molecular dynamics simulation techniques yield structural properties for the intrinsically disordered α-synuclein that are in agreement with various experiments via testing different force field parameters. Here, we extend our studies on the secondary structure properties of the disordered amyloid-β(1-40) peptide in aqueous solution. For these purposes, we conducted extensive replica exchange molecular dynamics simulations and obtained extensive molecular dynamics simulation trajectories from David E. Shaw group. Specifically, these molecular dynamics simulations were conducted using various force field parameters and obtained results are compared to our replica exchange molecular dynamics simulations and experiments. In this study, we calculated the secondary structure abundances and radius of gyration values for amyloid-β(1-40) that were simulated using varying force field parameter sets and different simulation techniques. In addition, the intrinsic disorder propensity, as well as sequence-based secondary structure predisposition of amyloid-β(1-40) and compared the findings with the results obtained from molecular simulations using various force field parameters and different simulation techniques. Our studies clearly show that the epitope region identification of amyloid-β(1-40) depends on the chosen simulation technique and chosen force field parameters. Based on comparison with experiments, we find that best computational results in agreement with experiments are obtained using the a99sb*-ildn, charmm36m, and a99sb-disp parameters for the amyloid-β(1-40) peptide in molecular dynamics simulations without parallel tempering or via replica exchange molecular dynamics simulations.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.