Replica Exchange Molecular Dynamics Simulations Research Articles

Melatonin Inhibits hIAPP Oligomerization by Preventing β-Sheet and Hydrogen Bond Formation of the Amyloidogenic Region Revealed by Replica-Exchange Molecular Dynamics Simulation.

The pathogenesis of type 2 diabetes (T2D) is highly related to the abnormal self-assembly of the human islet amyloid polypeptide (hIAPP) into amyloid aggregates. To inhibit hIAPP aggregation is considered a promising therapeutic strategy for T2D treatment. Melatonin (Mel) was reported to effectively impede the accumulation of hIAPP aggregates and dissolve preformed fibrils. However, the underlying mechanism at the atomic level remains elusive. Here, we performed replica-exchange molecular dynamics (REMD) simulations to investigate the inhibitory effect of Mel on hIAPP oligomerization by using hIAPP20–29 octamer as templates. The conformational ensemble shows that Mel molecules can significantly prevent the β-sheet and backbone hydrogen bond formation of hIAPP20–29 octamer and remodel hIAPP oligomers and transform them into less compact conformations with more disordered contents. The interaction analysis shows that the binding behavior of Mel is dominated by hydrogen bonding with a peptide backbone and strengthened by aromatic stacking and CH–π interactions with peptide sidechains. The strong hIAPP–Mel interaction disrupts the hIAPP20–29 association, which is supposed to inhibit amyloid aggregation and cytotoxicity. We also performed conventional MD simulations to investigate the influence and binding affinity of Mel on the preformed hIAPP1–37 fibrillar octamer. Mel was found to preferentially bind to the amyloidogenic region hIAPP20–29, whereas it has a slight influence on the structural stability of the preformed fibrils. Our findings illustrate a possible pathway by which Mel alleviates diabetes symptoms from the perspective of Mel inhibiting amyloid deposits. This work reveals the inhibitory mechanism of Mel against hIAPP20–29 oligomerization, which provides useful clues for the development of efficient anti-amyloid agents.

Fundamental Redesign of the TIGER2hs Kernel to Address Severe Parameter Sensitivity.

Replica exchange molecular dynamics simulations are one of the most popular approaches to enhance conformational sampling of molecular systems. Applications range from protein folding to protein-protein or other host-guest interactions, as well as binding free energy calculations. While these methods are computationally expensive, highly accurate results can be obtained. We recently developed TIGER2hs, an improved version of the temperature intervals with global exchange of replicas (TIGER2) algorithm. This method combines the replica-based enhanced sampling in an explicit solvent with a hybrid solvent energy evaluation. During the exchange attempts, bulk water is replaced by an implicit solvent model, allowing sampling with significantly less replicas than parallel tempering (REMD). This enables accurate enhanced sampling calculations with only a fraction of computational resources compared to REMD. Our latest results highlight several issues with sampling imbalance and parameter sensitivity within the original TIGER2 exchange algorithms that affect the overall state populations. A high sensitivity on replica number and maximum temperature is eliminated by changing to a pairwise exchange kernel (PE) without additional sorting. Simulations are controlled by adjusting the average temperature change per exchange ⟨ΔT/χ⟩ to below 30 K to mimic a controlled temperature mixing of replicas similar to REMD. Thus, this parameter provides an applicable property for selecting combinations of replica number and maximum temperature to adjust simulations for best accuracy, with flexible resource investment. This increases the robustness of the method and ensures results in excellent agreement with REMD, as demonstrated for three different peptides.

Tachykinin Neuropeptides and Amyloid β (25-35) Assembly: Friend or Foe?

Amyloid β (Aβ) protein is responsible for Alzheimer's disease, and one of its important fragments, Aβ(25-35), is found in the brain and has been shown to be neurotoxic. Tachykinin neuropeptides, including Neuromedin K (NK), Kassinin, and Substance P, have been reported to reduce Aβ(25-35)'s toxicity in cells even though they share similar primary structures with Aβ(25-35). Here, we seek to understand the molecular mechanisms of how these peptides interact with Aβ(25-35) and to shed light on why some peptides with similar primary structures are toxic and others nontoxic. We use both experimental and computational methods, including ion mobility mass spectrometry and enhanced-sampling replica-exchange molecular dynamics simulations, to study the aggregation pathways of Aβ(25-35), NK, Kassinin, Substance P, and mixtures of the latter three with Aβ(25-35). NK and Substance P were observed to remove the higher-order oligomers (i.e., hexamers and dodecamers) of Aβ(25-35), which are related to its toxicity, although Substance P did so more slowly. In contrast, Kassinin was found to promote the formation of these higher-order oligomers. This result conflicts with what is expected and is elaborated on in the text. We also observe that even though they have significant structural homology with Aβ(25-35), NK, Kassinin, and Substance P do not form hexamers with a β-sheet structure like Aβ(25-35). The hexamer structure of Aβ(25-35) has been identified as a cylindrin, and this structure has been strongly correlated to toxic species. The reasons why the three tachykinin peptides behave so differently when mixed with Aβ(25-35) are discussed.

