Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer.

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

Depletion of Mannose Receptor-Positive Tumor-associated Macrophages via a Peptide-targeted Star-shaped Polyglutamate Inhibits Breast Cancer Progression in Mice.

Although many studies have explored the depletion of tumor-associated macrophages (TAM) as a therapeutic strategy for solid tumors, currently available compounds suffer from poor efficacy and dose-limiting side effects. Here, we developed a novel TAM-depleting agent ("OximUNO") that specifically targets CD206+ TAMs and demonstrated efficacy in a triple-negative breast cancer (TNBC) mouse model. OximUNO comprises a star-shaped polyglutamate (St-PGA) decorated with the CD206-targeting peptide mUNO that carries the chemotherapeutic drug doxorubicin (DOX). In the TNBC model, a fluorescently labeled mUNO-decorated St-PGA homed to CD206+ TAMs within primary lesions and metastases. OximUNO exhibited no acute liver or kidney toxicity in vivo. Treatment with OximUNO reduced the progression of primary tumor lesions and pulmonary metastases, significantly diminished the number of CD206+ TAMs and increased the CD8/FOXP3 expression ratio (indicating immunomodulation). Our findings suggest the potential benefit of OximUNO as a TAM-depleting agent for TNBC treatment. Importantly, our studies also represent a novel design of a peptide-targeted St-PGA as a targeted therapeutic nanoconjugate. Significance: A peptide-targeted nanoformulation of DOX exclusively eliminates mannose receptor+ TAMs in breast cancer models, generating response without off-target effects (a drawback of many TAM-depleting agents under clinical study).

In Silico Selection and Evaluation of Pugnins with Antibacterial and Anticancer Activity Using Skin Transcriptome of Treefrog (Boana pugnax).

In order to combat bacterial and cancer resistance, we identified peptides (pugnins) with dual antibacterial l-anticancer activity from the Boana pugnax (B. pugnax) skin transcriptome through in silico analysis. Pugnins A and B were selected owing to their high similarity to the DS4.3 peptide, which served as a template for their alignment to the B. pugnax transcriptome, as well as their function as part of a voltage-dependent potassium channel protein. The secondary peptide structure stability in aqueous medium was evaluated as well, and after interaction with the Escherichia coli (E. coli) membrane model using molecular dynamics. These pugnins were synthesized via solid-phase synthesis strategy and verified by Reverse phase high-performance liquid chromatography (RP-HPLC) and mass spectrometry. Subsequently, their alpha-helix structure was determined by circular dichroism, after which antibacterial tests were then performed to evaluate their antimicrobial activity. Cytotoxicity tests against cancer cells also showed selectivity of pugnin A toward breast cancer (MFC7) cells, and pugnin B toward prostate cancer (PC3) cells. Alternatively, flow cytometry revealed necrotic cell damage with a major cytotoxic effect on human keratinocytes (HaCaT) control cells. Therefore, the pugnins found in the transcriptome of B. pugnax present dual antibacterial-anticancer activity with reduced selectivity to normal eukaryotic cells.

A Global Review on Short Peptides: Frontiers and Perspectives.

Peptides are fragments of proteins that carry out biological functions. They act as signaling entities via all domains of life and interfere with protein-protein interactions, which are indispensable in bio-processes. Short peptides include fundamental molecular information for a prelude to the symphony of life. They have aroused considerable interest due to their unique features and great promise in innovative bio-therapies. This work focusing on the current state-of-the-art short peptide-based therapeutical developments is the first global review written by researchers from all continents, as a celebration of 100 years of peptide therapeutics since the commencement of insulin therapy in the 1920s. Peptide "drugs" initially played only the role of hormone analogs to balance disorders. Nowadays, they achieve numerous biomedical tasks, can cross membranes, or reach intracellular targets. The role of peptides in bio-processes can hardly be mimicked by other chemical substances. The article is divided into independent sections, which are related to either the progress in short peptide-based theranostics or the problems posing challenge to bio-medicine. In particular, the SWOT analysis of short peptides, their relevance in therapies of diverse diseases, improvements in (bio)synthesis platforms, advanced nano-supramolecular technologies, aptamers, altered peptide ligands and in silico methodologies to overcome peptide limitations, modern smart bio-functional materials, vaccines, and drug/gene-targeted delivery systems are discussed.

Differential Stability of Aurein 1.2 Pores in Model Membranes of Two Probiotic Strains.

