Comparison of the transmembrane helices of bovine rhodopsin in the crystal structure and the C α template based on cryo-electron microscopy maps and sequence analysis of the G protein-coupled receptors

G protein-coupled receptors (GPCR) are activated by a diverse array of extracellular signals, ranging from light to polypeptide molecules. The receptors propagate these signals intracellularly using G protein secondary messenger pathways. A common feature in the architecture of these receptors is their seven transmembrane domains. The first crystal structure of a GPCR, bovine rhodopsin, has recently been solved at 2.8 Å. We compared the seven membrane-spanning helices (TMH) from the crystal structure of bovine rhodopsin with those from the low-resolution model of bovine rhodopsin based on the cryo-electron microscopy structure of frog rhodopsin developed by Dr Joyce Baldwin. The model developed by Baldwin used a consensus sequence approach to predict the rotational position of each helix with respect to the other six helices. Superposition of the entire helix bundle of the Baldwin model with the crystal structure gave a RMS difference (RMSD) of 3.2 Å for the 198 C f atoms which suggests a high level of similarity in the arrangement of the helices. Except for TMH IV (RMSD of 4.0 Å), the position of corresponding helices within the helix bundle overlapped well. The superposition of individual helices showed that the RMSD values over 3 Å in the global superposition were largely due to one or more of the following: (i) differences in the unraveling and kinks for these helices, (ii) translation of TMH perpendicular to the membrane and (iii) rotation of helices up to 31°, except for TMH IV in which an additional contribution to the RMSD came from the aforementioned observation. As other crystal structures of GPCRs become available, a comparison with the Baldwin consensus model may reveal larger differences than those observed here.     

Simulation of the packing of idealized transmembrane α-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala₂₀ helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and δ-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials. Received: 16 November 1998 / Revised version: 26 March 1999 / Accepted: 8 April 1999     

An automatic method for predicting transmembrane protein structures using cryo-EM and evolutionary data

The transmembrane (TM) domains of many integral membrane proteins are composed of alpha-helix bundles. Structure determination at high resolution (<4 A) of TM domains is still exceedingly difficult experimentally. Hence, some TM-protein structures have only been solved at intermediate (5-10 A) or low (>10 A) resolutions using, for example, cryo-electron microscopy (cryo-EM). These structures reveal the packing arrangement of the TM domain, but cannot be used to determine the positions of individual amino acids. The observation that typically, the lipid-exposed faces of TM proteins are evolutionarily more variable and less charged than their core provides a simple rule for orienting their constituent helices. Based on this rule, we developed score functions and automated methods for orienting TM helices, for which locations and tilt angles have been determined using, e.g., cryo-EM data. The method was parameterized with the aim of retrieving the native structure of bacteriorhodopsin among near- and far-from-native templates. It was then tested on proteins that differ from bacteriorhodopsin in their sequences, architectures, and functions, such as the acetylcholine receptor and rhodopsin. The predicted structures were within 1.5-3.5 A from the native state in all cases. We conclude that the computational method can be used in conjunction with cryo-EM data to obtain approximate model structures of TM domains of proteins for which a sufficiently heterogeneous set of homologs is available. We also show that in those proteins in which relatively short loops connect neighboring helices, the scoring functions can discriminate between near- and far-from-native conformations even without the constraints imposed on helix locations and tilt angles that are derived from cryo-EM.     

FoldGPCR: Structure prediction protocol for the transmembrane domain of G protein‐coupled receptors from class A

