The Effect of a Mutation in the Thyroid Stimulating Hormone Receptor (TSHR) on Development,Behaviour and TH Levels in Domesticated Chickens

The thyroid stimulating hormone receptor (TSHR) has been suggested to be a “domestication locus” in the chicken, due to a strong selective sweep over the gene found in domesticated chickens, differentiating them from their wild ancestor the Red Junglefowl (RJF). We investigated the effect of the mutation on development (incubation time), behaviour and thyroid hormone levels in intercross chickens homozygous for the mutation (d/d), wild type homozygotes (w/w) or heterozygotes (d/w). This allowed an assessment of the effect of genotype at this locus against a random mix of RJF and WL genotypes throughout the rest of the genome, controlling for family effects. The d/d genotype showed a longer incubation time, less fearful behaviours, lower number of aggressive behaviours and decreased levels of the thyroid hormone T4, in comparison to the w/w genotype. The difference between TSHR genotypes (d/d vs. w/w) in these respects mirrors the differences in development and behaviour between pure domesticated White Leghorns and pure RJF chickens. Higher individual T3 and T4 levels were associated with more aggression. Our study indicates that the TSHR mutation affects typical domestication traits, possibly through modifying plasma levels of thyroid hormones, and may therefore have been important during the evolution of the domestic chicken.     

鸡SNP多样性的比较研究与群体有效规模的估算

以红色原鸡(Red Jungle Fowl, RJF)、丝羽乌骨鸡(Taihe silk Chicken, TS)、隐性白洛克鸡(White Recessive Rock, WRR)为资源群, 在鸡一号染色体的Contig.060226.1上选取了一个200 kb的区域, 比较研究了3个群体的SNP (single nucleotide polymorphism)多样性、估算了鸡的初始群体有效规模大小(effective size of population, Ne)。红色原鸡、丝羽乌骨鸡、隐性白洛克鸡3个群体的平均杂合度分别为0.28533±0.03475、0.32926±0.03919、0.30168±0.04038。显著性检验差异不显著(P=0.2368>0.05)。根据Latter 和Nei的方法对鸡的群体有效规模进行了估算, 鸡的初始群体有效规模大小为20 000～150 000。鸡在驯养的早期阶段经历过严厉的瓶颈效应, 但瓶颈效应对鸡各品种SNP的多样性并未产生显著影响。笔者认为, 鸡在驯养的早期阶段群体有效规模足够大, 品种分化过程中群体迅速扩张, 品种间的广泛杂交(特别是和红色原鸡之间)以及鸡基因组的高重组率等因素是导致家鸡和原鸡以及各家鸡品种间SNP多样性没有显著差别的重要原因。     

Molecular Insights into the Transmembrane Domain of the Thyrotropin Receptor

The thyrotropin receptor (TSHR) is a G protein-coupled receptor (GPCR) that is member of the leucine-rich repeat subfamily (LGR). In the absence of crystal structure, the success of rational design of ligands targeting the receptor internal cavity depends on the quality of the TSHR models built. In this subfamily, transmembrane helices (TM) 2 and 5 are characterized by the absence of proline compared to most receptors, raising the question of the structural conformation of these helices. To gain insight into the structural properties of these helices, we carried out bioinformatics and experimental studies. Evolutionary analysis of the LGR family revealed a deletion in TM5 but provided no information on TM2. Wild type residues at positions 2.58, 2.59 or 2.60 in TM2 and/or at position 5.50 in TM5 were substituted to proline. Depending on the position of the proline substitution, different effects were observed on membrane expression, glycosylation, constitutive cAMP activity and responses to thyrotropin. Only proline substitution at position 2.59 maintained complex glycosylation and high membrane expression, supporting occurrence of a bulged TM2. The TSHR transmembrane domain was modeled by homology with the orexin 2 receptor, using a protocol that forced the deletion of one residue in the TM5 bulge of the template. The stability of the model was assessed by molecular dynamics simulations. TM5 straightened during the equilibration phase and was stable for the remainder of the simulations. Our data support a structural model of the TSHR transmembrane domain with a bulged TM2 and a straight TM5 that is specific of glycoprotein hormone receptors.     

