Gianvito Grasso

Josephin Domain Structural Conformations Explored by Metadynamics in Essential Coordinates

The Josephin Domain (JD), i.e. the N-terminal domain of Ataxin 3 (At3) protein, is an interesting example of competition between physiological function and aggregation risk. In fact, the fibrillogenesis of Ataxin 3, responsible for the spinocerebbellar ataxia 3, is strictly related to the JD thermodynamic stability. Whereas recent NMR studies have demonstrated that different JD conformations exist, the likelihood of JD achievable conformational states in solution is still an open issue. Marked differences in the available NMR models are located in the hairpin region, supporting the idea that JD has a flexible hairpin in dynamic equilibrium between open and closed states. In this work we have carried out an investigation on the JD conformational arrangement by means of both classical molecular dynamics (MD) and Metadynamics employing essential coordinates as collective variables. We provide a representation of the free energy landscape characterizing the transition pathway from a JD open-like structure to a closed-like conformation. Findings of our in silico study strongly point to the closed-like conformation as the most likely for a Josephin Domain in water.

Investigation of the Josephin Domain Protein-Protein Interaction by Molecular Dynamics

Spinocerebellar ataxia (SCA) 3, the most common form of SCA, is a neurodegenerative rare disease characterized by polyglutamine tract expansion and self-assembly of Ataxin3 (At3) misfolded proteins into highly organized fibrillar aggregates. The At3 N-terminal Josephin Domain (JD) has been suggested as being responsible for mediating the initial phase of the At3 double-step fibrillogenesis. Several issues concerning the residues involved in the JD’s aggregation and, more generally, the JD clumping mechanism have not been clarified yet. In this paper we present an investigation focusing on the JD protein-protein interaction by means of molecular modeling. Our results suggest possible aminoacids involved in JD contact together with local and non-local effects following JD dimerization. Surprisingly, JD conformational changes following the binding may involve ubiquitin binding sites and hairpin region even though they do not pertain to the JD interaction surfaces. Moreover, the JD binding event has been found to alter the hairpin open-like conformation toward a closed-like arrangement over the simulated timescale. Finally, our results suggest that the JD aggregation might be a multi-step process, with an initial fast JD-JD binding mainly driven by Arg101, followed by slower structural global rearrangements involving the exposure to the solvent of Leu84-Trp87, which might play a role in a second step of JD aggregation.

Cell Penetrating Peptide Adsorption on Magnetite and Silica Surfaces: A Computational Investigation.

Magnetic nanoparticles (MNPs) represent one of the most promising materials as they can act as a versatile platform in the field of bionanotechnology for enhanced imaging, diagnosis, and treatment of various diseases. Silica is the most common compound for preparing coated iron oxide NPs since it improves colloidal stability and the binding affinity for various organic molecules. Biomolecules such as cell penetrating peptides (CPPs) might be employed to decorate MNPs, combining their promising physicochemical properties with a cell penetrating ability. In this work, a computational investigation on adsorption of Antennapedia homeodomain-derived penetrating peptide (pAntp) on silica and magnetite (MAG) surfaces is presented. By employing umbrella sampling molecular dynamics, we provided a quantitative estimation of the pAntp-surface adsorption free energy to highlight the influence of surface hydroxylation state on the adsorption mechanism. The interaction between peptide and surface has shown to be mainly driven by electrostatics. In case of MAG surface, also an important contribution of van der Waals (VdW) attraction was observed. Our data suggest that a competitive mechanism between MNPs and cell membrane might partially inhibit the CPP to carry out its membrane penetrating function.

Conformational fluctuations of the AXH monomer of Ataxin-1

In this paper, we report the results of molecular dynamics simulations of AXH monomer of Ataxin‐1. The AXH domain plays a crucial role in Ataxin‐1 aggregation, which accompanies the initiation and progression of Spinocerebellar ataxia type 1. Our simulations involving both classical and replica exchange molecular dynamics, followed by principal component analysis of the trajectories obtained, reveal substantial conformational fluctuations of the protein structure, especially in the N‐terminal region. We show that these fluctuations can be generated by thermal noise since the free energy barriers between conformations are small enough for thermally stimulated transitions. In agreement with the previous experimental findings, our results can be considered as a basis for a future design of ataxin aggregation inhibitors that will require several key conformations identified in the present study as molecular targets for ligand binding. Proteins 2016; 84:52–59.

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

