Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.

Central European Institute of Technology Masaryk University Brno Czech Republic

National Centre for Biomolecular Research Faculty of Science Masaryk University Brno Czech Republic

National Centre for Biomolecular Research Faculty of Science Masaryk University Brno Czech Republic; Central European Institute of Technology Masaryk University Brno Czech Republic

See more in PubMed

Uversky V.N. Natively unfolded proteins: a point where biology waits for physics. Protein Sci. 2002;11:739–756. PubMed PMC

Dunker A.K., Brown C.J., Obradović Z. Intrinsic disorder and protein function. Biochemistry. 2002;41:6573–6582. PubMed

Tompa P. Unstructural biology coming of age. Curr. Opin. Struct. Biol. 2011;21:419–425. PubMed

Chi S.-W., Kim D.-H., Han K.H. Pre-structured motifs in the natively unstructured preS1 surface antigen of hepatitis B virus. Protein Sci. 2007;16:2108–2117. PubMed PMC

Fuxreiter M., Simon I., Tompa P. Preformed structural elements feature in partner recognition by intrinsically unstructured proteins. J. Mol. Biol. 2004;338:1015–1026. PubMed

Vacic V., Oldfield C.J., Dunker A.K. Characterization of molecular recognition features, MoRFs, and their binding partners. J. Proteome Res. 2007;6:2351–2366. PubMed PMC

Nováček J., Židek L., Sklenář V. Toward optimal-resolution NMR of intrinsically disordered proteins. J. Magn. Reson. 2014;241:41–52. PubMed

Nowakowski M., Saxena S., Koźmiński W. Applications of high dimensionality experiments to biomolecular NMR. Prog. Nucl. Magn. Reson. Spectrosc. 2015;90–91:49–73. PubMed

Papoian G.A. Proteins with weakly funneled energy landscapes challenge the classical structure-function paradigm. Proc. Natl. Acad. Sci. USA. 2008;105:14237–14238. PubMed PMC

Rabatinová A., Šanderová H., Krásný L. The δ subunit of RNA polymerase is required for rapid changes in gene expression and competitive fitness of the cell. J. Bacteriol. 2013;195:2603–2611. PubMed PMC

Papoušková V., Kadeřávek P., Žídek L. Structural study of the partially disordered full-length δ subunit of RNA polymerase from Bacillus subtilis. ChemBioChem. 2013;14:1772–1779. PubMed

Nagatsu T., Levitt M., Udenfriend S. Tyrosine hydroxylase. The initial step in norepinephrine biosynthesis. J. Biol. Chem. 1964;239:2910–2917. PubMed

Molinoff P.B., Axelrod J. Biochemistry of catecholamines. Annu. Rev. Biochem. 1971;40:465–500. PubMed

Louša P., Nedozrálová H., Hritz J. Phosphorylation of the regulatory domain of human tyrosine hydroxylase 1 monitored using non-uniformly sampled NMR. Biophys. Chem. 2017;223:25–29. PubMed

Jansen S., Melková K., Žídek L. Quantitative mapping of microtubule-associated protein 2c (MAP2c) phosphorylation and regulatory protein 14-3-3ζ-binding sites reveals key differences between MAP2c and its homolog Tau. J. Biol. Chem. 2017;292:6715–6727. PubMed PMC

Melková K., Zapletal V., Žídek L. Functionally specific binding regions of microtubule-associated protein 2c exhibit distinct conformations and dynamics. J. Biol. Chem. 2018;293:13297–13309. PubMed PMC

Melková K., Zapletal V., Žídek L. Structure and functions of microtubule associated proteins Tau and MAP2c: similarities and differences. Biomolecules. 2019;9:E105. PubMed PMC

Motáčková V., Nováček J., Sklenář V. Strategy for complete NMR assignment of disordered proteins with highly repetitive sequences based on resolution-enhanced 5D experiments. J. Biomol. NMR. 2010;48:169–177. PubMed PMC

Nováček J., Janda L., Sklenář V. Efficient protocol for backbone and side-chain assignments of large, intrinsically disordered proteins: transient secondary structure analysis of 49.2 kDa microtubule associated protein 2c. J. Biomol. NMR. 2013;56:291–301. PubMed

Ottiger M., Delaglio F., Bax A. Measurement of J and dipolar couplings from simplified two-dimensional NMR spectra. J. Magn. Reson. 1998;131:373–378. PubMed

Srb P., Nováček J., Žídek L. Triple resonance 15Ν NMR relaxation experiments for studies of intrinsically disordered proteins. J. Biomol. NMR. 2017;69:133–146. PubMed

