Fine structures of intrinsically disordered proteins

Abstract

We report simulation studies of 33 single intrinsically disordered proteins (IDPs) using coarse-grained bead-spring models where interactions among different amino acids are introduced through a hydropathy matrix and additional screened Coulomb interaction for the charged amino acid beads. Our simulation studies of two different hydropathy scales (HPS1, HPS2) [Dignon et al., PLoS Comput. Biol. 14, e1005941 (2018); Tesei et al. Proc. Natl. Acad. Sci. U. S. A. 118, e2111696118 (2021)] and the comparison with the existing experimental data indicate an optimal interaction parameter ϵ = 0.1 and 0.2 kcal/mol for the HPS1 and HPS2 hydropathy scales. We use these best-fit parameters to investigate both the universal aspects as well as the fine structures of the individual IDPs by introducing additional characteristics. (i) First, we investigate the polymer-specific scaling relations of the IDPs in comparison to the universal scaling relations [Bair et al., J. Chem. Phys. 158, 204902 (2023)] for the homopolymers. By studying the scaled end-to-end distances ⟨RN2⟩/(2Lℓp) and the scaled transverse fluctuations l̃⊥2=⟨l⊥2⟩/L, we demonstrate that IDPs are broadly characterized with a Flory exponent of ν ≃ 0.56 with the conclusion that conformations of the IDPs interpolate between Gaussian and self-avoiding random walk chains. Then, we introduce (ii) Wilson charge index (W) that captures the essential features of charge interactions and distribution in the sequence space and (iii) a skewness index (S) that captures the finer shape variation of the gyration radii distributions as a function of the net charge per residue and charge asymmetry parameter. Finally, our study of the (iv) variation of ⟨Rg⟩ as a function of salt concentration provides another important metric to bring out finer characteristics of the IDPs, which may carry relevant information for the origin of life.