Building protein diagrams on the web with the residue-based diagram editor RbDe

Mapping the Nonstructural Transmembrane Proteins of Severe Acute Respiratory Syndrome Coronavirus 2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) responsible for the disease coronavirus-19 disease (COVID-19) has wreaked havoc on the health and economy of humanity. In addition, the disease is observed in domestic and wild animals. The disease has impacted directly and indirectly every corner of the planet. Currently, there are no effective therapies for the treatment of COVID-19. Vaccination to protect against COVID-19 started in December 2020. SARS-CoV-2 is an enveloped virus with a single-stranded RNA genome of 29.8 kb. More than two-thirds of the genome comprise Orf1ab encoding 16 nonstructural proteins (nsps) followed by mRNAs encoding structural proteins, spike (S), envelop (E), membrane (M), and nucleocapsid (N). These genes are interspaced with several accessory genes (open reading frames [Orfs] 3a, 3b, 6, 7a, 7b, 8, 9b, 9c, and 10). The functions of these proteins are of particular interest for understanding the pathogenesis of SARS-CoV-2. Several of the nsps (nsp3, nsp4, and nsp6) and Orf3a are transmembrane proteins involved in regulating the host immunity, modifying host cell organelles for viral replication and escape and hence considered drug targets. In this paper, we report mapping the transmembrane structure of the nsps of SARS-CoV-2.

The Structure of the Membrane Protein of SARS-CoV-2 Resembles the Sugar Transporter SemiSWEET

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the disease COVID-19 that has decimated the health and economy of our planet. The virus causes the disease not only in people but also in companion and wild animals. People with diabetes are at risk of the disease. As yet we do not know why the virus has been highly successful in causing the pandemic within 3 months of its first report. The structural proteins of SARS include membrane glycoprotein (M), envelope protein (E), nucleocapsid protein (N), and the spike protein (S).

METHODS

The structure and function of the most abundant structural protein of SARS-CoV-2, the membrane (M) glycoprotein, is not fully understood. Using in silico analyses we determined the structure and potential function of the M protein.

RESULTS

The M protein of SARS-CoV-2 is 98.6% similar to the M protein of bat SARS-CoV, maintains 98.2% homology with pangolin SARS-CoV, and has 90% homology with the M protein of SARS-CoV; whereas, the similarity is only 38% with the M protein of MERS-CoV. In silico analyses showed that the M protein of SARS-CoV-2 has a triple helix bundle, forms a single 3-trans-membrane domain, and is homologous to the prokaryotic sugar transport protein SemiSWEET. SemiSWEETs are related to the PQ-loop family whose members function as cargo receptors in vesicle transport, mediate movement of basic amino acids across lysosomal membranes, and are also involved in phospholipase flippase function.

CONCLUSIONS

The advantage and role of the M protein having a sugar transporter-like structure is not clearly understood. The M protein of SARS-CoV-2 interacts with S, E, and N protein. The S protein of the virus is glycosylated. It could be hypothesized that the sugar transporter-like structure of the M protein influences glycosylation of the S protein. Endocytosis is critical for the internalization and maturation of RNA viruses, including SARS-CoV-2. Sucrose is involved in endosome and lysosome maturation and may also induce autophagy, pathways that help in the entry of the virus. Overall, it could be hypothesized that the SemiSWEET sugar transporter-like structure of the M protein may be involved in multiple functions that may aid in the rapid proliferation, replication, and immune evasion of the SARS-CoV-2 virus. Biological experiments would validate the presence and function of the SemiSWEET sugar transporter.

Tachykinin-1 receptor antagonism suppresses substance-P- and compound 48/80-induced mast cell activation from rat mast cells expressing functional mas-related GPCR B3

OBJECTIVE

Mice and rats are important animal models for mast cell (MC) study. However, rat Mas-related-GPCR-B3 receptor (MRGPRB3) has been less studied than its mouse counterpart. Therefore, we aimed to characterize rat MRGPRB3.

METHODS

Mrgprb3 mRNA expression was assessed in peritoneal cells (RPCs) and peritoneal MCs (RPMCs) of wild-type rats, RPCs of MC-deficient rats, and RBL-2H3 cells by reverse-transcriptase polymerase chain reaction (RT-PCR). RPMCs, MRGPRX2-transfected and non-transfected RBL-2H3 cells were activated by 15-30 min incubation with DNP-BSA, substance-P (SP), or compound-48/80. L732138 or CP96344 was used as a tachykinin/neurokinin-1-receptor antagonist. Histamine release from MCs was measured by HPLC fluorometry.

