Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

E(z), a depth-dependent potential for assessing the energies of insertion of amino acid side-chains into membranes: derivation and applications to determining the orientation of transmembrane and interfacial helices

We have developed an empirical residue-based potential (E(z) potential) for protein insertion in lipid membranes. Propensities for occurrence as a function of depth in the bilayer were calculated for the individual amino acid types from their distribution in known structures of helical membrane proteins. The propensities were then fit to continuous curves and converted to a potential using a reverse-Boltzman relationship. The E(z) potential demonstrated a good correlation with experimental data such as amino acid transfer free energy scales (water to membrane center and water to interface), and it incorporates transmembrane helices of varying composition in the membrane with trends similar to those obtained with translocon-mediated insertion experiments. The potential has a variety of applications in the analysis of natural membrane proteins as well as in the design of new ones. It can help in calculating the propensity of single helices to insert in the bilayer and estimate their tilt angle with respect to the bilayer normal. It can be utilized to discriminate amphiphilic helices that assume a parallel orientation at the membrane interface, such as those of membrane-active peptides. In membrane protein design applications, the potential allows an environment-dependent selection of amino acid identities.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

Marginally hydrophobic transmembrane α‐helices shaping membrane protein folding

Cells have developed an incredible machinery to facilitate the insertion of membrane proteins into the membrane. While we have a fairly good understanding of the mechanism and determinants of membrane integration, more data is needed to understand the insertion of membrane proteins with more complex insertion and folding pathways. This review will focus on marginally hydrophobic transmembrane helices and their influence on membrane protein folding. These weakly hydrophobic transmembrane segments are by themselves not recognized by the translocon and therefore rely on local sequence context for membrane integration. How can such segments reside within the membrane? We will discuss this in the light of features found in the protein itself as well as the environment it resides in. Several characteristics in proteins have been described to influence the insertion of marginally hydrophobic helices. Additionally, the influence of biological membranes is significant. To begin with, the actual cost for having polar groups within the membrane may not be as high as expected; the presence of proteins in the membrane as well as characteristics of some amino acids may enable a transmembrane helix to harbor a charged residue. The lipid environment has also been shown to directly influence the topology as well as membrane boundaries of transmembrane helices—implying a dynamic relationship between membrane proteins and their environment.     

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.     

How strongly do sequence conservation patterns and empirical scales correlate with exposure patterns of transmembrane helices of membrane proteins?

Given the difficulty in determining high-resolution structures of helical membrane proteins, sequence-based prediction methods can be useful in elucidating diverse physiological processes mediated by this important class of proteins. Predicting the angular orientations of transmembrane (TM) helices about the helix axes, based on the helix parameters from electron microscopy data, is a classical problem in this regard. This problem has triggered the development of a number of different empirical scales. Recently, sequence conservation patterns were also made use of for improved predictions. Empirical scales and sequence conservation patterns (collectively termed as "prediction scales") have also found frequent applications in other research areas of membrane proteins: for example, in structure modeling and in prediction of buried TM helices. This trend is expected to grow in the near future unless there are revolutionary developments in the experimental characterization of membrane proteins. Thus, it is timely and imperative to carry out a comprehensive benchmark test over the prediction scales proposed so far to determine their pros and cons. In the current analysis, we use exposure patterns of TM helices as a golden standard, because if one develops a prediction scale that correlates perfectly with exposure patterns of TM helices, it will enable one to predict buried residues (or buried faces) of TM helices with an accuracy of 100%. Our analysis reveals several important points. (1) It demonstrates that sequence conservation patterns are much more strongly correlated with exposure patterns of TM helices than empirical scales. (2) Scales that were specifically parameterized using structure data (structure-based scales) display stronger correlation than hydrophobicity-based scales, as expected. (3) A nonnegligible difference is observed among the structure-based scales in their correlational property, suggesting that not every learning algorithm is equally effective. (4) A straightforward framework of optimally combining sequence conservation patterns and empirical scales is proposed, which reveals that improvements gained from combining the two sources of information are not dramatic in almost all cases. In turn, this calls for the development of fundamentally different scales that capture the essentials of membrane protein folding for substantial improvements.     

