GPCR-I-TASSER: A Hybrid Approach to G Protein-Coupled Receptor Structure Modeling and the Application to the Human Genome

Genomic and phenotypic evidence support visual and olfactory shifts in primate evolution

Sensory trade-offs between vision and olfaction in the evolution and radiation of primates have long been debated. However, insights have been limited by a lack of sensory gene sequences and accompanying functional predictions. Here we conduct large-scale functional analyses of visual and olfactory receptors and related brain regions across extant primates. Our results reveal a visual shift from ultraviolet to violet colour sensitivity in early haplorrhine primates, followed by acceleration in the rhodopsin retinal release rates at the origin of anthropoids, both of which are expected to greatly enhance visual acuity under brighter light conditions. Additionally, we find that the sensitivity of olfactory receptors shifted from narrowly to broadly tuned early in anthropoid evolution. In contrast, strepsirrhines appear to have retained sensitive dim-light vision and underwent functional enhancement of narrowly tuned olfactory receptors. Our models indicate that this would have enhanced odorant discrimination and facilitated olfaction-mediated physiology and behaviour. These differences in tuning patterns of olfactory receptors between major primate lineages mirror well-established morphological differences in external anatomy and brain structures, revealing new mechanisms of olfactory adaptation and evolutionary plasticity. Our multisystem analyses reveal patterns of co-evolution in genomic, molecular and neuroanatomical traits that are consistent with a sensory 'reallocation' rather than strict trade-offs.

Recent Applications of In Silico Approaches for Studying Receptor Mutations Associated with Human Pathologies

In recent years, the advent of computational techniques to predict the potential activity of a drug interacting with a receptor or to predict the structure of unidentified proteins with aberrant characteristics has significantly impacted the field of drug design. We provide a comprehensive review of the current state of in silico approaches and software for investigating the effects of receptor mutations associated with human diseases, focusing on both frequent and rare mutations. The reported techniques include virtual screening, homology modeling, threading, docking, and molecular dynamics. This review clearly shows that it is common for successful studies to integrate different techniques in drug design, with docking and molecular dynamics being the most frequently used techniques. This trend reflects the current emphasis on developing novel therapies for diseases resulting from receptor mutations with the recently discovered AlphaFold algorithm as the driving force.

Resistance and virulence genes characteristic of a South Asia Clade (I) Candida auris strain isolated from blood in Beijing

INTRODUCTION

Candida auris is a globally disseminated invasive ascomycetous yeast, that imposes a substantial burden on healthcare systems. It has been documented to have spread to over 40 countries across six continents, necessitating in-depth comprehension through advanced techniques like Whole-Genome Sequencing.

METHOD

This study entailed the isolation and Whole-Genome Sequencing of a fluconazole-resistant C. auris strain (CA01) obtained from a patient's blood in Beijing. Genome analysis was conducted to classify the strain, and molecular docking was performed to understand the impact of mutations on drug resistance.

RESULTS

Genome analysis revealed that CA01 belongs to the South Asia Clade (I) and shares the closest genetic relationship with previously reported strains BJCA001 and BJCA002. Notably, unlike BJCA001, CA01 exhibits significant resistance to fluconazole primarily due to the A395T mutation in the ERG11 gene. Molecular docking studies demonstrated that this mutation leads to geometric changes in the active site where fluconazole binds, resulting in decreased binding affinity. Additionally, the present findings have identified several core virulence genes in C. auris, such as RBF1.

DISCUSSION

The findings from this study expand the understanding of the genetic diversity and adaptive mechanisms of C. auris within the South Asia Clade (I). The observed fluconazole resistance driven by the ERG11 mutation A395T highlights the need for heightened awareness and adaptation in clinical treatment strategies in China. This study provides critical insights into drug resistance and virulence profiles at a genetic level, which could guide future therapeutic and management strategies for C. auris infections.

Comprehensive understanding of the role of GPER in estrogen receptor-alpha negative breast cancer

G protein-coupled estrogen receptor (GPER) plays a prominent role in facilitating the rapid, non-genomic signaling of estrogens in breast cancer cells. Herein, a comprehensive overview of the role of GPER in ER-ɑ-negative breast cancer is provided. Activation of GPER affected proliferation, metastasis and epithelial mesenchymal transition in ER-ɑ negative breast cancer cells. Clinical studies have demonstrated that GPER positivity was strongly correlated with larger tumor size and advanced clinical stage, suggesting that GPER/ERK signaling may play a role in promoting tumor progression. Strong evidence existed that environmental contaminants like bisphenol A have a carcinogenic potential mediated by GPER activation. The complexity of the cross talk between GPER and other receptors including ER-β, ER-α36, Estrogen-related receptor α (ERRα) and androgen receptor has been discussed. The potential utility of small molecules and phytoestrogens targeting GPER, adds valuable insights into its therapeutic potential. This review holds promises in advancing our understanding of GPER role in ER-ɑ-negative breast cancer. Overall, the consequences of GPER activation are still an area of active research and the implication are not entirely clear.

Evolution of olfactory receptor superfamily in bats based on high throughput molecular modelling

The origin of flight and laryngeal echolocation in bats is likely to have been accompanied by evolutionary changes in other aspects of their sensory biology. Of all sensory modalities in bats, olfaction is perhaps the least well understood. Olfactory receptors (ORs) function in recognizing odour molecules, with crucial roles in evaluating food, as well as in processing social information. Here we compare OR repertoire sizes across taxa and apply a new pipeline that integrates comparative genome data with protein structure modelling and then we employ molecular docking techniques with small molecules to analyse OR functionality based on binding energies. Our results suggest a sharp contraction in odorant recognition of the functional OR repertoire during the origin of bats, consistent with a reduced dependence on olfaction. We also compared bat lineages with contrasting different ecological characteristics and found evidence of differences in OR gene expansion and contraction, and in the composition of ORs with different tuning breadths. The strongest binding energies of ORs in non-echolocating fruit-eating bats were seen to correspond to ester odorants, although we did not detect a quantitative advantage of functional OR repertoires in these bats compared with echolocating insectivorous species. Overall, our findings based on molecular modelling and computational docking suggest that bats have undergone olfactory evolution linked to dietary adaptation. Our results from extant and ancestral bats help to lay the groundwork for targeted experimental functional tests in the future.

Next-generation neuropeptide Y receptor small-molecule agonists inhibit mosquito-biting behavior

BACKGROUND

Female Aedes aegypti mosquitoes can spread disease-causing pathogens when they bite humans to obtain blood nutrients required for egg production. Following a complete blood meal, host-seeking is suppressed until eggs are laid. Neuropeptide Y-like receptor 7 (NPYLR7) plays a role in endogenous host-seeking suppression and previous work identified small-molecule NPYLR7 agonists that inhibit host-seeking and blood-feeding when fed to mosquitoes at high micromolar doses.

METHODS

Using structure-activity relationship analysis and structure-guided design we synthesized 128 compounds with similarity to known NPYLR7 agonists.

RESULTS

Although in vitro potency (EC50) was not strictly predictive of in vivo effect, we identified three compounds that reduced blood-feeding from a live host when fed to mosquitoes at a dose of 1 μM-a 100-fold improvement over the original reference compound.

CONCLUSIONS

Exogenous activation of NPYLR7 represents an innovative vector control strategy to block mosquito biting behavior and prevent mosquito-human host interactions that lead to pathogen transmission.

The G Protein Estrogen Receptor (GPER) is involved in the resistance to the CDK4/6 inhibitor palbociclib in breast cancer

BACKGROUND

