Next Article in Journal
Guided Bone Regeneration in a Periodontally Compromised Individual with Autogenous Tooth Bone Graft: A Radiomics Analysis
Previous Article in Journal
Bovine Fibroblast-Derived Extracellular Matrix Promotes the Growth and Preserves the Stemness of Bovine Stromal Cells during In Vitro Expansion
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JFB
Volume 14
Issue 4
10.3390/jfb14040219
jfb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Interaction Mechanism of Intramuscular Gene Delivery Materials with Cell Membranes

by
Zhanpeng Cui
1,
Yang Jiao
1,
Linyu Pu
2,
Jianlin Chen
3,
Ming Liu
4,
James Zhenggui Tang
5 and
Gang Wang
1,*
1
National Engineering Research Center for Biomaterials, College of Biomedical Engineering, Sichuan University, Chengdu 610064, China
2
School of Materials and Chemistry, Southwest University of Science and Technology, Mianyang 621010, China
3
School of Laboratory Medicine, Chengdu Medical College, Chengdu 610500, China
4
Department of Medical Oncology/Gastric Cancer Center, West China Hospital, Sichuan University, Chengdu 610041, China
5
Research Institute of Healthcare Science, Faculty of Science & Engineering, University of Wolverhampton, Wolverhampton WV1 1SB, UK
*
Author to whom correspondence should be addressed.
J. Funct. Biomater. 2023, 14(4), 219; https://doi.org/10.3390/jfb14040219
Submission received: 17 February 2023 / Revised: 25 March 2023 / Accepted: 8 April 2023 / Published: 14 April 2023
(This article belongs to the Section Biomaterials for Drug Delivery)
Browse Figures
Versions Notes

Abstract

:
It has been confirmed that skeletal muscle cells have the capability to receive foreign plasmid DNA (pDNA) and express functional proteins. This provides a promisingly applicable strategy for safe, convenient, and economical gene therapy. However, intramuscular pDNA delivery efficiency was not high enough for most therapeutic purposes. Some non-viral biomaterials, especially several amphiphilic triblock copolymers, have been shown to significantly improve intramuscular gene delivery efficiency, but the detailed process and mechanism are still not well understood. In this study, the molecular dynamics simulation method was applied to investigate the structure and energy changes of the material molecules, the cell membrane, and the DNA molecules at the atomic and molecular levels. From the results, the interaction process and mechanism of the material molecules with the cell membrane were revealed, and more importantly, the simulation results almost completely matched the previous experimental results. This study may help us design and optimize better intramuscular gene delivery materials for clinical applications.

1. Introduction

In 1990, Wolff et al. injected naked plasmid (pDNA) into the skeletal muscles of mice and found that skeletal muscle cells could not only take up the pDNA but also express the corresponding protein [1]. This discovery opened a new direction for gene therapy. Although successful trials were acquired in the treatment of lower extremity ischemia [2], further applications, such as the treatment of cancers and diabetes, struggled to achieve satisfactory outcomes since the expression levels of functional proteins were not high enough for therapeutic purpose, which was mainly due to the extremely low efficiency of pDNA delivery into skeletal muscle cells [3,4]. Viral vectors have high gene delivery efficiency, but they are prone to trigger inflammatory responses and immunological reactions, as well as posing a risk of oncogenesis [5,6,7,8]. It was reported that an overdose injection of viral vector triggered an excessively strong inflammatory response, and eventually resulted in the death of the patient [9]. Therefore, non-viral gene carriers were anticipated to play an important role in the transfer of pDNA into skeletal muscle cells. In 1996, Mumper et al. mixed polyvinyl pyrrolidone (PVP) and polyvinyl alcohol (PVA) with pDNA, injected the mixture into the tibialis anterior (TA) muscles of mice and acquired an improved expression level of exogenous genes, which was about 10 times that achieved by pure pDNA injection alone [10]. This finding confirmed that pDNA delivery efficiency to skeletal muscles could be improved with the help of non-viral materials, and these materials were safer than viral vectors.
Non-viral gene delivery materials have been developed thereafter. Several new materials have been reported, such as PEG-Plys block copolymer [11], PEG-conjugated dendritic molecule PAMAM [12], and triblock copolymer PEG-PLGA-PEG [13]. Later, as commonly used excipients, Pluronic compounds came of interest to researchers. Pluronic compounds have an amphiphilic hydrophilic-hydrophobic-hydrophilic or hydrophobic-hydrophilic-hydrophobic triblock structure. It has been reported that Pluronic L64, P85, F127, P104, and SP1017 were able to improve the efficiency of gene transfer into skeletal muscles [14,15,16,17]. In our previous studies, Pluronic L64 was shown to be one of the most efficient and stable materials for pDNA delivery to skeletal muscles [2], and higher gene delivery efficiency was acquired when Pluronic L64 was combined with a low-voltage electric field [18,19], and was even high enough for therapeutic trials on tumors and diabetes [20,21]. Further, inspired by the similarity of amphiphilic structure between Pluronic L64 and cell membranes, amphiphilic triblock copolymers of LPL [22], rL2PL2, and rL3PL3 were designed and synthesized [23]. The results showed that LPL and rL2PL2 had better performance than Pluronic L64, evident from intramuscular gene delivery/expression efficiency and biocompatibility. Based on these studies, it can be concluded that amphiphilic triblock copolymers are able to disturb the structural stability of cytomembranes, leading to increased permeability, and the electric neutrality of the material molecules can prevent them from being trapped in the extracellular matrix of skeletal muscle tissues [24]. Therefore, electroneutral amphiphilic molecules with both “hydrophilic-hydrophobic-hydrophilic” and “hydrophobic-hydrophilic-hydrophobic” structures are ideal models for molecule design to optimize intramuscular gene delivery materials. However, the materials of LPL, rL2PL2, and rL3PL3 were all designed based on bionics, and the interaction mechanism of the materials with cytomembranes was only revealed at the cellular level [25]. To optimize the molecule design, it is necessary to understand the interaction process and mechanism between the materials and cell membranes at the molecular or even atomic scale.
In the early 1980s, molecular dynamics (MD) simulation was first applied to study biological systems [26]. With the development of computer technology, MD simulation has been applied to many studies on phospholipid membranes, such as the damage of graphene on cell membranes [27], the damage of graphene oxide on red blood cell membranes [28], the mechanism of action of cell-penetrating peptides [29], and the mechanism of cell membrane toxicity of petroleum oil droplets [30]. Therefore, MD simulation provides a novel tool for this kind of research. MD simulation is based on Newtonian mechanics, and a large number of theoretical calculations could be used to describe the interaction between molecules to obtain the mechanism of action at the atomic scale [31,32]. Furthermore, the conformations produced in MD simulation can be displayed in three dimensions by visualization software, which makes the outputs of MD more intuitive and easier to understand [33].
Compared with traditional experimental methods, MD simulation can be used to study the interaction between molecules in a system by establishing an appropriate model. Consequently, it can effectively shorten the experimental period, improve experimental efficiency, and reduce the consumption of experimental resources. From the perspectives of gene delivery materials, the data files obtained by MD simulation can be used to reveal the detailed interaction process among material molecules, cell membranes, and DNA, and to further explore the interaction mechanisms by means of structure and energy analysis.
Although there were many literature reports on the MD simulation of phospholipid membranes [34,35,36,37,38], few reports were found on the MD simulation regarding the interaction of gene delivery materials, especially amphiphilic triblocks, with cell membranes, and many conjectures were still derived from experience and assumptions [2,22,23,29,39,40,41,42]. Therefore, we will use MD simulation, taking the four previously reported materials L64, LPL, rL2PL2, and rL3PL3 as samples, to explore the mechanisms of their interactions with cell membranes by dynamically analyzing the intermolecular interactions and energy changes. It is anticipated that understanding the interaction process and mechanism at molecular and atomic levels will provide new inspiration for designing better materials for intramuscular gene delivery.

