Structural and thermodynamic insight into the process of “weak” dimerization of the ErbB4 transmembrane domain by solution NMR

Substrate Selection Criteria in Regulated Intramembrane Proteolysis

Alzheimer's disease is the most common form of dementia encountered in an aging population. Characteristic amyloid deposits of Aβ peptides in the brain are generated through cleavage of amyloid precursor protein (APP) by γ-secretase, an intramembrane protease. Cryo-EM structures of substrate γ-secretase complexes revealed details of the process, but how substrates are recognized and enter the catalytic site is still largely ignored. γ-Secretase cleaves a diverse range of substrate sequences without a common consensus sequence, but strikingly, single point mutations within the transmembrane domain (TMD) of specific substrates may greatly affect cleavage efficiencies. Previously, conformational flexibility was hypothesized to be the main criterion for substrate selection. Here we review the 3D structure and dynamics of several γ-secretase substrate TMDs and compare them with mutants shown to affect the cleavage efficiency. In addition, we present structural and dynamic data on ITGB1, a known nonsubstrate of γ-secretase. A comparison of biophysical details between these TMDs and changes generated by introducing crucial mutations allowed us to unravel common principles that differ between substrates and nonsubstrates. We identified three motifs in the investigated substrates: a highly flexible transmembrane domain, a destabilization of the cleavage region, and a basic signature at the end of the transmembrane helix. None of these appears to be exclusive. While conformational flexibility on its own may increase cleavage efficiency in well-known substrates like APP or Notch1, our data suggest that the three motifs seem to be rather variably combined to determine whether a transmembrane helix is efficiently recognized as a γ-secretase substrate.

Transmembrane dimers of type 1 receptors sample alternate configurations: MD simulations using coarse grain Martini 3 versus AlphaFold2 Multimer

Structures and dynamics of transmembrane (TM) receptor regions are key to understanding their signaling mechanism across membranes. Here we examine configurations of TM region dimers, assembled using the recent Martini 3 force field for coarse-grain (CG) molecular dynamics simulations. At first glance, our results show only a reasonable agreement with ab initio predictions using PREDDIMER and AlphaFold2 Multimer and with nuclear magnetic resonance (NMR)-derived structures. 5 of 11 CG TM structures are similar to the NMR structures (within <3.5 Å root-mean-square deviation [RMSD]) compared with 10 and 9 using PREDDIMER and AlphaFold2, respectively (with 8 structures of the later within 1.5 Å). Surprisingly, AlphaFold2 predictions are closer to NMR structures when the 2001 instead of 2020 database is used for training. The CG simulations reveal that alternative configurations of TM dimers readily interconvert with a predominant population. The implications for transmembrane signaling are discussed, including for the development of peptide-based pharmaceuticals.

Thermodynamic analysis of the GASright transmembrane motif supports energetic model of dimerization

The GASright motif, best known as the fold of the glycophorin A transmembrane dimer, is one of the most common dimerization motifs in membrane proteins, characterized by its hallmark GxxxG-like sequence motifs (GxxxG, AxxxG, GxxxS, and similar). Structurally, GASright displays a right-handed crossing angle and short interhelical distance. Contact between the helical backbones favors the formation of networks of weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O=C). To understand the factors that modulate the stability of GASright, we previously presented a computational and experimental structure-based analysis of 26 predicted dimers. We found that the contributions of van der Waals packing and Cα-H hydrogen bonding to stability, as inferred from the structural models, correlated well with relative dimerization propensities estimated experimentally with the in vivo assay TOXCAT. Here we test this model with a quantitative thermodynamic analysis. We used Förster resonance energy transfer (FRET) to determine the free energy of dimerization of a representative subset of seven of the 26 original TOXCAT dimers using FRET. To overcome the technical issue arising from limited sampling of the dimerization isotherm, we introduced a globally fitting strategy across a set of constructs comprising a wide range of stabilities. This strategy yielded precise thermodynamic data that show strikingly good agreement between the original propensities and ΔG° of association in detergent, suggesting that TOXCAT is a thermodynamically driven process. From the correlation between TOXCAT and thermodynamic stability, the predicted free energy for all the 26 GASright dimers was calculated. These energies correlate with the in silico ΔE scores of dimerization that were computed on the basis of their predicted structure. These findings corroborate our original model with quantitative thermodynamic evidence, strengthening the hypothesis that van der Waals and Cα-H hydrogen bond interactions are the key modulators of GASright stability.

Genetic interaction between GABRA1 and ERBB4 variants in the pathogenesis of genetic generalized epilepsy

