Investigation of the Josephin Domain protein-protein interaction by molecular dynamics

The Josephin domain (JD) containing proteins are predicted to bind to the same interactors: Implications for spinocerebellar ataxia type 3 (SCA3) studies using Drosophila melanogaster mutants

Spinocerebellar ataxia type 3, also known as Machado-Joseph disease (SCA3/ MJD), is the most frequent polyglutamine (polyQ) neurodegenerative disorder. It is caused by a pathogenic expansion of the polyQ tract, located at the C-terminal region of the protein encoded by the ATXN3 gene. This gene codes for a deubiquitinating enzyme (DUB) that belongs to a gene family, that in humans is composed by three more genes (ATXN3L, JOSD1, and JOSD2), that define two gene lineages (the ATXN3 and the Josephins). These proteins have in common the N-terminal catalytic domain (Josephin domain, JD), that in Josephins is the only domain present. In ATXN3 knock-out mouse and nematode models, the SCA3 neurodegeneration phenotype is not, however, reproduced, suggesting that in the genome of these species there are other genes that are able to compensate for the lack of ATXN3. Moreover, in mutant Drosophila melanogaster, where the only JD protein is coded by a Josephin-like gene, expression of the expanded human ATXN3 gene reproduces multiple aspects of the SCA3 phenotype, in contrast with the results of the expression of the wild type human form. In order to explain these findings, phylogenetic, as well as, protein–protein docking inferences are here performed. Here we show multiple losses of JD containing genes across the animal kingdom, suggesting partial functional redundancy of these genes. Accordingly, we predict that the JD is essential for binding with ataxin-3 and proteins of the Josephin lineages, and that D. melanogaster mutants are a good model of SCA3 despite the absence of a gene from the ATXN3 lineage. The molecular recognition regions of the ataxin-3 binding and those predicted for the Josephins are, however, different. We also report different binding regions between the two ataxin-3 forms (wild-type (wt) and expanded (exp)). The interactors that show an increase in the interaction strength with exp ataxin-3, are enriched in extrinsic components of mitochondrial outer membrane and endoplasmatic reticulum membrane. On the other hand, the group of interactors that show a decrease in the interaction strength with exp ataxin-3 is significantly enriched in extrinsic component of cytoplasm.

Prediction of Protein–Protein Interactions Between Alsin DH/PH and Rac1 and Resulting Protein Dynamics

Alsin is a protein of 1,657 amino acids known for its crucial role in vesicular trafficking in neurons thanks to its ability to interact with two guanosine triphosphatases, Rac1 and Rab5. Evidence suggests that Rac1 can bind Alsin central region, composed by a Dbl Homology (DH) domain followed by a Pleckstrin Homology (PH) domain, leading to Alsin relocalization. However, Alsin three-dimensional structure and its relationship with known biological functions of this protein are still unknown. In this work, a homology model of the Alsin DH/PH domain was developed and studied through molecular dynamics both in the presence and in the absence of its binding partner, Rac1. Due to different conformations of DH domain, the presence of Rac1 seems to stabilize an open state of the protein, while the absence of its binding partner results in closed conformations. Furthermore, Rac1 interaction was able to reduce the fluctuations in the second conserved region of DH motif, which may be involved in the formation of a homodimer. Moreover, the dynamics of DH/PH was described through a Markov State Model to study the pathways linking the open and closed states. In conclusion, this work provided an all-atom model for the DH/PH domain of Alsin protein; moreover, molecular dynamics investigations suggested underlying molecular mechanisms in the signal transduction between Rac1 and Alsin, providing the basis for a deeper understanding of the whole structure–function relationship for Alsin protein.

Structural, energetic and lipophilic analysis of SARS-CoV-2 non-structural protein 9 (NSP9)

In SARS-CoV-2 replication complex, the Non-structural protein 9 (Nsp9) is an important RNA binding subunit in the RNA-synthesizing machinery. The dimeric forms of coronavirus Nsp9 increase their nucleic acid binding affinity and the N-finger motif appears to play a critical role in dimerization. Here, we present a structural, lipophilic and energetic study about the Nsp9 dimer of SARS-CoV-2 through computational methods that complement hydrophobicity scales of amino acids with molecular dynamics simulations. Additionally, we presented a virtual N-finger mutation to investigate whether this motif contributes to dimer stability. The results reveal for the native dimer that the N-finger contributes favorably through hydrogen bond interactions and two amino acids bellowing to the hydrophobic region, Leu45 and Leu106, are crucial in the formation of the cavity for potential drug binding. On the other hand, Gly100 and Gly104, are responsible for stabilizing the α-helices and making the dimer interface remain stable in both, native and mutant (without N-finger motif) systems. Besides, clustering results for the native dimer showed accessible cavities to drugs. In addition, the energetic and lipophilic analysis reveal that the higher binding energy in the native dimer can be deduced since it is more lipophilic than the mutant one, increasing non-polar interactions, which is in line with the result of MM-GBSA and SIE approaches where the van der Waals energy term has the greatest weight in the stability of the native dimer. Overall, we provide a detailed study on the Nsp9 dimer of SARS-CoV-2 that may aid in the development of new strategies for the treatment and prevention of COVID-19.

