Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone

The molecular interaction between common polymer chains and the cell membrane is unknown. Molecular dynamics simulations offer an emerging tool to characterise the nature of the interaction between common degradable polymer chains used in biomedical applications, such as polycaprolactone, and model cell membranes. Herein we characterise with all-atomistic and coarse-grained molecular dynamics simulations the interaction between single polycaprolactone chains of varying chain lengths with a phospholipid membrane. We find that the length of the polymer chain greatly affects the nature of interaction with the membrane, as well as the membrane properties. Furthermore, we next utilise advanced sampling techniques in molecular dynamics to characterise the two-dimensional free energy surface for the interaction of varying polymer chain lengths (short, intermediate, and long) with model cell membranes. We find that the free energy minimum shifts from the membrane-water interface to the hydrophobic core of the phospholipid membrane as a function of chain length. Finally, we perform coarse-grained molecular dynamics simulations of slightly larger membranes with polymers of the same length and characterise the results as compared with all-atomistic molecular dynamics simulations. These results can be used to design polymer chain lengths and chemistries to optimise their interaction with cell membranes at the molecular level.     

Analysis of brush-particle interactions using self-consistent-field theory

Non-specific protein adsorption can be reduced by attaching polymer chains by one end to a sorbent surface. End-grafted polymer modified surfaces have also found application in size-based chromatographic bioseparations. To better understand how to tailor surfaces for these applications, a numerical SCF model has been used to calculate theoretical results for the polymer density distribution of interacting polymer chains around a solute particle positioned at a fixed distance from a surface. In addition, the excess energy required to move the particle into the polymer chains (interaction energy) is calculated using a statistical mechanical treatment of the lattice model. The effect of system variables such as particle size, chain length, surface density and Flory interaction parameters on density distributions and interaction energies is also studied. Calculations for the interaction of a solute particle with a surface covered by many polymer chains (a brush) show that the polymer segments will fill in behind the particle quite rapidly as it moves toward the surface. When there is no strong energetic attraction between the polymer and solute we predict that the interaction energy will be purely repulsive upon compression due to losses in conformational entropy of the polymer chains. Above a critical chain length, which depends upon particle size, a maximum in the force required to move the particle toward the surface is observed due to an engulfment of the particle as chains attempt to access the free volume behind the particle. If an attraction exists between the polymer and solute, such that a minimum in the interaction energy is seen, the optimum conditions for solute repulsion occur at the highest surface density attainable. Long chain length can lead to increased solute concentration within the polymer layer due to the fact that an increased number of favourable polymer–solute contacts are able to occur than with short chains at a similar entropic penalty.     

Influence of polymer structure on adsorption behavior at solid–liquid interface by Monte Carlo simulation

Monte Carlo simulations for the adsorption of polymers including random copolymer, homopolymer, diblock copolymer and two kinds of triblock copolymers, respectively, in nonselective solvent at solid–liquid interface have been performed on a simple lattice model. The effect of polymer structure on adsorption properties was examined. In simulations, all polymeric molecules are modeled as self-avoiding linear chains composed of two segments A and B while A is attractive to the surface and B is non-attractive. It was found that for all polymers, the size distribution of various configurations is determined by the linked sequence of segments and the interaction energy between segment and surface. The results of simulation show that the adsorbed amount always increases with increasing bulk concentration but the adsorption layer thickness is mostly dependent on the adsorption energy at a fixed fraction of segments A. On the other hand, diblock copolymer has always the highest surface coverage and adsorbed amount, while random copolymers and homopolymers give generally the smallest surface coverage and adsorbed amount. It is shown that the sequence of polymer chains, i.e. molecular structure, is the most important factor in affecting adsorption properties at the same composition and interaction between segment and surface. The results also show that the adsorption behavior of random copolymers is remarkably different from that of block copolymers, but acting like homopolymer.     

Viscometric and spectroscopic studies on the solution behaviour of hydrophobically modified cellulosic polymers

The solution properties of hydroxyethyl cellulose (HEC) and hydrophobically modified hydroxyethyl cellulose (HM-HEC) have been investigated by means of viscometric and spectroscopic techniques involving free radical and fluorescent probes. The greater viscosity of HM-HEC solutions above a critical polymer concentration (C_p) of approximately 0·2% has been interpreted in terms of the formation of a three-dimensional network structure in which the polymer chains are effectively crosslinked by the intermolecular association of neighbouring hydrophobic side chains. C_p is considerably less than the predicted polymer coil overlap concentration (C*) of approximately 1%.

