Understanding the Cellular Uptake of pH-Responsive Zwitterionic Gold Nanoparticles: A Computer Simulation Study

Unveiling Interactions of Tumor-Targeting Nanoparticles with Lipid Bilayers Using a Titratable Martini Model

pH-responsive nanoparticles are ideal vehicles for drug delivery and are widely used in cell imaging in targeted therapy of cancer, which usually has a weakly acidic microenvironment. In this work, we constructed a titratable molecular model for nanoparticles grafted with ligands of pH-sensitive carboxylic acids and investigated the interactions between the nanoparticles and the lipid bilayer in varying pH environments. We mainly examined the effect of the grafting density of the pH-sensitive ligands of the nanoparticles on the interactions of the nanoparticles with the lipid bilayer. The results show that the nanoparticles can penetrate the lipid bilayer only when the pH value is lower than a critical value, which can be readily modulated to the specific pH value of the tumor microenvironment by changing the ligand grafting density. This work provides some insights into modulating the interactions between the pH-sensitive nanoparticles and cellular membranes to realize targeted drug delivery to tumors based on their specific pH environment.

Ligand Lipophilicity Determines Molecular Mechanisms of Nanoparticle Adsorption to Lipid Bilayers

The interactions of ligand-functionalized nanoparticles with the cell membrane affect cellular uptake, cytotoxicity, and related behaviors, but relating these interactions to ligand properties remains challenging. In this work, we perform coarse-grained molecular dynamics simulations to study how the adsorption of ligand-functionalized cationic gold nanoparticles (NPs) to a single-component lipid bilayer (as a model cell membrane) is influenced by ligand end group lipophilicity. A set of 2 nm diameter NPs, each coated with a monolayer of organic ligands that differ only in their end groups, was simulated to mimic NPs recently studied experimentally. Metadynamics calculations were performed to determine key features of the free energy landscape for adsorption as a function of the distance of the NP from the bilayer and the number of NP-lipid contacts. These simulations revealed that NP adsorption is thermodynamically favorable for all NPs due to the extraction of lipids from the bilayer and into the NP monolayer. To resolve ligand-dependent differences in adsorption behavior, string method calculations were performed to compute minimum free energy pathways for adsorption. These calculations revealed a surprising nonmonotonic dependence of the free energy barrier for adsorption on ligand end group lipophilicity. Large free energy barriers are predicted for the least lipophilic end groups because favorable NP-lipid contacts are initiated only through the unfavorable protrusion of lipid tail groups out of the bilayer. The smallest free energy barriers are predicted for end groups of intermediate lipophilicity which promote NP-lipid contacts by intercalating within the bilayer. Unexpectedly, large free energy barriers are also predicted for the most lipophilic end groups which remain sequestered within the ligand monolayer rather than intercalating within the bilayer. These trends are broadly in agreement with past experimental measurements and reveal how subtle variations in ligand lipophilicity dictate adsorption mechanisms and associated kinetics by influencing the interplay of lipid-ligand interactions.

Development of nano-delivery systems for loaded bioactive compounds: using molecular dynamics simulations

Over the past decade, a remarkable surge in the development of functional nano-delivery systems loaded with bioactive compounds for healthcare has been witnessed. Notably, the demanding requirements of high solubility, prolonged circulation, high tissue penetration capability, and strong targeting ability of nanocarriers have posed interdisciplinary research challenges to the community. While extensive experimental studies have been conducted to understand the construction of nano-delivery systems and their metabolic behavior in vivo, less is known about these molecular mechanisms and kinetic pathways during their metabolic process in vivo, and lacking effective means for high-throughput screening. Molecular dynamics (MD) simulation techniques provide a reliable tool for investigating the design of nano-delivery carriers encapsulating these functional ingredients, elucidating the synthesis, translocation, and delivery of nanocarriers. This review introduces the basic MD principles, discusses how to apply MD simulation to design nanocarriers, evaluates the ability of nanocarriers to adhere to or cross gastrointestinal mucosa, and regulates plasma proteins in vivo. Moreover, we presented the critical role of MD simulation in developing delivery systems for precise nutrition and prospects for the future. This review aims to provide insights into the implications of MD simulation techniques for designing and optimizing nano-delivery systems in the healthcare food industry.

How Does the Study MD of pH-Dependent Exposure of Nanoparticles Affect Cellular Uptake of Anticancer Drugs?

The lack of knowledge about the uptake of NPs by biological cells poses a significant problem for drug delivery. For this reason, designing an appropriate model is the main challenge for modelers. To address this problem, molecular modeling studies that can describe the mechanism of cellular uptake of drug-loaded nanoparticles have been conducted in recent decades. In this context, we developed three different models for the amphipathic nature of drug-loaded nanoparticles (MTX-SS-γ-PGA), whose cellular uptake mechanism was predicted by molecular dynamics studies. Many factors affect nanoparticle uptake, including nanoparticle physicochemical properties, protein-particle interactions, and subsequent agglomeration, diffusion, and sedimentation. Therefore, the scientific community needs to understand how these factors can be controlled and the NP uptake of nanoparticles. Based on these considerations, in this study, we investigated for the first time the effects of the selected physicochemical properties of the anticancer drug methotrexate (MTX) grafted with hydrophilic-γ-polyglutamic acid (MTX-SS-γ-PGA) on its cellular uptake at different pH values. To answer this question, we developed three theoretical models describing drug-loaded nanoparticles (MTX-SS-γ-PGA) at three different pH values, such as (1) pH 7.0 (the so-called neutral pH model), (2) pH 6.4 (the so-called tumor pH model), and (3) pH 2.0 (the so-called stomach pH model). Exceptionally, the electron density profile shows that the tumor model interacts more strongly with the head groups of the lipid bilayer than the other models due to charge fluctuations. Hydrogen bonding and RDF analyses provide information about the solution of the NPs with water and their interaction with the lipid bilayer. Finally, dipole moment and HOMO-LUMO analysis showed the free energy of the solution in the water phase and chemical reactivity, which are particularly useful for determining the cellular uptake of the NPs. The proposed study provides fundamental insights into molecular dynamics (MD) that will allow researchers to determine the influence of pH, structure, charge, and energetics of NPs on the cellular uptake of anticancer drugs. We believe that our current study will be useful in developing a new model for drug delivery to cancer cells with a much more efficient and less time-consuming model.

Intramembrane Nanoaggregates of Antimicrobial Peptides Play a Vital Role in Bacterial Killing

Recent developments in antimicrobial peptides (AMPs) have focused on the rational design of short sequences with less than 20 amino acids due to their relatively low synthesis costs and ease of correlation of the structure-function relationship. However, gaps remain in the understanding of how short cationic AMPs interact with the bacterial outer and inner membranes to affect their antimicrobial efficacy and dynamic killing. The membrane-lytic actions of two designed AMPs, G(IIKK)₃ I-NH₂ (G₃ ) and G(IIKK)₄ I-NH₂ (G₄ ), and previously-studied controls GLLDLLKLLLKAAG-NH₂ (LDKA, biomimetic) and GIGAVLKVLTTGLPALISWIKRKR-NH₂ (Melittin, natural) are examined. The mechanistic processes of membrane damage and the disruption strength of the four AMPs are characterized by molecular dynamics simulations and experimental measurements including neutron reflection and scattering. The results from the combined studies are characterized with distinctly different intramembrane nanoaggregates formed upon AMP-specific binding, reflecting clear influences of AMP sequence, charge and the chemistry of the inner and outer membranes. G₃ and G₄ display different nanoaggregation with the outer and inner membranes, and the smaller sizes and further extent of insertion of the intramembrane nanoaggregates into bacterial membranes correlate well with their greater antimicrobial efficacy and faster dynamic killing. This work demonstrates the crucial roles of intramembrane nanoaggregates in optimizing antimicrobial efficacy and dynamic killing.

Modeling Polyzwitterion-Based Drug Delivery Platforms: A Perspective of the Current State-of-the-Art and Beyond

Drug delivery platforms are anticipated to have biocompatible and bioinert surfaces. PEGylation of drug carriers is the most approved method since it improves water solubility and colloid stability and decreases the drug vehicles’ interactions with blood components. Although this approach extends their biocompatibility, biorecognition mechanisms prevent them from biodistribution and thus efficient drug transfer. Recent studies have shown (poly)zwitterions to be alternatives for PEG with superior biocompatibility. (Poly)zwitterions are super hydrophilic, mainly stimuli-responsive, easy to functionalize and they display an extremely low protein adsorption and long biodistribution time. These unique characteristics make them already promising candidates as drug delivery carriers. Furthermore, since they have highly dense charged groups with opposite signs, (poly)zwitterions are intensely hydrated under physiological conditions. This exceptional hydration potential makes them ideal for the design of therapeutic vehicles with antifouling capability, i.e., preventing undesired sorption of biologics from the human body in the drug delivery vehicle. Therefore, (poly)zwitterionic materials have been broadly applied in stimuli-responsive “intelligent” drug delivery systems as well as tumor-targeting carriers because of their excellent biocompatibility, low cytotoxicity, insignificant immunogenicity, high stability, and long circulation time. To tailor (poly)zwitterionic drug vehicles, an interpretation of the structural and stimuli-responsive behavior of this type of polymer is essential. To this end, a direct study of molecular-level interactions, orientations, configurations, and physicochemical properties of (poly)zwitterions is required, which can be achieved via molecular modeling, which has become an influential tool for discovering new materials and understanding diverse material phenomena. As the essential bridge between science and engineering, molecular simulations enable the fundamental understanding of the encapsulation and release behavior of intelligent drug-loaded (poly)zwitterion nanoparticles and can help us to systematically design their next generations. When combined with experiments, modeling can make quantitative predictions. This perspective article aims to illustrate key recent developments in (poly)zwitterion-based drug delivery systems. We summarize how to use predictive multiscale molecular modeling techniques to successfully boost the development of intelligent multifunctional (poly)zwitterions-based systems.

Atomistic Molecular Dynamics Simulation Study on the Interaction between Atomically Precise Thiolate-Protected Gold Nanoclusters and Phospholipid Membranes

The interaction of atomically precise monolayer thiolate (SR) protected gold nanoclusters (Au NCs) with the phospholipid membranes has been studied by the all-atom molecular dynamics (AAMD) simulations. The effect of cluster size, type, and the surface charge density of protection ligand was studied. The simulation results show gold nanoclusters with different size and surface modifications have much different transmembrane behaviors. The Au₂₅(SR)₁₈ cluster was found to possess the best affinity to the phospholipid membranes among six atomically accurate clusters Au₂₅(SR)₁₈, Au₃₆(SR)₂₄, Au₄₄(SR)₂₈, Au₆₈(SR)₃₂, Au₁₄₄(SR)₆₀, and Au₃₁₄(SR)₉₆. Using the Au₂₅ NC as a model, this work also found that the aggregation mode of the surface ligands and the surface charge density are the important factors affecting the interaction between the gold nanoclusters and the phospholipid membranes. Moreover, the balance of hydrophilic and hydrophobic ligands on the surface of Au NCs is beneficial to the high permeability.

The interplay between surface-functionalized gold nanoparticles and negatively charged lipid vesicles

The comprehensive understanding of the interactions between gold nanoparticles (AuNPs) and phospholipid vesicles has important implications in various biomedical applications; however, this is not yet well understood. Here, coarse-grained molecular dynamics (CGMD) simulations were performed to study the interactions between functionalized AuNPs and negatively charged lipid vesicles, and the effects of the surface chemistry and surface charge density (SCD) of AuNPs were analyzed. It is revealed that AuNPs with different surface ligands adhere to the membrane surface (anionic AuNPs) or get into the vesicle bilayer (hydrophobic and cationic AuNPs). Due to the loose arrangement of lipid molecules, AuNPs penetrate curved vesicle membranes more easily than planar lipid bilayers. Cationic AuNPs present three different interaction modes with the vesicle, namely insertion, partial penetration and complete penetration, which are decided by the SCD difference. Both hydrophobic interaction and electrostatic interaction play crucial roles in the interplay between cationic AuNPs and lipid vesicles. For the cationic AuNP with a low SCD, it gets into the lipid bilayer without membrane damage through the hydrophobic interaction, and it is finally stabilized in the hydrophobic interior of the vesicle membrane in a thermodynamically stable "snorkeling" configuration. For the cationic AuNP with a high SCD, it crosses the vesicle membrane and gets into the vesicle core through a membrane pore induced by strong electrostatic interaction. In this process, the membrane structure is destroyed. These findings provide a molecular-level understanding of the interplay between AuNPs and lipid vesicles, which may further expand the application of functional AuNPs in modern biomedicine.

