Molecular Dynamics Study to Investigate the Dimeric Structure of the Full-Length α-Synuclein in Aqueous Solution

In silico study on graphene quantum dots modified with various functional groups inhibiting α‑synuclein dimerization

HYPOTHESIS

Graphene quantum dots (GQDs) with various functional groups are hypothesized to inhibit the α-synuclein (αS) dimerization, a crucial step in Parkinson's disease pathogenesis. The potential of differently functionalized GQDs is systematically explored.

EXPERIMENTS

All-atom replica-exchange molecular dynamics simulations (accumulating to 75.6 μs) in explicit water were performed to study the dimerization of the αS non-amyloid component region and the influence of GQDs modified with various functional groups. Conformation ensemble, binding behavior, and free energy analysis were conducted.

FINDINGS

All studied GQDs inhibit β-sheet and backbone hydrogen bond formation in αS dimers, leading to looser oligomeric conformations. Charged GQDs severely impede the growth of extended β-sheets by providing extra contact surface. GQD binding primarily disrupts αS inter-peptide interactions through π-π stacking, CH-π interactions, and for charged GQDs, additionally through salt-bridge and hydrogen bonding interactions. GQD-COO- showed the most optimal inhibitory effect, binding mode, and intensity, which holds promise for the development of nanomedicines targeting amyloid aggregation in neurodegenerative diseases.

Structural and dynamic insights into α-synuclein dimer conformations

Parkinson disease is associated with the aggregation of the protein α-synuclein. While α-synuclein can exist in multiple oligomeric states, the dimer has been a subject of extensive debates. Here, using an array of biophysical approaches, we demonstrate that α-synuclein in vitro exhibits primarily a monomer-dimer equilibrium in nanomolar concentrations and up to a few micromolars. We then use spatial information from hetero-isotopic cross-linking mass spectrometry experiments as restrains in discrete molecular dynamics simulations to obtain the ensemble structure of dimeric species. Out of eight structural sub-populations of dimers, we identify one that is compact, stable, abundant, and exhibits partially exposed β-sheet structures. This compact dimer is the only one where the hydroxyls of tyrosine 39 are in proximity that may promote dityrosine covalent linkage upon hydroxyl radicalization, which is implicated in α-synuclein amyloid fibrils. We propose that this α-synuclein dimer features etiological relevance to Parkinson disease.

Effect of Electric Field on α-Synuclein Fibrils: Revealed by Molecular Dynamics Simulations

The self-association of amylogenic proteins to the fibril form is considered a pivotal factor in the pathogenesis of neurodegenerative diseases, including Parkinson’s disease (PD). PD causes unintended or uncontrollable movements in its common symptoms. α-synuclein is the major cause of PD development and thus has been the main target of numerous studies to suppress and sequester its expression or effectively degrade it. Nonetheless, to date, there are no efficient and proven ways to prevent pathological protein aggregation. Recent investigations proposed applying an external electric field to interrupt the fibrils. This method is a non-invasive approach that has a certain benefit over others. We performed molecular dynamics (MD) simulations by applying an electric field on highly toxic fibrils of α-synuclein to gain a molecular-level insight into fibril disruption mechanisms. The results revealed that the applied external electric field induces substantial changes in the conformation of the α-synuclein fibrils. Furthermore, we show the threshold value for electric field strength required to completely disrupt the α-synuclein fibrils by opening the hydrophobic core of the fibril. Thus, our findings might serve as a valuable foundation to better understand molecular-level mechanisms of the α-synuclein fibrils disaggregation process under an applied external electric field.

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with α-Substituent Electrophiles Targeting Cysteine