Korzhnev D.M., Billeter M., Orekhov V.Y. NMR studies of Brownian tumbling and internal motions in proteins. Prog. Nucl. Magn. Reson. Spectrosc. 2001;38:197–266.

Efron B. Bootstrap methods: another look at the jackknife. Ann. Stat. 1979;7:1–26.

Petoukhov M.V., Franke D., Svergun D.I. New developments in the ATSAS program package for small-angle scattering data analysis. J. Appl. Cryst. 2012;45:342–350. PubMed PMC

Riback J.A., Bowman M.A., Sosnick T.R. Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water. Science. 2017;358:238–241. PubMed PMC

Lindorff-Larsen K., Piana S., Shaw D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010;78:1950–1958. PubMed PMC

Piana S., Lindorff-Larsen K., Shaw D.E. How robust are protein folding simulations with respect to force field parameterization? Biophys. J. 2011;100:L47–L49. PubMed PMC

Huang J., Rauscher S., MacKerell A.D., Jr. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods. 2017;14:71–73. PubMed PMC

Jorgensen W.L. Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water. J. Am. Chem. Soc. 1981;103:335–340.

MacKerell A.D., Bashford D., Karplus M. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B. 1998;102:3586–3616. PubMed

Piana S., Donchev A.G., Shaw D.E. Water dispersion interactions strongly influence simulated structural properties of disordered protein states. J. Phys. Chem. B. 2015;119:5113–5123. PubMed

Hess B., Bekker H., Fraaije J.G.E.M. LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 1997;18:1463–1472.

Essmann U., Perera L., Pedersen L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995;103:8577–8593.

Berendsen H.J.C., Postma J.P.M., Haak J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984;81:3684–3690.

Bussi G., Donadio D., Parrinello M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007;126:014101. PubMed

Parrinello M., Rahman A. Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 1981;52:7182–7190.

Shen Y., Bax A. SPARTA+: a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural network. J. Biomol. NMR. 2010;48:13–22. PubMed PMC

Nielsen J.T., Mulder F.A.A. POTENCI: prediction of temperature, neighbor and pH-corrected chemical shifts for intrinsically disordered proteins. J. Biomol. NMR. 2018;70:141–165. PubMed

Zweckstetter M. NMR: prediction of molecular alignment from structure using the PALES software. Nat. Protoc. 2008;3:679–690. PubMed

Nodet G., Salmon L., Blackledge M. Quantitative description of backbone conformational sampling of unfolded proteins at amino acid resolution from NMR residual dipolar couplings. J. Am. Chem. Soc. 2009;131:17908–17918. PubMed

Salmon L., Nodet G., Blackledge M. NMR characterization of long-range order in intrinsically disordered proteins. J. Am. Chem. Soc. 2010;132:8407–8418. PubMed

Sezer D., Freed J.H., Roux B. Simulating electron spin resonance spectra of nitroxide spin labels from molecular dynamics and stochastic trajectories. J. Chem. Phys. 2008;128:165106. PubMed PMC

Salvi N., Abyzov A., Blackledge M. Multi-timescale dynamics in intrinsically disordered proteins from NMR relaxation and molecular simulation. J. Phys. Chem. Lett. 2016;7:2483–2489. PubMed

Urbańczyk M., Bernin D., Kazimierczuk K. Iterative thresholding algorithm for multiexponential decay applied to PGSE NMR data. Anal. Chem. 2013;85:1828–1833. PubMed

Robustelli P., Piana S., Shaw D.E. Developing a molecular dynamics force field for both folded and disordered protein states. Proc. Natl. Acad. Sci. USA. 2018;115:E4758–E4766. PubMed PMC

Bax A. Weak alignment offers new NMR opportunities to study protein structure and dynamics. Protein Sci. 2003;12:1–16. PubMed PMC

Fraczkiewicz R., Braun W. Exact and efficient analytical calculation of the accessible surface areas and their gradients for macromolecules. J. Comput. Chem. 1998;19:319–333.

Song D., Luo R., Chen H.-F. The IDP-specific force field ff14IDPSFF improves the conformer sampling of intrinsically disordered proteins. J. Chem. Inf. Model. 2017;57:1166–1178. PubMed PMC

Streamlining NMR Chemical Shift Predictions for Intrinsically Disordered Proteins: Design of Ensembles with Dimensionality Reduction and Clustering

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

Convergent views on disordered protein dynamics from NMR and computational approaches

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins