Studies of the binding and structure of adrenocorticotropin peptides in membrane mimics by NMR spectroscopy and pulsed-field gradient diffusion

Scaling Protein-Water Interactions in the Martini 3 Coarse-Grained Force Field to Simulate Transmembrane Helix Dimers in Different Lipid Environments

Martini 3, the latest version of the widely used Martini force field for coarse-grained molecular dynamics simulations, is a promising tool to investigate proteins in phospholipid bilayers. However, simulating other lipid environments, such as detergent micelles, presents challenges due to the absence of validated parameters for their constituent molecules. Here, we propose parameters for the micelle-forming surfactant, dodecylphosphocholine (DPC). These result in micelle assembly with aggregation numbers in agreement with the experimental values. However, we identified a lack of hydrophobic interactions between transmembrane helix protein dimers and the tails of DPC molecules, preventing insertion and stabilization of the protein in the micelles. This problem was also observed for protein insertion by self-assembling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or dipalmitoylphosphatidylcholine (DPPC) bilayers. We propose the reduction of the nonbonded interactions between protein and water beads by 10% as a simple and effective solution to this problem that enables protein encapsulation in phospholipid micelles and bilayers without altering protein dimerization or the bilayer structure.

Three-Dimensional Structure of the Antimicrobial Peptide Cecropin P1 in Dodecylphosphocholine Micelles and the Role of the C-Terminal Residues

Cecropin P1 (CP1) isolated from a large roundworm Ascaris suum, which is found in pig intestines, has been extensively studied as a model antimicrobial peptide (AMP). However, despite being a model AMP, its antibacterial mechanism is not well understood, particularly the function of its C-terminus. By using an Escherichia coli overexpression system with calmodulin as a fusion partner, we succeeded in the mass expression of recombinant peptides, avoiding toxicity to the host and degradation of CP1. The structure of the recombinant ¹⁵N- and ¹³C-labeled CP1 and its C-terminus truncated analogue in dodecylphosphocholine (DPC) micelles was determined by NMR. In this membrane-mimetic environment, CP1 formed an α-helix for almost its entire length, except for a short region at the C-terminus, and there was no evidence of a hinge, which is considered important for the expression of activity in other cecropins. Several NMR analyses showed that the entire length of CP1 was protected from water by micelles. Since the loss of the C-terminus of the analogue had little effect on the NMR structure or its interaction with the micelle, we investigated another role of the C-terminus of CP1 in its antimicrobial activity. The results showed that the C-terminal region affected the DNA-binding capacity of CP1, and this mechanism of action was also newly suggested that it contributed to the antimicrobial activity of CP1.

Finite-Size Effects and Optimal System Sizes in Simulations of Surfactant Micelle Self-Assembly

The spontaneous formation of micelles in aqueous solutions is governed by the amphipathic nature of surfactants and is practically interesting due to the regular use of micelles as membrane mimics, for the characterization of protein structure, and for drug design and delivery. We performed a systematic characterization of the finite-size effect observed in single-component dodecylphosphocholine (DPC) micelles with the coarse-grained MARTINI model. Of multiple coarse-grained solvent models investigated using large system sizes, the nonpolarizable solvent model was found to most accurately reproduce SANS spectra of 100 mM DPC in aqueous solution. We systematically investigated the finite-size effect at constant 100 mM concentration in 23 systems of sizes 40-150 DPC, confirming the finite-size effect to manifest as an oscillation in the mean micelle aggregation number about the thermodynamic aggregation number as the system size increases. The oscillations in aggregation number mostly diminish once the system supports the formation of three micelles. Similar oscillations were observed in the estimated critical micelle concentration with a mean value of 1.10 mM, which is in agreement with experiment to 0.1 mM. The accuracy of using a multiscale simulation approach to avoid finite-size effects in the micelle size distribution and SANS spectra using MARTINI and CHARMM36 was explored using multiple long time scale 500 DPC coarse-grained simulations, which were back-mapped to CHARMM36 all-atom systems. It was found that the MARTINI model generally occupies more volume than the all-atom model, leading to the formation of micelles that are of a reasonable radius of gyration but are smaller in aggregation number. The systematic characterization of the finite-size effect and exploration of multiscale modeling presented in this work provide guidance for the accurate modeling of micelles in simulations.

Determination of Structure and Micellar Interactions of Small Antimicrobial Peptides by Solution-State NMR

NMR spectroscopy is a well-established technique to determine the structure of peptides and small proteins in solution, also when bound to detergent micelles or phospholipid bicelles. The structure of the peptide alone is, however, not conveying the full picture, if the peptide is bound to a micelle, since it does not tell anything about the orientation of the peptide in the micelle. This article describes how to obtain that information together with information on peptide structure.

Diffusion NMR studies of macromolecular complex formation, crowding and confinement in soft materials

Label-free methods to obtain hydrodynamic size from diffusion measurements are desirable in environments that contain multiple macromolecular species at a high total concentration: one example is the crowded cellular environment. In complex, multi-species macromolecular environments - in this article, we feature aqueous systems involving polymers, surfactants and proteins - the link between dynamics and size is harder to unpack due to macromolecular crowding and confinement. In this review, we demonstrate that the pulsed-field gradient NMR technique, with its spectral separation of different chemical components, is ideal for studying the dynamics of the entire system simultaneously and without labelling, in a wide range of systems. The simultaneous measurement of the dynamics of multiple components allows for internal consistency checks and enables quantitative statements about the link between macromolecular dynamics, size, complex formation and crowding in soft materials.

