Characterization of the AXH domain of Ataxin-1 using enhanced sampling and functional mode analysis

Governing dynamics and preferential binding of the AXH domain influence the aggregation pathway of Ataxin-1

The present state of understanding the mechanism of Spinocerebellar Ataxia-1, a fatal neurodegenerative disease linked to the protein Ataxin-1 (ATXN1), is baffled by a set of self-contradictory, and hence, inconclusive observations. This fallacy poses a bottleneck to the effective designing of curable drugs as the field is currently missing the specific druggable site. To understand the fundamentals of pathogenesis, we tried to decipher the intricacies of the extremely complicated landscape by targeting the relevant species that supposedly dictate the structure-function paradigm. The atomic-level description and characterization of the dynamism of the systems reveal the existence of structural polymorphism in all the leading stakeholders of the overall system. The very existence of conformational heterogeneity in every species creates numerous possible combinations of favorable interactions because of the variability in segmental cross-talks and hence claims its role in the choice of routes between functional activity and dysfunctional disease-causing aggregation. Despite this emergent configurational diversity, there is a common mode of operative intermolecular forces that dictates the extent of stability of all the multimeric complexes due to the localized population of a specific type of residue. The present research proposes a dynamic switch mechanism between aggregability and functional activity, based on the logical interpretation of the estimated variables, which is practically dictated by the effective concentration of the interacting species involved in the cell.

Toward the design and development of peptidomimetic inhibitors of the Ataxin-1 aggregation pathway

Spinocerebellar ataxia type 1 is a degenerative disorder caused by polyglutamine expansions and aggregation of Ataxin-1. The interaction between Capicua (CIC) and the AXH domain of Ataxin-1 protein has been suggested as a possible driver of aggregation for the expanded Ataxin-1 protein and the subsequent onset of spinocerebellar ataxia 1. Experimental studies have demonstrated that short constructs of CIC may prevent such aggregation and suggested this as a possible candidate to inspire the rational design of peptidomimetics. In this work, molecular modeling techniques, namely the alchemical mutation and force field-based molecular dynamics, have been employed to propose a pipeline for the rational design of a CIC-inspired inhibitor of the ataxin-1 aggregation pathway. In particular, this study has shown that the alchemical mutation can estimate the affinity between AXH and CIC with good correlation with experimental data, while molecular dynamics shed light on molecular mechanisms that occur for stabilization of the interaction between the CIC-inspired construct and the AXH domain of Ataxin-1. This work lays the foundation for a rational methodology for the in silico screening and design of peptidomimetics, which can expedite and streamline experimental studies to identify strategies for inhibiting the ataxin-1 aggregation pathway.

In silico investigation of Alsin RLD conformational dynamics and phosphoinositides binding mechanism

Alsin is a protein known for its major role in neuronal homeostasis and whose mutation is associated with early-onset neurodegenerative diseases. It has been shown that its relocalization from the cytoplasm to the cell membrane is crucial to induce early endosomes maturation. In particular, evidences suggest that the N-terminal regulator of chromosome condensation 1 like domain (RLD) is necessary for membrane association thanks to its affinity to phosphoinositides, membrane lipids involved in the regulation of several signaling processes. Interestingly, this domain showed affinity towards phosphatidylinositol 3-phosphate [PI(3)P], which is highly expressed in endosomes membrane. However, Alsin structure has not been experimentally resolved yet and molecular mechanisms associated with its biological functions are mostly unknown. In this work, Alsin RLD has been investigated through computational molecular modeling techniques to analyze its conformational dynamics and obtain a representative 3D model of this domain. Moreover, a putative phosphoinositide binding site has been proposed and PI(3)P interaction mechanism studied. Results highlight the substantial conformational stability of Alsin RLD secondary structure and suggest the role of one highly flexible region in the phosphoinositides selectivity of this domain.

Extended diffusion theory: Recovering dynamics from biased/accelerated molecular simulations

