Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers

Essential collective dynamics analysis reveals nonlocal interactions of alpha-synuclein38-95 monomers with fibrillar seeds

Alpha-synuclein (α-syn) is a cytoplasmic protein commonly found in the nervous system. In solution, α-syn adopts disordered unfolded conformations, although it can also form α-helices upon binding to lipid membranes. Under conditions that are not yet fully understood, α-syn can misfold and aggregate, giving rise to β-sheet rich amyloid fibrils, which then tend to accumulate in degenerating neurons. This leads to Parkinson's disease (PD) and several other conditions collectively termed synucleinopathies. Development of disease-modifying treatments requires detailed understanding of structure and dynamics of α-syn's misfolded aggregates. We have employed 1000 ns long, all-atom molecular dynamics simulations to investigate the interaction of monomeric α-syn_38-95 fragments, which contain the most important amyloidogenic regions, with preformed fibrillar seeds composed of staggered, β-sheet rich α-syn chains of matching length. The simulations indicate that α-syn_38-95 monomers tend to form aggregates with the fibrillar seeds, although we have not observed alignment of the monomeric chains with β-strands of the fibril. To analyze the stability of these aggregates, we have employed the essential collective dynamics method, which allows making accurate assessment of dynamical coupling across individual atoms in macromolecules and supramolecular complexes. The analysis revealed extensive dynamical coupling across initially monomeric α-syn chains and the fibrillar seeds including distal regions thereof that did not contact the monomer directly. We have discussed structural origins of these long-range interactions, their impacts for the stability of α-syn aggregates, and potential implications for the development of anti-PD treatments.

Correlated conformational dynamics of the human GluN1-GluN2A type N-methyl-D-aspartate (NMDA) receptor

N-Methyl-D-aspartate receptors (NMDARs) are glutamate-gated ion channels found in the nerve cell membranes. As a result of overexcitation of NMDARs, neuronal death occurs and may lead to diseases such as epilepsy, stroke, Alzheimer's disease, and Parkinson's disease. In this study, human GluN1- GluN2A type NMDAR structure is modeled based on the X-ray structure of the Xenopus laevis template and missing loops are added by ab-initio loop modeling. The final structure is chosen according to two different model assessment scores. To be able to observe the structural changes upon ligand binding, glycine and glutamate molecules are docked into the corresponding binding sites of the receptor. Subsequently, molecular dynamics simulations of 1.3 μs are performed for both apo and ligand-bound structures. Structural parameters, which have been considered to show functionally important changes in previous NMDAR studies, are monitored as conformational rulers to understand the dynamics of the conformational changes. Moreover, principal component analysis (PCA) is performed for the equilibrated part of the simulations. From these analyses, the differences in between apo and ligand-bound simulations can be summarized as the following: The girdle right at the beginning of the pore loop, which connects M2 and M3 helices of the ion channel, partially opens. Ligands act like an adhesive for the ligand-binding domain (LBD) by keeping the bi-lobed structure together and consequently this is reflected to the overall dynamics of the protein as an increased correlation of the LBD with especially the amino-terminal domain (ATD) of the protein.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

The neural γ2α1β2α1β2 gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models-one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2-M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.

A novel Gerstmann-Sträussler-Scheinker disease mutation defines a precursor for amyloidogenic 8 kDa PrP fragments and reveals N-terminal structural changes shared by other GSS alleles

To explore pathogenesis in a young Gerstmann-Sträussler-Scheinker Disease (GSS) patient, the corresponding mutation, an eight-residue duplication in the hydrophobic region (HR), was inserted into the wild type mouse PrP gene. Transgenic (Tg) mouse lines expressing this mutation (Tg.HRdup) developed spontaneous neurologic syndromes and brain extracts hastened disease in low-expressor Tg.HRdup mice, suggesting de novo formation of prions. While Tg.HRdup mice exhibited spongiform change, PrP aggregates and the anticipated GSS hallmark of a proteinase K (PK)-resistant 8 kDa fragment deriving from the center of PrP, the LGGLGGYV insertion also imparted alterations in PrP's unstructured N-terminus, resulting in a 16 kDa species following thermolysin exposure. This species comprises a plausible precursor to the 8 kDa PK-resistant fragment and its detection in adolescent Tg.HRdup mice suggests that an early start to accumulation could account for early disease of the index case. A 16 kDa thermolysin-resistant signature was also found in GSS patients with P102L, A117V, H187R and F198S alleles and has coordinates similar to GSS stop codon mutations. Our data suggest a novel shared pathway of GSS pathogenesis that is fundamentally distinct from that producing structural alterations in the C-terminus of PrP, as observed in other prion diseases such as Creutzfeldt-Jakob Disease and scrapie.

Prion diseases can be sporadic, infectious or genetic. The central event of all prion diseases is the structural conversion of the cellular prion protein (PrP^C) to its disease associated conformer, PrP^Sc. Gerstmann-Sträussler-Scheinker Disease (GSS) is a genetic prion disease presenting as a multi-systemic neurological syndrome. A novel mutation, an eight amino acid insertion, was discovered in a young GSS patient. We created transgenic mice expressing this mutation and found that they recapitulate key features of the disease; namely PrP deposition in the brain and a low molecular weight proteinase K (PK) resistant internal PrP fragment. While structural investigations did not reveal a gross alteration in the conformation of this mutant PrP, the insertion lying at the boundary of the globular domain causes alterations in the unstructured amino terminal portion of the protein such that it becomes resistant to digestion by the enzyme thermolysin. We demonstrate by kinetic analysis and sequential digestion that this novel thermolysin resistant species is a precursor to the pathognomonic PK resistant fragment. Analysis of samples from other GSS patients revealed this same signature, suggesting a common molecular pathway.

