Membrane-induced folding and structure of membrane-bound annexin A1 N-terminal peptides: implications for annexin-induced membrane aggregation

Annexins constitute a family of calcium-dependent membrane-binding proteins and can be classified into two groups, depending on the length of the N-terminal domain unique for each individual annexin. The N-terminal domain of annexin A1 can adopt an α-helical conformation and has been implicated in mediating the membrane aggregation behavior of this protein. Although the calcium-independent interaction of the annexin A1 N-terminal domain has been known for some time, there was no structural information about the membrane interaction of this secondary membrane-binding site of annexin A1. This study used circular dichroism spectroscopy to show that a rat annexin A1 N-terminal peptide possesses random coil structure in aqueous buffer but an α-helical structure in the presence of small unilamellar vesicles. The binding of peptides to membranes was confirmed by surface pressure (Langmuir film balance) measurements using phosphatidylcholine/phosphatidylserine monolayers, which show a significant increase after injection of rat annexin A1 N-terminal peptides. Lamellar neutron diffraction with human and rat annexin A1 N-terminal peptides reveals an intercalation of the helical peptides with the phospholipid bilayer, with the helix axis lying parallel to the surface of membrane. Our findings confirm that phospholipid membranes assist the folding of the N-terminal peptides into α-helical structures and that this conformation enables favorable direct interactions with the membrane. The results are consistent with the hypothesis that the N-terminal domain of annexin A1 can serve as a secondary membrane binding site in the process of membrane aggregation by providing a peripheral membrane anchor.     

Contribution of Long-Range Interactions to the Secondary Structure of an Unfolded Globin

This work explores the effect of long-range tertiary contacts on the distribution of residual secondary structure in the unfolded state of an α-helical protein. N-terminal fragments of increasing length, in conjunction with multidimensional nuclear magnetic resonance, were employed. A protein representative of the ubiquitous globin fold was chosen as the model system. We found that, while most of the detectable α-helical population in the unfolded ensemble does not depend on the presence of the C-terminal region (corresponding to the native G and H helices), specific N-to-C long-range contacts between the H and A-B-C regions enhance the helical secondary structure content of the N terminus (A-B-C regions). The simple approach introduced here, based on the evaluation of N-terminal polypeptide fragments of increasing length, is of general applicability to identify the influence of long-range interactions in unfolded proteins.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

How Amphipols Embed Membrane Proteins: Global Solvent Accessibility and Interaction with a Flexible Protein Terminus

Amphipathic polymers called amphipols provide a valuable alternative to detergents for keeping integral membrane proteins soluble in aqueous buffers. Here, we characterize spatial contacts of amphipol A8-35 with membrane proteins from two architectural classes: The 8-stranded β-barrel outer membrane protein OmpX and the α-helical protein bacteriorhodopsin. OmpX is well structured in A8-35, with its barrel adopting a fold closely similar to that in dihexanoylphosphocholine micelles. The accessibility of A8-35-trapped OmpX by a water-soluble paramagnetic molecule is highly similar to that in detergent micelles and resembles the accessibility in the natural membrane. For the α-helical protein bacteriorhodopsin, previously shown to keep its fold and function in amphipols, NMR data show that the imidazole protons of a polyhistidine tag at the N-terminus of the protein are exchange protected in the presence of detergent and lipid bilayer nanodiscs, but not in amphipols, indicating the absence of an interaction in the latter case. Overall, A8-35 exhibits protein interaction properties somewhat different from detergents and lipid bilayer nanodiscs, while maintaining the structure of solubilized integral membrane proteins.     

Depth of α-synuclein in a bilayer determined by fluorescence, neutron reflectometry, and computation

α-Synuclein (α-syn) membrane interactions are implicated in the pathogenesis of Parkinson's disease. Fluorescence and neutron reflectometry (NR) measurements reveal that α-syn penetrates ～9–14 Å into the outer leaflet of the bilayer, with a substantial portion of the membrane-bound polypeptide extending into the aqueous solvent. For the first time, to our knowledge, we used NR to obtain direct quantitative evidence of α-syn-induced membrane thinning. To examine the effect of specific residues on membrane penetration depths, we used a series of W4-containing N-terminal peptides. We identified that the first 15 residues (P15) nearly recapitulate the features of the full-length protein (i.e., partition constants, molecular mobility, and insertion of the W4 side chain into the bilayer), and found that as few as the first four N-terminal residues are sufficient for vesicle binding. Although at least one imperfect amphipathic repeat sequence (KAKEGV) is required for α-helical formation, secondary structural formation has little effect on membrane affinity. To develop an N-terminal α-syn model for bilayer interactions, we performed molecular-dynamics simulations of the P15 peptide submerged in a bilayer. The simulation results are highly consistent with experimental data indicating a broad low-energy region (8.5–14.5 Å) for W4 insertion.     