Epilepsy is a complex neurological disease that can be caused by both genetic and environmental factors. Many studies have been conducted to investigate the genetic risk variants and molecular mechanisms of epilepsy. Disruption of excitation-inhibition balance (E/I balance) is one of the widely accepted disease mechanisms of epilepsy. The maintenance of E/I balance is an intricate process that is governed by multiple proteins. Using whole exome sequencing (WES), we identified a novel GABRA1 c.448G>A (p.E150K) variant and ERBB4 c.1972A>T (p.I658F, rs190654033) variant in a Malaysian Chinese family with genetic generalized epilepsy (GGE). The GGE may be triggered by dysregulation of E/I balance mechanism. Segregation of the variants in the family was verified by Sanger sequencing. All family members with GGE inherited both variants. However, family members who carried only one of the variants did not show any symptoms of GGE. Both the GABRA1 and ERBB4 variants were predicted damaging by MutationTaster and CADD, and protein structure analysis showed that the variants had resulted in the formation of additional hydrogen bonds in the mutant proteins. GABRA1 variant could reduce the efficiency of GABA_A receptors, and constitutively active ERBB4 receptors caused by the ERBB4 variant promote internalization of GABA_A receptors. The interaction between the two variants may cause a greater disruption in E/I balance, which is more likely to induce a seizure. Nevertheless, this disease model was derived from a single small family, further studies are still needed to confirm the verifiability of the purported disease model.

An extracellular receptor tyrosine kinase motif orchestrating intracellular STAT activation

The ErbB4 receptor isoforms JM-a and JM-b differ within their extracellular juxtamembrane (eJM) domains. Here, ErbB4 isoforms are used as a model to address the effect of structural variation in the eJM domain of receptor tyrosine kinases (RTK) on downstream signaling. A specific JM-a-like sequence motif is discovered, and its presence or absence (in JM-b-like RTKs) in the eJM domains of several RTKs is demonstrated to dictate selective STAT activation. STAT5a activation by RTKs including the JM-a like motif is shown to involve interaction with oligosaccharides of N-glycosylated cell surface proteins such as β1 integrin, whereas STAT5b activation by JM-b is dependent on TYK2. ErbB4 JM-a- and JM-b-like RTKs are shown to associate with specific signaling complexes at different cell surface compartments using analyses of RTK interactomes and super-resolution imaging. These findings provide evidence for a conserved mechanism linking a ubiquitous extracellular motif in RTKs with selective intracellular STAT signaling.

Lipid-Mediated Association of the Slg1 Transmembrane Domains in Yeast Plasma Membranes

Clustering of transmembrane proteins underlies a multitude of fundamental biological processes at the plasma membrane (PM) such as receptor activation, lateral domain formation, and mechanotransduction. The self-association of the respective transmembrane domains (TMDs) has also been suggested to be responsible for the micron-scaled patterns seen for integral membrane proteins in the budding yeast PM. However, the underlying interplay between the local lipid composition and the TMD identity is still not mechanistically understood. In this work, we combined coarse-grained molecular dynamics simulations of simplified bilayer systems with high-resolution live-cell microscopy to analyze the distribution of a representative helical yeast TMD from the PM sensor Slg1 within different lipid environments. In our simulations, we specifically evaluated the effects of acyl chain saturation and anionic lipid head groups on the association of two TMDs. We found that weak lipid-protein interactions significantly affect the configuration of TMD dimers and the free energy of association. Increased amounts of unsaturated phospholipids (PLs) strongly reduced the helix-helix interaction, while the presence of anionic phosphatidylserine (PS) hardly affected the dimer formation. We could experimentally confirm this surprising lack of effect of PS using the network factor, a mesoscopic measure of PM pattern formation in yeast cells. Simulations also showed that the formation of TMD dimers in turn increased the order parameter of the surrounding lipids and induced long-range perturbations in lipid organization. In summary, our results shed new light on the mechanisms of lipid-mediated dimerization of TMDs in complex lipid mixtures.

An Unbiased Functional Genetics Screen Identifies Rare Activating ERBB4 Mutations

Despite the relatively high frequency of somatic ERBB4 mutations in various cancer types, only a few activating ERBB4 mutations have been characterized, primarily due to lack of mutational hotspots in the ERBB4 gene. Here, we utilized our previously published pipeline, an in vitro screen for activating mutations, to perform an unbiased functional screen to identify potential activating ERBB4 mutations from a randomly mutated ERBB4 expression library. Ten potentially activating ERBB4 mutations were identified and subjected to validation by functional and structural analyses. Two of the 10 ERBB4 mutants, E715K and R687K, demonstrated hyperactivity in all tested cell models and promoted cellular growth under two-dimensional and three-dimensional culture conditions. ERBB4 E715K also promoted tumor growth in in vivo Ba/F3 cell mouse allografts. Importantly, all tested ERBB4 mutants were sensitive to the pan-ERBB tyrosine kinase inhibitors afatinib, neratinib, and dacomitinib. Our data indicate that rare ERBB4 mutations are potential candidates for ERBB4-targeted therapy with pan-ERBB inhibitors.

STATEMENT OF SIGNIFICANCE

ERBB4 is a member of the ERBB family of oncogenes that is frequently mutated in different cancer types but the functional impact of its somatic mutations remains unknown. Here, we have analyzed the function of over 8,000 randomly mutated ERBB4 variants in an unbiased functional genetics screen. The data indicate the presence of rare activating ERBB4 mutations in cancer, with potential to be targeted with clinically approved pan-ERBB inhibitors.

Protein yoga: Conformational versatility of the Hemolysin II C-terminal domain detailed by NMR structures for multiple states

