Structure Of Transmembrane Domain Research Articles

Discovery of small molecule binders of human FSHR(TMD) with novel structural scaffolds by integrating structural bioinformatics and machine learning algorithms

BackgroundThe activation of follicle stimulating hormone receptor (FSHR) by FSH and the consequent downstream signaling activities are crucial for reproductive health. The role of FSHR in tumor progression as well as osteoporosis advancement has also been well established. Currently, steroid preparations of estrogen and progesterone are being used for managing fertility, in spite of the harmful side effects, as there has not been much success in identification of effective FSHR modulators. Structure-based drug design initiatives for identification of potent and specific FSHR modulators have been impeded due to the non-availability of the complete crystal structure of hFSHR complexed with FSH. MethodsIn this study, we have modeled the 3D structure of transmembrane domain (TMD) of hFSHR and identified molecules that demonstrate good binding affinity by virtual screening of drug-like library of compounds. The 3D structural and pharmacophoric features of the binders and non-binders obtained from virtual screening were further used to develop Support Vector Machine based classifier for TMD binding. Based on the observations from docking and SVM classification, a small molecule was identified for extensive MD simulations and in vitro assays for FSHR modulatory activity. ResultsThe molecule selected based on docking score and SVM prediction was found to inhibit FSH-induced cAMP activity by 80% at 300 μM concentration. ConclusionThe study proposes 1,3-diphenyl-1H-pyrazole-5-carboxylate as a promising scaffold for the design of new and potent FSHR allosteric inhibitors.

The T Cell Antigen Receptor α Transmembrane Domain Coordinates Triggering through Regulation of Bilayer Immersion and CD3 Subunit Associations

Initial molecular details of cellular activation following αβT cell antigen receptor (TCR) ligation by peptide-major histocompatibility complexes (pMHC) remain unexplored. We determined the nuclear magnetic resonance (NMR) structure of the TCRα subunit transmembrane (TM) domain revealing a bipartite helix whose segmentation fosters dynamic movement. Positively charged TM residues Arg251 and Lys256 project from opposite faces of the helix, with Lys256 controlling immersion depth. Their modification caused stepwise reduction in TCR associations with CD3ζζ homodimers and CD3εγ plus CD3εδ heterodimers, respectively, leading to an activated transcriptome. Optical tweezers revealed that Arg251 and Lys256 mutations altered αβTCR-pMHC bond lifetimes, while mutations within interacting TCRα connecting peptide and CD3δ CxxC motif juxtamembrane elements selectively attenuated signal transduction. Our findings suggest that mechanical forces applied during pMHC ligation initiate Tcell activation via a dissociative mechanism, shifting disposition of those basic sidechains to rearrange TCR complex membrane topology and weaken TCRαβ and CD3 associations.

Dissecting conformational changes in APP's transmembrane domain linked to ε-efficiency in familial Alzheimer's disease.

The mechanism by which familial Alzheimer’s disease (FAD) mutations within the transmembrane domain (TMD) of the Amyloid Precursor Protein (APP) affect ε-endoproteolysis is only poorly understood. Thereby, mutations in the cleavage domain reduce ε-efficiency of γ-secretase cleavage and some even shift entry into production lines. Since cleavage occurs within the TMD, a relationship between processing and TMD structure and dynamics seems obvious. Using molecular dynamic simulations, we dissect the dynamic features of wild-type and seven FAD-mutants into local and global components. Mutations consistently enhance hydrogen-bond fluctuations upstream of the ε-cleavage sites but maintain strong helicity there. Dynamic perturbation-response scanning reveals that FAD-mutants target backbone motions utilized in the bound state. Those motions, obscured by large-scale motions in the pre-bound state, provide (i) a dynamic mechanism underlying the proposed coupling between binding and ε-cleavage, (ii) key sites consistent with experimentally determined docking sites, and (iii) the distinction between mutants and wild-type.

Bioinformatics analysis of ROP8 protein to improve vaccine design against Toxoplasma gondii

