Archaeal cyclopentane fragment in a surfactant's hydrophobic tail decreases the Krafft point.

Archaea are prokaryotic microorganisms famous for their ability to adapt to extreme environments, including low and high temperatures. Archaeal lipids often are macrocycles with two polar heads and a hydrophobic core that contains methyl groups and in-line cycles. Here we present the design of novel general-purpose surfactants that have inherited features of archaeal lipids. These are C12 and C14 carboxylic acids containing in-line cyclopentanes. The cyclopentanes disturb the chain packing, which results in remarkable expansion of the operational range of the surfactant into the low-temperature region. We report synthesis and properties of these novel archaea-like surfactants and details of their chain packing derived from thermodynamics model predictions, molecular dynamics simulations, and experimental data on CMC and Krafft points.

Cyclopentane rings in hydrophobic chains of a phospholipid enhance the bilayer stability to electric breakdown.

Archaeal lipids ensure unprecedented stability of archaea membranes in extreme environments. Here, we incorporate a characteristic structural feature of an archaeal lipid, the cyclopentane ring, into hydrocarbon chains of a short-chain (C12) phosphatidylcholine to explore whether the insertion would allow such a lipid (1,2-di-(3-(3-hexylcyclopentyl)-propanoate)-sn-glycero-3-phosphatidylcholine, diC12cp-PC) to form stable bilayers at room temperature. According to fluorescence-based assays, in water diC12cp-PC formed liquid-crystalline bilayers at room temperature. Liposomes produced from diC12cp-PC retained calcein for over a week when stored at +4 °C. diC12cp-PC could also form model bilayer lipid membranes that were by an order of magnitude more stable to electrical breakdown than egg PC membranes. Molecular dynamics simulation showed that the cyclopentane fragment fixes five carbon atoms (or four C-C bonds), which is compensated by the higher mobility of the rest of the chain. This was found to be the reason for the remarkable stability of the diC12cp-PC bilayer: restricted conformational mobility of a chain segment increases the membrane bending modulus (compared to a normal hydrocarbon chain of the same length). Here, higher stiffness practically does not affect the line tension of a membrane pore edge. Rather it makes it more difficult for diC12cp-PC to rearrange in order to line the edge of a hydrophilic pore; therefore, fewer pores are formed.

Structural basis of the signal transduction via transmembrane domain of the human growth hormone receptor.

BACKGROUND: Prior studies of the human growth hormone receptor (GHR) revealed a distinct role of spatial rearrangements of its dimeric transmembrane domain in signal transduction across membrane. Detailed structural information obtained in the present study allowed elucidating the bases of such rearrangement and provided novel insights into receptor functioning. METHODS: We investigated the dimerization of recombinant TMD fragment GHR254-294 by means of high-resolution NMR in DPC micelles and molecular dynamics in explicit POPC membrane. RESULTS: We resolved two distinct dimeric structures of GHR TMD coexisting in membrane-mimicking micellar environment and providing left- and right-handed helix-helix association via different dimerization motifs. Based on the available mutagenesis data, the conformations correspond to the dormant and active receptor states and are distinguished by cis-trans isomerization of Phe-Pro266 bond in the transmembrane helix entry. Molecular dynamic relaxations of the structures in lipid bilayer revealed the role of the proline residue in functionally significant rearrangements of the adjacent juxtamembrane region supporting alternation between protein-protein and protein-lipid interactions of this region that can be triggered by ligand binding. Also, the importance of juxtamembrane SS bonding for signal persistency, and somewhat unusual aspects of transmembrane region interaction with water molecules were demonstrated. CONCLUSIONS: Two alternative dimeric structures of GHR TMD attributed to dormant and active receptor states interchange via allosteric rearrangements of transmembrane helices and extracellular juxtamembrane regions that support coordination between protein-protein and protein-lipid interactions. GENERAL SIGNIFICANCE: This study provides a holistic vision of GHR signal transduction across the membrane emphasizing the role of protein-lipid interactions.