Both reversible noncovalent inhibitors and irreversible covalent inhibitors targeting tyrosine kinases have their disadvantages. The reversible covalent inhibitors with electrophilic group cyanoacrylamide as warheads reacting with cysteine residues could solve the dilemmas. However, there are still several unresolved issues regarding the electrophilic groups. In this manuscript, a series of EGFR inhibitors with double electron-withdrawing substituents introduced into the Cα position on the olefin bond were designed and synthesized. The binding structures and characteristics of inhibitors with the kinase in both the first noncovalent binding phase and the second covalent binding step were explored and combined with molecular docking and molecular dynamics simulations. Then, the reverse β-elimination reactions of the thiol-Michael adducts were investigated by applying density functional theory calculations. In addition, the effects of different electrophilic substituents of Cα on the binding between the inhibitors and kinase were elucidated. The results suggested that the electrophilicity and size of the electron-withdrawing groups play an important role in the specific interactions during the reaction. The compounds with the electron-withdrawing groups that had medium electrostatic and steric complementarity to the kinase active site could cooperatively stabilize the complexes and showed relatively good potent activities in the kinase assay experiment. The mechanical and structural information in this study could enhance our understanding of the functioning of the electron-withdrawing groups in the covalent inhibitors. The results might help to design efficient cysteine targeting inhibitors in the future.

Concerted enhanced-sampling simulations to elucidate the helix-fibril transition pathway of intrinsically disordered α-Synuclein

Fibril formation of α-synuclein is linked with Parkinson's disease. The intrinsically disordered nature of α-syn provides extensive conformational plasticity and becomes difficult to characterize its transition pathway from native monomeric to disease-associated fibril form. We implemented different simulation methods such as steered dynamics-umbrella sampling, and replica-exchange and conventional MD simulations to access various conformational states of α-syn. Nineteen distinct intermediate structures were identified by free energy landscape and cluster analysis. They were then sorted based on secondary structure and solvent exposure of fibril-core residues to illustrate the fibril dissociation pathway. The analysis showed that following the initial dissociation of the polypeptide chain from the fibril, α-syn might form either compact-conformations by long-range interactions or extended-conformations stabilized by local interactions. This leads α-syn to adapt two different pathways. The secondary structure, solvation, contact distance, interaction energies and backbone dihedrals of thirty-two selected residues were analyzed for all the 19 intermediates. The results suggested that formation of β-turns, reorganization of salt bridges, and dihedral changes in the hydrophobic regions are the major driving forces for helix-fibril transition. Structural features of the intermediates also correlated with the earlier experimental and computational studies. The study provides critical information on the fibrillation pathway of α-syn.

Engineering of a protein probe with multiple inputs and multiple outputs for evaluation of alpha synuclein aggregation states

The aggregation of α-synuclein (αS) into oligomers and fibrils is implicated in the pathology of Parkinson's Disease (PD). While a molecular probe for rapid and comprehensive evaluation of αS aggregation states is critical for a better understanding of PD pathology, identification of therapeutic candidates, and the development of early diagnostic strategies, no such probe has yet to be developed. A structurally flexible αS variant, PG65, was previously developed as a target binding-driven, conformation-switching molecular probe for rapid αS oligomer detection. Though informative, detection using PG65 provides no comprehensive assessment of the αS aggregation states. In the present study, we report engineering of a molecular probe, PG65-MIMO (a PG65 variant with Multiple-Inputs and Multiple-Outputs), that rapidly (within 2 hr) produces comprehensive information on αS aggregation states. PG65-MIMO generates distinct fluorescence responses to the three major αS conformers (monomers, oligomers, and fibrils). PG65-MIMO also displays unique fluorescent signals for αS oligomers, depending on the tris(2-carboxyethyl)phosphine (TCEP) concentration. Our results suggest that the TCEP dependent signaling of PG65-MIMO may be associated with its conformational states. Overall, our study illustrates engineering of an αS variant to create a molecular probe for handling multiple inputs and multiple outputs, addressing the technological gap in αS detection.

Dimerization of α-Synuclein Fragments Studied by Isothermal-Isobaric Replica-Permutation Molecular Dynamics Simulation