The application of DOSY NMR and molecular dynamics simulations to explore the mechanism(s) of micelle binding of antimicrobial peptides containing unnatural amino acids

Anionic and zwitterionic micelles are often used as simple models for the lipids found in bacterial and mammalian cell membranes to investigate antimicrobial peptide-lipid interactions. In our laboratory we have employed a variety of 1D, 2D, and diffusion ordered (DOSY) NMR experiments to investigate the interactions of antimicrobial peptides containing unnatural amino acids with SDS and DPC micelles. Complete assignment of the proton spectra of these peptides is prohibited by the incorporation of a high percentage of unnatural amino acids which don't contain amide protons into the backbone. However preliminary assignment of the TOCSY spectra of compound 23 in the presence of both micelles indicated multiple conformers are present as a result of binding to these micelles. Chemical Shift Indexing agreed with previously collected CD spectra that indicated on binding to SDS micelles compound 23 adopts a mixture of α-helical structures and on binding to DPC micelles this peptide adopts a mixture of helical and β-turn/sheet like structures. DOSY NMR experiments also indicated that the total positive charge and the relative placement of that charge at the N-terminus or C-terminus are important in determining the mole fraction of the peptide that will bind to the different micelles. DOSY and (1) H-NMR experiments indicated that the length of Spacer #1 plays a major role in defining the binding conformation of these analogs with SDS micelles. Results obtained from molecular simulations studies of the binding of compounds 23 and 36 with SDS micelles were consistent with the observed NMR results.

Photo-Fries rearrangement of aryl acetamides: regioselectivity induced by the aqueous micellar green environment

NMR spectroscopy shows the location of acetanilides within the shells and hydrophobic cores in micellar solutions. Irradiation of acetanilides in aqueous micellar solutions involves C–N homolytic cleavage to yield singlet radical pairs that selectively provide 2-aminoacetophenone derivatives.

CHARMM36 united atom chain model for lipids and surfactants

Molecular simulations of lipids and surfactants require accurate parameters to reproduce and predict experimental properties. Previously, a united atom (UA) chain model was developed for the CHARMM27/27r lipids (Hénin, J., et al. J. Phys. Chem. B. 2008, 112, 7008-7015) but suffers from the flaw that bilayer simulations using the model require an imposed surface area ensemble, which limits its use to pure bilayer systems. A UA-chain model has been developed based on the CHARMM36 (C36) all-atom lipid parameters, termed C36-UA, and agreed well with bulk, lipid membrane, and micelle formation of a surfactant. Molecular dynamics (MD) simulations of alkanes (heptane and pentadecane) were used to test the validity of C36-UA on density, heat of vaporization, and liquid self-diffusion constants. Then, simulations using C36-UA resulted in accurate properties (surface area per lipid, X-ray and neutron form factors, and chain order parameters) of various saturated- and unsaturated-chain bilayers. When mixed with the all-atom cholesterol model and tested with a series of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)/cholesterol mixtures, the C36-UA model performed well. Simulations of self-assembly of a surfactant (dodecylphosphocholine, DPC) using C36-UA suggest an aggregation number of 53 ± 11 DPC molecules at 0.45 M of DPC, which agrees well with experimental estimates. Therefore, the C36-UA force field offers a useful alternative to the all-atom C36 lipid force field by requiring less computational cost while still maintaining the same level of accuracy, which may prove useful for large systems with proteins.

Peptide-lipid interactions: experiments and applications

The interactions between peptides and lipids are of fundamental importance in the functioning of numerous membrane-mediated cellular processes including antimicrobial peptide action, hormone-receptor interactions, drug bioavailability across the blood-brain barrier and viral fusion processes. Moreover, a major goal of modern biotechnology is obtaining new potent pharmaceutical agents whose biological action is dependent on the binding of peptides to lipid-bilayers. Several issues need to be addressed such as secondary structure, orientation, oligomerization and localization inside the membrane. At the same time, the structural effects which the peptides cause on the lipid bilayer are important for the interactions and need to be elucidated. The structural characterization of membrane active peptides in membranes is a harsh experimental challenge. It is in fact accepted that no single experimental technique can give a complete structural picture of the interaction, but rather a combination of different techniques is necessary.

CHARMM-GUI micelle builder for pure/mixed micelle and protein/micelle complex systems

Micelle Builder in CHARMM-GUI, http://www.charmm-gui.org/input/micelle , is a web-based graphical user interface to build pure/mixed micelle and protein/micelle complex systems for molecular dynamics (MD) simulation. The robustness of Micelle Builder is tested by simulating four detergent-only homogeneous micelles of DHPC (dihexanoylphosphatidylcholine), DPC (dodecylphosphocholine), TPC (tetradecylphosphocholine), and SDS (sodium dodecyl sulfate) and comparing the calculated micelle properties with experiments and previous simulations. As a representative protein/micelle model, Pf1 coat protein is modeled and simulated in DHPC micelles with three different numbers of DHPC molecules. While the number of DHPC molecules in direct contact with Pf1 protein converges during the simulation, distinct behavior and geometry of micelles lead to different protein conformations in comparison to that in bilayers. It is our hope that CHARMM-GUI Micelle Builder can be used for simulation studies of various protein/micelle systems to better understand the protein structure and dynamics in micelles as well as distribution of detergents and their dynamics around proteins.

Understanding the interaction between valsartan and detergents by NMR techniques and molecular dynamics simulation

