Molecular dynamics simulations on the conformational transitions from the GA 98 (GA 88) to GB 98 (GB 88) proteins

Exploring the Effect of Enhanced Sampling on Protein Stability Prediction

Changes in protein stability due to side-chain mutations are evaluated by alchemical free-energy calculations based on classical molecular dynamics (MD) simulations in explicit solvent using the GROMOS force field. Three proteins of different complexity with a total number of 93 single-point mutations are analyzed, and the relative free-energy differences are discussed with respect to configurational sampling and (dis)agreement with experimental data. For the smallest protein studied, a 34-residue WW domain, the starting structure dependence of the alchemical free-energy changes, is discussed in detail. Deviations from previous simulations for the two other proteins are shown to result from insufficient sampling in the earlier studies. Hamiltonian replica exchange in combination with multiple starting structures and sufficient sampling time of more than 100 ns per intermediate alchemical state is required in some cases to reach convergence.

Insights into the aggregation mechanism of RNA recognition motif domains in TDP-43: a theoretical exploration

The transactive response DNA-binding protein 43 (TDP-43) is associated with several diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) due to pathogenic aggregations. In this work, we examined the dimer, tetramer and hexamer models built from the RRM domains of TDP-43 using molecular dynamics simulations in combination with the protein-protein docking. Our results showed that the formations of the dimer models are mainly achieved by the interactions of the RRM1 domains. The parallel β-sheet layers between the RRM1 domains provide most of the binding sites in these oligomer models, and thus play an important role in the aggregation process. The approaching of the parallel β-sheet layers from small oligomer models gradually expand to large ones through the allosteric communication between the α1/α2 helices of the RRM1 domains, which maintains the binding affinities and interactions in the larger oligomer models. Using the repeatable-superimposing method based on the tetramer models, we proposed a new aggregation mechanism of RRM domains in TDP-43, which could well characterize the formation of the large aggregation models with the repeated, helical and rope-like structures. These new insights help to understand the amyloid-like aggregation phenomena of TDP-43 protein in ALS and FTLD diseases.