Investigating the stabilisation of IFN-α2a by replica exchange molecular dynamics simulation.

Current biopharmaceuticaldrugs are mainlya class of peptides or proteins that play an essential role in the treatment of many diseases. Such peptides/proteins are usually thermally unstable and may lose their bioactivity when exposed to ambient conditions. Therefore, they are not suitable for long-term storage. Lyophilisation is the most common method to prolong shelf life of solid peptide/protein drugs; however, the freeze-drying process can lead to irreversible damage. In the present study, human interferon-alpha2a (IFN-α2a) was selected as a model protein drug; four disaccharides (β-lactose, β-maltose, sucrose, and trehalose) were selected as bioactive protectants. We investigated the effects of different protectants on IFN-α2a under various ambient conditions (vacuum, dry state, and aqueous solution) using replica exchange molecular dynamics simulation. The protective effect of β-maltose on IFN-α2a was the highest inaqueous solution and dry state, β-lactose showed a poor protective effect in all three conditions, the performance of sucrose was good in all conditions, and trehalose showed a better protective effect under vacuum conditions and in aqueous solution. Disaccharides form H-bonds with water, thereby preventing water from the tertiary structure of proteins. Trehalose forms strong H-bonds with water which explains its extraordinary stability.

Quantitative Agreement between Conformational Substates of Holo Calcium-Loaded Calmodulin Detected by Double Electron–Electron Resonance EPR and Predicted by Molecular Dynamics Simulations

Calcium-loaded calmodulin (CaM/4Ca2+) comprises two domains that undergo rigid body reorientation from a predominantly extended conformation to a compact one upon binding target peptides. A recent replica-exchange molecular dynamics (MD) simulation on holo CaM/4Ca2+ suggested the existence of distinct structural clusters (substates) along the path from extended to compact conformers in the absence of substrates. Here, we experimentally demonstrate the existence of CaM/4Ca2+ substates trapped in local minima by three freezing/annealing regimes (slow, 40 s; intermediate, 1.5 s; fast, 0.5 ms) using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure interdomain distances between nitroxide spin-labels positioned at A17C and A128C in the N- and C-terminal domains, respectively. The DEER echo curves were directly fit to population-optimized P(r) pairwise distance distributions calculated from the coordinates of the MD clusters and compact crystal structure. DEER data on fully deuterated CaM/4Ca2+ were acquired at multiple values of the second echo period (10-35 μs) and analyzed globally to eliminate instrumental and overfitting artifacts and ensure accurate populations, peak positions, and widths. The DEER data for all three freezing regimes are quantitatively accounted for within experimental error by 5-6 distinct conformers comprising a predominantly populated extended form (60-75%) and progressively more compact states whose populations decrease as the degree of compactness increases. The shortest interdomain separation is found in the compact crystal structure, which has an occupancy of 4-6%. Thus, CaM/4Ca2+ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

The Amyloidogenic Peptide Amyloid Beta(16–22) Displays Facet Dependent Conformation on Metal Surfaces

Currently, it is not understood how metal nanoparticles influence the formation of protein fibrils, although recent literature highlights that the shape and chemical composition of such nanoparticles can strongly influence the process. Understanding this process at a fundamental level can potentially unlock routes to the development of new therapeutics, as well as novel materials for technological applications. This requires a microscopic picture of the behaviour of amyloidogenic proteins on metal surfaces. Using replica exchange molecular dynamics simulations, we investigate the conformation of the model amyloidogenic peptide, Aβ(16–22), on different gold and silver surfaces. The conformation of the peptide on gold surfaces also shows a strong facet dependence, with fibril-like conformations being promoted in the 100 surface and inhibited on the 111 surface. A smaller degree of facet dependence is seen for silver with the peptide behaving similar on both of these. The difference in the facet dependence can be related to the difference between direct adsorption onto the gold 111 surface, with a preference towards indirect (water mediated) adsorption onto the other surfaces. This new information on the behaviour of an amyloidogenic peptide on metal surfaces can give insight into the size-dependent effect of nanoparticles on fibril formation and the use of surfaces to control fibrillation.

