Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 × 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the ‘quality’ of the starting structure seems mainly to influence the simulation stability of the transmembrane β-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal α-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a ∼10 ns timescale.     

Understanding peptide interactions with the lipid bilayer: a guide to membrane protein engineering

In recent years, studies of the interactions of designed hydrophobic and amphipathic polypeptides with biological membranes have progressed considerably. In-plane and transmembrane helical domains have been engineered as well as sequences that exhibit dynamic distributions of different topologies. These sequences not only help our understanding of the thermodynamic interaction contributions that determine peptide orientation, but also exhibit interesting biological functions as autonomous units. In addition, helices are considered to be folding intermediates of multi-spanning membrane proteins. The specificity and thermodynamic stability of transmembrane helix-helix interactions or the assembly of beta sheets at the membrane surface have, therefore, become a focus of ongoing research.     

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells’ membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane’s mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.     

Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

Molecular dynamics simulations of microstructure and mixing dynamics of cryoprotective solvents in water and in the presence of a lipid membrane

Molecular dynamics (MD) simulation is used to investigate the solubility behavior of cryoprotective (CP) solvents, such as DMSO, ethylene glycol (EG) and glycerol (GL), in pure water and in the presence of a lipid membrane. The MD study is focused on an equilibration timescale required for mixing large CP aggregates with aqueous and aqueous/lipid environments. The MD analysis demonstrates that DMSO mixes rapidly with water, so that all solute molecules are uniformly distributed in the equilibrium aqueous solution. Our investigation of the microstructure of binary EG/water and GL/water systems reveals that, despite the miscibility of both CP solvents with water, they are not ideally mixed in aqueous solutions at the molecular level. The MD simulations show that the mixing dynamics of the large CP cluster and surrounding water is found to be strongly dependent on nature of hydrophilic and hydrophobic interactions acting between cryoprotectant molecules. In particular, a spatial hydrogen-bond network formed between CP molecules plays an important role in the mixing dynamics between CP agents and water. A further analysis on the mixing behavior of the CP solvents with pure water and with aqueous solutions at a lipid membrane interface shows that, due to strong binding of the CP molecules to membrane surface, the equilibration process in the lipid environment becomes very slow, at least of the order of microseconds. The MD results are discussed in the context of the better understanding on the composition of the aqueous mixtures of the EG and GL solvents. Knowledge of the microstructure and the dynamics of these systems helps to develop better cryopreservation protocols and to propose more optimal cooling/warming regimes for cellular cryosolutions.     

Cellular uptake of transportan 10 and its analogs in live cells: Selectivity and structure-activity relationship studies

Transportan 10 (TP10) is an amphipathic cell-penetrating peptide with high translocation ability. In order to obtain more details of structure-activity relationship of TP10, we evaluated the effects of structure and charge on its translocation ability. Our results demonstrated that disrupting the helical structure or Arg substitution could remarkably decrease the cellular uptake of TP10. However, increasing the number of positive charge was an effective strategy to enhance translocation ability of TP10. Furthermore, the molecular dynamics simulation supported the results derived from experiments, suggesting that higher membrane disturbance leads to higher cellular uptake of peptides. In addition, our study also demonstrated TP10 and its analogs preferentially entered cancer cells rather than normal cells. The uptake selectivity toward cancer cells makes TP10 and its analogs as potent CPPs for drug delivery.     

Microscopic calculations of local lipid membrane permittivities and diffusion coefficients for application to electroporation analyses

Interaction of electric fields with biological systems has begun to receive considerable attention for applications that include field-assisted drug delivery, medical interventions, and genetic engineering. External fields induce the strongest effects at membranes with electroporation being a common feature. Membrane transport in this context of poration is often based on continuum approaches utilizing macroscopic parameters such as the permittivity, diffusion coefficients, and mobilities. In such modeling, field dependences, local inhomogeneities, and microscopic details are usually ignored. Here, a molecular dynamics (MD) scheme is used for a more rigorous and physically realistic evaluation of such parameters for potential application to electroporative transport model development. A suitable membrane structure containing a nanopore derived from MD analysis is used as the initial geometric configuration. Both static and frequency dependent diffusion coefficients have been evaluated. Permittivities are also calculated and shown to be dramatically non-uniform in the vicinity of membranes under high external fields. A positive feedback mechanism leading to enhanced membrane fields is discussed.     

Effects of deformability and thermal motion of lipid membrane on electroporation: by molecular dynamics simulations

Effects of mechanical properties and thermal motion of POPE lipid membrane on electroporation were studied by molecular dynamics simulations. Among simulations in which specific atoms of lipids were artificially constrained at their equilibrium positions using a spring with force constant of 2.0 kcal/(mol Ų) in the external electric field of 1.4 kcal/(mol Å e), only constraint on lateral motions of lipid tails prohibited electroporation while non-tail parts had little effects. When force constant decreased to 0.2 kcal/(mol Ų) in the position constraints on lipid tails in the external electric field of 2.0 kcal/(mol Å e), water molecules began to enter the membrane. Position constraints of lipid tails allow water to penetrate from both sides of membrane. Thermal motion of lipids can induce initial defects in the hydrophobic core of membrane, which are favorable nucleation sites for electroporation. Simulations at different temperatures revealed that as the temperature increases, the time taken to the initial pore formation will decrease.     

Direct AFM observation of saposin C-induced membrane domains in lipid bilayers: from simple to complex lipid mixtures

Saposin C (Sap C) is a small glycoprotein required by glucosylceramidase (GCase) for hydrolysis of glucosylceramide to ceramide and glucose in lysosomes. The molecular mechanism underlying Sap C stimulation of the enzyme activation is not fully understood. Here, atomic force microscopy (AFM) has been used to study Sap C-membrane interactions under physiological conditions. First, to establish how Sap C-membrane interactions affect membrane structure, lipid bilayers containing zwitterionic and anionic phospholipids were used. It was observed that Sap C induced two types of membrane restructuring effects, i.e., the formation of patch-like domains and membrane destabilization. Bilayers underwent extensive structural reorganization. To validate the biological importance of the membrane restructuring effects, interaction of Sap C with lipid bilayers composed of cholesterol, sphingomyelin, and zwitterionic and anionic phospholipids were studied. Although similar membrane restructuring effects were observed, Sap C-membrane interactions, in this case, were remarkably modulated and their effects were restricted to a limited area. As a result, nanometer-sized domains were formed. The establishment of a model membrane system will allow us to further study the dynamics, structure and mechanism of the Sap C-associated membrane domains and to examine the important role that these domains may play in enzyme activation.     

Interaction of S413-PV cell penetrating peptide with model membranes: relevance to peptide translocation across biological membranes.

Cell penetrating peptides (CPPs) have been successfully used to mediate the intracellular delivery of a wide variety of molecules of pharmacological interest both in vitro and in vivo, although the mechanisms by which the cellular uptake occurs remain unclear and controversial. Following our previous work demonstrating that the cellular uptake of the S4(13)-PV CPP occurs mainly through an endocytosis-independent mechanism, we performed a detailed biophysical characterization of the interaction of this peptide with model membranes. We demonstrate that the interactions of the S4(13)-PV peptide with membranes are essentially of electrostatic nature. As a consequence of its interaction with negatively charged model membranes, the S4(13)-PV peptide becomes buried into the lipid bilayer, which occurs concomitantly with significant peptide conformational changes that are consistent with the formation of a helical structure. Comparative studies using two related peptides demonstrate that the conformational changes and the extent of cell penetration are dependent on the peptide sequence, indicating that the helical structure acquired by the S4(13)-PV peptide is relevant for its nonendocytic uptake. Overall, our data suggest that the cellular uptake of the S4(13)-PV CPP is a consequence of its direct translocation through cell membranes, following conformational changes induced by peptide-membrane interactions.     

Improved protein free energy calculation by more accurate treatment of nonbonded energy: Application to chymotrypsin inhibitor 2, V57A

We developed a software package for improved free energy calculation, in which spherical solvent boundary potential, cell multipole method, and Nosé-Hoover equation are employed. The performance of the developed software package is demonstrated in the case of valine to alanine mutation of the 57th residue in chymotrypsin inhibitor 2. By using this package, we obtained results quantitatively comparable to experimental results. By the free energy component analysis, it is shown that leucine 51, arginine 65, arginine 67, and phenylalanine 69 residues contribute significantly to the total free energy shift, ΔΔG. Among them, contribution from the hydrophilic arginine 67 residue, which is in close contact with the mutation site, is the largest. Structure around the mutation site is largely changed by the mutation. The structure change is caused mainly by two effects, hydrophobic interaction and short-range interaction along the sequence. Effects of Nosé-Hoover algorithm and Kirkwood reaction field are also discussed. Proteins 30:388–400, 1998. © 1998 Wiley-Liss, Inc.