2. Computational Details

2.1. The Force Field Selection

The Charmm force field was selected for all simulations in this study. The force field version is Charmm36, which is an all-atom force field and can be used to study the interaction among nucleic acids, proteins, lipids, and other molecules [43,44]. The water model is TIP3P [45], which is suitable for the CHARMM force field.

2.2. The MD Simulation Software

The molecular dynamics simulation software used in this study was Gromacs 2018.4 [46,47,48]. The visualization software was VMD 1.9.3 [49].

2.3. Modeling of Material Molecules

In this study, four kinds of gene delivery polymer molecules, L64 and LPL (https://pubs.acs.org/doi/10.1021/am503808b, accessed on 15 August 2014), and rL2PL2 and rL3PL3 (https://pubs.acs.org/doi/10.1021/acsami.6b02592, accessed on 16 May 2016) were selected as the experimental objects. These four molecules are pure polymers. The structures of the molecules used in this paper are based on the actual experimental data from the previous studies of Liu S et al. [18] and Pu L et al. [22,23]. The experimental parameters of the gene delivery polymers are shown in Table 1.
The molecular structures are shown in Figure 1. The force field parameter files of the materials were obtained from the website CGenFF [50,51,52]. The total charge of each molecule was 0.
NVT stands for the canonical ensemble, a pre-equilibration step in which there are three constant quantities, namely, constant number of particles, constant volume, and constant temperature.
NPT, an isothermal isobaric ensemble, is the constant temperature and constant pressure ensemble, which is also a pre-equilibration step that includes three constant quantities, namely, constant particle number, constant temperature, and constant pressure. This approach can make the particle arrangement in the system approximate the actual particle arrangement, thereby improving the fidelity of the model and turning the system into a stable state.
After modeling the four material molecules, the energy in a vacuum was minimized in a grid with V = 30.0 nm × 30.0 nm × 30.0 nm, and the convergence limit was set to 100.0 kJ/mol·nm. Then, the NPT simulation at 500.0 ps was performed, and the pre-equilibration end conformation was taken as the initial conformation for the subsequent formal simulation.
The DNA molecule used in this study consisted of 14 consecutive A-T base pairs and was built using Avogadro 1.2.0 software [53,54]. The charge of the DNA molecule model in the environment of normal saline was −28 e−. The force field parameter files of the DNA were obtained from the website CGenFF [50,51,52].

2.4. Establishment of Cell Membranes

The cytoplasmic membrane was selected as a cell membrane containing 240 lipid molecules. It is composed of POPC (Figure 1E) and POPS (Figure 1F). The phase transition temperatures of POPC and POPS were set to 271.0 K and 300.0 K [55,56], respectively. The number ratio of the two lipid molecules, POPC:POPS was set to 4:1. The cell membrane was established by the Charmm-GUI server [57,58,59]. The thickness of the phospholipid membrane was about 4.0 nm. The membrane was located in the center of the grid (9.0 nm × 9.0 nm × 16.0 nm). The surface of the phospholipid membrane was perpendicular to the z-axis. The concentration of NaCl outside the membrane was 0.15 M.

2.5. Molecular Dynamics Simulation Ideas

In this study, we established 21 systems. The information on lattice size and the number of water molecules is shown in Supplementary Materials Table S1.
The deuterium order parameter can be written as:
S CD = 3 cos 2 θ i − 1 2
where θi is the angle between the vector connecting the i − 1 and i + 1 beads of the lipid tail and the bilayer normal z, and SCD is the average of the two chains of the same lipids in the entire bilayer and the designated simulation time. This allows a more quantitative assessment of the order of the phospholipid’s molecular orientation in membrane phospholipids. In this system, the order parameter can reflect the orderly arrangement of phospholipid molecules in the phospholipid membrane to quantitatively reflect the disturbance degree of the material to the cell membrane.
After the modeling was completed, a pre-equilibrium simulation was performed. The steps included energy minimization, NVT, NPT1, and NPT2. NVT and NPT1 limited the position of the phospholipid bilayer in the z-axis direction and the position of the material molecule in the three-dimensional direction. When using the modeling command to fill water molecules, the spatial orientation of all water molecules is consistent, that is, the high orientation of water molecules is consistent. This orientation certainly does not conform to the real situation. Therefore, before the simulation starts, we need first to carry out a pre-equilibrium step, that is, during the pre-equilibrium period, limit the orientation of material molecules and phospholipid molecules, and only let water molecules move in an irregular manner to approach the true situation of the irregular arrangement of water molecules in the natural state. The force constant of position limitation was 1000.0 kJ/mol·nm2, which was used to balance the orientation of the water molecules. The energy minimization adopted the steepest descent method. The energy convergence limit was set at 100 kJ/mol·nm, and the maximum number of calculation steps for energy minimization was set at 100,000.
The temperature control algorithm of NVT used the Berendsen method [60]. The pre-equilibration temperature was 310.0 K, and the duration was 1.0 ns. The temperature and pressure control algorithms of NPT1 also used the Berendsen method. The temperature and pressure control were set to 310.0 K and 1.0 bar, respectively, and the duration was 2.0 ns. NPT2 canceled the position restriction of the phospholipid membrane on the basis of NPT1 and continued to equilibrate for 2.0 ns.
For formal MD simulations, the temperature control algorithm used the V-rescale method and the voltage control algorithm used the Parrinello–Rahman method [61,62]. The simulation duration was 100 ns. The simulation step size was 0.001 ps (corresponding to Systems A-E, K-N). In systems S-V, the simulation duration was 300 ns. The process flow of this study is shown in Scheme 1.:
The potential of mean force (PMF) refers to the potential energy obtained by integrating the mean force. When the direction of the average force is one-dimensional, the PMF can also reflect changes in free energy, thereby determining the spontaneity of intermolecular interactions [63,64,65]. (corresponding to Systems F-J, O-R).
The initial modeling legend is shown in Supplementary Materials Figure S1.
In systems A, B, D, and E, in order to study the disturbing effect of the four material molecules on the cell membrane, we inserted four material molecules into the cell membrane at 1.0 nm, created a phospholipid membrane disturbing environment (corresponding to Figure 2, Figure 3A, Figure 4A and Figure 5), each system was repeated three times independently. The configurations were chosen every 0.1 ns during the molecular dynamics simulations.
In system C, we put LPL directly into the center of the phospholipid membrane and observe its disturbing effect (corresponding to LPL in), each system was repeated three times independently. The configurations were chosen every 0.1 ns during the molecular dynamics simulations.
In system F-J, we place the four material molecules above the phospholipid membrane, and the distance between them and the centroid of the phospholipid membrane is 4.0 nm. This initial model is used to calculate the change in free energy of the four material molecules, and the calculation methods were PMF and WHAM. For the PMF simulation, we apply a constant force of 1000.0 kJ/mol·nm2 perpendicular to the phospholipid membrane and set its maximum traction speed as 0.01 nm/ps. The traction distance is 9.0 nm, so we set the sampling window to 0.20 nm, and the sampling window is 50, to ensure that there is enough intersection between the two adjacent windows. The pre-equilibration time of each window was 300.0 ps, and the relaxation time was 8.0 ns. After the relaxation of all windows was completed, the weighted distribution histogram was drawn using the WHAM module of Gromacs [66]. The change curve of PMF free energy was obtained, and each system was repeated twice independently.
In system K-N, we put the material molecules and DNA into a lattice separately, and the distance between the DNA and material molecules is 0.3 nm—close but without contact—and each system was repeated three times independently. The configurations were chosen every 0.1 ns during the molecular dynamics simulations.
In system O-R, we put DNA and different phospholipid membranes into a system, where DNA molecules are placed above the phospholipid membrane. The distance from the centroid of the phospholipid membrane is 5.0 nm, and the double helix direction of DNA is parallel to the surface of the phospholipid membrane. This initial model is used to calculate the change in DNA transmembrane free energy, and the calculation methods were PMF and WHAM. For the PMF simulation, we apply a constant force of 1000.0 kJ/mol·nm2 perpendicular to the phospholipid membrane and set its maximum traction speed as 0.01 nm/ps. The traction distance is 10.0 nm, so we set the sampling window to 0.18 nm, and the sampling window is 60, to ensure that there is enough intersection between the two adjacent windows. The pre-equilibration time of each window was 300.0 ps, and the relaxation time was 8.0 ns. After the relaxation of all windows was completed, the weighted distribution histogram was drawn using the WHAM module of Gromacs [66], and each system was repeated twice independently.
In system S-U, we put L64, LPL, and rL2PL2 at 1.0 nm above the surface of the phospholipid membrane to observe how the material molecules enter the phospholipid membrane.

3. Results and Discussion

3.1. Interaction Process between Material Molecules and Phospholipid Membranes

In order to observe how the material molecules interacted with cell membranes, the material molecules were inserted about 1.0 nm into the membrane, and molecular dynamics simulation was performed for 100 ns. It can be found that when L64 was not repelled by the cell membrane, it flowed with the phospholipid membrane and gradually entered the centroid of the phospholipid membrane. The molecules became linearly arranged with time, and the hydrophobic segment flowed stably in the cell membrane while the hydrophilic end swung in the solvent (Figure 2).
When L64 molecules were located above the phospholipid membrane, they floated randomly outside the membrane and were not did not easily come into contact with the phospholipid membrane. In order to increase the probability of L64 entering the cell membrane, we set up a bilayer membranes model and put L64 molecules between the two membranes. In this system, L64 molecules were found to automatically enter the phospholipid membrane (Supplementary Materials Figure S2). At 139.2 ns, the hydrophobic part of L64 began to come into contact with the hydrophobic region of the membrane. Then, the L64 molecule gradually moved toward the center of the membrane, and the final conformation was consistent with that when L64 was inserted 1.0 nm into the phospholipid membrane. The overall extension of L64 became linear. Its hydrophobic segment flowed stably in the hydrophobic area of the phospholipid membrane, while the hydrophilic end swung randomly in the water. The results revealed that the entry of L64 molecules into the cell membrane mainly relied on the hydrophobic interaction.
As the LPL molecule has the “hydrophobic-hydrophilic-hydrophobic” structure, whereas molecules of L64, rL2PL2, and rL3PL3 are “hydrophilic-hydrophobic-hydrophilic” in arrangement, the LPL molecule was initially inserted 1.0 nm into the phospholipid membrane. About 14.8 ns later, one of the hydrophobic segments re-entered the water and floated there, while another hydrophobic segment end was still kept in the hydrophobic region of the membrane. This pattern was maintained until the end of the simulation (Figure 3A). In addition, we also observed how the LPL molecule gradually embedded in the phospholipid membrane when it was initially located outside the membrane (Figure 3B). At 275.2 ns, a hydrophobic segment of the LPL molecule protruded into the hydrophobic region of the membrane, forming a mosaic state. At 300 ns, the conformation was largely the same as that when LPL was inserted 1.0 nm into the membrane—namely, one hydrophobic segment was embedded in the membrane and the other floated irregularly in the water. The results revealed that the entry of the LPL molecule into the cell membrane also depended on the hydrophobic interaction, similar to the case with L64.
When rL2PL2 was inserted into the membrane, it showed a similar phenomenon to L64. From 55.0 ns, the hydrophobic segment of rL2PL2 was embedded in the membrane, and the hydrophilic ends floated on the membrane surface (Figure 4A). However, unlike L64, rL2PL2 molecules did not have a clear tendency to move to the center of the membrane but were embedded in a relatively fixed position on the membrane, indicating that rL2PL2 could interact with both the hydrophobic and hydrophilic regions of the phospholipid membranes. The hydrophilic ends of rL2PL2 were also in contact with the external aqueous solution. Meanwhile, the membrane entry process of rL2PL2 was also observed when rL2PL2 was put on the surface of the membrane (Figure 4B). Similar to L64, when the hydrophobic part of rL2PL2 contacted the hydrophobic region of the membrane, rL2PL2 could gradually enter the center of the membrane, until the hydrophobic part was inside and the hydrophilic ends were outside the membrane.
When the rL3PL3 molecule was inserted into the phospholipid membrane, the phenomenon was similar to that of the rL2PL2. But the difference was that rL3PL3 did not exhibit a mosaic state similar to that of rL2PL2 until 70.0 ns. The rL3PL3 molecule embedded in a shallower place of the membrane and its hydrophilic segments tended to be clumped together (Figure 5).
To sum up, no matter what material, as long as the hydrophobic part comes into contact with the hydrophobic region of the membrane, the material will enter the membrane spontaneously and will not be repelled from the cell membrane.

3.2. Transmembrane Free Energy Analysis of the Materials

In this step, the transmembrane free energy of the four materials was analyzed. After the attachment of the materials to the cell membrane, the free energy sharply decreased (Figure 6), indicating that the four materials could easily attach and enter the phospholipid membrane. In summary, the energy changes were rL3PL3 > rL2PL2 = LPL > L64. However, when the materials left the membrane, in all cases the energy slowly increased to a varying degree, implying that the four materials molecules all had a tendency to be detained in the center of the cell membrane, which may be due to electric interaction between the material and membrane. It is obvious that rL2PL2 caused a higher energy change than LPL while exiting the membrane, implying that the LPL molecule could exit the membrane more easily.
In summary, the four materials showed the same characteristics in that the free energy reduced while entering the membrane and increased while exiting the membrane. These results also explained that when the material molecules were inserted 1.0 nm into the membrane, they were not repelled but spontaneously moved toward the center of the membrane.
We found that after entering the cell membrane, L64, LPL, and rL2PL2 would change from their initial shape to a linear or dumbbell shape, while rL3PL3 still kept the initial clustering structure (Figure 2, Figure 3, Figure 4 and Figure 5). Meanwhile, the final morphologies of the four molecules were also confirmed in the PMF process (Figure S3). According to the molecular structure, we speculated that the conformational change may be caused by intramolecular hydrogen bonds. To address this, we analyzed the number of intramolecular hydrogen bonds in the four molecules. The results showed that the rL3PL3 molecule contained about 13 intramolecular hydrogen bonds, while other materials had fewer than 10 bonds (Figure 7). This meant that rL3PL3 had stronger intramolecular interaction than the other molecules, so the rL3PL3 molecule could not easily undergo a conformational change but was inclined to keep its initial conformation. The existence of excessive intramolecular hydrogen bonds would reduce the flexibility of the conformational change of the molecule, thus finally weakening the disturbance effects. The results of hydrogen bond analysis were also consistent with our previous research, namely, rL3PL3 had the worst gene delivery results among the four molecules [2,22,23].

3.3. Parametric Analysis of Phospholipid Membranes

To further understand the mechanism of these four materials in gene delivery, we analyzed the changes in deuterium order parameters and the cell membrane transverse diffusion coefficient and compared the perturbation effects of the two final conformations of LPL on cell membranes. The results showed that L64, LPL, and rL2PL2 had potent disturbance effects on the cell membrane, but rL3PL3 had just a little (Figure 8A). This may be the reason why rL3PL3 did not function as well as other materials in previous intramuscular gene delivery studies [2,22,23]. At the same time, there was little difference in the disturbance effects of “LPL in” and “LPL on” (Figure 8A,B), indicating that the contact area with the phospholipid membrane did not affect the disturbance effects of LPL.
We found that the lateral diffusion coefficients (Msd) of these materials did not decrease. Even the lowest rL2PL2 group had the same lateral diffusion coefficient as the “None” group (Figure 8C). This indicated that when these materials were used for gene delivery, they would not affect the fluidity of the cell membrane. These results were considered to explain the reason why the materials had almost no cell membrane toxicity [2,22,23].
In order to study the contact details between material molecules and cell membranes, we calculated the atom number in the phospholipid membrane in contact with the material molecules, which can reflect the contact area between the material molecules and the membrane. The results showed that for the “hydrophilic-hydrophobic-hydrophilic” molecules, the number of contacted atoms ran rL2PL2 > L64 > rL3PL3 (Figure 9). This was completely consistent with the disturbance effects (Figure 8A,B) and gene delivery efficiency data of the three molecules in our previous research [2,22,23]. Therefore, we found that for such materials, the number of contacted atoms between the material and the membrane could directly affect the disturbance effect and consequently, the intramuscular gene delivery efficiency. Conversely, as a “hydrophobic-hydrophilic-hydrophobic” triblock material, the “LPL on” had obviously fewer atomic contacts than “LPL in”, but both the “LPL on” and “LPL in” had the same disturbance effects (Figure 8A,B). This means that the delivery mechanism of LPL may be different from that of “hydrophilic-hydrophobic-hydrophilic” materials L64, rL2PL2, and rL3PL3.

3.4. The Interaction between DNA Molecule and Material Molecule

To analyze the interaction between a DNA molecule and a material molecule, a model system was built. In this model, the initial centroid distance was 3.0 nm and the lattice was 14.0 nm × 8.0 nm × 8.0 nm. The same pre-equilibration steps and parameters were used for the material and the membrane, and the simulation time was 100.0 ns. According to the real-time minimum distance map, no contact was displayed between the material molecule and the DNA molecule for most of the duration of the simulation (Figure 10). Therefore, it could be concluded that there was basically no interaction between the materials and DNA. These results were consistent with previous gel retardation assays [2,22,23].

3.5. Transmembrane Free Energy of DNA Molecule

In our previous experimental studies, it was confirmed that in the presence of material molecules, the delivery efficiency of DNA was higher than that of DNA alone [2,22,23]. Here, a free energy analysis was performed to explore the mechanism for the improved gene delivery efficiency. In this step, material molecules were placed in the membrane to perturb the membrane in advance for 100.0 ns, and then the transmembrane free energy of DNA was calculated. Only the molecules of L64 and LPL were investigated in this simulation since they had opposite amphiphilic structures.
For cell membranes without material perturbation, the transmembrane energy barrier to DNA was 97.3 kJ/mol. After L64 perturbation, the transmembrane energy barrier decreased to 85.6 kJ/mol, and after LPL perturbation, it decreased to 57.2 kJ/mol. However, when DNA encountered the LPL molecule in the membrane, the energy barrier became 85.2 kJ/mol (Figure 11). The results showed that after LPL perturbed the cell membrane, the DNA transmembrane energy barrier would significantly decrease (97.3 kJ/mol to 57.2 kJ/mol), which consequently made it easier for DNA to penetrate the membrane. It was also clear that in terms of lowering the DNA transmembrane energy barrier, LPL was better than L64 (57.2 kJ/mol vs. 85.6 kJ/mol). The results were still consistent with previous experimental results regarding intramuscular gene delivery [22]. In addition, as the DNA transmembrane energy barrier would be negatively influenced while encountering the LPL molecules in the membrane (57.2 kJ/mol vs. 85.2 kJ/mol), it would be better to reduce the retention of material molecules in the membrane.

4. Discussion

Traditional cationic polymer gene delivery materials, such as polyethyleneimine (PEI), usually electrostatically interact with DNA molecules and condense them into a more compact structure, called, for example, the PEI-DNA complex. When DNA is compressed, it can be protected from being degraded by nucleases outside and inside target cells, and there is an interaction between the DNA and material molecules throughout the action process [67,68]. In contrast, the nonionic triblock copolymer materials we studied here showed almost no interaction with DNA during the entire interaction, and contact between them was absent most of the time (Figure 10). These characteristics have been proved by DNA retardation assays in previous studies [2,22,23]. Therefore, this kind of material and DNA act independently and successfully during the transmembrane process.
Based on the simulation results, the interaction between material molecules and DNA molecules, the process of material molecules contacting and entering the membrane, and the change in DNA transmembrane free energy were systematically explored and displayed. The results showed that the four material molecules all had a tendency to interact and enter the membrane, and to be embedded in the phospholipid membrane at a certain time (Figure 2, Figure 3, Figure 4 and Figure 5). The analysis of the free energy showed that this process was accompanied by a decrease in free energy, revealing that the materials could enter the cell membrane spontaneously (Figure 6). When the hydrophilic ends of the materials faced the phospholipid membrane, the free energy was reduced to a lower value (Figure 6B), so, it could be inferred that the hydrophilic end of a material molecule should first contact and enter the phospholipid membrane. Subsequently, both the “hydrophilic-hydrophobic-hydrophilic” and “hydrophobic-hydrophilic-hydrophobic” molecules utilized their hydrophobic ends to contact with the hydrophobic area of the membrane and then enter the membrane through the hydrophobic interaction (Figure 2, Figure 3, Figure 4 and Figure 5). When the distance between the material molecule and the mass center of the membrane was zero, the free energy reached its lowest value (Figure 6), which also demonstrated that the free energy change process is spontaneous and energy-independent.
For the “hydrophilic-hydrophobic-hydrophilic” copolymer molecules, including L64, rL2PL2, and rL3PL3, it was found that their perturbation effects on the cell membrane depended on the number of contacted phospholipid atoms (Figure 9). The number of phospholipid atoms contacted by the material molecules ran rL2PL2 > L64 > rL3PL3 (Figure 9). Previous experiments revealed that rL2PL2 had the best in vivo gene delivery and expression outcomes, but L64 had almost the same results as those of rL3PL3 [23]. This may be attributed to the following reasons. L64 almost has no intramolecular hydrogen bonds (Figure 7), so it can be fully stretched in the cell membrane and, therefore, has a large contact area. However, the molecular weight and size of L64 are relatively small, thus limiting its perturbation effect on the membrane. Conversely, rL3PL3 is the largest molecule, but it has more internal hydrogen bonds than the other molecules (Figure 7), which may limit the ability of the molecule to change its confirmation, thus reducing the contact with membrane atoms and weakening the perturbation effect on the membrane. Meanwhile, the large size of the rL3PL3 molecule can also make it more difficult for it to enter the cell membrane, thus remaining at a shallower position than rL2PL2 (Figure 5).
From the intramuscular gene delivery and expression experiments, we already knew that LPL was more effective than L64, and rL2PL2 was better than L64 and rL3PL3 [22,23]. However, we still do not know which is better, LPL or rL2PL2. Further, we need to make clear which kind of molecular structure is better, “hydrophobic-hydrophilic-hydrophobic” or “hydrophilic-hydrophobic-hydrophilic”. This is of special importance for further intramuscular gene delivery material design and optimization. When the DNA molecule encountered the LPL molecule in the membrane, the DNA transmembrane energy barrier would significantly increase again (Figure 11E), which made it difficult for DNA to penetrate the membrane. So, if the material molecules found it easier to leave the cell membrane than to remain within, there should be less opportunity for the DNA molecules to encounter the material molecules. However, the free energy of the material molecules was found to rise again when the molecules left the membrane (Figure 6A), indicating that there is resistance to them leaving the membrane. So, if a material molecule had a lower free energy change while exiting the membrane, it would be easier for the material molecule to leave the membrane. Here, it was clear that LPL had an obviously lower free energy change than rL2PL2 (Figure 6A) while exiting the membrane, so, LPL found it easier than rL2PL2 to leave the membrane. Moreover, as LPL was found to embed in a relatively shallow position in the cell membrane (Figure 3A), it was able to produce a disturbance effect faster than rL2PL2. In addition, for the “hydrophobic-hydrophilic-hydrophobic” copolymer of LPL, whether it was “LPL in” or “LPL on” state, the perturbation effect was largely the same (Figure 8A,B). The number of contacted phospholipid atoms had little effect on the perturbation effect. Considering that LPL has an opposite molecular structure to that of the cell membrane and LPL finds it easier to leave the membrane, we speculate that LPL molecules are more likely to pass through the cell membrane, i.e., the disturbance results mainly rely on the transmembrane action. This is quite different from the case with “hydrophilic-hydrophobic-hydrophilic” molecules, which mainly depend on the intermembrane perturbation effect via interaction with contacted phospholipid atoms.

5. Conclusions

In this study, we revealed the details of several typical and efficient intramuscular gene delivery materials at the atomic level, including the process involved in the material molecules entering and exiting the phospholipid membrane, interactions between material molecules and the membrane, interactions between material molecules and DNA molecules, and also including the free energy changes in the materials, and the changes in the transmembrane energy barrier to DNA molecules. More importantly, the simulation results almost completely matched our previous experimental results, thus proving the rationality of the molecular design of these materials. Regardless of whether the structure was “hydrophilic-hydrophobic-hydrophilic” or “hydrophobic-hydrophilic-hydrophobic”, the material molecules approached the cell membrane with their hydrophilic ends at first and then gradually entered the interior of the membrane through the interaction between their hydrophobic segments with the hydrophobic region of the membrane. They had no adverse effect on the fluidity of the cell membranes, indicating their biocompatibility with skeletal muscle cells. By analyzing the energy change, the spontaneity of the material molecules entering the cell membrane was proved. Upon the perturbation effects of material molecules, the transmembrane free energy barrier of DNA was decreased, making it easier for DNA to pass through the cell membrane. The simulation results also showed that there was little-to-no interaction between DNA and material molecules.
In addition, this study revealed and explained that the “hydrophobic-hydrophilic-hydrophobic” structure of LPL had several advantages over the “hydrophilic-hydrophobic-hydrophilic” structure of rL2PL2. The results of this study will help us to design and optimize intramuscular gene delivery materials.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jfb14040219/s1, Table S1: Details of the various systems simulated in this study; Figure S1: Initial model of four-groups; (A) System A, B, D, E; (B) System F-J; (C) System K-N; (D) System O-R; Figure S2: The process of L64 entering the phospholipid membrane (the model of 2 membranes); Figure S3: Final conformations of four molecules in PMF across phospholipid membranes: (a) L64; (b)LPL; (c) rL2PL2; (d) rL3PL3.

Author Contributions

The manuscript was written through contributions of all authors. G.W. and Z.C. designed the project. Z.C. and Y.J. conducted the project with assistance from L.P., J.C. and M.L., Z.C. and Y.J. processed the data. Z.C., Y.J. and G.W. wrote the manuscript. G.W. and J.Z.T. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31971390), the 1.3.5 Project for Disciplines of Excellence, West China Hospital, Sichuan University (Grant No. ZYJC21043), and the International Cooperative Project of Sichuan Province on Science and Technology Innovation (No. 2021YFH0142).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data used in this article has been presented in the main text and Supporting Information.