The mechanism of proton translocation in respiratory complex I from molecular dynamics

Abstract

Respiratory complex I, the biggest enzyme of respiratory chain, plays a key role in energy production by the mitochondrial respiratory chain and has been implicated in many human neurodegenerative diseases. Recently, the crystal structure of respiratory complex I is reported. We perform 50?ns molecular dynamics simulations on the membrane domain of respiratory complex I under two hypothetical states (oxidized state and reduced state). We find that the density of water molecules in the trans-membrane domain under reduced state is bigger than that under oxidized state. The connecting elements (helix HL and β-hairpins-helix element) fluctuate stronger under reduced state than that under oxidized state, causing more internal water molecules and facilitating the proton conduction. The conformational changes of helix HL and the crucial charged residue Glu in TM5 play key roles in the mechanism of proton translocation. Our results illustrate the dynamic behavior and the potential mechanism of respiratory complex I, which provides the structural basis for drug design of respiratory complex I.     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.     

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.     

Transmembrane peptides used to investigate the homo-oligomeric interface and binding hotspot of latent membrane protein 1

Epstein-Barr virus (EBV), a human γ-herpesvirus, establishes lifelong infection by targeting the adaptive immune system of the host through memory B cells. Although normally benign, EBV contributes to lymphoid malignancies and lymphoproliferative syndromes in immunocompromised individuals. The viral oncoprotein latent membrane protein 1 (LMP-1) is essential for B lymphocyte immortalization by EBV. The constitutive signaling activity of LMP-1 is dependent on homo-oligomerization of its six-spanning hydrophobic transmembrane domain (TMD). However, the mechanism driving LMP-1 intermolecular interaction is poorly understood. Here, we show that the fifth transmembrane helix (TM5) of LMP-1 strongly self-associates, forming a homotrimeric complex mediated by a polar residue embedded in the membrane, D150. Replacement of this aspartic acid residue with alanine disrupts TM5 self-association in detergent micelles and bacterial cell membranes. A full-length LMP-1 variant harboring the D150A substitution is deficient in NFκB activation, supporting the key role of the fifth transmembrane helix in constitutive activation of signaling by this oncoprotein.     

Identification of secretin,vasoactive intestinal peptide and glucagon binding sites: from chimaeric receptors to point mutations

We have identified two basic residues that are important for the recognition of secretin and vasoactive intestinal peptide (VIP) by their respective receptors. These two peptides containing an Asp residue at position 3 interacted with an arginine residue in transmembrane helix 2 (TM2) of the receptor, and the lysine residue in extracellular loop 1 (ECL1) stabilized the active receptor conformation induced by the ligand. The glucagon receptor possesses a Lys instead of an Arg in TM2, and an Ile instead of Lys in ECL1; it markedly prefers a Gln side chain in position 3 of the ligand. Our results suggested that, in the wild-type receptor, the Ile side chain prevented access to the TM2 Lys side chain, but oriented the glucagon Gln(3) side chain to its proper binding site. In the double mutant, the ECL1 Lys allowed an interaction between negatively charged residues in position 3 of glucagon and the TM2 Arg, resulting in efficient receptor activation by [Asp(3)]glucagon as well as by glucagon.     

A Tyrosine Residue on the TSH Receptor Stabilizes Multimer Formation

Background

The thyrotropin stimulating hormone receptor (TSHR) is a G protein coupled receptor (GPCR) with a large ectodomain. The ligand, TSH, acting via this receptor regulates thyroid growth and thyroid hormone production and secretion. The TSH receptor (TSHR) undergoes complex post –translational modifications including intramolecular cleavage and receptor multimerization. Since monomeric and multimeric receptors coexist in cells, understanding the functional role of just the TSHR multimers is difficult. Therefore, to help understand the physiological significance of receptor multimerization, it will be necessary to abrogate multimer formation, which requires identifying the ectodomain and endodomain interaction sites on the TSHR. Here, we have examined the contribution of the ectodomain to constitutive multimerization of the TSHR and determined the possible residue(s) that may be involved in this interaction.

