title: Going Backward: A Flexible Geometric Approach to Reverse Transformation from Coarse Grained to Atomistic Models.
abstract: The conversion of coarse-grained to atomistic models is an important step in obtaining insight about atomistic scale processes from coarse-grained simulations. For this process, called backmapping or reverse transformation, several tools are available, but these commonly require libraries of molecule fragments or they are linked to a specific software package. In addition, the methods are usually restricted to specific molecules and to a specific force field. Here, we present an alternative method, consisting of geometric projection and subsequent force-field based relaxation. This method is designed to be simple and flexible, and offers a generic solution for resolution transformation. For simple systems, the conversion only requires a list of particle correspondences on the two levels of resolution. For special cases, such as nondefault protonation states of amino acids and virtual sites, a target particle list can be specified. The mapping uses simple building blocks, which list the particles on the different levels of resolution. For conversion to higher resolution, the initial model is relaxed with several short cycles of energy minimization and position-restrained MD. The reconstruction of an atomistic backbone from a coarse-grained model is done using a new dedicated algorithm. The method is generic and can be used to map between any two particle based representations, provided that a mapping can be written. The focus of this work is on the coarse-grained MARTINI force field, for which mapping definitions are written to allow conversion to and from the higher-resolution force fields GROMOS, CHARMM, and AMBER, and to and from a simplified three-bead lipid model. Together, these offer the possibility to simulate mesoscopic membrane structures, to be transformed to MARTINI and subsequently to an atomistic model for investigation of detailed interactions. The method was tested on a set of systems ranging from a simple, single-component bilayer to a large protein-membrane-solvent complex. The results demonstrate the efficiency and the efficacy of the new approach.
authors.affiliation: Department of Biological Sciences, University of Calgary, 2500 University Drive NW , Calgary, Alberta, Canada T2N 1N4.|Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 7, 9747 AG Groningen, The Netherlands.|Computational Biology, Department of Biology, University of Erlangen-Nürnberg , Staudtstr. 5, 91058 Erlangen, Germany.
authors.forename: Tsjerk A
authors.initials: TA
authors.lastname: Wassenaar
authors.affiliation: Computational Biology, Department of Biology, University of Erlangen-Nürnberg , Staudtstr. 5, 91058 Erlangen, Germany.
authors.forename: Kristyna
authors.initials: K
authors.lastname: Pluhackova
authors.affiliation: Computational Biology, Department of Biology, University of Erlangen-Nürnberg , Staudtstr. 5, 91058 Erlangen, Germany.
authors.forename: Rainer A
authors.initials: RA
authors.lastname: Böckmann
authors.affiliation: Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 7, 9747 AG Groningen, The Netherlands.
authors.forename: Siewert J
authors.initials: SJ
authors.lastname: Marrink
authors.affiliation: Department of Biological Sciences, University of Calgary, 2500 University Drive NW , Calgary, Alberta, Canada T2N 1N4.
authors.forename: D Peter
authors.initials: DP
authors.lastname: Tieleman