The cyclin D1-cyclin dependent kinases (CDK)4/6 inhibitor palbociclib in combination with endocrine therapy shows remarkable efficacy in the management of estrogen receptor (ER)-positive and HER2-negative advanced breast cancer (BC). Nevertheless, resistance to palbociclib frequently arises, highlighting the need to identify new targets toward more comprehensive therapeutic strategies in BC patients.

METHODS

BC cell lines resistant to palbociclib were generated and used as a model system. Gene silencing techniques and overexpression experiments, real-time PCR, immunoblotting and chromatin immunoprecipitation studies as well as cell viability, colony and 3D spheroid formation assays served to evaluate the involvement of the G protein-coupled estrogen receptor (GPER) in the resistance to palbociclib in BC cells. Molecular docking simulations were also performed to investigate the potential interaction of palbociclib with GPER. Furthermore, BC cells co-cultured with cancer-associated fibroblasts (CAFs) isolated from mammary carcinoma, were used to investigate whether GPER signaling may contribute to functional cell interactions within the tumor microenvironment toward palbociclib resistance. Finally, by bioinformatics analyses and k-means clustering on clinical and expression data of large cohorts of BC patients, the clinical significance of novel mediators of palbociclib resistance was explored.

RESULTS

Dissecting the molecular events that characterize ER-positive BC cells resistant to palbociclib, the down-regulation of ERα along with the up-regulation of GPER were found. To evaluate the molecular events involved in the up-regulation of GPER, we determined that the epidermal growth factor receptor (EGFR) interacts with the promoter region of GPER and stimulates its expression toward BC cells resistance to palbociclib treatment. Adding further cues to these data, we ascertained that palbociclib does induce pro-inflammatory transcriptional events via GPER signaling in CAFs. Of note, by performing co-culture assays we demonstrated that GPER contributes to the reduced sensitivity to palbociclib also facilitating the functional interaction between BC cells and main components of the tumor microenvironment named CAFs.

CONCLUSIONS

Overall, our results provide novel insights on the molecular events through which GPER may contribute to palbociclib resistance in BC cells. Additional investigations are warranted in order to assess whether targeting the GPER-mediated interactions between BC cells and CAFs may be useful in more comprehensive therapeutic approaches of BC resistant to palbociclib.

Estetrol/GPER/SERPINB2 transduction signaling inhibits the motility of triple-negative breast cancer cells

BACKGROUND

Estetrol (E4) is a natural estrogen produced by the fetal liver during pregnancy. Due to its favorable safety profile, E4 was recently approved as estrogenic component of a new combined oral contraceptive. E4 is a selective ligand of estrogen receptor (ER)α and ERβ, but its binding to the G Protein-Coupled Estrogen Receptor (GPER) has not been described to date. Therefore, we aimed to explore E4 action in GPER-positive Triple-Negative Breast Cancer (TNBC) cells.

METHODS

The potential interaction between E4 and GPER was investigated by molecular modeling and binding assays. The whole transcriptomic modulation triggered by E4 in TNBC cells via GPER was explored through high-throughput RNA sequencing analyses. Gene and protein expression evaluations as well as migration and invasion assays allowed us to explore the involvement of the GPER-mediated induction of the plasminogen activator inhibitor type 2 (SERPINB2) in the biological responses triggered by E4 in TNBC cells. Furthermore, bioinformatics analysis was aimed at recognizing the biological significance of SERPINB2 in ER-negative breast cancer patients.

RESULTS

After the molecular characterization of the E4 binding capacity to GPER, RNA-seq analysis revealed that the plasminogen activator inhibitor type 2 (SERPINB2) is one of the most up-regulated genes by E4 in a GPER-dependent manner. Worthy, we demonstrated that the GPER-mediated increase of SERPINB2 is engaged in the anti-migratory and anti-invasive effects elicited by E4 in TNBC cells. In accordance with these findings, a correlation between SERPINB2 levels and a good clinical outcome was found in ER-negative breast cancer patients.

CONCLUSIONS

Overall, our results provide new insights into the mechanisms through which E4 can halt migratory and invasive features of TNBC cells.

Next Generation Neuropeptide Y Receptor Small Molecule Agonists Inhibit Mosquito Biting Behavior

Female Aedes aegypti mosquitoes can spread disease-causing pathogens when they bite humans to obtain blood nutrients required for egg production. Following a complete blood meal, host-seeking is suppressed until eggs are laid. Neuropeptide Y-like Receptor 7 (NPYLR7) plays a role in endogenous host-seeking suppression and previous work identified small molecule NPYLR7 agonists that suppress host-seeking and blood feeding when fed to mosquitoes at high micromolar doses. Using structure activity relationship analysis and structure-guided design we synthesized 128 compounds with similarity to known NPYLR7 agonists. Although in vitro potency (EC50) was not strictly predictive of in vivo effect, we identified 3 compounds that suppressed blood feeding from a live host when fed to mosquitoes at a 1 μM dose, a 100-fold improvement over the original reference compound. Exogenous activation of NPYLR7 represents an innovative vector control strategy to block mosquito biting behavior and prevent mosquito/human host interactions that lead to pathogen transmission.

HIV-1 Vif protein sequence variations in South African people living with HIV and their influence on Vif-APOBEC3G interaction

PURPOSE

Despite extensive research, HIV-1 remains a global epidemic with variations in pathogenesis across regions and subtypes. The Viral Infectivity Factor (Vif) protein, which neutralizes the host protein APOBEC3G, has been implicated in differences in clinical outcomes among people living with HIV (PLHIV). Most studies on Vif sequence diversity have focused on subtype B, leaving gaps in understanding Vif variations in HIV-1C regions like South Africa. This study aimed to identify and compare Vif sequence diversity in a cohort of 51 South African PLHIV and other HIV-1C prevalent regions.

METHODS

Sanger sequencing was used for Vif analysis in the cohort, and additional sequences were obtained from the Los Alamos database. Molecular modeling and docking techniques were employed to study the influence of subtype-specific variants on Vif-APOBEC3G binding affinity.

RESULTS

The findings showed distinct genetic variations between Vif sequences from India and Uganda, while South African sequences had wider distribution and closer relatedness to both. Specific amino acid substitutions in Vif were associated with geographic groups. Molecular modeling and docking analyses consistently identified specific residues (ARGR19, LYS26, TYR30, TYR44, and TRP79) as primary contributors to intermolecular contacts between Vif and APOBEC3G, essential for their interaction. The Indian Vif variant exhibited the highest predicted binding affinity to APOBEC3G among the studied groups.

CONCLUSIONS

These results provide insights into Vif sequence diversity in HIV-1C prevalent regions and shed light on differential pathogenesis observed in different geographical areas. The identified Vif amino acid residues warrant further investigation for their diagnostic, prognostic, and therapeutic potential.

Gpr75 knockout mice display age-dependent cone photoreceptor cell loss

GPR75 is an orphan G protein-coupled receptor for which there is currently limited information and its function in physiology and disease is only recently beginning to emerge. This orphan receptor is expressed in the retina but its function in the eye is unknown. The earliest studies on GPR75 were conducted in the retina, where the receptor was first identified and cloned and mutations in the receptor were identified as a possible contributor to retinal degenerative disease. Despite these sporadic reports, the function of GPR75 in the retina and in retinal disease has yet to be explored. To assess whether GPR75 has a functional role in the retina, the retina of Gpr75 knockout mice was characterized. Knockout mice displayed a mild progressive retinal degeneration, which was accompanied by oxidative stress. The degeneration was because of the loss of both M-cone and S-cone photoreceptor cells. Housing mice under constant dark conditions reduced oxidative stress but did not prevent cone photoreceptor cell loss, indicating that oxidative stress is not a primary cause of the observed retinal degeneration. Studies here demonstrate an important role for GPR75 in maintaining the health of cone photoreceptor cells and that Gpr75 knockout mice can be used as a model to study cone photoreceptor cell loss.