The interaction of the polymers with an anionic surfactant, sodium dodecyl sulphate (SDS) has also been investigated. A mechanism involving the interaction of free surfactant with the regions of intermolecular hydrophobic association is suggested to account for the considerable differences in the rheological behaviour of the polymers in the presence of SDS.     

Interactions of mono- and divalent counterions with alkali- and enzyme-deesterified pectins in salt-free solutions

The free fractions of monovalent and divalent counterions were determined on salt-free solutions of pectins. The effects of charge density, distribution of the carboxyl groups, polymer concentration, and the nature of the counterion were investigated by determinating the calcium and sodium activity coefficients (with specific electrodes) and by measuring the transport parameters (by conductimetry). Poor agreement for calcium ions was found with the Manning theory. The strong binding of these ions to highly charged polymers, which is ascribed to a dimerization process was demonstrated in very dilute solutions.     

Membrane affinity and antibacterial properties of cationic polyelectrolytes with different hydrophobicity

The antibacterial behavior of cationic polyelectrolytes is studied using model membrane experiments and in vitro bacterial investigations. The molecular interaction with lipid films is evaluated by the degree of penetration of the polymers into Langmuir monolayers of neutral or negatively charged lipids. The polymer/lipid interaction results in structural changes of the penetrated lipid layer visualized using AFM. The polymers are found to be effective in inhibiting the proliferation of E. coli, B. subtilis and S. aureus. The influence of the chemical structure on the functional behavior is related to the conformational properties. An optimum structure is identified on the basis of antibacterial and hemolytic tests as well as membrane-destroying efficacy of the antimicrobial polymers.     

Biophysical characterization of quaternary pyridinium functionalized polynorbornenes for DNA complexation and their cellular interactions

Cationic polymers with hydrophobic side chains have gained great interest as DNA carriers since they form a compact complex with negatively charged DNA phosphate groups and interact with the cell membrane. Amphiphilic polyoxanorbornenes with different quaternary alkyl pyridinium side chains with ethyl‐p(OPy2) and hexyl units‐p(OPy6) bearing 10 kDa MWT were synthesized by living Ring‐Opening Metathesis Polymerization method. The physicochemical characteristics: critical micellar concentration, size distribution, surface charge, and condensation of polymer/DNA complex were investigated. Morphology of complexes was monitored by Atomic force microscopy. Cytotoxicity and interaction of these complexes with model lipid vesicles mimicking the cell membrane were examined. These polymers were enabled to form small sized complexes of DNA, which interact with model membrane vesicles. It was found that the nature of hydrophobicity of the homopolymers significantly impacts rates of DNA complexation and the surface charge of the resulting complexes. These results highlight the prospect of the further examinations of these polymers as gene carriers.     

Examining the contributions of lipid shape and headgroup charge on bilayer behavior

To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na⁺ ions relative to the water/membrane interface. For one PA system, the Na⁺ ions are localized around the phosphate group. In the second PA system, the Na⁺ ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na⁺ ions. We find that the Na⁺ ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A_lipid and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na⁺ ions located near the ester groups than between the two PA bilayers.     

Brownian dynamics simulation of the lateral distribution of charged membrane components

Brownian dynamics simulations were performed to study the contribution of electric interactions between charged membrane components to their lateral distribution in a two-dimensional viscous liquid (bilayer lipid membrane). The electrostatic interaction potential was derived from an analytical solution of the linearized Poisson-Boltzmann equation for point charges in an electrolyte solution — membrane — electrolyte solution system. Equilibrium as well as dynamic quantities were investigated. The lateral organization of membrane particles, modelled by mobile cylinders in a homogeneous membrane separating two electrolyte solutions was described by spatial distribution functions, diffusion coefficients and cluster statistics. Disorder, local order and crystal-like arrangements were observed as a function of the particle charge, the closest possible distances between the charges and the particle density. The simulations revealed that the system is very sensitive to the position of the charges with respect to the electrolyte solution — membrane interface. Electrostatic interactions of charges placed directly on the membrane surface were almost negligible, whereas deeper charges demonstrated pronounced interaction. Biologically relevant parameters corresponded at most to local and transient ordering. It was found that lateral electric forces can give rise to a preferred formation of clusters with an even number of constituents provided that the closest possible charge-charge distances are small. It is concluded that lateral electrostatic interactions can account for local particle aggregations, but their impact on the global arrangement and movement of membrane components is limited. Correspondence to: D. Walther     

Membrane fusion by proline-rich Rz1 lipoprotein, the bacteriophage lambda Rz1 gene product.