The C-terminal domain of Bacillus cereus hemolysin II (HlyIIC), stabilizes the trans-membrane-pore formed by the HlyII toxin and may aid in target cell recognition. Initial efforts to determine the NMR structure of HlyIIC were hampered by cis/trans isomerization about the single proline at position 405 that leads to doubling of NMR resonances. We used the mutant P405M-HlyIIC that eliminates the cis proline to determine the NMR structure of the domain, which revealed a novel fold. Here, we extend earlier studies to the NMR structure determination of the cis and trans states of WT-HlyIIC that exist simultaneously in solution. The primary structural differences between the cis and trans states are in the loop that contains P405, and structurally adjacent loops. Thermodynamic linkage analysis shows that at 25 C the cis proline, which already has a large fraction of 20% in the unfolded protein, increases to 50% in the folded state due to coupling with the global stability of the domain. The P405M or P405A substitutions eliminate heterogeneity due to proline isomerization but lead to the formation of a new dimeric species. The NMR structure of the dimer shows that it is formed through domain-swapping of strand β5, the last segment of secondary structure following P405. The presence of P405 in WT-HlyIIC strongly disfavors the dimer compared to the P405M-HlyIIC or P405A-HlyIIC mutants. The WT proline may thus act as a "gatekeeper," warding off aggregative misfolding.

Transmembrane Peptides as a New Strategy to Inhibit Neuraminidase-1 Activation

Sialidases, or neuraminidases, are involved in several human disorders such as neurodegenerative, infectious and cardiovascular diseases, and cancers. Accumulative data have shown that inhibition of neuraminidases, such as NEU1 sialidase, may be a promising pharmacological target, and selective inhibitors of NEU1 are therefore needed to better understand the biological functions of this sialidase. In the present study, we designed interfering peptides (IntPep) that target a transmembrane dimerization interface previously identified in human NEU1 that controls its membrane dimerization and sialidase activity. Two complementary strategies were used to deliver the IntPep into cells, either flanked to a TAT sequence or non-tagged for solubilization in detergent micelles. Combined with molecular dynamics simulations and heteronuclear nuclear magnetic resonance (NMR) studies in membrane-mimicking environments, our results show that these IntPep are able to interact with the dimerization interface of human NEU1, to disrupt membrane NEU1 dimerization and to strongly decrease its sialidase activity at the plasma membrane. In conclusion, we report here new selective inhibitors of human NEU1 of strong interest to elucidate the biological functions of this sialidase.

Structural basis of the transmembrane domain dimerization and rotation in the activation mechanism of the TRKA receptor by nerve growth factor

Tropomyosin-receptor kinases (TRKs) are essential for the development of the nervous system. The molecular mechanism of TRKA activation by its ligand nerve growth factor (NGF) is still unsolved. Recent results indicate that at endogenous levels most of TRKA is in a monomer-dimer equilibrium and that the binding of NGF induces an increase of the dimeric and oligomeric forms of this receptor. An unsolved issue is the role of the TRKA transmembrane domain (TMD) in the dimerization of TRKA and the structural details of the TMD in the active dimer receptor. Here, we found that the TRKA-TMD can form dimers, identified the structural determinants of the dimer interface in the active receptor, and validated this interface through site-directed mutagenesis together with functional and cell differentiation studies. Using in vivo cross-linking, we found that the extracellular juxtamembrane region is reordered after ligand binding. Replacement of some residues in the juxtamembrane region with cysteine resulted in ligand-independent active dimers and revealed the preferred dimer interface. Moreover, insertion of leucine residues into the TMD helix induced a ligand-independent TRKA activation, suggesting that a rotation of the TMD dimers underlies NGF-induced TRKA activation. Altogether, our findings indicate that the transmembrane and juxtamembrane regions of TRKA play key roles in its dimerization and activation by NGF.

Thyroid Peroxidase Activity is Inhibited by Phenolic Compounds-Impact of Interaction

The aim of this study was to estimate the mode of thyroid peroxidase (TPO) inhibition by polyphenols: Chlorogenic acid, rosmarinic acid, quercetin, and rutin. All the tested polyphenols inhibited TPO; the IC50 values ranged from 0.004 mM to 1.44 mM (for rosmarinic acid and rutin, respectively). All these pure phytochemical substances exhibited different modes of TPO inhibition. Rutin and rosmarinic acid showed competitive, quercetin-uncompetitive and chlorogenic acid-noncompetitive inhibition effect on TPO. Homology modeling was used to gain insight into the 3D structure of TPO and molecular docking was applied to study the interactions of the inhibitors with their target at the molecular level. Moreover, the type and strength of mutual interactions between the inhibitors (expressed as the combination index, CI) were analyzed. Slight synergism, antagonism, and moderate antagonism were found in the case of the combined addition of the pure polyphenols. Rutin and quercetin as well as rutin and rosmarinic acid acted additively (CI = 0.096 and 1.06, respectively), while rutin and chlorogenic acid demonstrated slight synergism (CI = 0.88) and rosmarinic acid with quercetin and rosmarinic acid with chlorogenic acid showed moderate antagonism (CI = 1.45 and 1.25, respectively). The mixture of chlorogenic acid and quercetin demonstrated antagonism (CI = 1.79). All the polyphenols showed in vitro antiradical ability against 2,2'-azinobis-(3-ethylbenzothiazoline-6-sulfonic acid), ABTS. The highest ability (expressed as IC50) was exhibited by rosmarinic acid (0.12 mM) and the lowest value was ascribed to quercetin (0.45 mM).

