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.

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.

Structures of MERS-CoV macro domain in aqueous solution with dynamics: Impacts of parallel tempering simulation techniques and CHARMM36m and AMBER99SB force field parameters.

A novel virus, severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), causing coronavirus disease 2019 (COVID-19) worldwide appeared in 2019. Detailed scientific knowledge of the members of the Coronaviridae family, including the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is currently lacking. Structural studies of the MERS-CoV proteins in the current literature are extremely limited. We present here detailed characterization of the structural properties of MERS-CoV macro domain in aqueous solution. Additionally, we studied the impacts of chosen force field parameters and parallel tempering simulation techniques on the predicted structural properties of MERS-CoV macro domain in aqueous solution. For this purpose, we conducted extensive Hamiltonian-replica exchange molecular dynamics simulations and Temperature-replica exchange molecular dynamics simulations using the CHARMM36m and AMBER99SB parameters for the macro domain. This study shows that the predicted secondary structure properties including their propensities depend on the chosen simulation technique and force field parameter. We perform structural clustering based on the radius of gyration and end-to-end distance of MERS-CoV macro domain in aqueous solution. We also report and analyze the residue-level intrinsic disorder features, flexibility and secondary structure. Furthermore, we study the propensities of this macro domain for protein-protein interactions and for the RNA and DNA binding. Overall, results are in agreement with available nuclear magnetic resonance spectroscopy findings and present more detailed insights into the structural properties of MERS CoV macro domain in aqueous solution. All in all, we present the structural properties of the aqueous MERS-CoV macro domain using different parallel tempering simulation techniques, force field parameters and bioinformatics tools.

Transition Metal Ion Interactions with Disordered Amyloid-ß Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies.

Monomers and oligomers of the amyloid-ß peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-ß aggregation. Due to the disordered nature of amyloid-ß, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-ß monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-ß obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-ß structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-ß and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-ß monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-ß obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-ß. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

Alanine Scanning Effects on the Biochemical and Biophysical Properties of Intrinsically Disordered Proteins: A Case Study of the Histidine to Alanine Mutations in Amyloid-ß42.

Alanine scanning is a tool in molecular biology that is commonly used to evaluate the contribution of a specific amino acid residue to the stability and function of a protein. Additionally, this tool is also used to understand whether the side chain of a specific amino acid residue plays a role in the protein's bioactivity. Furthermore, computational alanine scanning methods are utilized to predict the thermodynamic properties of proteins. These studies are utilized with the assumption that the biochemical and biophysical properties of a protein do not change with alanine scanning. Our study was dedicated to analyze the effect of alanine scanning on the biochemical and biophysical properties of intrinsically disordered proteins. To this end, we studied the impact of widely used histidine to alanine mutations in amyloid-ß (Aß). We found that the secondary and tertiary contacts, salt bridge formations, and thermodynamic properties, as well as disorder propensities and aggregation predisposition of Aß, are impacted by the single and triple point histidine to alanine mutations. Experimental and computational studies employing the alanine scanning technique for mutating histidine to alanine in the analysis of intrinsically disordered proteins have to consider these effects.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology.

Amyloid-ß and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-ß and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-ß and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-ß and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-ß and α-synuclein, such as ß-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

Structures prediction and replica exchange molecular dynamics simulations of α-synuclein: A case study for intrinsically disordered proteins.

In recent years, a variety of three-dimensional structure prediction tools, including AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold, have been employed in the investigation of intrinsically disordered proteins. However, a comprehensive validation of these tools specifically for intrinsically disordered proteins has yet to be conducted. In this study, we utilize AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold to predict the structure of a model intrinsically disordered α-synuclein protein. Additionally, extensive replica exchange molecular dynamics simulations of the intrinsically disordered protein are conducted. The resulting structures from both structure prediction tools and replica exchange molecular dynamics simulations are analyzed for radius of gyration, secondary and tertiary structure properties, as well as Cα and Hα chemical shift values. A comparison of the obtained results with experimental data reveals that replica exchange molecular dynamics simulations provide results in excellent agreement with experimental observations. However, none of the structure prediction tools utilized in this study can fully capture the structural characteristics of the model intrinsically disordered protein. This study shows that a cluster of ensembles are required for intrinsically disordered proteins. Artificial-intelligence based structure prediction tools such as AlphaFold3 and C-I-TASSER could benefit from stochastic sampling or Monte Carlo simulations for generating an ensemble of structures for intrinsically disordered proteins.

Liquid-Liquid Phase Separation Associated with Intrinsically Disordered Proteins: Experimental and Computational Tools.