The Effect of Loops on the Structural Organization of α-Helical Membrane Proteins

Loops connecting the transmembrane (TM) α-helices in membrane proteins are expected to affect the structural organization of the thereby connected helices and the helical bundles as a whole. This effect, which has been largely ignored previously, is studied here by analyzing the x-ray structures of 41 α-helical membrane proteins. First we define the loop flexibility ratio, R, and find that 53% of the loops are stretched, where a stretched loop constrains the distance between the two connected helices. The significance of this constraining effect is supported by experiments carried out with bacteriorhodopsin and rhodopsin, in which cutting or eliminating their (predominately stretched) loops has led to a decrease in protein stability, and for rhodopsin, in most cases, also to the destruction of the structure. We show that for nonstretched loops in the extramembranous regions, the fraction of hydrophobic residues is comparable to that for soluble proteins; furthermore (as is also the case for soluble proteins), the hydrophobic residues in these regions are preferentially buried. This is expected to lead to the compact structural organization of the loops, which is transferred to the TM helices, causing them to assemble. We argue that a soluble protein complexed with a membrane protein similarly promotes compactness; other properties of such complexes are also studied. We calculate complementary attractive interactions between helices, including hydrogen bonds and van der Waals interactions of sequential motifs, such as GXXXG. The relative and combined effects of all these factors on the association of the TM helices are discussed and protein structures with only a few of these factors are analyzed. Our study emphasizes the need for classifying membrane proteins into groups according to structural organization. This classification should be considered when procedures for structural analysis or prediction are developed and applied. Detailed analysis of each structure is provided at http://flan.blm.cs.cmu.edu/memloop/     

Dodecyl Maltoside Protects Membrane Proteins In Vacuo

Molecular dynamics simulations have been used to characterize the effects of transfer from aqueous solution to a vacuum to inform our understanding of mass spectrometry of membrane-protein-detergent complexes. We compared two membrane protein architectures (an α-helical bundle versus a β-barrel) and two different detergent types (phosphocholines versus an alkyl sugar) with respect to protein stability and detergent packing. The β-barrel membrane protein remained stable as a protein-detergent complex in vacuum. Zwitterionic detergents formed conformationally destabilizing interactions with an α-helical membrane protein after detergent micelle inversion driven by dehydration in vacuum. In contrast, a nonionic alkyl sugar detergent resisted micelle inversion, maintaining the solution-phase conformation of the protein. This helps to explain the relative stability of membrane proteins in the presence of alkyl sugar detergents such as dodecyl maltoside.     

Peptide Nanopores and Lipid Bilayers: Interactions by Coarse-Grained Molecular-Dynamics Simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (α-helical and β-barrel), and the seven different bilayer systems range in thickness from ∼28 to ∼43 Å. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.     

In silico chaperonin-like cycle helps folding of proteins for structure prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-α-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.     

How Many 3D Structures Do We Need to Train a Predictor？

It has been shown that the progress in the determination of membrane protein structure grows exponentially, with approximately the same growth rate as that of the water-soluble proteins. In order to investigate the effect of this, on the performance of prediction algorithms for both α-helical and β-barrel membrane proteins, we conducted a prospective study based on historical records. We trained separate hidden Markov models with different sized training sets and evaluated their performance on topology pred...     

Molecular dynamics simulation of Bombolitin II in the dipalmitoylphosphatidylcholine membrane bilayer

The orientation behavior of Bombolitin II (BLT2) in the dipalmitoylphosphatidylcholine membrane bilayer was investigated by using molecular-dynamics simulation. During the 20-ns simulation, the BLT2 began to tilt and finally reached the angle of 51° from the membrane-normal. The structure of the peptide formed the amphipathic α-helical structure during the entire simulation time. The peptide tilts with its hydrophobic side faced to the hydrophobic core of the bilayer. We analyzed the mechanism of the tilting behavior of the peptide associated with the membrane in detail. The analysis showed that the hydrogen-bond interaction and the electrostatic interaction were found to exist between Lys¹² and a lipid molecule. These interactions are considered to work as an important factor in tilting the peptide to the membrane-normal.     

Contact prediction for beta and alpha‐beta proteins using integer linear optimization and its impact on the first principles 3D structure prediction method ASTRO‐FOLD