Methodology/Principal Findings

We studied ectodomain multimer formation by expressing the extracellular domain of the TSHR linked to a glycophosphotidyl (GPI) anchor in both stable and transient expression systems. Using co-immunoprecipitation and FRET of tagged receptors, we established that the TSH receptor ectodomain was capable of multimerization even when totally devoid of the transmembrane domain. Further, we studied the effect of two residues that likely made critical contact points in this interaction. We showed that a conserved tyrosine residue (Y116) on the convex surface of the LRR3 was a critical residue in ectodomain multimer formation since mutation of this residue to serine totally abrogated ectodomain multimers. This abrogation was not seen with the mutation of cysteine 176 on the inner side of the LRR5, demonstrating that inter-receptor disulfide bonding was not involved in ectodomain multimer formation. Additionally, the Y116 mutation in the intact wild type receptor enhanced receptor degradation.

Conclusions/Significance

These data establish the TSH receptor ectodomain as one site of multimerization, independent of the transmembrane region, and that this interaction was primarily via a conserved tyrosine residue in LRR3.     

Roles of Subunit NuoK (ND4L) in the Energy-transducing Mechanism of Escherichia coli NDH-1 (NADH:Quinone Oxidoreductase)

The bacterial H⁺-translocating NADH:quinone oxidoreductase (NDH-1) catalyzes electron transfer from NADH to quinone coupled with proton pumping across the cytoplasmic membrane. The NuoK subunit (counterpart of the mitochondrial ND4L subunit) is one of the seven hydrophobic subunits in the membrane domain and bears three transmembrane segments (TM1–3). Two glutamic residues located in the adjacent transmembrane helices of NuoK are important for the energy coupled activity of NDH-1. In particular, mutation of the highly conserved carboxyl residue (_KGlu-36 in TM2) to Ala led to a complete loss of the NDH-1 activities. Mutation of the second conserved carboxyl residue (_KGlu-72 in TM3) moderately reduced the activities. To clarify the contribution of NuoK to the mechanism of proton translocation, we relocated these two conserved residues. When we shifted _KGlu-36 along TM2 to positions 32, 38, 39, and 40, the mutants largely retained energy transducing NDH-1 activities. According to the recent structural information, these positions are located in the vicinity of _KGlu-36, present in the same helix phase, in an immediately before and after helix turn. In an earlier study, a double mutation of two arginine residues located in a short cytoplasmic loop between TM1 and TM2 (loop-1) showed a drastic effect on energy transducing activities. Therefore, the importance of this cytosolic loop of NuoK (_KArg-25, _KArg-26, and _KAsn-27) for the energy transducing activities was extensively studied. The probable roles of subunit NuoK in the energy transducing mechanism of NDH-1 are discussed.     

Influence of proline residues in transmembrane helix packing

Integral membrane proteins often contain proline residues in their alpha-helical transmembrane (TM) fragments, which may strongly influence their folding and association. Pro-scanning mutagenesis of the helical domain of glycophorin A (GpA) showed that replacement of the residues located at the center abrogates helix packing while substitution of the residues forming the ending helical turns allows dimer formation. Synthetic TM peptides revealed that a point mutation of one of the residues of the dimerization motif (L75P) located at the N-terminal helical turn of the GpA TM fragment, adopts a secondary structure and oligomeric state similar to the wild-type sequence in detergents. In addition, both glycosylation mapping in biological membranes and molecular dynamics showed that the presence of a proline residue at the lipid/water interface has as an effect the extension of the helical end. Thus, helix packing can be an important factor that determines appearance of proline in TM helices. Membrane proteins might accumulate proline residues at the two ends of their TM segments in order to modulate the exposition of key amino acid residues at the interface for molecular recognition events while allowing stable association and native folding.     

