Lipid-protein nanodiscs: possible application in high-resolution NMR investigations of membrane proteins and membrane-active peptides

High-resolution NMR is shown to be applicable for investigation of membrane proteins and membrane-active peptides embedded into lipid-protein nanodiscs (LPNs). (15)N-Labeled K+-channel from Streptomyces lividans (KcsA) and the antibiotic antiamoebin I from Emericellopsis minima (Aam-I) were embedded in LPNs of different lipid composition. Formation of stable complexes undergoing isotropic motion in solution was confirmed by size-exclusion chromatography and (31)P-NMR spectroscopy. The 2D 1H-(15)N-correlation spectra were recorded for KcsA in the complex with LPN containing DMPC and for Aam-I in LPNs based on DOPG, DLPC, DMPC, and POPC. The spectra recorded were compared with those in detergent-containing micelles and small bicelles commonly used in high-resolution NMR spectroscopy of membrane proteins. The spectra recorded in LPN environments demonstrated similar signal dispersion but significantly increased (1)H(N) line width. The spectra of Aam-I embedded in LPNs containing phosphatidylcholine showed significant selective line broadening, thus suggesting exchange process(es) between several membrane-bound states of the peptide. (15)N relaxation rates were measured to obtain the effective rotational correlation time of the Aam-I molecule. The obtained value (approximately 40 nsec at 45 degrees C) is indicative of additional peptide motions within the Aam-I/LPN complex.

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.

Cell-free Production of the Extracellular Domain of the Nicotinic Acetylcholine Receptor

Collapse

Adaptation of a membrane-active peptide to heterogeneous environment. II. The role of mosaic nature of the membrane surface

In the first article of this series we demonstrated the importance of specific intrapeptide interactions and peptide-lipid contacts for the membrane binding of penetratin (pAntp). Here in focus was detailed characterization of spatial hydrophobic/hydrophilic properties of the bilayer surface and their influence on the binding mode of pAntp. From the hydrophobicity point of view, the solvent-accessible surfaces of lipid bilayers possess a distinctly "mosaic" character. This correlates well with the occurrence of dynamic clusters of hydrophobic surface area formed by hydrocarbon tails of phospholipids exposed on the interface. Such mosaic patterns are specific for lipid bilayers of particular composition. In an anionic membrane, they determine initial stages of pAntp adsorption, which strongly depends on the "complementarity" between polarity properties of the peptide and its local interfacial environment. If high complementarity is established, then pAntp penetrates deeply into the membrane without significant destabilization of its initial secondary structure. Alternatively, partial unfolding of pAntp takes place in order to compensate unfavorable peptide-membrane interactions upon embedding. Such effects explain complicated behavior of membrane-active peptides, especially if the target membrane surface is of distinctly mosaic nature, depending on the microscopic properties of the water-lipid interface, pAntp is capable of adopting different pathways to exercise its biological activity.

Spatial structure and pH-dependent conformational diversity of dimeric transmembrane domain of the receptor tyrosine kinase EphA1

Eph receptors are found in a wide variety of cells in developing and mature tissues and represent the largest family of receptor tyrosine kinases, regulating cell shape, movements, and attachment. The receptor tyrosine kinases conduct biochemical signals across plasma membrane via lateral dimerization in which their transmembrane domains play an important role. Structural-dynamic properties of the homodimeric transmembrane domain of the EphA1 receptor were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics in explicit lipid bilayer. EphA1 transmembrane segments associate in a right-handed parallel alpha-helical bundle, region (544-569)(2), through the N-terminal glycine zipper motif A(550)X(3)G(554)X(3)G(558). Under acidic conditions, the N terminus of the transmembrane helix is stabilized by an N-capping box formed by the uncharged carboxyl group of Glu(547), whereas its deprotonation results in a rearrangement of hydrogen bonds, fractional unfolding of the helix, and a realignment of the helix-helix packing with appearance of additional minor dimer conformation utilizing seemingly the C-terminal GG4-like dimerization motif A(560)X(3)G(564). This can be interpreted as the ability of the EphA1 receptor to adjust its response to ligand binding according to extracellular pH. The dependence of the pK(a) value of Glu(547) and the dimer conformational equilibrium on the lipid head charge suggests that both local environment and membrane surface potential can modulate dimerization and activation of the receptor. This makes the EphA1 receptor unique among the Eph family, implying its possible physiological role as an "extracellular pH sensor," and can have relevant physiological implications.