Aurein 1.2 is an antimicrobial peptide from the skin secretion of an Australian frog. In the previous experimental work, we reported a differential action of aurein 1.2 on two probiotic strains Lactobacillus delbrueckii subsp. bulgaricus (CIDCA 331) and Lactobacillus delbrueckii subsp. lactis (CIDCA 133). The differences found were attributed to the bilayer compositions. Cell cultures and CIDCA 331-derived liposomes showed higher susceptibility than the ones derived from the CIDCA 133 strain, leading to content leakage and structural disruption. Here, we used molecular dynamics simulations to explore these systems at the atomistic level. We hypothesize that if the antimicrobial peptides organized themselves to form a pore, it will be more stable in membranes that emulate the CIDCA 331 strain than in those of the CIDCA 133 strain. To test this hypothesis, we simulated preassembled aurein 1.2 pores embedded into bilayer models that emulate the two probiotic strains. It was found that the general behavior of the systems depends on the composition of the membrane rather than the preassemble system characteristics. Overall, it was observed that aurein 1.2 pores are more stable in the CIDCA 331 model membranes. This fact coincides with the high susceptibility of this strain against antimicrobial peptide. In contrast, in the case of the CIDCA 133 model membranes, peptides migrate to the water-lipid interphase, the pore shrinks, and the transport of water through the pore is reduced. The tendency of glycolipids to make hydrogen bonds with peptides destabilizes the pore structures. This feature is observed to a lesser extent in CIDCA 331 due to the presence of anionic lipids. Glycolipid transverse diffusion (flip-flop) between monolayers occurs in the pore surface region in all the cases considered. These findings expand our understanding of the antimicrobial peptide resistance properties of probiotic strains.

The Effects of Calcium on Lipid-Protein Interactions and Ion Flux in the Cx26 Connexon Embedded into a POPC Bilayer.

Gap junctions provide a communication pathway between adjacent cells. They are formed by paired connexons that reside in the plasma membrane of their respective cell and their activity can be modulated by the bilayer composition. In this work, we study the dynamic behavior of a Cx26 connexon embedded in a POPC lipid bilayer, studying: the membrane protein interactions and the ion flux though the connexon pore. We analyzed extensive atomistic molecular dynamics simulations for different conditions, with and without calcium ions. We found that lipid-protein interactions were mainly mediated by hydrogen bonds. Specific amino acids were identified forming hydrogen bonds with the POPC lipids (ARG98, ARG127, ARG165, ARG216, LYS22, LYS221, LYS223, LYS224, SER19, SER131, SER162, SER219, SER222, THR18 and TYR97, TYR155, TYR212, and TYR217). In the presence of calcium ions, we found subtle differences on the HB lifetimes. Finally, these MD simulations are able to identify and explain differential chlorine flux through the pore depending on the presence or absence of the calcium ions and its distribution within the pore.

Magnesium interactions with a CX26 connexon in lipid bilayers.

Following our previous work, where we described the interaction of calcium with the Cx26 hemichannel, we further explore the same system by atomistic molecular dynamics simulations considering a different di-cation, magnesium. Specifically, the interaction of magnesium di-cation with the previously reported calcium binding sites (ASP2, ASP117, ASP159, GLU114, GLU119, GLU120, and VAL226) was investigated to identify similarities and differences between them. In order to do so, four extensive simulations were carried out. Two of them considered a Cx26 hemichannel embedded on a POPC bilayer with one of the di-cations and a sodium-chlorine solution. For the remaining two, no di-cations were included and a sodium-chlorine or potassium-chlorine solution was considered. Potassium has a similar atomic mass to calcium, and sodium to magnesium, but they both differ in charge (1e and 2e respectively). Magnesium and calcium, even having the same charge, showed different affinity for the explored protein. From the calcium binding sites referred above, we found that the magnesium di-cations only binds strongly to the GLU114 site of one connexin. For the sodium and potassium simulations, no specific interactions with the protein were found. Altogether, these results suggest that mass and steric effects play an important role in determining cation binding to Cx26 hemichannels.

Combining nuclear magnetic resonance with molecular dynamics simulations to address sumatriptan interaction with model membranes.

