Tyrosine Regulates β-Sheet Structure Formation in Amyloid-β42: A New Clustering Algorithm for Disordered Proteins

Tyrosine - a structural glue for hierarchical protein assembly

Protein self-assembly, guided by the interplay of sequence- and environment-dependent liquid-liquid phase separation (LLPS), constitutes a fundamental process in the assembly of numerous intrinsically disordered proteins. Heuristic examination of these proteins has underscored the role of tyrosine residues, evident in their conservation and pivotal involvement in initiating LLPS and subsequent liquid-solid phase transitions (LSPT). The development of tyrosine-templated constructs, designed to mimic their natural counterparts, emerges as a promising strategy for creating adaptive, self-assembling systems with diverse applications. This review explores the central role of tyrosine in orchestrating protein self-assembly, delving into key interactions and examining its potential in innovative applications, including responsive biomaterials and bioengineering.

Tuning the viscoelastic properties of peptide coacervates by single amino acid mutations and salt kosmotropicity

Coacervation, or liquid-liquid phase separation (LLPS) of biomacromolecules, is increasingly recognized to play an important role both intracellularly and in the extracellular space. Central questions that remain to be addressed are the links between the material properties of coacervates (condensates) and both the primary and the secondary structures of their constitutive building blocks. Short LLPS-prone peptides, such as GY23 variants explored in this study, are ideal model systems to investigate these links because simple sequence modifications and the chemical environment strongly affect the viscoelastic properties of coacervates. Herein, a systematic investigation of the structure/property relationships of peptide coacervates was conducted using GY23 variants, combining biophysical characterization (plate rheology and surface force apparatus, SFA) with secondary structure investigations by infrared (IR) and circular dichroism (CD) spectroscopy. Mutating specific residues into either more hydrophobic or more hydrophilic residues strongly regulates the viscoelastic properties of GY23 coacervates. Furthermore, the ionic strength and kosmotropic characteristics (Hofmeister series) of the buffer in which LLPS is induced also significantly impact the properties of formed coacervates. Structural investigations by CD and IR indicate a direct correlation between variations in properties induced by endogenous (peptide sequence) or exogenous (ionic strength, kosmotropic characteristics, aging) factors and the β-sheet content within coacervates. These findings provide valuable insights to rationally design short peptide coacervates with programmable materials properties that are increasingly used in biomedical applications.

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.

A computational study of metal ions interaction with amyloid-β 1-42 peptide structure in hyperpyrexia: Implications for Alzheimer disease

Given the current context of the SARS-CoV-19 pandemic, among the interfering risky factors with the Aβ peptide aggregation in the brains of Alzheimer's disease (AD) patients can be hyperpyrexia and increased intracranial pressure (ICP). According to our hypothesis on the relationship between hyperpyrexia and cognitive decline in AD, two models of Aβ peptides were used in this study: the structure of AD amyloid beta-peptide and near-atomic resolution fibril structures of the Aβ peptide. Therefore, the binding templates were constructed for Aβ peptide regions able to bind 9 different metal ions. The fragment transformation method was used for the structural comparison between Aβ chains. Molecular dynamics simulation (MDS) was applied using the Nose-Poincare-Anderson equation to generate a theoretically correct NPT (isothermal-isobaric ensemble). The smallest dissimilarities were observed in the case of Cu⁺ binding potential followed by Co²⁺, both with similar variation. Structural changes have also occurred as a result of the dynamic simulation. All these changes suggest an aggravating factor in both hyperpyretic and AD conditions. Our findings suggest that elevated temperature and increased intracranial pressure rise the effect of peptide aggregation, by converting α-helix motif to β-sheet and random coil conformation, which are related to the formation of senile plaques in AD brains.

Rad5 HIRAN domain: Structural insights into its interaction with ssDNA through molecular modeling approaches