The Conformation of the Epidermal Growth Factor Receptor Transmembrane Domain Dimer Dynamically Adapts to the Local Membrane Environment.

The epidermal growth factor receptor (EGFR) family is an important class of receptor tyrosine kinases, mediating a variety of cellular responses in normal biological processes and in pathological states of multicellular organisms. Different modes of dimerization of the human EGFR transmembrane domain (TMD) in different membrane mimetics recently prompted us to propose a novel signal transduction mechanism based on protein-lipid interaction. However, the experimental evidence for it was originally obtained with slightly different TMD fragments used in the two different mimetics, compromising the validity of the comparison. To eliminate ambiguity, we determined the nuclear magnetic resonance (NMR) structure of the bicelle-incorporated dimer of the EGFR TMD fragment identical to the one previously used in micelles. The NMR results augmented by molecular dynamics simulations confirm the mutual influence of the TMD and lipid environment, as is required for the proposed lipid-mediated activation mechanism. They also reveal the possible functional relevance of a subtle interplay between the concurrent processes in the lipid and protein during signal transduction.

PREDDIMER: a web server for prediction of transmembrane helical dimers.

SUMMARY: Here we present PREDDIMER, a web tool for prediction of dimer structure of transmembrane (TM) helices. PREDDIMER allows (i) reconstruction of a number of dimer structures for given sequence(s) of TM protein fragments, (ii) ranking and filtering of predicted structures according to respective values of a scoring function, (iii) visualization of predicted 3D dimer structures and (iv) visualization of surface hydrophobicity of TM helices and their contacting (interface) regions represented as 2D maps. RESULTS: We implemented online the original PREDDIMER algorithm and benchmarked the server on 11 TM sequences, whose 3D dimer conformations were obtained previously by nuclear magnetic resonance spectroscopy. In the most of tested cases backbone root-mean-square deviations of closest predicted conformations from the experimental reference are below 3 Å. A randomization test displays good anticorrelation (-0.82) between values of the scoring function and statistical significance of the prediction 'by chance'. Going beyond a single dimer conformation, our web tool predicts an ensemble of possible conformations, which may be useful for explanation of a functioning of bitopic membrane proteins, e.g. receptor tyrosine kinases. AVAILABILITY AND IMPLEMENTATION: PREDDIMER can be accessed for free on the web at http://model.nmr.ru/preddimer/ CONTACT: newant@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Role of dimerization efficiency of transmembrane domains in activation of fibroblast growth factor receptor 3.

Mutations in transmembrane (TM) domains of receptor tyrosine kinases are shown to cause a number of inherited diseases and cancer development. Here, we use a combined molecular modeling approach to understand molecular mechanism of effect of G380R and A391E mutations on dimerization of TM domains of human fibroblast growth factor receptor 3 (FGFR3). According to results of Monte Carlo conformational search in the implicit membrane and further molecular dynamics simulations, TM dimer of this receptor is able to form a number of various conformations, which differ significantly by the free energy of association in a full-atom model bilayer. The aforementioned mutations affect dimerization efficiency of TM segments and lead to repopulation of conformational ensemble for the dimer. Particularly, both mutations do not change the dimerization free energy of the predominant (putative "non-active") symmetric conformation of TM dimer, while affect dimerization efficiency of its asymmetric ("intermediate") and alternative symmetric (putative "active") models. Results of our simulations provide novel atomistic prospective of the role of G380 and A391E mutations in dimerization of TM domains of FGFR3 and their consecutive contributions to the activation pathway of the receptor.

Multistate organization of transmembrane helical protein dimers governed by the host membrane.

Association of transmembrane (TM) helices taking place in the cell membrane has an important contribution to the biological function of bitopic proteins, among which receptor tyrosine kinases represent a typical example and a potent target for medical applications. Since this process depends on a complex interplay of different factors (primary structures of TM domains and juxtamembrane regions, composition and phase of the local membrane environment, etc.), it is still far from being fully understood. Here, we present a computational modeling framework, which we have applied to systematically analyze dimerization of 18 TM helical homo- and heterodimers of different bitopic proteins, including the family of epidermal growth factor receptors (ErbBs). For this purpose, we have developed a novel surface-based modeling approach, which not only is able to predict particular conformations of TM dimers in good agreement with experiment but also provides screening of their conformational heterogeneity. Using all-atom molecular dynamics simulations of several of the predicted dimers in different model membranes, we have elucidated a putative role of the environment in selection of particular conformations. Simulation results clearly show that each particular bilayer preferentially stabilizes one of possible dimer conformations, and that the energy gain depends on the interplay between structural properties of the protein and the membrane. Moreover, the character of protein-driven perturbations of the bilayer is reflected in the contribution of a particular membrane to the free energy gain. We have found that the approximated dimerization strength for ErbBs family can be related to their oncogenic ability.

Cytostatic effects of structurally different ginsenosides on yeast cells with altered sterol biosynthesis and transport.

Triterpene glycosides are a diverse group of plant secondary metabolites, consisting of a sterol-like aglycon and one or several sugar groups. A number of triterpene glycosides show membranolytic activity, and, therefore, are considered to be promising antimicrobial drugs. However, the interrelation between their structure, biological activities, and target membrane lipid composition remains elusive. Here we studied the antifungal effects of four Panax triterpene glycosides (ginsenosides) with sugar moieties at the C-3 (ginsenosides Rg3, Rh2), C-20 (compound K), and both (ginsenoside F2) positions in Saccharomyces cerevisiae mutants with altered sterol plasma membrane composition. We observed reduced cytostatic activity of the Rg3 and compound K in the UPC2-1 strain with high membrane sterol content. Moreover, LAM gene deletion reduced yeast resistance to Rg3 and digitonin, another saponin with glycosylated aglycon in the C-3 position. LAM genes encode plasma membrane-anchored StARkin superfamily-member sterol transporters. We also showed that the deletion of the ERG6 gene that inhibits ergosterol biosynthesis at the stage of zymosterol increased the cytostatic effects of Rg3 and Rh2, but not the other two tested ginsenosides. At the same time, in silico simulation revealed that the substitution of ergosterol with zymosterol in the membrane changes the spatial orientation of Rg3 and Rh2 in the membranes. These results imply that the plasma membrane sterol composition defines its interaction with triterpene glycoside depending on their glycoside group position. Our results also suggest that the biological role of membrane-anchored StARkin family protein is to protect eukaryotic cells from triterpenes glycosylated at the C-3 position.

Phospholipase A2 way to hydrolysis: Dint formation, hydrophobic mismatch, and lipid exclusion.

Phospholipase A2 (PLA2) exerts a wide range of biological effects and attracts a lot of attention of researchers. Two sites are involved in manifestation of PLA2 enzymatic activity: catalytic site responsible for substrate binding and fatty acid cleavage from the sn-2 position of a glycerophospholipid, and interface binding site (IBS) responsible for the protein binding to lipid membrane. IBS is formed by positively charged and hydrophobic amino acids on the outer surface of the protein molecule. Understanding the mechanism of PLA2 interaction with the lipid membrane is the most challenging step in biochemistry of this enzyme. We used a combination of experimental and computer simulation techniques to clarify molecular details of bee venom PLA2 interaction with lipid bilayers formed by palmitoyloleoylphosphatidylcholine or dipalmitoylphosphatidylcholine. We found that after initial enzyme contact with the membrane, a network of hydrogen bonds was formed. This led to deformation of the interacting leaflet and dint formation. The bilayer response to the deformation depended on its phase state. In a gel-phase bilayer, diffusion of lipids is restricted therefore chain melting occurred in both leaflets of the bilayer. In the case of a fluid-phase bilayer, lateral diffusion is possible, and lipid polar head groups were excluded from the contact area. As a result, the bilayer became thinner and a large hydrophobic area was formed. We assume that relative ability of a bilayer to come through lipid redistribution process defines the rate of initial stages of the catalysis.

All-d-Enantiomeric Peptide D3 Designed for Alzheimer's Disease Treatment Dynamically Interacts with Membrane-Bound Amyloid-ß Precursors.

Alzheimer's disease (AD) is a severe neurodegenerative pathology with no effective treatment known. Toxic amyloid-ß peptide (Aß) oligomers play a crucial role in AD pathogenesis. All-d-Enantiomeric peptide D3 and its derivatives were developed to disassemble and destroy cytotoxic Aß aggregates. One of the D3-like compounds is approaching phase II clinical trials; however, high-resolution details of its disease-preventing or pharmacological actions are not completely clear. We demonstrate that peptide D3 stabilizing Aß monomer dynamically interacts with the extracellular juxtamembrane region of a membrane-bound fragment of an amyloid precursor protein containing the Aß sequence. MD simulations based on NMR measurement results suggest that D3 targets the amyloidogenic region, not compromising its α-helicity and preventing intermolecular hydrogen bonding, thus creating prerequisites for inhibition of early steps of Aß conversion into ß-conformation and its toxic oligomerization. An enhanced understanding of the D3 action molecular mechanism facilitates development of effective AD treatment and prevention strategies.

Left-handed dimer of EphA2 transmembrane domain: Helix packing diversity among receptor tyrosine kinases.

The Eph receptor tyrosine kinases and their membrane-bound ephrin ligands control a diverse array of cell-cell interactions in the developing and adult organisms. During signal transduction across plasma membrane, Eph receptors, like other receptor tyrosine kinases, are involved in lateral dimerization and subsequent oligomerization presumably with proper assembly of their single-span transmembrane domains. Spatial structure of dimeric transmembrane domain of EphA2 receptor embedded into lipid bicelle was obtained by solution NMR, showing a left-handed parallel packing of the transmembrane helices (535-559)(2). The helices interact through the extended heptad repeat motif L(535)X(3)G(539)X(2)A(542)X(3)V(546)X(2)L(549) assisted by intermolecular stacking interactions of aromatic rings of (FF(557))(2), whereas the characteristic tandem GG4-like motif A(536)X(3)G(540)X(3)G(544) is not used, enabling another mode of helix-helix association. Importantly, a similar motif AX(3)GX(3)G as was found is responsible for right-handed dimerization of transmembrane domain of the EphA1 receptor. These findings serve as an instructive example of the diversity of transmembrane domain formation within the same family of protein kinases and seem to favor the assumption that the so-called rotation-coupled activation mechanism may take place during the Eph receptor signaling. A possible role of membrane lipid rafts in relation to Eph transmembrane domain oligomerization and Eph signal transduction was also discussed.

Probing temperature and capsaicin-induced activation of TRPV1 channel via computationally guided point mutations in its pore and TRP domains.

In a recent computational study, we revealed some mechanistic aspects of TRPV1 (transient receptor potential channel 1) thermal activation and gating and proposed a set of probable functionally important residues - "hot spots" that have not been characterized experimentally yet. In this work, we analyzed TRPV1 point mutants G643A, I679A + A680G, and K688G/P combining molecular modeling, biochemistry, and electrophysiology. The substitution G643A reduced maximal conductivity that resulted in a normal response to moderate stimuli, but a relatively weak response to more intensive activation. I679A + A680G channel was severely toxic for oocytes most probably due to abnormally increased basal activity of the channel ("always open" gates). The replacement K688G presumably facilitated movements of TRP domain and disturbed its coupling to the pore, thus leading to spontaneous activation and enhanced desensitization of the channel. Finally, mutation K688P was suggested to impair TRP domain directed movement, and the mutated channel showed ~100-fold less sensitivity to the capsaicin, enhanced desensitization and weaker activation by the heat. Our results provide a better understanding of TRPV1 thermal and capsaicin-induced activation and gating. These observations provide a structural basis for understanding some aspects of TRPV1 channel functioning and depict potentially pathogenic mutations.

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.

Adaptation of a membrane-active peptide to heterogeneous environment. I. Structural plasticity of the peptide.

Some membrane-active peptides undergo drastic changes of conformation and/or orientation on water-lipid interfaces. Among the most notable examples is penetratin (pAntp), a short cell-penetrating peptide. To delineate the driving forces behind pAntp-membrane interactions, we used, in this series of two papers, a combined modeling approach that includes: (1) molecular dynamics simulations of pAntp in zwitterionic and anionic lipid bilayers, (2) free energy perturbation calculations of model residue-residue contacts, and (3) detailed analysis of spatial hydrophobic/hydrophilic properties of the peptide/membrane systems. In this first article, we consider the role of conformational plasticity of the peptide in different membrane surroundings, as well as the ability of pAntp to form stable specific residue-residue interactions and make contacts with particular lipids. It was shown that pAntp displays a complicated conformational behavior. Basic and aromatic residues of the peptide form energetically favorable pairs in water and apolar environments, which facilitate membrane insertion of the peptide and stabilization of the membrane-bound state. These residues are also capable of "trapping" lipid heads, thereby affecting their dynamics and microscopic organization of the water-lipid interface. The latter effect is much more pronounced in anionic bilayers and might be related to the initial stage of peptide-induced destabilization of lipid bilayers.

Specific refolding pathway of viscumin A chain in membrane-like medium reveals a possible mechanism of toxin entry into cell.

How is a water-soluble globular protein able to spontaneously cross a cellular membrane? It is commonly accepted that it undergoes significant structural rearrangements on the lipid-water interface, thus acquiring membrane binding and penetration ability. In this study molecular dynamics (MD) simulations have been used to explore large-scale conformational changes of the globular viscumin A chain in a complex environment - comprising urea and chloroform/methanol (CHCl3/MeOH) mixture. Being well-packed in aqueous solution, viscumin A undergoes global structural rearrangements in both organic media. In urea, the protein is "swelling" and gradually loses its long-distance contacts, thus resembling the "molten globule" state. In CHCl3/MeOH, viscumin A is in effect turned "inside out". This is accompanied with strengthening of the secondary structure and surface exposure of hydrophobic epitopes originally buried inside the globule. Resulting solvent-adapted models were further subjected to Monte Carlo simulations with an implicit hydrophobic slab membrane. In contrast to only a few point surface contacts in water and two short regions with weak protein-lipid interactions in urea, MD-derived structures in CHCl3/MeOH reveal multiple determinants of membrane interaction. Consequently it is now possible to propose a specific pathway for the structural adaptation of viscumin A with respect to the cell membrane - a probable first step of its translocation into cytoplasmic targets.

Effects of Sterols on the Interaction of SDS, Benzalkonium Chloride, and A Novel Compound, Kor105, with Membranes.

Sterols change the biophysical properties of lipid membranes. Here, we analyzed how sterols affect the activity of widely used antimicrobial membrane-active compounds, sodium dodecyl sulfate (SDS) and benzalkonium chloride (BAC). We also tested a novel benzalkonium-like substance, Kor105. Our data suggest that benzalkonium and Kor105 disturb the ordering of the membrane lipid packaging, and this disturbance is dampened by cholesterol. The disturbance induced by Kor105 is stronger than that induced by BAC because of the higher rigidity of the Kor105 molecule due to a shorter linker between the phenyl group and quaternary nitrogen. On the contrary, individual SDS molecules do not cause the disturbance. Thus, in the tested range of concentrations, SDS-membrane interaction is not influenced by cholesterol. To study how sterols influence the biological effects of these chemicals, we used yeast strains lacking Lam1-4 proteins. These proteins transport sterols from the plasma membrane into the endoplasmic reticulum. We found that the mutants are resistant to BAC and Kor105 but hypersensitive to SDS. Together, our findings show that sterols influence the interaction of SDS versus benzalkonium chloride and Kor105 with the membranes in a completely different manner.

