Genetic Algorithm Optimization of Force Field Parameters: Application to a Coarse-Grained Model of RNA

Determining force field parameters for molecular dynamics simulations of reduced models of biomolecules is a long, troublesome, and exhaustive process that is often performed manually. To improve this parametrization procedure we apply a continuous-space Genetic Algorithm (GA). GA is implemented to optimize parameters of a coarse-grained potential energy function of ribonucleic acid (RNA) molecules. The parameters obtained using GA are correctly reproducing the dynamical behavior of an RNA helix and other RNA tertiary motifs. Therefore, GA can be a useful tool for force field parametrization of the effective potentials in coarse-grained molecular models.