Computer Simulations on a pH-Responsive Anticancer Drug Delivery System Using Zwitterion-Grafted Polyamidoamine Dendrimer Unimolecular Micelles

Unimolecular micelles have attracted wide attention in the field of drug delivery because of their thermodynamic stability and uniform size distribution. However, their drug loading/release mechanisms at the molecular level have been poorly understood. In this work, the stability and drug loading/release behaviors of unimolecular micelles formed using generation-5 polyamidoamine-graft-poly(carboxybetaine methacrylate) (PAMAM(G5)-PCBMA) were studied by dissipative particle dynamics simulations. In addition, the unimolecular micelles formed using generation-5 polyamidoamine-graft-poly(ethyleneglycol methacrylate) (PAMAM(G5)-PEGMA) were used as a comparison. The simulation results showed that PAMAM(G5)-PCBMA can spontaneously form core-shell unimolecular micelles. The PAMAM(G5) dendrimer constitutes a hydrophobic core to load the doxorubicin (DOX), while the zwitterionic PCBMA serves as a protective shell to improve the stability of the unimolecular micelle. The DOX can be encapsulated into the cavity of PAMAM(G5) at the physiological pH 7.4. The drug loading efficiency and drug loading content showed some regularities with the increase in the drug concentration. At the acidic pH 5.0, the loaded DOX can be released gradually from the hydrophobic core. The comparison of DOX-loaded morphologies between the PAMAM(G5)-PCBMA system and PAMAM(G5)-PEGMA system showed that the former has better monodisperse stability. This work could offer theoretical guidance for the design and development of promising unimolecular micelles for drug delivery.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Molecular mechanism of HIV-1 TAT peptide and its conjugated gold nanoparticles translocating across lipid membranes

The trans-acting activator of transcription (TAT) peptide, which is derived from human immunodeficiency virus-1 (HIV-1), has been widely used as an effective nanocarrier to transport extracellular substances into cells. However, the underlying translocation mechanism of TAT peptide across cell membranes still remains controversial. Besides, the molecular process of TAT peptide facilitating the transport of extracellular substances into cells is largely unknown. In this study, we explore the interactions of TAT peptides and their conjugated gold nanoparticles with lipid membranes by coarse-grained molecular dynamics simulations. It is found that the TAT peptides can hardly penetrate through the membrane at low peptide concentrations; after the concentration increases to a threshold value, they can cross the membrane through an induced nanopore due to the transmembrane electrostatic potential difference. The translocation of TAT peptides is mainly caused by the overall structural changes of membranes. Furthermore, we demonstrate that the translocation of gold nanoparticles (AuNPs) across the membrane is significantly affected by the number of grafted TAT peptides on the particle surface. The transmembrane efficiency of AuNPs may even be reduced when a small number of peptides modify them; whereas, when the number of grafted peptides increases to a certain value, the TAT-AuNP complex can translocate across the membrane in a pore-mediated way. Based on our findings, an effective strategy has been proposed to enhance the delivery efficiency of AuNPs. The present study can improve our understanding of the interactions between TAT peptides and cell membranes; it may also give some insightful suggestions on the design and development of nanocarriers with high efficiency for the delivery of nanoparticles and drugs.

pH-Responsive Zwitterionic Copolymer DHA-PBLG-PCB for Targeted Drug Delivery: A Computer Simulation Study