Molecular dynamics and free energy calculations of clozapine bound to D2 and H1 receptors reveal a cardiometabolic mitigated derivative

Most atypical antipsychotics derive from a high dropout of drug treatments due to adverse cardiometabolic side effects. These side effects are caused, in part, by the H1 receptor blockade. The current work sought a clozapine derivative with a reduced affinity for the H1 receptor while maintaining its therapeutic effect linked to D2 receptor binding. Explicit molecular dynamics simulations and end-point free energy calculations of clozapine in complex with the D2 and H1 receptors embedded in cholesterol-rich lipid bilayers were performed to analyze the intermolecular interactions and address the relevance of clozapine-functional groups. Based on that, free energy perturbation calculations were performed to measure the change in free energy of clozapine structural modifications. Our results indicate the best clozapine derivative is the iodine atom substitution for chlorine. The latter is mainly due to electrostatic interaction loss for the H1 receptor, while the halogen orientation out of the D2 active site reduces the impact on the affinity.Communicated by Ramaswamy H. Sarma.

Development of anti-PD-L1 antibody based on structure prediction of AlphaFold2

Accurate structural information plays a crucial role in comprehending biological processes and designing drugs. Indeed, the remarkable precision of the AlphaFold2 has facilitated significant advancements in predicting molecular structures, encompassing antibodies and antigens. This breakthrough has paved the way for rational drug design, ushering in new possibilities in the field of pharmaceutical development. Within this study, performing analysis and humanization guided by the structures predicted by AlphaFold2. Notably, the resulting humanized antibody, h3D5-hIgG1, demonstrated exceptional binding affinity to the PD-L1 protein. The KD value of parental antibody 3D5-hIgG1 was increased by nearly 7 times after humanization. Both h3D5-hIgG1 and 3D5-hIgG1 bound to cells expressing human PD-L1 with EC50 values of 5.13 and 9.92nM, respectively. Humanization resulted in a twofold increase in the binding capacity of the antibody, with h3D5-hIgG1 exhibiting superior performance compared to the parental antibody 3D5-hIgG1. Furthermore, h3D5-hIgG1 promoted cytokine secretion of T cells, and significantly suppressed MC38-hPD-L1 tumor growth. This study highlights the potential for artificial intelligence-assisted drug development, which is poised to become a prominent trend in the future.

Developing human olfactory network and exploring olfactory receptor-odorant interaction

The Olfactory receptor (OR)-odorant interactions are perplexed. ORs can bind to structurally diverse odorants associated with one or more odor percepts. Various attempts have been made to understand the intricacies of OR-odorant interaction. In this study, experimentally documented OR-odorant interactions are investigated comprehensively to; (a) suggest potential odor percepts for ORs based on the OR-OR network; (b) determine how odorants interacting with specific ORs differ in terms of inherent pharmacophoric features and molecular properties, (c) identify molecular interactions that explained OR-odorant interactions of selective ORs; and (d) predict the probable role of ORs other than olfaction. Human olfactory receptor network (hORnet) is developed to study possible odor percepts for ORs. We identified six molecular properties which showed variation and significant patterns to differentiate odorants binding with five ORs. The pharmacophore analysis revealed that odorants subset of five ORs follow similar pharmacophore hypothesis, (one hydrogen acceptor and two hydrophobic regions) but differ in terms of distance and orientation of pharmacophoric features. To ascertain the binding site residues and key interactions between the selected ORs and their interacting odorants, 3D-structure modelling, docking and molecular dynamics studies were carried out. Lastly, the potential role of ORs beyond olfaction is explored. A human OR-OR network was developed to suggest possible odor percepts for ORs using empirically proven OR-odorant interactions. We sought to find out significant characteristics, molecular properties, and molecular interactions that could explain OR-odorant interactions and add to the understanding of the complex issue of odor perception.Communicated by Ramaswamy H. Sarma.

Adaptations to a cold climate promoted social evolution in Asian colobine primates

The biological mechanisms that underpin primate social evolution remain poorly understood. Asian colobines display a range of social organizations, which makes them good models for investigating social evolution. By integrating ecological, geological, fossil, behavioral, and genomic analyses, we found that colobine primates that inhabit colder environments tend to live in larger, more complex groups. Specifically, glacial periods during the past 6 million years promoted the selection of genes involved in cold-related energy metabolism and neurohormonal regulation. More-efficient dopamine and oxytocin pathways developed in odd-nosed monkeys, which may have favored the prolongation of maternal care and lactation, increasing infant survival in cold environments. These adaptive changes appear to have strengthened interindividual affiliation, increased male-male tolerance, and facilitated the stepwise aggregation from independent one-male groups to large multilevel societies.

Helix encoder: a compound-protein interaction prediction model specifically designed for class A GPCRs

Class A G protein-coupled receptors (GPCRs) represent the largest class of GPCRs. They are essential targets of drug discovery and thus various computational approaches have been applied to predict their ligands. However, there are a large number of orphan receptors in class A GPCRs and it is difficult to use a general protein-specific supervised prediction scheme. Therefore, the compound-protein interaction (CPI) prediction approach has been considered one of the most suitable for class A GPCRs. However, the accuracy of CPI prediction is still insufficient. The current CPI prediction model generally employs the whole protein sequence as the input because it is difficult to identify the important regions in general proteins. In contrast, it is well-known that only a few transmembrane helices of class A GPCRs play a critical role in ligand binding. Therefore, using such domain knowledge, the CPI prediction performance could be improved by developing an encoding method that is specifically designed for this family. In this study, we developed a protein sequence encoder called the Helix encoder, which takes only a protein sequence of transmembrane regions of class A GPCRs as input. The performance evaluation showed that the proposed model achieved a higher prediction accuracy compared to a prediction model using the entire protein sequence. Additionally, our analysis indicated that several extracellular loops are also important for the prediction as mentioned in several biological researches.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

Structural analysis of human G-protein-coupled receptor 17 ligand binding sites

The human G protein coupled membrane receptor (GPR17), the sensor of brain damage, is identified as a biomarker for many neurological diseases. In human brain tissue, GPR17 exist in two isoforms, long and short. While cryo-electron microscopy technology has provided the structure of the long isoform of GPR17 with Gi complex, the structure of the short isoform and its activation mechanism remains unclear. Recently, we theoretically modeled the structure of the short isoform of GPR17 with Gi signaling protein and identified novel ligands. In the present work, we demonstrated the presence of two distinct ligand binding sites in the short isoform of GPR17. The molecular docking of GPR17 with endogenous (UDP) and synthetic ligands (T0510.3657, MDL29950) found the presence of two distinct binding pockets. Our observations revealed that endogenous ligand UDP can bind stronger in two different binding pockets as evidenced by glide and autodock vina scores, whereas the other two ligand's binding with GPR17 has less docking score. The analysis of receptor-UDP interactions shows complexes' stability in the lipid environment by 100 ns atomic molecular dynamics simulations. The amino acid residues VAL83, ARG87, and PHE111 constitute ligand binding site 1, whereas site 2 constitutes ASN67, ARG129, and LYS232. Root mean square fluctuation analysis showed the residues 83, 87, and 232 with higher fluctuations during molecular dynamics simulation in both binding pockets. Our findings imply that the residues of GPR17's two binding sites are crucial, and their interaction with UDP reveals the protein's hidden signaling and communication properties. Furthermore, this finding may assist in the development of targeted therapies for the treatment of neurological diseases.

