Interactions of calmodulin with death-associated protein kinase peptides: experimental and modeling studies

We have studied the interactions between calmodulin (CaM) and three target peptides from the death-associated protein kinase (DAPK) protein family using both experimental and modeling methods, aimed at determining the details of the underlying biological regulation mechanisms. Experimentally, calorimetric binding free energies were determined for the complexes of CaM with peptides representing the DAPK2 wild-type and S308D mutant, as well as DAPK1. The observed affinity of CaM was very similar for all three studied peptides. The DAPK2 and DAPK1 peptides differ significantly in sequence and total charge, while the DAPK2 S308D mutant is designed to model the effects of DAPK2 Ser308 phosphorylation. The crystal structure of the CaM-DAPK2 S308D mutant peptide is also reported. The structures of CaM-DAPK peptide complexes present a mode of CaM-kinase interaction, in which bulky hydrophobic residues at positions 10 and 14 are both bound to the same hydrophobic cleft. To explain the microscopic effects underlying these interactions, we performed free energy calculations based on the approximate MM-PBSA approach. For these highly charged systems, standard MM-PBSA calculations did not yield satisfactory results. We proposed a rational modification of the approach which led to reasonable predictions of binding free energies. All three complexes are strongly stabilized by two effects: electrostatic interactions and buried surface area. The strong favorable interactions are to a large part compensated by unfavorable entropic terms, in which vibrational entropy is the largest contributor. The electrostatic component of the binding free energy followed the trend of the overall peptide charge, with strongest interactions for DAPK1 and weakest for the DAPK2 mutant. The electrostatics was dominated by interactions of the positively charged residues of the peptide with the negatively charged residues of CaM. The nonpolar binding free energy was comparable for all three peptides, the largest contribution coming from the Trp305. About two-thirds of the buried surface area corresponds to nonpolar residues, showing that hydrophobic interactions play an important role in these CaM-peptide complexes. The simulation results agree with the experimental data in predicting a small effect of the S308D mutation on CaM interactions with DAPK2, suggesting that this mutation is not a good model for the S308 phosphorylation.     

Crystal structure of a myristoylated CAP-23/NAP-22 N-terminal domain complexed with Ca2+/calmodulin

A variety of viral and signal transduction proteins are known to be myristoylated. Although the role of myristoylation in protein-lipid interaction is well established, the involvement of myristoylation in protein-protein interactions is less well understood. CAP-23/NAP-22 is a brain-specific protein kinase C substrate protein that is involved in axon regeneration. Although the protein lacks any canonical calmodulin (CaM)-binding domain, it binds CaM with high affinity. The binding of CAP-23/NAP-22 to CaM is myristoylation dependent and the N-terminal myristoyl group is directly involved in the protein-protein interaction. Here we show the crystal structure of Ca2+-CaM bound to a myristoylated peptide corresponding to the N-terminal domain of CAP-23/NAP-22. The myristoyl moiety of the peptide goes through a hydrophobic tunnel created by the hydrophobic pockets in the N- and C-terminal domains of CaM. In addition to the myristoyl group, several amino-acid residues in the peptide are important for CaM binding. This is a novel mode of binding and is very different from the mechanism of binding in other CaM-target complexes.     

Retention of Conformational Entropy upon Calmodulin Binding to Target Peptides Is Driven by Transient Salt Bridges

Calmodulin (CaM) is a highly flexible calcium-binding protein that mediates signal transduction through an ability to differentially bind to highly variable binding sequences in target proteins. To identify how binding affects CaM motions, and its relationship to conformational entropy and target peptide sequence, we have employed fully atomistic, explicit solvent molecular dynamics simulations of unbound CaM and CaM bound to five different target peptides. The calculated CaM conformational binding entropies correlate with experimentally derived conformational entropies with a correlation coefficient R² of 0.95. Selected side-chain interactions with target peptides restrain interhelical loop motions, acting to tune the conformational entropy of the bound complex via widely distributed CaM motions. In the complex with the most conformational entropy retention (CaM in complex with the neuronal nitric oxide synthase binding sequence), Lys-148 at the C-terminus of CaM forms transient salt bridges alternating between Glu side chains in the N-domain, the central linker, and the binding target. Additional analyses of CaM structures, fluctuations, and CaM-target interactions illuminate the interplay between electrostatic, side chain, and backbone properties in the ability of CaM to recognize and discriminate against targets by tuning its conformational entropy, and suggest a need to consider conformational dynamics in optimizing binding affinities.     