The fusogenic properties of Rz1, the proline-rich lipoprotein that is the bacteriophage lambda Rz1 gene product, were studied. Light scattering was used to monitor Rz1-induced aggregation of artificial neutral (dipalmitoylphosphatidylcholine/cholesterol) and negatively charged (dipalmitoylphosphatidylcholine/cholesterol/dioleoylphosphatidylserin e) liposomes. Fluorescence assays [the resonance energy transfer between N-(7-nitro-2,1,3-benzoxadiazol-4-yl)phosphatidylethanolamine and N-(lissamine rhodamine B sulfonyl)dihexadecanol-sn-glycero-3-phosphoethanolamine lipid fluorescent probes, as well as fluorescent complex formation between terbium ions and dipicolinic acid encapsulated in two liposome populations and calcein fluorescence] were used to monitor Rz1-induced lipid mixing, contents mixing and leakage of neutral and negatively charged liposomes. The results demonstrated that Rz1 caused adhesion of neutral and negatively charged liposomes with concomitant lipid mixing; membrane distortion, leading to the fusion of liposomes and hence their internal content mixing; and local destruction of the membrane accompanied by leakage of the liposome contents. The use of artificial membranes showed that Rz1 induced the fusion of membranes devoid of any proteins. This might mean that the proline stretch of Rz1 allowed interaction with membrane lipids. It is suggested that Rz1-induced liposome fusion was mediated primarily by the generation of local perturbation in the bilayer lipid membrane and to a lesser extent by electrostatic forces.     

Synthesis and physicochemical and dynamic mechanical properties of a water-soluble chitosan derivative as a biomaterial

The physicochemical and rheological properties of a water-soluble chitosan (WSC) derivative were characterized in order to facilitate its use as a novel material for biomedical applications. The WSC was prepared by conjugating glycidyltrimethylammonium chloride (GTMAC) onto chitosan chains. Varying the molar ratio of GTMAC to chitosan from 3:1 to 6:1 produced WSCs with a degree of substitution (DS) that ranged from 56% to 74%. The WSC with the highest DS was soluble in water up to concentrations of 25 g/dL at room temperature. An increase in the polymer concentration gradually increased both the pH and conductivity of the WSC solutions. The rheological properties of the WSC solutions were found to be dependent on the salt and polymer concentrations as well as the DS value. In the absence of salt, the rheological behavior of the WSC was found to be typical of that for a polyelectrolyte in the dilute solution regime. However, the addition of salt decreased the viscosity of the polymer solution due to the reduction of electrostatic repulsions by the positively charged trimethylated ammonium groups of the WSC. In the concentrated regime, the viscosity of the WSCs was found to follow a power-law expression. The lowest DS WSC had the more favorable viscoelastic properties that were attributed to its high molecular weight, as confirmed by the stress relaxation spectra and intrinsic viscosity measurements. The effect of DS on the degree of interaction between WSC and the lipid egg phosphatidylcholine was investigated by FTIR analysis. Overall, the lower DS WSC had enhanced rheological properties and was capable of engaging in stronger intermolecular physical interactions.     

Interaction Between Equally Charged Membrane Surfaces Mediated by Positively and Negatively Charged Macro-Ions

In biological systems, charged membrane surfaces are surrounded by charged molecules such as electrolyte ions and proteins. Our recent experiments in the systems of giant phospholipid vesicles indicated that some of the blood plasma proteins (macro-ions) may promote adhesion between equally charged membrane surfaces. In this work, theory was put forward to describe an IgG antibody-mediated attractive interaction between negatively charged membrane surfaces which was observed in experiments on giant phospholipid vesicles with cardiolipin-containing membranes. The attractive interactions between negatively charged membrane surfaces in the presence of negatively and positively charged spherical macro-ions are explained using functional density theory and Monte Carlo simulations. Both, the rigorous solution of the variational problem within the functional density theory and the Monte Carlo simulations show that spatial and orientational ordering of macro-ions may give rise to an attractive interaction between negatively charged membrane surfaces. It is also shown that the distinctive spatial distribution of the charge within the macro-ions (proteins) is essential in this process.     

Grafting of polylysine with polyethylenoxide prevents demixing of O-pyromellitylgramicidin in lipid membranes