Molecular insight into mechanism of antimicrobial action of the beta-hairpin peptide arenicin: specific oligomerization in detergent micelles

Arenicins are 21-residue cationic antimicrobial peptides isolated from marine polychaeta Arenicola marina. The peptides exhibit potent broad-spectrum antimicrobial activity. In water solution arenicin-2 adopts a beta-hairpin conformation, stabilized by one disulfide and nine hydrogen bonds. To determine the propensity for the peptide oligomerization in membrane mimetic systems, the recombinant arenicin-2 was overexpressed as a fused form in Escherichia coli. The arenicin-2 oligomerization and intermolecular packing in membrane mimicking environment were investigated using high-resolution NMR spectroscopy. The present studies show that arenicin-2 preserves a beta-hairpin structure and forms asymmetric dimers upon incorporation into the dodecylphosphocholine micelle. Two monomers of arenicin-2 are aligned parallel to each other by the N-terminal strands of the beta-hairpin (CN upward arrow upward arrowNC type of association). Polyacrylamide gel electrophoresis analysis indicated that in environment of anionic SDS micelles the arenicin-2 might undergo further oligomerization and form tetramers. Our results afford further molecular insight into possible mechanism of antimicrobial action of arenicins.

Three-dimensional structure/hydrophobicity of latarcins specifies their mode of membrane activity

Latarcins, linear peptides from the Lachesana tarabaevi spider venom, exhibit a broad-spectrum antimicrobial activity, likely acting on the bacterial cytoplasmic membrane. We study their spatial structures and interaction with model membranes by a combination of experimental and theoretical methods to reveal the structure-activity relationship. In this work, a 26 amino acid peptide, Ltc1, was investigated. Its spatial structure in detergent micelles was determined by (1)H nuclear magnetic resonance (NMR) and refined by Monte Carlo simulations in an implicit water-octanol slab. The Ltc1 molecule was found to form a straight uninterrupted amphiphilic helix comprising 8-23 residues. A dye-leakage fluorescent assay and (31)P NMR spectroscopy established that the peptide does not induce the release of fluorescent marker nor deteriorate the bilayer structure of the membranes. The voltage-clamp technique showed that Ltc1 induces the current fluctuations through planar membranes when the sign of the applied potential coincides with the one across the bacterial inner membrane. This implies that Ltc1 acts on the membranes via a specific mechanism, which is different from the carpet mode demonstrated by another latarcin, Ltc2a, featuring a helix-hinge-helix structure with a hydrophobicity gradient along the peptide chain. In contrast, the hydrophobic surface of the Ltc1 helix is narrow-shaped and extends with no gradient along the axis. We have also disclosed a number of peptides, structurally homologous to Ltc1 and exhibiting similar membrane activity. This indicates that the hydrophobic pattern of the Ltc1 helix and related antimicrobial peptides specifies their activity mechanism. The latter assumes the formation of variable-sized lesions, which depend upon the potential across the membrane.

Divalent cation coordination and mode of membrane interaction in cyclotides: NMR spatial structure of ternary complex Kalata B7/Mn2+/DPC micelle

The cyclotides are the family of hydrophobic bioactive plant peptides, characterized by a circular protein backbone and three knot forming disulfide bonds. It is believed that membrane activity of the cyclotides underlines their antimicrobial, cytotoxic and hemolytic properties, but the specific interactions with divalent cations can be also involved. To assess the mode of membrane interaction and divalent cation coordination in cyclotides, the spatial structure of the Möbius cyclotide Kalata B7 from the African perennial plant Oldenlandia affinis was determined in the presence of anisotropic membrane mimetic (dodecylphosphocholine micelles). The model of peptide/cation/micelle complex was built using 5-doxylstearate and Mn2+ relaxation probes. Results show that the peptide binds to the micelle surface with relatively high affinity by two hydrophobic loops (loop 2 - Thr6-Leu7 and loop 5 - Trp19-Ile21). The partially hydrated divalent cation is coordinated by charged side-chain of Glu3, aromatic side chain of Tyr11 and free carbonyls of Thr4 and Thr9, and is located in direct contact with the polar head-groups of detergent. The comparison with data about other cyclotides indicates that divalent cation coordination is the invariant property of all cyclotides, but the mode of peptide/membrane interactions is varied. Probably, the specific cation/peptide interactions play a major, but yet not known, role in the biological activity of the cyclotides.

