AMBER-DYES in AMBER: Implementation of fluorophore and linker parameters into AmberTools

Molecular dynamics (MD) simulations of explicit representations of fluorescent dyes attached via a linker to a protein allow, e.g., probing commonly used approximations for dye localization and/or orientation or modeling Förster resonance energy transfer. However, setting up and performing such MD simulations with the AMBER suite of biomolecular simulation programs has remained challenging due to the unavailability of an easy-to-use set of parameters within AMBER. Here, we adapted the AMBER-DYES parameter set derived by Graen et al. [J. Chem. Theory Comput. 10, 5505 (2014)] into "AMBER-DYES in AMBER" to generate a force field applicable within AMBER for commonly used fluorescent dyes and linkers attached to a protein. In particular, the computationally efficient graphics processing unit (GPU) implementation of the AMBER MD engine can now be exploited to overcome sampling issues of dye movements. The implementation is compatible with state-of-the-art force fields such as GAFF, GAFF2, ff99SB, ff14SB, lipid17, and GLYCAM_06j, which allows simulating post-translationally modified proteins and/or protein-ligand complexes and/or proteins in membrane environments. It is applicable with frequently used water models such as TIP3P, TIP4P, TIP4P-Ew, and OPC. For ease of use, a LEaP-based workflow was created, which allows attaching (multiple) dye/linker combinations to a protein prior to further system preparation steps. Following the parameter development described by Graen et al. [J. Chem. Theory Comput. 10, 5505 (2014)] and the adaptation steps described here, AMBER-DYES in AMBER can be extended by additional linkers and fluorescent molecules.