Resolving the conformational dynamics of ErbB growth factor receptor dimers

The combinatorial dimerization of the ErbB growth factor receptors (ErbB¹- ErbB⁴) are critical for their function. Here, we have characterized the conformational dynamics of ErbB transmembrane homo-dimers and hetero-dimers by using a coarse-grain simulation framework. All dimers, except ErbB^4-4 and ErbB^1-4, exhibit at least two conformations. The reported NMR structures correspond to one of these conformations, representing the N-terminal active state in ErbB^1-1 (RH2), ErbB^2-2 (RH1) and ErbB^4-4 (RH) homo-dimers and the LH dimer in ErbB^3-3 homo-dimer, validating the computational approach. Further, we predict a right-handed ErbB^3-3 dimer conformer that warrants experimental testing. The five hetero-dimers that have not yet been experimentally resolved display prominent right-handed dimers associating by the SmXXXSm motif. Our results provide insights into the constitutive signaling of ErbB⁴ after cleavage of the extracellular region. The presence of the inactive-like dimer conformers leading to symmetric kinase domains gives clues on the autoinhibition of the receptor dimers. The dimer states characterized here represent an important step towards understanding the combinatorial cross associations in the ErbB family.

Probing the effect of membrane contents on transmembrane protein-protein interaction using solution NMR and computer simulations

The interaction between the secondary structure elements is the key process, determining the spatial structure and activity of a membrane protein. Transmembrane (TM) helix-helix interaction is known to be especially important for the function of so-called type I or bitopic membrane proteins. In the present work, we present the approach to study the helix-helix interaction in the TM domains of membrane proteins in various lipid environment using solution NMR spectroscopy and phospholipid bicelles. The technique is based on the ability of bicelles to form particles with the size, depending on the lipid/detergent ratio. To implement the approach, we report the experimental parameters of "ideal bicelle" models for four kinds of zwitterionic phospholipids, which can be also used in other structural studies. We show that size of bicelles and type of the rim-forming detergent do not affect substantially the spatial structure and stability of the model TM dimer. On the other hand, the effect of bilayer thickness on the free energy of the dimer is dramatic, while the structure of the protein is unchanged in various lipids with fatty chains having a length from 12 to 18 carbon atoms. The obtained data is analyzed using the computer simulations to find the physical origin of the observed effects.

Bacterial expression, purification and biophysical characterization of the smallest plant reticulon isoform, RTNLB13

Reticulons are a large family of integral membrane proteins that are ubiquitous in eukaryotes and play a key role in functional remodelling of the endoplasmic reticulum membrane. The reticulon family is especially large in plants, with the Arabidopsis thaliana genome containing twenty-one isoforms. Reticulons vary in length but all contain a conserved C-terminal reticulon homology domain (RHD) that associates with membranes. An understanding of the structure and membrane interactions of RHDs is key to unlocking their mechanism of function, however no three-dimensional structure has been solved. We believe that this is, in part, due to difficulties in obtaining reticulon proteins in yields sufficient for structural study. To address this, we report here the first bacterial overexpression, purification, and biophysical investigation of a reticulon protein from plants, the RTNLB13 protein from A. thaliana. RTNLB13 is the smallest plant reticulon and is made up of a single RHD. We used circular dichroism, SDS-PAGE and analytical ultracentrifugation to reveal that RTNLB13 is 45% α-helical in a number of detergent environments, monomeric at low concentrations, and capable of self-association at higher concentrations. We used solution-state NMR to screen the effect of detergent type on the fold of isotopically-enriched RTNLB13, and found that ∼60% of the expected protein peaks were broadened due to slow dynamics. This broadening points toward a large network of protein-membrane interactions throughout the sequence. We have interpreted our results in light of current literature and suggest a preliminary description of RTNLB13 structure and topology.

NMR relaxation parameters of methyl groups as a tool to map the interfaces of helix-helix interactions in membrane proteins

In the case of soluble proteins, chemical shift mapping is used to identify the intermolecular interfaces when the NOE-based calculations of spatial structure of the molecular assembly are impossible or impracticable. However, the reliability of the membrane protein interface mapping based on chemical shifts or other relevant parameters was never assessed. In the present work, we investigate the predictive power of various NMR parameters that can be used for mapping of helix-helix interfaces in dimeric TM domains. These parameters are studied on a dataset containing three structures of helical dimers obtained for two different proteins in various membrane mimetics. We conclude that the amide chemical shifts have very little predictive value, while the methyl chemical shifts could be used to predict interfaces, though with great care. We suggest an approach based on conversion of the carbon NMR relaxation parameters of methyl groups into parameters of motion, and one of such values, the characteristic time of methyl rotation, appears to be a reliable sensor of interhelix contacts in transmembrane domains. The carbon NMR relaxation parameters of methyl groups can be measured accurately and with high sensitivity and resolution, making the proposed parameter a useful tool for investigation of protein-protein interfaces even in large membrane proteins. An approach to build the models of transmembrane dimers based on perturbations of methyl parameters and TMDOCK software is suggested.