Role of the N-terminal region of the skeletal muscle myosin light chain kinase target sequence in its interaction with calmodulin.

The binding of calmodulin (CaM) to four synthetic peptide analogues of the skeletal muscle myosin light chain kinase (sk-MLCK) target sequence has been studied using 1H-NMR. The 18-residue peptide WFF is anchored to CaM via the interaction of the Trp 4 side chain with the C-domain and the Phe 17 side chain with the N-domain of the protein. A peptide corresponding to the first 10 residues (WF10) does not provide the second anchoring residue and is not long enough to span both domains of CaM. 1H-NMR spectroscopy indicates that the WF10 peptide interacts specifically with the C-domain of CaM, and the chemical shifts of the bound Trp side chain are very similar in the CaM:WF10 and CaM:WFF complexes. Binding of the C-domain of CaM to the strongly basic region around Trp 4 of this MLCK sequence may be an important step in target recognition. Comparison of 1H-NMR spectra of CaM bound to WFF, a Trp 4-->Phe analogue (FFF), or a Trp 4-->Phe/Phe 17-->Trp analogue (FFW) suggests that all three peptides bind to CaM in the same orientation, i.e., with the peptide side chain in position 4 interacting with the C-domain and the side chain in position 17 interacting with the N-domain. This indicates that a Trp residue in position 4 is not an absolute requirement for binding this target sequence and that interchanging the Trp 4 and Phe 17 residues does not reverse the orientation of the bound peptide, in confirmation of the deduction from previous indirect studies using circular dichroism (Findlay WA, Martin SR, Beckingham K, Bayley PM, 1995, Biochemistry 34:2087-2094). Molecular modeling/energy minimization studies indicate that only minor local changes in the protein structure are required to accommodate binding of the bulkier Trp 17 side chain of the FFW peptide to the N-domain of CaM.     

A novel target recognition revealed by calmodulin in complex with Ca2+-calmodulin-dependent kinase kinase.

The structure of calcium-bound calmodulin (Ca2+/CaM) complexed with a 26-residue peptide, corresponding to the CaM-binding domain of rat Ca2+/CaM-dependent protein kinase kinase (CaMKK), has been determined by NMR spectroscopy. In this complex, the CaMKK peptide forms a fold comprising an alpha-helix and a hairpin-like loop whose C-terminus folds back on itself. The binding orientation of this CaMKK peptide by the two CaM domains is opposite to that observed in all other CaM-target complexes determined so far. The N- and C-terminal hydrophobic pockets of Ca2+/CaM anchor Trp 444 and Phe 459 of the CaMKK peptide, respectively. This 14-residue separation between two key hydrophobic groups is also unique among previously determined CaM complexes. The present structure represents a new and distinct class of Ca2+/CaM target recognition that may be shared by other Ca2+/CaM-stimulated proteins.     

The heterodimer calmodulin: myosin light-chain kinase as a prototype vertebrate calcium signal transduction complex.

The heterodimer complex of calmodulin (CaM) and the protein kinase catalytic subunit of myosin light chain kinase from vertebrate smooth muscle and non-muscle tissues (sm/nmMLCK) is one of the most extensively characterized CaM-regulated enzyme complexes and it has an established in vivo role in the transduction of calcium signals into biological responses. We have used a combination of approaches to the study of CaM and sm/nmMLCK in order to derive initial insight into the key features of each protein and of the CaM-MLCK heterodimeric complex that are involved in protein-protein and calcium-protein recognition and regulation of enzyme activity. On-going studies are described here that include site-specific mutagenesis, fluorescence spectroscopy, enzymology and peptide analog analysis. These and previous results indicate that: (1), both electrostatic and hydrophobic features are important in the functionally correct interactions between CaM and MLCK; (2), even the interactions between CaM and peptide analogs of the CaM binding site of MLCK are heterogeneous and non-trivial in nature; (3), amino-acid residues that have been conserved in CaM across millions of years of evolution and that are conserved in CaMs with quantitative MLCK activator activity can be mutated without any detectable effect on activity and (4), structures different from the prototypical EF-hand domain of CaM can have similar calcium-binding activity in the presence of a CaM binding structure.     