Dynamical properties are of great importance in determining the behavior of synthetic and natural molecules, but capturing them by computational methods is a nontrivial task. Very often the time scales of the relevant phenomena are far beyond the typical time windows accessible by classical Molecular Dynamics (MD) simulations, currently limited to the order of microseconds on standard laboratory workstations. On the other hand, biased and accelerated simulations allow for fast and thorough exploration of the molecular conformational space, but they lose the dynamic information. The problem of recovering dynamics from biased/accelerated simulations is a very active field of research, but no totally robust/reliable solutions have been given yet. In this paper it is shown how the Smoluchowski equation, in the framework of Diffusion Theory (DT), can be used to bridge this gap, and dynamical properties, in the form of time correlation functions (TCFs), can be extracted also from such kind of simulations. DT is first extended (EDT) to express the mobility tensors entering the Smoluchowski operator in terms of a recently introduced unified and regularized Rotne-Prager-Yamakawa approximation, [P. J. Zuk, E. Wajnryb, K. A. Mizerski, P. Szymczak, J. Fluid. Mech. 2014, 741, R5, 1-13] also involving mixed rotation-translation contributions, and rotation-rotation terms beside the classical translation-translation ones, so far used in DT. Then, the method is applied to recover the dynamics of a nontrivial example of a peptide in explicit water from the first 200 ns of a Replica Exchange Molecular Dynamics simulation, which is a popular computational method that destroys the long time dynamics. EDT dynamics were found to favorably compare against those coming from a standard MD simulation of the same system, requiring a time window of 30 μs to converge. This result shows that EDT is a tool of practical value to recover the long time dynamics of systems in diffusive regimes from biased/accelerated simulations, to be exploited in those cases when direct evaluation by standard MD is unfeasible.

Molecular mechanisms of metal ions in regulating the catalytic efficiency of D-psicose 3-epimerase revealed by multiple short molecular dynamic simulations and free energy predictions

Communicated by Ramaswamy H. Sarma.

Aminoacid substitutions in the glycine zipper affect the conformational stability of amyloid beta fibrils

The aggregation of amyloid-beta peptides is associated with the pathogenesis of Alzheimer's disease. The hydrophobic core of the amyloid beta sequence contains a GxxxG repeated motif, called glycine zipper, which involves crucial residues for assuring stability and promoting the process of fibril formation. Mutations in this motif lead to a completely different oligomerization pathway and rate of fibril formation. In this work, we have tested G33L and G37L residue substitutions by molecular dynamics simulations. We found that both protein mutations may lead to remarkable changes in the fibril conformational stability. Results suggest the disruption of the glycine zipper as a possible strategy to reduce the aggregation propensity of amyloid beta peptides. On the basis of our data, further investigations may consider this key region as a binding site to design/discover novel effective inhibitors.Communicated by Ramaswamy H. Sarma.

The Extent of Human Apolipoprotein A-I Lipidation Strongly Affects the β-Amyloid Efflux Across the Blood-Brain Barrier in vitro

Much evidence suggests a protective role of high-density lipoprotein (HDL) and its major apolipoprotein apoA-I, in Alzheimer's disease (AD). The biogenesis of nascent HDL derived from a first lipidation of apoA-I, which is synthesized by the liver and intestine but not in the brain, in a process mediated by ABCA1. The maturation of nascent HDL in mature spherical HDL is due to a subsequent lipidation step, LCAT-mediated cholesterol esterification, and the change of apoA-I conformation. Therefore, different subclasses of apoA-I-HDL simultaneously exist in the blood circulation. Here, we investigated if and how the lipidation state affects the ability of apoA-I-HDL to target and modulate the cerebral β-amyloid (Aβ) content from the periphery, that is thus far unclear. In particular, different subclasses of HDL, each with different apoA-I lipidation state, were purified from human plasma and their ability to cross the blood-brain barrier (BBB), to interact with Aβ aggregates, and to affect Aβ efflux across the BBB was assessed in vitro using a transwell system. The results showed that discoidal HDL displayed a superior capability to promote Aβ efflux in vitro (9 × 10-5 cm/min), when compared to apoA-I in other lipidation states. In particular, no effect on Aβ efflux was detected when apoA-I was in mature spherical HDL, suggesting that apoA-I conformation, and lipidation could play a role in Aβ clearance from the brain. Finally, when apoA-I folded its structure in discoidal HDL, rather than in spherical ones, it was able to cross the BBB in vitro and strongly destabilize the conformation of Aβ fibrils by decreasing the order of the fibril structure (-24%) and the β-sheet content (-14%). These data suggest that the extent of apoA-I lipidation, and consequently its conformation, may represent crucial features that could exert their protective role in AD pathogenesis.

Destabilizing the AXH Tetramer by Mutations: Mechanisms and Potential Antiaggregation Strategies