Rhoptry proteins (ROPs) are involved in the different stages of Toxoplasma gondii (T. gondii) invasion and are also critical for survival within host cells. ROP8 is expressed in the early stages of infection and have a key role in the parasitophorous vacuole (PV) formation. In this paper, we have combined several bioinformatics online servers for immunogenicity prediction of ROP8 protein. In this study, several bioinformatics approaches were used to analyze the different aspects of ROP8 protein, including the physico-chemical properties, transmembrane domain, subcellular localization, secondary and tertiary structure, B and T-cell potential epitopes, and other important characteristics of this protein. The findings showed that ROP8 protein had 60 potential post-translational modification sites. Also, only one transmembrane domain was recognized for this protein. The secondary structure of ROP8 protein comprises 33.04% alpha-helix, 18.26% extended strand, and 48.70% random coil. Moreover, several potential B and T-cell epitopes were identified for ROP8. In addition, the obtained findings from antigenicity and allergenicity evaluation remarked that this protein is immunogenic and non-allergen. Based on the results of Ramachandran plot, 94.8%, 4.1%, and 1.1% of amino acid residues were incorporated in the favored, allowed, and outlier regions, respectively. This paper provides a foundation for further investigations, and laid a theoretical basis for the development of an appropriate vaccine against toxoplasmosis. More studies are needed experimentally using the ROP8 alone or in combination with other antigens in the future.

Structure of APP-C991–99 and implications for role of extra-membrane domains in function and oligomerization

The 99 amino acid C-terminal fragment of Amyloid Precursor Protein APP-C99 (C99) is cleaved by γ-secretase to form Aβ peptide, which plays a critical role in the etiology of Alzheimer's Disease (AD). The structure of C99 consists of a single transmembrane domain flanked by intra and intercellular domains. While the structure of the transmembrane domain has been well characterized, little is known about the structure of the flanking domains and their role in C99 processing by γ-secretase. To gain insight into the structure of full-length C99, REMD simulations were performed for monomeric C99 in model membranes of varying thickness. We find equilibrium ensembles of C99 from simulation agree with experimentally-inferred residue insertion depths and protein backbone chemical shifts. In thin membranes, the transmembrane domain structure is correlated with extra-membrane structural states and the extra-membrane domain structural states become less correlated to each other. Mean and variance of the transmembrane and G37G38 hinge angles are found to increase with thinning membrane. The N-terminus of C99 forms β-strands that may seed aggregation of Aβ on the membrane surface, promoting amyloid formation. In thicker membranes the N-terminus forms α-helices that interact with the nicastrin domain of γ-secretase. The C-terminus of C99 becomes more α-helical as the membrane thickens, forming structures that may be suitable for binding by cytoplasmic proteins, while C-terminal residues essential to cytotoxic function become α-helical as the membrane thins. The heterogeneous but discrete extra-membrane domain states analyzed here open the path to new investigations of the role of C99 structure and membrane in amyloidogenesis. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

In with the Good and out with the Bad – The Role of SR‐B1 in Lowering Blood Cholesterol Levels

According to the CDC, 1 in 3 adults in the US has high cholesterol, which is linked to cardiovascular disease. Cholesterol, a hydrophobic molecule, is transported in the blood by binding a lipoprotein. Scavenger receptor class B type 1 (SR‐B1), a transmembrane protein that binds high density lipoprotein (HDL) carrying cholesterol, helps transfer cholesterol from HDL into liver cells where it can be excreted in bile. This is important for lowering blood cholesterol levels and reducing the risk of heart disease and stroke. SR‐B1 belongs to a larger family of scavenger receptor proteins that includes lysosome membrane protein 2 (LIMP‐2), and SR‐B1 and LIMP‐2 share structural homology. The structure of the C‐terminal transmembrane domain of SR‐B1, determined through NMR, is composed of three helices with a leucine zipper motif. Mutagenesis studies implicate this leucine zipper motif in dimerization of SR‐B1 and proper receptor function. A crystal structure of LIMP‐2 revealed the existence of a large, primarily hydrophobic cavity that runs the entire length of the protein that suggests a role in the selective transfer of cholesterol into the liver cell. Using 3D printing technology, the Cedarburg SMART (Students Modeling A Research Topic) Team utilized the known structure of LIMP‐2 to investigate structure‐function relationships in SR‐B1. Determining important HDL/SR‐B1 interactions will increase knowledge of how SR‐B1 functions in transferring cholesterol from plasma HDL to the liver for excretion and can lead to the development of medical treatments to increase cholesterol excretion by the liver and prevent heart disease and stroke.Support or Funding InformationThe MSOE Center for BioMolecular Modeling would like to acknowledge and thank the National Institutes of Health Science Education Partnership Award (NIH‐SEPA 1R25OD010505‐01) and the National Institutes of Health Clinical and Translational Science Award (NIH‐CTSA UL1RR031973) for their support in funding the 2017–2018 SMART Team Team program.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