The Rad5 protein is an SWI/SNF family ubiquitin ligase that contains an N-terminal HIRAN domain and a RING C3HC4 motif. The HIRAN domain is critical for recognition of the stalled replication fork during the replication process and acts as a sensor to initiate the damaged DNA checkpoint. It is a conserved domain widely distributed in eukaryotic organisms and is present in several DNA-binding proteins from all kingdoms. Here we showed that distant species have important differences in key residues that affect affinity for ssDNA. Based on these findings, we hypothesized that different HIRAN domains might affect fork reversal and translesion synthesis through different metabolic processes. To address this question, we predicted the tertiary structure of both yeast and human HIRAN domains using molecular modeling. Structural dynamics experiments showed that the yeast HIRAN domain exhibited higher structural denaturation than its human homolog, although both domains became stable in the presence of ssDNA. Analysis of atomic contacts revealed that a greater number of interactions between the ssDNA nucleotides and the Rad5 domain are electrostatic. Taken together, these results provide new insights into the molecular mechanism of the HIRAN domain of Rad5 and may guide us to further elucidate differences in the ancient eukaryotes HIRAN sequences and their DNA affinity. Communicated by Ramaswamy H. Sarma.

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.

Interactions of the Aβ(1-42) Peptide with Boron Nitride Nanoparticles of Varying Curvature in an Aqueous Medium: Different Pathways to Inhibit β-Sheet Formation

The aggregation of amyloid β (Aβ) peptide triggered by its conformational changes leads to the commonly known neurodegenerative disease of Alzheimer's. It is believed that the formation of β sheets of the peptide plays a key role in its aggregation and subsequent fibrillization. In the current study, we have investigated the interactions of the Aβ(1-42) peptide with boron nitride nanoparticles and the effects of the latter on conformational transitions of the peptide through a series of molecular dynamics simulations. In particular, the effects of curvature of the nanoparticle surface are studied by considering boron nitride nanotubes (BNNTs) of varying diameter and also a planar boron nitride nanosheet (BNNS). Altogether, the current study involves the generation and analysis of 9.5 μs of dynamical trajectories of peptide-BNNT/BNNS pairs in an aqueous medium. It is found that BN nanoparticles of different curvatures that are studied in the present work inhibit the conformational transition of the peptide to its β-sheet form. However, such an inhibition effect follows different pathways for BN nanoparticles of different curvatures. For the BNNT with the highest surface curvature, i.e., (3,3) BNNT, the nanoparticle is found to inhibit β-sheet formation by stabilizing the helical structure of the peptide, whereas for planar BNNS, the β-sheet formation is prevented by making more favorable pathways available for transitions of the peptide to conformations of random coils and turns. The BNNTs with intermediate curvatures are found to exhibit diverse pathways of their interactions with the peptide, but in all cases, essentially no formation of the β sheet is found whereas substantial β-sheet formation is observed for Aβ(1-42) in water in the absence of any nanoparticle. The current study shows that BN nanoparticles have the potential to act as effective tools to prevent amyloid formation from Aβ peptides.

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.

Green Fluorescent Tripeptide Nanostructures: Synergetic Effects of Oxidation and Hierarchical Assembly

Peptide-based fluorescent materials hold promise for applications in energy harvesting and biomedicine. One remaining challenge is to overcome the barrier of fluorescence red shift, especially for peptide probes merely made of natural amino acids. Here we demonstrated an about 100 nm fluorescence red shift using a tripeptide GYK. Under UV illumination or in Fenton reactions, the hydroxyl free radical-based oxidation crosslinks GYK tripeptides into dimers, which can further act as building blocks to hierarchically assemble into nanostructures of different sizes and finally can shift the fluorescence from blue to green color (∼510 nm). Such assemblies can form core-shell-like nanostructures through further crosslinking of their surface with additional GYK monomers, which can not only make the nanostructures more robust but also efficiently improve their quantum yields. This research will deepen our understanding of bioluminescence, which sheds light on various biomedical applications of peptide-based fluorescent probes.

Folding and self-assembly of short intrinsically disordered peptides and protein regions

Proteins and peptide fragments are highly relevant building blocks in self-assembly for nanostructures with plenty of applications. Intrinsically disordered proteins (IDPs) and protein regions (IDRs) are defined by the absence of a well-defined secondary structure, yet IDPs/IDRs show a significant biological activity. Experimental techniques and computational modelling procedures for the characterization of IDPs/IDRs are discussed. Directed self-assembly of IDPs/IDRs allows reaching a large variety of nanostructures. Hybrid materials based on the derivatives of IDPs/IDRs show a promising performance as alternative biocides and nanodrugs. Cell mimicking, in vivo compartmentalization, and bone regeneration are demonstrated for IDPs/IDRs in biotechnological applications. The exciting possibilities of IDPs/IDRs in nanotechnology with relevant biological applications are shown.