Specificities of Modeling of Membrane Proteins Using Multi-Template Homology Modeling

Structures of membrane proteins are challenging to determine experimentally and currently represent only about 2% of the structures in the Protein Data Bank. Because of this disparity, methods for modeling membrane proteins are fewer and of lower quality than those for modeling soluble proteins. However, better expression, crystallization, and cryo-EM techniques have prompted a recent increase in experimental structures of membrane proteins, which can act as templates to predict the structure of closely related proteins through homology modeling. Because homology modeling relies on a structural template, it is easier and more accurate than fold recognition methods or de novo modeling, which are used when the sequence similarity between the query sequence and the sequence of related proteins in structural databases is below 25%. In homology modeling, a query sequence is mapped onto the coordinates of a single template and refined. With the increase in available templates, several templates often cover overlapping segments of the query sequence. Multi-template modeling can be used to identify the best template for local segments and join them into a single model. Here we provide a protocol for modeling membrane proteins from multiple templates in the Rosetta software suite. This approach takes advantage of several integrated frameworks, namely, RosettaScripts, RosettaCM, and RosettaMP with the membrane scoring function.

Molecular modeling and dynamics studies of the synthetic small molecule agonists with GPR17 and P2Y1 receptor

The human Guanine Protein coupled membrane Receptor 17 (hGPR17), an orphan receptor that activates uracil nucleotides and cysteinyl leukotrienes is considered as a crucial target for the neurodegenerative diseases. Yet, the detailed molecular interaction of potential synthetic ligands of GPR17 needs to be characterized. Here, we have studied a comparative analysis on the interaction specificity of GPR17-ligands with hGPR17 and human purinergic G protein-coupled receptor (hP2Y1) receptors. Previously, we have simulated the interaction stability of synthetic ligands such as T0510.3657, AC1MLNKK, and MDL29951 with hGPR17 and hP2Y1 receptor in the lipid environment. In the present work, we have comparatively studied the protein-ligand interaction of hGPR17-T0510.3657 and P2Y1-MRS2500. Sequence analysis and structural superimposition of hGPR17 and hP2Y1 receptor revealed the similarities in the structural arrangement with the local backbone root mean square deviation (RMSD) value of 1.16 Å and global backbone RMSD value of 5.30 Å. The comparative receptor-ligand interaction analysis between hGPR17 and hP2Y1 receptor exposed the distinct binding sites in terms of geometrical properties. Further, the molecular docking of T0510.3657 with the hP2Y1 receptor have shown non-specific interaction. The experimental validation also revealed that Gi-coupled activation of GPR17 by specific ligands leads to the adenylyl cyclase inhibition, while there is no inhibition upon hP2Y1 activation. Overall, the above findings suggest that T0510.3657-GPR17 binding specificity could be further explored for the treatment of numerous neuronal diseases. Communicated by Ramaswamy H. Sarma.

Multi-omics approach to identifying isoform variants as therapeutic targets in cancer patients

Cancer-specific alternatively spliced events (ASE) play a role in cancer pathogenesis and can be targeted by immunotherapy, oligonucleotide therapy, and small molecule inhibition. However, identifying actionable ASE targets remains challenging due to the uncertainty of its protein product, structure impact, and proteoform (protein isoform) function. Here we argue that an integrated multi-omics profiling strategy can overcome these challenges, allowing us to mine this untapped source of targets for therapeutic development. In this review, we will provide an overview of current multi-omics strategies in characterizing ASEs by utilizing the transcriptome, proteome, and state-of-art algorithms for protein structure prediction. We will discuss limitations and knowledge gaps associated with each technology and informatics analytics. Finally, we will discuss future directions that will enable the full integration of multi-omics data for ASE target discovery.

Multi-state modeling of G-protein coupled receptors at experimental accuracy

The family of G-protein coupled receptors (GPCRs) is one of the largest protein families in the human genome. GPCRs transduct chemical signals from extracellular to intracellular regions via a conformational switch between active and inactive states upon ligand binding. While experimental structures of GPCRs remain limited, high-accuracy computational predictions are now possible with AlphaFold2. However, AlphaFold2 only predicts one state and is biased toward either the active or inactive conformation depending on the GPCR class. Here, a multi-state prediction protocol is introduced that extends AlphaFold2 to predict either active or inactive states at very high accuracy using state-annotated templated GPCR databases. The predicted models accurately capture the main structural changes upon activation of the GPCR at the atomic level. For most of the benchmarked GPCRs (10 out of 15), models in the active and inactive states were closer to their corresponding activation state structures. Median RMSDs of the transmembrane regions were 1.12 Å and 1.41 Å for the active and inactive state models, respectively. The models were more suitable for protein-ligand docking than the original AlphaFold2 models and template-based models. Finally, our prediction protocol predicted accurate GPCR structures and GPCR-peptide complex structures in GPCR Dock 2021, a blind GPCR-ligand complex modeling competition. We expect that high accuracy GPCR models in both activation states will promote understanding in GPCR activation mechanisms and drug discovery for GPCRs. At the time, the new protocol paves the way towards capturing the dynamics of proteins at high-accuracy via machine-learning methods.

Effects of endocrine active contaminating pesticides on RACK1 expression and immunological consequences in THP-1 cells

We have previously demonstrated that RACK1, which expression is under steroid hormone control, plays an important role in the activation of immune cells and its expression can be useful to evaluate the immunotoxic profile of endocrine disrupting chemicals (EDCs). Hence, we investigated the effects of three contaminating and persistent pesticides: the fungicide vinclozolin (VIN), the herbicide atrazine (ATR) and the insecticide cypermethrin (CYP) on RACK1 expression and on innate immune response. VIN resulted in modest alteration of RACK1 while ATR and CYP reduced in a dose dependent manner RACK1 expression, ultimately leading to the decrease in lipopolysaccharide-induced IL-8 and TNF-α release and CD86 and CD54 surface marker expression. Moreover, our data indicate that, after exposure to EDCs, alterations of RACK1 expression can also occur with mechanisms not directly mediated by an interaction with a nuclear or membrane steroid receptors. Therefore, RACK1 could represent a useful EDCs screening tool to evaluate their immunotoxic potential and to dissect their mechanisms of action.

Structural dynamics of Smoothened (SMO) in the ciliary membrane and its interaction with membrane lipids