Both natural and synthetic polycations can induce demixing of negatively charged components in artificial and possibly in natural membranes. This process can result in formation of clusters (binding of several components to a polycation chain) and/or domains (aggregation of clusters and formation of a separate phase enriched in some particular component). In order to distinguish between these two phenomena, a model lipid membrane system containing ion channels, formed by a negatively charged peptide, O-pyromellitylgramicidin, and polycations of different structures was used. Microelectrophoresis of liposomes, changes in boundary potential of planar bilayers, the shape of compression curves and potentials of lipid and lipid/peptide monolayers were used to monitor the electrostatic factors in polymer adsorption to the membrane and peptide-polymer interactions. The synthesized PEO-grafted polylysine, PLL-PEO20000, did not induce peptide demixing monitored by stabilization of the gramicidin channels, in contrast to parent polylysine (PLL). Both polymers were shown to bind effectively to negatively charged liposomes and lipid monolayers, suggesting that the ineffectiveness of PLL-PEO20000 was not due to reduction of its binding. It was hypothesized that PLL-PEO20000 could not induce domain formation due to steric hindrance of long PEO chains preventing lateral fusion of clusters. Another copolymer, PLL-PEO4000, having four PEO chains of 4000 Da, exhibited intermediate effect between PLL and PLL-PEO20000, which shows the importance of the copolymer architecture for the effect on the lateral distribution of OPg channels. The model system can be relevant to regulation of lateral organization of ion channels and other components in natural membrane systems.     

Grafting of polylysine with polyethylenoxide prevents demixing of O-pyromellitylgramicidin in lipid membranes

Both natural and synthetic polycations can induce demixing of negatively charged components in artificial and possibly in natural membranes. This process can result in formation of clusters (binding of several components to a polycation chain) and/or domains (aggregation of clusters and formation of a separate phase enriched in some particular component). In order to distinguish between these two phenomena, a model lipid membrane system containing ion channels, formed by a negatively charged peptide, O-pyromellitylgramicidin, and polycations of different structures was used. Microelectrophoresis of liposomes, changes in boundary potential of planar bilayers, the shape of compression curves and potentials of lipid and lipid/peptide monolayers were used to monitor the electrostatic factors in polymer adsorption to the membrane and peptide-polymer interactions. The synthesized PEO-grafted polylysine, PLL-PEO20000, did not induce peptide demixing monitored by stabilization of the gramicidin channels, in contrast to parent polylysine (PLL). Both polymers were shown to bind effectively to negatively charged liposomes and lipid monolayers, suggesting that the ineffectiveness of PLL-PEO20000 was not due to reduction of its binding. It was hypothesized that PLL-PEO20000 could not induce domain formation due to steric hindrance of long PEO chains preventing lateral fusion of clusters. Another copolymer, PLL-PEO4000, having four PEO chains of 4000 Da, exhibited intermediate effect between PLL and PLL-PEO20000, which shows the importance of the copolymer architecture for the effect on the lateral distribution of OPg channels. The model system can be relevant to regulation of lateral organization of ion channels and other components in natural membrane systems.     

Molecular dynamics simulation of isothermal crystallisation of polymer chains around single polymer lamella

The isothermal crystallisation of polyethylene (PE) chains around single PE lamella in vacuum is investigated by molecular dynamic simulation. The crystallisation process is analysed in terms of the orientational order parameters, principal moments of inertia for the simulated systems. The effects of charge interactions between the polymer chains and lamella are discussed. It is found that the crystallisation process for uncharged systems can be divided into three stages: (1) adsorption, (2) orientation and (3) arrangement. The single polymer lamella changes a little during the three stages. PE chains are arranged parallel to the chain direction of the stems in the crystalline state. When considering the effect of charge interactions between the polymer chains and lamella, a different crystallisation process appears. The single polymer lamella is affected by the charged polymer chains.     

Structural analysis of telechelic polymer solution using dissipative particle dynamics simulations

A telechelic polymer is an amphiphilic polymer that can form micellar structures when dissolved in water. A telechelic polymer solution shows viscoelastic behaviour owing to the formation of characteristic networks, i.e. loops, bridges and dangling chains. For industrial purposes, telechelic polymers have many applications as thickening agents, such as in paints and cosmetics. Thus, it is desirable to predict and control the rheological properties of telechelic polymers. However, detailed studies at the molecular level have not yet been performed. In this study, I use the dissipative particle dynamics (DPD) method to investigate the relationship between the characteristic structural properties and the molecular structure in telechelic polymer solutions. I show that the morphology of telechelic polymer solutions depends on the concentration and chain length, the distribution of the end-to-end distance, the mean square end-to-end distance, the mean square radius of gyration and the time-averaged mean square displacement. Although an effect of entanglement is important for properties of polymer melts, the polymer chain composed of DPD particles cannot reproduce it. Therefore, I compare telechelic polymer solutions with and without the segmental repulsive potential (SRP), which can simulate the effect of entanglement in DPD simulations. The results indicate that it is necessary to include the SRP in DPD simulations to correctly analyse the behaviour of telechelic polymer solutions.     