Valsartan (VST) is one of the Angiotensin II receptor antagonists, which is widely used in clinical hypertension treatment. It is believed that VST incorporates into biological membranes before it binds to AT(1) receptor. Herein the interactions between VST and detergents, mimicking the membrane environment, were investigated by using nuclear magnetic resonance (NMR) techniques and molecular dynamics (MD) simulation. We observed that VST has two conformers (trans and cis) exchanging slowly in DPC (dodecyl-phosphocholine) micelles, a widely used detergent. The changes of chemical shifts, relaxation rates, and self-diffusion coefficients of VST protons indicate that both conformers have strong interactions with DPC. NOE cross peaks and MD simulation reveal that DPC interacts with VST not only through the hydrophobic lipid chain, but also the hydrophilic headgroup, locating VST at the charged headgroup and upper part of the micelles. Our results are in good agreement with the Raman spectroscopic studies of VST in the DPPC (dipalmitoyl-phosphatidylcholine) bilayers by Potamitis et al. (Biochim. Biophys. Acta. 2011). The concentration ratio of trans over cis conformers is 0.94, showing that two conformers have the same affinities with the detergent, which is significantly smaller than our previous results obtained in SDS (sodium dodecyl sulfate) micelles. MD simulation suggested that the cis conformer has slightly lower binding free energy than the trans conformer when interacting with DPC. The conformational change of VST was further investigated in two detergents, CTAB (hexadecyltrimethylammonium bromide) and Tween-20 (polysorbate 20). Ratios of conformer A and B in the presence of detergents are in the order of DPC, CTAB < Tween-20 < SDS, which is correlated with the charge characters of their head groups. NMR investigations and MD simulations indicate that the electrostatic interaction plays an essential role in the binding process of VST with detergents, and the hydrophobic interaction influences the packing of the drug in the micelles. These results may be of help in understanding delivery processes of sartan drugs in cell membranes.

Structure and orientation of the gH625-644 membrane interacting region of herpes simplex virus type 1 in a membrane mimetic system

Glycoprotein H (gH) of the herpes simplex virus type 1 is involved in the complex mechanism of membrane fusion of the viral envelope with host cells. The virus requires four glycoproteins (gB, gD, gH, gL) to execute fusion and the role played by gH remains mysterious. Mutational studies have revealed several regions of gH ectodomain required for fusion and identified the segment from amino acid 625 to 644 as the most fusogenic region. Here, we studied the behavior in a membrane-mimicking DPC micellar environment of a peptide encompassing this region (gH625-644) and determined its NMR solution structure and its orientation within the micelles.

Membrane catalysis of peptide-receptor binding

The membrane catalysis hypothesis states that a peptide ligand activates its target receptor after an initial interaction with the surrounding membrane. Upon membrane binding and interaction, the ligand is structured such that receptor binding and activation is encouraged. As evidence for this hypothesis, there are numerous studies concerning the conformation that peptides adopt in membrane mimetic environments. This mini-review analyzes the features of ligand peptides with an available high-resolution membrane-induced structure and a characterized membrane-binding region. At the peptide-membrane interface, both amphipathic helices and turn structures are commonly formed in peptide ligands and both hydrophobic and electrostatic interactions can be responsible for membrane binding. Apelin is the ligand to the G-protein coupled receptor (GPCR) named APJ, with various important physiological effects, which we have recently characterized both in solution and bound to anionic micelles. The structural changes that apelin undergoes when binding to micelles provide strong evidence for membrane catalysis of apelin-APJ interactions.

Solution structure of the HsapBK K+ channel voltage-sensor paddle sequence

Voltage-gated potassium channels open and close in response to changes in the membrane potential. In this study, we have determined the NMR solution structure of the putative S3b-S4 voltage-sensor paddle fragment, the part that moves to mediate voltage gating, of the HsapBK potassium channel in dodecylphosphocholine (DPC) micelles. This paper presents the first structure of the S3b-S4 fragment from a BK channel. Diffusion coefficients as determined from PFG NMR experiments showed that a well-defined complex between the peptide and DPC molecules was formed. The structure reveals a helix-turn-helix motif, which is in agreement with crystal structures of other voltage-gated potassium channels, thus indicating that it is feasible to study the isolated fragment. The paddle motifs generally contain several basic residues, implicated in the gating. The critical Arg residues in this structure all reside on the surface, which is in agreement with crystal structures of K(v) channels. Similarities in the structure of the S3b-S4 fragment in BK and K(v) channels as well as important differences are seen, which may be important for explaining the details in paddle movement within a bilayer.

Addition of a cholesterol group to an HIV-1 peptide fusion inhibitor dramatically increases its antiviral potency

Peptides derived from the heptad repeat 2 (HR2) region of the HIV fusogenic protein gp41 are potent inhibitors of viral infection, and one of them, enfuvirtide, is used for the treatment of therapy-experienced AIDS patients. The mechanism of action of these peptides is binding to a critical intermediate along the virus-cell fusion pathway, and accordingly, increasing the affinity for the intermediate yields more potent inhibitors. We took a different approach, namely to increase the potency of the HR2 peptide inhibitor C34 by targeting it to the cell compartment where fusion occurs, and we show here that a simple, yet powerful way to accomplish this is attachment of a cholesterol group. C34 derivatized with cholesterol (C34-Chol) shows dramatically increased antiviral potency on a panel of primary isolates, with IC(90) values 15- to 300-fold lower than enfuvirtide and the second-generation inhibitor T1249, making C34-Chol the most potent HIV fusion inhibitor to date. Consistent with its anticipated mechanism of action, the antiviral activity of C34-Chol is unusually persistent: washing target cells after incubation with C34-Chol, but before triggering fusion, increases IC(50) only 7-fold, relative to a 400-fold increase observed for C34. Moreover, derivatization with cholesterol extends the half-life of the peptide in vivo. In the mouse, s.c. administration of 3.5 mg/kg C34-Chol yields a plasma concentration 24 h after injection >300-fold higher than the measured IC(90) values. Because the fusion machinery targeted by C34-Chol is similar in several other enveloped viruses, we believe that these findings may be of general utility.