The AXH domain of protein Ataxin 1 is thought to play a key role in the misfolding and aggregation pathway responsible for Spinocerebellar ataxia 1. For this reason, a molecular level understanding of AXH oligomerization pathway is crucial to elucidate the aggregation mechanism, which is thought to trigger the disease. This study employs classical and enhanced molecular dynamics to identify the structural and energetic basis of AXH tetramer stability. Results of this work elucidate molecular mechanisms behind the destabilizing effect of protein mutations, which consequently affect the AXH tetramer assembly. Moreover, results of the study draw attention for the first time, to our knowledge, to the R638 protein residue, which is shown to play a key role in AXH tetramer stability. Therefore, R638 might be also implicated in the AXH oligomerization pathway and stands out as a target for future experimental studies focused on self-association mechanisms and fibril formation of full-length ATX1.

Protein Environment: A Crucial Triggering Factor in Josephin Domain Aggregation: The Role of 2,2,2-Trifluoroethanol

The protein ataxin-3 contains a polyglutamine stretch that triggers amyloid aggregation when it is expanded beyond a critical threshold. This results in the onset of the spinocerebellar ataxia type 3. The protein consists of the globular N-terminal Josephin domain and a disordered C-terminal tail where the polyglutamine stretch is located. Expanded ataxin-3 aggregates via a two-stage mechanism: first, Josephin domain self-association, then polyQ fibrillation. This highlights the intrinsic amyloidogenic potential of Josephin domain. Therefore, much effort has been put into investigating its aggregation mechanism(s). A key issue regards the conformational requirements for triggering amyloid aggregation, as it is believed that, generally, misfolding should precede aggregation. Here, we have assayed the effect of 2,2,2-trifluoroethanol, a co-solvent capable of stabilizing secondary structures, especially α-helices. By combining biophysical methods and molecular dynamics, we demonstrated that both secondary and tertiary JD structures are virtually unchanged in the presence of up to 5% 2,2,2-trifluoroethanol. Despite the preservation of JD structure, 1% of 2,2,2-trifluoroethanol suffices to exacerbate the intrinsic aggregation propensity of this domain, by slightly decreasing its conformational stability. These results indicate that in the case of JD, conformational fluctuations might suffice to promote a transition towards an aggregated state without the need for extensive unfolding, and highlights the important role played by the environment on the aggregation of this globular domain.

Self-Assembled Ligands Targeting TLR7: A Molecular Level Investigation

Toll-like receptors (TLRs) are pattern recognition transmembrane proteins that play an important role in innate immunity. In particular, TLR7 plays a role in detecting nucleic acids derived from viruses and bacteria. The huge number of pathologies in which TLR7 is involved has led to an increasing interest in developing new compounds targeting this protein. Several conjugation strategies were proposed for TLR7 agonists to increase the potency while maintaining a low toxicity. In this work, we focus the attention on two promising classes of TLR7 compounds derived from the same pharmacophore conjugated with phospholipid and polyethylene glycol (PEG). A multidisciplinary investigation has been carried out by molecular dynamics (MD), dynamic light scattering (DLS), electron paramagnetic resonance (EPR), and cytotoxicity assessment. DLS and MD indicated how only the phospholipid conjugation provides the compound abilities to self-assemble in an orderly fashion with a maximal pharmacophore exposition to the solvent. Further EPR and cytotoxicity experiments highlighted that phospholipid compounds organize in stable aggregates and well interact with TLR7, whereas PEG conjugation was characterized by poorly stable aggregates at the cells surface. The methodological framework proposed in this study may be used to investigate, at a molecular level, the interactions generally occurring between aggregated ligands, to be used as drugs, and protein receptors.

Free energy landscape of siRNA-polycation complexation: Elucidating the effect of molecular geometry, polymer flexibility, and charge neutralization

The success of medical threatments with DNA and silencing interference RNA is strongly related to the design of efficient delivery technologies. Cationic polymers represent an attractive strategy to serve as nucleic-acid carriers with the envisioned advantages of efficient complexation, low cost, ease of production, well-defined size, and low polydispersity index. However, the balance between efficacy and toxicity (safety) of these polymers is a challenge and in need of improvement. With the aim of designing more effective polycationic-based gene carriers, many parameters such as carrier morphology, size, molecular weight, surface chemistry, and flexibility/rigidity ratio need to be taken into consideration. In the present work, the binding mechanism of three cationic polymers (polyarginine, polylysine and polyethyleneimine) to a model siRNA target is computationally investigated at the atomistic level. In order to better understand the polycationic carrier-siRNA interactions, replica exchange molecular dynamic simulations were carried out to provide an exhaustive exploration of all the possible binding sites, taking fully into account the siRNA flexibility together with the presence of explicit solvent and ions. Moreover, well-tempered metadynamics simulations were employed to elucidate how molecular geometry, polycation flexibility, and charge neutralization affect the siRNA-polycations free energy landscape in term of low-energy binding modes and unbinding free energy barriers. Significant differences among polymer binding modes have been detected, revealing the advantageous binding properties of polyarginine and polylysine compared to polyethyleneimine.

