The optimisation of the helix/helix interaction of a transmembrane dimer is improved by the IMPALA restraint field

The interaction and orientation of Peptide KL4 in model membranes

We report on the orientation and location of synthetic pulmonary surfactant peptide KL₄, (KLLLL)₄K, in model lipid membranes. The partitioning depths of selectively deuterated leucine residues within KL₄ were determined in DPPC:POPG (4:1) and POPC:POPG (4:1) bilayers by oriented neutron diffraction. These measurements were combined with an NMR-generated model of the peptide structure to determine the orientation and partitioning of the peptide at the lipid-water interface. The results demonstrate KL₄ adopting an orientation that interacts with a single membrane leaflet. These observations are consistent with past ²H NMR and EPR studies (Antharam et al., 2009; Turner et al., 2014).

Multiscale modelling of the extracellular matrix

The extracellular matrix is a complex three-dimensional network of molecules that provides cells with a complex microenvironment. The major constituents of the extracellular matrix such as collagen, elastin and associated proteins form supramolecular assemblies contributing to its physicochemical properties and organization. The structure of proteins and their supramolecular assemblies such as fibrils have been studied at the atomic level (e.g., by X-ray crystallography, Nuclear Magnetic Resonance and cryo-Electron Microscopy) or at the microscopic scale. However, many protein complexes are too large to be studied at the atomic level and too small to be studied by microscopy. Most extracellular matrix components fall into this intermediate scale, so-called the mesoscopic scale, preventing their detailed characterization. Simulation and modelling are some of the few powerful and promising approaches that can deepen our understanding of mesoscale systems. We have developed a set of modelling tools to study the self-organization of the extracellular matrix and large motion of macromolecules at the mesoscale level by taking advantage of the dynamics of articulated rigid bodies as a mean to study a larger range of motions at the cost of atomic resolution.

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 new in-silico method for determination of helical transmembrane domains based on the PepLook scan: application to IL-2Rβ and IL-2Rγc receptor chains

BACKGROUND

Modeling of transmembrane domains (TMDs) requires correct prediction of interfacial residues for in-silico modeling and membrane insertion studies. This implies the defining of a target sequence long enough to contain interfacial residues. However, too long sequences induce artifactual polymorphism: within tested modeling methods, the longer the target sequence, the more variable the secondary structure, as though the procedure were stopped before the end of the calculation (which may in fact be unreachable). Moreover, delimitation of these TMDs can produce variable results with sequence based two-dimensional prediction methods, especially for sequences showing polymorphism. To solve this problem, we developed a new modeling procedure using the PepLook method. We scanned the sequences by modeling peptides from the target sequence with a window of 19 residues.

RESULTS

Using sequences whose NMR-structures are already known (GpA, EphA1 and Erb2-HER2), we first determined that the hydrophobic to hydrophilic accessible surface area ratio (ASAr) was the best criterion for delimiting the TMD sequence. The length of the helical structure and the Impala method further supported the determination of the TMD limits. This method was applied to the IL-2Rβ and IL-2Rγ TMD sequences of Homo sapiens, Rattus norvegicus, Mus musculus and Bos taurus.

CONCLUSIONS

We succeeded in reducing the variation in the TMD limits to only 2 residues and in gaining structural information.

[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.

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.

Energetics of the native and non-native states of the glycophorin transmembrane helix dimer

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures.

Prediction of membrane protein orientation in lipid bilayers: a theoretical approach

Over the past few years, several three-dimensional (3-D) structures of membrane proteins have been described with increasing accuracy, but their relationship with membranes are still not well understood. Recently, we have developed an empirical method, Integral Membrane Protein and Lipid Association (IMPALA), to predict the insertion of molecules (lipids, drugs) into lipid bilayers (Proteins 30 (1998) 357). The IMPALA uses a Monte Carlo minimisation procedure to calculate the depth and the angle of insertion of membrane-interacting molecules taking into account the restraints dictated by a lipid bilayer. In this paper, we use IMPALA to test the insertion of 23 integral membranous proteins (IMPs) and 2 soluble proteins into membranes. Four IMP are studied in detail: OmpA, maltoporin, MsCl channel and bacteriorhodopsin. The 3-D structures of the proteins are kept constant and the insertion into membrane is monitored by minimising the value of the restraint representing the sum of two terms, one for lipid perturbation and the other for hydrophobicity. The two soluble proteins are rejected from the membrane whereas, under the same conditions, all the membrane proteins remain inside, if the solvent accessible surface of the amino acids located inside the pore of porins is ignored. The results give the tilt angle of the IMP helices or strands with respect to the membrane surface and the depth of the protein mass centre insertion. We conclude that the restraint terms of IMPALA could be used to study the insertion of model structures or complexes of proteins within membranes.