Effect of external pulling forces on the length distribution of peptides

Background: The distribution of the length of a polypeptide, or that of the distance between any two of its atoms, is an important property as it can be analytically or numerically estimated for a number of polymer models. Importantly, it is directly measurable through a number of different experimental techniques. Length distributions can be straightforwardly assessed from molecular dynamics simulation; however, true convergence through full accurate coverage of the length range is difficult to achieve. Methods: The application of external constant force combined with the weighted-histogram analysis method (WHAM) is used to enhance sampling of unlikely 'long' or 'short' conformations and obtain the potential of mean force, while also collecting dynamic properties of the chain under variable tension. Results: We demonstrate the utility of constant force to enhance the sampling efficiency and obtain experimentally measurable quantities on a series of short peptides, including charge-rich sequences that are known to be highly helical but whose properties are distinct from those of helical peptides undergoing helix-coil transitions. Conclusions: Force-enhanced sampling enhances the range and accuracy of the length-based potential of mean force of the peptide, in particular those sequences that contain increased numbers of charged residues. General significance: This approach allows users to simultaneously probe the force-dependent behaviour of peptides directly, enhance the range and accuracy of the length-based PMF of the peptide and also test the convergence of simulations by comparing the overlap of PMF profiles from different constant forces. This article is part of a special issue entitled Recent developments of molecular dynamics.