Elucidation of the molecular interactions that enable stable assembly and structural diversity in multicomponent immune receptors

Multicomponent immune receptors are essential complexes in which distinct ligand-recognition and signaling subunits are held together by interactions between acidic and basic residues of their transmembrane helices. A 2:1 acidic-to-basic motif in the transmembrane domains of the subunits is necessary and sufficient to assemble these receptor complexes. Here, we study a prototype for these receptors, a DAP12-NKG2C 2:1 heterotrimeric complex, in which the two DAP12 subunits each contribute a single transmembrane Asp residue, and the NKG2C subunit contributes a Lys to form the complex. DAP12 can also associate with 20 other subunits using a similar motif. Here, we use molecular-dynamics simulations to understand the basis for the high affinity and diversity of interactions in this group of receptors. Simulations of the transmembrane helices with differing protonation states of the Asp-Asp-Lys triad identified a structurally stable interaction in which a singly-protonated Asp-Asp pair forms a hydrogen-bonded carboxyl-carboxylate clamp that clasps onto a charged Lys side chain. This polar motif was also supported by density functional theory and a Protein Data Bank-wide search. In contrast, the helices are dynamic at sites distal to the stable carboxyl-carboxylate clamp motif. Such a locally stable but globally dynamic structure is well suited to accommodate the sequence and structural variations in the transmembrane helices of multicomponent receptors, which mix and match subunits to create combinatorial functional diversity from a limited number of subunits. It also supports a signaling mechanism based on multisubunit clustering rather than propagation of rigid conformational changes through the membrane.

Coarse Grained Molecular Dynamic Simulations for the Study of TNF Receptor Family Members' Transmembrane Organization

The Tumor Necrosis Factor (TNF) and the TNF receptor (TNFR) superfamilies are composed of 19 ligands and 30 receptors, respectively. The oligomeric properties of ligands, both membrane bound and soluble, has been studied most. However, less is known about the oligomeric properties of TNFRs. Earlier reports identified the extracellular, membrane-distal, cysteine-rich domain as a pre-ligand assembly domain which stabilizes receptor dimers and/or trimers in the absence of ligand. Nevertheless, recent reports based on structural nuclear magnetic resonance (NMR) highlight the intrinsic role of the transmembrane domains to form dimers (p75NTR), trimers (Fas), or dimers of trimers (DR5). Thus, understanding the structural basis of transmembrane oligomerization may shed light on the mechanism for signal transduction and the impact of disease-associated mutations in this region. To this end, here we used an in silico coarse grained molecular dynamics approach with Martini force field to study TNFR transmembrane homotypic interactions. We have first validated this approach studying the three TNFR described by NMR (p75NTR, Fas, and DR5). We have simulated membrane patches composed of 36 helices of the same receptor equidistantly distributed in order to get unbiassed information on spontaneous proteins assemblies. Good agreement was found in the specific residues involved in homotypic interactions and we were able to observe dimers, trimers, and higher-order oligomers corresponding to those reported in NMR experiments. We have, applied this approach to study the assembly of disease-related mutations being able to assess their impact on oligomerization stability. In conclusion, our results showed the usefulness of coarse grained simulations with Martini force field to study in an unbiased manner higher order transmembrane oligomerization.

Transmembrane Helix Assembly by Max-Min Ant System Algorithm

Because of the rapid progress in biochemical and structural studies of membrane proteins, considerable attention has been given on developing efficient computational methods for solving low-to-medium resolution structures using sparse structural data. In this study, we demonstrate a novel algorithm, max-min ant system (MMAS), designed to find an assembly of α-helical transmembrane proteins using a rigid helix arrangement guided by distance constraints. The new algorithm generates a large variety with finite number of orientations of transmembrane helix bundle and finds the solution that is matched with the provided distance constraints based on the behavior of ants to search for the shortest possible path between their nest and the food source. To demonstrate the efficiency of the novel search algorithm, MMAS is applied to determine the transmembrane packing of KcsA and MscL ion channels from a limited distance information extracted from the crystal structures, and the packing of KvAP voltage sensor domain using a set of 10 experimentally determined constraints, and the results are compared with those of two popular used stochastic methods, simulated annealing Monte Carlo method and genetic algorithm.