Building reliable structural models of G protein‐coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template‐based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1–3 Å Cα RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter‐residue contacts between the template and target, and exploits an all‐atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor‐ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human β₂‐adrenergic receptor (β₂AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X‐ray crystal structure of the β₂AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix‐helix packing interactions that distinguish rhodopsin from β₂AR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Thermodynamic and kinetic stabilities of transmembrane helix bundles as revealed by single-pair FRET analysis: Effects of the number of membrane-spanning segments and cholesterol

The tertiary structures and conformational dynamics of transmembrane (TM) helical proteins are maintained by the interhelical interaction network in membranes, although it is complicated to analyze the underlying driving forces because the amino acid sequences can involve multiple and various types of interactions. To obtain insights into basal and common effects of the number of membrane-spanning segments and membrane cholesterol, we measured stabilities of helix bundles composed of simple TM helices (AALALAA)₃ (1TM) and (AALALAA)₃-G₅-(AALALAA)₃ (2TM). Association–dissociation dynamics for 1TM–1TM, 1TM–2TM, and 2TM–2TM pairs were monitored to compare stabilities of 2-, 3-, and 4-helical bundles, respectively, with single-pair fluorescence resonance energy transfer (sp-FRET) in liposome membranes. Both thermodynamic and kinetic stabilities of the helix bundles increased with a greater number of membrane-spanning segments in POPC. The presence of 30 mol% cholesterol strongly enhanced the formation of 1TM–1TM and 1TM–2TM bundles (~ ? 9 kJ mol^?1), whereas it only weakly stabilized the 2TM–2TM bundle (~ ? 3 kJ mol^?1). Fourier transform infrared-polarized attenuated total reflection (ATR-FTIR) spectroscopy revealed an ~30° tilt of the helix axis relative to bilayer normal for the 1TM–2TM pair in the presence of cholesterol, suggesting the formation of a tilted helix bundle to release high lateral pressure at the center of cholesterol-containing membranes. These results demonstrate that the number of membrane-spanning segments affects the stability and structure of the helix bundle, and their cholesterol-dependences. Such information is useful to understand the basics of folding and assembly of multispanning TM proteins.     

On the spatial disposition of the fifth transmembrane helix and the structural integrity of the transmembrane binding site in the opioid and ORL1 G protein-coupled receptor family

Evidence from statistical cluster analyses of a multiple sequence alignment of G protein-coupled receptor seven-helix folds supports the existence of structurally conserved transmembrane (TM) ligand binding sites in the opioid/opioid receptor-like (ORL1) and amine receptor families. Based on the expectation that functionally conserved regions in homologous proteins will display locally higher levels of sequence identity compared with global sequence similarities that pertain to the overall fold, this approach may have wider applications in functional genomics to annotate sequence data. Binding sites in models of the kappa-opioid receptor seven-helix bundle built from the rhodopsin templates of Baldwin et al. (1997) [J. Mol. Biol., 272, 144-164] and Herzyk and Hubbard (1998) [J. Mol. Biol., 281, 742-751] are compared. The Herzyk and Hubbard template is found to be in better accord with experimental studies of amine, opioid and rhodopsin receptors owing to the reduced physical separation of the extracellular parts of TM helices V and VI and differences in the rotational orientation of the N-terminal of helix V that reveal side chain accessibilities in the Baldwin et al. structure to be out of phase with relative alkylation rates of engineered cysteine residues in the TM binding site of the alpha(2A)-adrenergic receptor. TM helix V in the Baldwin et al. template has been remodelled with a different proline kink to satisfy experimental constraints. A recent proposal that rotation of helix V is associated with receptor activation is critically discussed.     

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.     

Three-dimensional model for meta-II rhodopsin,an activated G-protein-coupled receptor

A novel approach that iteratively combined the results of energy calculations and experimental data was used to generate a three-dimensional (3D) model of the photoactivated state (R*) of bovine rhodopsin (Rh). The approach started with simplified energy calculations in an effort to find a set of sterically and energetically reasonable options for transmembrane (TM) helix arrangements with all-trans-retinal. Various 3D models of TM helix packing found by computations were then compared to limited site-directed spin-label experimental data regarding the transition of the TM helices of Rh in the inactive state (R) to those in the R* state to identify the most plausible model of the TM helical bundle. At the next step, all non-TM structural elements, such as the non-TM helix 8, the N- and C-terminal fragments, and the loops connecting TM helices, were reconstructed, and after the entire R* structure had been relaxed, all other currently available additional experimental data, both mutational and spectroscopic, on the structure of the meta-II state of rhodopsin were used to validate the resulting 3D model.     

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.     

Structures of the transmembrane helices of the G-protein coupled receptor, rhodopsin.

An hypothesis is tested that individual peptides corresponding to the transmembrane helices of the membrane protein, rhodopsin, would form helices in solution similar to those in the native protein. Peptides containing the sequences of helices 1, 4 and 5 of rhodopsin were synthesized. Two peptides, with overlapping sequences at their termini, were synthesized to cover each of the helices. The peptides from helix 1 and helix 4 were helical throughout most of their length. The N- and C-termini of all the peptides were disordered and proline caused opening of the helical structure in both helix 1 and helix 4. The peptides from helix 5 were helical in the middle segment of each peptide, with larger disordered regions in the N- and C-termini than for helices 1 and 4. These observations show that there is a strong helical propensity in the amino acid sequences corresponding to the transmembrane domain of this G-protein coupled receptor. In the case of the peptides from helix 4, it was possible to superimpose the structures of the overlapping sequences to produce a construct covering the whole of the sequence of helix 4 of rhodopsin. As similar superposition for the peptides from helix 1 also produced a construct, but somewhat less successfully because of the disordering in the region of sequence overlap. This latter problem was more severe for helix 5 and therefore a single peptide was synthesized for the entire sequence of this helix, and its structure determined. It proved to be helical throughout. Comparison of all these structures with the recent crystal structure of rhodopsin revealed that the peptide structures mimicked the structures seen in the whole protein. Thus similar studies of peptides may provide useful information on the secondary structure of other transmembrane proteins built around helical bundles.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Improving prediction of helix–helix packing in membrane proteins using predicted contact numbers as restraints

One of the challenging problems in tertiary structure prediction of helical membrane proteins (HMPs) is the determination of rotation of α‐helices around the helix normal. Incorrect prediction of helix rotations substantially disrupts native residue–residue contacts while inducing only a relatively small effect on the overall fold. We previously developed a method for predicting residue contact numbers (CNs), which measure the local packing density of residues within the protein tertiary structure. In this study, we tested the idea of incorporating predicted CNs as restraints to guide the sampling of helix rotation. For a benchmark set of 15 HMPs with simple to rather complicated folds, the average contact recovery (CR) of best‐sampled models was improved for all targets, the likelihood of sampling models with CR greater than 20% was increased for 13 targets, and the average RMSD100 of best‐sampled models was improved for 12 targets. This study demonstrated that explicit incorporation of CNs as restraints improves the prediction of helix–helix packing. Proteins 2017; 85:1212–1221. © 2017 Wiley Periodicals, Inc.     

High-resolution prediction of protein helix positions and orientations

We have developed a new method for predicting helix positions in globular proteins that is intended primarily for comparative modeling and other applications where high precision is required. Unlike helix packing algorithms designed for ab initio folding, we assume that knowledge is available about the qualitative placement of all helices. However, even among homologous proteins, the corresponding helices can demonstrate substantial differences in positions and orientations, and for this reason, improperly positioned helices can contribute significantly to the overall backbone root-mean-square deviation (RMSD) of comparative models. A helix packing algorithm for use in comparative modeling must obtain high precision to be useful, and for this reason we utilize an all-atom protein force field (OPLS) and a Generalized Born continuum solvent model. To reduce the computational expense associated with using a detailed, physics-based energy function, we have developed new hierarchical and multiscale algorithms for sampling the helices and flanking loops. We validate the method using a test suite of 33 cases, which are drawn from a diverse set of high-resolution crystal structures. The helix positions are reproduced with an average backbone RMSD of 0.6 A, while the average backbone RMSD of the complete loop-helix-loop region (i.e., the helix with the surrounding loops, which are also repredicted) is 1.3 A.     

Transmembrane helix dimerization: beyond the search for sequence motifs

Studies of the dimerization of transmembrane (TM) helices have been ongoing for many years now, and have provided clues to the fundamental principles behind membrane protein (MP) folding. Our understanding of TM helix dimerization has been dominated by the idea that sequence motifs, simple recognizable amino acid sequences that drive lateral interaction, can be used to explain and predict the lateral interactions between TM helices in membrane proteins. But as more and more unique interacting helices are characterized, it is becoming clear that the sequence motif paradigm is incomplete. Experimental evidence suggests that the search for sequence motifs, as mediators of TM helix dimerization, cannot solve the membrane protein folding problem alone. Here we review the current understanding in the field, as it has evolved from the paradigm of sequence motifs into a view in which the interactions between TM helices are much more complex. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Structure and function in rhodopsin: mapping light-dependent changes in distance between residue 316 in helix 8 and residues in the sequence 60-75, covering the cytoplasmic end of helices TM1 and TM2 and their connection loop CL1.

Double-spin-labeled mutants of rhodopsin were prepared containing a nitroxide side chain at position 316 in the cytoplasmic surface helix H8, and a second nitroxide in the sequence of residues 60-75, which includes the cytoplasmic loop CL1 and cytoplasmic ends of helices TM1 and TM2. Magnetic dipole-dipole interactions between the spins were analyzed to provide interspin distance distributions in both the dark and photoactivated states of rhodopsin. In the dark state in solutions of dodecyl maltoside, the interspin distances are found to be consistent with structural models of the nitroxide side chain and rhodopsin, both derived from crystallography. Photoactivation of rhodopsin shows a pattern of increases in internitroxide distance between the reference, position 316 in H8, and residues in CL1 and TM2 that suggests an outward displacement of TM2 relative to H8 by approximately 3 A.     

Comparative analysis of the heptahelical transmembrane bundles of G protein-coupled receptors

Background

G protein-coupled receptors represent a large family of eukaryotic membrane proteins, and are involved in almost all physiological processes in humans. Recent advances in the crystallographic study of these receptors enable a detailed comparative analysis of the commonly shared heptahelical transmembrane bundle. Systematic comparison of the bundles from a variety of receptors is indispensable for understanding not only of the structural diversification optimized for the binding of respective ligands but also of the structural conservation required for the common mechanism of activation accompanying the interaction changes among the seven helices.

Methodology/Principal Findings

We have examined the bundles of 94 polypeptide chains from almost all available structures of 11 receptors, which we classified into either inactivated chain or activated chain, based on the type of bound ligand. For the inactivated chains, superposition of 200 residue bundles by secondary structure matching demonstrated that the bound ligands share a laterally limited cavity in the extracellular section of the bundle. Furthermore, a distinct feature was found for helix III of bovine rhodopsin, which might have evolved to lower its activity in the presence of 11-cis-retinal, to a level that other receptors could hardly achieve with any currently available ligands.

Conclusions/Significance

Systematic analysis described here would be valuable for understanding of the rearrangement of seven helices which depends on the ligand specificity and activation state of the receptors.     

Conformational flexibility of the N-terminal domain of apolipoprotein a-I bound to spherical lipid particles

Lipid binding of human apolipoprotein A-I (apoA-I) occurs initially through the C-terminal alpha-helices followed by conformational reorganization of the N-terminal helix bundle. This led us to hypothesize that apoA-I has multiple lipid-bound conformations, in which the N-terminal helix bundle adopts either open or closed conformations anchored by the C-terminal domain. To investigate such possible conformations of apoA-I at the surface of a spherical lipid particle, site-specific labeling of the N- and C-terminal helices in apoA-I by N-(1-pyrene)maleimide was employed after substitution of a Cys residue for Val-53 or Phe-229. Neither mutagenesis nor the pyrene labeling caused discernible changes in the lipid-free structure and lipid interaction of apoA-I. Taking advantage of a significant increase in fluorescence when a pyrene-labeled helix is in contact with the lipid surface, we monitored the behaviors of the N- and C-terminal helices upon binding of apoA-I to egg PC small unilamellar vesicles. Comparison of the binding isotherms for pyrene-labeled apoA-I as well as a C-terminal helical peptide suggests that an increase in surface concentration of apoA-I causes dissociation of the N-terminal helix from the surface leaving the C-terminal helix attached. Consistent with this, isothermal titration calorimetry measurements showed that the enthalpy of apoA-I binding to the lipid surface under near saturated conditions is much less exothermic than that for binding at a low surface concentration, indicating the N-terminal helix bundle is out of contact with lipid at high apoA-I surface concentrations. Interestingly, the presence of cholesterol significantly induces the open conformation of the helix bundle. These results provide insight into the multiple lipid-bound conformations that the N-terminal helix bundle of apoA-I can adopt on a lipid or lipoprotein particle, depending upon the availability of space on the surface and the surface composition.     

Computational de novo design of a four-helix bundle protein—DND_4HB

The de novo design of proteins is a rigorous test of our understanding of the key determinants of protein structure. The helix bundle is an interesting de novo design model system due to the diverse topologies that can be generated from a few simple α-helices. Previously, noncomputational studies demonstrated that connecting amphipathic helices together with short loops can sometimes generate helix bundle proteins, regardless of the bundle''s exact sequence. However, using such methods, the precise positions of helices and side chains cannot be predetermined. Since protein function depends on exact positioning of residues, we examined if sequence design tools in the program Rosetta could be used to design a four-helix bundle with a predetermined structure. Helix position was specified using a folding procedure that constrained the design model to a defined topology, and iterative rounds of rotamer-based sequence design and backbone refinement were used to identify a low energy sequence for characterization. The designed protein, DND_4HB, unfolds cooperatively (T_m >90°C) and a NMR solution structure shows that it adopts the target helical bundle topology. Helices 2, 3, and 4 agree very closely with the design model (backbone RMSD = 1.11 Å) and >90% of the core side chain χ1 and χ2 angles are correctly predicted. Helix 1 lies in the target groove against the other helices, but is displaced 3 Å along the bundle axis. This result highlights the potential of computational design to create bundles with atomic-level precision, but also points at remaining challenges for achieving specific positioning between amphipathic helices.     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.