Combination of Cα-H Hydrogen Bonds and van der Waals Packing Modulates the Stability of GxxxG-Mediated Dimers in Membranes

The GxxxG motif is frequently found at the dimerization interface of a transmembrane structural motif called GAS_right, which is characterized by a short interhelical distance and a right-handed crossing angle between the helices. In GAS_right dimers, such as glycophorin A (GpA), BNIP3, and members of the ErbB family, the backbones of the helices are in contact, and they invariably display networks of 4 to 8 weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O═C hydrogen bonds). These networks of weak hydrogen bonds at the helix-helix interface are presumably stabilizing, but their energetic contribution to dimerization has yet to be determined experimentally. Here, we present a computational and experimental structure-based analysis of GAS_right dimers of different predicted stabilities, which show that a combination of van der Waals packing and Cα-H hydrogen bonding predicts the experimental trend of dimerization propensities. This finding provides experimental support for the hypothesis that the networks of Cα-H hydrogen bonds are major contributors to the free energy of association of GxxxG-mediated dimers. The structural comparison between groups of GAS_right dimers of different stabilities reveals distinct sequence as well as conformational preferences. Stability correlates with shorter interhelical distances, narrower crossing angles, better packing, and the formation of larger networks of Cα-H hydrogen bonds. The identification of these structural rules provides insight on how nature could modulate stability in GAS_right and finely tune dimerization to support biological function.

Piecing it together: Unraveling the elusive structure-function relationship in single-pass membrane receptors

The challenge of crystallizing single-pass plasma membrane receptors has remained an obstacle to understanding the structural mechanisms that connect extracellular ligand binding to cytosolic activation. For example, the complex interplay between receptor oligomerization and conformational dynamics has been, historically, only inferred from static structures of isolated receptor domains. A fundamental challenge in the field of membrane receptor biology, then, has been to integrate experimentally observable dynamics of full-length receptors (e.g. diffusion and conformational flexibility) into static structural models of the disparate domains. In certain receptor families, e.g. the ErbB receptors, structures have led somewhat linearly to a putative model of activation. In other families, e.g. the tumor necrosis factor (TNF) receptors, structures have produced divergent hypothetical mechanisms of activation and transduction. Here, we discuss in detail these and other related receptors, with the goal of illuminating the current challenges and opportunities in building comprehensive models of single-pass receptor activation. The deepening understanding of these receptors has recently been accelerated by new experimental and computational tools that offer orthogonal perspectives on both structure and dynamics. As such, this review aims to contextualize those technological developments as we highlight the elegant and complex conformational communication between receptor domains. This article is part of a Special Issue entitled: Interactions between membrane receptors in cellular membranes edited by Kalina Hristova.

New Insights into Molecular Organization of Human Neuraminidase-1: Transmembrane Topology and Dimerization Ability

Neuraminidase 1 (NEU1) is a lysosomal sialidase catalyzing the removal of terminal sialic acids from sialyloconjugates. A plasma membrane-bound NEU1 modulating a plethora of receptors by desialylation, has been consistently documented from the last ten years. Despite a growing interest of the scientific community to NEU1, its membrane organization is not understood and current structural and biochemical data cannot account for such membrane localization. By combining molecular biology and biochemical analyses with structural biophysics and computational approaches, we identified here two regions in human NEU1 - segments 139-159 (TM1) and 316-333 (TM2) - as potential transmembrane (TM) domains. In membrane mimicking environments, the corresponding peptides form stable α-helices and TM2 is suited for self-association. This was confirmed with full-size NEU1 by co-immunoprecipitations from membrane preparations and split-ubiquitin yeast two hybrids. The TM2 region was shown to be critical for dimerization since introduction of point mutations within TM2 leads to disruption of NEU1 dimerization and decrease of sialidase activity in membrane. In conclusion, these results bring new insights in the molecular organization of membrane-bound NEU1 and demonstrate, for the first time, the presence of two potential TM domains that may anchor NEU1 in the membrane, control its dimerization and sialidase activity.

Helix-helix interactions in membrane domains of bitopic proteins: Specificity and role of lipid environment

Interaction between transmembrane helices often determines biological activity of membrane proteins. Bitopic proteins, a broad subclass of membrane proteins, form dimers containing two membrane-spanning helices. Some aspects of their structure-function relationship cannot be fully understood without considering the protein-lipid interaction, which can determine the protein conformational ensemble. Experimental and computer modeling data concerning transmembrane parts of bitopic proteins are reviewed in the present paper. They highlight the importance of lipid-protein interactions and resolve certain paradoxes in the behavior of such proteins. Besides, some properties of membrane organization provided a clue to understanding of allosteric interactions between distant parts of proteins. Interactions of these kinds appear to underlie a signaling mechanism, which could be widely employed in the functioning of many membrane proteins. Treatment of membrane proteins as parts of integrated fine-tuned proteolipid system promises new insights into biological function mechanisms and approaches to drug design. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Impact of Differential Detergent Interactions on Transmembrane Helix Dimerization Affinities

