Lipidic moieties were described by common Martini parameters. RLs CG topologies (not available in the Martini FF) were built starting from those of monogalactosyldiacylglycerol (MGDG) and digalactosyldiacylglycerol (DGDG) in their anionic form (López et al., 2013 (
link); van Eerden et al., 2015 (
link)). Out of many different models available, differing in the mapping, Martini bead types, rhamnose/s bead bonds, dihedral angles as well as constraint and/or exclusion definitions, we selected the one shown in
Figure 2 on the basis of the comparison of results obtained by all-atom (AA) and CG simulations. In particular, a RL (mono- or di-RL) was simulated in both water (10 ns) and POPC membrane (128 total lipids, 500 ns) using both C36 (2 fs time steps) and Martini 2 (20 fs time step) FFs.
First, CG distributions of bonds, angles and dihedral angles were compared to the atomistic ones (for the comparison we used the centre of mass of atoms corresponding to each CG bead). The choice of the best model required several tests until a good reproduction of CG bonds, angles and dihedral angles was achieved (see
Supplementary Figures S1, S2 for the optimised mono- and di-RLs topologies, shown in
Supplementary Tables S1, S2).
To ensure that the chosen CG topologies could also reproduce membrane properties, we monitored both area per lipid and membrane thickness of (a) RL-containing POPC membrane as well as for (b) a more complex model containing POPC/POPG/ergosterol (53/23/25) and 1 RL/25 lipids (40/60, mono- and di-RL) (see
Supplementary Figure S3). Both models allowed us to compare the results with previous works from our group (Monnier et al., 2019 (
link)). These latter models were run with 256 lipids, 4 mono-RLs and 6 di-RLs for 500 ns (AA) and with 864 lipids, 9 mono-RLs and 27 di-RLs for 10 µs (CG).
For model (a), membrane properties were nicely reproduced with stable trajectories (20 fs time step). Good results were also achieved for model (b) (see
Supplementary Figure S3) when reducing the time step (14 fs), as previously reported for other complex glycolipids (López et al., 2013 (
link)).
The final parameters used for mono- and di-RLs topologies are listed in
Supplementary Tables S1, S2. They correspond to those of the model which gave the best results in terms of system stability, bond, angle and dihedral angle distributions (
Supplementary Figure S1), and structural membrane properties (
Supplementary Figure S3). Mappings and Martini bead assignments are shown in
Figure 2A.
Rodríguez-Moraga N., Ramos-Martín F., Buchoux S., Rippa S., D’Amelio N, & Sarazin C. (2023). The effect of rhamnolipids on fungal membrane models as described by their interactions with phospholipids and sterols: An in silico study. Frontiers in Chemistry, 11, 1124129.