Understanding peptide interactions with the lipid bilayer: a guide to membrane protein engineering

In recent years, studies of the interactions of designed hydrophobic and amphipathic polypeptides with biological membranes have progressed considerably. In-plane and transmembrane helical domains have been engineered as well as sequences that exhibit dynamic distributions of different topologies. These sequences not only help our understanding of the thermodynamic interaction contributions that determine peptide orientation, but also exhibit interesting biological functions as autonomous units. In addition, helices are considered to be folding intermediates of multi-spanning membrane proteins. The specificity and thermodynamic stability of transmembrane helix-helix interactions or the assembly of beta sheets at the membrane surface have, therefore, become a focus of ongoing research.     

A highly unusual palindromic transmembrane helical hairpin formed by SARS coronavirus E protein

The agent responsible for the recent severe acute respiratory syndrome (SARS) outbreak is a previously unidentified coronavirus. While there is a wealth of epidemiological studies, little if any molecular characterization of SARS coronavirus (SCoV) proteins has been carried out. Here we describe the molecular characterization of SCoV E protein, a critical component of the virus responsible for virion envelope morphogenesis. We conclusively show that SCoV E protein contains an unusually short, palindromic transmembrane helical hairpin around a previously unidentified pseudo-center of symmetry, a structural feature which seems to be unique to SCoV. The hairpin deforms lipid bilayers by way of increasing their curvature, providing for the first time a molecular explanation of E protein's pivotal role in viral budding. The molecular understanding of this critical component of SCoV may represent the beginning of a concerted effort aimed at inhibiting its function, and consequently, viral infectivity.     

Molecular modeling of c-erbB2 receptor dimerization: Coiled-coil structure of wild and oncogenic transmembrane domains—Stabilization by interhelical hydrogen bonds in the oncogenic form

Dimerization models of c-erbB2 transmembrane domains (Leu651-Ile675) are studied by molecular mechanics and molecular dynamics simulations. Both wild and Glu mutated transmembrane helices exhibit the same relative orientation for favorable associations and dimerize preferentially in left-handed coiled-coil structures. The mutation point 659 belongs to the interfacing residues, and in the transforming domain, symmetric hydrogen bonds between Glu carboxylic groups stabilize the dimeric structure. The same helix packing found for the wild dimers, except side-chain—side-chain hydrogen bonds, suggests that the transmembrane domains dimerize according to similar process. Structural and energetical characterization of the models are presented. © 1997 John Wiley & Sons, Inc. Biopoly 42: 157–168, 1997     

Modeling of peptides and proteins in a membrane environment: II. Structural and energetic aspects of glycophorin A in a lipid bilayer

The conformational space of a hydrophobic peptide fragment of glycophorin A in a lipid membrane was studied with the Monte Carlo method using the solvation model described in the first communication of this series. The simulation was performed for various starting orientations of the peptide relative the membrane bilayer: outside, inside, partially immersed, and transbilayer. We showed that the membrane substantially stabilizes the α-helical conformation of the central hydrophobic part of the glycophorin A molecule, which for the most part is immersed in the apolar core of the bilayer. For various conformational states, energy values were calculated and the orientations of the peptide relative to the membrane were characterized. Depending on the thickness of the bilayer, either an entirely α-helical conformation in transbilayer orientation or a conformation with a kink in the central part of the helix with theN- andC-termini exposed on one side of the membrane corresponds to the minimal-energy structure. The transmembrane orientation of glycophorin A is energetically advantageous when the membrane thickness is close to the length of its hydrophobic helical portion, which is consistent with the effect ofhydrophobic match observed experimentally. The prospects for further refinement of the model are discussed. For communication I, see [1].