Production of MD trajectories: All-atom MD simulations were performed on a local High Performance Cluster (HPC) called Chemfarm with nodes of 2 GPUs and 36 CPUs using the CHARMM-36 force field of GROMACS 2020.2 package. All simulations comprise trajectories of 50ns, with atomic coordinates recorded every 2ps. Two sets of triplicates were performed for binary systems (containing two types of lipids), whereas one set of triplicates was performed for each pure system (containing a single lipid-type). For each replica, new initial random velocities were generated, and for each set of replicas new initial atom positions were set as well. For the binary systems, five different concentration ratios were performed: The following MD protocols were used: the integration time step was 2fs; the isobaric-isothermal (NPT) ensemble was employed; the pressure was set to 1 bar using isotropic coupling to the Parrinello-Rahman barostat with a time constant of 5ps and an isothermal compressibility of 4.5*10 -5 bar -1 ; the temperature was kept constant at 300K using the Nosé-Hoover thermostat with a time constant of 1ps. Bonds with hydrogen atoms were constrained using the Linear Constraint Solver (LINCS).