The Smoothened receptor (SMO, a 7 pass transmembrane domain, Class F GPCR family protein) plays a crucial role in the Hedgehog (HH) signaling pathway, which is involved in embryonic development and is implicated in various types of cancer throughout the animal kingdom. In the absence of HH signaling, SMO is inhibited by Patched 1 (PTC1; a 12 pass transmembrane domain protein), which is localized in the primary cilia. HH binding leads to the dislocation of PTC1 from the cilia, thus making way for SMO to localize in the primary cilia, as an essential prerequisite for its activation. We have carried out MARTINI coarse-grained molecular dynamics simulations of SMO in POPC and in ciliary membrane models, respectively, to study the interactions of SMO with cholesterol and other lipid molecules in the ciliary membrane, and to gain molecular-level insights into the role of the primary cilia in shaping the functional dynamics of SMO. We are able to identify the interaction of membrane cholesterols with definite sites and domains within SMO and relate them with known cholesterol-binding sequence and structure motifs. We show that cholesterol interactions with the transmembrane domain TMD, unlike those with the cysteine-rich domain (CRD) and the intracellular domain (ICD), are through residues belonging to known cholesterol-binding motifs. Notably, a few persistent interactions of cholesterol with lower TM cholesterol-binding domains are governed by the presence of multiple cholesterol-binding motifs. These analyses have also helped to identify and define a strict cholesterol consensus motif (CCM), which may well steer cholesterol into the hitherto identified binding sites within the TMD of SMO. We have also reported the interaction of phosphatidylinositol 4-phosphate with the intracellular region of transmembrane (TM) helices (TM1, TM3, TM4, and TM5), intracellular loop1, helix8, and Arg/Lys clusters of the ICD. Structural analysis of SMO domains shows significant changes in the CRD and ICD, during the course of the simulation. Further detailed analysis of the dynamics of the TMD reveals the movements of TM5, TM6, and TM7, linked with the helix8, which are possibly involved in shaping the conformational disposition of the ICD. The movement of these TM helices could possibly be a consequence of interactions involving the extracellular domain and extracellular loops. In addition, our analysis also shows that phosphatidylinositol-4-phosphate (PI4P), along with some ICD cholesterols, are implicated in anchoring SMO in the membrane.

Computational insights on molecular interactions of acifran with GPR109A and GPR109B

Acifran is a well-known agonist of G-protein-coupled receptor protein, namely GPR109A. Acifran is primarily used in the treatment of dyslipidemia, myocardial infractions, and atherosclerosis in humans due to its lower vascular and metabolic side effects. However, experimental and computational studies on interaction sites of acifran with GPR proteins (GPR109A and GPR109B) are lacking. Our computational studies using docking and molecular dynamics simulation revealed that acifran binds distinctly to both GPR109A and GPR109B, but with lower affinity to the latter. The weak binding of acifran-GPR109B is mainly due to the presence of residues S91 and N94 in ECL1 and I178 amino acid in ECL2 region of GPR109B, whereas R111 and R251 residues in TMH3 and TMH6 are crucial for GPR109A-acifran complex stability. Additionally, molecular mechanics/Poisson-Boltzmann solvent accessible surface area (MM/PBSA) analysis revealed that both GPR109A- and GPR109B-acifran complexes are energetically stable with lower calculated binding free energy values for the latter. Energy-minimized structures of GPR109A-acifran and GPR109B-acifran complex.

Computational and experimental approaches to probe GPCR activation and signaling

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Hyaluronic acid-GPRC5C signalling promotes dormancy in haematopoietic stem cells

Bone marrow haematopoietic stem cells (HSCs) are vital for lifelong maintenance of healthy haematopoiesis. In inbred mice housed in gnotobiotic facilities, the top of the haematopoietic hierarchy is occupied by dormant HSCs, which reversibly exit quiescence during stress. Whether HSC dormancy exists in humans remains debatable. Here, using single-cell RNA sequencing, we show a continuous landscape of highly purified human bone marrow HSCs displaying varying degrees of dormancy. We identify the orphan receptor GPRC5C, which enriches for dormant human HSCs. GPRC5C is also essential for HSC function, as demonstrated by genetic loss- and gain-of-function analyses. Through structural modelling and biochemical assays, we show that hyaluronic acid, a bone marrow extracellular matrix component, preserves dormancy through GPRC5C. We identify the hyaluronic acid-GPRC5C signalling axis controlling the state of dormancy in mouse and human HSCs.

Additive and epistatic effects influence spectral tuning in molluscan retinochrome opsin

The relationship between genotype and phenotype is nontrivial due to often complex molecular pathways that make it difficult to unambiguously relate phenotypes to specific genotypes. Photopigments, an opsin apoprotein bound to a light-absorbing chromophore, present an opportunity to directly relate the amino acid sequence to an absorbance peak phenotype (λmax). We examined this relationship by conducting a series of site-directed mutagenesis experiments of retinochrome, a non-visual opsin, from two closely related species: the common bay scallop, Argopecten irradians, and the king scallop, Pecten maximus. Using protein folding models, we identified three amino acid sites of likely functional importance and expressed mutated retinochrome proteins in vitro. Our results show that the mutation of amino acids lining the opsin binding pocket are responsible for fine spectral tuning, or small changes in the λmax of these light sensitive proteins Mutations resulted in a blue or red shift as predicted, but with dissimilar magnitudes. Shifts ranged from a 16 nm blue shift to a 12 nm red shift from the wild-type λmax. These mutations do not show an additive effect, but rather suggests the presence of epistatic interactions. This work highlights the importance of binding pocket shape in the evolution of spectral tuning and builds on our ability to relate genotypic changes to phenotypes in an emerging model for opsin functional analysis.

Effect of Enterohepatic Circulation on the Accumulation of Per- and Polyfluoroalkyl Substances: Evidence from Experimental and Computational Studies

The pharmacokinetic characteristics of per- and polyfluoroalkyl substances (PFAS) affect their distribution and bioaccumulation in biological systems. The enterohepatic circulation leads to reabsorption of certain chemicals from bile back into blood and the liver and thus influences their elimination, yet its influence on PFAS bioaccumulation remains unclear. We explored the role of enterohepatic circulation in PFAS bioaccumulation by examining tissue distribution of various PFAS in wild fish and a rat model. Computational models were used to determine the reabsorbed fractions of PFAS by calculating binding affinities of PFAS for key transporter proteins of enterohepatic circulation. The results indicated that higher concentrations were observed in blood, the liver, and bile compared to other tissues for some PFAS in fish. Furthermore, exposure to a PFAS mixture on the rat model showed that the reabsorption phenomenon appeared during 8-12 h for most long-chain PFAS. Molecular docking calculations suggest that PFAS can bind to key transporter proteins via electrostatic and hydrophobic interactions. Further regression analysis adds support to the hypothesis that binding affinity of the apical sodium-dependent bile acid transporter is the most important variable to predict the human half-lives of PFAS. This study demonstrated the critical role of enterohepatic circulation in reabsorption, distribution, and accumulation of PFAS.

Cryptotanshinone inhibits proliferation and induces apoptosis of breast cancer MCF-7 cells via GPER mediated PI3K/AKT signaling pathway

G protein-coupled estrogen receptor (GPER) was reported to be a potential target in the breast cancer therapy. This study aimed to illuminate the function of GPER and its mediated PI3K/AKT pathway in cryptotanshinone (CPT) inducing cell apoptosis and antiproliferation effect on GPER positive breast cancer MCF-7 cells. Cell proliferation was tested by MTT assay. Apoptosis rates were tested by Annexin V-FITC/PI double staining and the cell cycle was researched by flow cytometry. Autodock vina was applied to make molecular docking between CPT or estradiol and GPER. siRNA technique and GPER specific agonist G-1 or antagonist G-15 were applied to verify the mediated function of GPER. Apoptosis and cell cycle related proteins, as well as the key proteins on PI3K/AKT signaling pathway were detected by western blot. The results indicated that CPT could exert antiproliferation effects by arresting cell cycle in G2/M phase and downregulating the expression of cyclin D, cyclin B and cyclin A. Besides, apoptosis induced by CPT was observed. CPT might be a novel GPER binding compounds. Significantly, suppression of PI3K/AKT signal transduction by CPT was further increased by G-1 and decreased by G-15. The study revealed that the effect of antiproliferation and apoptosis treating with CPT on MCF-7 cells might be through the downregulation of PI3K/AKT pathway mediated by activated GPER.

