A modified Jarzynski free-energy estimator to eliminate non-conservative forces and its application in nanoparticle-membrane interactions

Computational methods to understand interactions in bio-complex systems are however limited to time-scales typically much shorter than in Nature. For example, on the nanoscale level, interactions between nanoparticles (NPs)/molecules/peptides and membranes are central in complex biomolecular processes such as membrane-coated NPs or cellular uptake. This can be remedied by the application of e.g. Jarzynski's equality where thermodynamic properties are extracted from non-equilibrium simulations. Although, the out of equilibrium work leads to non-conservative forces. We here propose a correction Pair Forces method, that removes these forces. Our proposed method is based on the calculation of pulling forces in backward and forward directions for the Jarzynski free-energy estimator using steered molecular dynamics simulation. Our results show that this leads to much improvement for NP-membrane translocation free energies. Although here we have demonstrated the application of the method in molecular dynamics simulation, it could be applied for experimental approaches

Electronic polarization effects on membrane translocation of anti-cancer drugs

Free-energy calculations are crucial for investigating biomolecular interactions. However, in theoretical studies, the neglect of electronic polarization can reduce predictive capabilities, specifically for free-energy calculations. To effectively mimick polarization, we explore a Charge Switching (CS) model, aiming to narrow the gap between computational and experimental results. The model requires quantum-level partial charge calculations of the molecule in different environments, combined with atomistic MD simulations. Studying three different anti-cancer drug molecules with three different phospholipid membranes, we show that the method significantly improves agreement with available experimental data. In contrast, using conventional fixed charge atomistic methods, qualitative discrepancies with experiments are observed, and we show that neglecting polarization may lead to an unphysical free energy sign inversion. While the CS method is here applied to anti-cancer drug-membrane translocation, it could be used more generally to study processes considering solvent effects