Acknowledgments

The authors would thank the National Natural Science Foundation of China (31971390), the 1.3.5 Project for Disciplines of Excellence, West China Hospital, Sichuan University (Grant No. ZYJC21043), and the International Cooperative Project of Sichuan Province on Science and Technology Innovation (No. 2021YFH0142).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wolff, J.A.; Malone, R.W.; Williams, P.; Chong, W.; Acsadi, G.; Jani, A.; Felgner, P.L. Direct Gene Transfer into Mouse Muscle In Vivo. Science 1990, 247, 1465–1468. [Google Scholar] [CrossRef] [PubMed]
  2. Song, H.; Liu, S.; Li, C.; Geng, Y.; Wang, G.; Gu, Z. Pluronic L64-mediated stable HIF-1α expression in muscle for therapeutic angiogenesis in mouse ischemic limb. Int. J. Nanomed. 2014, 9, 3439–3452. [Google Scholar] [CrossRef] [Green Version]
  3. Gao, X.; Kim, K.-S.; Liu, D. Nonviral Gene Delivery: What We Know and What Is Next. AAPS J. 2007, 9, E92–E104. [Google Scholar] [CrossRef] [PubMed]
  4. Bosman, F.T.; Stamenkovic, I. Functional structure and composition of the extracellular matrix. J. Pathol. 2003, 200, 423–428. [Google Scholar] [CrossRef] [PubMed]
  5. Conwell, C.-C.; Huang, L. Recent advances in non-viral gene delivery. Adv. Genet. 2005, 53, 3–18. [Google Scholar] [CrossRef]
  6. Wang, W.; Li, W.; Ma, N.; Steinhoff, G. Non-viral gene delivery methods. Curr. Pharm. Biotechnol. 2013, 14, 46–60. [Google Scholar] [CrossRef]
  7. Aied, A.; Greiser, U.; Pandit, A.; Wang, W. Polymer gene delivery: Overcoming the obstacles. Drug Discov. Today. 2013, 18, 1090–1098. [Google Scholar] [CrossRef]
  8. Ramamoorth, M.; Narvekar, A. Non-viral vectors in gene therapy-an overview. J. Clin. Diagn. Res. 2015, 9, GE01–GE06. [Google Scholar] [CrossRef]
  9. Wirth, T.; Parker, N.; Ylä-Herttuala, S. History of gene therapy. Gene 2013, 525, 162–169. [Google Scholar] [CrossRef]
  10. Mumper, R.J.; Duguid, J.G.; Anwer, K.; Barron, M.K.; Nitta, H.; Rolland, A.P. Polyvinyl derivatives as novel interactive polymers for controlled gene delivery to muscle. Pharm. Res. 1996, 13, 701–709. [Google Scholar] [CrossRef]
  11. Osada, K.; Shiotani, T.; Tockary, T.A.; Kobayashi, D.; Oshima, H.; Ikeda, S.; Christie, R.J.; Itaka, K.; Kataoka, K. Enhanced gene expression promoted by the quantized folding of pDNA within polyplex micelles. Biomaterials 2012, 33, 325–332. [Google Scholar] [CrossRef]
  12. Qi, R.; Gao, Y.; Tang, Y.; He, R.-R.; Liu, T.-L.; He, Y.; Sun, S.; Li, B.-Y.; Li, Y.-B.; Liu, G. PEG-conjugated PAMAM Dendrimers Mediate Efficient Intramuscular Gene Expression. AAPS J. 2009, 11, 395–405. [Google Scholar] [CrossRef] [Green Version]
  13. Chang, C.-W.; Choi, D.; Kim, W.J.; Yockman, J.W.; Christen, L.V.; Kim, Y.-H.; Kim, S.W. Non-ionic amphiphilic biodegradable PEG–PLGA–PEG copolymer enhances gene delivery efficiency in rat skeletal muscle. J. Control. Release 2007, 118, 245–253. [Google Scholar] [CrossRef]
  14. Kabanov, A.V.; Batrakova, E.V.; Alakhov, V.Y. Pluronic® block copolymers for overcoming drug resistance in cancer. Adv. Drug Deliv. Rev. 2002, 54, 759–779. [Google Scholar] [CrossRef]
  15. Pitard, B.; Pollard, H.; Agbulut, O.; Lambert, O.; Vilquin, J.T.; Cherel, Y.; Abadie, J.; Samuel, J.L.; Rigaud, J.L.; Menoret, S.; et al. A nonionic amphiphile agent promotes gene delivery in vivo to skeletal and cardiac muscles. Hum. Gene Ther. 2002, 13, 1767–1775. [Google Scholar] [CrossRef]
  16. Lavigne, M.D.; Pohlschmidt, M.; Novo, J.F.; Higgins, B.; Alakhov, V.; Lochmuller, H.; Sakuraba, H.; Goldspink, G.; MacDermot, K.; Górecki, D.C. Promoter dependence of plasmid-pluronics targeted alpha galactosidase A expression in skeletal muscle of Fabry mice. Mol. Ther. 2005, 12, 985–990. [Google Scholar] [CrossRef]
  17. Yang, Z.; Zhu, J.; Sriadibhatla, S.; Gebhart, C.; Alakhov, V.; Kabanov, A. Promoter- and strain-selective enhancement of gene expression in a mouse skeletal muscle by a polymer excipient Pluronic P85. J. Control. Release 2005, 108, 496–512. [Google Scholar] [CrossRef]
  18. Lui, S.; Ma, L.; Tan, R.; Lu, Q.; Geng, Y.; Wang, G.; Gu, Z. Safe and efficient local gene delivery into skeletal muscle via a combination of Pluronic L64 and modified electrotransfer. Gene Ther. 2014, 21, 558–565. [Google Scholar] [CrossRef] [Green Version]
  19. He, Y.; Liu, Y.; Sun, Z.; Han, F.; Tang, J.Z.; Gao, R.; Wang, G. The proper strategy to compress and protect plasmid DNA in the Pluronic L64-electropulse system for enhanced intramuscular gene delivery. Regen. Biomater. 2019, 6, 289–298. [Google Scholar] [CrossRef]
  20. Ma, L.; Wang, G.; Liu, S.; Bi, F.; Liu, M.; Wang, G. Intramuscular expression of plasmid-encoded FVII-Fc immunoconjugate for tumor immunotherapy by targeting tumoral blood vessels and cells. Front. Oncol. 2021, 11, 638591. [Google Scholar] [CrossRef]
  21. Deng, L.; Yang, P.; Li, C.; Xie, L.; Lu, W.; Zhang, Y.; Liu, M.; Wang, G. Prolonged control of insulin-dependent diabetes via intramuscular expression of plasmid-encoded single-strand insulin analogue. Genes Dis. 2022; in press. [Google Scholar] [CrossRef]
  22. Pu, L.; Geng, Y.; Liu, S.; Chen, J.; Luo, K.; Wang, G.; Gu, Z. Electroneutralized Amphiphilic Triblock Copolymer with a Peptide Dendron for Efficient Muscular Gene Delivery. ACS Appl. Mater. Interfaces 2014, 6, 15344–15351. [Google Scholar] [CrossRef] [PubMed]
  23. Pu, L.; Wang, J.; Li, N.; Chai, Q.; Irache, J.M.; Wang, G.; Tang, J.Z.; Gu, Z. Synthesis of Electroneutralized Amphiphilic Copolymers with Peptide Dendrons for Intramuscular Gene Delivery. ACS Appl. Mater. Interfaces 2016, 8, 13724–13734. [Google Scholar] [CrossRef] [PubMed]
  24. Cui, Z.; Jiao, Y.; Pu, L.; Tang, J.Z.; Wang, G. The Progress of Non-Viral Materials and Methods for Gene Delivery to Skeletal Muscle. Pharmaceutics 2022, 14, 2428. [Google Scholar] [CrossRef]
  25. Chen, J.; Luo, J.; Zhao, Y.; Pu, L.; Lu, X.; Gao, R.; Wang, G.; Gu, Z. Increase in Transgene Expression by Pluronic L64-Mediated Endosomal/Lysosomal Escape through Its Membrane-Disturbing Action. ACS Appl. Mater. Interfaces 2015, 7, 7282–7293. [Google Scholar] [CrossRef]
  26. Salsbury, F.R., Jr. Molecular dynamics simulations of protein dynamics and their relevance to drug discovery. Curr. Opin. Pharmacol. 2010, 10, 738–744. [Google Scholar] [CrossRef] [Green Version]
  27. Dallavalle, M.; Calvaresi, M.; Bottoni, A.; Melle-Franco, M.; Zerbetto, F. Graphene can wreak havoc with cell membranes. ACS Appl. Mater. Interfaces 2015, 7, 4406–4414. [Google Scholar] [CrossRef]
  28. Liu, L.; Zhang, S.; Zhao, L.; Gu, Z.; Duan, G.; Zhou, B.; Yang, Z.; Zhou, R. Superior Compatibility of C2N with Human Red Blood Cell Membranes and the Underlying Mechanism. Small 2018, 14, e180539. [Google Scholar] [CrossRef]
  29. Gao, X.; Hong, S.; Liu, Z.; Yue, T.; Dobnikar, J.; Zhang, X. Membrane potential drives direct translocation of cell-penetrating peptides. Nanoscale 2019, 11, 1949–1958. [Google Scholar] [CrossRef]
  30. Tian, L.; Wu, G. Microsecond molecular dynamics simulation of the adsorption and penetration of oil droplets on cellular membrane. J. Hazard. Mater. 2020, 397, 122683. [Google Scholar] [CrossRef]
  31. Harvey, M.J.; Fabritiis, G.D. High-throughput molecular dynamics the powerful new tool for drug discovery. Drug Discov. Today 2012, 17, 1059–1062. [Google Scholar] [CrossRef]
  32. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef] [Green Version]
  33. Hildebrand, P.W.; Rose, A.S.; Tiemann, J.K.S. Bringing Molecular Dynamics Simulation Data into View. Trends Biochem. Sci. 2019, 44, 902–913. [Google Scholar] [CrossRef] [Green Version]
  34. He, X.; Lin, M.; Lu, T.; Qu, Z.; Xu, F. Molecular analysis of interactions between a PAMAM dendrimer-paclitaxel conjugate and a biomembrane. Phys. Chem. Chem. Phys. 2015, 17, 29507–29517. [Google Scholar] [CrossRef]
  35. Voulgarakis, N.-K.; Rasmussen, K.-Ø.; Welch, P.-M. Dendrimers as synthetic gene vectors: Cell membrane attachment. J. Chem. Phys. 2009, 130, 155101. [Google Scholar] [CrossRef]
  36. Drenscko, M.; Loverde, S.-M. Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone. Mol. Simul. 2019, 45, 859–867. [Google Scholar] [CrossRef]
  37. Sun, C.; Tang, T. Study on the role of polyethylenimine as gene delivery carrier using molecular dynamics simulations. J. Adhes. Sci. Technol. 2012, 28, 399–416. [Google Scholar] [CrossRef]
  38. Arvayo-Zatarain, J.A.; Favela-Rosales, F.; Contreras-Aburto, C.; Urrutia-Bañuelos, E.; Maldonado, A. Molecular dynamics simulation study of the effect of halothane on mixed DPPC/DPPE phospholipid membranes. J. Mol. Model. 2019, 25, 4. [Google Scholar] [CrossRef]
  39. Sasaki, E.; Hayashi, Y.; Kimura, Y.; Sashida, S.; Hamano, N.; Nirasawa, K.; Hamada, K.; Katagiri, F.; Kikkawa, Y.; Sakai, T.; et al. Alpha-dystroglycan binding peptide A2G80-modified stealth liposomes as a muscle-targeting carrier for Duchenne muscular dystrophy. J. Control. Release 2021, 329, 1037–1045. [Google Scholar] [CrossRef]
  40. Nirasawa, K.; Hamada, K.; Naraki, Y.; Kikkawa, Y.; Sasaki, E.; Endo-Takahashi, Y.; Hamano, N.; Katagiri, F.; Nomizu, M.; Negishi, Y. Development of A2G80 peptide-gene complex for targeted delivery to muscle cells. J. Control. Release 2021, 329, 988–996. [Google Scholar] [CrossRef]
  41. Pomel, C.; Leborgne, C.; Cheradame, H.; Scherman, D.; Kichler, A.; Guegan, P. Synthesis and evaluation of amphiphilic poly(tetrahydrofuran-b-ethylene oxide) copolymers for DNA delivery into skeletal muscle. Pharm. Res. 2008, 25, 2963–2971. [Google Scholar] [CrossRef] [PubMed]
  42. Rasolonjatovo, B.; Illy, N.; Bennevault, V.; Mathé, J.; Midoux, P.; Le Gall, T.; Haudebourg, T.; Montier, T.; Lehn, P.; Pitard, B.; et al. Temperature-Sensitive Amphiphilic Non-Ionic Triblock Copolymers for Enhanced In Vivo Skeletal Muscle Transfection. Macromol. Biosci. 2020, 20, e1900276. [Google Scholar] [CrossRef] [PubMed]
  43. Brooks, B.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar] [CrossRef]
  44. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.W.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D., Jr.; Pastor, R.W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  46. Spoel, D.V.D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  47. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  48. Kutzner, C.; Páll, S.; Fechner, M.; Esztermann, A.; Groot, B.L.D.; Grubmüller, H. More Bang for Your Buck: Improved use of GPU Nodes for GROMACS 2018. J. Comput. Chem. 2019, 40, 2418–2431. [Google Scholar] [CrossRef] [Green Version]
  49. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  50. MacKerell, A.D., Jr.; Vanommeslaeghe, K. Automation of the CHARMM General Force Field (CGenFF) I: Bond Perception and Atom Typing. J. Chem. Inf. Model. 2012, 52, 3144–3154. [Google Scholar] [CrossRef]
  51. Vanommeslaeghe, K.; Raman, E.P.; MacKerell, A.D., Jr. Automation of the CHARMM General Force Field (CGenFF) II: Assignment of Bonded Parameters and Partial Atomic Charges. J. Chem. Inf. Model. 2012, 52, 3155–3168. [Google Scholar] [CrossRef] [Green Version]
  52. Yu, W.; He, X.; Vanommeslaeghe, K.; MacKerell, A.D., Jr. Extension of the CHARMM General Force Field to Sulfonyl-Containing Compounds and Its Utility in Biomolecular Simulations. J. Comput. Chem. 2012, 33, 2451–2468. [Google Scholar] [CrossRef] [Green Version]
  53. Avogadro: An Open-Source Molecular Builder and Visualization Tool, Version 1.2.0; Available online: http://avogadro.cc/ (accessed on 13 August 2012).
  54. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeersch, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [Green Version]
  55. Tieleman, D.P.; Forrest, L.R.; Sansom, S.P.; Berendsen, H.J.C. Lipid Properties and the Orientation of Aromatic Residues in OmpF, Influenza M2, and Alamethicin Systems: Molecular Dynamics Simulations. Biochemistry 1998, 37, 17554–17561. [Google Scholar] [CrossRef]
  56. Mukhopadhyay, P.; Monticelli, L.; Tieleman, D.P. Molecular Dynamics Simulation of a Palmitoyl-Oleoyl Phosphatidylserine Bilayer with Na+ Counterions and NaCl. Biophys. J. 2004, 86, 1601–1609. [Google Scholar] [CrossRef] [Green Version]
  57. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  58. Jo, S.; Lim, J.B.; Klauda, J.B.; Im, W. CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes. Biophys. J. 2009, 97, 50–58. [Google Scholar] [CrossRef] [Green Version]
  59. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.C.; Dávila-Contreras, E.M.; Qi, Y.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder Toward Realistic Biological Membrane Simulations. J. Comput. Chem. 2014, 35, 1997–2004. [Google Scholar] [CrossRef] [Green Version]
  60. Berendsen, H.J.C.; Postma, J.P.M.; Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  61. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [Green Version]
  62. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  63. Roux, B. The calculation of the potential of mean force using computer simulations. Comput. Phys. Comm. 1995, 91, 275–282. [Google Scholar] [CrossRef]
  64. Kumar, S.; Rosenerg, J.M.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A. Multidimensional Free-Energy Calculations Using the Weighted Histogram Analysis Method. J. Comput. Chem. 1995, 16, 1339–1350. [Google Scholar] [CrossRef]
  65. Kästner, J.; Thiel, W. Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: “Umbrella integration”. J. Chem. Phys. 2005, 123, 144104. [Google Scholar] [CrossRef]
  66. Hub, J.S.; Groot, B.L.; Spoel, D. g_whams-A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. J. Chem. Theory Comput. 2010, 6, 3713–3720. [Google Scholar] [CrossRef] [Green Version]
  67. Sirsi, S.R.; Schray, R.C.; Guan, X.; Lykens, N.M.; Williams, J.H.; Erney, M.L.; Lutz, G.J. Functionalized PEG–PEI Copolymers Complexed to Exon-Skipping Oligonucleotides Improve Dystrophin Expression in mdx Mice. Hum. Gene Ther. 2008, 19, 795–806. [Google Scholar] [CrossRef]
  68. Lee, J.-L.; Lo, C.-W.; Ka, S.-M.; Chen, A.; Chen, W.-S. Prolonging the expression duration of ultrasound-mediated gene transfection using PEI nanoparticles. J. Control. Release 2012, 160, 64–71. [Google Scholar] [CrossRef]