Distribution and Structure Analysis of Fibril-Forming Peptides Focusing on Concentration Dependency.

We focus on the concentration dependency of fibril-forming peptides, which have the potential of aggregation by themselves. In this study, we performed replica-exchange molecular dynamics simulations of Lys-Phe-Phe-Glu (KFFE) fragments, which are known to form fibrils in experiments under different concentration environments. The analysis by static structure factors suggested that the density fluctuation of the KFFE fragments becomes large as the concentration increases. It was also found that the number of β-structures and oligomers also increases under a high concentration environment. Hence, a high concentration environment of fibril-forming peptides is likely to cause protein aggregation.

Effect of Cholesterol Molecules on Aβ1-42 Wild-Type and Mutants Trimers.

Alzheimer’s disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain, and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with a solute tempering simulation of an Aβ42 trimer with cholesterol and compare it with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex, involving hydrophobic residues L17–F20 and L30–M35 with a non-negligible contribution of loop residues D22–K28 and N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations moderately impact cholesterol-binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into aggregates but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates; this provides a mechanistic link between cholesterol and Alzheimer’s disease.

Influence of Aqueous Arginine Solution on Regulating Conformational Stability and Hydration Properties of the Secondary Structural Segments of a Protein at Elevated Temperatures: A Molecular Dynamics Study.

The effects of aqueous arginine solution on the conformational stability of the secondary structural segments of a globular protein, ubiquitin, and the structure and dynamics of the surrounding water and arginine were examined by performing atomistic molecular dynamics (MD) simulations. Attempts have been made to identify the osmolytic efficacy of arginine solution, and its influence in guiding the hydration properties of the protein at an elevated temperature of 450 K. The similar properties of the protein in pure water at elevated temperatures were computed and compared. Replica exchange MD simulation was performed to explore the arginine solution's sensitivity in stabilizing the protein conformations for a wide range of temperatures (300-450 K). It was observed that although all the helices and strands of the protein undergo unfolding at elevated temperature in pure water, they exhibited native-like conformational dynamics in the presence of arginine at both ambient and elevated temperatures. We find that the higher free energy barrier between the folded native and unfolded states of the protein primarily arises from the structural transformation of α-helix, relative to the strands. Our study revealed that the water structure around the secondary segments depends on the nature of amino acid compositions of the helices and strands. The reorientation of water dipoles around the helices and strands was found hindered due to the presence of arginine in the solution; such hindrance reduces the possibility of exchange of hydrogen bonds that formed between the secondary segments of protein and water (PW), and as a result, PW hydrogen bonds take longer time to relax than in pure water. On the other hand, the origin of slow relaxation of protein-arginine (PA) hydrogen bonds was identified to be due to the presence of different types of protein-bound arginine molecules, where arginine interacts with the secondary structural segments of the protein through multiple/bifurcated hydrogen bonds. These protein-bound arginine formed different kinds of bridged PA hydrogen bonds between amino acid residues of the same secondary segments or among multiple bonds and helped protein to conserve its native folded form firmly.

Counterintuitive Electrostatics upon Metal Ion Coordination to a Receptor with Two Homotopic Binding Sites.

The consecutive binding of two potassium ions to a bis(18-crown-6) analogue of Tröger’s base (BCETB) in water was studied by isothermal titration calorimetry using four different salts, KCl, KI, KSCN, and K2SO4. A counterintuitive result was observed: the enthalpy change associated with the binding of the second ion is more negative than that of the first (ΔHbind,2° < ΔHbind,1°). This remarkable finding is supported by continuum electrostatic theory as well as by atomic scale replica exchange molecular dynamics simulations, where the latter robustly reproduces experimental trends for all simulated salts, KCl, KI, and KSCN, using multiple force fields. While an enthalpic K+–K+attraction in water poses a small, but fundamentally important, contribution to the overall interaction, the probability of the collapsed conformation (COL) of BCETB, where both crown ether moieties (CEs) of BCETB are bent in toward the cavity, was found to increase successively upon binding of the first and second potassium ions. The promotion of the COL conformation reveals favorable intrinsic interactions between the potassium coordinated CEs, which further contribute to the observation that ΔHbind,2° < ΔHbind,1°. While the observed trend is independent of the counterion, the origin of the significantly larger magnitude of the difference ΔHbind,2° – ΔHbind,1° observed experimentally for KSCN was studied in light of the weaker hydration of the thiocyanate anion, resulting in an enrichment of thiocyanate ions close to BCETB compared to the other studied counterions.