Interactions between transmembrane (TM) helices play a critical role in the fundamental processes required for cells to communicate and exchange materials with their surroundings. Our understanding of the factors that promote TM helix interactions has greatly benefited from our ability to study these interactions in the solution phase through the use of membrane-mimetic micelles. However, less is known about the potential influence of juxtamembrane regions flanking the interacting TM helices that may modulate dimerization affinities, even when the interacting surface itself is not altered. To investigate this question, we used solution NMR to quantitate the dimerization affinity of the major coat protein from the M13 bacteriophage in sodium dodecyl sulfate (SDS), a well-characterized model of a single-spanning self-associating TM protein. Here, we showed that a shorter construct lacking the N-terminal amphipathic helix has a higher dimerization affinity relative to that of the full-length protein, with no change in the helical structure between the monomeric and dimeric states in both cases. Although this translated into a 0.6 kcal/mol difference in free energy when the SDS solvent was approximated as a continuous phase, there were deviations from this model at high protein to detergent ratios. Instead, the equilibria were better fit to a model that treats the empty micelle as an active participant in the reaction, giving rise to standard free energies of association that were the same for both full-length and TM-segment constructs. According to this model, the higher apparent affinity of the shorter peptide could be completely explained by the enhanced detergent binding by the monomer relative to that bound by the dimer. Therefore, differential detergent binding between the monomeric and dimeric states provides a mechanism by which TM helix interactions can be modulated by noninteracting juxtamembrane regions.

Understanding single-pass transmembrane receptor signaling from a structural viewpoint-what are we missing?

Single-pass transmembrane receptors are involved in essential processes of both physiological and pathological nature and represent more than 1300 proteins in the human genome. Despite the high biological relevance of these receptors, the mechanisms of the signal transductions they facilitate are incompletely understood. One major obstacle is the lack of structures of the transmembrane domains that connect the extracellular ligand-binding domains to the intracellular signaling platforms. Over a period of almost 20 years since the first structure was reported, only 21 of these receptors have become represented by a transmembrane domain structure. This scarceness stands in strong contrast to the significance of these transmembrane α-helices for receptor functionality. In this review, we explore the properties and qualities of the current set of structures, as well as the methodological difficulties associated with their characterization and the challenges left to be overcome. Without an increased and focused effort to bring this class of proteins on par with the remaining membrane protein field, a serious lag in their biological understanding looms. Design of pharmaceutical agents, prediction of mutational affects in relation to disease, and deciphering of functional mechanisms require high-resolution structural information, especially when dealing with a domain carrying so much functionality in so few residues.

Characterization of Small Isotropic Bicelles with Various Compositions

Structural studies of membrane proteins are of great importance and interest, with solution and solid state NMR spectroscopy being very promising tools for that task. However, such investigations are hindered by a number of obstacles, and in the first place by the fact that membrane proteins need an adequate environment that models the cell membrane. One of the most widely used and prospective membrane mimetics is isotropic bicelles. While large anisotropic bicelles are well-studied, the field of small bicelles contains a lot of "white spots". The present work reports the radii of particles and concentration of the detergents in the monomeric state in solutions of isotropic bicelles, formed by 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC), 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS), 3-[(3-cholamidopropyl)dimethylammonio]-2-hydroxy-1-propanesulfonate (CHAPSO), and sodium cholate, as a function of lipid/detergent ratio and temperature. These parameters were measured using (1)H NMR diffusion spectroscopy for the bicelles composed of lipids with saturated fatty chains of different length and lipids, containing unsaturated fatty acid residue. The influence of a model transmembrane protein (membrane domain of rat TrkA) on the properties of bicelles and the effect of the bicelle size and composition on the properties of the transmembrane protein were investigated with heteronuclear NMR and nuclear Overhauser effect spectroscopy. We show that isotropic bicelles that are applicable for solution NMR spectroscopy behave as predicted by the theoretical models and are likely to be bicelles rather than mixed micelles. Using the obtained data, we propose a simple approach to control the size of bicelles at low concentrations. On the basis of our results, we compared different rim-forming agents and selected CHAPS as a detergent of choice for structural studies in bicelles, if the deuteration of the detergent is not required.

Structural Basis of p75 Transmembrane Domain Dimerization

Dimerization of single span transmembrane receptors underlies their mechanism of activation. p75 neurotrophin receptor plays an important role in the nervous system, but the understanding of p75 activation mechanism is still incomplete. The transmembrane (TM) domain of p75 stabilizes the receptor dimers through a disulfide bond, essential for the NGF signaling. Here we solved by NMR the three-dimensional structure of the p75-TM-WT and the functionally inactive p75-TM-C257A dimers. Upon reconstitution in lipid micelles, p75-TM-WT forms the disulfide-linked dimers spontaneously. Under reducing conditions, p75-TM-WT is in a monomer-dimer equilibrium with the Cys(257) residue located on the dimer interface. In contrast, p75-TM-C257A forms dimers through the AXXXG motif on the opposite face of the α-helix. Biochemical and cross-linking experiments indicate that AXXXG motif is not on the dimer interface of p75-TM-WT, suggesting that the conformation of p75-TM-C257A may be not functionally relevant. However, rather than mediating p75 homodimerization, mutagenesis of the AXXXG motif reveals its functional role in the regulated intramembrane proteolysis of p75 catalyzed by the γ-secretase complex. Our structural data provide an insight into the key role of the Cys(257) in stabilization of the weak transmembrane dimer in a conformation required for the NGF signaling.