In this work, the self-assembled behaviors of zwitterionic copolymer docosahexaenoic acid- b-poly(γ-benzyl-l-glutamate)- b-poly(carboxybetaine methacrylate) (DHA-PBLG-PCB) and the loading and release mechanism of the anticancer drug doxorubicin (DOX) was investigated via computer simulations. The effects of polymer concentration, drug content, and pH on polymeric micelles were explored by dissipative particle dynamics (DPD) simulations. Simulation results show that DHA-PBLG₁₅-PCB₁₀ can self-assemble into core-shell micelles; in addition, the drug-loaded micelles have a pH-responsive feature. DOX can be encapsulated into the core-shell micelle under normal physiological pH conditions, whereas it can be released under acidic pH conditions. The self-assembled behaviors of copolymer DHA-PBLG-PEG were also studied to have a comparison with those of DHA-PBLG-PCB. The DHA-PBLG₁₅-PCB₁₀ system has a stable structure and it has a great potential to serve as drug delivery vehicles for targeted drug delivery.

Energy landscape for the insertion of amphiphilic nanoparticles into lipid membranes: A computational study

Amphiphilic, monolayer-protected gold nanoparticles (NPs) have been shown to enter cells via a non-endocytic, non-disruptive pathway that could be valuable for biomedical applications. The same NPs were also found to insert into a series of model cell membranes as a precursor to cellular uptake, but the insertion mechanism remains unclear. Previous simulations have demonstrated that an amphiphilic NP can insert into a single leaflet of a planar lipid bilayer, but in this configuration all charged end groups are localized to one side of the bilayer and it is unknown if further insertion is thermodynamically favorable. Here, we use atomistic molecular dynamics simulations to show that an amphiphilic NP can reach the bilayer midplane non-disruptively if charged ligands iteratively "flip" across the bilayer. Ligand flipping is a favorable process that relaxes bilayer curvature, decreases the nonpolar solvent-accessible surface area of the NP monolayer, and increases attractive ligand-lipid electrostatic interactions. Analysis of end group hydration further indicates that iterative ligand flipping can occur on experimentally relevant timescales. Supported by these results, we present a complete energy landscape for the non-disruptive insertion of amphiphilic NPs into lipid bilayers. These findings will help guide the design of NPs to enhance bilayer insertion and non-endocytic cellular uptake, and also provide physical insight into a possible pathway for the translocation of charged biomacromolecules.

Surface Charge-Dependent Cellular Uptake of Polystyrene Nanoparticles

The evaluation of the role of physicochemical properties in the toxicity of nanoparticles is important for the understanding of toxicity mechanisms and for controlling the behavior of nanoparticles. The surface charge of nanoparticles is suggested as one of the key parameters which decide their biological impact. In this study, we synthesized fluorophore-conjugated polystyrene nanoparticles (F-PLNPs), with seven different types of surface functional groups that were all based on an identical core, to evaluate the role of surface charge in the cellular uptake of nanoparticles. Phagocytic differentiated THP-1 cells or non-phagocytic A549 cells were incubated with F-PLNP for 4 h, and their cellular uptake was quantified by fluorescence intensity and confocal microscopy. The amount of internalized F-PLNPs showed a good positive correlation with the zeta potential of F-PLNPs in both cell lines (Pearson’s r = 0.7021 and 0.7852 for zeta potential vs. cellular uptake in THP-1 cells and nonphagocytic A549 cells, respectively). This result implies that surface charge is the major parameter determining cellular uptake efficiency, although other factors such as aggregation/agglomeration, protein corona formation, and compositional elements can also influence the cellular uptake partly or indirectly.

Progress in ligand design for monolayer-protected nanoparticles for nanobio interfaces

Ligand-functionalized inorganic nanoparticles, also known as monolayer-protected nanoparticles, offer great potential as vehicles for in vivo delivery of drugs, genes, and other therapeutics. These nanoparticles offer highly customizable chemistries independent of the size, shape, and functionality imparted by the inorganic core. Their success as drug delivery agents depends on their interaction with three major classes of biomolecules: nucleic acids, proteins, and membranes. Here, the authors discuss recent advances and open questions in the field of nanoparticle ligand design for nanomedicine, with a focus on atomic-scale interactions with biomolecules. While the importance of charge and hydrophobicity of ligands for biocompatibility and cell internalization has been demonstrated, ligand length, flexibility, branchedness, and other properties also influence the properties of nanoparticles. However, a comprehensive understanding of ligand design principles lies in the cost associated with synthesizing and characterizing diverse ligand chemistries and the ability to carefully assess the structural integrity of biomolecules upon interactions with nanoparticles.