Steered molecular dynamics simulations reveal a self-protecting configuration of nanoparticles during membrane penetration

Computational assessment of lipid facilitated membrane permeation of vancomycin using force-probe molecular dynamic simulation

In this study the efficacy of different edible lipids for drug permeation enhancement of vancomycin through biological membrane was investigated using molecular dynamic simulation. In this regard, at first the ability of the lipids for complex formation with the drug was evaluated for number of most common edible lipids including tripalmitin (TPA), trimyristin (TMY), labrafil (LAB), glycerol monostearate (GMS), glycerol monooleate (GMO), Distearoylphosphorylethanolamine (DSPE), dipalmitoylphosphatidylethanolamine (DPPE), Dipalmitoylphosphatidylcholine (DPPC), cholesterol (CL), stearic acid (SA), palmitic acid (PA) and oleic acid (OA). Then the complexes were pulled thorough a bilayer membrane while the changes in force were probed. The results showed that besides the SA, PA and OA the other examined lipids were able to perform a perfect molecular complex with the drug. Also the results of pulling simulation revealed that the least of force was needed for drug transmittance through the membrane when it was covered by LAB, TMY and DSPE. These results indicated that these lipids can be the excellent materials of choice as permeation enhancer for preparing a proper oral formulation of vancomycin.Communicated by Ramaswamy H. Sarma.

Computer-aided nanodrug discovery: recent progress and future prospects

Nanodrugs, which utilise nanomaterials in disease prevention and therapy, have attracted considerable interest since their initial conceptualisation in the 1990s. Substantial efforts have been made to develop nanodrugs for overcoming the limitations of conventional drugs, such as low targeting efficacy, high dosage and toxicity, and potential drug resistance. Despite the significant progress that has been made in nanodrug discovery, the precise design or screening of nanomaterials with desired biomedical functions prior to experimentation remains a significant challenge. This is particularly the case with regard to personalised precision nanodrugs, which require the simultaneous optimisation of the structures, compositions, and surface functionalities of nanodrugs. The development of powerful computer clusters and algorithms has made it possible to overcome this challenge through in silico methods, which provide a comprehensive understanding of the medical functions of nanodrugs in relation to their physicochemical properties. In addition, machine learning techniques have been widely employed in nanodrug research, significantly accelerating the understanding of bio-nano interactions and the development of nanodrugs. This review will present a summary of the computational advances in nanodrug discovery, focusing on the understanding of how the key interfacial interactions, namely, surface adsorption, supramolecular recognition, surface catalysis, and chemical conversion, affect the therapeutic efficacy of nanodrugs. Furthermore, this review will discuss the challenges and opportunities in computer-aided nanodrug discovery, with particular emphasis on the integrated "computation + machine learning + experimentation" strategy that can potentially accelerate the discovery of precision nanodrugs.

Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways

Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.

Investigation of the ionic conditions in SiRNA-mediated delivery through its carriers in the cell membrane: a molecular dynamic simulation

SiRNA is a new generation of drug molecules and a new approach for treating a variety of diseases such as cancer and viral infections. SiRNA delivery to cells and translocation into cytoplasm are the main challenges in the clinical application of siRNA. Lipid carriers are one of the most successful carriers for siRNA delivery. In this study, we investigated the interaction of siRNA with a zwitterionic bilayer and how ion concentration and lipid conjugation can affect it. The divalent cation such as Mg²⁺ ions could promote the siRNA adsorption on the bilayer surface. The cation ions can bind to the head groups of lipids and the grooves of siRNA molecules and form bridges between the siRNA and bilayer surface. Our findings demonstrated the bridges formed by divalent ions could facilitate the attachment of siRNA to the membrane surface. We showed that the divalent cations can regulate the bridging-driven membrane attachment and it seems the result of this modulation can be used for designing biomimetic devices. In the following, we examined the effect of cations on the interaction between siRNA modified by cholesterol and the membrane surface. Our MD simulations showed that in the presence of Mg²⁺, the electrostatic and vdW energy between the membrane and siRNA were higher compared to those in the presence of NA⁺. We showed that the electrostatic interaction between membrane and siRNA cannot be facilitated only by cholesterol conjugated. Indeed, cations are essential to create coulomb repulsion and enable membrane attachment. This study provides important insight into liposome carriers for siRNA delivery and could help us in the development of siRNA-based therapeutics. Due to the coronavirus pandemic outbreak, these results may shed light on the new approach for treating these diseases and their molecular details.

The penetration of human defensin 5 (HD5) through bacterial outer membrane: simulation studies