Implicit solvent simulations of DPC micelle formation

The formation of micelles by dodecylphosphocholine (DPC) is modeled by treating the surfactants in atomic detail and the solvent implicitly, in the spirit of the EEF1 solvation model for proteins. The solvation parameters of the DPC atoms are carried over from those of similar atoms in proteins. A slight adjustment of the parameters for the headgroup was found necessary for obtaining an aggregation number consistent with experiment. Molecular dynamics simulations of 960 DPC molecules at different concentrations are used to obtain the aggregation number, the micelle size distribution, and the CMC. At 20 mM concentration we obtain an aggregation number of 53-56 and a CMC of 1.25 mM, values close to the experimental ones. At 100 mM the aggregation number increases to 90. Simulations of individual micelles of varying size show that the effective energy per surfactant molecule is initially a decreasing function of aggregation number but stabilizes at about 60 molecules. The van der Waals term and the desolvation of nonpolar groups contribute to micellization, whereas the desolvation of polar groups opposes it. From the difference between the effective energy and the free energy (calculated from the CMC), the translational and rotational entropy contributions to the free energy are estimated at about 7 kcal/mol per monomer. The micelles obtained here are more irregular than those obtained in explicit water simulations. This modeling approach allows the study of larger surfactant aggregates for longer times and the extraction of thermodynamic in addition to structural information.

A Model for Cell-Surface-Exposed Carbohydrate Moieties Suitable for Structural Studies by NMR Spectroscopy

In the present study a synthetic glycolipid system is presented that can be readily incorporated into phospholipid micelles and that allows the study of cell-surface-exposed carbohydrate units by high-resolution NMR techniques. Here, we present an efficient route for the synthesis of glycolipid compounds that contain mannose, mannobiose, or mannotriose coupled either directly to an alkyl chain or through a poly(ethylene glycol) linker. Furthermore, we have validated our model system by measuring the binding of cyanovirin N (CV-N), a cyanobacterial protein that binds with nanomolar affinity to the terminal arms of high-mannose structures of the HIV surface-envelope glycoprotein gp120, to glycolipids the carbohydrate portions of which comprise the corresponding high-mannose moieties. From the results of chemical-shift mapping with uniformly (15)N-labelled CV-N, we conclude that binding to the protein occurs at sites similar to those involved in binding the nonconjugated carbohydrates. We characterized the insertion of the glycolipids into dodecylphosphocholine (DPC) micelles by measuring translational diffusion, and we observed that the diffusion constants of the glycolipids were very similar to those of the DPC micelles themselves, but significantly deviated from those of the free glycolipids. We also present experimental proof that the glycolipids remain inserted in the micelles while binding to CV-N. Finally, by addition of a ligand that had a higher affinity to CV-N but which was not attached did not couple to a lipid anchor, CV-N could be released from the glycolipid and, hence, from the micelle-associated state.

Triggering and visualizing the aggregation and fusion of lipid membranes in microfluidic chambers

We present a method that makes it possible to trigger, observe, and quantify membrane aggregation and fusion of giant liposomes in microfluidic chambers. Using electroformation from spin-coated films of lipids on transparent indium tin oxide electrodes, we formed two-dimensional networks of closely packed, surface-attached giant liposomes. We investigated the effects of fusogenic agents by simply flowing these molecules into the chambers and analyzing the resulting shape changes of more than 100 liposomes in parallel. We used this setup to quantify membrane fusion by several well-studied mechanisms, including fusion triggered by Ca2+, polyethylene glycol, and biospecific tethering. Directly observing many liposomes simultaneously proved particularly useful for studying fusion events in the presence of low concentrations of fusogenic agents, when fusion was rare and probabilistic. We applied this microfluidic fusion assay to investigate a novel 30-mer peptide derived from a recently identified human receptor protein, B5, that is important for membrane fusion during the entry of herpes simplex virus into host cells. This peptide triggered fusion of liposomes at an approximately 6 times higher probability than control peptides and caused irreversible interactions between adjacent membranes; it was, however, less fusogenic than Ca2+ at comparable concentrations. Closely packed, surface-attached giant liposomes in microfluidic chambers offer a method to observe membrane aggregation and fusion in parallel without requiring the use of micromanipulators. This technique makes it possible to characterize rapidly novel fusogenic agents under well-defined conditions.

A pulsed field gradient NMR diffusion investigation of enkephalin peptide-sodium dodecyl sulfate micelle association

