Comparing native and irradiated E. coli lactose repressor-operator complex by molecular dynamics simulation.

The function of the E. coli lactose operon requires the binding of the tetrameric repressor protein to the operator DNA. We have previously shown that gamma-irradiation destabilises the repressor-operator complex because the repressor gradually loses its DNA-binding ability (Radiat Res 170:604-612, 2008). It was suggested that the observed oxidation of tyrosine residues and the concomitant structural changes of irradiated headpieces (DNA-binding domains of repressor monomers) could be responsible for the inactivation. To unravel the mechanisms that lead to repressor-operator complex destabilisation when tyrosine oxidation occurs, we have compared by molecular dynamic simulations two complexes: (1) the native complex formed by two headpieces and the operator DNA, and (2) the damaged complex, in which all tyrosines are replaced by their oxidation product 3,4-dihydroxyphenylalanine (DOPA). On a 20 ns time scale, MD results show effects consistent with complex destabilisation: increased flexibility, increased DNA bending, modification of the hydrogen bond network, and decrease of the positive electrostatic potential at the protein surface and of the global energy of DNA-protein interactions.

Transmembrane helix packing of ErbB/Neu receptor in membrane environment: a molecular dynamics study.

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.

Inhibition of PlexA1-mediated brain tumor growth and tumor-associated angiogenesis using a transmembrane domain targeting peptide.

The neuropilin-plexin receptor complex regulates tumor cell migration and proliferation and thus is an interesting therapeutic target. High expression of neuropilin-1 is indeed associated with a bad prognosis in glioma patients. Q-RTPCR and tissue-array analyses showed here that Plexin-A1 is highly expressed in glioblastoma and that the highest level of expression correlates with the worse survival of patients. We next identified a developmental and tumor-associated pro-angiogenic role of Plexin-A1. Hence, by using molecular simulations and a two-hybrid like assay in parallel with biochemical and cellular assays we developed a specific Plexin-A1 peptidic antagonist disrupting transmembrane domain-mediated oligomerization of the receptor and subsequent signaling and functional activity. We found that this peptide exhibits anti-tumor activity in vivo on different human glioblastoma models including glioma cancer stem cells. Thus, screening Plexin-A1 expression and targeting Plexin-A1 in glioblastoma patients exhibit diagnostic and therapeutic value.

Transmembrane peptides from tyrosine kinase receptor. Mutation-related behavior in a lipid bilayer investigated by molecular dynamics simulations.

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer.

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.

Transmembrane recognition of the semaphorin co-receptors neuropilin 1 and plexin A1: coarse-grained simulations.

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.

Transmembrane domain targeting peptide antagonizing ErbB2/Neu inhibits breast tumor growth and metastasis.

Breast cancer is still a deadly disease despite major achievements in targeted therapies designed to block ligands or ligand-binding subunits of major tyrosine kinase receptors. Relapse is significant and metastases deleterious, which demands novel strategies for fighting this disease. Here, we report a proof-of-concept experiment demonstrating that small peptides interfering with the transmembrane domain of the tyrosine kinase epidermal growth factor receptor ErbB2 exhibit anticancer properties when used at micromolar dosages in a genetically engineered mouse model of breast cancer. Different assays demonstrate the specificity of the ErbB2-targeting peptide, which induces long-term reduction of ErbB2 phosphorylation and Akt signaling consistent with reduced tumor cell proliferation and increased survival. Microcomputed tomography analysis established the antimetastatic activity of the peptide and its impact on primary tumor growth. This reveals the interior of the cell membrane as an unexplored dimension for drug design.

Ligand entry pathways in the ligand binding domain of PPARγ receptor.

To address the question of ligand entry process, we report targeted molecular dynamics simulations of the entry of the flexible ionic ligand GW0072 in the ligand binding domain of the nuclear receptor PPARγ. Starting with the ligand outside the receptor the simulations led to a ligand docked inside the binding pocket resulting in a structure very close to the holo-form of the complex. The results showed that entry process is guided by hydrophobic interactions and that entry pathways are very similar to exit pathways. We suggest that TMD method may help in discriminating between ligands generated by in silico docking.

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.

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

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer.

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.