Identification and Structural Characterization of a New Three-Finger Toxin Hemachatoxin from Hemachatus haemachatus Venom

Snake venoms are rich sources of biologically active proteins and polypeptides. Three-finger toxins are non-enzymatic proteins present in elapid (cobras, kraits, mambas and sea snakes) and colubrid venoms. These proteins contain four conserved disulfide bonds in the core to maintain the three-finger folds. Although all three-finger toxins have similar fold, their biological activities are different. A new three-finger toxin (hemachatoxin) was isolated from Hemachatus haemachatus (Ringhals cobra) venom. Its amino acid sequence was elucidated, and crystal structure was determined at 2.43 Å resolution. The overall fold is similar to other three-finger toxins. The structure and sequence analysis revealed that the fold is maintained by four highly conserved disulfide bonds. It exhibited highest similarity to particularly P-type cardiotoxins that are known to associate and perturb the membrane surface with their lipid binding sites. Also, the increased B value of hemachotoxin loop II suggests that loop II is flexible and may remain flexible until its interaction with membrane phospholipids. Based on the analysis, we predict hemachatoxin to be cardiotoxic/cytotoxic and our future experiments will be directed to characterize the activity of hemachatoxin.     

Membrane insertion of the heptameric staphylococcal alpha-toxin pore. A domino-like structural transition that is allosterically modulated by the target cell membrane

Staphylococcal alpha-toxin forms heptameric pores on eukaryotic cells. After binding to the cell membrane in its monomeric form, the toxin first assembles into a heptameric pre-pore. Subsequently, the pre-pore transforms into the final pore by membrane insertion of an amphipathic beta-barrel, which comprises the "central loop" domains of all heptamer subunits. The process of membrane insertion was analyzed here using a set of functionally altered toxin mutants. The results show that insertion may be initiated within an individual protomer when its NH2 terminus activates its central loop. The activated state is then shared with the central loops of the residual heptamer subunits, which results in cooperative membrane penetration. This cooperation of the central loops commences while these are still remote from the lipid bilayer. Nevertheless, it is subject to modulation by the target membrane, which therefore acts across a distance much like an allosteric effector. However, while allosteric transitions usually are reversible, membrane insertion of alpha-toxin is an irreversible event, and we show here that it can proceed to completion in a domino-like fashion when triggered by as little as a single foreign atom within the entire heptamer.     

Insertion of bitopic membrane proteins into the inner membrane of mitochondria involves an export step from the matrix

The mitochondrial inner membrane contains a large number of polytopic proteins that are derived from prokaryotic ancestors of mitochondria. Little is known about the intramitochondrial sorting of these proteins. We chose two proteins of known topology as examples to study the pathway of insertion into the inner membrane; Mrs2 and Yta10 are bitopic proteins that expose negatively charged loops of different complexity into the intermembrane space. Here we show that both Mrs2 and Yta10 transiently accumulate as sorting intermediates in the matrix before they integrate into the inner membrane. The sorting pathway of both proteins can be separated into two sequential reactions: (i) import into the matrix and (ii) insertion from the matrix into the inner membrane. The latter process was found to depend on the membrane potential and, in this respect, is similar to the insertion of membrane proteins in bacteria. A comparison of the charge distribution of intermembrane space loops in a variety of mitochondrial inner membrane proteins suggests that this mode of "conservative sorting" might be the typical insertion route for polytopic inner membrane proteins that originated from bacterial ancestors.     

Membrane Insertion of Marginally Hydrophobic Transmembrane Helices Depends on Sequence Context

In mammalian cells, most integral membrane proteins are initially inserted into the endoplasmic reticulum membrane by the so-called Sec61 translocon. However, recent predictions suggest that many transmembrane helices (TMHs) in multispanning membrane proteins are not sufficiently hydrophobic to be recognized as such by the translocon. In this study, we have screened 16 marginally hydrophobic TMHs from membrane proteins of known three-dimensional structure. Indeed, most of these TMHs do not insert efficiently into the endoplasmic reticulum membrane by themselves. To test if loops or TMHs immediately upstream or downstream of a marginally hydrophobic helix might influence the insertion efficiency, insertion of marginally hydrophobic helices was also studied in the presence of their neighboring loops and helices. The results show that flanking loops and nearest-neighbor TMHs are sufficient to ensure the insertion of many marginally hydrophobic helices. However, for at least two of the marginally hydrophobic helices, the local interactions are not enough, indicating that post-insertional rearrangements are involved in the folding of these proteins.     

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.     

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.     

Glycosphingolipid-facilitated membrane insertion and internalization of cobra cardiotoxin. The sulfatide.cardiotoxin complex structure in a membrane-like environment suggests a lipid-dependent cell-penetrating mechanism for membrane binding polypeptides

Cobra cardiotoxins, a family of basic polypeptides having lipid- and heparin-binding capacities similar to the cell-penetrating peptides, induce severe tissue necrosis and systolic heart arrest in snakebite victims. Whereas cardiotoxins are specifically retained on the cell surface via heparan sulfate-mediated processes, their lipid binding ability appears to be responsible, at least in part, for cardiotoxin-induced membrane leakage and cell death. Although the exact role of lipids involved in toxin-mediated cytotoxicity remains largely unknown, monoclonal anti-sulfatide antibody O4 has recently been shown to inhibit the action of CTX A3, the major cardiotoxin from Taiwan cobra venom, on cardiomyocytes by preventing cardiotoxin-induced membrane leakage and CTX A3 internalization into mitochondria. Here, we show that anti-sulfatide acts by blocking the binding of CTX A3 to the sulfatides in the plasma membrane to prevent sulfatide-dependent CTX A3 membrane pore formation and internalization. We also describe the crystal structure of a CTX A3-sulfatide complex in a membrane-like environment at 2.3 angstroms resolution. The unexpected orientation of the sulfatide fatty chains in the structure allows prediction of the mode of toxin insertion into the plasma membrane. CTX A3 recognizes both the headgroup and the ceramide interfacial region of sulfatide to induce a lipid conformational change that may play a key role in CTX A3 oligomerization and cellular internalization. This proposed lipid-mediated toxin translocation mechanism may also shed light on the cellular uptake mechanism of the amphiphilic cell-penetrating peptides known to involve multiple internalization pathways.     

Membrane insertion pathway of annexin B12: thermodynamic and kinetic characterization by fluorescence correlation spectroscopy and fluorescence quenching

Experimental determination of the free energy stabilizing the structure of membrane proteins in their native lipid environment is undermined by the lack of appropriate methods and suitable model systems. Annexin B12 (ANX) is a soluble protein which reversibly inserts into lipid membranes under mildly acidic conditions, which makes it a good experimental model for thermodynamic studies of folding and stability of membrane proteins. Here we apply fluorescence correlation spectroscopy for quantitative analysis of ANX partitioning into large unilamellar vesicles containing either 25% or 75% anionic lipids. Membrane binding of ANX results in changes of autocorrelation time and amplitude, both of which are used in quantitative analysis. The thermodynamic scheme describing acid-induced membrane interactions of ANX considers two independent processes: pH-dependent formation of a membrane-competent form near the membrane interface and its insertion into the lipid bilayer. Our novel fluorescence lifetime topology method demonstrates that the insertion proceeds via an interfacial refolded intermediate state, which can be stabilized by anionic lipids. Lipid titration measurements are used to determine the free energy of both transmembrane insertion and interfacial penetration, which are found to be similar, approximately -10-12 kcal/mol. The formation of the membrane-competent form, examined in a lipid saturation experiment, was found to depend on the local proton concentration near the membrane interface, occurring with pK = 4.3 and involving the protonation of two residues. Our results demonstrate that fluorescence correlation spectroscopy is a convenient tool for the quantitative characterization of the energetics of transmembrane insertion and that pH-triggered ANX insertion is a useful model for studying the thermodynamic stability of membrane proteins.     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

Crystallographic and molecular dynamics analysis of loop motions unmasking the peptidoglycan-binding site in stator protein MotB of flagellar motor

Background

The C-terminal domain of MotB (MotB-C) shows high sequence similarity to outer membrane protein A and related peptidoglycan (PG)-binding proteins. It is believed to anchor the power-generating MotA/MotB stator unit of the bacterial flagellar motor to the peptidoglycan layer of the cell wall. We previously reported the first crystal structure of this domain and made a puzzling observation that all conserved residues that are thought to be essential for PG recognition are buried and inaccessible in the crystal structure. In this study, we tested a hypothesis that peptidoglycan binding is preceded by, or accompanied by, some structural reorganization that exposes the key conserved residues.

Methodology/Principal Findings

We determined the structure of a new crystalline form (Form B) of Helicobacter pylori MotB-C. Comparisons with the existing Form A revealed conformational variations in the petal-like loops around the carbohydrate binding site near one end of the β-sheet. These variations are thought to reflect natural flexibility at this site required for insertion into the peptidoglycan mesh. In order to understand the nature of this flexibility we have performed molecular dynamics simulations of the MotB-C dimer. The results are consistent with the crystallographic data and provide evidence that the three loops move in a concerted fashion, exposing conserved MotB residues that have previously been implicated in binding of the peptide moiety of peptidoglycan.

Conclusion/Significance

Our structural analysis provides a new insight into the mechanism by which MotB inserts into the peptidoglycan mesh, thus anchoring the power-generating complex to the cell wall.     

Molecular dynamics simulations of discoidal bilayers assembled from truncated human lipoproteins

Human apolipoprotein A-1 (apo A-1) is the major protein component of high-density lipoproteins. The apo A-1 lipid-binding domain was used as a template for the synthesis of amphipathic helical proteins termed membrane scaffold proteins, employed to self-assemble soluble monodisperse discoidal particles called Nanodiscs. In these particles, membrane scaffold proteins surround a lipid bilayer in a belt-like fashion forming bilayer disks of discrete size and composition. Here we investigate the structure of Nanodiscs through molecular dynamics simulations in which Nanodiscs were built from scaffold proteins of various lengths. The simulations showed planar or deformed Nanodiscs depending on optimal length and alignment of the scaffold proteins. Based on mean surface area per lipid calculations, comparison of small-angle x-ray scattering curves, and the relatively planar shape of Nanodiscs made from truncated scaffold proteins, one can conclude that the first 17 to 18 residues of the 200-residue apo A-1 lipid-binding domain are not involved in formation of the protein "belts" surrounding the lipid bilayer. To determine whether the addition of an integral membrane protein has an effect on the overall structure of a Nanodisc, bacteriorhodopsin was embedded into a Nanodisc and simulated using molecular dynamics, revealing a planar disk with a slightly rectangular shape.     

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.     

Predicting transmembrane beta-barrels and interstrand residue interactions from sequence