Pulsed field gradient NMR (PFG-NMR) diffusion experiments were used to investigate the binding of leucine and methionine enkephalin peptides to anionic sodium dodecyl sulfate (SDS) micelles. The study was undertaken to investigate the mechanism of interaction between enkephalin peptides and SDS micelles and to determine if NMR-derived association constants, K(eq), can predict the elution order in electrokinetic chromatography (EKC). In EKC, peptides are separated on the basis of their interactions with micelles. The Leu-enkephalin peptide-micelle association constant increased from 130 +/- 8 to 1459 +/- 57 and 1744 +/- 64 M(-1), respectively, when an Arg or Lys was added to the C-terminus. The association constant of Leu-enkephalinamide was approximately equal to that of Leu-enkephalin-Arg. Substitution of Phe4 with a Trp or Gly2 with an Ala in the Leu-enkephalin peptides also increased the micelle binding affinity. These results confirm that the interaction of Leu-enkephalin peptides with SDS micelles is largely electrostatic and that the non-polar amino acid side chains interact with the hydrophobic micelle core. The peptide-micelle association constants for the cationic Met-enkephalin peptides were also greater than their zwitterionic counterparts. For example, the Met-enkephalin K(eq) value was 162 +/- 9 M(-1), while the association constants for Met-enkephalin-Arg, Met-enkephalin-Lys, and Met-enkephalinamide were, respectively, 674 +/- 31, 426 +/- 23, and 453 +/- 27 M(-1). In both Met-enkephalin and Met-enkephalinamide, replacing Gly2 with an Ala did not significantly increase the association constant. These results confirm that with the Met-enkephalin peptides, there was little or no interaction of the amino acid side chains with the micelle core. For both the Leu-enkephalin and Met-enkephalin peptides, the association constants were consistent with EKC results, in that the peptides with smaller K(eq) values were found to elute before those with larger association constants.

NMR Conformational Analyses on (des-bromo) Neuropeptide B [1–23] and Neuropeptide W [1–23]: The Importance of α-helices, a Cation-π Interaction and a β-Turn

The preferred conformations of the orphan G-protein coupled receptor agonists (des-bromo) neuropeptide B [1-23] and neuropeptide W [1-23], referred to as NPB and NPW, have been determined by (1)H NMR, CD, and molecular modeling. The sequences of NPB and NPW are WYKPAAGHSSYSVGRAAGLLSGL and WYKHVASPRYHTVGRAAGLLMGL, respectively. These are hypothalamic peptides that exert their biological actions on GPR7 and GPR8 receptors. Micellar solutions using the membrane mimetic, sodium dodecylsulphate-d(25) (SDS), were used to mimic a physiological environment for the peptides. The secondary structure of NPB consists of a type II beta-turn involving residues Lys(3) to Ala(6). The C-terminal region of NPB exists in a conformational equilibrium between different secondary structures, including an alpha-helix from residues Arg(15) to Ser(21), and a 3(10)-helix from residues Ser(12) to Ser(21). The N-terminus of NPW exhibits a cation-pi interaction between the Lys(3) side chain and the quadrupole moment of the Trp(1) indole group. At the C-terminus of NPW, a well-defined alpha-helical conformation exists from Arg(15) to Met(21). As NPB and NPW have 91% sequence homology from residues Val(13) to Leu(23), with only residue 21 differing between the two peptides, the similar C-terminal secondary structures of these two peptides are consistent with the sequences. This is supported by the similar CD spectra. The different secondary structures at the N-termini for NPB and NPW point to the importance of the N-terminus in receptor binding. This is consistent with the work of Fujii et al. [J. Biol. Chem. 277, 34010-34016 (2002)] who observed that iodination of the NPB Tyr(2) resulted in decreased agonistic activity at GPR7. In addition, Tanaka et al. [Proc. Natl. Acad. Sci. USA 100, 6251-6256 (2003)] showed that deletion of Trp(1) from NPB or NPW drastically decreased activity at GPR7 for NPB and GPR7 and GPR8 for NPW. Therefore, we postulate that the N-terminus is involved in membrane recognition and receptor binding.

The interactions of the HIV gp41 fusion peptides with zwitterionic membrane mimics determined by NMR spectroscopy

The wild-type (wt) N-terminal 23-residue fusion peptide (FP) of the human immunodeficiency virus (HIV) fusion protein gp41 and its V2E mutant have been studied by nuclear magnetic resonance (NMR) spectroscopy in dodecylphosphocholine (DPC) micelles as membrane mimics. A number of NMR techniques have been used. Pulsed field-gradient diffusion measurements in DPC and in 4:1 DPC/sodium dodecylsulfate mixed micelles showed that there is no major difference between the partition coefficients of the fusogenic wt peptide and the V2E mutant in these micelles, indicating that there is no correlation between the activity of the fusion peptides and their membrane affinities. The nuclear Overhauser enhancement (NOE) patterns and the chemical shift index for these two peptides indicated that both FP are in an alpha helical conformation between the Ile4 to Leu12 or to Ala15 region. Simulated annealing showed that the helical region extends from Ile4 to Met19. The two FPs share similar conformational characteristics, indicating that the conformation of the FP is not an important factor determining its activity. The spin-label studies, utilizing spin labels 5- and 16-doxystearic acids in the DPC micelles, provided clear indication that the wt FP inserts its N-terminus into the micelles while the V2E mutant does not insert into the micelles. The conclusion from the spin-label results is corroborated by deuterium amide proton exchange experiments. The correlation between the oblique insertion of the FP and its fusogenic activity is in excellent agreement with results from our molecular dynamics simulation and from other previous studies.

Interaction and structural study of kinin peptide bradykinin and ganglioside monosialylated 1 micelle

Partitioning of small proteins and peptides from the aqueous to membrane phase is often coupled with folding. In this work we examine the binding and folding of the kinin peptide, bradykinin (BK), in the presence of the ganglioside monosialylated 1 (GM1) micelle. Using two-dimensional NMR techniques, we have shown that at low concentration, GM1 micelle is able to induce a turn conformation to BK. A pulsed-field gradient diffusion NMR study indicated that the peptide partitions into the GM1 micelle with a DeltaG(part) of -3.14 +/- 0.03 kcal/mol. A saturation transfer difference (STD) NMR study indicated that the binding is mostly through hydrophobic residues.

