Sandipan Mohanty

Structural Reorganisation and Potential Toxicity of Oligomeric Species Formed during the Assembly of Amyloid Fibrils

Increasing evidence indicates that oligomeric protein assemblies may represent the molecular species responsible for cytotoxicity in a range of neurological disorders including Alzheimer and Parkinson diseases. We use all-atom computer simulations to reveal that the process of oligomerization can be divided into two steps. The first is characterised by a hydrophobic coalescence resulting in the formation of molten oligomers in which hydrophobic residues are sequestered away from the solvent. In the second step, the oligomers undergo a process of reorganisation driven by interchain hydrogen bonding interactions that induce the formation of β sheet rich assemblies in which hydrophobic groups can become exposed. Our results show that the process of aggregation into either ordered or amorphous species is largely determined by a competition between the hydrophobicity of the amino acid sequence and the tendency of polypeptide chains to form arrays of hydrogen bonds. We discuss how the increase in solvent-exposed hydrophobic surface resulting from such a competition offers an explanation for recent observations concerning the cytotoxicity of oligomeric species formed prior to mature amyloid fibrils.

Oligomerization of Amyloid Aβ16–22 Peptides Using Hydrogen Bonds and Hydrophobicity Forces

The 16-22 amino-acid fragment of the beta-amyloid peptide associated with the Alzheimers disease, Abeta, is capable of forming amyloid fibrils. Here we study the aggregation mechanism of Abeta16-22 peptides by unbiased thermodynamic simulations at the atomic level for systems of one, three, and six Abeta16-22 peptides. We find that the isolated Abeta16-22 peptide is mainly a random coil in the sense that both the alpha-helix and beta-strand contents are low, whereas the three- and six-chain systems form aggregated structures with a high beta-sheet content. Furthermore, in agreement with experiments on Abeta16-22 fibrils, we find that large parallel beta-sheets are unlikely to form. For the six-chain system, the aggregated structures can have many different shapes, but certain particularly stable shapes can be identified.

PROFASI: A Monte Carlo simulation package for protein folding and aggregation†

We present a flexible and efficient program package written in C++, PROFASI, for simulating protein folding and aggregation. The systems are modeled using an all‐atom description of the protein chains with only torsional degrees of freedom, and implicit water. The program package has a modular structure that makes the interaction potential easy to modify. The currently implemented potential is able to fold several peptides with about 20 residues, and has also been used to study aggregation and force‐induced unfolding. The simulation methods implemented in PROFASI are Monte Carlo‐based and include a semilocal move and simulated tempering. Adding new updates is easy. The code runs fast in both single‐ and multi‐chain applications, as is illustrated by several examples.

Formation and Growth of Oligomers: A Monte Carlo Study of an Amyloid Tau Fragment

Small oligomers formed early in the process of amyloid fibril formation may be the major toxic species in Alzheimers disease. We investigate the early stages of amyloid aggregation for the tau fragment AcPHF6 (Ac-VQIVYK-NH2) using an implicit solvent all-atom model and extensive Monte Carlo simulations of 12, 24, and 36 chains. A variety of small metastable aggregates form and dissolve until an aggregate of a critical size and conformation arises. However, the stable oligomers, which are β-sheet-rich and feature many hydrophobic contacts, are not always growth-ready. The simulations indicate instead that these supercritical oligomers spend a lengthy period in equilibrium in which considerable reorganization takes place accompanied by exchange of chains with the solution. Growth competence of the stable oligomers correlates with the alignment of the strands in the β-sheets. The larger aggregates seen in our simulations are all composed of two twisted β-sheets, packed against each other with hydrophobic side chains at the sheet–sheet interface. These β-sandwiches show similarities with the proposed steric zipper structure for PHF6 fibrils but have a mixed parallel/antiparallel β-strand organization as opposed to the parallel organization found in experiments on fibrils. Interestingly, we find that the fraction of parallel β-sheet structure increases with aggregate size. We speculate that the reorganization of the β-sheets into parallel ones is an important rate-limiting step in the formation of PHF6 fibrils.

Folding Thermodynamics of Peptides