A mechanism for agonist activation of the glucagon-like peptide-1 (GLP-1) receptor through modelling & molecular dynamics

The receptor for glucagon-like peptide 1 (GLP-1R) is a validated drug target for the treatment of type 2 diabetes and obesity. Recently the first three structures of GLP-1R were published – an X-ray structure of the apo transmembrane domain in the inactive conformation; an X-ray structure of the full-length receptor bound to a truncated peptide agonist; and a cryo-EM structure of the full-length receptor bound with GLP-1 and coupled to the G protein Gs. Since the inactive structure was incomplete, and the two active-state structures shared significant differences, we utilised all available knowledge to build hybrid models of the full length active and inactive state receptors. The two models were simulated using molecular dynamics and the output trajectories analysed and compared to reveal insights into the mechanism for agonist-mediated receptor activation. His-7, Glu-9 and Asp-15 of GLP-1 act together to destabilise transmembrane helix 6 and extracellular loop 3 in order to generate an active conformation of GLP-1R.

Partial gene identification and protein analysis of β-galactosidase from bird cherry-oat aphid

Bird cherry-oat aphid, Rhopalosiphum padi L. (Hemiptera: Aphididae) is a key pest of wheat and barley in the wide area of the world. Beta-galactosidase (EC 3.2.1.23) is an important digestive enzyme of R. padi which catalyses the hydrolysis of terminal, non-reducing beta-D-galactose residues in beta-D-galactosides. In this study a part of β-galactosidase (βgal) gene was isolated from R. padi (designated as Rp-βgal-JN835184), which contained 186 nucleotide coding and 62 amino acids. Homology search was done by BLAST in GenBank database. The identified protein sequence was aligned with 33 βgal protein sequences from various insects, two animals including human and mouse, two samples from rice and thale cress and two bacteria samples including Synechococcus sp. and Escherschia coli. Rp-βgal-JN835184 was analysed by computational tools to predict the protein properties, such as molecular mass, isoelectric point, signal peptide, transmembrane domain, secondary and spatial structure. Protein structure analysis revealed the deduced Rp-βgal-JN835184 had extensive homology with insect βgals and other organisms. Phylogenetic analysis of Rp-βgal-JN835184 showed that it has a close relationship with the bacteria, Synechococcus sp. Accordingly, βgals should be functional proteins involved in the biosynthesis of lactose and are derived from a common ancestor.

Functional characterization and substrate specificity analysis of Δ6-desaturase from marine microalga Isochrysis sp.

To express a Δ6-desaturase gene and produce gamma-linolenic acid (GLA) and stearidonic acid (SDA) in prokaryotic expression system (Escherichia coli), and analyze its substrate specificity in the omega-3 fatty acid biosynthetic pathway. Full-length ORF (1448bp) of Δ6Des-Iso was isolated from Isochrysis sp. and characterized using multiple sequence alignment, phylogenetic analysis, transmembrane domain, and protein tertiary structure. Δ6Des-Iso is a front-end desaturase consisting of three conserved histidine domains and a cytochrome b5 domain. Δ6Des-Iso was cloned and expressed in E. coli with the production of GLA and SDA. Recombinant E. coli utilized 27 and 8% of exogenously supplied alpha-linolenic acid (ALA) and linoleic acid (LA) to produce 6.3% of SDA and 2.3% of GLA, respectively, suggesting that isolated Δ6Des-Iso is specific to the omega-3 pathway. For the first time production of GLA and SDA in a prokaryotic system was achieved.

An Exhaustive Search Algorithm to Aid NMR-Based Structure Determination of Rotationally Symmetric Transmembrane Oligomers