RESULTS

Mrgprb3 mRNA expression was found in all cells, with the highest level in wild-type RPCs. All cells responded to DNP-BSA, but only MRGPRX2-transfected-RBL-2H3 cells and RPMCs responded to all activators. L732138 (0.1-10 μM) and CP96344 (1-100 μM) suppressed SP (10 μM)-induced RPMC activation. L732138 inhibition was dose independent, whereas CP96344 inhibition occurred in a dose-dependent manner. Additionally, only CP96344 suppressed SP (100 μM)- and compound-48/80 (10 μg/mL)-induced RPMC activation.

CONCLUSIONS

RPMCs expressing functional MRGPRB3 response upon MRGPRX2 ligands to regulated MC-mediated activities. It`s provide novel insights for future pseudo-allergic studies in rodents.

Structural Mapping and Functional Characterization of Zebrafish Class B G-Protein Coupled Receptor (GPCR) with Dual Ligand Selectivity towards GLP-1 and Glucagon

GLP-1 and glucagon regulate glucose metabolism through a network of metabolic pathways initiated upon binding to their specific receptors that belong to class B G-protein coupled receptors (GPCRs). The therapeutic potential of glucagon is currently being evaluated, while GLP-1 is already used in the treatment of type 2 diabetes and obesity. Development of a second generation of GLP-1 based therapeutics depends on a molecular and structural understanding of the interactions between the GLP-1 receptor (GLP-1R) and its ligand GLP-1. There is considerable sequence conservation between GLP-1 and glucagon and between the hGLP-1R and human glucagon receptor (hGCGR), yet each receptor recognizes only its own specific ligand. Glucagon receptors in fish and frogs also exhibit ligand selectivity only towards glucagon and not GLP-1. Based on competitive binding experiments and assays of increase in intracellular cAMP, we demonstrate here that a GPCR in zebrafish (Danio rerio) exhibits dual ligand selectivity towards GLP-1 and glucagon, a characteristic not found in mammals. Further, many structural features found in hGLP-1R and hGCGR are also found in this zebrafish GPCR (zfGPCR). We show this by mapping of its sequence and structural features onto the hGLP-1R and hGCGR based on their partial and complementary crystal structures. Thus, we propose that zfGPCR represents a dual GLP-1R/GCGR. The main differences between the three receptors are in their stalk regions that connect their N-terminal extracellular domains (NECDs) with their transmembrane domains and the absence of loop 3 in the NECD in zfGLP-1R/GCGR. These observations suggest that the interactions between GLP-1 and glucagon with loop 3 and the stalk regions may induce different conformational changes in hGLP-1R and hGCGR upon ligand binding and activation that lead to selective recognition of their native ligands.

The sodium/multivitamin transporter: a multipotent system with therapeutic implications

The Na(+)/multivitamin transporter (SMVT) is a member of the solute:sodium symporter family that catalyzes the Na(+)-dependent uptake of the structurally diverse water-soluble vitamins pantothenic acid (vitamin B5) and biotin (vitamin H), α-lipoic acid-a vitamin-like substance with strong antioxidant properties-and iodide. The organic substrates of SMVT play central roles in the cellular metabolism and are, therefore, essential for normal human health and development. For example, biotin deficiency leads to growth retardation, dermatological disorders, and neurological disorders. Animal studies have shown that biotin deficiency during pregnancy is directly correlated to embryonic growth retardation, congenital malformation, and death of the embryo. This chapter focuses on the structural and functional features of the human isoform of SMVT (hSMVT); the discovery of which was greatly facilitated by the cloning and expression of hSMVT in tractable expression systems. Special emphasis will be given to mechanistic implications of the transport process of hSMVT that will inform our understanding of the molecular determinants of hSMVT-mediated transport in dynamic context to alleviate the development and optimization of hSMVT as a multipotent platform for drug delivery.

Molecular evolutionary characterization of a V1R subfamily unique to strepsirrhine primates

Vomeronasal receptor genes have frequently been invoked as integral to the establishment and maintenance of species boundaries among mammals due to the elaborate one-to-one correspondence between semiochemical signals and neuronal sensory inputs. Here, we report the most extensive sample of vomeronasal receptor class 1 (V1R) sequences ever generated for a diverse yet phylogenetically coherent group of mammals, the tooth-combed primates (suborder Strepsirrhini). Phylogenetic analysis confirms our intensive sampling from a single V1R subfamily, apparently unique to the strepsirrhine primates. We designate this subfamily as V1Rstrep. The subfamily retains extensive repertoires of gene copies that descend from an ancestral gene duplication that appears to have occurred prior to the diversification of all lemuriform primates excluding the basal genus Daubentonia (the aye-aye). We refer to the descendent clades as V1Rstrep-α and V1Rstrep-β. Comparison of the two clades reveals different amino acid compositions corresponding to the predicted ligand-binding site and thus potentially to altered functional profiles between the two. In agreement with previous studies of the mouse lemur (genus, Microcebus), the majority of V1Rstrep gene copies appear to be intact and under strong positive selection, particularly within transmembrane regions. Finally, despite the surprisingly high number of gene copies identified in this study, it is nonetheless probable that V1R diversity remains underestimated in these nonmodel primates and that complete characterization will be limited until high-coverage assembled genomes are available.

Protter: interactive protein feature visualization and integration with experimental proteomic data

SUMMARY

The ability to integrate and visualize experimental proteomic evidence in the context of rich protein feature annotations represents an unmet need of the proteomics community. Here we present Protter, a web-based tool that supports interactive protein data analysis and hypothesis generation by visualizing both annotated sequence features and experimental proteomic data in the context of protein topology. Protter supports numerous proteomic file formats and automatically integrates a variety of reference protein annotation sources, which can be readily extended via modular plug-ins. A built-in export function produces publication-quality customized protein illustrations, also for large datasets. Visualizations of surfaceome datasets show the specific utility of Protter for the integrated visual analysis of membrane proteins and peptide selection for targeted proteomics.

AVAILABILITY AND IMPLEMENTATION

The Protter web application is available at http://wlab.ethz.ch/protter. Source code and installation instructions are available at http://ulo.github.io/Protter/.

CONTACT

wbernd@ethz.ch

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cholecystokinin receptors in Atlantic salmon: molecular cloning, gene expression, and structural basis

The peptide hormone cholecystokinin (CCK) exerts a wide range of digestive and CNS-related physiological signaling via CCK receptors in brain and gut. There is very limited information available on these receptors in Atlantic salmon. The aim of this study was to characterize CCK receptors in gut and brain of salmon. We have identified and cloned one CCK-1 receptor and duplicates of CCK-2 receptor in salmon. The phylogenetic analysis indicates the existence of one common ancestor gene for all CCK receptors. CCK-1R mRNA is highly expressed in pancreas followed by midgut, hindgut, gallbladder, and stomach indicating an involvement in pancreatic regulation and gallbladder contractions. CCK-2R1/gastrin mRNA is expressed at high levels in midgut and at relatively low levels in stomach, gallbladder, and pancreas. We postulate CCK-2R1/gastrin receptor to have gastrin-related functions because of its distribution and abundance in gastro-intestinal (GI) tissues. CCK-2R2 is relatively abundant in brain but has low expression levels in gut tissues supporting the hypothesis for involvement in the gut-brain signaling. Major functional motifs and ligand interaction sites in salmon are conserved with that of mammals. This information will be instrumental for comparative studies and further targeting receptor activation and selectivity of biological responses of CCK in salmon.

Residues within the transmembrane domain of the glucagon-like peptide-1 receptor involved in ligand binding and receptor activation: modelling the ligand-bound receptor

The C-terminal regions of glucagon-like peptide-1 (GLP-1) bind to the N terminus of the GLP-1 receptor (GLP-1R), facilitating interaction of the ligand N terminus with the receptor transmembrane domain. In contrast, the agonist exendin-4 relies less on the transmembrane domain, and truncated antagonist analogs (e.g. exendin 9-39) may interact solely with the receptor N terminus. Here we used mutagenesis to explore the role of residues highly conserved in the predicted transmembrane helices of mammalian GLP-1Rs and conserved in family B G protein coupled receptors in ligand binding and GLP-1R activation. By iteration using information from the mutagenesis, along with the available crystal structure of the receptor N terminus and a model of the active opsin transmembrane domain, we developed a structural receptor model with GLP-1 bound and used this to better understand consequences of mutations. Mutation at Y152 [transmembrane helix (TM) 1], R190 (TM2), Y235 (TM3), H363 (TM6), and E364 (TM6) produced similar reductions in affinity for GLP-1 and exendin 9-39. In contrast, other mutations either preferentially [K197 (TM2), Q234 (TM3), and W284 (extracellular loop 2)] or solely [D198 (TM2) and R310 (TM5)] reduced GLP-1 affinity. Reduced agonist affinity was always associated with reduced potency. However, reductions in potency exceeded reductions in agonist affinity for K197A, W284A, and R310A, while H363A was uncoupled from cAMP generation, highlighting critical roles of these residues in translating binding to activation. Data show important roles in ligand binding and receptor activation of conserved residues within the transmembrane domain of the GLP-1R. The receptor structural model provides insight into the roles of these residues.

Evidence for positive selection in the extracellular domain of human cytomegalovirus encoded G protein-coupled receptor US28

The human cytomegalovirus (HCMV)-encoded chemokine receptor US28 is also a seven-transmembrane G-protein coupled receptor, whose signaling pathway is known for its involvement in host immune system evasion. HCMV infection can result in serious disease in immunocompromised individuals and is also linked to atherosclerosis and cardiovascular disease. Identifying amino acid residues that play a crucial role in successful viral adaptation in response to the host's immune defense is critical for effective drug design. In this study maximum likelihood-based codon substitution analyses were carried out to determine whether any codon of US28 has evolved adaptively. If the rate of nonsynonymous (dn) to the rate of synonymous (ds) nucleotide substitutions (ω = dn/ds) is greater than one, the codon is said to be under positive selection, indicating adaptive evolution. Although the overall ω for US28 gene was 0.154, indicating that most codon sites were subject to strong purifying selection, five codon sites are under strong positive selection. Three (E18D/L, D19A/E/G, and R267K/Q) of these positively selected sites are located in extracellular domains, the domains that play a crucial role for successful viral adaptation in response to the host's immune defense. The C-terminal (R329Q/W) and the fifth transmembrane domain (V190I), each have one positively selected site. These results suggest that relative to the extracellular domains, amino acid residues present in intracellular domains are more selectively constrained. A few amino acid residues in extracellular domains of US28 evolved more rapidly, presumably due to positive selection pressure resulting from ligand-binding and pathogen interactions of extracellular domains.

Local and global dynamics of the G protein-coupled receptor CXCR1

The local and global dynamics of the chemokine receptor CXCR1 are characterized using a combination of solution NMR and solid-state NMR experiments. In isotropic bicelles (q = 0.1), only 13% of the expected number of backbone amide resonances is observed in (1)H/(15)N HSQC solution NMR spectra of uniformly (15)N-labeled samples; extensive deuteration and the use of TROSY made little difference in the 800 MHz spectra. The limited number of observed amide signals is ascribed to mobile backbone sites and assigned to specific residues in the protein; 19 of the signals are from residues at the N-terminus and 25 from residues at the C-terminus. The solution NMR spectra display no evidence of local backbone motions from residues in the transmembrane helices or interhelical loops of CXCR1. This finding is reinforced by comparisons of solid-state NMR spectra of both magnetically aligned and unoriented bilayers containing either full-length or doubly N- and C-terminal truncated CXCR1 constructs. CXCR1 undergoes rapid rotational diffusion about the normal of liquid crystalline phospholipid bilayers; reductions in the frequency span and a change to axial symmetry are observed for both carbonyl carbon and amide nitrogen chemical shift powder patterns of unoriented samples containing (13)C- and (15)N-labeled CXCR1. In contrast, when the phospholipids are in the gel phase, CXCR1 does not undergo rapid global reorientation on the 10(4) Hz time scale defined by the carbonyl carbon and amide nitrogen chemical shift powder patterns.

Oligomeric size of the m2 muscarinic receptor in live cells as determined by quantitative fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET), measured by fluorescence intensity-based microscopy and fluorescence lifetime imaging, has been used to estimate the size of oligomers formed by the M(2) muscarinic cholinergic receptor. The approach is based on the relationship between the apparent FRET efficiency within an oligomer of specified size (n) and the pairwise FRET efficiency between a single donor and a single acceptor (E). The M(2) receptor was fused at the N terminus to enhanced green or yellow fluorescent protein and expressed in Chinese hamster ovary cells. Emission spectra were analyzed by spectral deconvolution, and apparent efficiencies were estimated by donor-dequenching and acceptor-sensitized emission at different ratios of enhanced yellow fluorescent protein-M(2) receptor to enhanced green fluorescent protein-M(2) receptor. The data were interpreted in terms of a model that considers all combinations of donor and acceptor within a specified oligomer to obtain fitted values of E as follows: n = 2, 0.495 +/- 0.019; n = 4, 0.202 +/- 0.010; n = 6, 0.128 +/- 0.006; n = 8, 0.093 +/- 0.005. The pairwise FRET efficiency determined independently by fluorescence lifetime imaging was 0.20-0.24, identifying the M(2) receptor as a tetramer. The strategy described here yields an explicit estimate of oligomeric size on the basis of fluorescence properties alone. Its broader application could resolve the general question of whether G protein-coupled receptors exist as dimers or larger oligomers. The size of an oligomer has functional implications, and such information can be expected to contribute to an understanding of the signaling process.

Evidence of positive selection at codon sites localized in extracellular domains of mammalian CC motif chemokine receptor proteins

Background

CC chemokine receptor proteins (CCR1 through CCR10) are seven-transmembrane G-protein coupled receptors whose signaling pathways are known for their important roles coordinating immune system responses through targeted trafficking of white blood cells. In addition, some of these receptors have been identified as fusion proteins for viral pathogens: for example, HIV-1 strains utilize CCR5, CCR2 and CCR3 proteins to obtain cellular entry in humans. The extracellular domains of these receptor proteins are involved in ligand-binding specificity as well as pathogen recognition interactions.

In mammals, the majority of chemokine receptor genes are clustered together; in humans, seven of the ten genes are clustered in the 3p21-24 chromosome region. Gene conversion events, or exchange of DNA sequence between genes, have been reported in chemokine receptor paralogs in various mammalian lineages, especially between the cytogenetically closely located pairs CCR2/5 and CCR1/3. Datasets of mammalian orthologs for each gene were analyzed separately to minimize the potential confounding impact of analyzing highly similar sequences resulting from gene conversion events.

Molecular evolution approaches and the software package Phylogenetic Analyses by Maximum Likelihood (PAML) were utilized to investigate the signature of selection that has acted on the mammalian CC chemokine receptor (CCR) gene family. The results of neutral vs. adaptive evolution (positive selection) hypothesis testing using Site Models are reported. In general, positive selection is defined by a ratio of nonsynonymous/synonymous nucleotide changes (dN/dS, or ω) >1.

Results

Of the ten mammalian CC motif chemokine receptor sequence datasets analyzed, only CCR2 and CCR3 contain amino acid codon sites that exhibit evidence of positive selection using site based hypothesis testing in PAML. Nineteen of the twenty codon sites putatively indentified as likely to be under positive selection code for amino acid residues located in extracellular domains of the receptor protein products.

Conclusions

These results suggest that amino acid residues present in intracellular and membrane-bound domains are more selectively constrained for functional signal transduction and homo- or heterodimerization, whereas amino acid residues in extracellular domains of these receptor proteins evolve more quickly, perhaps due to heightened selective pressure resulting from ligand-binding and pathogen interactions of extracellular domains.

Specific domains in anterior pharynx-defective 1 determine its intramembrane interactions with nicastrin and presenilin

γ-Secretase, a multisubunit transmembrane protease comprised of presenilin, nicastrin, presenilin enhancer 2, and anterior pharynx-defective one, participates in the regulated intramembrane proteolysis of Type I membrane proteins including the amyloid precursor protein (APP). Although Aph-1 is thought to play a structural role in the assembly of γ-secretase complex and several transmembrane domains (TMDs) of Aph-1 have been shown to be critical for its function, the importance of the other domains of Aph-1 remains elusive. We screened a series of Aph-1 mutants and focused on nine mutations distributed in six different TMDs of human APH-1aS, assessing their ability to complement mouse embryonic fibroblasts lacking Aph-1. We showed that mutations in TMD4 (G126) and TMD5 (H171) of Aph-1aS prevented the formation of the Nct/Aph-1 subcomplex. Importantly, although mutations in TMD3 (Q83/E84/R85) and TMD6 (H197) of APH-1aS did not affect Nct/Aph-1 subcomplex formation, both mutations prevented further association/endoproteolysis of PS1. We propose a model that identifies critical TMDs of Aph-1 for associations with Nct and PS for the stepwise assembly of γ-secretase components.

Familial glucocorticoid deficiency type 1 due to a novel compound heterozygous MC2R mutation

OBJECTIVE

Description of the clinical, biochemical and genetic features of a Polish patient with familial glucocorticoid deficiency.

METHODS

Detailed clinical investigation, hormonal analysis and sequencing of the coding region of the melanocortin 2 receptor (MC2R) gene in this patient.

RESULTS

We report on a 3-month-old boy with familial glucocorticoid deficiency who presented at the age of 3 months with skin hyperpigmentation, muscle weakness, mild jaundice and constipation. Hormonal analyses revealed high ACTH and TSH serum concentrations, low serum cortisol concentration along with normal blood electrolytes. On hydrocortisone supplementation, the disease symptoms disappeared and the child recovered completely. His physical and mental development progresses normally. Genetic analysis disclosed a novel compound heterozygous MC2R mutation p.Leu46fs and p.Val49Met.

CONCLUSION

The heterozygous p.Leu46fs mutation adds to the small number of MC2R nonsense mutations and is the first frameshift mutation within the first transmembrane domain of the receptor. According to molecular modeling the Val49Met mutation results in a structural change of the first transmembrane domain and in a potential novel interaction of the transmembrane domains I and VII.

High-resolution NMR spectroscopy of a GPCR in aligned bicelles

Solid-state NMR spectra with single-site resolution of CXCR1, a G protein-coupled receptor (GPCR), were obtained in magnetically aligned phospholipid bicelles. These results demonstrate that GPCRs in phospholipid bilayers are suitable samples for structure determination by solid-state NMR. The spectra also enable studies of drug-receptor interactions.

Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo

Drosophila olfactory sensory neurons (OSNs) each express two odorant receptors (ORs): a divergent member of the OR family and the highly conserved, broadly expressed receptor OR83b. OR83b is essential for olfaction in vivo and enhances OR function in vitro, but the molecular mechanism by which it acts is unknown. Here we demonstrate that OR83b heterodimerizes with conventional ORs early in the endomembrane system in OSNs, couples these complexes to the conserved ciliary trafficking pathway, and is essential to maintain the OR/OR83b complex within the sensory cilia, where odor signal transduction occurs. The OR/OR83b complex is necessary and sufficient to promote functional reconstitution of odor-evoked signaling in sensory neurons that normally respond only to carbon dioxide. Unexpectedly, unlike all known vertebrate and nematode chemosensory receptors, we find that Drosophila ORs and OR83b adopt a novel membrane topology with their N-termini and the most conserved loops in the cytoplasm. These loops mediate direct association of ORs with OR83b. Our results reveal that OR83b is a universal and integral part of the functional OR in Drosophila. This atypical heteromeric and topological design appears to be an insect-specific solution for odor recognition, making the OR/OR83b complex an attractive target for the development of highly selective insect repellents to disrupt olfactory-mediated host-seeking behaviors of insect disease vectors.

This study reveals a novel membrane topology for olfactory receptors in Drosophila and details the molecular mechanisms of receptor localization at the sensory cilia.

Adaptive evolution of MRGX2, a human sensory neuron specific gene involved in nociception

MRGX2, a G-protein-coupled receptor, is specifically expressed in the sensory neurons of the human peripheral nervous system and involved in nociception. Here, we studied DNA polymorphism patterns and evolution of the MRGX2 gene in world-wide human populations and the representative nonhuman primate species. Our results demonstrated that MRGX2 had undergone adaptive changes in the path of human evolution, which were likely caused by Darwinian positive selection. The patterns of DNA sequence polymorphisms in human populations showed an excess of derived substitutions, which against the expectation of neutral evolution, implying that the adaptive evolution of MRGX2 in humans was a relatively recent event. The reconstructed secondary structure of the human MRGX2 revealed that three of the four human-specific amino acid substitutions were located in the extra-cellular domains. Such critical substitutions may alter the interactions between MRGX2 protein and its ligand, thus, potentially led to adaptive changes of the pain-perception-related nervous system during human evolution.