Aggregates and fibrils of intrinsically disordered α-synuclein are associated with Parkinson's disease. Within a non-amyloid β component (NAC) spanning from the 61st to the 95th residue of α-synuclein, an 11-residue segment called NACore (₆₈GAVVTGVTAVA₇₈) is an essential region for both fibril formation and cytotoxicity. Although NACore peptides alone are known to form aggregates and amyloid fibrils, the mechanisms of aggregation and fibrillation remain unknown. This study investigated the dimerization process of NACore peptides as the initial stage of the aggregation and fibrillation processes. We performed an isothermal-isobaric replica-permutation molecular dynamics simulation, which is one of the efficient sampling methods, for the two NACore peptides in explicit water over 96 μs. The simulation succeeded in sampling a variety of dimer structures. An analysis of secondary structure revealed that most of the NACore dimers form intermolecular β-bridges. In particular, more antiparallel β-bridges were observed than parallel β-bridges. We also found that intramolecular secondary structures such as α-helix and antiparallel β-bridge are stabilized in the pre-dimer state. However, we identified that the intermolecular β-bridges tend to form directly between residues with no specific structure rather than via the intramolecular β-bridges. This is because the NACore peptides still have a low propensity to form the intramolecular secondary structures even though they are stabilized in the pre-dimer state.

The role of water in the primary nucleation of protein amyloid aggregation

The understanding of the complex conformational landscape of amyloid aggregation and its modulation by relevant physicochemical and cellular factors is a prerequisite for elucidating some of the molecular basis of pathology in amyloid related diseases, and for developing and evaluating effective disease-specific therapeutics to reduce or eliminate the underlying sources of toxicity in these diseases. Interactions of proteins with solvating water have been long considered to be fundamental in mediating their function and folding; however, the relevance of water in the process of protein amyloid aggregation has been largely overlooked. Here, we provide a perspective on the role water plays in triggering primary amyloid nucleation of intrinsically disordered proteins (IDPs) based on recent experimental evidences. The initiation of amyloid aggregation likely results from the synergistic effect between both protein intermolecular interactions and the properties of the water hydration layer of the protein surface. While the self-assembly of both hydrophobic and hydrophilic IDPs would be thermodynamically favoured due to large water entropy contributions, large desolvation energy barriers are expected, particularly for the nucleation of hydrophilic IDPs. Under highly hydrating conditions, primary nucleation is slow, being facilitated by the presence of nucleation-active surfaces (heterogeneous nucleation). Under conditions of poor water activity, such as those found in the interior of protein droplets generated by liquid-liquid phase separation, however, the desolvation energy barrier is significantly reduced, and nucleation can occur very rapidly in the bulk of the solution (homogeneous nucleation), giving rise to structurally distinct amyloid polymorphs. Water, therefore, plays a key role in modulating the transition free energy of amyloid nucleation, thus governing the initiation of the process, and dictating the type of preferred primary nucleation and the type of amyloid polymorph generated, which could vary depending on the particular microenvironment that the protein molecules encounter in the cell.

α-Synuclein in traumatic and vascular diseases of the central nervous system

α-Synuclein (α-Syn) is a small, soluble, disordered protein that is widely expressed in the nervous system. Although its physiological functions are not yet fully understood, it is mainly involved in synaptic vesicle transport, neurotransmitter synthesis and release, cell membrane homeostasis, lipid synthesis, mitochondrial and lysosomal activities, and heavy metal removal. The complex and inconsistent pathological manifestations of α-Syn are attributed to its structural instability, mutational complexity, misfolding, and diverse posttranslational modifications. These effects trigger mitochondrial dysfunction, oxidative stress, and neuroinflammatory responses, resulting in neuronal death and neurodegeneration. Several recent studies have discovered the pathogenic roles of α-Syn in traumatic and vascular central nervous system diseases, such as traumatic spinal cord injury, brain injury, and stroke, and in aggravating the processes of neurodegeneration. This review aims to highlight the structural and pathophysiological changes in α-Syn and its mechanism of action in traumatic and vascular diseases of the central nervous system.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

Effects of natural compounds on conformational properties and hairpin formation of amyloid-β42 monomer: docking and molecular dynamics simulation study

