β-sheet-like formation during the mechanical unfolding of prion protein

Interfacial interactions between spider silk protein and cellulose studied by molecular dynamics simulation

CONTEXT

Due to their excellent biocompatibility and degradability, cellulose/spider silk protein composites hold a significant value in biomedical applications such as tissue engineering, drug delivery, and medical dressings. The interfacial interactions between cellulose and spider silk protein affect the properties of the composite. Therefore, it is important to understand the interfacial interactions between spider silk protein and cellulose to guide the design and optimization of composites. The study of the adsorption of protein on specific surfaces of cellulose crystal can be very complex using experimental methods. Molecular dynamics simulations allow the exploration of various physical and chemical changes at the atomic level of the material and enable an atomic description of the interactions between cellulose crystal planes and spider silk protein. In this study, molecular dynamics simulations were employed to investigate the interfacial interactions between spider silk protein (NTD) and cellulose surfaces. Findings of RMSD, RMSF, and secondary structure showed that the structure of NTD proteins remained unchanged during the adsorption process. Cellulose contact numbers and hydrogen bonding trends on different crystalline surfaces suggest that van der Waals forces and hydrogen bonding interactions drive the binding of proteins to cellulose. These findings reveal the interaction between cellulose and protein at the molecular level and provide theoretical guidance for the design and synthesis of cellulose/spider silk protein composites.

METHODS

MD simulations were all performed using the GROMACS-5.1 software package and run with CHARMM36 carbohydrate force field. Molecular dynamics simulations were performed for 500 ns for the simulated system.

Unraveling the Histidine Tautomerism Effect on the Initial Stages of Prion Misfolding: New Insights from a Computational Perspective

The aggregation and structural conversion of normal prion peptide (PrP^C) into the pathogenic scrapie form (PrP^Sc), which can act as a seed to enhance prion amyloid fiber formation, is believed to be a crucial event in prionopathies. Previous research suggests that the prion monomer may play an important role in oligomer generation during disease pathogenesis. In the present study, extensive replica-exchange molecular dynamics (REMD) simulations were conducted to explore the conformational characteristics of the huPrP (125-160) monomer under the histidine tautomerism effect. Investigating the structural characteristics and fibrilization process is challenging because two histidine tautomers [N_ε2-H (ε) and N_δ1-H (δ)] can occur in the open neutral state. Molecular dynamics (MD) simulation outcomes have shown that the toxic εδ and δδ isomer (containing several and broader local minima) had the highest α-helix structures, with contents of 21.11% and 21.01%, respectively, and may have a strong influence on the organizational behavior of a monomeric prion. The amino acids aspartate 20 (D20)-asparagine 29 (N29) and isoleucine 15 (I15)-histidine 16 (H16), D20-arginine 27 (R27) as well as N29 formed α-helix with the highest probabilities in the δδ and εδ isomer, accordingly. On the basis of our findings, we propose the histidine tautomerization hypothesis as a new prion accumulation mechanism, which may exist to induce the formation of prion accumulates. Overall, our tautomerism hypothesis constitutes a promising perspective for enhancing understanding of prion disease pathobiology and may help in the design of a good inhibitor.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Atomistic Simulation of the Rate-Dependent Ductile-to-Brittle Failure Transition in Bicrystalline Metal Nanowires

The mechanical properties and plastic deformation mechanisms of metal nanowires have been studied intensely for many years. One of the important yet unresolved challenges in this field is to bridge the gap in properties and deformation mechanisms reported for slow strain rate experiments (∼10^-2 s^-1), and high strain rate molecular dynamics (MD) simulations (∼10⁸ s^-1) such that a complete understanding of strain rate effects on mechanical deformation and plasticity can be obtained. In this work, we use long time scale atomistic modeling based on potential energy surface exploration to elucidate the atomistic mechanisms governing a strain-rate-dependent incipient plasticity and yielding transition for face centered cubic (FCC) copper and silver nanowires. The transition occurs for both metals with both pristine and rough surfaces for all computationally accessible diameters (<10 nm). We find that the yield transition is induced by a transition in the incipient plastic event from Shockley partials nucleated on primary slip systems at MD strain rates to the nucleation of planar defects on non-Schmid slip planes at experimental strain rates, where multiple twin boundaries and planar stacking faults appear in copper and silver, respectively. Finally, we demonstrate that, at experimental strain rates, a ductile-to-brittle transition in failure mode similar to previous experimental studies on bicrystalline silver nanowires is observed, which is driven by differences in dislocation activity and grain boundary mobility as compared to the high strain rate case.

Mechanical Deformation Mechanisms and Properties of Prion Fibrils Probed by Atomistic Simulations

Prion fibrils, which are a hallmark for neurodegenerative diseases, have recently been found to exhibit the structural diversity that governs disease pathology. Despite our recent finding concerning the role of the disease-specific structure of prion fibrils in determining their elastic properties, the mechanical deformation mechanisms and fracture properties of prion fibrils depending on their structures have not been fully characterized. In this work, we have studied the tensile deformation mechanisms of prion and non-prion amyloid fibrils by using steered molecular dynamics simulations. Our simulation results show that the elastic modulus of prion fibril, which is formed based on left-handed β-helical structure, is larger than that of non-prion fibril constructed based on right-handed β-helix. However, the mechanical toughness of prion fibril is found to be less than that of non-prion fibril, which indicates that infectious prion fibril is more fragile than non-infectious (non-prion) fibril. Our study sheds light on the role of the helical structure of amyloid fibrils, which is related to prion infectivity, in determining their mechanical deformation mechanisms and properties.

N-Terminal Acetylation Preserves α-Synuclein from Oligomerization by Blocking Intermolecular Hydrogen Bonds

The abnormal aggregation of α-synuclein (α-Syn) is closely associated with Parkinson's disease. Different post-translational modifications of α-Syn have been identified and contribute distinctly in α-Syn aggregation and cytotoxicity. Recently, α-Syn was reported to be N-terminally acetylated in cells, yet the functional implication of this modification, especially in α-Syn oligomerization, remains unclear. By using a solid-state nanopore system, we found that N-terminal acetylation can significantly decrease α-Syn oligomerization. Replica-exchange molecular dynamics simulations further revealed that addition of an acetyl group at the N-terminus disrupts intermolecular hydrogen bonds, which slows down the initial α-Syn oligomerization. Our finding highlights the essential role of N-terminal acetylation of α-Syn in preserving its native conformation against pathological aggregation.