The phenomenon of Liquid-Liquid Phase Separation (LLPS) serves as a vital mechanism for the spatial organization of biomolecules, significantly influencing the elementary processes within the cellular milieu. Intrinsically disordered proteins, or proteins endowed with intrinsically disordered regions, are pivotal in driving this biophysical process, thereby dictating the formation of non-membranous cellular compartments. Compelling evidence has linked aberrations in LLPS to the pathogenesis of various neurodegenerative diseases, underscored by the disordered proteins' proclivity to form pathological aggregates. This study meticulously evaluates the arsenal of contemporary experimental and computational methodologies dedicated to the examination of intrinsically disordered proteins within the context of LLPS. Through a discerning discourse on the capabilities and constraints of these investigative techniques, we unravel the intricate contributions of these ubiquitous proteins to LLPS and neurodegeneration. Moreover, we project a future trajectory for the field, contemplating on innovative research tools and their potential to elucidate the underlying mechanisms of LLPS, with the ultimate goal of fostering new therapeutic avenues for combating neurodegenerative disorders.

Current Stage and Future Perspectives for Homology Modeling, Molecular Dynamics Simulations, Machine Learning with Molecular Dynamics, and Quantum Computing for Intrinsically Disordered Proteins and Proteins with Intrinsically Disordered Regions.

The structural ensembles of intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) cannot be easily characterized using conventional experimental techniques. Computational techniques complement experiments and provide useful insights into the structural ensembles of IDPs and proteins with IDRs. Herein, we discuss computational techniques such as homology modeling, molecular dynamics simulations, machine learning with molecular dynamics, and quantum computing that can be applied to the studies of IDPs and hybrid proteins with IDRs. We also provide useful future perspectives for computational techniques that can be applied to IDPs and hybrid proteins containing ordered domains and IDRs.

Structural Properties of Rat Intestinal Fatty Acid-Binding Protein with its Dynamics: Insights into Intrinsic Disorder.

BACKGROUND: The rat intestinal fatty acid-binding protein (I-FABP) is expressed in the small intestine and is involved in the absorption and transport of dietary fatty acids. It is used as a marker for intestinal injury and is associated with various gastrointestinal disorders. I-FABP has been studied extensively using conventional experimental and computational techniques. However, the detection of intrinsically disordered regions requires the application of special sampling molecular dynamics simulations along with certain bioinformatics because conventional computational and experimental studies face challenges in identifying the features of intrinsic disorder. METHOD: Replica exchange molecular dynamics simulations were conducted along with bioinformatics studies to gain deeper insights into the structural properties of I-FABP. Specifically, the Cα and Hα chemical shift values werecalculated, and the findings were compared to the experiments. Furthermore, secondary and tertiary structure properties were also calculated, and the protein was clustered using k-means clustering. The end-to-end distance and radius of gyration values were reported for the protein in an aqueous solution medium. In addition, its disorder tendency was studied using various bioinformatics tools. RESULTS AND CONCLUSION: It was reported that I-FABP is a flexible protein with regions that demonstrate intrinsic disorder characteristics. This flexibility and intrinsic disorder characteristics of I-- FABP may be related to its nature in ligand binding processes.

Paving the Way for Synthetic Intrinsically Disordered Polymers for Soft Robotics.

Nature is full of examples of processes that, through evolution, have been perfected over the ages to effectively use matter and sustain life. Here, we present our strategies for designing intrinsically disordered smart polymers for soft robotics applications that are bio-inspired by intrinsically disordered proteins. Bio-inspired intrinsically disordered smart and soft polymers designed using our deep understanding of intrinsically disordered proteins have the potential to open new avenues in soft robotics. Together with other desirable traits, such as robustness, dynamic self-organization, and self-healing abilities, these systems possess ideal characteristics that human-made formations strive for but often fail to achieve. Our main aim is to develop materials for soft robotics applications bio-inspired by intrinsically disordered proteins to address what we see as the largest current barriers in the practical deployment of future soft robotics in various areas, including defense. Much of the current literature has focused on the de novo synthesis of tailor-made polymers to perform specific functions. With bio-inspired polymers, the complexity of protein folding mechanisms has limited the ability of researchers to reliably engineer specific structures. Unlike existing studies, our work is focused on utilizing the high flexibility of intrinsically disordered proteins and their self-organization characteristics using synthetic quasi-foldamers.

Frontotemporal Dementia-Related V57E Mutation Impairs Mitochondrial Function and Alters the Structural Properties of CHCHD10.

The V57E pathological variant of the mitochondrial coiled-coil-helix-coiled-coil-helix domain-containing protein 10 (CHCHD10) plays a role in frontotemporal dementia. The wild-type and V57E mutant CHCHD10 proteins contain intrinsically disordered regions, and therefore, these regions hampered structural characterization of these proteins using conventional experimental tools. For the first time in the literature, we represent that the V57E mutation is pathogenic to mitochondria as it increases mitochondrial superoxide and impairs mitochondrial respiration. In addition, we represent here the structural ensemble properties of the V57E mutant CHCHD10 and describe the impacts of V57E mutation on the structural ensembles of wild-type CHCHD10 in aqueous solution. We conducted experimental and computational studies for this research. Namely, MitoSOX Red staining and Seahorse Mito Stress experiments, atomic force microscopy measurements, bioinformatics, homology modeling, and multiple-run molecular dynamics simulation computational studies were conducted. Our experiments show that the V57E mutation results in mitochondrial dysfunction, and our computational studies present that the structural ensemble properties of wild-type CHCHD10 are impacted by the frontotemporal dementia-associated V57E genetic mutation.

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.