Gαs proteinC-terminal α-helix at the interface: does the plasma membrane play a critical role in the Gαs protein functionality?

The heterotrimeric guanine nucleotide-binding regulatory proteins (G proteins, Galphabetagamma) mediate the signalling process of a large number of receptors, known as G protein-coupled receptors. The C-terminal domain of the heterotrimeric G protein alpha-subunit plays a key role in the selective activation of G proteins by their cognate receptors. The interaction of this domain can take place at the end of a cascade including several successive conformational modifications. Galpha(s)(350-394) is the 45-mer peptide corresponding to the C-terminal region of the Galpha(s) subunit. In the crystal structure of the Galpha(s) subunit it encompasses the alpha4/beta6 loop, the beta6 beta-sheet segment and the alpha5 helix region. Following a previous study based on the synthesis, biological activity and conformational analysis of shorter peptides belonging to the same Galpha(s) region, Galpha(s)(350-394) was synthesized and investigated. The present study outlines the central role played by the residues involved in the alpha4/beta6 loop and beta6/alpha5 loops in the stabilization of the C-terminal Galpha(s)alpha-helix. H(2)O/(2)H(2)O exchange experiments, and NMR diffusion experiments show interesting evidence concerning the interaction between the SDS micelles and the polypeptide. These data prompt intriguing speculations on the role of the intracellular environment/cellular membrane interface in the stabilization and functionality of the C-terminal Galpha(s) region.

Aqueous and Micelle-bound Structural Characterization of the Epidermal Growth Factor Receptor Juxtamembrane Domain Containing Basolateral Sorting Motifs

The EGF receptor is the prototype for four highly related receptors constituting the ErbB family. The EGF receptor is normally targeted to the basolateral membrane in polarized epithelial cells, where it relays information from underlying tissues. Two basolateral sorting signals have been mapped to the cytoplasmic juxtamembrane region of the receptor, a dominant signal comprised of a polyproline core (667-PXXP) and a preceding basic residue (Arg662), and a consensus leucine-based signal (658-LL) responsible for residual sorting when the 667-PXXP signal is absent or defective. The goal of this study was to define the structure of these signals, and gain some insights into how these structures might be regulated by cellular microenvironment. Structural information was acquired for two peptides corresponding to EGF receptor residues Arg645 and Ala674 in aqueous solution or in the presence of membrane-mimicking dodecylphosphocholine micelles, using a variety of NMR and CD spectroscopic methods. Chemical shift data indicate that the 667-PXXP signal does not bind to the micelles and is in random coil state in both aqueous solution and a micellar environment, raising the possibility that 667-PXXP switches to an ordered structure during interaction with the basolateral sorting machinery. In contrast, the adjacent region including 658-LL does bind to micelles mediated by a highly positively charged region located between Arg645 and Arg656. The micelle-bound region also includes Thr654, a known substrate for PKC. This suggests a distinct mode of regulation for this signal involving membrane association and/or phosphorylation.

A multidimensional 1H NMR investigation of the conformation of methionine-enkephalin in fast-tumbling bicelles

Enkephalins are pentapeptides found in the central nervous system. It is believed that these neuropeptides interact with the nerve cell membrane to adopt a conformation suitable for their binding to an opiate receptor. In this work, we have determined the three-dimensional structure of methionine-enkephalin (Menk) in fast-tumbling bicelles using multidimensional (1)H NMR. Bicelles were selected as model membranes because both their bilayer organization and composition resemble those of natural biomembranes. The effect of the membrane composition on the peptide conformation was explored using both zwitterionic (PC bicelles) and negatively charged bicelles (Bic/PG). Pulsed field gradient experiments allowed the determination of the proportion of Menk bound to the model membranes. Approximately 60% of the water-soluble enkephalin was found to associate to the bicellar systems. Structure calculations from torsion angle and NOE-based distance constraints suggest the presence of both micro - and delta-selective conformers of Menk in each system and slightly different conformers in PC bicelles and Bic/PG. As opposed to previous studies of enkephalins in membrane mimetic systems, our results show that these opiate peptides could adopt several conformations in a membrane environment, which is consistent with the flexibility and poor selectivity of enkephalins.

PFG-NMR Investigations of the Binding of Cationic Neuropeptides to Anionic and Zwitterionic Micelles

The mechanism by which peptides bind to micelles is believed to be a two-phase process, involving (i). initial electrostatic interactions between the peptide and micelle surface, followed by (ii). hydrophobic interactions between peptide side chains and the micelle core. To better characterize the electrostatic portion of this process, a series of pulse field gradient nuclear magnetic resonance (PFG-NMR) spectroscopic experiments were conducted on a group of neuropeptides with varying net cationic charges (+1 to +3) and charge location to determine both their diffusion coefficients and partition coefficients when in the presence of detergent micelles. Two types of micelles were chosen for the study, namely anionic sodium dodecylsulfate (SDS) and zwitterionic dodecylphosphocholine (DPC) micelles. Results obtained from this investigation indicate that in the case of the anionic SDS micelles, peptides with a larger net positive charge bind to a greater extent than those with a lesser net positive charge (bradykinin > substance P > neurokinin A > Met-enkephalin). In contrast, when in the presence of zwitterionic DPC micelles, the degree of mixed-charge nature of the peptide affects binding (neurokinin A > substance P > Met-enkephalin > bradykinin). Partition coefficients between the peptides and the micelles follow similar trends for both micelle types. Diffusion coefficients for the peptides in SDS micelles, when ranked from largest to smallest, follow a trend where increasing net positive charge results in the smallest diffusion coefficient: Met-enkephalin > neurokinin A > bradykinin > substance P. Diffusion coefficients when in the presence of DPC micelles, when ranked from largest to smallest, follow a trend where the presence of negatively-charged side chains results in the smallest diffusion coefficient: bradykinin > Met-enkephalin > substance P > neurokinin A.