Characterization of Heparin's Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy.

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin's solution conformational dynamics and structure-function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule's underlying conformational ensemble, which then provides important information necessary for understanding structure-function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin's conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. The anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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.

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.

Unraveling the Influence of K280 Acetylation on the Conformational Features of Tau Core Fragment: A Molecular Dynamics Simulation Study.

Abnormal aggregation of the microtubule-associated protein Tau is closely associated with tauopathies, including Alzheimer’s disease and chronic traumatic encephalopathy. The hexapeptide 275VQIINK280 (PHF6*), a fibril-nucleating core motif of Tau, has been shown to play a vital role in the aggregation of Tau. Mounting experiment evidence demonstrated the acetylation of a single-lysine residue K280 in the PHF6* was a critical event for the formation of pathological Tau amyloid deposits. However, the underlying mechanisms by which K280 acetylation affects Tau aggregation at the atomic level remain elusive. In this work, we performed replica exchange molecular dynamics simulations to investigate the influence of acetylation of K280 on the aggregation of PHF6*. Our simulations show that acetylation of K280 not only enhances the self-assembly capability of PHF6* peptides but also increases the β-sheet structure propensity of the PHF6*. The inter-molecular interactions among PHF6* peptides are strengthened by the acetylation of K280, resulting in an increased ordered β-sheet-rich conformations of the PHF6* assemblies along with a decrease of the structural diversity. The residue-pairwise contact frequency analysis shows that K280 acetylation increases the interactions among the hydrophobic chemical groups from PHF6* peptides, which promotes the aggregation of PHF6*. This study offers mechanistic insights into the effects of acetylation on the aggregation of PHF6*, which will be helpful for an in-depth understanding of the relationship between acetylation and Tau aggregation at the molecular level.

Interaction of KRSR Peptide with Titanium Dioxide Anatase (100) Surface: A Molecular Dynamics Simulation Study.

Due to its tensile strength and excellent biocompatibility, titanium (Ti) is commonly used as an implant material in medicine and dentistry. The success of dental implants depends on the formation of a contact between the oxidized surface of Ti implant and the surrounding bone tissue. The adsorption of proteins and peptides to the implant surface allows the bone-forming osteoblast cells to adhere to such modified surfaces. Recently, it has been observed that tetrapeptide KRSR (Lys-Arg-Ser-Arg) functionalization could promote osteoblast adhesion to implant surfaces. This may facilitate the establishment of an efficient bone-to implant contact and improve implant stability during the healing process. GROMACS, a molecular dynamics software package was used to perform a 200 ns simulation of adsorption of the KRSR peptide to the TiO2 (anatase) surface in an aqueous environment. The molecule conformations were mapped with Replica Exchange Molecular Dynamics (REMD) simulations to assess the possible peptide conformations on the anatase surface, and the umbrella sampling method was used to calculate the binding energy of the most common conformation. The simulations have shown that the KRSR peptide migrates and attaches to the surface in a stable position. The dominant amino acid residue interacting with the TiO2 surface was the N-terminal charged lysine (K) residue. REMD indicated that there is a distinct conformation that is taken by the KRSR peptide. In this conformation the surface interacts only with the lysine residue while the ser (S) and arg (R) residues interact with water molecules farther from the surface. The binding free energy of the most common conformation of KRSR peptide to the anatase (100) surface was ΔG = −8.817 kcal/mol. Our result suggests that the N-terminal lysine residue plays an important role in the adhesion of KRSR to the TiO2 surface and may influence the osseointegration of dental implants.

Replica exchange molecular dynamics simulations reveal self-association sites in M-crystallin caused by mutations provide insights of cataract

Crystallins are ubiquitous, however, prevalence is seen in eye lens. Eye lens crystallins are long-lived and structural intactness is required for maintaining lens transparency and protein solubility. Mutations in crystallins often lead to cataract. In this study, we performed mutations at specific sites of M-crystallin, a close homologue of eye lens crystallin and studied by using replica exchange molecular dynamics simulation with generalized Born implicit solvent model. Mutations were made on the Ca2+ binding residues (K34D and S77D) and in the hydrophobic core (W45R) which is known to cause congenital cataract in homologous γD-crystallin. The chosen mutations caused large motion of the N-terminal Greek key, concomitantly broke the interlocking Greek keys interactions and perturbed the compact core resulting in several folded and partially unfolded states. Partially unfolded states exposed large hydrophobic patches that could act as precursors for self-aggregation. Accumulation of such aggregates is the potential cause of cataract in homologous eye lens crystallins.

