This paper presents a new model for polymersomes developed using a poly(ethylene oxide)-poly(butadiene) diblock copolymer bilayer. The model is based on a coarse-grained approach using the MARTINI force field. Since no MARTINI parameters exist for poly(butadiene), we have refined these parameters using quantum mechanical calculations and molecular dynamics simulations. The model has been validated using extensive molecular dynamics simulations in systems with several hundred polymer units and reaching up to 6 μs. These simulations show that the copolymer coarse grain model self-assemble into bilayers and that NPT and NPT ensemble runs reproduce key structural and mechanical experimental properties for different copolymer length chains with a similar hydrophilic weight fraction.

No abstract

No abstract

In this work, we developed a coarse-grained model of sumatriptan suitable for extensive molecular dynamics simulations. First, we confirmed the interfacial distribution of this drug in bilayers through cryogenic transmission electron microscopy and small-angle X-ray scattering techniques, as was predicted by our previous atomistic simulations. Based on these simulations, we developed a coarse-grained model for sumatriptan able to reproduce its overall molecular behavior, captured by atomistic simulations and experiments. We then tested the sumatriptan model in a micellar environment along with experimental characterization of sumatriptan-loaded micelles. The simulation results showed good agreement with photon correlation spectroscopy and electrophoretic mobility experiments performed in this work. The particle size of the obtained micelles was comparable with the simulated ones; meanwhile, zeta-potential results suggest adsorption of the drug on the micellar surface. This model is a step forward in the search for a suitable drug-delivery system for sumatriptan.

In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.

Molecular dynamics simulation is a very powerful tool to understand biomolecular processes. In this chapter, we go over different applications of this methodology to drug delivery systems (DDS) carried out in the group. DDS-a formulation or a device that enables the introduction of a therapeutic substance in the body and improves its efficacy and safety by controlling the rate, time, and place of release of drugs-are an important component of drug development and therapeutics. Biocompatible nanoparticles are materials in the nanoscale that emerged as important players, improving efficacy of approved drugs, for example. The molecular understanding of the encapsulation process could be very helpful to guide the nanocarrier for a specific system. Here we discuss different applications of drug delivery carriers, such as liposomes, polymeric micelles, and polymersomes using atomistic and coarse grain (CG) molecular dynamics simulations.

Connexinophaties are a collective of diseases related to connexin channels and hemichannels. In particular many Cx26 alterations are strongly associated to human deafness. Calcium plays an important role on this structures regulation. Here, using calcium as a probe, extensive atomistic Molecular Dynamics simulations were performed on the Cx26 hemichannel embedded in a lipid bilayer. Exploring different initial conditions and calcium concentration, simulation reached ~4 μs. Several analysis were carried out in order to reveal the calcium distribution and localization, such as electron density profiles, density maps and distance time evolution, which is directly associated to the interaction energy. Specific amino acid interactions with calcium and their stability were capture within this context. Few of these sites such as, GLU42, GLU47, GLY45 and ASP50, were already suggested in the literature. Besides, we identified novel calcium biding sites: ASP2, ASP117, ASP159, GLU114, GLU119, GLU120 and VAL226. To the best of our knowledge, this is the first time that these sites are reported within this context. Furthermore, since various pathologies involving the Cx26 hemichannel are associated with pathogenic variants in the corresponding CJB2 gene, using ClinVar, we were able to spatially associate the 3D positions of the identified calcium binding sites within the framework of this work with reported pathogenic variants in the CJB2 gene. This study presents a first step on finding associations between molecular features and pathological variants of the Cx26 hemichannel.

In this work, we analyzed the behavior of Pluronic F127 through molecular dynamics simulations at the coarse-grain level, focusing on the micellar and lamellar phases. To this aim, two initial polymer conformations were considered, S-shape and U-shape, for both simulated phases. Through the simulations, we were able to examine the structural and mechanical properties that are difficult to access through experiments. Since no transition between S and U shapes was observed in our simulations, we inferred that all single co-polymers had memory of their initial configuration. Nevertheless, most copolymers had a more complex amorphous structure, where hydrophilic beads were part of the lamellar-like core. Finally, an overall comparison of the micellar a lamellar phases showed that the lamellar thickness was in the same order of magnitude as the micelle diameter (approx. 30 nm). Therefore, high micelle concentration could lead to lamellar formation. With this new information, we could understand lamellae as orderly packed micelles.