Peptide conformational changes induced by tryptophan-phosphocholine interactions in a micelle

Sodium dodecylsulfate (SDS) and dodecylphosphocholine (DPC) micelles are often used to mimic the membrane- or receptor-bound states of peptides in NMR studies. From the present examination of a 26-residue analog of exendin-4 (TrEX4) by NMR and CD in water, aqueous 30% trifluoroethanol (TFE), and bound to both SDS and DPC micelles, it is clear that these two lipid micelles can yield very different peptide structures. The Trp-cage fold (also observed in 30% TFE) is present when TrEX4 is bound to SDS micelles; however, tertiary structure is absent in the presence of DPC micelles. The loss of tertiary structure is attributed to an energetically favorable interaction (estimated as 2-3 kcal/mol) of the tryptophan side chain with the phosphocholine head groups. These dramatic structural differences suggest that care must be taken when using either SDS or DPC to mimic the membrane- or receptor-bound states.

Molecular dynamics simulations of adrenocorticotropin (1-24) peptide in a solvated dodecylphosphocholine (DPC) micelle and in a dimyistoylphosphatidylcholine (DMPC) bilayer

The structure and interactions of the 1-24 fragment of the adrenocorticotropin hormone, ACTH (1-24), with membrane have been studied by molecular dynamics (MD) simulation in an NPT ensembles in two explicit membrane mimics, a dodecylphosphocholine (DPC) micelle and a dimyristoylphosphatidylcholine (DMPC) bilayer. The starting configuration of the peptide/lipid systems had the 1-10 segment of the peptide lying on the surface of the model membrane, the same as the equilibrated structure (by MD) of ACTH (1-10) in a DPC micelle. The simulations showed that the peptide adopts the surface-binding mode and essentially the same structure in both systems. Thus the results of this work lend support to the assumption that micelles are reasonable mimics for biological membranes for the study of peptide binding. The 1-10 segment is slightly tilted from the parallel orientation to the interface and interacts strongly with the membrane surface while the more polar 11-24 segment shows little tendency to interact with the membrane surface, preferring to reside primarily in the aqueous phase. Furthermore, the 1-10 segment of the peptide binds to the DPC micelle in essentially the same way as ACTH (1-10). Thus the MD results are in excellent agreement with the model of interaction of ACTH (1-24) with membrane derived from NMR experiments. The secondary structure and the hydration of the peptide and the interactions of specific residues with the lipid head groups have also been analyzed.

Molecular dynamics simulation of adrenocorticotropin (1-10) peptide in a solvated dodecylphosphocholine micelle

Adrenocorticotropin (ACTH) (1-10), an adrenocorticotropin hormone fragment, has been studied by molecular dynamics (MD) simulation in an NPT ensemble in an explicit dodecylphosphocholine (DPC) micelle. Two starting configurations of the peptide/micelle system, corresponding to the insertion and surface-binding modes, were used. A common equilibrated configuration, in which the peptide lies parallel to the micellar surface, was reached from both simulations. In the initial part of the simulations, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated before the peptide reached an equilibrium configuration with respect to the micelle. Analyses of the trajectories from the subsequent free (unrestrained) MD simulation showed that ACTH (1-10) does not conform strictly to a helical structure. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide protons of the peptide and the micellar head groups. However, the extent of the latter interaction is less pronounced than in the negatively charged SDS micelle. The final structure enhances the amphipathic nature of the peptide, facilitating better interactions at the water-hydrophobic interface. The primary hydrophobic interactions with the micelle came from the side chains of Met4, Phe7, and Trp9. All peptide bonds were either hydrated or were involved in intramolecular hydrogen bonding. The interactions with the DPC micelle, the conformation of the bound peptide, and the dynamics of the peptide, as revealed by the time correlation functions of the N-H bonds, were compared with those of the ACTH (1-10)/SDS system studied previously by MD simulations.

Effects of the incorporation of CHAPS into SDS micelles on neuropeptide-micelle binding: separation of the role of electrostatic interactions from hydrophobic interactions

It is well known that neuropeptides interact with lipid vesicles in a manner similar to biological membranes, with electrostatic interactions between the two providing a mechanism for concentrating the peptide at the vesicle's surface, followed by hydrophobic interactions between the peptide and the core of the vesicle that induce and stabilize secondary structure motifs. In an effort to understand these interactions to a greater extent, our group has developed a series of anionic micelles (SDS) containing various concentrations of the bile salt CHAPS, which is used as a model for cholesterol. The incorporation of CHAPS into the hydrophobic core of these micelles should alter the degree to which the neuropeptide can insert itself, affecting structure. These interactions were investigated using two-dimensional NMR, pulse-field gradient (PFG) NMR, and molecular modeling experiments. The results of this study clearly indicate that electrostatic and hydrophobic interactions between the micelle and neuropeptide are completely independent of one another. Increasing the concentration of CHAPS to 15 mM in the micelles blocks the insertion of the hydrophobic side chains of the neuropeptide into the hydrophobic core of the micelles. The electrostatic interactions as determined by diffusion measurements are not affected by the presence of increasing CHAPS concentration. Our observations are consistent with the predictions of Seelig (A. Seelig and J. Seelig, "Interaction of Drugs and Peptides with the Lipid Membrane," in Structure and Function of 7TM Receptors, T. W. Schwartz, S. A. Hjorth, and T. S. Kastrup, Eds., Munksgaard: Location, 1996).

