The architecture of transmembrane and cytoplasmic juxtamembrane regions of Toll-like receptors

Toll-like receptors (TLRs) are the important participants of the innate immune response. Their spatial organization is well studied for the ligand-binding domains, while a lot of questions remain unanswered for the membrane and cytoplasmic regions of the proteins. Here we use solution NMR spectroscopy and computer simulations to investigate the spatial structures of transmembrane and cytoplasmic juxtamembrane regions of TLR2, TLR3, TLR5, and TLR9. According to our data, all the proteins reveal the presence of a previously unreported structural element, the cytoplasmic hydrophobic juxtamembrane α-helix. As indicated by the functional tests in living cells and bioinformatic analysis, this helix is important for receptor activation and plays a role, more complicated than a linker, connecting the transmembrane and cytoplasmic parts of the proteins.

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.

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.

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.

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.

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.

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.

BILMIX: a new approach to restore the size polydispersity and electron density profiles of lipid bilayers from liposomes using small-angle X-ray scattering data

Small-angle X-ray scattering (SAXS) is one of the major tools for the study of model membranes, but interpretation of the scattering data remains non-trivial. Current approaches allow the extraction of some structural parameters and the electron density profile of lipid bilayers. Here it is demonstrated that parametric modelling can be employed to determine the polydispersity of spherical or ellipsoidal vesicles and describe the electron density profile across the lipid bilayer. This approach is implemented in the computer program BILMIX. BILMIX delivers a description of the electron density of a lipid bilayer from SAXS data and simultaneously generates the corresponding size distribution of the unilamellar lipid vesicles.

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.

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 (CHCl₃/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 CHCl₃/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 CHCl₃/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.

Probing the effect of membrane contents on transmembrane protein-protein interaction using solution NMR and computer simulations

The interaction between the secondary structure elements is the key process, determining the spatial structure and activity of a membrane protein. Transmembrane (TM) helix-helix interaction is known to be especially important for the function of so-called type I or bitopic membrane proteins. In the present work, we present the approach to study the helix-helix interaction in the TM domains of membrane proteins in various lipid environment using solution NMR spectroscopy and phospholipid bicelles. The technique is based on the ability of bicelles to form particles with the size, depending on the lipid/detergent ratio. To implement the approach, we report the experimental parameters of "ideal bicelle" models for four kinds of zwitterionic phospholipids, which can be also used in other structural studies. We show that size of bicelles and type of the rim-forming detergent do not affect substantially the spatial structure and stability of the model TM dimer. On the other hand, the effect of bilayer thickness on the free energy of the dimer is dramatic, while the structure of the protein is unchanged in various lipids with fatty chains having a length from 12 to 18 carbon atoms. The obtained data is analyzed using the computer simulations to find the physical origin of the observed effects.

Pore formation in lipid membrane II: Energy landscape under external stress

Lipid membranes are extremely stable envelopes allowing cells to survive in various environments and to maintain desired internal composition. Membrane permeation through formation of transversal pores requires substantial external stress. Practically, pores are usually formed by application of lateral tension or transmembrane voltage. Using the same approach as was used for obtaining continuous trajectory of pore formation in the stress-less membrane in the previous article, we now consider the process of pore formation under the external stress. The waiting time to pore formation proved a non-monotonous function of the lateral tension, dropping from infinity at zero tension to a minimum at the tension of several millinewtons per meter. Transmembrane voltage, on the contrary, caused the waiting time to decrease monotonously. Analysis of pore formation trajectories for several lipid species with different spontaneous curvatures and elastic moduli under various external conditions provided instrumental insights into the mechanisms underlying some experimentally observed phenomena.

Pore formation in lipid membrane I: Continuous reversible trajectory from intact bilayer through hydrophobic defect to transversal pore

Lipid membranes serve as effective barriers allowing cells to maintain internal composition differing from that of extracellular medium. Membrane permeation, both natural and artificial, can take place via appearance of transversal pores. The rearrangements of lipids leading to pore formation in the intact membrane are not yet understood in details. We applied continuum elasticity theory to obtain continuous trajectory of pore formation and closure, and analyzed molecular dynamics trajectories of pre-formed pore reseal. We hypothesized that a transversal pore is preceded by a hydrophobic defect: intermediate structure spanning through the membrane, the side walls of which are partially aligned by lipid tails. This prediction was confirmed by our molecular dynamics simulations. Conversion of the hydrophobic defect into the hydrophilic pore required surmounting some energy barrier. A metastable state was found for the hydrophilic pore at the radius of a few nanometers. The dependence of the energy on radius was approximately quadratic for hydrophobic defect and small hydrophilic pore, while for large radii it depended on the radius linearly. The pore energy related to its perimeter, line tension, thus depends of the pore radius. Calculated values of the line tension for large pores were in quantitative agreement with available experimental data.

HER2 Transmembrane Domain Dimerization Coupled with Self-Association of Membrane-Embedded Cytoplasmic Juxtamembrane Regions

Receptor tyrosine kinases of the human epidermal growth factor receptor (HER or ErbB) family transduce biochemical signals across plasma membrane, playing a significant role in vital cellular processes and in various cancers. Inactive HER/ErbB receptors exist in equilibrium between the monomeric and unspecified pre-dimerized states. After ligand binding, the receptors are involved in strong lateral dimerization with proper assembly of their extracellular ligand-binding, single-span transmembrane, and cytoplasmic kinase domains. The dimeric conformation of the HER2 transmembrane domain that is believed to support the cytoplasmic kinase domain configuration corresponding to the receptor active state was previously described in lipid bicelles. Here we used high-resolution NMR spectroscopy in another membrane-mimicking micellar environment and identified an alternative HER2 transmembrane domain dimerization coupled with self-association of membrane-embedded cytoplasmic juxtamembrane region. Such a dimerization mode appears to be capable of effectively inhibiting the receptor kinase activity. This finding refines the molecular mechanism regarding the signal propagation steps from the extracellular to cytoplasmic domains of HER/ErbB receptors.

Liquid but durable: molecular dynamics simulations explain the unique properties of archaeal-like membranes

Archaeal plasma membranes appear to be extremely durable and almost impermeable to water and ions, in contrast to the membranes of Bacteria and Eucaryota. Additionally, they remain liquid within a temperature range of 0-100°C. These are the properties that have most likely determined the evolutionary fate of Archaea, and it may be possible for bionanotechnology to adopt these from nature. In this work, we use molecular dynamics simulations to assess at the atomistic level the structure and dynamics of a series of model archaeal membranes with lipids that have tetraether chemical nature and "branched" hydrophobic tails. We conclude that the branched structure defines dense packing and low water permeability of archaeal-like membranes, while at the same time ensuring a liquid-crystalline state, which is vital for living cells. This makes tetraether lipid systems promising in bionanotechnology and material science, namely for design of new and unique membrane nanosystems.

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.

Structural, dynamic, and functional aspects of helix association in membranes: a computational view

This review surveys recent achievements of molecular computer modeling in understanding spatial structure, dynamics, and mechanisms of functioning of transmembrane α-helical dimers in membranes. The factors driving self-association of hydrophobic helices in the membrane milieu are considered with examples of their applications to biologically relevant problems. The emphasis is made on the recent results, which help to understand important aspects of structure-function relations for these systems and their biological activity. Limitations and shortcomings of the methods, along with their perspectives in design of new membrane active agents, are discussed.

COMPUTER SIMULATIONS OF MEMBRANE-LYTIC PEPTIDES: PERSPECTIVES IN DRUG DESIGN

Structure activity relationships were investigated for membrane-lytic peptides (MLP) Ltc1 and Ltc2a from the latarcin family. The peptides were studied via long-term molecular dynamics (MD) simulations in different membrane environments (detergent micelles, mixed lipid bilayers mimiking eukaryotic and bacterial membranes). The calculated structure of Ltc2a in sodium dodecyl sulfate micelle agrees well with the data obtained by 1H-NMR spectroscopy. This validates the applied modeling approach. The binding mode of MLPs is governed by several factors: (i) the membrane surface curvature; (ii) the conformational plasticity and hydrophobic organization of the peptide, which depend on the arrangement of charged, non-polar and helix-breaking residues in the amino acid sequence. In contrast to Ltc1, insertion of Ltc2a into model membranes induces significant changes in dynamic behavior of lipids in the contact region. Such a prominent membrane destabilization correlates with high membrane-lytic activity of Ltc2a. In all cases the "membrane response" has a local character and is caused by formation of specific peptide-lipid contacts. Results of MD simulations of Ltc2a in model membranes were used to develop a number of its analogs with predefined activity.

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.

Structure elucidation of dimeric transmembrane domains of bitopic proteins

The interaction between transmembrane helices is of great interest because it directly determines biological activity of a membrane protein. Either destroying or enhancing such interactions can result in many diseases related to dysfunction of different tissues in human body. One much studied form of membrane proteins known as bitopic protein is a dimer containing two membrane-spanning helices associating laterally. Establishing structure-function relationship as well as rational design of new types of drugs targeting membrane proteins requires precise structural information about this class of objects. At present time, to investigate spatial structure and internal dynamics of such transmembrane helical dimers, several strategies were developed based mainly on a combination of NMR spectroscopy, optical spectroscopy, protein engineering and molecular modeling. These approaches were successfully applied to homo- and heterodimeric transmembrane fragments of several bitopic proteins, which play important roles in normal and in pathological conditions of human organism.

Computer simulations and modeling-assisted ToxR screening in deciphering 3D structures of transmembrane alpha-helical dimers: ephrin receptor A1

Membrane-spanning segments of numerous proteins (e.g. receptor tyrosine kinases) represent a novel class of pharmacologically important targets, whose activity can be modulated by specially designed artificial peptides, the so-called interceptors. Rational construction of such peptides requires understanding of the main factors driving peptide-peptide association in lipid membranes. Here we present a new method for rapid prediction of the spatial structure of transmembrane (TM) helix-helix complexes. It is based on computer simulations in membrane-like media and subsequent refinement/validation of the results using experimental studies of TM helix dimerization in a bacterial membrane by means of the ToxR system. The approach was applied to TM fragments of the ephrin receptor A1 (EphA1). A set of spatial structures of the dimer was proposed based on Monte Carlo simulations in an implicit membrane followed by molecular dynamics relaxation in an explicit lipid bilayer. The resulting models were employed for rational design of wild-type and mutant genetic constructions for ToxR assays. The computational and the experimental data are self-consistent and provide an unambiguous spatial model of the TM dimer of EphA1. The results of this work can be further used to develop new biologically active 'peptide interceptors' specifically targeting membrane domains of proteins.

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.

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.

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.

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.

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.

Spatial Structure and Activity Mechanism of a Novel Spider Antimicrobial Peptide,

Latarcins (Ltc), linear peptides (ca. 25 amino acid long) isolated from the venom of the Lachesana tarabaevi spider, exhibit a broad-spectrum antibacterial activity, most likely acting on the bacterial plasmatic membrane. We study the structure-activity relationships in the series of these compounds. At the first stage, we investigated the spatial structure of one of the peptides, Ltc2a, and its mode of membrane perturbation. This was done by a combination of experimental and theoretical methods. The approach includes (i) structural study of the peptide by CD spectroscopy in phospholipid liposomes and by (1)H NMR in detergent micelles, (ii) determination of the effect on the liposomes by a dye leakage fluorescent assay and (31)P NMR spectroscopy, (iii) refinement of the NMR-derived spatial structure via Monte Carlo simulations in an implicit water-octanol slab, and (iv) calculation of the molecular hydrophobicity potential. The molecule of Ltc2a was found to consist of two helical regions (residues 3-9 and 13-21) connected via a poorly ordered fragment. The effect of the peptide on the liposomes suggests the carpet mechanism of the membrane deterioration. This is also supported by the analysis of hydrophobic/hydrophilic characteristics of Ltc2a and homologous antimicrobial peptides. These peptides exhibiting a helix-hinge-helix structural motif are characterized by a distinct and feebly marked amphiphilicity of their N- and C-terminal helices, respectively, and by a hydrophobicity gradient along the peptide chain. The approach we suggested may be useful in studying not only other latarcins but also a wider class of membrane-active peptides.

Association of transmembrane helices: what determines assembling of a dimer?

Self-association of two hydrophobic alpha-helices is studied via unrestrained Monte Carlo (MC) simulations in a hydrophobic slab described by an effective potential. The system under study represents two transmembrane (TM) segments of human glycophorin A (GpA), which form homo-dimers in membranes. The influence of TM electrostatic potential, thickness and hydrophobicity degree of lipid bilayer is investigated. It is shown that the membrane environment stabilizes alpha-helical conformation of GpA monomers, induces their TM insertion and facilitates inter-helical contacts. Head-to-head orientation of the helices is promoted by the voltage difference across the membrane. Subsequent "fine-tuned" assembling of the dimer is mediated by van der Waals interactions. Only the models of dimer, calculated in a hydrophobic slab with applied voltage agree with experimental data, while simulations in vacuo or without TM voltage fail to give reasonable results. The moderate structural heterogeneity of GpA dimers (existence of several groups of states with close energies) is proposed to reflect their equilibrium dynamics in membrane-mimics. The calculations performed for GpA mutants G83A and G86L permit rationalization of mutagenesis data for them. The results of Monte Carlo simulations critically depend on the parameters of the membrane model: adequate description of helix association is achieved in the water-cyclohexane-water system with the membrane thickness 30-34 A, while in membranes with different hydrophobicities and thickness unrealistic conformations of the dimer are found. The computational approach permits efficient prediction of TM helical oligomers based solely on the sequences of interacting peptides.

Effect of Lipid Composition on the “Membrane Response” Induced by a Fusion Peptide

To understand the initial stages of membrane destabilization induced by viral proteins, the factors important for binding of fusion peptides to cell membranes must be identified. In this study, effects of lipid composition on the mode of peptides' binding to membranes are explored via molecular dynamics (MD) simulations of the peptide E5, a water-soluble analogue of influenza hemagglutinin fusion peptide, in two full-atom hydrated lipid bilayers composed of dimyristoyl- and dipalmitoylphosphatidylcholine (DMPC and DPPC, respectively). The results show that, although the peptide has a common folding motif in both systems, it possesses different modes of binding. The peptide inserts obliquely into the DMPC membrane mainly with its N-terminal alpha helix, while in DPPC, the helix lies on the lipid/water interface, almost parallel to the membrane surface. The peptide seriously affects structural and dynamical parameters of surrounding lipids. Thus, it induces local thinning of both bilayers and disordering of acyl chains of lipids in close proximity to the binding site. The "membrane response" significantly depends upon lipid composition: distortions of DMPC bilayer are more pronounced than those in DPPC. Implications of the observed effects to molecular events on initial stages of membrane destabilization induced by fusion peptides are discussed.

Development of the force field parameters for phosphoimidazole and phosphohistidine

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.

Pyrenemethyl ara-uridine-2'-carbamate: a strong interstrand excimer in the major groove of a DNA duplex

The synthesis of new nucleoside derivatives, ara-uridine-2'-carbamates, and their incorporation into synthetic DNA oligomers is described. The modification directs ligands into the major groove of duplex DNA and somewhat destabilizes the duplexes of modified oligonucleotides with complementary DNA or RNA. In the case of pyrenemethyl carbamate modification in DNA-DNA duplexes, the destabilization is considerably reduced. The pyrenemethyl derivative also shows remarkable spectral properties: a "reversed" absorbance change for pyrene at 350 nm in the course of denaturation of the DNA duplex, as compared to the change seen in the nucleotide absorbance at 260 nm. This derivatization also causes pronounced sequence-dependent excimer formation in the major groove.