Transmembrane beta-barrel (TMB) proteins are embedded in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. The cellular location and functional diversity of beta-barrel outer membrane proteins (omps) makes them an important protein class. At the present time, very few nonhomologous TMB structures have been determined by X-ray diffraction because of the experimental difficulty encountered in crystallizing transmembrane proteins. A novel method using pairwise interstrand residue statistical potentials derived from globular (nonouter membrane) proteins is introduced to predict the supersecondary structure of transmembrane beta-barrel proteins. The algorithm transFold employs a generalized hidden Markov model (i.e., multitape S-attribute grammar) to describe potential beta-barrel supersecondary structures and then computes by dynamic programming the minimum free energy beta-barrel structure. Hence, the approach can be viewed as a "wrapping" component that may capture folding processes with an initiation stage followed by progressive interaction of the sequence with the already-formed motifs. This approach differs significantly from others, which use traditional machine learning to solve this problem, because it does not require a training phase on known TMB structures and is the first to explicitly capture and predict long-range interactions. TransFold outperforms previous programs for predicting TMBs on smaller (相似文献   

Interaction of Cardiotoxin A5 with Membrane: Role of Conformational Heterogeneity and Hydrophobic Properties

The hypothesis that local conformational differences of the snake venom cardiotoxins (cytotoxins, CT) may play a significant role in their interaction with membrane was tested by molecular modeling of the behavior of the CT A5 from the venom of Naja atra in water and at the water–membrane interface. Two models of the CT A5 spatial structure are known: the first was obtained by X-ray analysis and the second, by NMR studies in solution. A molecular dynamics (MD) analysis demonstrated that loop II of the toxin has a fixed -like shape in water, which does not depend on its initial structure. An interaction of the experimentally derived (X-ray and NMR) conformations and MD simulated conformations of CT A5 with the lipid bilayer was studied by the Monte Carlo method using the previously developed model of the implicit membrane. It is found that: (1) unlike the previously studied CT2 from the venom of cobra Naja oxiana, CT A5 has only loops I and II bound to the membrane with the involvement of a lesser number of hydrophobic residues. (2) A long hydrophobic area is formed on the surface of CT A5 due to the -like shape of loop II and the arrangement of loop I in proximity to loop II. This hydrophobic area favors the toxin embedding into the lipid bilayer. (3) The toxin retains its conformation upon interaction with the membrane. (4) The CT A5 molecule has close values of the potential energy in the membrane and in aqueous environment, which suggests dynamic character of the binding. The results of the molecular modeling indicate a definite configuration of loops I and II and, consequently, a specific character of distribution of polar and apolar properties on the toxin surface, which turns out to be the most energetically favorable.     

Membrane-bound orientation and position of the synaptotagmin C2B domain determined by site-directed spin labeling

Site-directed spin labeling is used to determine the orientation and depth of insertion of the second C2 domain from synaptotagmin I (C2B) into membrane vesicles composed of phosphatidylcholine (PC) and phosphatidylserine (PS). EPR line shapes of spin-labeled mutants located with the Ca(2+)-binding loops of C2B broaden in the presence of Ca(2+) and PC/PS vesicles, indicating that these loops undergo a Ca(2+)-dependent insertion into the membrane interface. Power saturation of the EPR spectra provides a position for each spin-labeled site along the bilayer normal, and these EPR-derived distance constraints, along with a high-resolution structure of the C2B domain, are used to generate a model for the domain orientation and position at the membrane interface. Our data show that the isolated C2B domain from synaptotagmin I penetrates PC/PS membranes, and that the backbone of Ca(2+)-binding loops 1 and 3 is inserted below the level of a plane defined by the lipid phosphates. The side chains of several loop residues are within the bilayer interior, and both Ca(2+)-binding sites are positioned near a plane defined by the lipid phosphates. A Tb(3+)-based fluorescence assay is used to compare the membrane affinity of the C2B domain to that of the first synaptotagmin C2 domain (C2A). Both C2A and C2B bind PC/PS (75:25) membrane vesicles with a micromolar lipid affinity in the presence of metal ion. These results indicate that C2A and C2B have a similar membrane affinity and position when bound to PC/PS (75:25) membrane vesicles. EPR spectroscopy indicates that the C2B domain has different interactions with PC/PS membranes containing 1 mol % phosphatidylinositol 4,5-bisphosphate.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

Surfactant Bilayers Maintain Transmembrane Protein Activity

In vitro studies of membrane proteins are of interest only if their structure and function are significantly preserved. One approach is to insert them into the lipid bilayers of highly viscous cubic phases rendering the insertion and manipulation of proteins difficult. Less viscous lipid sponge phases are sometimes used, but their relatively narrow domain of existence can be easily disrupted by protein insertion. We present here a sponge phase consisting of nonionic surfactant bilayers. Its extended domain of existence and its low viscosity allow easy insertion and manipulation of membrane proteins. We show for the first time, to our knowledge, that transmembrane proteins, such as bacteriorhodopsin, sarcoplasmic reticulum Ca²⁺ATPase (SERCA1a), and its associated enzymes, are fully active in a surfactant phase.     

Surfactant Bilayers Maintain Transmembrane Protein Activity

In vitro studies of membrane proteins are of interest only if their structure and function are significantly preserved. One approach is to insert them into the lipid bilayers of highly viscous cubic phases rendering the insertion and manipulation of proteins difficult. Less viscous lipid sponge phases are sometimes used, but their relatively narrow domain of existence can be easily disrupted by protein insertion. We present here a sponge phase consisting of nonionic surfactant bilayers. Its extended domain of existence and its low viscosity allow easy insertion and manipulation of membrane proteins. We show for the first time, to our knowledge, that transmembrane proteins, such as bacteriorhodopsin, sarcoplasmic reticulum Ca²⁺ATPase (SERCA1a), and its associated enzymes, are fully active in a surfactant phase.     

A generalized born implicit-membrane representation compared to experimental insertion free energies

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.     

Membrane Topology and Insertion of Membrane Proteins: Search for Topogenic Signals

Integral membrane proteins are found in all cellular membranes and carry out many of the functions that are essential to life. The membrane-embedded domains of integral membrane proteins are structurally quite simple, allowing the use of various prediction methods and biochemical methods to obtain structural information about membrane proteins. A critical step in the biosynthetic pathway leading to the folded protein in the membrane is its insertion into the lipid bilayer. Understanding of the fundamentals of the insertion and folding processes will significantly improve the methods used to predict the three-dimensional membrane protein structure from the amino acid sequence. In the first part of this review, biochemical approaches to elucidate membrane protein topology are reviewed and evaluated, and in the second part, the use of similar techniques to study membrane protein insertion is discussed. The latter studies search for signals in the polypeptide chain that direct the insertion process. Knowledge of the topogenic signals in the nascent chain of a membrane protein is essential for the evaluation of membrane topology studies.