NMR studies of adrenocorticotropin hormone peptides in sodium dodecylsulfate and dodecylphosphocholine micelles: proline isomerism and interactions of the peptides with micelles

Three adrenocorticotropin hormone (ACTH) fragments (1-10, 1-24, and 11-24) have been studied in water and in sodium dodecylsulfate (SDS) and dodecylphosphocholine (DPC) micelles by nuclear magnetic resonance spectroscopy. The trans-cis isomerism at all three proline sites (at positions 12, 19, and 24) was found in the 11-24 segment of the peptide. The population of the cis isomers changes with the environment of the peptide. Specifically, the presence of the DPC micelle does not affect the trans-cis equilibrium in the 11-24 segment from that in water. In contrast, the presence of the SDS micelles decreases the population of the cis isomer at Pro(24), but increases its population at Pro(12) and Pro(19). The effect of SDS micelles on the trans-cis equilibrium at these proline sites was discussed. Intermolecular nuclear Overhauser effect (NOE) correlations between the ACTH peptides and the micelles were observed. These correlations occurred only in the 1-10 segment of the peptides, and the hydrophobic side chains contributed most to the intermolecular NOE. The intermolecular NOE pattern corroborates the suggestion that the 1-10 segment of the ACTH peptides bind to these micelles via a surface-binding mode, with most of the interactions coming from the insertion of the hydrophobic side chains.

The structure and dynamics of ACTH (1-10) on the surface of a sodium dodecylsulfate (SDS) micelle: a molecular dynamics simulation study

ACTH (1-10), an adrenocorticotropin hormone fragment, was studied by molecular dynamics (MD) simulation in the NPT ensemble in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated during the equilibration stage of the simulation. The analyses of the trajectories from the subsequent unrestrained MD showed that ACTH (1-10) does not conform to a helical structure at the micelle-water interface; however, the structure is amphipathic. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide hydrogens of the peptide and the micelle head-groups. ACTH (1-10) was found to lie on the surface of the SDS micelle. Most of the hydrophobic interactions came from the side-chains of Met-4, Phe-7 and Trp-9. The peptide bonds were either hydrated or involved in intramolecular hydrogen bonding. Decreased hydration for the backbone of His-6 and Phe-7 was due to intermolecular hydrogen bonding with the SDS head-groups. The time correlation functions of the N-H bonds of the peptide in water and in the micelle showed that the motions of the peptide, except for the N- and C-termini, are significantly reduced when partitioned in the micelle.

The study of the conformation and interaction of two tachykinin peptides in membrane mimicking systems by NMR spectroscopy and pulsed field gradient diffusion

Pulsed-field gradient diffusion has been used to study the binding of two tachykinin peptides, [Tyr(8)]-substance P (SP) and [Tyr(0)]-neurokinin A (NKA) to two membrane-mimicking micelles, dodecylphosphocholine, and sodium dodecylsulfate. The structure of these peptides bound to the micelles have also been studied by using two-dimensional nmr and restrained simulated annealing calculations. No major difference in the structures of each peptide in the two micellar media was found. The difference between the micelle-bound structure of [Tyr(8)]SP and that of SP was also minor. The longer helical conformation on the C-terminus for [Tyr(0)]NKA was observed, compared with that for NKA. The relationship between the difference in the biological potencies of [Tyr(8)]SP and SP and the differences in their structure, especially the interaction of the side chains of the two aromatic residues, and the difference in their binding affinities to membrane was discussed. In addition, differences between the result of restrained molecular dynamics simulations of [Tyr(8)]SP in the presence of an explicit micelle and the present results were observed and discussed.

Molecular dynamics study of substance P peptides in a biphasic membrane mimic

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.

A membrane setting for the sorting motifs present in the adenovirus E3-13.7 protein which down-regulates the epidermal growth factor receptor

The adenovirus E3-13.7 protein interferes with endosomal protein sorting to down-regulate the epidermal growth factor receptor and related tyrosine kinase receptors. The cytoplasmic C terminus of this protein contains three protein sorting motifs which are related to the function of E3-13.7. In this study, the structure of a 23-residue polypeptide corresponding to this domain was examined using solution NMR and CD spectroscopic methods. The peptide was observed to exist in a mostly random structural state in aqueous solution but underwent high affinity association with dodecylphosphocholine micelles, where it adopted an ordered structure. The affinity of this peptide for the micellar surface and the structure of the bound peptide were independent of pH variation, surface charge, or attachment of a myristoyl anchor to the N-terminal. Studies with phospholipid vesicles suggested that the micellar structural results can be extrapolated to a true lipid bilayer. On the micellar surface all three sorting motifs are closely associated with the water/apolar interface: 72-YLRH and 87-LL lie within interfacial amphipathic helices, while 76-HPQY is non-helical and dimples just above the surface. These results contribute to the development of an understanding of the basis for specificity in recognition of sorting motifs by components of the cellular protein trafficking machinery.