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.     

Molecular design of modified polyacrylamide for the salt tolerance

In our work, three kinds of functional monomers were selected to modify polyacrylamide (PAM) or partially hydrolyzed polyacrylamide (HPAM) by molecular dynamics simulation so as to achieve the stronger salt-tolerance of modified HM-HPAM. The radius of gyration (R (g)), the hydrodynamic radius (R (H)), the effective length (L (ef)) and the intrinsic viscosity ([η]) for modified PAM or HPAM were studied in aqueous solutions with different ionic strength at 298?K. The results showed that modified HM-HPAM has a stronger salt tolerance and the salt tolerance increases gradually from HM-HPAM1 to HM-HPAM3 because the monomers with different steric hindrance would reduce the curliness of molecular chains and, consequently, improve the salt tolerance. So, introducing the steric hindrance monomer into polymer will increase the salt tolerance of the polymer and it is indicated that the simulated results agree with the experimental results very well. Furthermore, the radial distribution function (RDF) has been used to investigate the effect of NaCl on the hydration of the -COO- groups of the HM-HPAM from microscopic view.     

Mobility of adsorbed proteins: a Brownian dynamics study

We simulate the adsorption of lysozyme on a solid surface, using Brownian dynamics simulations. A protein molecule is represented as a uniformly charged sphere and interacts with other molecules through screened Coulombic and double-layer forces. The simulation starts from an empty surface and attempts are made to introduce additional proteins at a fixed time interval that is inversely proportional to the bulk protein concentration. We examine the effect of ionic strength and bulk protein concentration on the adsorption kinetics over a range of surface coverages. The structure of the adsorbed layer is examined through snapshots of the configurations and quantitatively with the radial distribution function. We extract the surface diffusion coefficient from the mean square displacement. At high ionic strengths the Coulombic interaction is effectively shielded, leading to increased surface coverage. This effect is quantified with an effective particle radius. Clustering of the adsorbed molecules is promoted by high ionic strength and low bulk concentrations. We find that lateral protein mobility decreases with increasing surface coverage. The observed trends are consistent with previous theoretical and experimental studies.     

A Robust Water Potential Parameterisation

Abstract

We compare molecular dynamics simulation results for the properties of liquid water predicted by four novel water potential models. These models are designed as a combination of parameters taken from the dedicated but brittle TIP3P water potential, and the more flexible but less accurate parameterisations such as the Dreiding and Universal force fields. We find that a hybrid of Dreiding and TIP3P delivers the best results, yielding a density, diffusion coefficient and radial distribution function in good agreement with experiment, performing in some respects even better than the dedicated reference TIP3P model. Another Dreiding based force field predicts semi-quantitative results for the water structure and dynamics while the Universal force field based models are incapable of simulating a condensed phase of water at all, continuing to expand indefinitely. These observations are useful for selecting and designing robust water force field parameterisations that can be used for general simulation purposes.     

Effect of Temperature and Capsid Tail on the Packing and Ejection of Viral DNA

We use a simulation technique based on molecular dynamics and stochastic rotation model to present the effect of temperature and capsid tail on the packaging and ejection processes of semiflexible polymers. We consider two types of solvents, a good solvent, where the polymer is neutral and repulsion interactions among its various sections are favored, and one where the polymer is charged, giving rise to extra electrostatic reaction. For tailless capsids, we find that packing a neutral polymer is slightly slower at higher temperatures whereas its ejection is slightly slower at lower temperatures. We find the same trend for a charged polymer but the effect is noticeably larger. At a high enough temperature, we notice that packing a charged polymer can be stopped. On the other hand, at fixed temperature and regardless whether the polymer is charged, packing is much easier for a capsid with a tail whereas ejection is much slower. The effect of including the tail on the dynamics of a charged polymer, in particular, is rather significant: more packing fraction is facilitated at higher temperatures due to more ordered polymer configuration inside the capsid. In contrast, during ejection the tail traps the last remaining beads for quite some time before allowing full ejection. We interpret these results in terms of entropic and electrostatic forces.     

Mechanism of fiber assembly: treatment of Aβ peptide aggregation with a coarse-grained united-residue force field

The growth mechanism of β-amyloid (Aβ) peptide fibrils was studied by a physics-based coarse-grained united-residue model and molecular dynamics (MD) simulations. To identify the mechanism of monomer addition to an Aβ_1-40 fibril, we placed an unstructured monomer at a distance of 20 Å from a fibril template and allowed it to interact freely with the latter. The monomer was not biased towards fibril conformation by either the force field or the MD algorithm. With the use of a coarse-grained model with replica-exchange molecular dynamics, a longer timescale was accessible, making it possible to observe how the monomers probe different binding modes during their search for the fibril conformation. Although different assembly pathways were seen, they all follow a dock-lock mechanism with two distinct locking stages, consistent with experimental data on fibril elongation. Whereas these experiments have not been able to characterize the conformations populating the different stages, we have been able to describe these different stages explicitly by following free monomers as they dock onto a fibril template and to adopt the fibril conformation (i.e., we describe fibril elongation step by step at the molecular level). During the first stage of the assembly (“docking”), the monomer tries different conformations. After docking, the monomer is locked into the fibril through two different locking stages. In the first stage, the monomer forms hydrogen bonds with the fibril template along one of the strands in a two-stranded β-hairpin; in the second stage, hydrogen bonds are formed along the second strand, locking the monomer into the fibril structure. The data reveal a free-energy barrier separating the two locking stages. The importance of hydrophobic interactions and hydrogen bonds in the stability of the Aβ fibril structure was examined by carrying out additional canonical MD simulations of oligomers with different numbers of chains (4-16 chains), with the fibril structure as the initial conformation. The data confirm that the structures are stabilized largely by hydrophobic interactions and show that intermolecular hydrogen bonds are highly stable and contribute to the stability of the oligomers as well.     

The importance of coulombic end effects: experimental characterization of the effects of oligonucleotide flanking charges on the strength and salt dependence of oligocation (L8+) binding to single-stranded DNA oligomers.

Binding constants Kobs, expressed per site and evaluated in the limit of zero binding density, are quantified as functions of salt (sodium acetate) concentration for the interactions of the oligopeptide ligand KWK6NH2 (designated L8+, with ZL = 8 charges) with three single-stranded DNA oligomers (ss dT-mers, with |ZD| = 15, 39, and 69 charges). These results provide the first systematic experimental information about the effect of changing |ZD| on the strength and salt dependence of oligocation-oligonucleotide binding interactions. In a comparative study of L8+ binding to poly dT and to a short dT oligomer (|ZD| = 10),. Proc. Natl. Acad. Sci. USA. 93:2511-2516) demonstrated the profound thermodynamic effects of phosphate charges that flank isolated nonspecific L8+ binding sites on DNA. Here we find that both Kobs and the magnitude of its power dependence on salt activity (|SaKobs|) increase monotonically with increasing |ZD|. The dependences of Kobs and SaKobs on |ZD| are interpreted by introducing a simple two-state thermodynamic model for Coulombic end effects, which accounts for our finding that when L8+ binds to sufficiently long dT-mers, both DeltaGobso = -RT ln Kobs and SaKobs approach the values characteristic of binding to poly-dT as linear functions of the reciprocal of the number of potential oligocation binding sites on the DNA lattice. Analysis of our L8+-dT-mer binding data in terms of this model indicates that the axial range of the Coulombic end effect for ss DNA extends over approximately 10 phosphate charges. We conclude that Coulombic interactions cause an oligocation (with ZL < |ZD|) to bind preferentially to interior rather than terminal binding sites on oligoanionic or polyanionic DNA, and we quantify the strong increase of this preference with decreasing salt concentration. Coulombic end effects must be considered when oligonucleotides are used as models for polyanionic DNA in thermodynamic studies of the binding of charged ligands, including proteins.     

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     

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl₂, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5–12) carrying three positive charges on lysine⁺ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.     

Particle Transport through Hydrogels Is Charge Asymmetric

Transport processes within biological polymer networks, including mucus and the extracellular matrix, play an important role in the human body, where they serve as a filter for the exchange of molecules and nanoparticles. Such polymer networks are complex and heterogeneous hydrogel environments that regulate diffusive processes through finely tuned particle-network interactions. In this work, we present experimental and theoretical studies to examine the role of electrostatics on the basic mechanisms governing the diffusion of charged probe molecules inside model polymer networks. Translational diffusion coefficients are determined by fluorescence correlation spectroscopy measurements for probe molecules in uncharged as well as cationic and anionic polymer solutions. We show that particle transport in the charged hydrogels is highly asymmetric, with diffusion slowed down much more by electrostatic attraction than by repulsion, and that the filtering capability of the gel is sensitive to the solution ionic strength. Brownian dynamics simulations of a simple model are used to examine key parameters, including interaction strength and interaction range within the model networks. Simulations, which are in quantitative agreement with our experiments, reveal the charge asymmetry to be due to the sticking of particles at the vertices of the oppositely charged polymer networks.     

Structural organization of a chain molecule with specific charge distribution: A molecular dynamics study

We report the results of molecular dynamics simulations of a charged bead-monomer chain molecule with charge distribution adopted from immunoglobulin-binding domain B1 of protein-g. The beads of the model are connected by invariable bonds and interact with each other via the Coulomb potential. To study the low-temperature conformational space of the designed model we use standard canonical, microcanonical and multicanonical molecular dynamics simulations. We find that at low temperature T = T c the chain undergoes a continuous freezing transition into one of many low-energy conformations. Below T c the molecule is a compact globule composed of an inner core, containing mostly charged monomers, and an outer corona, consisting of all the rest neutral units. All frozen conformations have almost equal potential energy but differ in structure. The potential energy surface of the model does not posses a pronounced ground-state minimum--an essential feature of protein-like heteropolymers.     

Quasi-elastic laser light scattering from solutions and gels of hemoglobin S.

Quasi-elastic light scattering has been used to examine solutions and gels of deoxyhemoglobin S. The autocorrelation function is found to decay with a characteristic exponential relaxation which can be ascribed to the diffusion of monomer (64,000 molecular weight) hemoglobin S molecules. In the absence of polymers, the relaxation time is in good agreement with previous measurements of the diffusion coefficient for solutions of normal human hemoglobin. In the presence of the polymer phase, a large (greater than 200-fold) increase in the scattered intensity is observed but no contribution to the decay of the autocorrelation function from the motion of the aligned polymer phase can be detected. Heterodyning between the time-independent scattering amplitude from the polymers and the time-dependent scattering of the diffusing monomers results in a twofold increase in the relaxation time arising from monomer diffusion.     

Calcium Enhances Binding of Aβ Monomer to DMPC Lipid Bilayer

Using isobaric-isothermal replica-exchange molecular dynamics and the all-atom explicit-solvent model, we studied the equilibrium binding of Aβ monomers to a zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer coincubated with calcium ions. Using our previous replica-exchange molecular dynamics calcium-free simulations as a control, we reached three conclusions. First, calcium ions change the tertiary structure of the bound Aβ monomer by destabilizing several long-range intrapeptide interactions, particularly the salt bridge Asp²³-Lys²⁸. Second, calcium strengthens Aβ peptide binding to the DMPC bilayer by enhancing electrostatic interactions between charged amino acids and lipid polar headgroups. As a result, Aβ monomer penetrates deeper into the bilayer, making disorder in proximal lipids and bilayer thinning more pronounced. Third, because calcium ions demonstrate strong affinity to negatively charged amino acids, a considerable influx of calcium into the area proximal to the bound Aβ monomer is observed. Consequently, the localizations of negatively charged amino acids and calcium ions in the Aβ binding footprint overlap. Based on our data, we propose a mechanism by which calcium ions strengthen Aβ-bilayer interactions. This mechanism involves two factors: 1) calcium ions make the DMPC bilayer partially cationic and thus attractive to the anionic Aβ peptide; and 2) destabilization of the Asp²³-Lys²⁸ salt bridge makes Lys²⁸ available for interactions with the bilayer. Finally, we conclude that a single Aβ monomer does not promote permeation of calcium ions through the zwitterionic bilayer.     

Evaluation of the energetic contribution of interhelical Coulombic interactions for coiled coil helix orientation specificity.

Coiled coils are formed by two or more alpha-helices that align in a parallel or an antiparallel relative orientation. The factors that determine a preference for a given relative helix orientation are incompletely understood. The helix orientation preference for the designed coiled coil, Acid-a1-Base-a1, was measured previously. This model system therefore provides a means for the experimental determination of the energetic contribution of a variety of interactions to helix orientation specificity.The antiparallel preference for Acid-a1-Base-a1 is imparted by a single buried polar interaction. Interhelical Coulombic interactions between residues at the e and g positions have been proposed to influence helix orientation preference. In the Acid-a1-Base-a1 heterodimer, potentially attractive Coulombic interactions are expected in both orientations. To determine the energetic consequences of Coulombic interactions for helix orientation preference, we have positioned a single charged residue in each peptide such that exclusively favorable interhelical Coulombic interactions can occur only in the parallel orientation. In contrast, two potentially repulsive interactions are expected in the antiparallel orientation. Because the buried polar interaction can occur only in the antiparallel orientation, interhelical Coulombic interactions favor the parallel orientation and the potential to form a buried polar interaction favors the antiparallel orientation. We find no clear preference for an antiparallel orientation in the resulting heterodimer, Acid-Ke-Base-Eg, suggesting that interhelical Coulombic interactions and a buried polar interaction are of approximately equal importance for helix orientation specificity. Stability measurements indicate that maintenance of all favorable electrostatic interactions and/or avoidance of two potentially repulsive interactions contributes approximately 2.1 kcal/mol to helix orientation preference.     

αB-crystallin polydispersity is a consequence of unbiased quaternary dynamics

The inherent heterogeneity of many protein assemblies complicates characterization of their structure and dynamics, as most biophysical techniques require homogeneous preparations of isolated components. For this reason, quantitative studies of the molecular chaperone αB-crystallin, which populates a range of interconverting oligomeric states, have been difficult, and the physicochemical basis for its polydispersity has remained unknown. Here, we perform mass spectrometry experiments to study αB-crystallin and extract detailed information as to its oligomeric distribution and exchange of subunits under a range of conditions. This allows a determination of the thermodynamic and kinetic parameters that govern the polydisperse ensemble and enables the construction of a simple energy profile for oligomerization. We find that the quaternary structure and dynamics of the protein can be explained using a simple model with just two oligomer-independent interactions (i.e., interactions that are energetically identical in all oligomers from 10mers to 40mers) between constituent monomers. As such, the distribution of oligomers is governed purely by the dynamics of individual monomers. This provides a new means for understanding the polydispersity of αB-crystallin and a framework for interrogating other heterogeneous protein assemblies.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

Molecular dynamic simulation study of cholesterol and conjugated double bonds in lipid bilayers

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is believed to have anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be associated with anti-obesity effects. In this paper we extend earlier molecular dynamics (MD) simulations of pure CLA–phosphatidylcholine bilayers to investigate the comparative effects of cholesterol on bilayers composed of the two respective isomers. Simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed in which, for each isomer, the simulated bilayers contained 10% and 30% cholesterol (Chol). From MD trajectories we calculate and compare structural properties of the bilayers, including areas per molecule, thickness of bilayers, tilt angle of cholesterols, order parameter profiles, and one and two-dimensional radial distribution function (RDF), as functions of Chol concentration. While the structural effect of cholesterol is approximately the same for both isomers, we find differences at an atomistic level in order parameter profiles and in two-dimensional radial distribution functions.