Machine learning-based prediction of drug and ligand binding in BCL-2 variants through molecular dynamics

Venetoclax is a BH3 (BCL-2 Homology 3) mimetic used to treat leukemia and lymphoma by inhibiting the anti-apoptotic BCL-2 protein thereby promoting apoptosis of cancerous cells. Acquired resistance to Venetoclax via specific variants in BCL-2 is a major problem for the successful treatment of cancer patients. Replica exchange molecular dynamics (REMD) simulations combined with machine learning were used to define the average structure of variants in aqueous solution to predict changes in drug and ligand binding in BCL-2 variants. The variant structures all show shifts in residue positions that occlude the binding groove, and these are the primary contributors to drug resistance. Correspondingly, we established a method that can predict the severity of a variant as measured by the inhibitory constant (Ki) of Venetoclax by measuring the structure deviations to the binding cleft. In addition, we also applied machine learning to the phi and psi angles of the amino acid backbone to the ensemble of conformations that demonstrated a generalizable method for drug resistant predictions of BCL-2 proteins that elucidates changes where detailed understanding of the structure-function relationship is less clear.

Integrative structural dynamics probing of the conformational heterogeneity in synaptosomal-associated protein 25

SUMMARYSNAP-25 (synaptosomal-associated protein of 25 kDa) is a prototypical intrinsically disordered protein (IDP) that is unstructured by itself but forms coiled-coil helices in the SNARE complex. With high conformational heterogeneity, detailed structural dynamics of unbound SNAP-25 remain elusive. Here, we report an integrative method to probe the structural dynamics of SNAP-25 by combining replica-exchange discrete molecular dynamics (rxDMD) simulations and label-based experiments at ensemble and single-molecule levels. The rxDMD simulations systematically characterize the coil-to-molten globular transition and reconstruct structural ensemble consistent with prior ensemble experiments. Label-based experiments using Förster resonance energy transfer and double electron-electron resonance further probe the conformational dynamics of SNAP-25. Agreements between simulations and experiments under both ensemble and single-molecule conditions allow us to assign specific helix-coil transitions in SNAP-25 that occur in submillisecond timescales and potentially play a vital role in forming the SNARE complex. We expect that this integrative approach may help further our understanding of IDPs.

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure.

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1-1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Moreover, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

A Comprehensive Insight into the Mechanisms of Dopamine in Disrupting Aβ Protofibrils and Inhibiting Aβ Aggregation.

Fibrillary aggregates of amyloid-β (Aβ) are the pathological hallmark of Alzheimer's disease (AD). Clearing Aβ deposition or inhibiting Aβ aggregation is a promising approach to treat AD. Experimental studies reported that dopamine (DA), an important neurotransmitter, can inhibit Aβ aggregation and disrupt Aβ fibrils in a dose-dependent manner. However, the underlying molecular mechanisms still remain mostly elusive. Herein, we investigated the effect of DA on Aβ42 protofibrils at three different DA-to-Aβ molar ratios (1:1, 2:1, and 10:1) using all-atom molecular dynamics simulations. Our simulations demonstrate that protonated DA at a DA-to-Aβ ratio of 2:1 exhibits stronger Aβ protofibril disruptive capacity than that at a molar-ratio of 1:1 by mostly disrupting the F4-L34-V36 hydrophobic core. When the ratio of DA-to-Aβ increases to 10:1, DA has a high probability to bind to the outer surface of protofibril and has negligible effect on the protofibril structure. Interestingly, at the same DA-to-Aβ ratio (10:1), a mixture of protonated (DA+) and deprotonated (DA0) DA molecules significantly disrupts Aβ protofibrils by the binding of DA0 to the F4-L34-V36 hydrophobic core. Replica-exchange molecular dynamics simulations of Aβ42 dimer show that DA+ inhibits the formation of β-sheets, K28-A42/K28-D23 salt-bridges, and interpeptide hydrophobic interactions and results in disordered coil-rich Aβ dimers, which would inhibit the subsequent fibrillization of Aβ. Further analyses reveal that DA disrupts Aβ protofibril and prevents Aβ dimerization mostly through π-π stacking interactions with residues F4, H6, and H13, hydrogen bonding interactions with negatively charged residues D7, E11, E22 and D23, and cation-π interactions with residues R5. This study provides a complete picture of the molecular mechanisms of DA in disrupting Aβ protofibril and inhibiting Aβ aggregation, which could be helpful for the design of potent drug candidates for the treatment/intervention of AD.