Analyzing lipid distributions and curvature in molecular dynamics simulations of complex membranes.

We describe methods to analyze lipid distributions and curvature in membranes with complex lipid mixtures and embedded membrane proteins. We discuss issues involved in these analyses, available tools to calculate curvature preferences of lipids and proteins, and focus on tools developed in our group for visual analysis of lipid-protein interactions and the analysis of membrane curvature.

An implementation of the Martini coarse-grained force field in OpenMM.

We describe a complete implementation of Martini 2 and Martini 3 in the OpenMM molecular dynamics software package. Martini is a widely used coarse-grained force field with applications in biomolecular simulation, materials, and broader areas of chemistry. It is implemented as a force field but makes extensive use of facilities unique to the GROMACS software, including virtual sites and bonded terms that are not commonly used in standard atomistic force fields. OpenMM is a flexible molecular dynamics package widely used for methods development and is competitive in speed on GPUs with other commonly used packages. OpenMM has facilities to easily implement new force field terms, external forces and fields, and other nonstandard features, which we use to implement all force field terms used in Martini 2 and Martini 3. This allows Martini simulations, starting with GROMACS topology files that are processed by custom scripts, with all the added flexibility of OpenMM. We provide a GitHub repository with test cases, compare accuracy and performance between GROMACS and OpenMM, and discuss the limitations of our implementation in terms of direct comparison with GROMACS. We describe a use case that implements the Modeling Employing Limited Data method to apply experimental constraints in a Martini simulation to efficiently determine the structure of a protein complex. We also discuss issues and a potential solution with the Martini 2 topology for cholesterol.

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers.

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial ß-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

Molecular modelling of TLR agonist Pam3CSK4 entrapment in PLA nanoparticles as a tool to explain loading efficiency and functionality.

Designing potent and safe-of-use therapies against cancers and infections remains challenging despite the emergence of novel molecule classes like checkpoint inhibitors or Toll-Like-Receptor ligands. The latest therapeutic perspectives under development for immune modulator administration exploits vectorization, and biodegradable delivery systems are one of the most promising vehicles. Nanoparticles based on Poly (D,L) Lactic Acid (PLA) as polymer for formulation are widely investigated due to its bioresorbable, biocompatible and low immunogen properties. We propose a PLA-based nanoparticle delivery system to vectorize Pam3CSK4, a lipopeptide TLR1/2 ligand and a potent activator of the proinflammatory transcription factor NF-κB that shows a self-assembling behavior from 30 µg/mL onwards. We demonstrate successful encapsulation of Pam3CSK4 in PLA nanoparticles by nanoprecipitation in a 40-180 µg/mL concentration range, with 99% of entrapment efficiency. By molecular modelling, we characterize drug/carrier interactions and conclude that Pam3CSK4 forms clusters onto the nanoparticle and is not encapsulated into the hydrophobic core. In silico predictions provide nanoprecipitation optimization and the mechanistic understanding of the particle dynamics. The loaded-Pam3CSK4 maintains bioactivity on TLR2, confirmed by in vitro experiments using reporter cell line HEK-Blue hTLR2. Our presented data and results are convincing evidence that Pam3CSK4-loaded in PLA nanoparticles represent a promising immune modulating system.