Nuclear magnetic resonance (NMR) has been an important source of structural restraints for solving structures of oligomeric transmembrane domains (TMDs) of cell surface receptors and viral membrane proteins. In NMR studies, oligomers are assembled using inter-protomer distance restraints. But, for oligomers that are higher than dimer, these distance restraints all have two-fold directional ambiguity, and resolving such ambiguity often requires time-consuming trial-and-error calculations using restrained molecular dynamics (MD) with simulated annealing (SA). We report an Exhaustive Search algorithm for Symmetric Oligomer (ExSSO), which can perform near-complete search of the symmetric conformational space in a very short time. In this approach, the predetermined protomer model is subject to full angular and spatial search within the symmetry space. This approach, which can be applied to any rotationally symmetric oligomers, was validated using the structures of the Fas death receptor, the HIV-1 gp41 fusion protein, the influenza proton channel, and the MCU pore. The algorithm is able to generate approximate oligomer solutions quickly as initial inputs for further refinement using the MD/SA method.

Assessment of the transmembrane domain structures in GPCR Dock 2013 models.

The community-wide blind prediction of G-protein coupled receptor (GPCR) structures and ligand docking has been conducted three times and the quality of the models was primarily assessed by the accuracy of ligand binding modes. The seven transmembrane (TM) helices of the receptors were taken as a whole; thus the model quality within the 7TM domains has not been evaluated. Here we evaluate the 7TM domain structures in the models submitted for the last round of prediction - GPCR Dock 2013. Applying the 7 × 7 RMSD matrix analysis described in our prior work, we show that the models vary widely in prediction accuracy of the 7TM structures, exhibiting diverse structural differences from the targets. For the prediction of the 5-hydroxytryptamine receptors, the top 7TM models are rather close to the targets, which however are not ranked top by ligand-docking. On the other hand, notable deviations of the TMs are found in in the previously identified top docking models that closely resemble other receptors. We further reveal reasons of success and failure in ligand docking for the models. This current assessment not only complements the previous assessment, but also provides important insights into the current status of GPCR modeling and ligand docking.

The L530R variation associated with recurrent kidney stones impairs the structure and function of TRPV5

TRPV5 is a Ca2+-selective channel that plays a key role in the reabsorption of Ca2+ ions in the kidney. Recently, a rare L530R variation (rs757494578) of TRPV5 was found to be associated with recurrent kidney stones in a founder population. However, it was unclear to what extent this variation alters the structure and function of TRPV5. To evaluate the function and expression of the TRPV5 variant, Ca2+ uptake in Xenopus oocytes and western blot analysis were performed. The L530R variation abolished the Ca2+ uptake activity of TRPV5 in Xenopus oocytes. The variant protein was expressed with drastic reduction in complex glycosylation. To assess the structural effects of this L530R variation, TRPV5 was modeled based on the crystal structure of TRPV6 and molecular dynamics simulations were carried out. Simulation results showed that the L530R variation disrupts the hydrophobic interaction between L530 and L502, damaging the secondary structure of transmembrane domain 5. The variation also alters its interaction with membrane lipid molecules. Compared to the electroneutral L530, the positively charged R530 residue shifts the surface electrostatic potential towards positive. R530 is attracted to the negatively charged phosphate group rather than the hydrophobic carbon atoms of membrane lipids. This shifts the pore helix where R530 is located and the D542 residue in the Ca2+-selective filter towards the surface of the membrane. These alterations may lead to misfolding of TRPV5, reduction in translocation of the channel to the plasma membrane and/or impaired Ca2+ transport function of the channel, and ultimately disrupt TRPV5-mediated Ca2+ reabsorption.

Spatial structure of TLR4 transmembrane domain in bicelles provides the insight into the receptor activation mechanism

Toll-like receptors (TLRs) play a key role in the innate and adaptive immune systems. While a lot of structural data is available for the extracellular and cytoplasmic domains of TLRs, and a model of the dimeric full-length TLR3 receptor in the active state was build, the conformation of the transmembrane (TM) domain and juxtamembrane regions in TLR dimers is still unclear. In the present work, we study the transmembrane and juxtamembrane parts of human TLR4 receptor using solution NMR spectroscopy in a variety of membrane mimetics, including phospholipid bicelles. We show that the juxtamembrane hydrophobic region of TLR4 includes a part of long TM α-helix. We report the dimerization interface of the TM domain and claim that long TM domains with transmembrane charged aminoacids is a common feature of human toll-like receptors. This fact is analyzed from the viewpoint of protein activation mechanism, and a model of full-length TLR4 receptor in the dimeric state has been proposed.