The β₄₂ amyloid peptides (Aβ) are identified as a candidate target for Alzheimer's drugs. Phenolic compounds can bind to the Aβ and inhibit amyloid formation. However, the inhibitory mechanism of phenolic compounds remains unclear. In this study, the molecular dynamic simulation and docking program were used to characterize the molecular details of inhibitory mechanism of the phenolic compounds. Our Results show that the phenolic compounds can bind to hydrophobic region in Aβ₄₂ monomer and alter hydrophobic interactions network at Aβ₄₂ which play a key role in β-sheet formation. The cluster analysis and interactions network analysis were used to probe conformational changes in Aβ₄₂. In most populated clusters of Aβ₄₂-phenolic compounds complexes, the sheet structures were not observed or reduced. It seems that the binding of phenolic compounds can induce larger conformational diversity for amyloid peptide and changes conformational properties of amyloid peptide. The phenolic compounds can deform β-Hairpin structure of Aβ by destabilizing salt bridges E22-K28 and D23-K28 which can alter the conformation of Aβ₄₂ in aqueous solution. These findings are in accordance with experimental results, to some extent give a molecular level interpretation for the inhibitory mechanism of phenolic compounds .In addition, this study may add important new details to the inhibitory mechanism of Alzheimer's drug.Communicated by Ramaswamy H. Sarma.

Probing the Basis of α-Synuclein Aggregation by Comparing Simulations to Single-Molecule Experiments

Intrinsically disordered proteins often play an important role in protein aggregation. However, it is challenging to determine the structures and interactions that drive the early stages of aggregation because they are transient and obscured in a heterogeneous mixture of disordered states. Even computational methods are limited because the lack of ordered structure makes it difficult to ensure that the relevant conformations are sampled. We address these challenges by integrating atomistic simulations with high-resolution single-molecule measurements reported previously, using the measurements to help discern which parts of the disordered ensemble of structures in the simulations are most probable while using the simulations to identify residues and interactions that are important for oligomer stability. This approach was applied to α-synuclein, an intrinsically disordered protein that aggregates in the context of Parkinson's disease. We simulated single-molecule pulling experiments on dimers, the minimal oligomer, and compared them to force spectroscopy measurements. Force-extension curves were simulated starting from a set of 66 structures with substantial structured content selected from the ensemble of dimer structures generated at zero force via Monte Carlo simulations. The pattern of contour length changes as the structures unfolded through intermediate states was compared to the results from optical trapping measurements on the same dimer to discern likely structures occurring in the measurements. Simulated pulling curves were generally consistent with experimental data but with a larger number of transient intermediates. We identified an ensemble of β-rich dimer structures consistent with the experimental data from which dimer interfaces could be deduced. These results suggest specific druggable targets in the structural motifs of α-synuclein that may help prevent the earliest steps of oligomerization.

Sulfonated Compounds Bind with Prostatic Acid Phosphatase (PAP248-286) to Inhibit the Formation of Amyloid Fibrils

The peptide segment of prostatic acid phosphatase (PAP248-286) aggregates to form SEVI (semen-derived enhancer of virus infection) amyloid fibrils. These are characteristic seminal amyloids that have the ability to promote the effect of HIV infection. In this paper, we explore the binding of sulfonated compounds with PAP248-286 through an in silico study. Three derivatives of suramin, NF110, NF279, and NF340, are selected. All of these sulfonated molecules bind to PAP248-286 and alter the conformation of the peptide, even though they have various structures, sizes, and configurations. The compounds bind with PAP248-286 through multiple interactions, such as hydrogen-bonding interactions, hydrophobic interactions, π-π stacking interactions, and electrostatic interactions. However, NF110, which has an X-shaped configuration, has the highest binding affinity of the three derivatives investigated. We also perform surface plasmon resonance and a Congo red assay to validate the results. The interactions between PAP248-286 and the sulfonated compounds are proposed to depend on the orientations of the sulfonate groups and the specific configurations of the compounds instead of the number of sulfonate groups. NF110 molecules occupy the exposed binding sites of PAP248-286, blocking interactions between the peptides. Therefore, these compounds are important in inhibiting the aggregation of PAP248-286. Herein, we provide useful information to develop new efficient microbicides to antagonize seminal amyloid fibrils and to block HIV transmission.