Role of lipid charge in organization of water/lipid bilayer interface: insights via computer simulations

Anionic unsaturated lipid bilayers represent suitable model systems that mimic real cell membranes: they are fluid and possess a negative surface charge. Understanding of detailed molecular organization of water-lipid interfaces in such systems may provide an important insight into the mechanisms of proteins' binding to membranes. Molecular dynamics (MD) of full-atom hydrated lipid bilayers is one of the most powerful tools to address this problem in silico. Unfortunately, wide application of computational methods for such systems is limited by serious technical problems. They are mainly related to correct treatment of long-range electrostatic effects. In this study a physically reliable model of an anionic unsaturated bilayer of 1,2-dioleoyl-sn-glycero-3-phosphoserine (DOPS) was elaborated and subjected to long-term MD simulations. Electrostatic interactions were treated with two different algorithms: spherical cutoff function and particle-mesh Ewald summation (PME). To understand the role of lipid charge in the system behavior, similar calculations were also carried out for zwitterionic bilayer composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). It was shown that, for the charged DOPS bilayer, the PME protocol performs much better than the cutoff scheme. In the last case a number of artifacts in the structural organization of the bilayer were observed. All of them were attributed to inadequate treatment of electrostatic interactions of lipid headgroups with counterions. Electrostatic properties, along with structural and dynamic parameters, of both lipid bilayers were investigated. Comparative analysis of the MD data reveals that the water-lipid interface of the DOPC bilayer is looser than that for DOPS. This makes possible deeper penetration of water molecules inside the zwitterionic (DOPC) bilayer, where they strongly interact with carbonyls of lipids. This can lead to thickening of the membrane interface in zwitterionic as compared to negatively charged bilayers.

The Membrane-proximal Fusion Domain of HIV-1 GP41 Reveals Sequence-specific and Fine-tuning Mechanism of Membrane Binding

The membrane interface-partitioning region preceding the transmembrane anchor of the human immunodeficiency virus type 1 (HIV-1) gp41 envelope protein is one of the sites responsible for virus binding to its host cell membrane and subsequent fusion events. Here, we used molecular modeling techniques to assess membrane interactions, structure, and hydrophobic properties of the fusion-active peptide representing this region, several of its homologs from different HIV-1 strains, as well as a peptide - defective gp41 phenotype - unable to mediate cell-cell fusion and virus entry. It is shown that the wild-type peptides bind to the water-membrane interface in alpha-helical conformation, while the mutant adopts partly destabilized helix-break-helix structure on the membrane surface. The wild-type peptides reveal specific "tilted oblique-oriented" pattern of hydrophobicity on their surfaces - the property specific for fusion regions of other viruses. Fusion peptides penetrate into the membrane with their N-termini and reveal "fine-tuning" interactions with membrane and water environments: the shift of this balance (e.g., due to point mutations) may dramatically change the mode of membrane binding, and therefore, may cause loss of fusion activity. The modeling results agree well with experimental data and provide a strategy to delineate fusogenic regions in amino acid sequences of viral proteins.

Recombinant expression, synthesis, purification, and solution structure of arenicin

Arenicins are 21-residue cationic antimicrobial peptides, isolated from marine polychaeta Arenicola marina. In order to determine a high-resolution three-dimensional structure of arenicin-2, the recombinant peptide was overexpressed as a fused form in Escherichia coli. Both arenicin isoforms were synthesized using the Fmoc-based solid-phase strategy. Recombinant and synthetic arenicins were purified, and their antimicrobial and spectroscopic properties were analyzed. NMR investigation shows that in water solution arenicin-2 displays a prolonged beta-hairpin, formed by two antiparallel beta-strands and stabilized by one disulfide and nine hydrogen bonds. A significant right-handed twist in the beta-sheet is deprived the peptide surface of amphipathicity. CD spectroscopic analysis indicates that arenicin-2 binds to the SDS and DPC micelles, and conformation of the peptide is significantly changed upon binding. Arenicin strongly binds to anionic lipid (POPE/POPG) vesicles in contrast with zwitterionic (POPC) ones. These results suggest that arenicins are membrane active peptides and point to possible mechanism of their selectivity toward bacterial cells.

Antiamoebin I in Methanol Solution: Rapid Exchange between Right-Handed and Left-Handed310-Helical Conformations

Antiamoebin I (Aam-I) is a membrane-active peptaibol antibiotic isolated from fungal species belonging to the genera Cephalosporium, Emericellopsis, Gliocladium, and Stilbella. Antiamoebin I has the amino acid sequence: Ac-Phe(1)-Aib-Aib-Aib-Iva-Gly-Leu-Aib(8)-Aib-Hyp-Gln-Iva-Hyp-Aib-Pro-Phl(16). By using the uniformly (13)C,(15)N-labeled sample of Aam-I, the set of conformationally dependent J couplings and (3h)J(NC) couplings through H-bonds were measured. Analysis of these data along with the data on magnetic nonequivalence of the (13)C(beta) nuclei (Deltadelta((13)C(beta))) in Aib and Iva residues allowed us to draw the univocal conclusion that the N-terminal part (Phe(1)-Gly(6)) of Aam-I in MeOH solution is in fast exchange between the right-handed and left-handed 3(10)-helical conformations, with an approximately equal population of both states. An additional conformational exchange process was found at the Aib(8) residue. The (15)N-NMR-relaxation and CD-spectroscopy measurements confirmed these findings. Molecular modeling and Monte Carlo simulations revealed that both exchange processes are correlated and coupled with significant hinge-bending motions around the Aib(8) residue. Our results explain relatively low activity of Aam-I with respect to other 15-amino acid residue peptaibols (for example, zervamicin) in functional and biological tests. The high dynamic 'propensity' possibly prevents both initial binding of the antiamoebin to the membrane and subsequent formation of stable ionic channels according to the barrel-stave mechanism.

Docking of ATP to Ca-ATPase: considering protein domain motions

Most standard molecular docking algorithms take into account only ligand flexibility, while numerous studies demonstrate that receptor flexibility may be also important. While some efficient methods have been proposed to take into account local flexibility of protein side chains, the influence of large-scale domain motions on the docking results still represents a challenge for computational methods. In this work we compared the results of ATP docking to different models of Ca-ATPase: crystallographic apo- and holo-forms of the enzyme as well as "flexible" target models generated via molecular dynamics (MD) simulations in water. MD simulations were performed for two different apo-forms and one holo-form of Ca2+-ATPase and reveal large-scale domain motions of type "closure", which is consistent with experimental structures. Docking to a set of MD-conformers yielded correct solutions with ATP bound in both domains regardless of the starting Ca2+-ATPase structure. Also, special attention was paid to proper ranking of docking solutions and some particular features of different scoring functions and their applicability for the model of "flexible" receptor. Particularly, the results of docking ATP were ranked by a scoring criterion specially designed to estimate ATP-protein interactions. This criterion includes stacking and hydrophobic interactions characteristic of ATP-protein complexes. The performance of this ligand-specific scoring function was considerably better than that of a standard scoring function used in the docking algorithm.

Unique dimeric structure of BNip3 transmembrane domain suggests membrane permeabilization as a cell death trigger

BNip3 is a prominent representative of apoptotic Bcl-2 proteins with rather unique properties initiating an atypical programmed cell death pathway resembling both necrosis and apoptosis. Many Bcl-2 family proteins modulate the permeability state of the outer mitochondrial membrane by forming homo- and hetero-oligomers. The structure and dynamics of the homodimeric transmembrane domain of BNip3 were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics energy relaxation in an explicit lipid bilayer. The right-handed parallel helix-helix structure of the domain with a hydrogen bond-rich His-Ser node in the middle of the membrane, accessibility of the node for water, and continuous hydrophilic track across the membrane suggest that the domain can provide an ion-conducting pathway through the membrane. Incorporation of the BNip3 transmembrane domain into an artificial lipid bilayer resulted in pH-dependent conductivity increase. A possible biological implication of the findings in relation to triggering necrosis-like cell death by BNip3 is discussed.

Propagation of Dynamic Changes in Barnase Upon Binding of Barstar: An NMR and Computational Study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.

Molecular lipophilicity in protein modeling and drug design

Hydrophobic interactions play a key role in the folding and maintenance of the 3-dimensional structure of proteins, as well as in the binding of ligands (e.g. drugs) to protein targets. Therefore, quantitative assessment of spatial hydrophobic (lipophilic) properties of these molecules is indispensable for the development of efficient computational methods in drug design. One possible solution to the problem lies in application of a concept of the 3-dimensional molecular hydrophobicity potential (MHP). The formalism of MHP utilizes a set of atomic physicochemical parameters evaluated from octanol-water partition coefficients (log P) of numerous chemical compounds. It permits detailed assessment of the hydrophobic and/or hydrophilic properties of various parts of molecules and may be useful in analysis of protein-protein and protein-ligand interactions. This review surveys recent applications of MHP-based techniques to a number of biologically relevant tasks. Among them are: (i) Detailed assessment of hydrophobic/hydrophilic organization of proteins; (ii) Application of this data to the modeling of structure, dynamics, and function of globular and membrane proteins, membrane-active peptides, etc. (iii) Employment of the MHP-based criteria in docking simulations for ligands binding to receptors. It is demonstrated that the application of the MHP-based techniques in combination with other molecular modeling tools (e.g. Monte Carlo and molecular dynamics simulations, docking, etc.) permits significant improvement to the standard computational approaches, provides additional important insights into the intimate molecular mechanisms driving protein assembling in water and in biological membranes, and helps in the computer-aided drug discovery process.

Method To Assess Packing Quality of Transmembrane α-Helices in Proteins. 1. Parametrization Using Structural Data

Integral membrane proteins (MPs) are pharmaceutical targets of exceptional importance. Modern methods of three-dimensional protein structure determination often fail to supply the fast growing field of structure-based drug design with the requested MPs' structures. That is why computational modeling techniques gain a special importance for these objects. Among the principal difficulties limiting application of these methods is the low quality of the MPs' models built in silico. In this series of two papers we present a computational approach to the assessment of the packing "quality" of transmembrane (TM) alpha-helical domains in proteins. The method is based on the concept of protein environment classes, whereby each amino acid residue is described in terms of its environment polarity and accessibility to the membrane. In the first paper we analyze a nonredundant set of 26 TM alpha-helical domains and compute the residues' propensities to five predefined classes of membrane-protein environments. Here we evaluate the proposed approach only by various test sets, cross-validation protocols and ability of the method to delimit the crystal structure of visual rhodopsin, and a number of its erroneous theoretical models. More advanced validation of the method is given in the second article of this series. We assume that the developed "membrane score" method will be helpful in optimizing computer models of TM domains of MPs, especially G-protein coupled receptors.

Method To Assess Packing Quality of Transmembrane α-Helices in Proteins. 2. Validation by “Correct vs Misleading” Test

We describe a set of tests designed to check the ability of the new "membrane score" method (see the first paper of this series) to assess the packing quality of transmembrane (TM) alpha-helical domains in proteins. The following issues were addressed: (1) Whether there is a relation between the score (S(mem)) of a model and its closeness to the "nativelike" conformation? (2) Is it possible to recognize a correct model among misfolded and erroneous ones? (3) To what extent the score of a homology-built model is sensitive to errors in sequence alignment? To answer the first question, two test cases were considered: (i) Several models of bovine aquaporin-1 (target protein) were built on the structural templates provided by its homologs with known X-ray structure. (ii) Side chains in the spatial models of visual rhodopsin and cytochrome c oxidase were rebuilt based on the backbone scaffolds taken from their crystal structures, and the resulting models were iteratively fitted into the full-atom X-ray conformations. It was shown that the higher the S(mem) value of a model is, the lower its root-mean-square deviation is from the "correct" (crystal) structure of a target. Furthermore, the "membrane score" method successfully identifies the rhodopsin crystal structure in an ensemble of "rotamer-type" decoys, thus providing the way to optimize mutual orientations of alpha-helices in models of TM domains. Finally, being applied to a set of homology models of rhodopsin built on its crystal structure with systematically shifted alignment, the approach demonstrates a prominent ability to detect alignment errors. We therefore assume that the "membrane score" method will be helpful in optimization of in silico models of TM domains in proteins, especially those in GPCRs.

