Passive Internalization of Bioactive β-Casein Peptides into Phospholipid (POPC) Bilayers. Free Energy Landscapes from Unbiased Equilibrium MD Simulations at μs-Time Scale

Food Biophysics - Absorption of bioactive peptides in the intestinal epithelium take place in the apical or the basolateral tight junctions of the cells. Depending on the peptide size and...     

Membrane adsorption,folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1–2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug‐specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans‐membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α‐helical MPs as well as β‐barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge‐based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade. Proteins 2015; 83:1–24. © 2014 Wiley Periodicals, Inc.     

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

Conformational States of Melittin at a Bilayer Interface

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.Unstructured peptides fold into the membrane interface because partitioned hydrogen-bonded peptide bonds are energetically favorable compared to free peptide bonds (1–3). This folding process is central to the mechanisms of antimicrobial and cell-penetrating peptides, as well as to lipid interactions and stabilities of larger membrane proteins (4). The energetics of peptide partitioning into membrane interfaces can be described by a thermodynamic cycle (Fig. 1). State A is a theoretical state representing the fully unfolded peptide in water, B is the unfolded peptide in the membrane interface, C is the peptide in water, and D is the folded peptide in the membrane. The population of peptides in solution (State C) is best described as an ensemble of folded and unfolded conformations, whereas the population of peptides in State D generally is assumed to have a single, well-defined helicity, as shown in Fig. 1 A (5). Given that, in principle, folding in solution and in the membrane interface should follow the same basic rules, peptides in state D could reasonably be assumed to also be an ensemble. A fundamental question (5) is therefore whether peptides in state D can be correctly described as having a single helicity. Because differentiating an ensemble of conformations and a single conformation may be an impossible experimental task (5), molecular-dynamics (MD) simulations provide a unique high-resolution view of the phenomenon.Open in a separate windowFigure 1Thermodynamic cycles for peptide partitioning into a membrane interface. States A and B correspond to the fully unfolded peptide in solution and membrane interface, respectively. The folded peptide in solution is best described as an ensemble of unfolded and folded conformations (State C). State D is generally assumed to be one of peptides with a narrow range of conformations, but the state could actually be an ensemble of states as in the case of State C.Melittin is a 26-residue, amphipathic peptide that partitions strongly into membrane interfaces and therefore has become a model system for describing folding energetics (3,6–8). Here, we describe the structural dynamics of melittin in a dioleoylphosphatidylcholine (DOPC) bilayer by means of two extensive MD simulations using two different force fields.We extended a 12-ns equilibrated melittin-DOPC system (9) by 17 μs using the Anton specialized hardware (10) with the CHARMM22/36 protein/lipid force field and CMAP correction (11,12) (see Fig. S1 and Fig. S2 in the Supporting Material). To explore force-field effects, a similar system was simulated for 2 μs using the OPLS force field (13) (see Methods in the Supporting Material). In agreement with x-ray diffraction measurements on melittin in DOPC multilayers (14), melittin partitioned spontaneously into the lipid headgroups at a position below the phosphate groups at similar depth as glycerol/carbonyl groups (Fig. 2).Open in a separate windowFigure 2Melittin partitioned into the polar headgroup region of the lipid bilayer. (A) Snapshot of the simulation cell showing two melittin molecules (MLT1 and MLT2, in yellow) at the lipid-water interface. (B) Density cross-section of the simulation cell extracted from the 17-μs simulation. The peptides are typically located below the lipid phosphate (PO4) groups, in a similar depth as the glycerol/carbonyl (G/C) groups.To describe the secondary structure for each residue, we defined helicity by backbone dihedral angles (φ, ψ) within 30° from the ideal α-helical values (–57°, –47°). The per-residue helicity in the CHARMM simulation displays excellent agreement with amide exchange rates from NMR measurements that show a proline residue to separate two helical segments, which are unfolded below Ala⁵ and above Arg²² (15) (Fig. 3 A). In contrast, the OPLS simulation failed to reproduce the per-residue helicity except for a short central segment (see Fig. S3).Open in a separate windowFigure 3Helicity and conformational distribution of melittin as determined via MD simulation. (A) Helicity per residue for MLT1 and MLT2. (B) Corresponding evolution of the helicity. (C) Conformational distributions over the entire 17-μs simulation.Circular dichroism experiments typically report an average helicity of ∼70% for melittin at membrane interfaces (3,6,16,17), but other methods yield average helicities as high as 85% (15,18). Our CHARMM simulations are generally consistent with the experimental results, especially amide-exchange measurements (15); melittin helicity averaged to 78% for MLT1, whereas MLT2 transitioned from 75% to 89% helicity at t ≈ 8 μs, with an overall average helicity of 82% (Fig. 3 B). However, in the OPLS simulation, melittin steadily unfolds over the first 1.3 μs, after which the peptide remains only partly folded, with an average helicity of 33% (see Fig. S3). Similar force-field-related differences in peptide helicity were recently reported, albeit at shorter timescales (19). Although suitable NMR data are not presently available, we have computed NMR quadrupolar splittings for future reference (see Fig. S4).To answer the question asked in this article—whether the conformational space of folded melittin in the membrane interface can be described by a narrow distribution—the helicity distributions for the equilibrated trajectories are shown in Fig. 3 C. Whereas MLT1 in the CHARMM simulation produces a single, narrow distribution of the helicity, MLT2 has a bimodal distribution as a consequence of the folding event at t ≈ 8 μs (Fig. 3 C). We note that CHARMM force fields have a propensity for helix-formation and this transition might therefore be an artifact. We performed a cluster analysis to describe the structure of the peptide in the membrane interface. The four most populated conformations in the CHARMM simulation are shown in Fig. 4. The dominant conformation for both peptides was a helix kinked at G12 and unfolded at the last 5–6 residues of the C-terminus. The folding transition of MLT2 into a complete helix is visible by the 48% occupancy of a fully folded helix.Open in a separate windowFigure 4Conformational clusters of the two melittin peptides (MLT1 and MLT2) from the 17-μs CHARMM simulation in DOPC. Clustering is based on Cα-RMSD with a cutoff criterion of 2 Å.We conclude that the general assumption when calculating folding energetics holds: Folded melittin partitioned into membrane interfaces can be described by a narrow distribution of conformations. Furthermore, extended (several microsecond) simulations are needed to differentiate force-field effects. Although the CHARMM and OPLS simulations would seem to agree for the first few hundred nanoseconds, the structural conclusions differ drastically with longer trajectories, with CHARMM parameters being more consistent with experiments. However, as implied by the difference in substate distributions between MLT1 and MLT2, 17 μs might not be sufficient to observe the fully equilibrated partitioning process. The abrupt change in MLT2 might indicate that the helicity will increase to greater than experimentally observed in a sufficiently long simulation. On the other hand, it could be nothing more than a transient fluctuation. Increased sampling will provide further indicators of convergence of the helix partitioning process.     

Monte Carlo folding of trans-membrane helical peptides in an implicit generalized Born membrane

An efficient Monte Carlo (MC) algorithm using concerted backbone rotations is combined with a recently developed implicit membrane model to simulate the folding of the hydrophobic transmembrane domain M2TM of the M2 protein from influenza A virus and Sarcolipin at atomic resolution. The implicit membrane environment is based on generalized Born theory and has been calibrated against experimental data. The MC sampling has previously been used to fold several small polypeptides and been shown to be equivalent to molecular dynamics (MD). In combination with a replica exchange algorithm, M2TM is found to form continuous membrane spanning helical conformations for low temperature replicas. Sarcolipin is only partially helical, in agreement with the experimental NMR structures in lipid bilayers and detergent micelles. Higher temperature replicas exhibit a rapidly decreasing helicity, in agreement with expected thermodynamic behavior. To exclude the possibility of an erroneous helical bias in the simulations, the model is tested by sampling a synthetic Alanine-rich polypeptide of known helicity. The results demonstrate there is no overstabilization of helical conformations, indicating that the implicit model captures the essential components of the native membrane environment for M2TM and Sarcolipin.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

How reliable are molecular dynamics simulations of membrane active antimicrobial peptides?

Membrane-active antimicrobial peptides (AMPs) are challenging to study experimentally, but relatively easy to investigate using molecular dynamics (MD) computer simulations. For this reason, a large number of MD studies of AMPs have been reported over recent years. Yet relatively little effort has focused on the validity of such simulations. Are these results reliable, and do they agree with what is known experimentally? And how much meaningful information can be obtained? To answer these questions, we demonstrate here some of the requirements and limitations of running MD simulations for several common AMPs: PGLa, melittin, maculatin and BP100. The two most important findings are: (a) simulation results depend strongly on force field parameters, making experimental verification of the simulations obligatory, and (b) slow orientational and conformational fluctuations mean that much longer sampling timescales (multi-μs) are needed if quantitative agreement between simulation averages and experimental data is to be achieved. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.