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.

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. METHODS: 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 IFABP may be related to its nature in ligand binding processes.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

Experiments face challenges in the analysis of intrinsically disordered proteins in solution due to fast conformational changes and enhanced aggregation propensity. Computational studies complement experiments, being widely used in the analyses of intrinsically disordered proteins, especially those positioned at the centers of neurodegenerative diseases. However, recent investigations - including our own - revealed that computer simulations face significant challenges and limitations themselves. In this review, we introduced and discussed some of the scientific challenges and limitations of computational studies conducted on intrinsically disordered proteins. We also outlined the importance of future developments in the areas of computational chemistry and computational physics that would be needed for generating more accurate data for intrinsically disordered proteins from computer simulations. Additional theoretical strategies that can be developed are discussed herein.

Epitope region identification challenges of intrinsically disordered proteins in neurodegenerative diseases: Secondary structure dependence of α-synuclein on simulation techniques and force field parameters.

Due to fast aggregation processes of many disordered proteins in neurodegenerative diseases, it is difficult to study their epitope regions at the monomeric and oligomeric levels. Computer simulations complement experiments and have been used to identify the epitope regions of proteins. Residues that adopt ß-sheet conformation play a central role in the oligomerization and aggregation mechanisms of such proteins, including α-synuclein, which is at the center of Parkinson's and Alzheimer's diseases. In this study, we simulated the monomeric α-synuclein protein in an aqueous environment to evaluate its secondary structure properties, including ß-sheet propensity, and radius of gyration by replica exchange molecular dynamics simulations. We also obtained the molecular dynamics simulation trajectories of α-synuclein that were conducted using various force field parameters by the David E. Shaw group. Using these trajectories, we calculated the impacts of force field parameters on α-synuclein secondary structure properties and radius of gyration values and obtained results are compared with our data from REMD simulations. This study shows that the chosen force field parameters and computer simulation techniques effect the predicted secondary structure properties and radius of gyration values of α-synuclein in water. Herewith, we illustrate the challenges in epitope region identification of intrinsically disordered proteins in neurodegenerative diseases by current computer simulations.