A simplified interaction potential for protein folding studies at the atomic level is discussed and tested on a set of peptides with approximately 20 residues each. The test set contains both alpha-helical (Trp cage, F(s)) and beta-sheet (GB1p, GB1m2, GB1m3, Betanova, LLM) peptides. The model, which is entirely sequence-based, is able to fold these different peptides for one and the same choice of model parameters. Furthermore, the melting behavior of the peptides is in good quantitative agreement with experimental data. Apparent folded populations obtained using different observables are compared, and are found to be very different for some of the peptides (e.g., Betanova). In other cases (in particular, GB1m2 and GB1m3), the different estimates agree reasonably well, indicating a more two-state-like melting behavior.

An effective all-atom potential for proteins

We describe and test an implicit solvent all-atom potential for simulations of protein folding and aggregation. The potential is developed through studies of structural and thermodynamic properties of 17 peptides with diverse secondary structure. Results obtained using the final form of the potential are presented for all these peptides. The same model, with unchanged parameters, is furthermore applied to a heterodimeric coiled-coil system, a mixed α/β protein and a three-helix-bundle protein, with very good results. The computational efficiency of the potential makes it possible to investigate the free-energy landscape of these 49–67-residue systems with high statistical accuracy, using only modest computational resources by todays standards. PACS Codes: 87.14.E-, 87.15.A-, 87.15.Cc

Simulation of Top7-CFr: a transient helix extension guides folding.

Protein structures often feature β-sheets in which adjacent β-strands have large sequence separation. How the folding process orchestrates the formation and correct arrangement of these strands is not comprehensively understood. Particularly challenging are proteins in which β-strands at the N and C termini are neighbors in a β-sheet. The N-terminal β-strand is synthesized early on, but it can not bind to the C terminus before the chain is fully synthesized. During this time, there is a danger that the β-strand at the N terminus interacts with nearby molecules, leading to potentially harmful aggregates of incompletely folded proteins. Simulations of the C-terminal fragment of Top7 show that this risk of misfolding and aggregation can be avoided by a “caching” mechanism that relies on the “chameleon” behavior of certain segments.

Distinct phases of free α‐synuclein—A Monte Carlo study

The α‐synuclein protein (αS), implicated in Parkinsons disease, shows conformational versatility. It aggregates into β‐sheet‐rich fibrils, occurs in helical membrane‐bound forms, is disordered as a free monomer, and has recently been suggested to have a folded helical tetramer as its main physiological form. Here, we use implicit solvent all‐atom Monte Carlo methods to explore the conformational ensemble sampled by the free αS monomer. We analyze secondary structure propensities, size, and topological properties and compare with existing experimental data. Our study suggests that free αS has two distinct phases. One phase has the expected disordered character. The other phase also shows large conformational variability. However, in this phase, the β‐strand content is substantial, and the backbone fold shows statistical similarities with that in αS fibrils. Presence of this phase is consistent with data from low‐temperature experiments. Conversion of disordered αS to this fibril‐like form requires the crossing of a rather large apparent free‐energy barrier. Proteins

SMMP v. 3.0 - Simulating proteins and protein interactions in Python and Fortran

We describe a revised and updated version of the program package SMMP. SMMP is an open-source FORTRAN package for molecular simulation of proteins within the standard geometry model. It is designed as a simple and inexpensive tool for researchers and students to become familiar with protein simulation techniques. SMMP 3.0 sports a revised API increasing its flexibility, an implementation of the Lund force field, multi-molecule simulations, a parallel implementation of the energy function, Python bindings, and more.

Equilibrium simulation of trp-cage in the presence of protein crowders

While steric crowders tend to stabilize globular proteins, it has been found that protein crowders can have an either stabilizing or destabilizing effect, where a destabilization may arise from nonspecific attractive interactions between the test protein and the crowders. Here, we use Monte Carlo replica-exchange methods to explore the equilibrium behavior of the miniprotein trp-cage in the presence of protein crowders. Our results suggest that the surrounding crowders prevent trp-cage from adopting its global native fold, while giving rise to a stabilization of its main secondary-structure element, an α-helix. With the crowding agent used (bovine pancreatic trypsin inhibitor), the trp-cage-crowder interactions are found to be specific, involving a few key residues, most of which are prolines. The effects of these crowders are contrasted with those of hard-sphere crowders.