13Molecular dynamics (MD) has greatly contributed to understanding and predicting the way 14 proteins fold. However, the time-scale and complexity of folding are not accessible via 15 classical MD. Furthermore, efficient folding pipelines involving enhanced MD techniques 16 are not routinely accessible. We aimed to determine whether perturbing the electrostatic 17 component of the MD force field can help expedite folding simulations. We developed 18 charge-perturbation dynamics (CPD), an MD-based simulation approach that involves 19 periodically perturbing the atomic charges to values non-native to the MD force field. CPD 20 obtains suitable sampling via multiple iterations in which a classical MD segment (with 21 native charges) is followed by a very short segment of perturbed MD (using the same force 22 field and conditions, but with non-native charges); subsequently, partially folded 23 intermediates are refined via a longer segment of classical MD. Among the partially folded 24 structures from low-energy regions of the free-energy landscape sampled, the lowest-25 energy conformer with high root-mean-square deviation to the starting structure and low 26 radius of gyration is defined as the folded structure. Upon benchmark testing, we found 27 that medium-length peptides such as an alanine-based pentadecapeptide, an amyloid-β 28 peptide, and the tryptophan-cage mini-protein can fold starting from their extended linear 29 structure in under 45 ns of CPD (total simulation time), versus over 100 ns of classical 30 MD. CPD not only achieved folding close to the desired conformation but also sampled 31 key intermediates along the folding pathway without prior knowledge of the folding 32 mechanism or final folded structure. Our findings confirmed that perturbing the 33 electrostatic component of the classical MD force field can help expedite folding 34 3 simulations without changing the MD algorithm or using expensive computing 35 architectures. CPD can be employed to probe the folding dynamics of known, putative, or 36 planned peptides, as well as to improve sampling in more advanced simulations or to guide 37 further experiments. 38 39 4 Author summary 40 Folding represents the process by which proteins assemble into biologically active 41 conformations. While computational techniques such as molecular dynamics (MD) have 42 provided invaluable insight into protein folding, efficient folding pipelines are not 43 routinely accessible. In MD, the behavior of the studied molecule is simulated under the 44 concerted action of multiple forces described by mathematical functions employing 45 optimized parameters. Using non-native parameters effectively perturbs the MD force 46field. We show that this can be exploited to help expedite folding simulations. Specifically,
47we developed charge-perturbation dynamics (CPD), an MD-based simulation approach 48 that involves periodically perturbing the force field by using non-native atomic charges.
49For folding medium-length peptides such as the tryptophan-cage mini-protein s...