Mutant Ataxin-3-Containing Aggregates (MATAGGs) in Spinocerebellar Ataxia Type 3: Dynamics of the Disorder

Spinocerebellar ataxia type 3 (SCA3) is the most common type of SCA worldwide caused by abnormal polyglutamine expansion in the coding region of the ataxin-3 gene. Ataxin-3 is a multi-faceted protein involved in various cellular processes such as deubiquitination, cytoskeletal organisation, and transcriptional regulation. The presence of an expanded poly(Q) stretch leads to altered processing and misfolding of the protein culminating in the production of insoluble protein aggregates in the cell. Various post-translational modifications affect ataxin-3 fibrillation and aggregation. This review provides an exhaustive assessment of the various pathogenic mechanisms undertaken by the mutant ataxin-3-containing aggregates (MATAGGs) for disease induction and neurodegeneration. This includes in-depth discussion on MATAGG dynamics including their formation, role in neuronal pathogenesis, and the debate over the toxic v/s protective nature of the MATAGGs in disease progression. Additionally, the currently available therapeutic strategies against SCA3 have been reviewed. The shift in the focus of such strategies, from targeting the steps that lead to or reduce aggregate formation to targeting the expression of mutant ataxin-3 itself via RNA-based therapeutics, has also been presented. We also discuss the intriguing promise that various growth and neurotrophic factors, especially the insulin pathway, hold in the modulation of SCA3 progression. These emerging areas show the newer directions through which SCA3 can be targeted including various preclinical and clinical trials. All these advances made in the last three decades since the discovery of the ataxin-3 gene have been critically reviewed here.

Insights into the interaction dynamics between volatile anesthetics and tubulin through computational molecular modelling

General anesthetics, able to reversibly suppress all conscious brain activity, have baffled medical science for decades, and little is known about their exact molecular mechanism of action. Given the recent scientific interest in the exploration of microtubules as putative functional targets of anesthetics, and the involvement thereof in neurodegenerative disorders, the present work focuses on the investigation of the interaction between human tubulin and four volatile anesthetics: ethylene, desflurane, halothane and methoxyflurane. Interaction sites on different tubulin isotypes are predicted through docking, along with an estimate of the binding affinity ranking. The analysis is expanded by Molecular Dynamics simulations, where the dimers are allowed to freely interact with anesthetics in the surrounding medium. This allowed for the determination of interaction hotspots on tubulin dimers, which could be linked to different functional consequences on the microtubule architecture, and confirmed the weak, Van der Waals-type interaction, occurring within hydrophobic pockets on the dimer. Both docking and MD simulations highlighted significantly weaker interactions of ethylene, consistent with its far lower potency as a general anesthetic. Overall, simulations suggest a transient interaction between anesthetics and microtubules in general anesthesia, and contact probability analysis shows interaction strengths consistent with the potencies of the four compounds.Communicated by Ramaswamy H. Sarma.

A new method for protein characterization and classification using geometrical features for 3D face analysis: An example of tubulin structures

This article reports on the results of research aimed to translate biometric 3D face recognition concepts and algorithms into the field of protein biophysics in order to precisely and rapidly classify morphological features of protein surfaces. Both human faces and protein surfaces are free-forms and some descriptors used in differential geometry can be used to describe them applying the principles of feature extraction developed for computer vision and pattern recognition. The first part of this study focused on building the protein dataset using a simulation tool and performing feature extraction using novel geometrical descriptors. The second part tested the method on two examples, first involved a classification of tubulin isotypes and the second compared tubulin with the FtsZ protein, which is its bacterial analog. An additional test involved several unrelated proteins. Different classification methodologies have been used: a classic approach with a support vector machine (SVM) classifier and an unsupervised learning with a k-means approach. The best result was obtained with SVM and the radial basis function kernel. The results are significant and competitive with the state-of-the-art protein classification methods. This leads to a new methodological direction in protein structure analysis.

Highlighting the effect of amyloid beta assemblies on the mechanical properties and conformational stability of cell membrane

Alzheimer disease (AD) is the most common cause of dementia, characterized by a progressive decline in cognitive function due to the abnormal aggregation and deposition of Amyloid beta (Aβ) fibrils in the brain of patients. In this context, the molecular mechanisms of protein misfolding and aggregation that are known to induce significant biophysical alterations in cells, including destabilization of plasma membranes, remain partially unclear. Physical interaction between the Aβ assemblies and the membrane leads to the disruption of the cell membrane in multiple ways including, surface carpeting, generation of transmembrane channels and detergent-like membrane dissolution. Understanding the impact of amyloidogenic protein in different stages of aggregation with the plasma membrane, plays a crucial role to fully elucidate the pathological mechanisms of AD. Within this framework, computer simulations represent a powerful tool able to shed lights on the interactions governing the structural influence of Aβ proteins on biological membrane. In this study, molecular dynamics (MD) simulations have been performed in order to characterize how POPC bilayer conformational and mechanical properties are affected by the interaction with Aβ_11-42 peptide, oligomer and fibril.