The goal of this work is to obtain a complete map on the interactions between sumatriptan, an amphiphilic ionizable anti-migraine drug, with lipid bilayers. To this end, we combined two physico-chemical techniques - nuclear magnetic resonance and molecular dynamics simulations - to obtain a detailed picture at different pH values. Both approaches were used considering the strength and constraints of each one. NMR experiments were performed at pH 7.4 where at least 95% of the drug molecules are in their protonated state. From NMR, sumatriptan shows partition on the interfacial region of model membranes (near the head groups and intercalating between adjacent lipids), inducing changes in chemical environment and affecting lipid dynamics of liposomes, in a dose dependent manner. Due to the experimental instability of lipid bilayers at high pH, we took advantage of the molecular dynamics power to emulate different pH values, to simulate sumatriptan in bilayers including at fully uncharged state. Simulations show that the neutral species have preferential orientation within the bilayer interface while the distribution of protonated drugs is independent on the initial conditions. In summary, several properties depicted the interfacial partition of the anti-migraine drug at the water-lipid interface at different conditions. Both techniques were found complementary to shed light on the structural and dynamics of sumatriptan-lipid bilayer interactions.

Rational design of polymer-lipid nanoparticles for docetaxel delivery.

In this work, a stable nanocarrier for the anti-cancer drug docetaxel was rational designed. The nanocarrier was developed based on the solid lipid nanoparticle preparation process aiming to minimize the total amount of excipients used in the final formulations. A particular interest was put on the effects of the polymers in the final composition. In this direction, two poloxoamers -Pluronic F127 and F68- were selected. Some poloxamers are well known to be inhibitors of the P-glycoprotein efflux pump. Additionally, their poly-ethylene-oxide blocks can help them to escape the immune system, making the poloxamers appealing to be present in a nanoparticle designed for the treatment of cancer. Within this context, a factorial experiment design was used to achieve the most suitable formulations, and also to identify the effects of each component on the final (optimized) systems. Two final formulations were chosen with sizes < 250 nm and PDI < 0.2. Then, using dynamic light scattering and nanotracking techniques, the stability of the formulations was assessed during six months. Structural studies were carried on trough different techniques: DSC, x-ray diffraction, FTIR-AR and Molecular Dynamics. The encapsulation efficiency of the anticancer drug docetaxel (> 90%) and its release dynamics from formulations were measured, showing that the polymer-lipid nanoparticle is suitable as a drug delivery system for the treatment of cancer.

Calcium interactions with Cx26 hemmichannel: Spatial association between MD simulations biding sites and variant pathogenicity.

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.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study.

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.

V-Shaped Molecular Configuration of Wax Esters of Jojoba Oil in a Langmuir Film Model.

The aim of the present work was to understand the interfacial properties of a complex mixture of wax esters (WEs) obtained from Jojoba oil (JO). Previously, on the basis of molecular area measurements, a hairpin structure was proposed as the hypothetical configuration of WEs, allowing their organization as compressible monolayers at the air-water interface. In the present work, we contributed with further experimental evidence by combining surface pressure (π), surface potential (Δ V), and PM-IRRAS measurements of JO monolayers and molecular dynamic simulations (MD) on a modified JO model. WEs were self-assembled in Langmuir films. Compression isotherms exhibited πlift-off at 100 Å2/molecule mean molecular area ( Alift-off) and a collapse point at πc ≈ 2.2 mN/m and Ac ≈ 77 Å2/molecule. The Δ V profile reflected two dipolar reorganizations, with one of them at A > Alift-off due to the release of loosely bound water molecules and another one at Ac < A < Alift-off possibly due to reorientations of a more tightly bound water population. This was consistent with the maximal SP value that was calculated according to a model that considered two populations of oriented water and was very close to the experimental value. The orientation of the ester group that was assumed in that calculation was coherent with the PM-IRRAS behavior of the carbonyl group with the C═O oriented toward the water and the C-O oriented parallel to the surface and was in accordance with their orientational angles (∼45 and ∼90°, respectively) determined by MD simulations. Taken together, the present results confirm a V shape rather than a hairpin configuration of WEs at the air-water interface.

A sumatriptan coarse-grained model to explore different environments: interplay with experimental techniques.

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.

Differential Interaction of Antimicrobial Peptides with Lipid Structures Studied by Coarse-Grained Molecular Dynamics Simulations.