Munc13-Like skMLCK Variants Cannot Mimic the Unique Calmodulin Binding Mode of Munc13 as Evidenced by Chemical Cross-Linking and Mass Spectrometry

Among the neuronal binding partners of calmodulin (CaM) are Munc13 proteins as essential presynaptic regulators that play a key role in synaptic vesicle priming and are crucial for presynaptic short-term plasticity. Recent NMR structural investigations of a CaM/Munc13-1 peptide complex have revealed an extended structure, which contrasts the compact structures of most classical CaM/target complexes. This unusual binding mode is thought to be related to the presence of an additional hydrophobic anchor residue at position 26 of the CaM binding motif of Munc13-1, resulting in a novel 1-5-8-26 motif. Here, we addressed the question whether the 1-5-8-26 CaM binding motif is a Munc13-related feature or whether it can be induced in other CaM targets by altering the motif''s core residues. For this purpose, we chose skeletal muscle myosin light chain kinase (skMLCK) with a classical 1-5-8-14 CaM binding motif and constructed three skMLCK peptide variants mimicking Munc13-1, in which the hydrophobic anchor amino acid at position 14 was moved to position 26. Chemical cross-linking between CaM and skMLCK peptide variants combined with high-resolution mass spectrometry yielded insights into the peptides'' binding modes. This structural comparison together with complementary binding data from surface plasmon resonance experiments revealed that skMLCK variants with an artificial 1-5-8-26 motif cannot mimic CaM binding of Munc13-1. Apparently, additional features apart from the spacing of the hydrophobic anchor residues are required to define the functional 1-5-8-26 motif of Munc13-1. We conclude that Munc13 proteins display a unique CaM binding behavior to fulfill their role as efficient presynaptic calcium sensors over broad range of Ca²⁺ concentrations.     

Calmodulin’s flexibility allows for promiscuity in its interactions with target proteins and peptides

The small bilobal calcium regulatory protein calmodulin (CaM) activates numerous target enzymes in response to transient changes in intracellular calcium concentrations. Binding of calcium to the two helix-loop-helix calcium-binding motifs in each of the globular domains induces conformational changes that expose a methionine-rich hydrophobic patch on the surface of each domain of the protein, which it uses to bind to peptide sequences in its target enzymes. Although these CaM-binding domains typically have little sequence identity, the positions of several bulky hydrophobic residues are often conserved, allowing for classification of CaM-binding domains into recognition motifs, such as the 1–14 and 1–10 motifs. For calcium-independent binding of CaM, a third motif known as the IQ motif is also common. Many CaM-peptide complexes have globular conformations, where CaM’s central linker connecting the two domains unwinds, allowing the protein to wrap around a single predominantly α-helical target peptide sequence. However, novel structures have recently been reported where the conformation of CaM is highly dissimilar to these globular complexes, in some instances with less than a full compliment of bound calcium ions, as well as novel stoichiometries. Furthermore, many divergent CaM isoforms from yeast and plant species have been discovered with unique calcium-binding and enzymatic activation characteristics compared to the single CaM isoform found in mammals.     

Structurally homologous binding of plant calmodulin isoforms to the calmodulin-binding domain of vacuolar calcium-ATPase

The discovery that plants contain multiple calmodulin (CaM) isoforms having variable sequence identity to mammalian CaM has sparked a flurry of new questions regarding the intracellular role of Ca(2+) regulation in plants. To date, the majority of research in this field has focused on the differential enzymatic regulation of various mammalian CaM-dependent enzymes by the different plant CaM isoforms. However, there is comparatively little information on the structural recognition of target enzymes found exclusively in plant cells. Here we have used a variety of spectroscopic techniques, including nuclear magnetic resonance, circular dichroism, and fluorescence spectroscopy, to study the interactions of the most conserved and most divergent CaM isoforms from soybean, SCaM-1, and SCaM-4, respectively, with a synthetic peptide derived from the CaM-binding domain of cauliflower vacuolar calcium-ATPase. Despite their sequence divergence, both SCaM-1 and SCaM-4 interact with the calcium-ATPase peptide in a similar calcium-dependent, stoichiometric manner, adopting an antiparallel binding orientation with an alpha-helical peptide. The single Trp residue is bound in a solvent-inaccessible hydrophobic pocket on the C-terminal domain of either protein. Thermodynamic analysis of these interactions using isothermal titration calorimetry demonstrates that the formation of each calcium-SCaM-calcium-ATPase peptide complex is driven by favorable binding enthalpy and is very similar to the binding of mammalian CaM to the CaM-binding domains of myosin light chain kinases and calmodulin-dependent protein kinase I.     