Complementarity of hydrophobic properties in ATP-protein binding: a new criterion to rank docking solutions

ATP is an important substrate of numerous biochemical reactions in living cells. Molecular recognition of this ligand by proteins is very important for understanding enzymatic mechanisms. Considerable insight into the problem may be gained via molecular docking simulations. At the same time, standard docking protocols are often insufficient to predict correct conformations for protein-ATP complexes. Thus, in most cases the native-like solutions can be found among the docking poses, but current scoring functions have only limited ability to discriminate them from false positives. To improve the selection of correct docking solutions obtained with the GOLD software, we developed a new ranking criterion specific for ATP-protein binding. The method is based on detailed analysis of the intermolecular interactions in 40 high-resolution 3D structures of ATP-protein complexes (the training set). We found that the most important factors governing this recognition are hydrogen-bonding, stacking between adenine and aromatic protein residues, and hydrophobic contacts between adenine and protein residues. To address the latter, we applied the formalism of 3D molecular hydrophobicity potential. The results obtained were used to construct an ATP-oriented scoring criterion as a linear combination of the terms describing these intermolecular interactions. The criterion was then validated using the test set of 10 additional ATP-protein complexes. As compared with the standard scoring functions, the new ranking criterion significantly improved the selection of correct docking solutions in both sets and allowed considerable enrichment at the top of the list containing docking poses with correct solutions.

A novel method for packing quality assessment of transmembrane alpha-helical domains in proteins

Here we present a novel method for assessment of packing quality for transmembrane (TM) domains of alpha-helical membrane proteins (MPs), based on analysis of available high-resolution experimental structures of MPs. The presented concept of protein-membrane environment classes permits quantitative description of packing characteristics in terms of membrane accessibility and polarity of the nearest protein groups. We demonstrate that the method allows identification of native-like conformations among the large set of theoretical MP models. The developed "membrane scoring function" will be of use for optimization of TM domain packing in theoretical models of MPs, first of all G-protein coupled receptors.

Conformation and mode of membrane interaction in cyclotides. Spatial structure of kalata B1 bound to a dodecylphosphocholine micelle

Cyclotides are a family of bioactive plant peptides that are characterized by a circular protein backbone and three conserved tightly packed disulfide bonds. The antimicrobial and hemolytic properties of cyclotides, along with the relative hydrophobicity of the peptides, point to the biological membrane as a target for cyclotides. To assess the membrane-induced conformation and orientation of cyclotides, the interaction of the Möbius cyclotide, kalata B1, from the African perennial plant Oldenlandia affinis, with dodecylphosphocholine micelles was studied using NMR spectroscopy. Under conditions where the cyclotide formed a well-defined complex with micelles, the spatial structure of kalata B1 was calculated from NOE and J couplings data, and the model for the peptide-micelle complex was built using 5- and 16-doxylstearate relaxation probes. The binding of divalent cations to the peptide-micelle complex was quantified by Mn2+ titration. The results show that the peptide binds to the micelle surface, with relatively high affinity, via two hydrophobic loops (loop 5, Trp19-Val21; and loop6, Leu27-Val29). The charged residues (Glu3 and Arg24), along with the cation-binding site (near Glu3) are segregated on the other side of the molecule and in contact with polar head groups of detergent. The spatial structure of kalata B1 is only slightly changed during incorporation into micelles and represents a distorted triple-stranded beta-sheet cross-linked by a cystine knot. Detailed structural analysis and comparison with other knottins revealed structural conservation of the two-disulfide motif in cyclic and acyclic peptides. The results thus obtained provide the first model for interaction of cyclotides with membranes and permit consideration of the cyclotides as membrane-active cationic antimicrobial peptides.