Two entry tunnels in mouse TAAR9 suggest the possibility of multi-entry tunnels in olfactory receptors

Orthosteric binding sites of olfactory receptors have been well understood for ligand-receptor interactions. However, a lack of explanation for subtle differences in ligand profile of olfactory receptors even with similar orthosteric binding sites promotes more exploration into the entry tunnels of the receptors. An important question regarding entry tunnels is the number of entry tunnels, which was previously believed to be one. Here, we used TAAR9 that recognizes important biogenic amines such as cadaverine, spermine, and spermidine as a model for entry tunnel study. We identified two entry tunnels in TAAR9 and described the residues that form the tunnels. In addition, we found two vestibular binding pockets, each located in one tunnel. We further confirmed the function of two tunnels through site-directed mutagenesis. Our study challenged the existing views regarding the number of entry tunnels in the subfamily of olfactory receptors and demonstrated the possible mechanism how the entry tunnels function in odorant recognition.

Targeting GPCRs and Their Signaling as a Therapeutic Option in Melanoma

Simple Summary

Sixteen G-protein-coupled receptors (GPCRs) have been involved in melanogenesis or melanomagenesis. Here, we review these GPCRs, their associated signaling, and therapies.

Abstract

G-protein-coupled receptors (GPCRs) serve prominent roles in melanocyte lineage physiology, with an impact at all stages of development, as well as on mature melanocyte functions. GPCR ligands are present in the skin and regulate melanocyte homeostasis, including pigmentation. The role of GPCRs in the regulation of pigmentation and, consequently, protection against external aggression, such as ultraviolet radiation, has long been established. However, evidence of new functions of GPCRs directly in melanomagenesis has been highlighted in recent years. GPCRs are coupled, through their intracellular domains, to heterotrimeric G-proteins, which induce cellular signaling through various pathways. Such signaling modulates numerous essential cellular processes that occur during melanomagenesis, including proliferation and migration. GPCR-associated signaling in melanoma can be activated by the binding of paracrine factors to their receptors or directly by activating mutations. In this review, we present melanoma-associated alterations of GPCRs and their downstream signaling and discuss the various preclinical models used to evaluate new therapeutic approaches against GPCR activity in melanoma. Recent striking advances in our understanding of the structure, function, and regulation of GPCRs will undoubtedly broaden melanoma treatment options in the future.

Evolution of brain-expressed biogenic amine receptors into olfactory trace amine-associated receptors

The family of trace amine-associated receptors (TAARs) is distantly related to G protein-coupled biogenic aminergic receptors. TAARs are found in the brain as well as in the olfactory epithelium where they detect biogenic amines. However, the functional relationship of receptors from distinct TAAR subfamilies and in different species is still uncertain. Here, we perform a thorough phylogenetic analysis of 702 TAAR-like (TARL) and TAAR sequences from 48 species. We show that a clade of Tarl genes has greatly expanded in lampreys, whereas the other Tarl clade consists of only one or two orthologs in jawed vertebrates and is lost in amniotes. We also identify two small clades of Taar genes in sharks related to the remaining Taar genes in bony vertebrates, which are divided into four major clades. We further identify ligands for 61 orphan TARLs and TAARs from sea lamprey, shark, ray-finned fishes, and mammals, as well as novel ligands for two 5-hydroxytryptamine receptor 4 orthologs, a serotonin receptor subtype closely related to TAARs. Our results reveal a pattern of functional convergence and segregation: TARLs from sea lamprey and bony vertebrate olfactory TAARs underwent independent expansions to function as chemosensory receptors, whereas TARLs from jawed vertebrates retain ancestral response profiles and may have similar functions to TAAR1 in the brain. Overall, our data provide a comprehensive understanding of the evolution and ligand recognition profiles of TAARs and TARLs.

Binding site identification of G protein-coupled receptors through a 3D Zernike polynomials-based method: application to C. elegans olfactory receptors

Studying the binding processes of G protein-coupled receptors (GPCRs) proteins is of particular interest both to better understand the molecular mechanisms that regulate the signaling between the extracellular and intracellular environment and for drug design purposes. In this study, we propose a new computational approach for the identification of the binding site for a specific ligand on a GPCR. The method is based on the Zernike polynomials and performs the ligand-GPCR association through a shape complementarity analysis of the local molecular surfaces. The method is parameter-free and it can distinguish, working on hundreds of experimentally GPCR-ligand complexes, binding pockets from randomly sampled regions on the receptor surface, obtaining an Area Under ROC curve of 0.77. Given its importance both as a model organism and in terms of applications, we thus investigated the olfactory receptors of the C. elegans, building a list of associations between 21 GPCRs belonging to its olfactory neurons and a set of possible ligands. Thus, we can not only carry out rapid and efficient screenings of drugs proposed for GPCRs, key targets in many pathologies, but also we laid the groundwork for computational mutagenesis processes, aimed at increasing or decreasing the binding affinity between ligands and receptors.

Search for Structural Basis of Interactions of Biogenic Amines with Human TAAR1 and TAAR6 Receptors

The identification and characterization of ligand-receptor binding sites are important for drug development. Trace amine-associated receptors (TAARs, members of the class A GPCR family) can interact with different biogenic amines and their metabolites, but the structural basis for their recognition by the TAARs is not well understood. In this work, we have revealed for the first time a group of conserved motifs (fingerprints) characterizing TAARs and studied the docking of aromatic (β-phenylethylamine, tyramine) and aliphatic (putrescine and cadaverine) ligands, including gamma-aminobutyric acid, with human TAAR1 and TAAR6 receptors. We have identified orthosteric binding sites for TAAR1 (Asp68, Asp102, Asp284) and TAAR6 (Asp78, Asp112, Asp202). By analyzing the binding results of 7500 structures, we determined that putrescine and cadaverine bind to TAAR1 at one site, Asp68 + Asp102, and to TAAR6 at two sites, Asp78 + Asp112 and Asp112 + Asp202. Tyramine binds to TAAR6 at the same two sites as putrescine and cadaverine and does not bind to TAAR1 at the selected Asp residues. β-Phenylethylamine and gamma-aminobutyric acid do not bind to the TAAR1 and TAAR6 receptors at the selected Asp residues. The search for ligands targeting allosteric and orthosteric sites of TAARs has excellent pharmaceutical potential.

Structural and Energetic Insights Into the Interaction of Niacin With the GPR109A Receptor

The transmembrane G-protein coupled receptor GPR109A has been previously shown to function as a receptor for niacin in mediating antilipolytic effects. Although administration of high doses of niacin has shown beneficial effects on lipid metabolism, however, it is often accompanied by disturbing side effects such as flushing, liver damage, glucose intolerance, or gastrointestinal problems. Thus, it is important to understand niacin-GPR109A interactions, which can be beneficial for the development of alternate drugs having antilipolytic effects with less or no side effects. To get into the structural insights on niacin binding to GPR109A, we have performed 100 nanoseconds long all-atom MD simulations of five niacin-GPR109A complexes (automatically docked pose 0, and randomly placed niacin in poses 1 to 4 in the receptor crevice) and analyzed using binding free energy calculations and H-bond analysis. Steered MD simulations were used to get an average force for niacin translocation between the bulk and the external crevice of the wild type and mutant (N86Y, W91 S, S178I, and triple mutant of all three residues) GPR109A receptors, as well as GPR109B (as a control that does not bind niacin). The H-bond analysis revealed that TMH3 residue R111 interacts with niacin in a total of 4 (poses 0 to 3) complexes, while residues C177, S178, and S179 contact niacin in complex pose 4, and all these complexes were energetically stable. According to steered MD simulations, all the GPR109A mutants and GPR109B required greater force than that of wild-type GPR109A to translocate in the external crevice, suggesting increased sterical obstacles. Thus, the residues N86 (at the junction of TMH2/ECL2), W91 (ECL2), R111 (TMH3), and ECL3 residues (C177, S178, S179) play an important role for optimal routing of niacin entry and to bind GPR109A.

