Helix Interactions in Membranes: Lessons from Unrestrained Monte Carlo Simulations

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

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

A helix heterodimer in a lipid bilayer: prediction of the structure of an integrin transmembrane domain via multiscale simulations

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling.

Structure elucidation of dimeric transmembrane domains of bitopic proteins

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

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

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

[Prediction of the spatial structure of proteins: emphasis on membrane targets]

Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.

Method To Assess Packing Quality of Transmembrane α-Helices in Proteins. 1. Parametrization Using Structural Data

Integral membrane proteins (MPs) are pharmaceutical targets of exceptional importance. Modern methods of three-dimensional protein structure determination often fail to supply the fast growing field of structure-based drug design with the requested MPs' structures. That is why computational modeling techniques gain a special importance for these objects. Among the principal difficulties limiting application of these methods is the low quality of the MPs' models built in silico. In this series of two papers we present a computational approach to the assessment of the packing "quality" of transmembrane (TM) alpha-helical domains in proteins. The method is based on the concept of protein environment classes, whereby each amino acid residue is described in terms of its environment polarity and accessibility to the membrane. In the first paper we analyze a nonredundant set of 26 TM alpha-helical domains and compute the residues' propensities to five predefined classes of membrane-protein environments. Here we evaluate the proposed approach only by various test sets, cross-validation protocols and ability of the method to delimit the crystal structure of visual rhodopsin, and a number of its erroneous theoretical models. More advanced validation of the method is given in the second article of this series. We assume that the developed "membrane score" method will be helpful in optimizing computer models of TM domains of MPs, especially G-protein coupled receptors.

A novel method for packing quality assessment of transmembrane alpha-helical domains in proteins

Here we present a novel method for assessment of packing quality for transmembrane (TM) domains of alpha-helical membrane proteins (MPs), based on analysis of available high-resolution experimental structures of MPs. The presented concept of protein-membrane environment classes permits quantitative description of packing characteristics in terms of membrane accessibility and polarity of the nearest protein groups. We demonstrate that the method allows identification of native-like conformations among the large set of theoretical MP models. The developed "membrane scoring function" will be of use for optimization of TM domain packing in theoretical models of MPs, first of all G-protein coupled receptors.

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

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