Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: difference between atomistic and coarse-grained simulations

Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT(OPLS)). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT(OPLS) PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT(CHARMM)) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16-22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT(OPLS) PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA)3 peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.

Prediction, refinement, and persistency of transmembrane helix dimers in lipid bilayers using implicit and explicit solvent/lipid representations: microsecond molecular dynamics simulations of ErbB1/B2 and EphA1

All-atom simulations are carried out on ErbB1/B2 and EphA1 transmembrane helix dimers in lipid bilayers starting from their solution/DMPC bicelle NMR structures. Over the course of microsecond trajectories, the structures remain in close proximity to the initial configuration and satisfy the majority of experimental tertiary contact restraints. These results further validate CHARMM protein/lipid force fields and simulation protocols on Anton. Separately, dimer conformations are generated using replica exchange in conjunction with an implicit solvent and lipid representation. The implicit model requires further improvement, and this study investigates whether lengthy all-atom molecular dynamics simulations can alleviate the shortcomings of the initial conditions. The simulations correct many of the deficiencies. For example, excessive helix twisting is eliminated over a period of hundreds of nanoseconds. The helix tilt, crossing angles, and dimer contacts approximate those of the NMR-derived structure, although the detailed contact surface remains off-set for one of two helices in both systems. Hence, even microsecond simulations are not long enough for extensive helix rotations. The alternate structures can be rationalized with reference to interaction motifs and may represent still sought after receptor states that are important in ErbB1/B2 and EphA1 signaling.

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.

Transmembrane domains interactions within the membrane milieu: principles, advances and challenges

Protein-protein interactions within the membrane are involved in many vital cellular processes. Consequently, deficient oligomerization is associated with known diseases. The interactions can be partially or fully mediated by transmembrane domains (TMD). However, in contrast to soluble regions, our knowledge of the factors that control oligomerization and recognition between the membrane-embedded domains is very limited. Due to the unique chemical and physical properties of the membrane environment, rules that apply to interactions between soluble segments are not necessarily valid within the membrane. This review summarizes our knowledge on the sequences mediating TMD-TMD interactions which include conserved motifs such as the GxxxG, QxxS, glycine and leucine zippers, and others. The review discusses the specific role of polar, charged and aromatic amino acids in the interface of the interacting TMD helices. Strategies to determine the strength, dynamics and specificities of these interactions by experimental (ToxR, TOXCAT, GALLEX and FRET) or various computational approaches (molecular dynamic simulation and bioinformatics) are summarized. Importantly, the contribution of the membrane environment to the TMD-TMD interaction is also presented. Studies utilizing exogenously added TMD peptides have been shown to influence in vivo the dimerization of intact membrane proteins involved in various diseases. The chirality independent TMD-TMD interactions allows for the design of novel short d- and l-amino acids containing TMD peptides with advanced properties. Overall these studies shed light on the role of specific amino acids in mediating the assembly of the TMDs within the membrane environment and their contribution to protein function. This article is part of a Special Issue entitled: Protein Folding in Membranes.

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.

Mechanistic and signaling analysis of Muc4-ErbB2 signaling module: new insights into the mechanism of ligand-independent ErbB2 activity

The membrane mucin Muc4 is aberrantly expressed in numerous epithelial carcinomas and is currently used as a cancer diagnostic and prognostic tool. Muc4 can also potentiate signal transduction by modulating differential ErbB2 phosphorylation in the absence and in the presence of the ErbB3 soluble ligand heregulin (HRG-beta1). These features of Muc4 suggest that Muc4 is not merely a cancer marker, but an oncogenic factor with a unique-binding/activation relationship with the receptor ErbB2. In the present study, we examined the signaling mechanisms that are associated with the Muc4-ErbB2 module by analyzing ErbB2 differential signaling in response to Muc4 expression. Our study was carried out in the A375 human melanoma and BT-474 breast cancer cell lines as our model systems. Quantitative and comparative signaling modulations were evaluated by immunoblot using phospho-specific antibodies, and densitometry analysis. Signaling complex components were identified by chemical cross-linking, fractionation by gel filtration, immunoprecipitation, and immunoblotting. Activated downstream signaling pathways were analyzed by an antibody microarray screen and immunoblot analyses. Our results indicate that Muc4 modulates ErbB2 signaling potential significantly by stabilizing and directly interacting with the ErbB2-ErbB3 heterodimer. Further analyses indicate that Muc4 promotes ErbB2 autocatalysis, but it has no effect on ErbB3 phosphorylation, although the chemical cross-linking data indicated that the signaling module is composed of Muc4, ErbB2, and ErbB3. Our microarray analysis indicates that Muc4 expression promotes cell migration by increasing the phosphorylation of the focal adhesion kinase and also through an increase in the levels of beta-catenin.