Intrinsically Disordered Synthetic Polymers in Biomedical Applications.

In biology and medicine, intrinsically disordered synthetic polymers bio-mimicking intrinsically disordered proteins, which lack stable three-dimensional structures, possess high structural/conformational flexibility. They are prone to self-organization and can be extremely useful in various biomedical applications. Among such applications, intrinsically disordered synthetic polymers can have potential usage in drug delivery, organ transplantation, artificial organ design, and immune compatibility. The designing of new syntheses and characterization mechanisms is currently required to provide the lacking intrinsically disordered synthetic polymers for biomedical applications bio-mimicked using intrinsically disordered proteins. Here, we present our strategies for designing intrinsically disordered synthetic polymers for biomedical applications based on bio-mimicking intrinsically disordered proteins.

Intrinsically disordered proteins and proteins with intrinsically disordered regions in neurodegenerative diseases.

Many different intrinsically disordered proteins and proteins with intrinsically disordered regions are associated with neurodegenerative diseases. These types of proteins including amyloid-ß, tau, α-synuclein, CHCHD2, CHCHD10, and G-protein coupled receptors are increasingly becoming evaluated as potential drug targets in the pharmaceutical-based treatment approaches. Here, we focus on the neurobiology of this class of proteins, which lie at the center of numerous neurodegenerative diseases, such as Alzheimer's and Parkinson's diseases, Huntington's disease, amyotrophic lateral sclerosis, frontotemporal dementia, Charcot-Marie-Tooth diseases, spinal muscular atrophy, and mitochondrial myopathy. Furthermore, we discuss the current treatment design strategies involving intrinsically disordered proteins and proteins with intrinsically disordered regions in neurodegenerative diseases. In addition, we emphasize that although the G-protein coupled receptors are traditionally investigated using structural biology-based models and approaches, current studies show that these receptors are proteins with intrinsically disordered regions and therefore they require new ways for their analysis.

From Quantum Mechanics, Classical Mechanics, and Bioinformatics to Artificial Intelligence Studies in Neurodegenerative Diseases.

The amyloid ß-protein is an intrinsically disordered protein that has the potential to assemble into myriad structures, including oligomers and fibrils. These structures are neurotoxic and are thought to initiate a cascade of events leading to Alzheimer's disease. Understanding this pathogenetic process and elucidating targets for drug therapy depends on elucidation of the structural dynamics of Aß assembly. In this chapter, we describe work packages required to determine the three-dimensional structures of Aß and of smaller bioactive fragments thereof, which may be important in AD pathogenesis. These packages include density functional theory, Car-Parrinello molecular dynamics simulations, temperature-dependent replica exchange molecular dynamics simulations, disorder predictors based on bioinformatics, and neural network deep learning.

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.

Insights into the structural properties of SARS-CoV-2 main protease.

SARS-CoV-2 is the infectious agent responsible for the coronavirus disease since 2019, which is the viral pneumonia pandemic worldwide. The structural knowledge on SARS-CoV-2 is rather limited. These limitations are also applicable to one of the most attractive drug targets of SARS-CoV-2 proteins - namely, main protease Mpro, also known as 3C-like protease (3CLpro). This protein is crucial for the processing of the viral polyproteins and plays crucial roles in interfering viral replication and transcription. In fact, although the crystal structure of this protein with an inhibitor was solved, Mpro conformational dynamics in aqueous solution is usually studied by molecular dynamics simulations without special sampling techniques. We conducted replica exchange molecular dynamics simulations on Mpro in water and report the dynamic structures of Mpro in an aqueous environment including root mean square fluctuations, secondary structure properties, radius of gyration, and end-to-end distances, chemical shift values, intrinsic disorder characteristics of Mpro and its active sites with a set of computational tools. The active sites we found coincide with the currently known sites and include a new interface for interaction with a protein partner.

Secondary structure dependence on simulation techniques and force field parameters: from disordered to ordered proteins.

Computer simulations are used for identifying the secondary structure properties of ordered and disordered proteins. However, our recent studies showed that the chosen computer simulation protocol, simulation technique, and force field parameter set for a disordered protein impact its predicted secondary structure properties. Here, we compare the outcome from computer simulations utilizing molecular dynamics simulations without parallel tempering techniques using various force field parameter sets and temperature-replica exchange molecular dynamics simulations both for a model ordered and two model disordered proteins. Specifically, the model ordered protein is the third IgG-binding domain of Protein G (GB3) and the two model disordered proteins are amyloid-ß(1-40) and α-synuclein in water. Our findings clearly indicate that temperature-replica exchange molecular dynamics simulations and molecular dynamics simulations without special sampling techniques yield similar results for the ordered GB3 protein whereas such agreement between simulation techniques using various force field parameter sets could not be obtained for disordered proteins. These findings clearly indicate that a consensus has to be reached via further development in computer simulation technique and force field parameter sets for disordered proteins.

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.