Alternative packing of EGFR transmembrane domain suggests that protein-lipid interactions underlie signal conduction across membrane

The human epidermal growth factor receptor (EGFR) of HER/ErbB receptor tyrosine kinase family mediates a broad spectrum of cellular responses transducing biochemical signals via lateral dimerization in plasma membrane, while inactive receptors can exist in both monomeric and dimeric forms. Recently, the dimeric conformation of the helical single-span transmembrane domains of HER/ErbB employing the relatively polar N-terminal motifs in a fashion permitting proper kinase activation was experimentally determined. Here we describe the EGFR transmembrane domain dimerization via an alternative weakly polar C-terminal motif A(661)xxxG(665) presumably corresponding to the inactive receptor state. During association, the EGFR transmembrane helices undergo a structural adjustment with adaptation of inter-molecular polar and hydrophobic interactions depending upon the surrounding membrane properties that directly affect the transmembrane helix packing. This might imply that signal transduction through membrane and allosteric regulation are inclusively mediated by coupled protein-protein and protein-lipid interactions, elucidating paradoxically loose linkage between ligand binding and kinase activation.

Role of GxxxG Motifs in Transmembrane Domain Interactions

Transmembrane (TM) helices of integral membrane proteins can facilitate strong and specific noncovalent protein-protein interactions. Mutagenesis and structural analyses have revealed numerous examples in which the interaction between TM helices of single-pass membrane proteins is dependent on a GxxxG or (small)xxx(small) motif. It is therefore tempting to use the presence of these simple motifs as an indicator of TM helix interactions. In this Current Topic review, we point out that these motifs are quite common, with more than 50% of single-pass TM domains containing a (small)xxx(small) motif. However, the actual interaction strength of motif-containing helices depends strongly on sequence context and membrane properties. In addition, recent studies have revealed several GxxxG-containing TM domains that interact via alternative interfaces involving hydrophobic, polar, aromatic, or even ionizable residues that do not form recognizable motifs. In multipass membrane proteins, GxxxG motifs can be important for protein folding, and not just oligomerization. Our current knowledge thus suggests that the presence of a GxxxG motif alone is a weak predictor of protein dimerization in the membrane.

Evolutionary-guided de novo structure prediction of self-associated transmembrane helical proteins with near-atomic accuracy

How specific protein associations regulate the function of membrane receptors remains poorly understood. Conformational flexibility currently hinders the structure determination of several classes of membrane receptors and associated oligomers. Here we develop EFDOCK-TM, a general method to predict self-associated transmembrane protein helical (TMH) structures from sequence guided by co-evolutionary information. We show that accurate intermolecular contacts can be identified using a combination of protein sequence covariation and TMH binding surfaces predicted from sequence. When applied to diverse TMH oligomers, including receptors characterized in multiple conformational and functional states, the method reaches unprecedented near-atomic accuracy for most targets. Blind predictions of structurally uncharacterized receptor tyrosine kinase TMH oligomers provide a plausible hypothesis on the molecular mechanisms of disease-associated point mutations and binding surfaces for the rational design of selective inhibitors. The method sets the stage for uncovering novel determinants of molecular recognition and signalling in single-spanning eukaryotic membrane receptors.

Intramembrane binding of VE-cadherin to VEGFR2 and VEGFR3 assembles the endothelial mechanosensory complex

VE-cadherin plays a critical role in endothelial shear stress mechanotransduction by interacting with VEGFRs through their transmembrane domains.

Endothelial responses to fluid shear stress are essential for vascular development and physiology, and determine the formation of atherosclerotic plaques at regions of disturbed flow. Previous work identified VE-cadherin as an essential component, along with PECAM-1 and VEGFR2, of a complex that mediates flow signaling. However, VE-cadherin’s precise role is poorly understood. We now show that the transmembrane domain of VE-cadherin mediates an essential adapter function by binding directly to the transmembrane domain of VEGFR2, as well as VEGFR3, which we now identify as another component of the junctional mechanosensory complex. VEGFR2 and VEGFR3 signal redundantly downstream of VE-cadherin. Furthermore, VEGFR3 expression is observed in the aortic endothelium, where it contributes to flow responses in vivo. In summary, this study identifies a novel adapter function for VE-cadherin mediated by transmembrane domain association with VEGFRs.

Structural and functional characterization of alternative transmembrane domain conformations in VEGF receptor 2 activation