Enhancing the binding of the β-sheet breaker peptide LPFFD to the amyloid-β fibrils by aromatic modifications: A molecular dynamics simulation study

Alzheimer's is a fatal neurodegenerative disease for which there is no cure at present. The disease is characterized by the presence of plaques in the brains of a patient, which are composed mainly of aggregates of the amyloid-β peptide in the form of β-sheet fibrils. Here, we investigated the possibility of exploiting the superior binding ability of aromatic amino acids to a particular model of the amyloid-β fibrils. which is a difficult target for drug design. The β-sheet breaker peptide LPFFD was modified with aromatic amino acids and its binding to these fibrils was studied. We found that the orientation and the electrostatic complementarity of the modified peptide with respect to the fibrils played a crucial role in determining whether its binding was improved by the aromatic amino acids. The modified LPFFD peptides were able to bind to those fibril residues. which are important in the aggregation of amyloid-β peptides and thus can potentially inhibit the further aggregation of the amyloid-beta peptides by blocking their interactions. We found that the tryptophan modified LPFFD peptides had the best binding affinities. In most cases, the aromatic amino acids in the N-terminus of the modified peptides made more contacts with the fibrils than those in the C-terminus. We also found that increasing the aromatic content did not significantly improve the binding of the LPFFD peptide to the fibrils. Our study can serve as a basis for the design of novel peptide-based drugs for Alzheimer's disease in which aromatic interactions play an important role.

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.

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.

Effect of AmyTrap, an amyloid-β binding drug, on Aβ induced mitochondrial dysfunction and tau phosphorylation in cultured neuroblastoma cells

Alzheimer's Disease (AD) is the most common cause of dementia, affecting 25 million people worldwide. Accumulation of Amyloid-β (Aβ) in the mitochondria has been shown to adversely affect key enzymes including pyruvate dehydrogenase (PDH), succinate dehydrogenase (SDH), oxoglutarate dehydrogenase (OGDH). Accumulation of Aβ is also believed to increase Tau expression and pathology. Tau, in its toxic state, results in synaptic damage causing memory and cognitive dysfunction. We are developing a drug to treat AD namely AmyTrap. The active pharmacological ingredient is a retro inverso, Aβ-binding peptide which sequesters Aβ. We wanted to examine the effect of AmyTrap peptide on Aβ-induced mitochondrial dysfunction and Tau phosphorylation. Therefore, we treated SH-SY5Y neuroblastoma cells with wild-type Aβ, a mutant AβY10A, AmyTrap peptide (RI-peptide), or Aβ and RI-peptide for 72 h. The mutant AβY10A is known to impact the self-aggregating property of Aβ as this Tyr10 is essential for self-aggregation. As expected, AβY10A reversed PDH, OGDH and SDH dysfunction to near normal levels. Further, AβY10A successfully reversed Tau phosphorylation, suggesting that Tyr10 is also associated with Aβ-induced cytotoxicity. RI-peptide was able to significantly reverse SDH dysfunction with limited effect on PDH and Tau phosphorylation. The findings are suggestive that the Tyr10 on Aβ plays a critical role in the self-aggregation. Further studies are warranted to expand these findings.

Intrinsically disordered proteins in various hypotheses on the pathogenesis of Alzheimer's and Parkinson's diseases

Amyloid-β (Aβ) and α-synuclein (αS) are two intrinsically disordered proteins (IDPs) at the centers of the pathogenesis of Alzheimer's and Parkinson's diseases, respectively. Different hypotheses have been proposed for explanation of the molecular mechanisms of the pathogenesis of these two diseases, with these two IDPs being involved in many of these hypotheses. Currently, we do not know, which of these hypothesis is more accurate. Experiments face challenges due to the rapid conformational changes, fast aggregation processes, solvent and paramagnetic effects in studying these two IDPs in detail. Furthermore, pathological modifications impact their structures and energetics. Theoretical studies using computational chemistry and computational biology have been utilized to understand the structures and energetics of Aβ and αS. In this chapter, we introduce Aβ and αS in light of various hypotheses, and discuss different experimental and theoretical techniques that are used to study these two proteins along with their weaknesses and strengths. We suggest that a promising solution for studying Aβ and αS at the center of varying hypotheses could be provided by developing new techniques that link quantum mechanics, statistical mechanics, thermodynamics, bioinformatics to machine learning. Such new developments could also lead to development in experimental techniques.

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.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

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.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.