Critical Role of Water Molecules in Proton Translocation by the Membrane-Bound Transhydrogenase

The nicotinamide nucleotide transhydrogenase (TH) is an integral membrane enzyme that uses the proton-motive force to drive hydride transfer from NADH to NADP+ in bacteria and eukaryotes. Here we solved a 2.2-Å crystal structure of the TH transmembrane domain (Thermus thermophilus) at pH 6.5. This structure exhibits conformational changes of helix positions from a previous structure solved at pH 8.5, and reveals internal water molecules interacting with residues implicated in proton translocation. Together with molecular dynamics simulations, we show that transient water flows across anarrow pore and a hydrophobic "dry" region inthe middle of the membrane channel, with key residues His42α2 (chain A) being protonated and Thr214β (chain B) displaying a conformational change, respectively, to gate the channel access to both cytoplasmic and periplasmic chambers. Mutation of Thr214β to Ala deactivated the enzyme. These data provide new insights into the gating mechanism of proton translocation in TH.

Human GLP-1 receptor transmembrane domain structure in complex with allosteric modulators.

The glucagon-like peptide-1 receptor (GLP-1R) and the glucagon receptor (GCGR) are members of the secretin-like class B family of G-protein-coupled receptors (GPCRs) and have opposing physiological roles in insulin release and glucose homeostasis. The treatment of type 2 diabetes requires positive modulation of GLP-1R to inhibit glucagon secretion and stimulate insulin secretion in a glucose-dependent manner. Here we report crystal structures of the human GLP-1R transmembrane domain in complex with two different negative allosteric modulators, PF-06372222 and NNC0640, at 2.7 and 3.0 Å resolution, respectively. The structures reveal a common binding pocket for negative allosteric modulators, present in both GLP-1R and GCGR and located outside helices V-VII near the intracellular half of the receptor. The receptor is in an inactive conformation with compounds that restrict movement of the intracellular tip of helix VI, a movement that is generally associated with activation mechanisms in class A GPCRs. Molecular modelling and mutagenesis studies indicate that agonist positive allosteric modulators target the same general region, but in a distinct sub-pocket at the interface between helices V and VI, which may facilitate the formation of an intracellular binding site that enhances G-protein coupling.

Ligand-guided homology modelling of the GABAB2 subunit of the GABAB receptor.

γ-aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the central nervous system, and disturbances in the GABAergic system have been implicated in numerous neurological and neuropsychiatric diseases. The GABAB receptor is a heterodimeric class C G protein-coupled receptor (GPCR) consisting of GABAB1a/b and GABAB2 subunits. Two GABAB receptor ligand binding sites have been described, namely the orthosteric GABA binding site located in the extracellular GABAB1 Venus fly trap domain and the allosteric binding site found in the GABAB2 transmembrane domain. To date, the only experimentally solved three-dimensional structures of the GABAB receptor are of the Venus fly trap domain. GABAB receptor allosteric modulators, however, show great therapeutic potential, and elucidating the structure of the GABAB2 transmembrane domain may lead to development of novel drugs and increased understanding of the allosteric mechanism of action. Despite the lack of x-ray crystal structures of the GABAB2 transmembrane domain, multiple crystal structures belonging to other classes of GPCRs than class A have been released within the last years. More closely related template structures are now available for homology modelling of the GABAB receptor. Here, multiple homology models of the GABAB2 subunit of the GABAB receptor have been constructed using templates from class A, B and C GPCRs, and docking of five clusters of positive allosteric modulators and decoys has been undertaken to select models that enrich the active compounds. Using this ligand-guided approach, eight GABAB2 homology models have been chosen as possible structural representatives of the transmembrane domain of the GABAB2 subunit. To the best of our knowledge, the present study is the first to describe homology modelling of the transmembrane domain of the GABAB2 subunit and the docking of positive allosteric modulators in the receptor.

Solid-State Structure of Abeta (Aβ) in Alzheimer's Disease.