Transmembrane signaling by receptor tyrosine kinases (RTKs) entails ligand-mediated dimerization and structural rearrangement of the extracellular domains. RTK activation also depends on the specific orientation of the transmembrane domain (TMD) helices, as suggested by pathogenic, constitutively active RTK mutants. Such mutant TMDs carry polar amino acids promoting stable transmembrane helix dimerization, which is essential for kinase activation. We investigated the effect of polar amino acids introduced into the TMD of vascular endothelial growth factor receptor 2, regulating blood vessel homeostasis. Two mutants showed constitutive kinase activity, suggesting that precise TMD orientation is mandatory for kinase activation. Nuclear magnetic resonance spectroscopy revealed that TMD helices in activated constructs were rotated by 180° relative to the interface of the wild-type conformation, confirming that ligand-mediated receptor activation indeed results from transmembrane helix rearrangement. A molecular dynamics simulation confirmed the transmembrane helix arrangement of wild-type and mutant TMDs revealed by nuclear magnetic resonance spectroscopy.

Tie2 and Eph receptor tyrosine kinase activation and signaling

The Eph and Tie cell surface receptors mediate a variety of signaling events during development and in the adult organism. As other receptor tyrosine kinases, they are activated on binding of extracellular ligands and their catalytic activity is tightly regulated on multiple levels. The Eph and Tie receptors display some unique characteristics, including the requirement of ligand-induced receptor clustering for efficient signaling. Interestingly, both Ephs and Ties can mediate different, even opposite, biological effects depending on the specific ligand eliciting the response and on the cellular context. Here we discuss the structural features of these receptors, their interactions with various ligands, as well as functional implications for downstream signaling initiation. The Eph/ephrin structures are already well reviewed and we only provide a brief overview on the initial binding events. We go into more detail discussing the Tie-angiopoietin structures and recognition.

A frequent, GxxxG-mediated, transmembrane association motif is optimized for the formation of interhelical Cα-H hydrogen bonds

Carbon hydrogen bonds between Cα-H donors and carbonyl acceptors are frequently observed between transmembrane helices (Cα-H···O=C). Networks of these interactions occur often at helix-helix interfaces mediated by GxxxG and similar patterns. Cα-H hydrogen bonds have been hypothesized to be important in membrane protein folding and association, but evidence that they are major determinants of helix association is still lacking. Here we present a comprehensive geometric analysis of homodimeric helices that demonstrates the existence of a single region in conformational space with high propensity for Cα-H···O=C hydrogen bond formation. This region corresponds to the most frequent motif for parallel dimers, GASright, whose best-known example is glycophorin A. The finding suggests a causal link between the high frequency of occurrence of GASright and its propensity for carbon hydrogen bond formation. Investigation of the sequence dependency of the motif determined that Gly residues are required at specific positions where only Gly can act as a donor with its "side chain" Hα. Gly also reduces the steric barrier for non-Gly amino acids at other positions to act as Cα donors, promoting the formation of cooperative hydrogen bonding networks. These findings offer a structural rationale for the occurrence of GxxxG patterns at the GASright interface. The analysis identified the conformational space and the sequence requirement of Cα-H···O=C mediated motifs; we took advantage of these results to develop a structural prediction method. The resulting program, CATM, predicts ab initio the known high-resolution structures of homodimeric GASright motifs at near-atomic level.

Primary and secondary dimer interfaces of the fibroblast growth factor receptor 3 transmembrane domain: characterization via multiscale molecular dynamics simulations

Receptor tyrosine kinases are single-pass membrane proteins that form dimers within the membrane. The interactions of their transmembrane domains (TMDs) play a key role in dimerization and signaling. Fibroblast growth factor receptor 3 (FGFR3) is of interest as a G380R mutation in its TMD is the underlying cause of ~99% of the cases of achondroplasia, the most common form of human dwarfism. The structural consequences of this mutation remain uncertain: the mutation shifts the position of the TMD relative to the lipid bilayer but does not alter the association free energy. We have combined coarse-grained and all-atom molecular dynamics simulations to study the dimerization of wild-type, heterodimer, and mutant FGFR3 TMDs. The simulations reveal that the helices pack together in the dimer to form a flexible interface. The primary packing mode is mediated by a Gx3G motif. There is also a secondary dimer interface that is more highly populated in heterodimer and mutant configurations that may feature in the molecular mechanism of pathology. Both coarse-grained and atomistic simulations reveal a significant shift of the G380R mutant dimer TMD relative to the bilayer to allow interactions of the arginine side chain with lipid headgroup phosphates.

Structure of FGFR3 transmembrane domain dimer: implications for signaling and human pathologies

Fibroblast growth factor receptor 3 (FGFR3) transduces biochemical signals via lateral dimerization in the plasma membrane, and plays an important role in human development and disease. Eight different pathogenic mutations, implicated in cancers and growth disorders, have been identified in the FGFR3 transmembrane segment. Here, we describe the dimerization of the FGFR3 transmembrane domain in membrane-mimicking DPC/SDS (9/1) micelles. In the solved NMR structure, the two transmembrane helices pack into a symmetric left-handed dimer, with intermolecular stacking interactions occurring in the dimer central region. Some pathogenic mutations fall within the helix-helix interface, whereas others are located within a putative alternative interface. This implies that although the observed dimer structure is important for FGFR3 signaling, the mechanism of FGFR3-mediated transduction across the membrane is complex. We propose an FGFR3 signaling mechanism that is based on the solved structure, available structures of isolated soluble FGFR domains, and published biochemical and biophysical data.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.