CGCompiler: Automated coarse-grained molecule parameterization via noise-resistant mixed-variable optimization

13 June 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Coarse-grained force-fields (CG FF) such as the Martini model entail a predefined, fixed set of Lennard-Jones parameters (building blocks) to model virtually all possible non-bonded interactions between chemically relevant molecules. Owing to its universality and transferability, the building block coarse-grained approach has gained a tremendous popularity over the last decade. The parameterization of molecules can be highly complex and often involves the selection and fine tuning of a large number of parameters (e.g., bead types and bond lengths) to optimally match multiple relevant targets simultaneously. The parameterization of a molecule within the building block CG approach is a mixed-variable optimization problem: The non-bonded interactions are discrete variables whereas the bonded interactions are continuous variables. Here, we pioneer the utility of mixed-variable particle swarm optimization in automatically parameterizing molecules within the Martini 3 coarse-grained force-field by matching both structural (e.g., RDFs) as well as thermodynamic data (phase-transition temperatures). For sake of demonstration, we parameterize the linker of the lipid sphingomyelin. The important advantage of our approach is that both bonded- and non-bonded interactions are simultaneously optimized while conserving the search efficiency of vector guided particle swarm optimization (PSO) methods over other metaheuristic search methods such as genetic algorithms. In addition, we explore noise-mitigation strategies in matching the phase transition temperatures of lipid membranes, where nucleation and concomitant hysteresis introduces a dominant noise term within the objective function. We propose that noise-resistant mixed-variable PSO methods can both improve as well as automate parameterization of molecules within building block CG FFs, such as Martini.

Keywords

coarse grained
molecular dynamics
molecule parameterization
optimization
sphingomyelin
Martini

Supplementary materials

Title
Description
Actions
Title
Supporting Information
Description
Supporting Information including i) methods to estimate melting temperature, ii) additional plots regarding noise, iii) sphingomyelin-cholesterol 2d-center-of-mass radial distribution functions , and iv) DPSM Gromacs topology
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.