Alzheimer's disease (AD) has become the most common neurodegenerative disease. The deposition of amyloid fibrils in the brain is one of the characteristics of AD. The fibrils are composed of amyloid-β peptide (Aβ). Aβ is produced through a series event of protease cleavage of a transmembrane protein called β-amyloid precursor protein (APP) which is commonly expressed in the brain. The production of Aβ and its propensity to aggregation to form oligomers and fibrils are believed to initiate a sequence of events that lead to AD dementia. The production of Aβ is influenced by the transmembrane domain (TM) structure of APP. The structure variety of different Aβ assemblies including oligomers and fibrils may result in different neurotoxicity to the brain. Therefore, enormous work has been carried out to study the structure of APP TM and various Aβ assemblies. Solid-state NMR has advantages in studying immobile protein structures with large molecular weight. In this review, solid-state NMR structure of APP TM and different Aβ assemblies will be discussed, especially various Aβ amyloid fibril structures. This structural information greatly enhanced our understanding in AD, providing fundamental knowledge that will help in finding a treatment for AD.

7×7 RMSD matrix: A new method for quantitative comparison of the transmembrane domain structures in the G-protein coupled receptors.

The G-protein coupled receptors (GPCRs) share a conserved heptahelical fold in the transmembrane (TM) region, but the exact arrangements of the seven TM helices vary with receptors and their activation states. The differences or the changes have been observed in the experimentally solved structures, but have not been systematically and quantitatively investigated due to lack of suitable methods. In this work, we describe a novel method, called 7×7 RMSD matrix that is proposed specifically for comparing the characteristic 7TM bundle structures of GPCRs. Compared to the commonly used overall TM bundle RMSD as a single parameter, a 7×7 RMSD matrix contains 49 parameters, which reveal changes of the relative orientations of the seven TMs. We demonstrate the novelty and advantages of this method by tackling two problems that are challenging for the existing methods. With this method, we are able to identify and quantify the helix movements in the activated receptor structures and reveal structural conservation and divergence as well as the structural relationships of different GPCRs in terms of the relative orientations of the seven TMs.

NMR Structure of the C-Terminal Transmembrane Domain of the HDL Receptor, SR-BI, and a Functionally Relevant Leucine Zipper Motif

The interaction of high-density lipoprotein (HDL) with its receptor, scavenger receptor BI (SR-BI), is critical for lowering plasma cholesterol levels and reducing the risk for cardiovascular disease. The HDL/SR-BI complex facilitates delivery of cholesterol into cells and is likely mediated by receptor dimerization. This work describes the use of nuclear magnetic resonance (NMR) spectroscopy to generate the first high-resolution structure of the C-terminal transmembrane domain of SR-BI. This region of SR-BI harbors a leucine zipper dimerization motif, which when mutated impairs the ability of the receptor to bind HDL and mediate cholesterol delivery. These losses in function correlate with the inability of SR-BI to form dimers. We also identify juxtamembrane regions of the extracellular domain of SR-BI that may interact with the lipid surface to facilitate cholesterol transport functions of the receptor.

In Situ Solid-State Nmr Spectroscopy of APP Transmembrane Domain (TM) Structure and Dimerization in Native E. Coli Membranes

Cleavage of the amyloid precursor protein (APP) by α-, β- and γ-secretases leads to generation of amyloid peptide (Aβ), which is the primary constituent of Alzheimer's disease (AD) plaques in the brain. Three consecutive GXXXG motifs and one GXXXA motif on the transmembrane domain (TM) of APP have been reported to be one or two important sites for mediating APP homodimerization in different experimental conditions, the difference are likely due to the membrane mimetic. We and several other groups recently demonstrated the feasibility to directly characterize recombinant proteins in native Escherichia coli (E. coli) membranes using magic-angle spinning (MAS) solid-state nuclear magnetic resonance (ssNMR)(2, 3). Here, we will present several strategies to further improve spectral sensitivity and resolution, and suppress background signals of E. coli proteins and lipids, including preparation of E. coli inner membranes, reverse 13C- and 15N-isotope labeling, and selective signal filtering with the frequency-selective REDOR pulse sequence. These improvements have allowed us to characterize the structure and dimerization of APP TM and its variants in a native membrane environment.