Human α-defensin 5 (HD5) is one of cationic antimicrobial peptides which plays a crucial role in an innate immune system in human body. HD5 shows the killing activity against a broad spectrum of pathogenic bacteria by making a pore in a bacterial membrane and penetrating into a cytosol. Nonetheless, its pore-forming mechanisms remain unclear. Thus, in this work, the constant-velocity steered molecular dynamics (SMD) simulation was used to simulate the permeation of a dimeric HD5 into a gram-negative lipopolysaccharide (LPS) membrane model. Arginine-rich HD5 is found to strongly interact with a LPS surface. Upon arrival, arginines on HD5 interact with lipid A head groups (a top part of LPS) and then drag these charged moieties down into a hydrophobic core resulting in the formation of water-filled pore. Although all arginines are found to interact with a membrane, Arg13 and Arg32 appear to play a dominant role in the HD5 adsorption on a gram-negative membrane. Furthermore, one chain of a dimeric HD5 is required for HD5 adhesion. The interactions of arginine-lipid A head groups play a major role in adhering a cationic HD5 on a membrane surface and retarding a HD5 passage in the meantime.

Metadynamics-Based Approaches for Modeling the Hypoxia-Inducible Factor 2α Ligand Binding Process

Several methods based on enhanced-sampling molecular dynamics have been proposed for studying ligand binding processes. Here, we developed a protocol that combines the advantages of steered molecular dynamics (SMD) and metadynamics. While SMD is proposed for investigating possible unbinding pathways of the ligand and identifying the preferred one, metadynamics, with the path collective variable (PCV) formalism, is suggested to explore the binding processes along the pathway defined on the basis of SMD, by using only two CVs. We applied our approach to the study of binding of two known ligands to the hypoxia-inducible factor 2α, where the buried binding cavity makes simulation of the process a challenging task. Our approach allowed identification of the preferred entrance pathway for each ligand, highlighted the features of the bound and intermediate states in the free-energy surface, and provided a binding affinity scale in agreement with experimental data. Therefore, it seems to be a suitable tool for elucidating ligand binding processes of similar complex systems.

Membrane lipids destabilize short interfering ribonucleic acid (siRNA)/polyethylenimine nanoparticles

Cell entry of polymeric nanoparticles (NPs) bearing polynucleotides is an important stage for successful gene delivery. In this work, we addressed the influence of cell membrane lipids on the integrity and configurational changes of NPs composed of short interfering ribonucleic acid (siRNA) and polyethylenimine. We focused on NPs derived from two different PEIs, unmodified low molecular weight PEI and linoleic acid (LA)-substituted PEI, and their interactions with two membrane lipids (zwitterionic 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) and anionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)). Our experiments showed that POPS liposomes interacted strongly with both types of NPs, which caused partial dissociation of the NPs. POPC liposomes, however, did not induce any dissociation. Consistent with the experiments, steered molecular dynamics simulations showed a stronger interaction between the NPs and the POPS membrane than between the NPs and the POPC membrane. Lipid substitution on the PEIs enhanced the stability of the NPs during membrane crossing; lipid association between PEIs of the LA-bearing NPs as well as parallel orientation of the siRNAs provided protection against their dissociation (unlike NPs from native PEI). Our observations provide valuable insight into the integrity and structural changes of PEI/siRNA NPs during membrane crossing which will help in the design of more effective carriers for nucleic acid delivery.

pH-Dependent aggregation and pH-independent cell membrane adhesion of monolayer-protected mixed charged gold nanoparticles

Design of pH-responsive monolayer-protected gold nanoparticles (AuNPs) that are mixed charged, with the ability to switch their net surface charge, based on the stimuli of environmental pH is a promising technique in nanomedicine. However, understanding of pH-responsive mixed charged AuNP behavior in terms of their stability and cellular interaction are still limited. In this work, we study the aggregation of pH-responsive AuNPs and their interaction with model lipid bilayers by adopting Martini coarse-grained (CG) molecular dynamics simulations. The surface of these AuNPs is decorated by both positively and negatively charged ligands. The AuNP is positively charged at low pH values due to protonation of negatively charged ligands. Its net charge is lowered at higher pH by increasing the ratio of deprotonated negatively charged ligands. We find that the AuNPs are severely aggregated at moderate pH value, where each AuNP has an overall neutral charge, whereas they are stable and dispersed at both low and high pH values. Further free energy analysis reveals that the energy barrier at a larger separation distance than the location of the hydrophobic driving force potential well, plays a key role that determines the stability of monolayer-protected AuNPs at different pH values. This energy barrier is dramatically decreased at moderate pH value, leading to severe aggregation of AuNPs. By investigating the interaction between AuNPs and model lipid bilayers, we find that all the AuNPs adhere onto the lipid bilayer, independent of the pH value. Moreover, the lipids present originally in the bilayer are extracted by the AuNPs through a process of protrusion and upward climbing. The extraction of lipids can cause dehydration and disruption of the bilayers when multiple AuNPs are adhered. Free energy analysis reveals that the penetration of AuNPs will induce a dramatic free energy increase because of deformation of the ligands with hydrophilic functional end groups. We have systematically studied the stability of pH-responsive AuNPs and their interactions with lipid bilayers by simulation, which might pave the way for the design of pH-responsive monolayer protected AuNPs for biomedical applications.