Identification of amino acids essential for calmodulin binding and activation of smooth muscle myosin light chain kinase.

Smooth muscle myosin light chain kinase (smMLCK) is a Ca(2+)-calmodulin (CaM)-dependent enzyme that phosphorylates the 20-kDa light chains of myosin. In a previous study (Bagchi, I.C., Kemp, B.E., and Means, A.R. (1989) J. Biol. Chem. 264, 15843-15849), we expressed in bacteria a 40-kDa fragment of smMLCK that displayed Ca(2+)-CaM-regulated catalytic activity. Initial mutagenesis experiments indicated that Gly811 and Arg812 were important for CaM-dependent activation of this 40-kDa enzyme. We have now carried out site-directed mutagenesis within the CaM-binding domain (Ser787 to Leu813) of this enzyme to identify amino acids that are critical for CaM binding and activation. Our studies reveal that the individual mutation of several hydrophobic amino acid residues such as Leu813, Ile810, and Trp800 and the glycine residue Gly804 also resulted in a severe decrease in or complete loss of CaM binding and activation of smMLCK. The hydrophobic residue (Trp800) and the basic residue (Arg812), both of which are mandatory for CaM binding to smMLCK, occur in analogous positions within the CaM-binding domain of a number of CaM-regulated enzymes. We conclude from these results that CaM binding by smMLCK is determined by an interplay of specific hydrophobic and electrostatic interactions which appear to be conserved among various target enzymes of CaM.     

Calcium-dependent and -independent binding of soybean calmodulin isoforms to the calmodulin binding domain of tobacco MAPK phosphatase-1

The recent finding of an interaction between calmodulin (CaM) and the tobacco mitogen-activated protein kinase phosphatase-1 (NtMKP1) establishes an important connection between Ca(2+) signaling and the MAPK cascade, two of the most important signaling pathways in plant cells. Here we have used different biophysical techniques, including fluorescence and NMR spectroscopy as well as microcalorimetry, to characterize the binding of soybean CaM isoforms, SCaM-1 and -4, to synthetic peptides derived from the CaM binding domain of NtMKP1. We find that the actual CaM binding region is shorter than what had previously been suggested. Moreover, the peptide binds to the SCaM C-terminal domain even in the absence of free Ca(2+) with the single Trp residue of the NtMKP1 peptides buried in a solvent-inaccessible hydrophobic region. In the presence of Ca(2+), the peptides bind first to the C-terminal lobe of the SCaMs with a nanomolar affinity, and at higher peptide concentrations, a second peptide binds to the N-terminal domain with lower affinity. Thermodynamic analysis demonstrates that the formation of the peptide-bound complex with the Ca(2+)-loaded SCaMs is driven by favorable binding enthalpy due to a combination of hydrophobic and electrostatic interactions. Experiments with CaM proteolytic fragments showed that the two domains bind the peptide in an independent manner. To our knowledge, this is the first report providing direct evidence for sequential binding of two identical peptides of a target protein to CaM. Discussion of the potential biological role of this interaction motif is also provided.     

Solution NMR Structure of the Ca2+-bound N-terminal Domain of CaBP7: A REGULATOR OF GOLGI TRAFFICKING*

Calcium-binding protein 7 (CaBP7) is a member of the calmodulin (CaM) superfamily that harbors two high affinity EF-hand motifs and a C-terminal transmembrane domain. CaBP7 has been previously shown to interact with and modulate phosphatidylinositol 4-kinase III-β (PI4KIIIβ) activity in in vitro assays and affects vesicle transport in neurons when overexpressed. Here we show that the N-terminal domain (NTD) of CaBP7 is sufficient to mediate the interaction of CaBP7 with PI4KIIIβ. CaBP7 NTD encompasses the two high affinity Ca²⁺ binding sites, and structural characterization through multiangle light scattering, circular dichroism, and NMR reveals unique properties for this domain. CaBP7 NTD binds specifically to Ca²⁺ but not Mg²⁺ and undergoes significant conformational changes in both secondary and tertiary structure upon Ca²⁺ binding. The Ca²⁺-bound form of CaBP7 NTD is monomeric and exhibits an open conformation similar to that of CaM. Ca²⁺-bound CaBP7 NTD has a solvent-exposed hydrophobic surface that is more expansive than observed in CaM or CaBP1. Within this hydrophobic pocket, there is a significant reduction in the number of methionine residues that are conserved in CaM and CaBP1 and shown to be important for target recognition. In CaBP7 NTD, these residues are replaced with isoleucine and leucine residues with branched side chains that are intrinsically more rigid than the flexible methionine side chain. We propose that these differences in surface hydrophobicity, charge, and methionine content may be important in determining highly specific interactions of CaBP7 with target proteins, such as PI4KIIIβ.     