Understanding the dynamics of monomeric, dimeric, and tetrameric α-synuclein structures in water

Human α‐synuclein (αS) is an intrinsically disordered protein associated with Parkinson's disease. Molecular mechanisms of corruptive misfolding and aggregation of αS resulting in the disease, as well as the structure and other properties of the corresponding oligomers are not entirely understood yet, preventing the development of efficient therapies. In this study, we investigate the folding dynamics of initially unfolded hypothetical αS constructs in water using all‐atom molecular dynamics simulations. We also employ the novel essential collective dynamics method to analyze the results obtained from the simulations. Our comparative analysis of monomeric, dimeric, and tetrameric αS models reveals pronounced differences in their structure and stability, emphasizing the importance of small oligomers, particularly dimers, in the process of misfolding.

Molecular mechanisms in the activation of abscisic acid receptor PYR1

The pyrabactin resistance 1 (PYR1)/PYR1-like (PYL)/regulatory component of abscisic acid (ABA) response (RCAR) proteins comprise a well characterized family of ABA receptors. Recent investigations have revealed two subsets of these receptors that, in the absence of ABA, either form inactive homodimers (PYR1 and PYLs 1–3) or mediate basal inhibition of downstream target type 2C protein phosphatases (PP2Cs; PYLs 4–10) respectively in vitro. Addition of ABA has been shown to release the apo-homodimers yielding ABA-bound monomeric holo-receptors that can interact with PP2Cs; highlighting a competitive-interaction process. Interaction selectivity has been shown to be mediated by subtle structural variations of primary sequence and ligand binding effects. Now, the dynamical contributions of ligand binding on interaction selectivity are investigated through extensive molecular dynamics (MD) simulations of apo and holo-PYR1 in monomeric and dimeric form as well as in complex with a PP2C, homology to ABA insensitive 1 (HAB1). Robust comparative interpretations were enabled by a novel essential collective dynamics approach. In agreement with recent experimental findings, our analysis indicates that ABA-bound PYR1 should efficiently bind to HAB1. However, both ABA-bound and ABA-extracted PYR1-HAB1 constructs have demonstrated notable similarities in their dynamics, suggesting that apo-PYR1 should also be able to make a substantial interaction with PP2Cs, albeit likely with slower complex formation kinetics. Further analysis indicates that both ABA-bound and ABA-free PYR1 in complex with HAB1 exhibit a higher intra-molecular structural stability and stronger inter-molecular dynamic correlations, in comparison with either holo- or apo-PYR1 dimers, supporting a model that includes apo-PYR1 in complex with HAB1. This possibility of a conditional functional apo-PYR1-PP2C complex was validated in vitro. These findings are generally consistent with the competitive-interaction model for PYR1 but highlight dynamical contributions of the PYR1 structure in mediating interaction selectivity suggesting added degrees of complexity in the regulation of the competitive-inhibition.

Protein pyrabactin resistance 1 (PYR1) belongs to a group of PYR1-like (PYL) proteins that regulate plant development and responses to conditions of drought and salinity. Recent studies have reported characterization of their molecular structures as well as elucidation of important aspects of their function; highlighting their roles as receptors for the stress responsive phytohormone, abscisic acid (ABA). However details of the molecular mechanisms regulating their receptor signalling remain enigmatic. In this work, we use molecular dynamics simulations complemented by a sophisticated statistical-mechanical analysis to investigate structural and dynamical properties of PYR1 protein and how its interaction with ABA modifies receptor-protein complex formation. Our results provide detailed insight into how the PYR1-mediated inactivation of its downstream phosphatase target is regulated by homodimer formation and yield new hypotheses, supported by in vitro experiments, for further investigation. Ultimately, this knowledge provides insight into how plants respond to stress, with potential applications in the development of crops with improved growth characteristics and higher stress tolerance.

Behavior of solvent-exposed hydrophobic groove in the anti-apoptotic Bcl-XL protein: clues for its ability to bind diverse BH3 ligands from MD simulations

Bcl-XL is a member of Bcl-2 family of proteins involved in the regulation of intrinsic pathway of apoptosis. Its overexpression in many human cancers makes it an important target for anti-cancer drugs. Bcl-XL interacts with the BH3 domain of several pro-apoptotic Bcl-2 partners. This helical bundle protein has a pronounced hydrophobic groove which acts as a binding region for the BH3 domains. Eight independent molecular dynamics simulations of the apo/holo forms of Bcl-XL were carried out to investigate the behavior of solvent-exposed hydrophobic groove. The simulations used either a twin-range cut-off or particle mesh Ewald (PME) scheme to treat long-range interactions. Destabilization of the BH3 domain-containing helix H2 was observed in all four twin-range cut-off simulations. Most of the other major helices remained stable. The unwinding of H2 can be related to the ability of Bcl-XL to bind diverse BH3 ligands. The loss of helical character can also be linked to the formation of homo- or hetero-dimers in Bcl-2 proteins. Several experimental studies have suggested that exposure of BH3 domain is a crucial event before they form dimers. Thus unwinding of H2 seems to be functionally very important. The four PME simulations, however, revealed a stable helix H2. It is possible that the H2 unfolding might occur in PME simulations at longer time scales. Hydrophobic residues in the hydrophobic groove are involved in stable interactions among themselves. The solvent accessible surface areas of bulky hydrophobic residues in the groove are significantly buried by the loop LB connecting the helix H2 and subsequent helix. These observations help to understand how the hydrophobic patch in Bcl-XL remains stable in the solvent-exposed state. We suggest that both the destabilization of helix H2 and the conformational heterogeneity of loop LB are important factors for binding of diverse ligands in the hydrophobic groove of Bcl-XL.