Single-spanning transmembrane domains in cell growth and cell-cell interactions: More than meets the eye?

As a whole, integral membrane proteins represent about one third of sequenced genomes, and more than 50% of currently available drugs target membrane proteins, often cell surface receptors. Some membrane protein classes, with a defined number of transmembrane (TM) helices, are receiving much attention because of their great functional and pharmacological importance, such as G protein-coupled receptors possessing 7 TM segments. Although they represent roughly half of all membrane proteins, bitopic proteins (with only 1 TM helix) have so far been less well characterized. Though they include many essential families of receptors, such as adhesion molecules and receptor tyrosine kinases, many of which are excellent targets for biopharmaceuticals (peptides, antibodies, et al.). A growing body of evidence suggests a major role for interactions between TM domains of these receptors in signaling, through homo and heteromeric associations, conformational changes, assembly of signaling platforms, etc. Significantly, mutations within single domains are frequent in human disease, such as cancer or developmental disorders. This review attempts to give an overview of current knowledge about these interactions, from structural data to therapeutic perspectives, focusing on bitopic proteins involved in cell signaling.

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

Essential role for the putative S6 inner pore region in the activation gating of the human TRPA1 channel

The ankyrin transient receptor potential channel TRPA1 is a sensory neuron-specific channel that is gated by various proalgesic agents such as allyl isothiocyanate (AITC), deep cooling or highly depolarizing voltages. How these disparate stimuli converge on the channel protein to open/close its ion-conducting pore is unknown. We identify several residues within the S6 inner pore-forming region of human TRPA1 that contribute to AITC and voltage-dependent gating. Alanine substitution in the conserved mid-S6 proline (P949A) strongly affected the activation/deactivation and ion permeation. The P949A was functionally restored by substitution with a glycine but not by the introduction of a proline at positions -1, -2 or +1, which indicates that P949 is structurally required for the normal functioning of the TRPA1 channel. Mutation N954A generated a constitutively open phenotype, suggesting a role in stabilizing the closed conformation. Alanine substitutions in the distal GXXXG motif decreased the relative permeability of the channel for Ca(2+) and strongly affected its activation/deactivation properties, indicating that the distal G962 stabilizes the open conformation. G958, on the other hand, provides additional tuning leading to decreased channel activity. Together these findings provide functional support for the critical role of the putative inner pore region in controlling the conformational changes that determine the transitions between the open and close states of the TRPA1 channel.

Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils

The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.

Insight into the recognition patterns of the ErbB receptor family transmembrane domains: heterodimerization models through molecular dynamics search

ErbB receptors undergo a complex interaction network defining hierarchical and competition relationships. Dimerization is driven entirely by receptor-receptor interactions and the transmembrane domains play a role in modulating the specificity and the selection of the partners during signal transduction. To shed light on the role of the GxxxG-like dimerization motifs in the formation of ErbB transmembrane heterodimers, we propose structural models resulting from conformational search method combined with molecular dynamics simulations. Left-handed structures of the transmembrane heterodimers are found preponderant over right-handed structures. All heterotypic heterodimers undergo two modes of association either via the N-terminal motif or the C-terminal motif. The transmembrane domain of ErbB3 impairs this C-terminal motif but also associates with the other partners owing to the presence of Gly residues. The two dimerization modes involve different orientations of the two helices. Thus, a molecular-switch model allowing the transition between the two dimerizing states may apply to the heterodimers and could help interpret receptor competition for the formation of homodimers and heterodimers. The comparison between experimental and theoretical results on the dimerization hierarchy of the transmembrane domains is not straightforward. However, we demonstrate that the intrinsic properties of the transmembrane sequences are an important component in heterodimer formation and that the ErbB2 and ErbB3 transmembrane domains have a strong power for heterodimerization as observed experimentally.