Expanding the clinical phenotype associated with NIPAL4 mutation: Study of a Tunisian consanguineous family with erythrokeratodermia variabilis-Like Autosomal Recessive Congenital Ichthyosis

Erythrokeratodermia variabilis (EKV) is a rare disorder of cornification usually associated with dominant mutations in the GJB3 and GJB4 genes encoding connexins (Cx)31 and 30.3. Genetic heterogeneity of EKV has already been suggested. We investigated at the clinical and genetic level a consanguineous Tunisian family with 2 sisters presenting an autosomal recessive form of EKV to better characterize this disease. Mutational analysis initially screened the connexin genes and Whole-exome sequencing (WES) was performed to identify the molecular aetiology of the particular EKV phenotype in the proband. Migratory shaped erythematous areas are the initial presenting sign followed by relatively stable hyperkeratotic plaques are the two predominates characteristics in both patients. However, remarkable variability of morphological and dominating features of the disease were observed between patients. In particular, the younger sister (proband) exhibited ichthyosiform-like appearance suggesting Autosomal Recessive Congenital Ichthyosis (ARCI) condition. No causative mutations were detected in the GJB3 and GJB4 genes. WES results revealed a novel missense homozygous mutation in NIPAL4 gene (c.835C>G, p.Pro279Ala) in both patients. This variant is predicted to be likely pathogenic. In addition, in silico analysis of the mutated 3D domain structure predicted that this variant would result in NIPA4 protein destabilization and Mg²⁺ transport perturbation, pointing out the potential role of NIPAL4 gene in the development and maintenance of the barrier function of the epidermis. Taken togheter, these results expand the clinical phenotype associated with NIPAL4 mutation and reinforce our hypothesis of NIPAL4 as the main candidate gene for the EKV-like ARCI phenotype.

Locating ligand binding sites in G-protein coupled receptors using combined information from docking and sequence conservation

GPCRs (G-protein coupled receptors) are the largest family of drug targets and share a conserved structure. Binding sites are unknown for many important GPCR ligands due to the difficulties of GPCR recombinant expression, biochemistry, and crystallography. We describe our approach, ConDockSite, for predicting ligand binding sites in class A GPCRs using combined information from surface conservation and docking, starting from crystal structures or homology models. We demonstrate the effectiveness of ConDockSite on crystallized class A GPCRs such as the beta2 adrenergic and A2A adenosine receptors. We also demonstrate that ConDockSite successfully predicts ligand binding sites from high-quality homology models. Finally, we apply ConDockSite to predict the ligand binding sites on a structurally uncharacterized GPCR, GPER, the G-protein coupled estrogen receptor. Most of the sites predicted by ConDockSite match those found in other independent modeling studies. ConDockSite predicts that four ligands bind to a common location on GPER at a site deep in the receptor cleft. Incorporating sequence conservation information in ConDockSite overcomes errors introduced from physics-based scoring functions and homology modeling.

Effects of Bisphenols on RACK1 Expression and Their Immunological Implications in THP-1 Cells

Receptor for activated C kinase 1 (RACK1) has an important role in immune activation, and is regulated through a balance between glucocorticoid and androgen levels. We have previously demonstrated that RACK1 expression can serve as a marker for evaluation of immunotoxic profiles of hormone-active substances, such as endocrine-disrupting chemicals (EDCs). In this study, we investigated the effects of three bisphenols (BPA, BPAF, BPS) on RACK1 expression and on the innate immune responses in the THP-1 human promyelocytic cell line, a validated model for this investigation. BPA and BPAF reduced RACK1 promoter transcriptional activity, mRNA expression, and protein levels. However, BPS had the opposite effect. As expected, these results on RACK1 were paralleled by lipopolysaccharide (LPS)-induced interleukin-8 (IL-8) and tumor necrosis factor-α (TNFα) production. Since BPA and BPAF induced RACK1 expression in the presence of glucocorticoid receptor (GR) antagonist mifepristone, a role of G-protein-coupled estrogen receptor (GPER) has been considered due to their known estrogenic profile. Therefore, additional molecular effects of BPA and BPAF were unmasked after treatment with different inhibitors of well-known pivotal players of GPER-mediated signaling. BPA exerted its effects on RACK1 via NF-κB, as shown using the NF-κB inhibitor BAY11-7085 and NF-κB-specific luciferase reporter assay. Conversely, BPAF induced RACK1 up-regulation via androgen receptor (AR) activation, as confirmed by treatment with AR antagonist flutamide. Indeed, a biased agonism profile for BPA and BPAF for GPER was suggested based on their different binding modes revealed by our molecular docking. Altogether, our data suggest that RACK1 could represent an important target of EDCs and serves as a screening tool for their immunotoxic potential. Furthermore, RACK1 can be exploited to unmask multiple molecular interactions of hormone-active substances to better dissect out their mechanisms of action.

Naturally Occurring Genetic Variants in the Oxytocin Receptor Alter Receptor Signaling Profiles

The hormone oxytocin is commonly administered during childbirth to initiate and strengthen uterine contractions and prevent postpartum hemorrhage. However, patients have wide variation in the oxytocin dose required for a clinical response. To begin to uncover the mechanisms underlying this variability, we screened the 11 most prevalent missense genetic variants in the oxytocin receptor (OXTR) gene. We found that five variants, V45L, P108A, L206V, V281M, and E339K, significantly altered oxytocin-induced Ca²⁺ signaling or β-arrestin recruitment and proceeded to assess the effects of these variants on OXTR trafficking to the cell membrane, desensitization, and internalization. The variants P108A and L206V increased OXTR localization to the cell membrane, whereas V281M and E339K caused OXTR to be retained inside the cell. We examined how the variants altered the balance between OXTR activation and desensitization, which is critical for appropriate oxytocin dosing. The E339K variant impaired OXTR activation, internalization, and desensitization to roughly equal extents. In contrast, V281M decreased OXTR activation but had no effect on internalization and desensitization. V45L and P108A did not alter OXTR activation but did impair β-arrestin recruitment, internalization, and desensitization. Molecular dynamics simulations predicted that V45L and P108A prevent extension of the first intracellular loop of OXTR, thus inhibiting β-arrestin binding. Overall, our data suggest mechanisms by which OXTR genetic variants could alter clinical response to oxytocin.

Convergent olfactory trace amine-associated receptors detect biogenic polyamines with distinct motifs via a conserved binding site

