Influence of force field choice on the conformational landscape of rat and human islet amyloid polypeptide.

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein (IDP) whose abnormal aggregation into toxic soluble oligomers and insoluble amyloid fibrils is a pathological feature in type-2 diabetes. Rat IAPP (rIAPP) differs from hIAPP by only six amino acids yet has a reduced tendency to aggregate or form fibrils. The structures of the monomeric forms of IAPP are difficult to characterize due to their intrinsically disordered nature. Molecular dynamics simulations can provide a detailed characterization of the monomeric forms of rIAPP and hIAPP in near-physiological conditions. In this work, the conformational landscapes of rIAPP and hIAPP as a function of secondary structure content were predicted using well-tempered bias exchange metadynamics simulations. Several combinations of commonly used biomolecular force fields and water models were tested. The predicted conformational preferences of both rIAPP and hIAPP are typical of IDPs, exhibiting dominant random coil structures but showing a low propensity for transient α-helical conformations. Predicted nuclear magnetic resonance Cα chemical shifts reveal different preferences with each force field towards certain conformations, with AMBERff99SBnmr2/TIP4Pd showing the best agreement with the experiment. Comparisons of secondary structure content demonstrate residue-specific differences between hIAPP and rIAPP that may reflect their different aggregation propensities.