Insights into the Effect of the G245S Single Point Mutation on the Structure of p53 and the Binding of the Protein to DNA

The transcription factor p53 is a potent tumor suppressor dubbed as the "guardian of the genome" because of its ability to orchestrate protective biological outputs in response to a variety of oncogenic stresses. Mutation and thus inactivation of p53 can be found in 50% of human tumors. The majority are missense mutations located in the DNA binding region. Among them, G245S is known to be a structural hotspot mutation. To understand the behaviors and differences between the wild-type and mutant, both a dimer of the wild type p53 (wt-p53) and its G245S mutant (G245S-mp53), complexed with DNA, were simulated using molecular dynamics for more than 1 μs. wt-p53 and G245S-mp53 apo monomers were simulated for 1 μs as well. Conformational analyses and binding energy evaluations performed underline important differences and therefore provide insights to understand the G245S-mp53 loss of function. Our results indicate that the G245S mutation destabilizes several structural regions in the protein that are crucial for DNA binding when found in its apo form and highlight differences in the mutant-DNA complex structure compared to the wt protein. These findings not only provide means that can be applied to other p53 mutants but also serve as structural basis for further studies aimed at the development of cancer therapies based on restoring the function of p53.

Chaperones in Polyglutamine Aggregation: Beyond the Q-Stretch

Expanded polyglutamine (polyQ) stretches in at least nine unrelated proteins lead to inherited neuronal dysfunction and degeneration. The expansion size in all diseases correlates with age at onset (AO) of disease and with polyQ protein aggregation, indicating that the expanded polyQ stretch is the main driving force for the disease onset. Interestingly, there is marked interpatient variability in expansion thresholds for a given disease. Between different polyQ diseases the repeat length vs. AO also indicates the existence of modulatory effects on aggregation of the upstream and downstream amino acid sequences flanking the Q expansion. This can be either due to intrinsic modulation of aggregation by the flanking regions, or due to differential interaction with other proteins, such as the components of the cellular protein quality control network. Indeed, several lines of evidence suggest that molecular chaperones have impact on the handling of different polyQ proteins. Here, we review factors differentially influencing polyQ aggregation: the Q-stretch itself, modulatory flanking sequences, interaction partners, cleavage of polyQ-containing proteins, and post-translational modifications, with a special focus on the role of molecular chaperones. By discussing typical examples of how these factors influence aggregation, we provide more insight on the variability of AO between different diseases as well as within the same polyQ disorder, on the molecular level.

Thermodynamic and kinetic stability of the Josephin Domain closed arrangement: evidences from replica exchange molecular dynamics

BACKGROUND

Molecular phenomena driving pathological aggregation in neurodegenerative diseases are not completely understood yet. Peculiar is the case of Spinocerebellar Ataxia 3 (SCA3) where the conformational properties of the AT-3 N-terminal region, also called Josephin Domain (JD), play a key role in the first step of aggregation, having the JD an amyloidogenic propensity itself. For this reason, unraveling the intimate relationship between JD structural features and aggregation tendency may lead to a step forward in understanding the pathology and rationally design a cure. In this connection, computational modeling has demonstrated to be helpful in exploring the protein molecular dynamics and mechanism of action.

RESULTS

Conformational dynamics of the JD is here finely investigated by replica exchange molecular dynamics simulations able to sample the microsecond time scale and to provide both a thermodynamic and kinetic description of the protein conformational changes. Accessible structural conformations of the JD have been identified in: open, intermediate and closed like arrangement. Data indicated the closed JD arrangement as the most likely protein arrangement. The protein transition from closed toward intermediate/open states was characterized by a rate constant higher than 700 ns. This result also explains the inability of classical molecular dynamics to explore transitions from closed to open JD configuration on a time scale of hundreds of nanoseconds.

CONCLUSION

This work provides the first kinetic estimation of the JD transition pathway from open-like to closed-like arrangement and vice-versa, indicating the closed-like arrangement as the most likely configuration for a JD in water environment. More widely, the importance of our results is also underscored considering that the ability to provide a kinetic description of the protein conformational changes is a scientific challenge for both experimental and theoretical approaches to date.

REVIEWERS

This article was reviewed by Oliviero Carugo, Bojan Zagrovic.