Familial L723P Mutation Can Shift the Distribution between the Alternative APP Transmembrane Domain Cleavage Cascades by Local Unfolding of the Ε-Cleavage Site Suggesting a Straightforward Mechanism of Alzheimer's Disease Pathogenesis.

Alzheimer's disease is an age-related pathology associated with accumulation of amyloid-ß peptides, products of enzymatic cleavage of amyloid-ß precursor protein (APP) by secretases. Several familial mutations causing early onset of the disease have been identified in the APP transmembrane (TM) domain. The mutations influence production of amyloid-ß, but the molecular mechanisms of this effect are unclear. The "Australian" (L723P) mutation located in the C-termini of APP TM domain is associated with autosomal-dominant, early onset Alzheimer's disease. Herein, we describe the impact of familial L723P mutation on the structural-dynamic behavior of APP TM domain studied by high-resolution NMR in membrane-mimicking micelles and augmented by molecular dynamics simulations in explicit lipid bilayer. We found L723P mutation to cause local unfolding of the C-terminal turn of the APP TM domain helix and increase its accessibility to water required for cleavage of the protein backbone by γ-secretase in the ε-site, thus switching between alternative ("pathogenic" and "non-pathogenic") cleavage cascades. These findings suggest a straightforward mechanism of the pathogenesis associated with this mutation, and are of generic import for understanding the molecular-level events associated with APP sequential proteolysis resulting in accumulation of the pathogenic forms of amyloid-ß. Moreover, age-related onset of Alzheimer's disease can be explained by a similar mechanism, where the effect of mutation is emulated by the impact of local environmental factors, such as oxidative stress and/or membrane lipid composition. Knowledge of the mechanisms regulating generation of amyloidogenic peptides of different lengths is essential for development of novel treatment strategies of the Alzheimer's disease.

Improving therapeutic potential of antibacterial spider venom peptides: coarse-grain molecular dynamics guided approach.

AIM: Spider venom is a rich source of antibacterial peptides, whose hemolytic activity is often excessive. METHODOLOGY: How to get rid of it? Using latarcins from Lachesana tarabaevi and oxyopinin Oxt 4a from Oxyopes takobius spider venoms we performed coarse-grained molecular dynamics simulations of these peptides in the presence of lipid bilayers, mimicking erythrocyte membranes. This identified hemolytically active fragments within Oxt 4a and latarcins. Then, we synthesized five 20-residue peptides, containing different parts of the Oxt 4a and latarcin-1 sequence, carrying mutations within the identified regions. CONCLUSION: The antibacterial and hemolytic tests suggested that the three of the synthesized peptides demonstrated substantial decrease in hemolytic activity, retaining, or even exceeding antibacterial potential of the parent peptides.

Impact of membrane partitioning on the spatial structure of an S-type cobra cytotoxin.

Cobra cytotoxins (CTs) belong to the three-fingered protein family. They are classified into S- and P-types, the latter exhibiting higher membrane-perturbing capacity. In this work, we investigated the interaction of CTs with phospholipid bilayers, using coarse-grained (CG) and full-atom (FA) molecular dynamics (MD). The object of this work is a CT of an S-type, cytotoxin I (CT1) from N.oxiana venom. Its spatial structure in aqueous solution and in the micelles of dodecylphosphocholine (DPC) were determined by 1H-NMR spectroscopy. Then, via CG- and FA MD-computations, we evaluated partitioning of CT1 molecule into palmitoyloleoylphosphatidylcholine (POPC) membrane, using the toxin spatial models, obtained either in aqueous solution, or detergent micelle. The latter model exhibits minimal structural changes upon partitioning into the membrane, while the former deviates from the starting conformation, loosing the tightly bound water molecule in the loop-2. These data show that the structural changes elicited by CT1 molecule upon incorporation into DPC micelle take place likely in the lipid membrane, although the mode of the interaction of this toxin with DPC micelle (with the tips of the all three loops) is different from its mode in POPC membrane (primarily with the tip of the loop-1 and both the tips of the loop-1 and loop-2).

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.