Jfb 14 00219 g001 550
Figure 1. Structures of various gene delivery polymers: (A) Pluronic L64; (B) LPL; (C) rL2PL2; (D) rL3PL3; and cell membrane molecules: (E) POPC; (F) POPS.
Figure 1. Structures of various gene delivery polymers: (A) Pluronic L64; (B) LPL; (C) rL2PL2; (D) rL3PL3; and cell membrane molecules: (E) POPC; (F) POPS.
Jfb 14 00219 g001
Jfb 14 00219 sch001 550
Scheme 1. The process flow of the study simulations.
Scheme 1. The process flow of the study simulations.
Jfb 14 00219 sch001
Jfb 14 00219 g002 550
Figure 2. The process of L64 entering the phospholipid membrane. L64 was inserted about 1.0 nm into the phospholipid membrane (initial state). Red: hydrophobic segment; blue: hydrophilic segments.
Figure 2. The process of L64 entering the phospholipid membrane. L64 was inserted about 1.0 nm into the phospholipid membrane (initial state). Red: hydrophobic segment; blue: hydrophilic segments.
Jfb 14 00219 g002
Jfb 14 00219 g003 550
Figure 3. The process of LPL entering the phospholipid membrane. (A) LPL was inserted 1.0 nm into the phospholipid membrane (initial state). (B) LPL was out of the membrane (initial state). Red: hydrophobic segments; blue: hydrophilic segments.
Figure 3. The process of LPL entering the phospholipid membrane. (A) LPL was inserted 1.0 nm into the phospholipid membrane (initial state). (B) LPL was out of the membrane (initial state). Red: hydrophobic segments; blue: hydrophilic segments.
Jfb 14 00219 g003
Jfb 14 00219 g004 550
Figure 4. The process of rL2PL2 entering the phospholipid membrane. (A) rL2PL2 was inserted 1.0 nm into the phospholipid membrane (initial state). (B) rL2PL2 was out of the membrane (initial state). Red: hydrophobic segment; blue: hydrophilic segments.
Figure 4. The process of rL2PL2 entering the phospholipid membrane. (A) rL2PL2 was inserted 1.0 nm into the phospholipid membrane (initial state). (B) rL2PL2 was out of the membrane (initial state). Red: hydrophobic segment; blue: hydrophilic segments.
Jfb 14 00219 g004
Jfb 14 00219 g005 550
Figure 5. The process of rL3PL3 entering the phospholipid membrane. rL3PL3 was inserted 1.0 nm into the phospholipid membrane. Red: hydrophobic segment; blue: hydrophilic segments.
Figure 5. The process of rL3PL3 entering the phospholipid membrane. rL3PL3 was inserted 1.0 nm into the phospholipid membrane. Red: hydrophobic segment; blue: hydrophilic segments.
Jfb 14 00219 g005
Jfb 14 00219 g006 550
Figure 6. Changes in transmembrane free energy of the four molecules. The phospholipid membrane was located at coordinates 2.0–6.0 nm. (A) L64, LPL rL2PL2, and rL3PL3. (B) rL3PL3 and rL3PL3-2. rL3PL3: hydrophobic end of rL3PL3 facing the membrane; rL3PL3−2: hydrophilic ends of rL3PL3 facing the membrane.
Figure 6. Changes in transmembrane free energy of the four molecules. The phospholipid membrane was located at coordinates 2.0–6.0 nm. (A) L64, LPL rL2PL2, and rL3PL3. (B) rL3PL3 and rL3PL3-2. rL3PL3: hydrophobic end of rL3PL3 facing the membrane; rL3PL3−2: hydrophilic ends of rL3PL3 facing the membrane.
Jfb 14 00219 g006
Jfb 14 00219 g007 550
Figure 7. Analysis of the number of intramolecular hydrogen bonds.
Figure 7. Analysis of the number of intramolecular hydrogen bonds.
Jfb 14 00219 g007
Jfb 14 00219 g008 550
Figure 8. Parameter analysis of the phospholipid membrane after interaction with the four molecules respectively. LPL in: the entire LPL molecule embedded in the membrane. LPL on: half of the LPL embedded in the membrane. Lipid: the membrane alone, without any gene delivery molecule. (A) The deuterated order parameter of the sn-1 chain; (B) the deuterated order parameter of the sn-2 chain; (C) the lateral diffusion coefficient Msd of the membrane. None: the membrane without any gene delivery molecule.
Figure 8. Parameter analysis of the phospholipid membrane after interaction with the four molecules respectively. LPL in: the entire LPL molecule embedded in the membrane. LPL on: half of the LPL embedded in the membrane. Lipid: the membrane alone, without any gene delivery molecule. (A) The deuterated order parameter of the sn-1 chain; (B) the deuterated order parameter of the sn-2 chain; (C) the lateral diffusion coefficient Msd of the membrane. None: the membrane without any gene delivery molecule.
Jfb 14 00219 g008
Jfb 14 00219 g009 550
Figure 9. Atom numbers in the phospholipid membrane in contact with different material molecules. LPL in: the whole LPL molecule was embedded in the membrane. LPL on: LPL was only half embedded in the membrane.
Figure 9. Atom numbers in the phospholipid membrane in contact with different material molecules. LPL in: the whole LPL molecule was embedded in the membrane. LPL on: LPL was only half embedded in the membrane.
Jfb 14 00219 g009
Jfb 14 00219 g010 550
Figure 10. Minimum distance change between the DNA molecule and the material molecule.
Figure 10. Minimum distance change between the DNA molecule and the material molecule.
Jfb 14 00219 g010
Jfb 14 00219 g011 550
Figure 11. Changes in DNA transmembrane free energy. The phospholipid membrane was located at coordinates 4.0–8.0 nm. (A) across pure membranes; (B) across L64-perturbed membranes (DNAs did not encounter L64 molecules); (C) across LPL-perturbed membranes (DNA encountered LPL molecules); (D) across LPL-perturbed membranes (DNA did not encounter LPL molecules); (E) Changes in DNA transmembrane free energy of (A–D).
Figure 11. Changes in DNA transmembrane free energy. The phospholipid membrane was located at coordinates 4.0–8.0 nm. (A) across pure membranes; (B) across L64-perturbed membranes (DNAs did not encounter L64 molecules); (C) across LPL-perturbed membranes (DNA encountered LPL molecules); (D) across LPL-perturbed membranes (DNA did not encounter LPL molecules); (E) Changes in DNA transmembrane free energy of (A–D).
Jfb 14 00219 g011
Table
Table 1. Experimental parameters of gene delivery polymers.
Table 1. Experimental parameters of gene delivery polymers.
Molar MassCritical Micelle Concentration (mg/mL)Work Concentration (wt%)
L6429000.430.1
LPL37750.0890.0075
rL2PL240250.0560.005
rL3PL359800.2510.02
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cui, Z.; Jiao, Y.; Pu, L.; Chen, J.; Liu, M.; Tang, J.Z.; Wang, G. The Interaction Mechanism of Intramuscular Gene Delivery Materials with Cell Membranes. J. Funct. Biomater. 2023, 14, 219. https://doi.org/10.3390/jfb14040219

AMA Style

Cui Z, Jiao Y, Pu L, Chen J, Liu M, Tang JZ, Wang G. The Interaction Mechanism of Intramuscular Gene Delivery Materials with Cell Membranes. Journal of Functional Biomaterials. 2023; 14(4):219. https://doi.org/10.3390/jfb14040219

Chicago/Turabian Style

Cui, Zhanpeng, Yang Jiao, Linyu Pu, Jianlin Chen, Ming Liu, James Zhenggui Tang, and Gang Wang. 2023. "The Interaction Mechanism of Intramuscular Gene Delivery Materials with Cell Membranes" Journal of Functional Biomaterials 14, no. 4: 219. https://doi.org/10.3390/jfb14040219

APA Style

Cui, Z., Jiao, Y., Pu, L., Chen, J., Liu, M., Tang, J. Z., & Wang, G. (2023). The Interaction Mechanism of Intramuscular Gene Delivery Materials with Cell Membranes. Journal of Functional Biomaterials, 14(4), 219. https://doi.org/10.3390/jfb14040219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop