The amyloid-β (Aβ) monomer, an intrinsically disordered peptide, is produced by the cleavage of the amyloid precursor protein, leading to Aβ-40 and Aβ-42 as major products. These two isoforms generate pathological aggregates, whose accumulation correlates with Alzheimer’s disease (AD). Experiments have shown that even though the natural abundance of Aβ-42 is smaller than that for Aβ-40, the Aβ-42 is more aggregation-prone compared to Aβ-40. Moreover, several single-point mutations are associated with early onset forms of AD. This work analyzes coarse-grained associative-memory, water-mediated, structure and energy model (AWSEM) simulations of normal Aβ-40 and Aβ-42 monomers, along with six single-point mutations associated with early onset disease. We analyzed the simulations using the energy landscape visualization method (ELViM), a reaction-coordinate-free approach suited to explore the frustrated energy landscapes of intrinsically disordered proteins. ELViM is shown to distinguish the monomer ensembles of variants that rapidly form fibers from those that do not form fibers as readily. It also delineates the amino acid contacts characterizing each ensemble. The results shed light on the potential of ELViM to probe intrinsically disordered proteins.

Ankyrin (ANK) repeat proteins are coded by tandem occurrences of patterns with around 33 amino acids. They often mediate protein–protein interactions in a diversity of biological systems. These proteins have an elongated non-globular shape and often display complex folding mechanisms. This work investigates the energy landscape of representative proteins of this class made up of 3, 4 and 6 ANK repeats using the energy-landscape visualisation method (ELViM). By combining biased and unbiased coarse-grained molecular dynamics AWSEM simulations that sample conformations along the folding trajectories with the ELViM structure-based phase space, one finds a three-dimensional representation of the globally funnelled energy surface. In this representation, it is possible to delineate distinct folding pathways. We show that ELViMs can project, in a natural way, the intricacies of the highly dimensional energy landscapes encoded by the highly symmetric ankyrin repeat proteins into useful low-dimensional representations. These projections can discriminate between multiplicities of specific parallel folding mechanisms that otherwise can be hidden in oversimplified depictions.

Rational design of proteins via mutagenesis is crucial for several biotechnological applications. A significant challenge of the computational strategies used to predict optimized mutations is to understand the influence of each amino acid during the folding process. In the present work, chymotrypsin inhibitor 2 (CI2) and several of its designed mutants have been simulated using a non-native hydrophobic and electrostatic potential as a structure-based Cα model. Through these simulations, we could identify the most critical folding stage to accelerate CI2 and also the charged residues responsible for providing its thermostability. The replacement of ionizable residues for hydrophobic ones tended to promote the formation of the CI2 secondary structure in the early transition state, which speeds up folding. However, this same replacement destabilized the native structure, and there was a decrease in the protein thermostability. Such a simple method proved to be capable of providing valuable information about thermodynamics and kinetics of CI2 and its mutations, thus being a fast alternative to the study of rational protein design.

The Cu/Zn Human Superoxide Dismutase (SOD1) is a dimeric metalloenzyme whose genetic mutations are directly related to amyotrophic lateral sclerosis (ALS), so understanding its folding mechanism is of fundamental importance. Currently, the SOD1 dimer formation is studied via molecular dynamics simulations using a simplified structure-based model and an all-atom model. Results from the simplified model reveal a mechanism dependent on distances between monomers, which are limited by constraints to mimic concentration dependence. The stability of intermediates (during the int state) is significantly affected by this distance, as well as by the presence of two folded monomers prior to dimer formation. The kinetics of interface formation are also highly dependent on the separation distance. The folding temperature of the dimer is about 4.2% higher than that of the monomer, a value not too different from experimental data. All-atom simulations on the apo dimer give binding free energy between monomers similar to experimental values. An intermediate state is evident for the apo form at a separation distance between monomers slightly larger than the native distance which has little formed interface between monomers. We have shown that this intermediate is stabilized by non-native intra- and intercontacts.