In this work; we investigated the differential interaction of amphiphilic antimicrobial peptides with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid structures by means of extensive molecular dynamics simulations. By using a coarse-grained (CG) model within the MARTINI force field; we simulated the peptide-lipid system from three different initial configurations: (a) peptides in water in the presence of a pre-equilibrated lipid bilayer; (b) peptides inside the hydrophobic core of the membrane; and (c) random configurations that allow self-assembled molecular structures. This last approach allowed us to sample the structural space of the systems and consider cooperative effects. The peptides used in our simulations are aurein 1.2 and maculatin 1.1; two well-known antimicrobial peptides from the Australian tree frogs; and molecules that present different membrane-perturbing behaviors. Our results showed differential behaviors for each type of peptide seen in a different organization that could guide a molecular interpretation of the experimental data. While both peptides are capable of forming membrane aggregates; the aurein 1.2 ones have a pore-like structure and exhibit a higher level of organization than those conformed by maculatin 1.1. Furthermore; maculatin 1.1 has a strong tendency to form clusters and induce curvature at low peptide-lipid ratios. The exploration of the possible lipid-peptide structures; as the one carried out here; could be a good tool for recognizing specific configurations that should be further studied with more sophisticated methodologies.

Diblock copolymer bilayers as model for polymersomes: A coarse grain approach.

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 NPNγT ensemble runs reproduce key structural and mechanical experimental properties for different copolymer length chains with a similar hydrophilic weight fraction.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations.

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

Triptan partition in model membranes.

In this work, we report a molecular dynamics simulations study of protonated triptans, sumatriptan and naratriptan, in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC). The simulations were carried out at two concentrations for each drug. Our results show partition between the lipid head-water interphase and water phase for both triptans, with increasing access to the water phase with increasing concentrations. The triptans were stabilized at the interphase through different specific interactions with the POPC bilayer such as hydrogen bonds, salt bridges, and cation-π. Besides, sumatriptan and naratriptan protonated molecules have no access to the hydrophobic region of the bilayer at the studied conditions. Similar results were found for both drugs, however protonated naratriptan shows slightly higher affinity for the water phase. This behavior was attributed to the bulky lateral amino group in its structure under the studied conditions (drugs were originally placed at the water phase). This work represents a first insight to the comprehensive understanding of triptan partition in model membranes.

Concentration effects of sumatriptan on the properties of model membranes by molecular dynamics simulations.

In this work, we report a molecular dynamics simulations study of protonated sumatriptan (pSMT) in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline at the fluid lamellar phase. The simulations were carried out at three different drug/lipid stoichiometries, 1:75, 1:10 and 1:3, under NPT conditions. Our results show partition of pSMT between the lipid head-water interphase and water phase. The main interactions that stabilized the systems were hydrogen bonds, salt bridges and cation-π. Besides, pSMT molecules have no access to the hydrophobic region of the bilayer at the studied concentrations. From an atomistic point of view, this work could contribute to the discussion of drug-membrane interactions regarding the limitation of sumatriptan to cross the blood-brain barrier.

Non-inclusion complexes between riboflavin and cyclodextrins.

OBJECTIVES: To investigate the molecular interaction between ß-cyclodextrin (ßCD) or hydroxypropyl-ß-cyclodextrin (HPßCD) and riboflavin (RF), and to test the anticancer potential of these formulations. METHODS: The physicochemical characterization of the association between RF and CDs was performed by UV-vis absorption, fluorescence, differential scanning calorimetry and NMR techniques. Molecular dynamics simulation was used to shed light on the mechanism of interaction of RF and CDs. Additionally, in-vitro cell culture tests were performed to evaluate the cytotoxicity of the RF-CD complexes against prostate cancer cells. KEY FINDINGS: Neither ßCD nor HPßCD led to substantial changes in the physicochemical properties of RF (with the exception of solubility). Additionally, rotating frame Overhauser effect spectroscopy experiments detected no spatial correlations between hydrogens from the internal cavity of CDs and RF, while molecular dynamics simulations revealed 'out-of-ring' RF-CD interactions. Notwithstanding, both RF-ßCD and RF-HPßCD complexes were cytotoxic to PC3 prostate cancer cells. CONCLUSIONS: The interaction between RF and either ßCD or HPßCD, at low concentrations, seems to be made through hydrogen bonding between the flavonoid and the external rim of both CDs. Regardless of the mechanism of complexation, our findings indicate that RF-CD complexes significantly increase RF solubility and potentiate its antitumour effect.