Hydroxamic acid-containing hydrogels for nonabsorbed iron chelation therapy: synthesis, characterization, and biological evaluation

Iron overload is a severe clinical condition and can be largely prevented by the use of iron-specific chelating agents. A successful iron chelator needs to be orally active, nontoxic, and selective. In this study, hydrogels containing pendant hydroxamic acid groups have been synthesized as potential nonabsorbed chelators for iron in the gastrointestinal tract. The synthetic method employed to introduce hydroxamic acid groups to polymer chains involved reaction of polymer gels based on N-acryloxysuccinimide, acryloyl chloride, and (2-hydroxyethyl)acrylate monomers with hydroxylamine. These hydroxamic acid-functionalized polymer gels swell favorably in water and effectively sequester iron. In vitro iron-binding properties of these hydrogels were evaluated from their binding isotherms by use of iron(II) alone and in the presence of other competing metal ions. These polymers bind iron over a broad pH range. The iron-binding properties of the polymers were found to depend on the concentration of hydroxamate groups on polymer chains. The in vivo iron-binding efficacy of the polymers was evaluated in rat as the animal model. The polymers prevented an increase in serum hemoglobin and hematocrit levels in the animals, thus suggesting the prevention of systemic absorption of dietary iron from the gastrointestinal tract. The animals also maintained normal body weight during the treatment period, indicating the absence of any apparent toxicity associated with these polymers.     

Physicochemical characterization of carboxymethyl lipid A derivatives in relation to biological activity

Lipopolysaccharide (LPS) from the outer membrane of Gram-negative bacteria belongs to the most potent activators of the mammalian immune system. Its lipid moiety, lipid A, the 'endotoxic principle' of LPS, carries two negatively charged phosphate groups and six acyl chain residues in a defined asymmetric distribution (corresponding to synthetic compound 506). Tetraacyl lipid A (precursor IVa or synthetic 406), which lacks the two hydroxylated acyl chains, is agonistically completely inactive, but is a strong antagonist to bioactive LPS when administered to the cells before LPS addition. The two negative charges of lipid A, represented by the two phosphate groups, are essential for agonistic as well as for antagonistic activity and no highly active lipid A are known with negative charges other than phosphate groups. We hypothesized that the phosphate groups could be substituted by other negatively charged groups without changing the endotoxic properties of lipid A. To test this hypothesis, we synthesized carboxymethyl (CM) derivatives of hexaacyl lipid A (CM-506 and Bis-CM-506) and of tetraacyl lipid A (Bis-CM-406) and correlated their physicochemical with their endotoxic properties. We found that, similarly to compounds 506 and 406, also for their carboxymethyl derivatives a particular molecular ('endotoxic') conformation and with that, a particular aggregate structure is a prerequisite for high cytokine-inducing capacity and antagonistic activity, respectively. In other parameters such as acyl chain melting behaviour, antibody binding, activity in the Limulus lysate assay, and partially the binding of 3-deoxy-D-manno-oct-2-ulosonic acid transferase, strong deviations from the properties of the phosphorylated compounds were observed. These data allow a better understanding of endotoxic activity and its structural prerequisites.     

Theory for carrier-mediated zero-current conductance of bilayers extended to allow for nonequilibrium of interfacial reactions,spatially dependent mobilities and barrier shape

Summary A generalized form of the electrodiffusion equation, allowing for any shape of symmetrical energy barrier and any spatial dependence of the diffusion coefficient, is used to deduce theoretically the carrier-mediated conductance for thin (e.g., bilayer) membranes in the limit of low applied current. Both the Nernst-Planck and the Eyring single-barrier treatments are special cases of this more general approach, which allows for the effect of non-uniform properties of the lipid and non-uniform profiles of the forces acting within the membrane interior. Two independent mechanisms for ions to cross the membrane-solution interfaces are considered; namely, (1) the reaction at the interface between ions from solution and carriers from the membrane, and (2) the partition across the interfaces of complexes already formed in the solution. The rates of these reactions are taken into account using the rate equations of chemical kinetics; and the Poisson-Boltzmann equation is integrated in the aqueous solutions to evaluate the effect of charged polar head groups of the lipid. The analysis leads to an expression for the conductance, which, in the approximation of constant field, is an explicit function of such experimentally variable parameters as the concentrations and types of permeant ions and carriers in the aqueous phases, the total ionic strength and the nature of the polar head groups of the lipid. The functional relationship observable in an unknown membrane can, in principle, enable one to deduce such information as the mechanism of ion permeation across the interfaces, the magnitude of the surface charge, and the degree of ion-carrier complexation in the aqueous solutions.