Ataxin‐1 is the protein responsible for the Spinocerebellar ataxia type 1, an incurable neurodegenerative disease caused by polyglutamine expansion. The AXH domain plays a pivotal role in physiological functions of Ataxin‐1. In Spinocerebellar ataxia 1, the AXH domain is involved in the misfolding and aggregation pathway. Here molecular modeling is applied to investigate the protein‐protein interactions contributing to the AXH dimer stability. Particular attention is focused on: (i) the characterization of AXH monomer‐monomer interface; (ii) the molecular description of the AXH monomer‐monomer interaction dynamics. Technically, an approach based on functional mode analysis, here applied to replica exchange molecular dynamics trajectories, was employed. The findings of this study are consistent with previous experimental results and elucidate the pivotal role of the I580 residue in mediating the AXH monomer‐monomer interaction dynamics. Proteins 2016; 84:666–673.

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.

Multivalent interacting glycodendrimer to prevent amyloid-peptide fibril formation induced by Cu(II): A multidisciplinary approach

Amyloid peptide fibrillogenesis induced by Cu(II) ions is a key event in the pathogenesis of Alzheimer’s disease. Dendrimers have been found to be active in preventing fibril formation. Therefore, they hold promise for the treatment of Alzheimer’s disease. In this study, the fibrillation mechanism of amyloid peptide Aβ 1-40 was studied by adding Cu(II) in the absence and presence of 4th generation poly(propyleneimine) glycodendrimer functionalized with sulfate groups, using dynamic light scattering (DLS), circular dichroism (CD), fluorescence, electron paramagnetic resonance (EPR) and molecular modeling (MD). The glycodendrimer was non-toxic to mHippoE-18 embryonic mouse hippocampal cells, selected as a nerve cell model, and decreased the toxicity of peptide aggregates formed after the addition of Cu(II). The binary systems of Cu(II)–glycodendrimer, Cu(II)–peptide, and glycodendrimer–peptide were first characterized. At the lowest Cu(II)/glycodendrimer molar ratios, Cu(II) was complexed by the internal-dendrimer nitrogen sites. After saturation of these sites, Cu(II) binding with sulfate groups occurred. Stable Cu(II)–peptide complexes formed within 5 min and were responsible for a transition from an α helix to a β-sheet conformation of Aβ 1-40. Glycodendrimer–peptide interactions provoked the stabilization of the α-helix, as demonstrated in the absence of Cu(II) by the Thioflavin T assay, and in the presence of Cu(II) by CD, EPR, and MD. Formation of fibrils is differentially modulated by glycodendrimer and Cu(II) concentrations for a fixed amount of Aβ 1-40. Therefore, this multidisciplinary study facilitated the recognition of optimal experimental conditions that allow the glycodendrimer to avoid the fibril formation induced by Cu(II).

Conformational Dynamics and Stability of U-Shaped and S-Shaped Amyloid β Assemblies

Alzheimer’s disease is the most fatal neurodegenerative disorder characterized by the aggregation and deposition of Amyloid β (Aβ) oligomers in the brain of patients. Two principal variants of Aβ exist in humans: Aβ1–40 and Aβ1–42. The former is the most abundant in the plaques, while the latter is the most toxic species and forms fibrils more rapidly. Interestingly, fibrils of Aβ1–40 peptides can only assume U-shaped conformations while Aβ1–42 can also arrange as S-shaped three-stranded chains, as recently discovered. As alterations in protein conformational arrangement correlate with cell toxicity and speed of disease progression, it is important to characterize, at molecular level, the conformational dynamics of amyloid fibrils. In this work, Replica Exchange Molecular Dynamics simulations were carried out to compare the conformational dynamics of U-shaped and S-shaped Aβ17–42 small fibrils. Our computational results provide support for the stability of the recently proposed S-shaped model due to the maximized interactions involving the C-terminal residues. On the other hand, the U-shaped motif is characterized by significant distortions resulting in a more disordered assembly. Outcomes of our work suggest that the molecular architecture of the protein aggregates might play a pivotal role in formation and conformational stability of the resulting fibrils.

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.

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.