An integer linear optimization model is presented to predict residue contacts in β, α + β, and α/β proteins. The total energy of a protein is expressed as sum of a C^α? C^α distance dependent contact energy contribution and a hydrophobic contribution. The model selects contact that assign lowest energy to the protein structure as satisfying a set of constraints that are included to enforce certain physically observed topological information. A new method based on hydrophobicity is proposed to find the β‐sheet alignments. These β‐sheet alignments are used as constraints for contacts between residues of β‐sheets. This model was tested on three independent protein test sets and CASP8 test proteins consisting of β, α + β, α/β proteins and it was found to perform very well. The average accuracy of the predictions (separated by at least six residues) was ～61%. The average true positive and false positive distances were also calculated for each of the test sets and they are 7.58 Å and 15.88 Å, respectively. Residue contact prediction can be directly used to facilitate the protein tertiary structure prediction. This proposed residue contact prediction model is incorporated into the first principles protein tertiary structure prediction approach, ASTRO‐FOLD. The effectiveness of the contact prediction model was further demonstrated by the improvement in the quality of the protein structure ensemble generated using the predicted residue contacts for a test set of 10 proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

The Relation between α-Helical Conformation and Amyloidogenicity

Amyloid fibrils are stable aggregates of misfolded proteins and polypeptides that are insoluble and resistant to protease activity. Abnormal formation of amyloid fibrils in vivo may lead to neurodegenerative disorders and other systemic amyloidosis, such as Alzheimer’s, Parkinson’s, and atherosclerosis. Because of their clinical importance, amyloids are under intense scientific research. It is believed that short polypeptide segments within proteins are responsible for the transformation of correctly folded proteins into parts of larger amyloid fibrils and that this transition is modulated by environmental factors, such as pH, salt concentration, interaction with the cell membrane, and interaction with metal ions. Most studies on amyloids focus on the amyloidogenic sequences. The focus of this study is on the structure of the amyloidogenic α-helical segments because the α-helical secondary structure has been recognized to be a key player in different stages of the amyloidogenesis process. We have previously shown that the α-helical conformation may be expressed by two parameters (θ and ρ) that form orthogonal coordinates based on the Ramachandran dihedrals (φ and ψ) and provide an illuminating interpretation of the α-helical conformation. By performing statistical analysis on α-helical conformations found in the Protein Data Bank, an apparent relation between α-helical conformation, as expressed by θ and ρ, and amyloidogenicity is revealed. Remarkably, random amino acid sequences, whose helical structures were obtained from the most probable dihedral angles, revealed the same dependency of amyloidogenicity, suggesting the importance of α-helical structure as opposed to sequence.     

Novel approach for alpha-helical topology prediction in globular proteins: generation of interhelical restraints

The protein folding problem represents one of the most challenging problems in computational biology. Distance constraints and topology predictions can be highly useful for the folding problem in reducing the conformational space that must be searched by deterministic algorithms to find a protein structure of minimum conformational energy. We present a novel optimization framework for predicting topological contacts and generating interhelical distance restraints between hydrophobic residues in alpha-helical globular proteins. It should be emphasized that since the model does not make assumptions about the form of the helices, it is applicable to all alpha-helical proteins, including helices with kinks and irregular helices. This model aims at enhancing the ASTRO-FOLD protein folding approach of Klepeis and Floudas (Journal of Computational Chemistry 2003;24:191-208), which finds the structure of global minimum conformational energy via a constrained nonlinear optimization problem. The proposed topology prediction model was evaluated on 26 alpha-helical proteins ranging from 2 to 8 helices and 35 to 159 residues, and the best identified average interhelical distances corresponding to the predicted contacts fell below 11 A in all 26 of these systems. Given the positive results of applying the model to several protein systems, the importance of interhelical hydrophobic-to-hydrophobic contacts in determining the folding of alpha-helical globular proteins is highlighted.     

CASP11 – An Evaluation of a Modular BCL::Fold-Based Protein Structure Prediction Pipeline

In silico prediction of a protein’s tertiary structure remains an unsolved problem. The community-wide Critical Assessment of Protein Structure Prediction (CASP) experiment provides a double-blind study to evaluate improvements in protein structure prediction algorithms. We developed a protein structure prediction pipeline employing a three-stage approach, consisting of low-resolution topology search, high-resolution refinement, and molecular dynamics simulation to predict the tertiary structure of proteins from the primary structure alone or including distance restraints either from predicted residue-residue contacts, nuclear magnetic resonance (NMR) nuclear overhauser effect (NOE) experiments, or mass spectroscopy (MS) cross-linking (XL) data. The protein structure prediction pipeline was evaluated in the CASP11 experiment on twenty regular protein targets as well as thirty-three ‘assisted’ protein targets, which also had distance restraints available. Although the low-resolution topology search module was able to sample models with a global distance test total score (GDT_TS) value greater than 30% for twelve out of twenty proteins, frequently it was not possible to select the most accurate models for refinement, resulting in a general decay of model quality over the course of the prediction pipeline. In this study, we provide a detailed overall analysis, study one target protein in more detail as it travels through the protein structure prediction pipeline, and evaluate the impact of limited experimental data.     