Cooperativity of hydrophobic anchor interactions: evidence for epitope selection by MHC class II as a folding process

Peptide binding to MHC class II (MHCII) molecules is stabilized by hydrophobic anchoring and hydrogen bond formation. We view peptide binding as a process in which the peptide folds into the binding groove and to some extent the groove folds around the peptide. Our previous observation of cooperativity when analyzing binding properties of peptides modified at side chains with medium to high solvent accessibility is compatible with such a view. However, a large component of peptide binding is mediated by residues with strong hydrophobic interactions that bind to their respective pockets. If these reflect initial nucleation events they may be upstream of the folding process and not show cooperativity. To test whether the folding hypothesis extends to these anchor interactions, we measured dissociation and affinity to HLA-DR1 of an influenza hemagglutinin-derived peptide with multiple substitutions at major anchor residues. Our results show both negative and positive cooperative effects between hydrophobic pocket interactions. Cooperativity was also observed between hydrophobic pockets and positions with intermediate solvent accessibility, indicating that hydrophobic interactions participate in the overall folding process. These findings point out that predicting the binding potential of epitopes cannot assume additive and independent contributions of the interactions between major MHCII pockets and corresponding peptide side chains.     

Tryptophan fluorescence of calmodulin binding domain peptides interacting with calmodulin containing unnatural methionine analogues

The interactions between the abundant methionine residues of the calcium regulatory protein calmodulin (CaM) and several of its binding targets were probed using fluorescence spectroscopy. Tryptophan steady-state fluorescence from peptides encompassing the CaM-binding domains of the target proteins myosin light chain kinase (MLCK), cyclic nucleotide phosphodiesterase (PDE) and caldesmon site A and B (CaD A, CaD B), and the model peptide melittin showed Ca(2+)-dependent blue-shifts in their maximum emission wavelength when complexed with wild-type CaM. Blue-shifts were also observed for complexes in which the CaM methionine residues were replaced by selenomethionine, norleucine and ethionine, and when a quadruple methionine to leucine C-terminal mutant of CaM was studied. Quenching of the tryptophan fluorescence intensity was observed with selenomethionine, but not with norleucine or ethionine substituted protein. Fluorescence quenching studies with added potassium iodide (KI) demonstrate that the non-native proteins limit the solvent accessibility of the Trp in the MLCK peptide to levels close to that of the wild-type CaM-MLCK interaction. Our results show that the methionine residues from CaM are highly sensitive to the target peptide in question, confirming the importance of their role in binding interactions. In addition, we provide evidence that the nature of binding in the CaM-CaD B complex is unique compared with the other complexes studied, as the Trp residue of this peptide remains partially solvent exposed upon binding to CaM.     

First observation of solution structures of bradykinin-penta-O-galloyl-D-glucopyranose complexes as determined by NMR and simulated annealing

Polyphenols (tannins) are known for their high propensity to precipitate proteins. They bind most strongly to proteins with a high proline content. Understanding the mechanism of this association is of prime interest because this interaction might induce protein conformational changes that may modify their biological activity. To investigate the interaction, an NMR study was carried out on the binding of a representative polyphenol, penta-O-galloyl-D-glucopyranose, to a nonapeptide hormone, bradykinin (BDK), where proline accounts for 30% of residues. Series of 1D and 2D-NMR experiments were performed. For the first time, a three-dimensional structure of complexes was determined using 2D-NMR experiments and molecular modeling. These structure calculations are a potent tool to understand how the association arises. They clearly show that the interaction is a complex phenomenon where several parameters are involved. The PGG/BDK complexes are formed by multiple weak interactions between peptide side chains and galloyl rings. Proline and arginine are good anchoring points and the glycine gives a certain flexibility in the peptide backbone that allows the polyphenol to approach and interact. Therefore, it is not only the hydrophobic stackings between galloyl rings and proline and hydrogen bonding involving arginine and aromatic rings which are important. The residue sequence and the side chain steric bulk also intervene.     