Insights into the transmembrane helix associations of kit ligand by molecular dynamics simulation and TOXCAT

Kit ligand (KITL) plays important roles in cell proliferation, differentiation, and survival via interaction with its receptor Kit. The previous studies demonstrated that KITL formed a noncovalent homodimer through transmembrane (TM) domain; however, the undergoing mechanism of transmembrane association that determines KITL TM dimerization is still not clear. Herein, molecular dynamics (MD) simulation strategy and TOXCAT assay were combined to characterize the dimerization interface and structure of KITL TM in details. KITL TM formed a more energetically favorable noncovalent dimer through a conserved SxxxGxxxG motif in the MD simulation. Furthermore, the TOXCAT results demonstrated that KITL TM self‐associated strongly in the bilayer membrane environment. Mutating any one of the small residues Ser11, Gly15 or Gly19 to Ile disrupted KITL TM dimerization dramatically, which further validated our MD simulation results. In addition, our results showed that Tyr22 could help to stabilize the TM interactions via interacting with the phosphoric group in the bilayer membrane. Pro7 did not induce helix kinks or swivel angles in KITL TM, but it was related with the pitch of the turn around this residue so as to affect the dimer formation. Combining the results of computer modeling and experimental mutagenesis studies on the KITL TM provide new insights for the transmembrane helix association of KITL dimerization. Proteins 2017; 85:1362–1370. © 2017 Wiley Periodicals, Inc.     

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.     

Incorporation of transmembrane peptides from the vacuolar H(+)-ATPase in phospholipid membranes: spin-label electron paramagnetic resonance and polarized infrared spectroscopy

Peptides were designed that are based on candidate transmembrane sequences of the V o-sector from the vacuolar H (+)-ATPase of Saccharomyces cerevisiae. Spin-label EPR studies of lipid-protein interactions were used to characterize the state of oligomerization, and polarized IR spectroscopy was used to determine the secondary structure and orientation, of these peptides in lipid bilayer membranes. Peptides corresponding to the second and fourth transmembrane domains (TM2 and TM4) of proteolipid subunit c (Vma3p) and of the putative seventh transmembrane domain (TM7) of subunit a (Vph1p) are wholly, or predominantly, alpha-helical in membranes of dioleoyl phosphatidylcholine. All three peptides self-assemble into oligomers of different sizes, in which the helices are differently inclined with respect to the membrane normal. The coassembly of rotor (Vma3p TM4) and stator (Vph1p TM7) peptides, which respectively contain the glutamate and arginine residues essential to proton transport by the rotary ATPase mechanism, is demonstrated from changes in the lipid interaction stoichiometry and helix orientation. Concanamycin, a potent V-ATPase inhibitor, and a 5-(2-indolyl)-2,4-pentadienoyl inhibitor that exhibits selectivity for the osteoclast subtype, interact with the membrane-incorporated Vma3p TM4 peptide, as evidenced by changes in helix orientation; concanamycin additionally interacts with Vph1p TM7, suggesting that both stator and rotor elements contribute to the inhibitor site within the membrane. Comparison of the peptide behavior in lipid bilayers is made with membranous subunit c assemblies of the 16-kDa proteolipid from Nephrops norvegicus, which can substitute functionally for Vma3p in S. cerevisiae.     

Structural and functional characterization of transmembrane segment VII of the Na+/H+ exchanger isoform 1