Helix-helix interactions and their impact on protein motifs and assemblies

Protein secondary structure elements are arranged in distinct structural motifs such as four-α-helix bundle, 8α/8β TIM-barrel, Rossmann dinucleotide binding fold, assembly of a helical rod. Each structural motif is characterized by a particular type of helix-helix interactions. A unique pattern of contacts is formed by interacting helices of the structural motif. In each type of fold, edges of the helix surface, which participate in the formation of helix-helix contacts with preceding and following helices, differ. This work shows that circular arrangements of the four, eight, and sixteen α-helices, which are found in the four-α-helical motif, TIM-barrel 8α/8β fold, and helical rod of 16.3¯ helices per turn correspondingly, can be associated with the mutual positioning of the edges of the helix surfaces. Edges (i, i+1)−(i+1, i+2) of the helix surface are central for the interhelical contacts in a four-α-helix bundle. Edges (i, i+1)−(i+2, i+3) are involved in the assembly of four-α-helix subunits into helical rod of a tobacco mosaic virus and a three-helix fragment of a Rossmann fold. In 8α/8β TIM-barrel fold, edges (i, i+1)−(i+5, i+6) are involved in the octagon arrangement. Approximation of a cross section of each motif with a polygon (n-gon, n=4, 8, 16) shows that a good correlation exists between polygon interior angles and angles formed by the edges of helix surfaces.     

Exploring peptide-membrane interactions with coarse-grained MD simulations

The interaction of α-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an α-helix within a membrane, or that change the probability that the α-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an α-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.     

Prediction of helix-helix contacts and interacting helices in polytopic membrane proteins using neural networks

Despite rapidly increasing numbers of available 3D structures, membrane proteins still account for less than 1% of all structures in the Protein Data Bank. Recent high-resolution structures indicate a clearly broader structural diversity of membrane proteins than initially anticipated, motivating the development of reliable structure prediction methods specifically tailored for this class of molecules. One important prediction target capturing all major aspects of a protein's 3D structure is its contact map. Our analysis shows that computational methods trained to predict residue contacts in globular proteins perform poorly when applied to membrane proteins. We have recently published a method to identify interacting alpha-helices in membrane proteins based on the analysis of coevolving residues in predicted transmembrane regions. Here, we present a substantially improved algorithm for the same problem, which uses a newly developed neural network approach to predict helix-helix contacts. In addition to the input features commonly used for contact prediction of soluble proteins, such as windowed residue profiles and residue distance in the sequence, our network also incorporates features that apply to membrane proteins only, such as residue position within the transmembrane segment and its orientation toward the lipophilic environment. The obtained neural network can predict contacts between residues in transmembrane segments with nearly 26% accuracy. It is therefore the first published contact predictor developed specifically for membrane proteins performing with equal accuracy to state-of-the-art contact predictors available for soluble proteins. The predicted helix-helix contacts were employed in a second step to identify interacting helices. For our dataset consisting of 62 membrane proteins of solved structure, we gained an accuracy of 78.1%. Because the reliable prediction of helix interaction patterns is an important step in the classification and prediction of membrane protein folds, our method will be a helpful tool in compiling a structural census of membrane proteins.     

Probing the Transmembrane Structure and Dynamics of Microsomal NADPH-cytochrome P450 oxidoreductase by Solid-State NMR

NADPH-cytochrome P450 oxidoreductase (CYPOR) is an essential redox partner of the cytochrome P450 (cyt P450) superfamily of metabolic enzymes. In the endoplasmic reticulum of liver cells, such enzymes metabolize ∼75% of the pharmaceuticals in use today. It is known that the transmembrane domain of CYPOR plays a crucial role in aiding the formation of a complex between CYPOR and cyt P450. Here we present the transmembrane structure, topology, and dynamics of the FMN binding domain of CYPOR in a native membrane-like environment. Our solid-state NMR results reveal that the N-terminal transmembrane domain of CYPOR adopts an α-helical conformation in the lipid membrane environment. Most notably, we also show that the transmembrane helix is tilted ∼13° from the lipid bilayer normal, and exhibits motions on a submillisecond timescale including rotational diffusion of the whole helix and fluctuation of the helical director axis. The approaches and the information reported in this study would enable further investigations on the structure and dynamics of the full-length NADPH-cytochrome P450 oxidoreductase and its interaction with other membrane proteins in a membrane environment.