Biogenic amines activate G-protein-coupled receptors (GPCRs) in the central nervous system in vertebrate animals. Several biogenic amines, when excreted, stimulate trace amine-associated receptors (TAARs), a group of GPCRs in the main olfactory epithelium, and elicit innate behaviors. How TAARs recognize amines with varying numbers of amino groups is largely unknown. We reasoned that a comparison between lamprey and mammalian olfactory TAARs, which are thought to have evolved independently but show convergent responses to polyamines, may reveal structural determinants of amine recognition. Here, we demonstrate that sea lamprey TAAR365 (sTAAR365) responds strongly to biogenic polyamines cadaverine, putrescine, and spermine, and shares a similar response profile as a mammalian TAAR (mTAAR9). Docking and site-directed mutagenesis analyses show that both sTAAR365 and mTAAR9 recognize the two amino groups of cadaverine with the conserved Asp^3.32 and Tyr^6.51 residues. sTAAR365, which has remarkable sensitivity for cadaverine (EC₅₀ = 4 nM), uses an extra residue, Thr^7.42, to stabilize ligand binding. These cadaverine recognition sites also interact with amines with four and three amino groups (spermine and spermidine, respectively). Glu^7.36 of sTAAR365 cooperates with Asp^3.32 and Thr^7.42 to recognize spermine, whereas mTAAR9 recognizes spermidine through an additional aromatic residue, Tyr^7.43. These results suggest a conserved mechanism whereby independently evolved TAAR receptors recognize amines with two, three, or four amino groups using the same recognition sites, at which sTAAR365 and mTAAR9 evolved distinct motifs. These motifs interact directly with the amino groups of the polyamines, a class of potent and ecologically important odorants, mediating olfactory signaling.

Nonantimicrobial Actions of Macrolides: Overview and Perspectives for Future Development

Macrolides are among the most widely prescribed broad spectrum antibacterials, particularly for respiratory infections. It is now recognized that these drugs, in particular azithromycin, also exert time-dependent immunomodulatory actions that contribute to their therapeutic benefit in both infectious and other chronic inflammatory diseases. Their increased chronic use in airway inflammation and, more recently, of azithromycin in COVID-19, however, has led to a rise in bacterial resistance. An additional crucial aspect of chronic airway inflammation, such as chronic obstructive pulmonary disease, as well as other inflammatory disorders, is the loss of epithelial barrier protection against pathogens and pollutants. In recent years, azithromycin has been shown with time to enhance the barrier properties of airway epithelial cells, an action that makes an important contribution to its therapeutic efficacy. In this article, we review the background and evidence for various immunomodulatory and time-dependent actions of macrolides on inflammatory processes and on the epithelium and highlight novel nonantibacterial macrolides that are being studied for immunomodulatory and barrier-strengthening properties to circumvent the risk of bacterial resistance that occurs with macrolide antibacterials. We also briefly review the clinical effects of macrolides in respiratory and other inflammatory diseases associated with epithelial injury and propose that the beneficial epithelial effects of nonantibacterial azithromycin derivatives in chronic inflammation, even given prophylactically, are likely to gain increasing attention in the future. SIGNIFICANCE STATEMENT: Based on its immunomodulatory properties and ability to enhance the protective role of the lung epithelium against pathogens, azithromycin has proven superior to other macrolides in treating chronic respiratory inflammation. A nonantibiotic azithromycin derivative is likely to offer prophylactic benefits against inflammation and epithelial damage of differing causes while preserving the use of macrolides as antibiotics.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Local Adaptation of Bitter Taste and Ecological Speciation in a Wild Mammal

Sensory systems are attractive evolutionary models to address how organisms adapt to local environments that can cause ecological speciation. However, tests of these evolutionary models have focused on visual, auditory, and olfactory senses. Here, we show local adaptation of bitter taste receptor genes in two neighboring populations of a wild mammal-the blind mole rat Spalax galili-that show ecological speciation in divergent soil environments. We found that basalt-type bitter receptors showed higher response intensity and sensitivity compared with chalk-type ones using both genetic and cell-based functional analyses. Such functional changes could help animals adapted to basalt soil select plants with less bitterness from diverse local foods, whereas a weaker reception to bitter taste may allow consumption of a greater range of plants for animals inhabiting chalk soil with a scarcity of food supply. Our study shows divergent selection on food resources through local adaptation of bitter receptors, and suggests that taste plays an important yet underappreciated role in speciation.

Empirical Potential Energy Function Toward ab Initio Folding G Protein-Coupled Receptors

Approximately 40-50 percent of all drugs targets are G protein-coupled receptors (GPCRs). Three-dimensional structure of GPCRs is important to probe their biophysical and biochemical functions and their pharmaceutical applications. Lacking reliable and high quality free function is one of the ugent problems of computational predicting the three-dimensional structure in this community. We proposed a GPCR-specified energy function composed of four novel empirical potential energy terms: a two-dimensional contact energy force field, knowledge-based helix pair connection distance energy term, knowledge-based helix pair angle restraint energy term and a disulfide bond energy term. To validate the energy function, we employed an ab initio GPCR three-dimensional structure predictor to test if the energy function improved the accuracy of prediction. We evaluated 28 solved GPCRs and found that 21(75 percent) targets were correctly folded (TM-score>0.5). Also, the average TM-score using the energy function was 0.54, which was improved 134 percent than the TM-score 0.23 for MODELLER energy function and 170 percent than the TM-score 0.20 for Rosetta membrane energy function. The results confirmed that our empirical potential energy function toward ab initio folding is competitive to state-of-the-art solutions for structural prediction of GPCRs.

Theoretical study of the interactions between peptide tyrosine tyrosine [PYY (1-36)], a newly identified modulator in type 2 diabetes pathophysiology, with receptors NPY1R and NPY4R

PURPOSE

Diabetes mellitus is a common condition in the clinically obese. Bariatric surgery is one of the ways to put type 2 diabetes in remission. Recent findings propose the appetite-regulator peptide tyrosine tyrosine (PYY) as a therapeutic option for patients with type 2 diabetes. This novel gut hormone restores impaired insulin and glucagon secretion in pancreatic islets and is implicated in type 2 diabetes reversal after bariatric surgery. The current study elucidates the interactions between PYY and the NPY1R and NPY4R receptors using computational methods.

METHODS

Protein structure prediction, molecular docking simulation, and molecular dynamics (MD) simulation were performed to elucidate the interactions of PYY with NPY1R and NPY4R.

RESULTS

The predicted binding models of PYY-NPY receptors are in agreement with those described in the literature, although different interaction partners are presented for the C-terminal tail of PYY. Non-polar interactions are predicted to drive the formation of the protein complex. The calculated binding energies show that PYY has higher affinity for NPY4R (ΔG_GBSA = -65.08 and ΔG_PBSA = -87.62 kcal/mol) than for NPY1R (ΔG_GBSA = -23.11 and ΔG_PBSA = -50.56 kcal/mol).

CONCLUSIONS

Based on the constructed models, the binding conformations obtained from docking and MD simulation for both the PYY-NPY1R and PYY-NPY4R complexes provide a detailed map of possible interactions. The calculated binding energies show a higher affinity of PYY for NPY4R. These findings may help to understand the mechanisms behind the improvement of diabetes following bariatric surgery.

Highly accurate protein structure prediction for the human proteome

Protein structures can provide invaluable information, both for reasoning about biological processes and for enabling interventions such as structure-based drug development or targeted mutagenesis. After decades of effort, 17% of the total residues in human protein sequences are covered by an experimentally determined structure1. Here we markedly expand the structural coverage of the proteome by applying the state-of-the-art machine learning method, AlphaFold2, at a scale that covers almost the entire human proteome (98.5% of human proteins). The resulting dataset covers 58% of residues with a confident prediction, of which a subset (36% of all residues) have very high confidence. We introduce several metrics developed by building on the AlphaFold model and use them to interpret the dataset, identifying strong multi-domain predictions as well as regions that are likely to be disordered. Finally, we provide some case studies to illustrate how high-quality predictions could be used to generate biological hypotheses. We are making our predictions freely available to the community and anticipate that routine large-scale and high-accuracy structure prediction will become an important tool that will allow new questions to be addressed from a structural perspective.