The Na(+)/H(+) exchanger isoform 1 is an integral membrane protein that regulates intracellular pH by exchanging one intracellular H(+) for one extracellular Na(+). It is composed of an N-terminal membrane domain of 12 transmembrane segments and an intracellular C-terminal regulatory domain. We characterized the structural and functional aspects of the critical transmembrane segment VII (TM VII, residues 251-273) by using alanine scanning mutagenesis and high resolution NMR. Each residue of TM VII was mutated to alanine, the full-length protein expressed, and its activity characterized. TM VII was sensitive to mutation. Mutations at 13 of 22 residues resulted in severely reduced activity, whereas other mutants exhibited varying degrees of decreases in activity. The impaired activities sometimes resulted from low expression and/or low surface targeting. Three of the alanine scanning mutant proteins displayed increased, and two displayed decreased resistance to the Na(+)/H(+) exchanger isoform 1 inhibitor EMD87580. The structure of a peptide of TM VII was determined by using high resolution NMR in dodecylphosphocholine micelles. TM VII is predominantly alpha-helical, with a break in the helix at the functionally critical residues Gly(261)-Glu(262). The relative positions and orientations of the N- and C-terminal helical segments are seen to vary about this extended segment in the ensemble of NMR structures. Our results show that TM VII is a critical transmembrane segment structured as an interrupted helix, with several residues that are essential to both protein function and sensitivity to inhibition.     

Repulsive separation of the cytoplasmic ends of transmembrane helices 3 and 6 is linked to receptor activation in a novel thyrotropin receptor mutant (M626I)

Ligand-dependent activation of G protein-coupled receptors (GPCRs) involves repositioning of the juxtacytoplasmic ends of transmembrane helices TM3 and TM6. This concept, inferred from site-directed spin labeling studies, is supported by chemical cross-linking of the cytoplasmic ends of TM3 and TM6 blocking GPCR activation. Here we report a novel constitutive active mutation (M626I) in TM6 of the TSH receptor (TSHR), identified in affected members of a family with nonautoimmune hyperthyroidism. The specific constitutive activity of M626I, measured by its basal cAMP generation corrected for cell surface expression, was 13-fold higher than that of wild-type TSHR. Homology modeling of the TSHR serpentine domain based on the rhodopsin crystal structure suggests that M626 faces the side chain of I515 of TM3 near the membrane-cytoplasmic junction. Steric hindrance of the introduced isoleucine by I515 is consistent with the fact that shorter or more flexible side chains at position 626 did not increase constitutivity. Furthermore, a reciprocal mutation at position 515 (I515M), when introduced into the M626I background, acts as revertant mutation by allowing accommodation of the isoleucine sidechain at position 626 and fully restoring the constitutive activity to the level of wild-type TSHR. Thus, repulsive separation of the juxtacytoplasmic TM6 and TM3 in the M626I model conclusively demonstrates a direct link between the opening of this cytoplasmic face of the receptor structure and G protein coupling.     

Structural analysis of a HAMP domain: the linker region of the phototransducer in complex with sensory rhodopsin II

Sensory rhodopsin II, the photophobic receptor from Natronomonas pharaonis (NpSRII)5, forms a 2:2 complex with its cognate transducer (N. pharaonis halobacterial transducer of rhodopsins II (NpHtrII)) in lipid membranes. Light activation of NpSRII leads to a displacement of helix F, which in turn triggers a rotation/screw-like motion of TM2 in NpHtrII. This conformational change is thought to be transmitted through the membrane adjacent conserved signal transduction domain in histidine kinases, adenylyl cyclases, methyl-accepting chemotaxis proteins, and phosphatases (HAMP domain) to the cytoplasmic signaling domain of the transducer. The architecture and function of the HAMP domain are still unknown. In order to obtain information on the structure and dynamics of this region, EPR experiments on a truncated transducer (NpHtrII(157)) and NpSRII, site-directed spin-labeled and reconstituted into purple membrane lipids, have been carried out. A nitroxide scanning involving residues in the transducer helix TM2, in the predicted AS-1 region, and at selected positions in the following connector and AS-2 regions of the HAMP domain has been performed. Accessibility and dynamics data allowed us to identify a helical region up to residue Ala(94) in the AS-1 amphipathic sequence, followed by a highly dynamic domain protruding into the water phase. Additionally, transducer-transducer and transducer-receptor proximity relations revealed the overall architecture of the AS-1 sequences in the 2:2 complex, which are suggested to form a molten globular type of a coiled-coil bundle.