Molecular dynamics simulations to the bidirectional adhesion signaling pathway of integrin αV β3

The bidirectional force transmission process of integrin through the cell membrane is still not well understood. Several possible mechanisms have been discussed in literature on the basis of experimental data, and in this study, we investigate these mechanisms by free and steered molecular dynamics simulations. For the first time, constant velocity pulling on the complete integrin molecule inside a dipalmitoyl-phosphatidylcholine membrane is conducted. From the results, the most likely mechanism for inside-out and outside-in signaling is the switchblade model with further separation of the transmembrane helices.

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.

Elucidating the role of surface chemistry on cationic phosphorus dendrimer-siRNA complexation

In the field of dendrimers targeting small interfering RNA (siRNA) delivery, dendrimer structural properties, such as the flexibility/rigidity ratio, play a crucial role in the efficiency of complexation. However, advances in organic chemistry have enabled the development of dendrimers that differ only by a single atom on their surface terminals. This is the case for cationic phosphorus dendrimers functionalized with either pyrrolidinium (DP) or morpholinium (DM) terminal groups. This small change was shown to strongly affect the dendrimer-siRNA complexation, leading to more efficient anti-inflammatory effects in the case of DP. Reasons for this different behavior can hardly be inferred only by biological in vitro and in vivo experiments due to the high number of variables and complexity of the investigated biological system. However, an understanding of how small chemical surface changes may completely modify the overall dendrimer-siRNA complexation is a significant breakthrough towards the design of efficient dendrimers for nucleic acid delivery. Herein, we present experimental and computational approaches based on isothermal titration calorimetry and molecular dynamics simulations to elucidate the molecular reasons behind different efficiencies and activities of DP and DM. Results of the present research highlight how chemical surface modifications may drive the overall dendrimer-siRNA affinity by influencing enthalpic and entropic contributions of binding free energy. Moreover, this study elucidates molecular reasons related to complexation stoichiometry that may be crucial in determining the dendrimer complexation efficiency.

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.

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.

Amyloidogenicity at a Distance: How Distal Protein Regions Modulate Aggregation in Disease

The misfolding of proteins to form amyloid is a key pathological feature of several progressive, and currently incurable, diseases. A mechanistic understanding of the pathway from soluble, native protein to insoluble amyloid is crucial for therapeutic design, and recent efforts have helped to elucidate the key molecular events that trigger protein misfolding. Generally, either global or local structural perturbations occur early in amyloidogenesis to expose aggregation-prone regions of the protein that can then self-associate to form toxic oligomers. Surprisingly, these initiating structural changes are often caused or influenced by protein regions distal to the classically amyloidogenic sequences. Understanding the importance of these distal regions in the pathogenic process has highlighted many remaining knowledge gaps regarding the precise molecular events that occur in classic aggregation pathways. In this review, we discuss how these distal regions can influence aggregation in disease and the recent technical and conceptual advances that have allowed this insight.

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.

Underlying mechanisms of cyclic peptide inhibitors interrupting the interaction of CK2α/CK2β: comparative molecular dynamics simulation studies

Protein-protein interactions (PPIs) are fundamental to all biological processes. Recently, the CK2β-derived cyclic peptide Pc has been demonstrated to efficiently antagonize the CK2α/CK2β interaction and strongly affect the phosphorylation of CK2β-dependent CK2 substrate specificity. The binding affinity of Pc to CK2α is destroyed to different extents by two single-point mutations of Tyr188 to Ala (Y188A) and Phe190 to Ala (F190A), which exert negative effects on the inhibitory activity (IC50) of Pc against the CK2α/CK2β interaction from 3.0 μM to 54.0 μM and ≫100 μM, respectively. However, the structural influences of Y188A and F190A mutations on the CK2α-Pc complex remain unclear. In this study, comparative molecular dynamics (MD) simulations, principal component analysis (PCA), domain cross-correlation map (DCCM) analysis and energy calculations were performed on wild type (WT), Y188A mutant, and F190A mutant systems. The results revealed that ordered communications between hydrophobic and polar interactions were essential for CK2α-Pc binding in the WT system. In addition to the loss of the hydrogen bond between Gln36 of CK2α and Gly189 of Pc in the two mutants, the improper recognition mechanisms occurred through different pathways. These pathways included the weakened hydrophobic interactions in the Y188A mutant as well as decreased polar and hydrophobic interactions in the F190A mutant. The energy analysis results qualitatively elucidated the instability of the two mutants and energetic contributions of the key residues. This study not only revealed the structural mechanisms for the decreased binding affinity of Y188A and F190A mutant CK2α-Pc complexes, but also provided valuable clues for the rational design of CK2α/CK2β subunit interaction inhibitors with high affinity and specificity.

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.

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.

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.