First observation of solution structures of bradykinin–penta-O-galloyl-d-glucopyranose complexes as determined by NMR and simulated annealing

Polyphenols (tannins) are known for their high propensity to precipitate proteins. They bind most strongly to proteins with a high proline content. Understanding the mechanism of this association is of prime interest because this interaction might induce protein conformational changes that may modify their biological activity. To investigate the interaction, an NMR study was carried out on the binding of a representative polyphenol, penta-O-galloyl-d-glucopyranose, to a nonapeptide hormone, bradykinin (BDK), where proline accounts for 30% of residues. Series of 1D and 2D-NMR experiments were performed. For the first time, a three-dimensional structure of complexes was determined using 2D-NMR experiments and molecular modeling. These structure calculations are a potent tool to understand how the association arises. They clearly show that the interaction is a complex phenomenon where several parameters are involved. The PGG/BDK complexes are formed by multiple weak interactions between peptide side chains and galloyl rings. Proline and arginine are good anchoring points and the glycine gives a certain flexibility in the peptide backbone that allows the polyphenol to approach and interact. Therefore, it is not only the hydrophobic stackings between galloyl rings and proline and hydrogen bonding involving arginine and aromatic rings which are important. The residue sequence and the side chain steric bulk also intervene.     

Cooperativity during the formation of peptide/MHC class II complexes

To generate an effective immune response, class II major histocompatibility complex molecules (MHCII) must present a diverse array of peptide ligands for recognition by T lymphocytes. Peptide/MHCII complexes are stabilized by hydrophobic anchoring of peptide side chains to pockets in the MHCII protein and the formation of hydrogen bonds to the peptide backbone. Many current models of peptide/MHCII association assume an additive and independent contribution of the interactions between major MHCII pockets and corresponding side chains in the peptide. However, significant conformational rearrangements occur in both the peptide and MHCII during binding. Therefore, we hypothesize that peptide binding to MHCII could be viewed as a folding process in which both molecules cooperate to produce the final conformation. To directly test this hypothesis, we adapt a serial mutagenesis strategy to study cooperativity in the interaction of the human MHCII HLA-DR1 and a peptide derived from influenza hemagglutinin. Substitutions in either the peptide or HLA-DR1 that are predicted to interfere with hydrogen bond formation show cooperative effects on complex stability and affinity. Substitution of a peptide side chain that provides a hydrophobic contact also contributes to the cooperative effect, suggesting a role for all energetic sources in the folding process. We propose that cooperativity throughout the peptide-binding groove reflects the folding of segments of the MHCII molecule into helices around the peptide with a concomitant folding of the peptide into a polyproline helix. The implications of cooperativity for peptide/MHCII structure and epitope selection are discussed.     

Molecular dynamics study of peptide-bilayer adsorption

Two 6-ns simulations of the somatostatin analog sandostatin and a 1-palmityl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer are presented. In the first simulation, the peptide was placed in a region of bulk water density and allowed to spontaneously move toward and bind to the bilayer surface. An attractive force between the peptide and bilayer drove the binding process, which was opposed by a significant frictional force caused by the solvent (water). During the approach of the peptide toward the bilayer the area of the interacting surface between the species was inversely proportional to the distance between them, supporting the application of such a relationship in continuum calculations of peptide-bilayer binding free energies. In the second simulation, the N-terminus of the surface-bound peptide was deprotonated. Consistent with experiment, this strengthened interactions between the peptide and the bilayer. Details of both peptide-bilayer complexes, including the orientation, percent buried surface area, and orientation of the lipid headgroups are in good agreement with those obtained from experiment. The location of the different side chains in the bilayer is in direct correlation with an interfacial hydrophobicity scale developed using model peptides. The aromatic side chains of the Phe and Trp residues all lie flat with respect to the bilayer surface in both complexes. Changes in lipid and water ordering due to peptide binding suggest a possible domination of lipophobic over hydrophobic effects, as proposed by other workers. Where appropriate, peptide and lipid properties in the bound states are compared with separate simulations of sandostatin and the bilayer in water, respectively, so as to monitor the response of the system to the binding process.