Abstract:Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 μs (five replicates of 200 ns, 300 ns, 1 μs, or 2 μs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, β amyloid (Aβ)-40, Aβ-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of β sheets for β-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of β-amyloid 42 starting from two different initial structures (Aβ421Z0Q and Aβ42model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 μs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00762?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00762</a>.Predicted measurement comparison (Table S1); secondary chemical shift conversion (Figure S1); time evolution of the cumulative Rg and the end-to-end distance (Figure S2); comparison of secondary chemical shifts between MD simulations and NMR experiments (Figure S3); cluster analyses (Figure S4); calculated NMR observables in simulations of Ab40 starting from the PDB structure and the extended structure (Figure S5); secondary structure propensity analysis of the trajectories (Figure S6); RMSD analyses of all of the trajectories under the tested force fields (Figure S7) (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00762/suppl_file/ci0c00762_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.

All-atom Molecular Dynamics Simulations of the Projection Domain of the Intrinsically Disordered htau40 Protein

Probing Conformational Dynamics of Tau Protein by Hydrogen/Deuterium Exchange Mass Spectrometry

Combining molecular dynamics simulations with small-angle X-ray and neutron scattering data to study multi-domain proteins in solution.

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau

Dancing with Strings: the Conformational Dynamics of Vqixxk Motifs Within Tau Protein in Monomer, Fibril and Hyper-Phosphorylated Filament States

Microsecond Dynamics During the Binding-induced Folding of an Intrinsically Disordered Protein

Conformational Dependence of Chemical Shifts in the Proline Rich Region of TAU Protein

nP-collabs: Investigating counterion mediated bridges in the multiply phosphorylated tau-R2 repeat

Residual structure in the repeat domain of tau: echoes of microtubule binding and paired helical filament formation

Molecular Dynamics Simulation of Protein Conformation Transition in Solution

n P-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat

Kinetic Ensemble of Tau Protein through the Markov State Model and Deep Learning Analysis

Resolution Exchange with Tunneling for Enhanced Sampling of Protein Landscapes

Structural Insight into Tau Protein'S Paradox of Intrinsically Disordered Behavior, Self-Acetylation Activity, and Aggregation

Translational diffusion of hydration water correlates with functional motions in folded and intrinsically disordered proteins

SOP-MULTI: A self-organized polymer based coarse-grained model for multi-domain and intrinsically disordered proteins with conformation ensemble consistent with experimental scattering data

Construction of an Intermediate-Resolution Lattice Model and Re-Examination of the Helix-Coil Transition: a Dynamic Monte Carlo Simulation.

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

All-atom Replica Exchange Molecular Simulation of Protein BBL

Temperature-Dependent Folding Pathways Of Pin1 Ww Domain: An All-Atom Molecular Dynamics Simulation Of A G(O)Over-Bar Model

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins