Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

New developments in first-principles excited-state dynamics simulations: unveiling the solvent specificity of excited anionic cluster relaxation and electron solvation

Charge-transfer-to-solvent excited iodide–polar solvent molecule clusters, [I^? (Solv)_n]^*, have attracted substantial interest over the past 20 years as they can undergo intriguing relaxation processes leading ultimately to the formation of gas-phase molecular analogues of the solvated electron. In this review article, we present a comprehensive overview of the development and application of state-of-the-art first-principles molecular dynamics simulation approaches to understand and interpret the results of femtosecond photoelectron spectroscopy experiments on [I^? (Solv)_n]^* relaxation, which point to a high degree of solvent specificity in the electron solvation dynamics. The intricate molecular details of the [I^? (Solv)_n]^* relaxation process are presented, and by contrasting the relaxation mechanisms of clusters with several different solvents (water, methanol and acetonitrile), the molecular basis of the solvent specificity of electron solvation in [I^? (Solv)_n]^* is uncovered, leading to a more refined view of the manifestation of electron solvation in small gas-phase clusters.     

Free energy surfaces of miniproteins with a betabetaalpha motif: replica exchange molecular dynamics simulation with an implicit solvation model

Designed miniproteins with a betabetaalpha motif, such as BBA5, 1FSD, and 1PSV can serve as a benchmark set to test the validity of all-atom force fields with computer simulation, because they contain all the basic structural elements in protein folding. Unfortunately, it was found that the standard all-atom force fields with the generalized Born (GB) implicit solvation model tend to produce distorted free energy surfaces for the betabetaalpha proteins, not only because energetically those proteins need to be described by more balanced weights of the alpha- and beta-strands, but also because the GB implicit solvation model suffers from overestimated salt bridge effects. In an attempt to resolve these problems, we have modified one of the standard all-atom force fields in conjunction with the GB model, such that each native state of the betabetaalpha proteins is in its free energy minimum state with reasonable energy barriers separating local minima. With this modified energy model, the free energy contour map in each protein was constructed from the replica exchange molecular dynamics REMD simulation. The resulting free energy surfaces are significantly improved in comparison with previous simulation results and consistent with general views on small protein folding behaviors with realistic topology and energetics of all three proteins.     

Membrane partitioning of peptide aggregates: coarse-grained molecular dynamics simulations

Abstract

Coarse-grained molecular dynamics (CGMD) simulation technique (MARTINI force field) is applied to monitor the aggregation of helical peptides representing the transmembrane sequence and its extension of bone marrow stromal cell antigen 2 (BST-2). One of the peptides is coupled with a protein transducing domain (PTD) of nine arginine residues (R9) at its N-terminal side as well as a peptide, pep11**, which has been shown to bind to human papilloma virus 16 (HPV16) E6 oncoprotein. A short hydrophobic stretch of the transmembrane domain (TMD) of BST-2 aggregates the fastest and inserts into a lipid membrane. An aggregate of R9-pep11** attaches to the membrane via simultaneous contact of many arginine residues. Monomers from the aggregates of the shortest of the hydrophobic TMDs dissolve into the opposing leaflet when the aggregate spans the bilayer. A ‘flipping’ of the individual monomeric peptides is not observed.

Communicated by Ramaswamy H. Sarma     

Implicit solvent model studies of the interactions of the influenza hemagglutinin fusion peptide with lipid bilayers

The "fusion peptide," a segment of approximately 20 residues of the influenza hemagglutinin (HA), is necessary and sufficient for HA-induced membrane fusion. We used mean-field calculations of the free energy of peptide-membrane association (DeltaG(tot)) to deduce the most probable orientation of the fusion peptide in the membrane. The main contributions to DeltaG(tot) are probably from the electrostatic (DeltaG(el)) and nonpolar (DeltaG(np)) components of the solvation free energy; these were calculated using continuum solvent models. The peptide was described in atomic detail and was modeled as an alpha-helix based on spectroscopic data. The membrane's hydrocarbon region was described as a structureless slab of nonpolar medium embedded in water. All the helix-membrane configurations, which were lower in DeltaG(tot) than the isolated helix in the aqueous phase, were in the same (wide) basin in configurational space. In each, the helix was horizontally adsorbed at the water-bilayer interface with its principal axis parallel to the membrane plane, its hydrophobic face dissolved in the bilayer, and its polar face in the water. The associated DeltaG(tot) value was approximately -8 to -10 kcal/mol (depending on the rotameric state of one of the phenylalanine residues). In contrast, the DeltaG(tot) values associated with experimentally observed oblique orientations were found to be near zero, suggesting they are marginally stable at best. The theoretical model did not take into account the interactions of the polar headgroups with the peptide and peptide-induced membrane deformation effects. Either or both may overcompensate for the DeltaG(tot) difference between the horizontal and oblique orientations.     

Antimicrobial peptide RP-1 structure and interactions with anionic versus zwitterionic micelles

Topologically, platelet factor-4 kinocidins consist of distinct N-terminal extended, C-terminal helical, and interposing gamma-core structural domains. The C-terminal alpha-helices autonomously confer direct microbicidal activity, and the synthetic antimicrobial peptide RP-1 is modeled upon these domains. In this study, the structure of RP-1 was assessed using several complementary techniques. The high-resolution structure of RP-1 was determined by NMR in anionic sodium dodecyl sulfate (SDS) and zwitterionic dodecylphosphocholine (DPC) micelles, which approximate prokaryotic and eukaryotic membranes, respectively. NMR data indicate the peptide assumes an amphipathic alpha-helical backbone conformation in both micelle environments. However, small differences were observed in the side-chain orientations of lysine, tyrosine, and phenylalanine residues in SDS versus DPC environments. NMR experiments with a paramagnetic probe indicated differences in positioning of the peptide within the two micelle types. Molecular dynamics (MD) simulations of the peptide in both micelle types were also performed to add insight into the peptide/micelle interactions and to assess the validity of this technique to predict the structure of peptides in complex with micelles. MD independently predicted RP-1 to interact only peripherally with the DPC micelle, leaving its spherical shape intact. In contrast, RP-1 entered deeply into and significantly distorted the SDS micelle. Overall, the experimental and MD results support a preferential specificity of RP-1 for anionic membranes over zwitterionic membranes. This specificity likely derives from differences in RP-1 interaction with distinct lipid systems, including subtle differences in side chain orientations, rather than gross changes in RP-1 structure in the two lipid environments.     

Interaction of an artificial antimicrobial peptide with lipid membranes

Antimicrobial peptides constitute an important part of the innate immune defense and are promising new candidates for antibiotics. Naturally occurring antimicrobial peptides often possess hemolytic activity and are not suitable as drugs. Therefore, a range of new synthetic antimicrobial peptides have been developed in recent years with promising properties. But their mechanism of action is in most cases not fully understood. One of these peptides, called V4, is a cyclized 19 amino acid peptide whose amino acid sequence has been modeled upon the hydrophobic/cationic binding pattern found in Factor C of the horseshoe crab (Carcinoscorpius rotundicauda). In this work we used a combination of biophysical techniques to elucidate the mechanism of action of V4. Langmuir-Blodgett trough, atomic force microscopy, Fluorescence Correlation Spectroscopy, Dual Polarization Interference, and confocal microscopy experiments show how the hydrophobic and cationic properties of V4 lead to a) selective binding of the peptide to anionic lipids (POPG) versus zwitterionic lipids (POPC), b) aggregation of vesicles, and above a certain concentration threshold to c) integration of the peptide into the bilayer and finally d) to the disruption of the bilayer structure. The understanding of the mechanism of action of this peptide in relation to the properties of its constituent amino acids is a first step in designing better peptides in the future.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Computational study of solution behavior of magainin 2 monomers

Antimicrobial peptides (AMPs) play crucial role as mediators of the primary host defense against microbial invasion. They are considered a promising alternative to antibiotics for multidrug resistant bacterial strains. For complete understanding of the antimicrobial defense mechanism, a detailed knowledge of the dynamics of peptide-membrane interactions, including atomistic studies on AMPs geometry and both peptide and membrane structural changes during the whole process is a prerequisite. We aim at clarifying the conformation dynamics of small linear AMPs in solution as a first step of in silico protocol for establishing a correspondence between certain amino-acid sequence motifs, secondary-structure elements, conformational dynamics in solution and the intensity and mode of interaction with the bacterial membrane. To this end, we use molecular dynamics simulations augmented by well-tempered metadynamics to study the free-energy landscape of two AMPs with close primary structure and different antibacterial activity – the native magainin 2 (MG2) and an analog (MG2m, with substitutions F5Y and F16W) in aqueous solution. We observe that upon solvation, the initial α-helical structures change differently. The native form remains structured, with three shorter α-helical motifs, connected by random coils, while the synthetic analog tends predominantly to a disordered conformation. Our results indicate the importance of the side-chains at positions 5 and 16 for maintaining the solvated peptide conformation. They also provide a modeling background for recent experimental observations, relating the higher α-helical content in solution (peptide pre-folding) in the case of small linear AMPs to a lower antibacterial activity.     

Protein-drug interactions: characterization of inhibitor binding in complexes of DHFR with trimethoprim and related derivatives

Structural and thermodynamic interactions for the binding of trimethoprim and related congeners to the binary complex of dihydrofolate reductase (from chicken) and NADPH are explored using free energy simulation methods. Good agreement between structures from experimental X-ray refinement and molecular dynamics simulations is found for the complexes. Agreement with thermodynamic measurements is found as well. Our thermodynamic calculations suggest that entropic contributions and desolvation thermodynamics can play a crucial role in overall binding, and that extreme care must be taken in the use of simple model building to rationalize or predict protein-drug binding.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Thermodynamic and kinetic characterization of a beta-hairpin peptide in solution: an extended phase space sampling by molecular dynamics simulations in explicit water

The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Interactions studies of CYP2D6 with quercetin and hyperoside by spectral analysis and molecular dynamics simulations

Some ingredients from herbal medicine can significantly affect the activity of CYP2D6, thus leading to serious interactions between herbs and drugs. Quercetin and hyperoside are active ingredients widely found in vegetables, fruits, and herbal medicines. Quercetin and hyperoside have many biological activities. In this work, the characteristic bindings of CYP2D6 with quercetin/hyperoside are revealed by multi-spectroscopy analysis, molecular docking, and molecular dynamics simulations. The fluorescence of CYP2D6 is statically quenched by quercetin and hyperoside. The binding constant (K_a) values of CYP2D6–quercetin/hyperoside range from 10⁴ L mol⁻¹, which indicates that these two flavonoids bind moderately to CYP2D6. Meanwhile, quercetin has a stronger quenching ability to CYP2D6 than that of hyperoside. The secondary structure of CYP2D6 is obviously changed by binding with quercetin/hyperoside. The docking results reveal that the quercetin/hyperoside enters the active site of CYP2D6 near heme and binds to CYP2D6 by hydrogen bonds and van der Waals forces. The molecular dynamics simulation results indicate that the binding of quercetin/hyperoside can stabilize the two complexes, enhance the flexibility of CYP2D6 backbone atoms, and make a more unfolded and looser structure of CYP2D6.     

Folding is not required for bilayer insertion: replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.     

Thermodynamic and structural analysis of interactions between peptide ligands and SEB

Staphylococcal enterotoxin B (SEB) is an exotoxin produced by Staphylococcus aureus and commonly associated with food poisoning. In this study, SEB‐binding peptides were identified by screening a phage displayed peptide library. The binding of peptides to SEB was tested with isothermal titration calorimetry (ITC) and of the five selected peptides, three showed affinity to SEB, with one measured to have the highest affinity constant (10⁵ M^?1). ITC revealed that the interaction of peptide ligands with SEB was driven entropically and the binding was dominated by hydrophobic interactions. Circular dichroism (CD) measurements and molecular dynamics (MD) simulations, together, give a structural insight into the interaction of peptides with SEB. While SEB binding peptides showed random coil structure before binding, after complex formation they had more ordered structures. The peptide with highest affinity to SEB showed stable conformation during MD simulation. Taken together, our approach about thermodynamic and structural characterization of peptide ligands can be used to develop aptamers, with high affinity and selectivity, for biosensor applications. Copyright © 2009 John Wiley & Sons, Ltd.     

Sampling efficiency in explicit and implicit membrane environments studied by peptide folding simulations

Efficient exploration of the conformational space of peptides embedded in biological membranes is vital to extract converged thermodynamic and kinetic data from simulation studies. Recently developed implicit membrane models promise vast increases in sampling efficiency compared to explicit membrane simulations, allowing for ab initio structure prediction and functional studies. In this study, a previously developed implicit membrane model, based on the generalized Born method, is compared to an explicit di‐palmitoyl‐phosphatidyl‐choline lipid bilayer and an octane slab membrane mimic. The complete folding process of a synthetic 16‐residue peptide is compared using these three setups. Since the comparison requires the entire folding pathway to be captured, individual simulations ranged up to 3 μs of MD. A quantitative sampling comparison using a wide range of performance metrics reveals that the implicit membrane model is at least 2 orders of magnitude more efficient than the simplest explicit setups. Proteins 2009. © 2008 Wiley‐Liss, Inc.