Modeling, molecular dynamics simulation, and mutation validation for structure of cannabinoid receptor 2 based on known crystal structures of GPCRs

Single nucleotide polymorphisms in the cannabinoid CB2 receptor: Molecular pharmacology and disease associations

Preclinical evidence implicating cannabinoid receptor 2 (CB2) in various diseases has led researchers to question whether CB2 genetics influence aetiology or progression. Associations between conditions and genetic loci are often studied via single nucleotide polymorphism (SNP) prevalence in case versus control populations. In the CNR2 coding exon, ~36 SNPs have high overall population prevalence (minor allele frequencies [MAF] ~37%), including non-synonymous SNP (ns-SNP) rs2501432 encoding CB2 63Q/R. Interspersed are ~27 lower frequency SNPs, four being ns-SNPs. CNR2 introns also harbour numerous SNPs. This review summarises CB2 ns-SNP molecular pharmacology and evaluates evidence from ~70 studies investigating CB2 genetic variants with proposed linkage to disease. Although CNR2 genetic variation has been associated with a wide variety of conditions, including osteoporosis, immune-related disorders, and mental illnesses, further work is required to robustly validate CNR2 disease links and clarify specific mechanisms linking CNR2 genetic variation to disease pathophysiology and potential drug responses.

Cannabinoid CB2 receptor orthologues; in vitro function and perspectives for preclinical to clinical translation

Cannabinoid CB2 receptor agonists are in development as therapeutic agents, including for immune modulation and pain relief. Despite promising results in rodent preclinical studies, efficacy in human clinical trials has been marginal to date. Fundamental differences in ligand engagement and signalling responses between the human CB2 receptor and preclinical model species orthologues may contribute to mismatches in functional outcomes. This is a tangible possibility for the CB2 receptor in that there is a relatively large degree of primary amino acid sequence divergence between human and rodent. Here, we summarise CB2 receptor gene and protein structure, assess comparative molecular pharmacology between CB2 receptor orthologues, and review the current status of preclinical to clinical translation for drugs targeted at the CB2 receptor, focusing on comparisons between human, mouse and rat receptors. We hope that raising wider awareness of, and proposing strategies to address, this additional challenge in drug development will assist in ongoing efforts toward successful therapeutic translation of drugs targeted at the CB2 receptor. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Metabolism-guided development of Ko143 analogs as ABCG2 inhibitors

ATP-binding cassette subfamily G member 2 (ABCG2), an efflux transporter, is involved in multiple pathological processes. Ko143 is a potent ABCG2 inhibitor; however, it is quickly metabolized through carboxylesterase 1-mediated hydrolysis of its t-butyl ester moiety. The current work aimed to develop more metabolically stable ABCG2 inhibitors. Novel Ko143 analogs were designed and synthesized by replacing the unstable t-butyl ester moiety in Ko143 with an amide group. The synthesized Ko143 analogs were evaluated for their ABCG2 inhibitory activity, binding mode with ABCG2, cytotoxicity, and metabolic stability. We found that the amide modification of Ko143 led to metabolically stable ABCG2 inhibitors. Among these Ko143 analogs, K2 and K34 are promising candidates with favorable oral pharmacokinetic profiles in mice. In summary, we synthesized novel Ko143 analogs with improved metabolic stability, which can potentially be used as lead compounds for the future development of ABCG2 inhibitors.

Anandamide reduces the migration of lymphocytes to the intestine by CB2 activation and reduces TNF-α in the target organs, protecting mice from graft-versus-host disease

Graft-versus-host disease (GVHD) is a serious inflammatory illness that often occurs as a secondary complication of bone marrow transplantation. Current therapies have limited effectiveness and fail to achieve a balance between inflammation and the graft-versus-tumor effect. In this study, we investigate the effects of the endocannabinoid anandamide on the complex pathology of GVHD. We assess the effects of an irreversible inhibitor of fatty acid amine hydrolase or exogenous anandamide and find that they increase survival and reduce clinical signs in GVHD mice. In the intestine of GVHD mice, treatment with exogenous anandamide also leads to a reduction in the number of CD3+, CD3+CD4+, and CD3+CD8+ cells, which reduces the activation of CD3+CD4+ and CD3+CD8+ cells, as assessed by enhanced CD28 expression, a T cell co-stimulatory molecule. Exogenous AEA was also able to reduce TNF-α and increase IL-10 in the intestine of GVHD mice. In the liver, exogenous AEA reduces injury, TNF-α levels, and the number of CD3+CD8+ cells. Interestingly, anandamide reduces Mac-1α, which lowers the adhesion of transplanted cells in mesenteric veins. These effects are mimicked by JWH133-a CB2 selective agonist-and abolished by treatment with a CB2 antagonist. Furthermore, the effects caused by anandamide treatment on survival were related to the CB2 receptor, as the CB2 antagonist abolished it. This study shows the critical role of the CB2 receptor in the modulation of the inflammatory response of GVHD by treatment with anandamide, the most prominent endocannabinoid.

Goods and Bads of the Endocannabinoid System as a Therapeutic Target: Lessons Learned after 30 Years

The cannabis derivative marijuana is the most widely used recreational drug in the Western world and is consumed by an estimated 83 million individuals (∼3% of the world population). In recent years, there has been a marked transformation in society regarding the risk perception of cannabis, driven by its legalization and medical use in many states in the United States and worldwide. Compelling research evidence and the Food and Drug Administration cannabis-derived cannabidiol approval for severe childhood epilepsy have confirmed the large therapeutic potential of cannabidiol itself, Δ9-tetrahydrocannabinol and other plant-derived cannabinoids (phytocannabinoids). Of note, our body has a complex endocannabinoid system (ECS)-made of receptors, metabolic enzymes, and transporters-that is also regulated by phytocannabinoids. The first endocannabinoid to be discovered 30 years ago was anandamide (N-arachidonoyl-ethanolamine); since then, distinct elements of the ECS have been the target of drug design programs aimed at curing (or at least slowing down) a number of human diseases, both in the central nervous system and at the periphery. Here a critical review of our knowledge of the goods and bads of the ECS as a therapeutic target is presented to define the benefits of ECS-active phytocannabinoids and ECS-oriented synthetic drugs for human health. SIGNIFICANCE STATEMENT: The endocannabinoid system plays important roles virtually everywhere in our body and is either involved in mediating key processes of central and peripheral diseases or represents a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of the components of this complex system, and in particular of key receptors (like cannabinoid receptors 1 and 2) and metabolic enzymes (like fatty acid amide hydrolase and monoacylglycerol lipase), will advance our understanding of endocannabinoid signaling and activity at molecular, cellular, and system levels, providing new opportunities to treat patients.

Understanding the Dynamics of the Structural States of Cannabinoid Receptors and the Role of Different Modulators

The cannabinoid receptors CB₁R and CB₂R are members of the G protein-coupled receptor (GPCR) family. These receptors have recently come to light as possible therapeutic targets for conditions affecting the central nervous system. However, because CB₁R is known to have psychoactive side effects, its potential as a drug target is constrained. Therefore, targeting CB₂R has become the primary focus of recent research. Using various molecular modeling studies, we analyzed the active, inactive, and intermediate states of both CBRs in this study. We conducted in-depth research on the binding properties of various groups of cannabinoid modulators, including agonists, antagonists, and inverse agonists, with all of the different conformational states of the CBRs. The binding effects of these modulators were studied on various CB structural features, including the movement of the transmembrane helices, the volume of the binding cavity, the internal fluids, and the important GPCR properties. Then, using in vitro experiments and computational modeling, we investigated how vitamin E functions as a lipid modulator to influence THC binding. This comparative examination of modulator binding to CBRs provides significant insight into the mechanisms of structural alterations and ligand affinity, which can directly help in the rational design of selective modulators that target either CB₁R or CB₂R.

Large interfacial relocation in RBD-ACE2 complex may explain fast-spreading property of Omicron

The Omicron variant of SARS-CoV-2 emerged in South African in late 2021. This variant has a large number of mutations, and regarded as fastest-spreading Covid variant. The spike RBD region of SARS-CoV-2 and its interaction with human ACE2 play fundamental role in viral infection and transmission. To explore the reason of fast-spreading properties of Omicron variant, we have modeled the interactions of Omicron RBD and human ACE2 using docking and molecular dynamics simulations. Results show that RBD-ACE2 binding site may drastically relocate with an enlarged interface. The predicted interface has large negative binding energies and shows stable conformation in molecular dynamics simulations. It was found that the interfacial area in Omicron RBD-ACE2 complex is increased up to 40% in comparison to wild-type Sars-Cov-2. Moreover, the number of hydrogen bonds significantly increased up to 80%. The key interacting residues become also very different in Omicron variant. The new binding interface can significantly accommodate R403, as a key RBD residue, near ACE2 surface which leads to two new strong salt bridges. The exploration of the new binding interface can help to understand the reasons of high transmission rate of Omicron.

Structure-based virtual screening discovers novel kidney-type glutaminase inhibitors

Glutaminase (GLS) serves a critical bioenergetic role for malignant tumor growth and has become a valuable therapeutic target for cancer treatment. Herein, we performed a structure-based virtual screening to discover novel GLS inhibitors and provide information for developing new GLS inhibitors. We identified critical pharmacological interactions in the GLS1 binding site by analyzing the known GLS1 inhibitors and selected potential inhibitors based on their docking score and pharmacological interactions. The inhibitory effects of compounds were further confirmed by enzymatic and cell viability assays. We treated colorectal cancer and triple-negative breast cancer cells with the selected candidates and measured the inhibitory efficacy of hit compounds on cell viability. In total, we identified three GLS1 inhibitors. The compounds identified from our structure-based virtual screening methodology exhibited great anticancer potential as a lead targeting glutamine metabolism.

Discovery of new chemotypes of dual 5-HT2A/D2 receptor antagonists with a strategy of drug design methodologies

Aim: Through the application of structure- and ligand-based methods, the authors aimed to create an integrative approach to developing a computational protocol for the rational drug design of potent dual 5-HT_2A/D₂ receptor antagonists without off-target activities on H₁ receptors. Materials & methods: Molecular dynamics and virtual docking methods were used to identify key interactions of the structurally diverse antagonists in the binding sites of the studied targets, and to generate their bioactive conformations for further 3D-quantitative structure-activity relationship modeling. Results & conclusion: Toward the goal of finding multi-potent drugs with a more effective and safer profile, the obtained results led to the design of a new set of dual antagonists and opened a new perspective on the therapy for complex brain diseases.

Covalent cannabinoid receptor ligands - structural insight and selectivity challenges

X-ray crystallography and cryogenic electronic microscopy have provided significant advancement in the knowledge of GPCR structure and have allowed the rational design of GPCR ligands. The class A GPCRs cannabinoid receptor type 1 and type 2 are implicated in many pathophysiological processes and thus rational design of drug and tool compounds is of great interest. Recent structural insight into cannabinoid receptors has already led to a greater understanding of ligand binding sites and receptor residues that likely contribute to ligand selectivity. Herein, classes of heterocyclic covalent cannabinoid receptor ligands are reviewed in light of the recent advances in structural knowledge of cannabinoid receptors, with particular discussion regarding covalent ligand selectivity and rationale design.

Novel Noninvasive Approaches to the Treatment of Obesity: From Pharmacotherapy to Gene Therapy

Recent insights into the pathophysiologic underlying mechanisms of obesity have led to the discovery of several promising drug targets and novel therapeutic strategies to address the global obesity epidemic and its comorbidities. Current pharmacologic options for obesity management are largely limited in number and of modest efficacy/safety profile. Therefore, the need for safe and more efficacious new agents is urgent. Drugs that are currently under investigation modulate targets across a broad range of systems and tissues, including the central nervous system, gastrointestinal hormones, adipose tissue, kidney, liver, and skeletal muscle. Beyond pharmacotherapeutics, other potential antiobesity strategies are being explored, including novel drug delivery systems, vaccines, modulation of the gut microbiome, and gene therapy. The present review summarizes the pathophysiology of energy homeostasis and highlights pathways being explored in the effort to develop novel antiobesity medications and interventions but does not cover devices and bariatric methods. Emerging pharmacologic agents and alternative approaches targeting these pathways and relevant research in both animals and humans are presented in detail. Special emphasis is given to treatment options at the end of the development pipeline and closer to the clinic (ie, compounds that have a higher chance to be added to our therapeutic armamentarium in the near future). Ultimately, advancements in our understanding of the pathophysiology and interindividual variation of obesity may lead to multimodal and personalized approaches to obesity treatment that will result in safe, effective, and sustainable weight loss until the root causes of the problem are identified and addressed.

How Do Modulators Affect the Orthosteric and Allosteric Binding Pockets?

Allosteric modulators (AMs) that bind allosteric sites can exhibit greater selectivity than the orthosteric ligands and can either enhance agonist-induced receptor activity (termed positive allosteric modulator or PAM), inhibit agonist-induced activity (negative AM or NAM), or have no effect on activity (silent AM or SAM). Until now, it is not clear what the exact effects of AMs are on the orthosteric active site or the allosteric binding pocket(s). In the present work, we collected both the three-dimensional (3D) structures of receptor-orthosteric ligand and receptor-orthosteric ligand-AM complexes of a specific target protein. Using our novel algorithm toolset, molecular complex characterizing system (MCCS), we were able to quantify the key residues in both the orthosteric and allosteric binding sites along with potential changes of the binding pockets. After analyzing 21 pairs of 3D crystal or cryo-electron microscopy (cryo-EM) complexes, including 4 pairs of GPCRs, 5 pairs of ion channels, 11 pairs of enzymes, and 1 pair of transcription factors, we found that the binding of AMs had little impact on both the orthosteric and allosteric binding pockets. In return, given the accurately predicted allosteric binding pocket(s) of a drug target of medicinal interest, we can confidently conduct the virtual screening or lead optimization without concern that the huge conformational change of the pocket could lead to the low accuracy of virtual screening.

Distinct Binding Mechanisms for Allosteric Sodium Ion in Cannabinoid Receptors

The therapeutic potential of cannabinoid receptors is not fully explored due to psychoactive side effects and lack of selectivity associated with orthosteric ligands. Allosteric modulators have the potential to become selective therapeutics for cannabinoid receptors. Biochemical experiments have shown the effects of the allosteric Na⁺ binding on cannabinoid receptor activity. However, the Na⁺ coordination site and binding pathway are still unknown. Here, we perform molecular dynamic simulations to explore Na⁺ binding in the cannabinoid receptors, CB₁ and CB₂. Simulations reveal that Na⁺ binds to the primary binding site from different extracellular sites for CB₁ and CB₂. A distinct secondary Na⁺ coordination site is identified in CB₁ that is not present in CB₂. Furthermore, simulations also show that intracellular Na⁺ could bind to the Na⁺ binding site in CB₁. Constructed Markov state models show that the standard free energy of Na⁺ binding is similar to the previously calculated free energy for other class A GPCRs.

In Silico Prediction and Validation of CB2 Allosteric Binding Sites to Aid the Design of Allosteric Modulators

Although the 3D structures of active and inactive cannabinoid receptors type 2 (CB2) are available, neither the X-ray crystal nor the cryo-EM structure of CB2-orthosteric ligand-modulator has been resolved, prohibiting the drug discovery and development of CB2 allosteric modulators (AMs). In the present work, we mainly focused on investigating the potential allosteric binding site(s) of CB2. We applied different algorithms or tools to predict the potential allosteric binding sites of CB2 with the existing agonists. Seven potential allosteric sites can be observed for either CB2-CP55940 or CB2-WIN 55,212-2 complex, among which sites B, C, G and K are supported by the reported 3D structures of Class A GPCRs coupled with AMs. Applying our novel algorithm toolset-MCCS, we docked three known AMs of CB2 including Ec2la (C-2), trans-β-caryophyllene (TBC) and cannabidiol (CBD) to each site for further comparisons and quantified the potential binding residues in each allosteric binding site. Sequentially, we selected the most promising binding pose of C-2 in five allosteric sites to conduct the molecular dynamics (MD) simulations. Based on the results of docking studies and MD simulations, we suggest that site H is the most promising allosteric binding site. We plan to conduct bio-assay validations in the future.

New insights into the function of Fascin in actin bundling: A combined theoretical and experimental study

Fascin, one of actin bundling proteins, plays an important role in the cross-linking of actin filaments (F-actin). Phosphorylation of Fascin is an important posttranslational modification to affect its structure and function. For example, a phosphomimetic mutation of Fascin-S39D decrease its bundling ability with F-actin significantly. In this paper, we studied the actin-bundling activity of Fascin by using molecular dynamics (MD) simulations and biochemical methods. All single-site mutations from serine/threonine to aspartic acid were mimicked by MD simulations. For five mutants (S146D, S156D, S218D, T239D and S259D), the mutated residues in domain 2 of Fascin were found to form salt-bridge interactions with an adjacent residue, indicating that mutations of these residues could potentially reduce actin-bundling activity. Further, F-actin-bundling assays and immunofluorescence technique showed S146D and T239D to have a strong effect on Fascin bundling with F-actin. Finally, we show that single-site mutations do not change the general shape of Fascin, but local structures near the mutated residues in Fascin-S146D and T239D become unstable, thereby affecting the ability of Fascin to bind with F-actin. These findings suggest that targeting domain 2 of Fascin would be very useful for the drug design. In addition, our study indicates that MD simulation is a useful method to screening which residues on Fascin are important.

MCCS, a novel characterization method for protein-ligand complex

Delineating the fingerprint or feature vector of a receptor/protein will facilitate the structural and biological studies, as well as the rational design and development of drugs with high affinities and selectivity. However, protein is complicated by its different functional regions that can bind to some of its protein partner(s), substrate(s), orthosteric ligand(s) or allosteric modulator(s) where cogent methods like molecular fingerprints do not work well. We here elaborate a scoring-function-based computing protocol Molecular Complex Characterizing System to help characterize the binding feature of protein-ligand complexes. Based on the reported receptor-ligand interactions, we first quantitate the energy contribution of each individual residue which may be an alternative of MD-based energy decomposition. We then construct a vector for the energy contribution to represent the pattern of the ligand recognition at a receptor and qualitatively analyze the matching level with other receptors. Finally, the energy contribution vector is explored for extensive use in similarity and clustering. The present work provides a new approach to cluster proteins, a perspective counterpart for determining the protein characteristics in the binding, and an advanced screening technique where molecular docking is applicable.

Medicinal Chemistry approach, pharmacology and neuroprotective benefits of CB2R modulators in neurodegenerative diseases

In the last decades, cannabinoid receptor 2 (CB₂R) has continued to receive attention as a key therapeutic target in neuroprotection. Indeed, several findings highlight the neuroprotective effects of CB₂R through suppression of both neuronal excitability and reactive microglia. Additionally, CB₂R seems to be a more promising target than cannabinoid receptor 1 (CB₁R) thanks to the lack of central side effects, its lower expression levels in the central nervous system (CNS), and its inducibility, since its expression enhances quickly in the brain following pathological conditions. This review aims to provide a thorough overview of the main natural and synthetic selective CB₂R modulators, their chemical classification and their potential therapeutic usefulness in neuroprotection, a crucial aspect for the treatment of neurodegenerative diseases.

Molecular Dynamics Simulations for the Determination of the Characteristic Structural Differences between Inactive and Active States of Wild Type and Mutants of the Orexin2 Receptor

The orexin2 receptor (OX2R), which is classified as a class A G protein-coupled receptor (GPCR), is the target of our study. We performed over 20 several-microsecond-scale molecular dynamics simulations of the wild type and mutants of OX2R to extract the characteristics of the structural changes taking place in the active state. We introduced mutations that exhibited the stable inactive state and the constitutively active state in class A GPCRs. In these simulations, significant characteristic structural changes were observed in the V309^6.40Y mutant, which corresponded to a constitutively active mutant. These conformational changes include the outward movement of the transmembrane helix 6 (TM6) and the inward movement of TM7, which are common structural changes in the activation of GPCRs. In addition, we extracted a suitable index for the quantitative evaluation of the active and inactive states of GPCRs, namely, the inter-atomic distance of Cα atoms between x(3.46) and Y(7.53). The structures of the inactive and active states solved by X-ray crystallography and cryo-electron microscopy can be classified using the inter-atomic distance. Furthermore, we clarified that the inward movement of TM7 requires the swapping of M305^6.36 on TM6 and L367^7.56 on TM7. Finally, we discussed the structural advantages of TM7 inward movement for GPCR activation.

Rising of LOXHD1 as a signature causative gene of down-sloping hearing loss in people in their teens and 20s

BACKGROUND

Down-sloping sensorineural hearing loss (SNHL) in people in their teens and 20s hampers efficient learning and communication and in-depth social interactions. Nonetheless, its aetiology remains largely unclear, with the exception of some potential causative genes, none of which stands out especially in people in their teens and 20s. Here, we examined the role and genotype-phenotype correlation of lipoxygenase homology domain 1 (LOXHD1) in down-sloping SNHL through a cohort study.

METHODS

Based on the Seoul National University Bundang Hospital (SNUBH) genetic deafness cohort, in which the patients show varying degrees of deafness and different onset ages (n=1055), we have established the 'SNUBH Teenager-Young Adult Down-sloping SNHL' cohort (10-35 years old) (n=47), all of whom underwent exome sequencing. Three-dimensional molecular modelling, minigene splicing assay and short tandem repeat marker genotyping were performed, and medical records were reviewed.

RESULTS

LOXHD1 accounted for 33.3% of all genetically diagnosed cases of down-sloping SNHL (n=18) and 12.8% of cases in the whole down-sloping SNHL cohort (n=47) of young adults. We identified a potential common founder allele, as well as an interesting genotype-phenotype correlation. We also showed that transcript 6 is necessary and probably sufficient for normal hearing.

CONCLUSIONS

LOXHD1 exceeds other genes in its contribution to down-sloping SNHL in young adults, rising as a signature causative gene, and shows a potential but interesting genotype-phenotype correlation.

Variously substituted 2-oxopyridine derivatives: Extending the structure-activity relationships for allosteric modulation of the cannabinoid CB2 receptor

We previously reported the 2-oxopyridine-3-carboxamide derivative EC21a as the first small synthetic CB2R positive allosteric modulator which displayed antinociceptive activity in vivo in an experimental mouse model of neuropathic pain. Herein, we extended the structure-activity relationships of EC21a through structural modifications regarding the p-fluoro benzyl moiety at position 1 and the amide group in position 3 of the central core. The characterization in vitro was assessed through radioligand binding experiments and functional assays (GTPγS, cAMP, βarrestin2). Among the new compounds, the derivatives A1 (SV-10a) and A5 (SB-13a) characterized respectively by fluorine atom or by chlorine atom in ortho position of the benzylic group at position 1 and by a cycloheptane-carboxamide at position 3 of the central core, showed positive allosteric behavior on CB2R. They enhanced the efficacy of CP55,940 in [³⁵S]GTPγS assay, and modulated CP55,940-dependent βarrestin2 recruitment and cAMP inhibition. The obtained results extend our knowledge of the structural requirements for interaction with the allosteric site of CB2R.

The discovery of a potent and selective pyrazolo-[2,3-e]-[1,2,4]-triazine cannabinoid type 2 receptor agonist

The development of selective CB₂ receptor agonists is a promising therapeutic approach for the treatment of inflammatory diseases, without CB₁ receptor mediated psychoactive side effects. Preliminary structure-activity relationship studies on pyrazoylidene benzamide agonists revealed the -ylidene benzamide moiety was crucial for functional activity at the CB₂ receptor. A small library of compounds with varying linkage moieties between the pyrazole and substituted phenyl group has culminated in the discovery of a potent and selective pyrazolo-[2,3-e]-[1,2,4]-triazine agonist 19 (CB₂R EC₅₀ = 19 nM, CB₁R EC₅₀ > 10 μM). Docking studies have revealed key structural features of the linkage group that are important for potent functional activity.

Cannabinoids and Cannabinoid Receptors: The Story so Far

Like most modern molecular biology and natural product chemistry, understanding cannabinoid pharmacology centers around molecular interactions, in this case, between the cannabinoids and their putative targets, the G-protein coupled receptors (GPCRs) cannabinoid receptor 1 (CB₁) and cannabinoid receptor 2 (CB₂). Understanding the complex structure and interplay between the partners in this molecular dance is required to understand the mechanism of action of synthetic, endogenous, and phytochemical cannabinoids. This review, with 91 references, surveys our understanding of the structural biology of the cannabinoids and their target receptors including both a critical comparison of the extant crystal structures and the computationally derived homology models, as well as an in-depth discussion about the binding modes of the major cannabinoids. The aim is to assist in situating structural biochemists, synthetic chemists, and molecular biologists who are new to the field of cannabis research.

•

Cannabinoid research has greatly expanded

•

Structural biology and computational chemistry jointly provide mechanistic insight

•

Structural data are being generated at an exponentially increasing rate

•

Phytocannabinoid targeting of other GPCR receptors deserves investigation

Cutaneous endocannabinoid system: Does it have a role on skin wound healing bearing fibrosis?

INTRODUCTION

Recently, the endocannabinoid system has been identified in skin and it has been linked with the formation of skin fibrosis and wound healing. We aimed to find and analyse reported data on compounds acting in the endocannabinoid system with significant effect in skin fibrosis.

METHODS

A literature search on PUBMED was conducted for studies published in English until February 2020 on cannabinoids and skin fibrosis. The initial search was performed with terms: "cannabinoid" AND "skin". This search retrieved 296 publications from which 18 directly related to skin fibrosis or wound healing process were included in this review.

RESULTS

CB1 receptor inactivation and CB2 receptor activation show anti-fibrotic effects on cellular and animal experimental models of cutaneous fibrosis. CB2 receptor activation also promotes re-epithelization. Other cannabinoid related receptors, like adenosine A2A receptors and PPAR-γ, are also involved. Their activation lead to a pro-fibrotic and anti-fibrotic effect, respectively.

CONCLUSION

Several molecular drug targets for endocannabinoid system were identified in skin. It may be a promising approach for the treatment of excessive skin fibrosis disorders.

Allosteric modulators targeting cannabinoid cb1 and cb2 receptors: implications for drug discovery

Allosteric modulators of cannabinoid receptors hold great therapeutic potential, as they do not possess intrinsic efficacy, but instead enhance or diminish the receptor's response of orthosteric ligands allowing for the tempering of cannabinoid receptor signaling without the desensitization, tolerance and dependence. Allosteric modulators of cannabinoid receptors have numerous advantages over the orthosteric ligands such as higher receptor type selectivity, probe dependence and biased signaling, so they have a great potential to separate the therapeutic benefits from side effects own of orthosteric ligands. This review aims to give an overview of the CB1 and CB2 receptor allosteric modulators highlighting the structure-activity relationship and pharmacological profile of each classes, and their future promise.

New Insights in Cannabinoid Receptor Structure and Signaling

BACKGROUND

Cannabinoid has long been used for medicinal purposes. Cannabinoid signaling has been considered the therapeutic target for treating pain, addiction, obesity, inflammation, and other diseases. Recent studies have suggested that in addition to CB1 and CB2, there are non-CB1 and non-CB2 cannabinoid-related orphan GPCRs including GPR18, GPR55, and GPR119. In addition, CB1 and CB2 display allosteric binding and biased signaling, revealing correlations between biased signaling and functional outcomes. Interestingly, new investigations have indicated that CB1 is functionally present within the mitochondria of striated and heart muscles directly regulating intramitochondrial signaling and respiration.

CONCLUSION

In this review, we summarize the recent progress in cannabinoid-related orphan GPCRs, CB1/CB2 structure, Gi/Gs coupling, allosteric ligands and biased signaling, and mitochondria-localized CB1, and discuss the future promise of this research.

Paeonol alleviates primary dysmenorrhea in mice via activating CB2R in the uterus

BACKGROUND AND PURPOSE

Primary dysmenorrhea is the most common gynaecologic problem in menstruating women and is characterized by spasmodic uterine contraction and pain symptoms associated with inflammatory disturbances. Paeonol is an active phytochemical component that has shown anti-inflammatory and analgesic effects in several animal models. The aim of this study was to explore whether paeonol is effective against dysmenorrhea and to investigate the potential mechanism of cannabinoid receptor signalling.

EXPERIMENTAL APPROACH

Dysmenorrhea was established by injecting oestradiol benzoate into female mice. The effects of paeonol on writhing time and latency, uterine pathology and inflammatory mediators were explored. Isolated uterine smooth muscle was used to evaluate the direct effect of paeonol on uterine contraction.

KEY RESULTS

The oral administration of paeonol reduced dysmenorrhea pain and PGE2 and TNF-α expression in the uterine tissues of mice, and paeonol was found to be distributed in lesions of the uterus. Paeonol almost completely inhibited oxytocin-, high potassium- and Ca²⁺-induced contractions in isolated uteri. Antagonists of CB2R (AM630) and the MAPK pathway (U0126), but not of CB1R (AM251), reversed the inhibitory effect of paeonol on uterine contraction. Paeonol significantly blocked L-type Ca²⁺ channels and calcium influx in uterine smooth muscle cells via CB2R. Molecular docking results showed that paeonol fits well with the binding site of CB2R.

CONCLUSIONS AND IMPLICATIONS

Paeonol partially acts through CB2R to restrain calcium influx and uterine contraction to alleviate dysmenorrhea in mice. These results suggest that paeonol has therapeutic potential for the treatment of dysmenorrhea.

Insight into Ginkgo biloba L. Extract on the Improved Spatial Learning and Memory by Chemogenomics Knowledgebase, Molecular Docking, Molecular Dynamics Simulation, and Bioassay Validations

Epilepsy is a common cause of serious cognitive disorders and is known to have impact on patients' memory and executive functions. Therefore, the development of antiepileptic drugs for the improvement of spatial learning and memory in patients with epileptic cognitive dysfunction is important. In the present work, we systematically predicted and analyzed the potential effects of Ginkgo terpene trilactones (GTTL) on cognition and pathologic changes utilizing in silico and in vivo approaches. Based on our established chemogenomics knowledgebase, we first conducted the network systems pharmacology analysis to predict that ginkgolide A/B/C may target 5-HT 1A, 5-HT 1B, and 5-HT 2B. The detailed interactions were then further validated by molecular docking and molecular dynamics (MD) simulations. In addition, status epilepticus (SE) was induced by lithium-pilocarpine injection in adult Wistar male rats, and the results of enzyme-linked immunosorbent assay (ELISA) demonstrated that administration with GTTL can increase the expression of brain-derived neurotrophic factor (BDNF) when compared to the model group. Interestingly, recent studies suggest that the occurrence of a reciprocal involvement of 5-HT receptor activation along with the hippocampal BDNF-increased expression can significantly ameliorate neurologic changes and reverse behavioral deficits in status epilepticus rats while improving cognitive function and alleviating neuronal injury. Therefore, we evaluated the effects of GTTL (bilobalide, ginkgolide A, ginkgolide B, and ginkgolide C) on synergistic antiepileptic effect. Our experimental data showed that the spatial learning and memory abilities (e.g., electroencephalography analysis and Morris water maze test for behavioral assessment) of rats administrated with GTTL were significantly improved under the middle dose (80 mg/kg, GTTL) and high dose (160 mg/kg, GTTL). Moreover, the number of neurons in the hippocampus of the GTTL group increased when compared to the model group. Our studies showed that GTTL not only protected rat cerebral hippocampal neurons against epilepsy but also improved the learning and memory ability. Therefore, GTTL may be a potential drug candidate for the prevention and/or treatment of epilepsy.

Structure-Based Design of Novel Biphenyl Amide Antagonists of Human Transient Receptor Potential Cation Channel Subfamily M Member 8 Channels with Potential Implications in the Treatment of Sensory Neuropathies

Structure-activity relationship studies of a reported menthol-based transient receptor potential cation channel subfamily M member 8 channel (TRPM8) antagonist, guided by computational simulations and structure-based design, uncovers a novel series of TRPM8 antagonists with >10-fold selectivity versus related TRP subtypes. Spiro[4.5]decan-8-yl analogue 14 inhibits icilin-evoked Ca²⁺ entry in HEK-293 cells stably expressing human TRPM8 (hTRPM8) with an IC₅₀ of 2.4 ± 1.0 nM, while in whole-cell patch-clamp recordings this analogue inhibits menthol-evoked currents with a hTRPM8 IC₅₀ of 64 ± 2 nM. Molecular dynamics (MD) simulations of compound 14 in our homology model of hTRPM8 suggest that this antagonist forms extensive hydrophobic contacts within the orthosteric site. In the wet dog shakes (WDS) assay, compound 14 dose-dependently blocks icilin-triggered shaking behaviors in mice. Upon local administration, compound 14 dose dependently inhibits cold allodynia evoked by the chemotherapy oxaliplatin in a murine model of peripheral neuropathy at microgram doses. Our findings suggest that 14 and other biphenyl amide analogues within our series can find utility as potent antagonist chemical probes derived from (-)-menthol as well as small molecule therapeutic scaffolds for chemotherapy-induced peripheral neuropathy (CIPN) and other sensory neuropathies.

Cryo-EM Structure of the Human Cannabinoid Receptor CB2-Gi Signaling Complex

Drugs selectively targeting CB2 hold promise for treating neurodegenerative disorders, inflammation, and pain while avoiding psychotropic side effects mediated by CB1. The mechanisms underlying CB2 activation and signaling are poorly understood but critical for drug design. Here we report the cryo-EM structure of the human CB2-G_i signaling complex bound to the agonist WIN 55,212-2. The 3D structure reveals the binding mode of WIN 55,212-2 and structural determinants for distinguishing CB2 agonists from antagonists, which are supported by a pair of rationally designed agonist and antagonist. Further structural analyses with computational docking results uncover the differences between CB2 and CB1 in receptor activation, ligand recognition, and G_i coupling. These findings are expected to facilitate rational structure-based discovery of drugs targeting the cannabinoid system.

Molecular dynamics simulations for genetic interpretation in protein coding regions: where we are, where to go and when

The increasing ease with which massive genetic information can be obtained from patients or healthy individuals has stimulated the development of interpretive bioinformatics tools as aids in clinical practice. Most such tools analyze evolutionary information and simple physical-chemical properties to predict whether replacement of one amino acid residue with another will be tolerated or cause disease. Those approaches achieve up to 80-85% accuracy as binary classifiers (neutral/pathogenic). As such accuracy is insufficient for medical decision to be based on, and it does not appear to be increasing, more precise methods, such as full-atom molecular dynamics (MD) simulations in explicit solvent, are also discussed. Then, to describe the goal of interpreting human genetic variations at large scale through MD simulations, we restrictively refer to all possible protein variants carrying single-amino-acid substitutions arising from single-nucleotide variations as the human variome. We calculate its size and develop a simple model that allows calculating the simulation time needed to have a 0.99 probability of observing unfolding events of any unstable variant. The knowledge of that time enables performing a binary classification of the variants (stable-potentially neutral/unstable-pathogenic). Our model indicates that the human variome cannot be simulated with present computing capabilities. However, if they continue to increase as per Moore's law, it could be simulated (at 65°C) spending only 3 years in the task if we started in 2031. The simulation of individual protein variomes is achievable in short times starting at present. International coordination seems appropriate to embark upon massive MD simulations of protein variants.

Structural insight into the serotonin (5-HT) receptor family by molecular docking, molecular dynamics simulation and systems pharmacology analysis

Serotonin (5-HT) receptors are proteins involved in various neurological and biological processes, such as aggression, anxiety, appetite, cognition, learning, memory, mood, sleep, and thermoregulation. They are commonly associated with drug abuse and addiction due to their importance as targets for various pharmaceutical and recreational drugs. However, due to a high sequence similarity/identity among 5-HT receptors and the unavailability of the 3D structure of the different 5-HT receptor, no report was available so far regarding the systematical comparison of the key and selective residues involved in the binding pocket, making it difficult to design subtype-selective serotonergic drugs. In this work, we first built and validated three-dimensional models for all 5-HT receptors based on the existing crystal structures of 5-HT_1B, 5-HT_2B, and 5-HT_2C. Then, we performed molecular docking studies between 5-HT receptors agonists/inhibitors and our 3D models. The results from docking were consistent with the known binding affinities of each model. Sequentially, we compared the binding pose and selective residues among 5-HT receptors. Our results showed that the affinity variation could be potentially attributed to the selective residues located in the binding pockets. Moreover, we performed MD simulations for 12 5-HT receptors complexed with ligands; the results were consistent with our docking results and the reported data. Finally, we carried out off-target prediction and blood-brain barrier (BBB) prediction for Captagon using our established hallucinogen-related chemogenomics knowledgebase and in-house computational tools, with the hope to provide more information regarding the use of Captagon. We showed that 5-HT_2C, 5-HT_5A, and 5-HT₇ were the most promising targets for Captagon before metabolism. Overall, our findings can provide insights into future drug discovery and design of medications with high specificity to the individual 5-HT receptor to decrease the risk of addiction and prevent drug abuse.

Discovery of novel natural compound inhibitors targeting estrogen receptor α by an integrated virtual screening strategy

Estrogen receptor (ER) is a nuclear hormone receptor and plays an important role in mediating the cellular effects of estrogen. ER can be classified into two receptors: estrogen receptor alpha (ERα) and beta (ERβ), and the former is expressed in 50~80% of breast tumors and has been extensively investigated in breast cancer for decades. Excessive exposure to estrogen can obviously stimulate the growth of breast cancers primarily mediated by ERα, and thus anti-estrogen therapies by small molecules are of concern to clinicians and pharmaceutical industry in the treatment of ERα-positive breast cancers. Although a series of estrogen receptor modulators have been developed, these drugs can lead to resistance and side effects. Therefore, the development of small molecule inhibitors with high target specificity has been intensified. In this pursuit, an integrated computer-aided virtual screening technique, including molecular docking and pharmacophore model screening, was used to screen traditional Chinese medicine (TCM) databases. The compounds with high docking scores and fit values were subjected to ADME (adsorption, distribution, metabolism, excretion) and toxicity prediction, and ten hits were identified as potential inhibitors targeting ERα. Molecular docking was used to investigate the binding modes between ERα and three most potent hits, and molecular dynamic simulations were chosen to explore the stability of these complexes. The rank of the predicted binding free energies evaluated by MM/GBSA is consistent with the docking score. These novel scaffolds discovered in the present study can be used as critical starting point in the drug discovery process for treating ERα-positive breast cancer. Graphical abstract .

Computational Systems Pharmacology-Target Mapping for Fentanyl-Laced Cocaine Overdose

The United States of America is fighting against one of its worst-ever drug crises. Over 900 people a week die from opioid- or heroin-related overdoses, while millions more suffer from opioid prescription addiction. Recently, drug overdoses caused by fentanyl-laced cocaine specifically are on the rise. Due to drug synergy and an increase in side effects, polydrug addiction can cause more risk than addiction to a single drug. In the present work, we systematically analyzed the overdose and addiction mechanism of cocaine and fentanyl. First, we applied our established chemogenomics knowledgebase and machine-learning-based methods to map out the potential and known proteins, transporters, and metabolic enzymes and the potential therapeutic target(s) for cocaine and fentanyl. Sequentially, we looked into the detail of (1) the addiction to cocaine and fentanyl by binding to the dopamine transporter and the μ opioid receptor (DAT and μOR, respectively), (2) the potential drug-drug interaction of cocaine and fentanyl via p-glycoprotein (P-gp) efflux, (3) the metabolism of cocaine and fentanyl in CYP3A4, and (4) the physiologically based pharmacokinetic (PBPK) model for two drugs and their drug-drug interaction at the absorption, distribution, metabolism, and excretion (ADME) level. Finally, we looked into the detail of JWH133, an agonist of cannabinoid 2-receptor (CB2) with potential as a therapy for cocaine and fentanyl overdose. All these results provide a better understanding of fentanyl and cocaine polydrug addiction and future drug abuse prevention.

DAKB-GPCRs: An Integrated Computational Platform for Drug Abuse Related GPCRs

Drug abuse (DA) or drug addiction is a complicated brain disorder which is commonly considered as neurobiological impairments caused by both genetic factors and environmental effects. Among DA-related targets, G protein-coupled receptors (GPCRs) play an important role in DA therapy. However, only 52 GPCRs have been published with crystal structures in the recent two decades. In the effort to overcome the limitations of crystal structure and conformational diversity of GPCRs, we built homology models and performed conformational searches by molecular dynamics (MD) simulation. To accelerate and facilitate the drug abuse research, we construct a DA-related GPCR-specific chemogenomics knowledgebase (KB) (DAKB-GPCRs) for its research that can be implemented with our established and novel chemogenomics tools as well as algorithms for data analysis and visualization. Our established TargetHunter and HTDocking tools, as well as our novel tools that include target classification and Spider Plot, are compiled into the platform. Our DAKB-GPCRs provides the following results for a query compound: (1) blood-brain barrier (BBB) plot via our BBB predictor, (2) docking scores via HTDocking, (3) similarity score via TargetHunter, (4) target classification via machine learning methods that utilize both docking scores and similarity scores, and (5) a drug-target interaction network via Spider Plot.

Insight of Captagon Abuse by Chemogenomics Knowledgebase-guided Systems Pharmacology Target Mapping Analyses

Captagon, known by its genetic name Fenethylline, is an addictive drug that complicates the War on Drugs. Captagon has a strong CNS stimulating effect than its primary metabolite, Amphetamine. However, multi-targets issues associated with the drug and metabolites as well as its underlying mechanisms have not been fully defined. In the present work, we applied our established drug-abuse chemogenomics-knowledgebase systems pharmacology approach to conduct targets/off-targets mapping (SP-Targets) investigation of Captagon and its metabolites for hallucination addiction, and also analyzed the cell signaling pathways for both Amphetamine and Theophylline with data mining of available literature. Of note, Amphetamine, an agonist for trace amine-associated receptor 1 (TAAR1) with enhancing dopamine signaling (increase of irritability, aggression, etc.), is the main cause of Captagon addiction; Theophylline, an antagonist that blocks adenosine receptors (e.g. A2aR) in the brain responsible for restlessness and painlessness, may attenuate the behavioral sensitization caused by Amphetamine. We uncovered that Theophylline's metabolism and elimination could be retarded due to competition and/or blockage of the CYP2D6 enzyme by Amphetamine; We also found that the synergies between these two metabolites cause Captagon's psychoactive effects to act faster and far more potently than those of Amphetamine alone. We carried out further molecular docking modeling and molecular dynamics simulation to explore the molecular interactions between Amphetamine and Theophylline and their important GPCRs targets, including TAAR1 and adenosine receptors. All of the systems pharmacology analyses and results will shed light insight into a better understanding of Captagon addiction and future drug abuse prevention.

Targeted inhibition of the type 2 cannabinoid receptor is a novel approach to reduce renal fibrosis

The cannabinoid receptor type 2 (CB2) is a G protein-coupled seven transmembrane receptor that transmits endogenous cannabinoid signaling. The role of CB2 in the pathogenesis of kidney injury and fibrosis remains poorly understood. Here we demonstrate that CB2 was induced, predominantly in kidney tubular epithelium, in various models of kidney disease induced by unilateral ureteral obstruction, adriamycin or ischemia/reperfusion injury. In vitro, forced expression of CB2 or treatment with a CB2 agonist was sufficient to trigger matrix gene expression, whereas knockdown of CB2 by siRNA abolished transforming growth factor-β1-induced signaling and fibrogenic responses in kidney tubular cells. CB2 also mediated fibroblasts and macrophage activation in vitro. Mice with genetic ablation of CB2 were protected against kidney injury after ureteral obstruction, validating a pathogenic role of CB2 in renal fibrosis in vivo. By using in silico screening and medicinal chemistry modifications, we discovered a novel compound, XL-001, that bound to CB2 with high affinity and selectivity and acted as an inverse agonist. Incubation with XL-001 inhibited in a dose-dependent fashion the fibrogenic response induced by CB2 overexpression, CB2 agonist or transforming growth factor-β1. In vivo, intraperitoneal injections of XL-001 ameliorated kidney injury, fibrosis and inflammation in both the obstruction and ischemia/reperfusion models. Delayed administration of XL-001 was also effective in ameliorating kidney fibrosis and inflammation. Thus, CB2 is a pathogenic mediator in kidney fibrosis and targeted inhibition with the novel inverse agonist XL-001 may provide a strategy in the fight against fibrotic kidney diseases.

An insight into paracetamol and its metabolites using molecular docking and molecular dynamics simulation

Paracetamol is a relatively safe analgesia/antipyretic drug without the risks of addiction, dependence, tolerance, and withdrawal when used alone. However, when administrated in an opioid/paracetamol combination product, which often contains a large quantity of paracetamol, it can be potentially dangerous due to the risk of hepatotoxicity. Paracetamol is known to be metabolized into N-(4-hydroxyphenyl)-arachidonamide (AM404) via fatty acid amide hydrolase (FAAH) and into N-acetyl-p-benzoquinone imine (NAPQI) via cytochrome P450 (CYP) enzymes. However, the underlying mechanism of paracetamol is still unclear. In addition, paracetamol has the potential to interact with other drugs that are also involved with CYP family enzymes (inducer/inhibitor/substrate), an example being illicit drugs. In our present work, we looked into the relationship between paracetamol and its metabolites (AM404 and NAPQI) using molecular docking and molecular dynamics (MD) simulations. We first carried out a series of molecular docking studies between paracetamol/AM404/NAQPI and their reported targets, including CYP 2E1, FAAH, TRPA1, CB1, and TRPV1. Subsequently, we performed MD simulations and energy decomposition for CB1-AM404, TRPV1-AM404, and TRPV1-NAPQI for further investigation of the dynamics interactions. Finally, we summarized and discussed the reported drug-drug interactions between paracetamol and central nervous system drugs, especially illicit drugs. Overall, we are able to provide new insights into the structural and functional roles of paracetamol and its metabolites that can inform the potential prevention and treatment of paracetamol overdose. Graphical abstract Paracetamol and its metabolites.

Emerging strategies targeting CB2 cannabinoid receptor: Biased agonism and allosterism

During these last years, the CB₂ cannabinoid receptor has emerged as a potential anti-inflammatory target in diseases such as multiple sclerosis, amyotrophic lateral sclerosis, Huntington's disease, ischemic stroke, autoimmune diseases, osteoporosis, and cancer. However, the development of clinically useful CB₂ agonists reveals to be very challenging. Allosterism and biased-signaling mechanisms at CB₂ receptor may offer new avenues for the development of improved CB₂ receptor-targeted therapies. Although there has been some exploration of CB₁ receptor activation by new CB₁ allosteric or biased-signaling ligands, the CB₂ receptor is still at initial stages in this domain. In an effort to understand the molecular basis behind these pharmacological approaches, we have analyzed and summarized the structural data reported so far at CB₂ receptor.

Identification of the First Synthetic Allosteric Modulator of the CB2 Receptors and Evidence of Its Efficacy for Neuropathic Pain Relief

The direct activation of cannabinoid receptors (CBRs) results in several beneficial effects; therefore several CBRs ligands have been synthesized and tested in vitro and in vivo. However, none of them reached an advanced phase of clinical development due mainly to side effects on the CNS. Medicinal chemistry approaches are now engaged to develop allosteric modulators that might offer a novel therapeutic approach to achieve potential therapeutic benefits avoiding inherent side effects of orthosteric ligands. Here we identify the first ever synthesized positive allosteric modulator (PAM) that targets CB₂Rs. The evidence for this was obtained using [³H]CP55940 and [³⁵S]GTPγS binding assays. This finding will be useful for the characterization of allosteric binding site(s) on CB₂Rs which will be essential for the further development of CB₂R allosteric modulators. Moreover, the new CB₂R PAM displayed antinociceptive activity in vivo in an experimental mouse model of neuropathic pain, raising the possibility that it might be a good candidate for clinical development.

The cyclohexene derivative MC-3129 exhibits antileukemic activity via RhoA/ROCK1/PTEN/PI3K/Akt pathway-mediated mitochondrial translocation of cofilin

The effects of MC-3129, a synthetic cyclohexene derivative, on cell viability and apoptosis have been investigated in human leukemia cells. Exposure of leukemia cells to MC-3129 led to the inhibition of cell viability and induction of apoptosis through the dephosphorylation and mitochondrial translocation of cofilin. A mechanistic study revealed that interruption of the RhoA/ROCK1/PTEN/PI3K/Akt signaling pathway plays a crucial role in the MC-3129-mediated dephosphorylation and mitochondrial translocation of cofilin and induction of apoptosis. Our in vivo study also showed that the MC-3129-mediated inhibition of the tumor growth in a mouse leukemia xenograft model is associated with the interruption of ROCK1/PTEN/PI3K/Akt signaling and apoptosis. Molecular docking suggested that MC-3129 might activate the RhoA/ROCK1 pathway by targeting LPAR2. Collectively, these findings suggest a hierarchical model, in which the induction of apoptosis by MC-3129 primarily results from the activation of RhoA/ROCK1/PTEN and inactivation of PI3K/Akt, leading to the dephosphorylation and mitochondrial translocation of cofilin, and culminating in cytochrome c release, caspase activation, and apoptosis. Our study reveals a novel role for RhoA/ROCK1/PTEN/PI3K/Akt signaling in the regulation of mitochondrial translocation of cofilin and apoptosis and suggests MC-3129 as a potential drug for the treatment of human leukemia.

Insight into tetrodotoxin blockade and resistance mechanisms of Nav 1.2 sodium channel by theoretical approaches

Na_v 1.2, a member of voltage-gated sodium channels (Na_v s) that are responsible for the generation and propagation of action potentials along the cell membrane, and play a vital role in the process of information transmission within the nervous system and muscle contraction, is preferentially expressed in the central nervous system. As a potent and selective blocker of Na_v s, tetrodotoxin (TTX) has been extensively studied in biological and chemical sciences, whereas the detailed mechanism by which it blocks nine Na_v 1 channel subtypes remain elusive. Despite the high structural similarity, the TTX metabolite 4,9-anhydro-TTX is 161 times less effective toward the mammalian Na_v 1.2, which puzzled us to ask a question why such a subtle structural variation results in the largely binding affinity difference. In the current work, an integrated computational strategy, including homology modeling, induced fit docking, explicit-solvent MD simulations, and free energy calculations, was employed to investigate the binding mechanism and conformational determinants of TTX analogs. Based on the computational results, the H-bond interactions between C4-OH and C9-OH of TTX and the outer ring carboxylates of the selectivity-filter residues, and the cation-π interaction between the primary amine of guanidinium of TTX and Phe385 determine the difference of their binding affinities. Moreover, the computationally simulations were carried out for the D384N and E945K mutants of hNa_v 1.2-TTX, and the rank of the predicted binding free energies is in accordance with the experimental data. These observations provide a valuable model to design potent and selective neurotoxins of Na_v 1.2 and shed light on the blocking mechanism of TTX to sodium channels.

Exploring the mechanism of F282L mutation-caused constitutive activity of GPCR by a computational study

G-protein-coupled receptors (GPCRs) are important drug targets and generally activated by ligands. However, some experiments found that GPCRs also give rise to constitutive activity through some mutations (viz., CAM), which are usually associated with different kinds of diseases. However, the mechanisms of CAMs and their roles in interactions with drug-ligands are unclear in experiments. Herein, we used microsecond molecular dynamics simulations to study the effect of one important F282L mutation on β2AR in order to address the questions above. With the aid of principle component and correlation analysis, our results revealed that the F282L mutation could increase the instability of the overall structure, increase the dramatic fluctuations of NPxxY and extracellular loops, and decrease restraint of the helices through weakening interhelical H-bonding and correlations between residues, which could partly contribute to the constitutive activity reported by the experiments. The observations from the protein structure network (PSN) analysis indicate that the mutant exhibits less information flow than the wild β2AR and weakens the role of TM5 and TM6 in the signal transmission, but it enhances the impact of TM3 on the orthosteric pathway and TM4 on the allosteric one. In addition, the results from the virtual screening reveal that the mutant prefers to select agonists rather than antagonists, similar to the active state but opposite of the inactive state, further confirming that the F282L mutation advances the activation of β2AR. Our observations provide valuable information for understanding the mechanism of the mutation-caused constitutive activity of GPCR and related drug-design.

Allosteric Modulation of Intact γ-Secretase Structural Dynamics

As a protease complex involved in the cleavage of amyloid precursor proteins that lead to the formation of amyloid β fibrils implicated in Alzheimer's disease, γ-secretase is an important target for developing therapeutics against Alzheimer's disease. γ-secretase is composed of four subunits: nicastrin (NCT) in the extracellular (EC) domain, presenilin-1 (PS1), anterior pharynx defective 1, and presenilin enhancer 2 in the transmembrane (TM) domain. NCT and PS1 play important roles in binding amyloid β precursor proteins and modulating PS1 catalytic activity. Yet, the molecular mechanisms of coupling between substrate/modulator binding and catalytic activity remain to be elucidated. Recent determination of intact human γ-secretase cryo-electron microscopy structure has opened the way for a detailed investigation of the structural dynamics of this complex. Our analysis, based on a membrane-coupled anisotropic network model, reveals two types of NCT motions, bending and twisting, with respect to PS1. These underlie the fluctuations between the "open" and "closed" states of the lid-like NCT with respect to a hydrophilic loop 1 (HL1) on PS1, thus allowing or blocking access of the substrate peptide (EC portion) to HL1 and to the neighboring helix TM2. In addition to this alternating access mechanism, fluctuations in the volume of the PS1 central cavity facilitate the exposure of the catalytic site for substrate cleavage. Druggability simulations show that γ-secretase presents several hot spots for either orthosteric or allosteric inhibition of catalytic activity, consistent with experimental data. In particular, a hinge region at the interface between the EC and TM domains, near the interlobe groove of NCT, emerges as an allo-targeting site that would impact the coupling between HL1/TM2 and the catalytic pocket, opening, to our knowledge, new avenues for structure-based design of novel allosteric modulators of γ-secretase protease activity.

Exploring the binding mechanisms of PDE5 with chromeno[2,3-c]pyrrol-9(2H)-one by theoretical approaches

Cyclic nucleotide phosphodiesterase type 5 (PDE5), exclusively specific for the cyclic guanosine monophosphate (cGMP), is an important drug target for the treatment of erectile dysfunction and pulmonary arterial hypertension (PAH). Although many PDE5 inhibitors have been approved, such as sildenafil, vardenafil, tadalafil and so on, extensive studies have reported some side effects, such as vision disturbance and hearing loss as a result of the amino acid sequence and the secondary structural similarity of other PDEs to the catalytic domain of PDE5. In this study, multiple docking strategies, molecular dynamics (MD) simulations, free energy calculations and decomposition were employed to explore the structural determinants of PDE5 with a series of chromeno[2,3-c]pyrrol-9(2H)-one derivatives. First, reliable docking results were obtained using quantum mechanics/molecular mechanics (QM/MM) docking. Then, MD simulations and MM/GBSA free energy calculations were used to explore the dynamic binding process and characterize the binding modes of the inhibitors with different activities. The predicted binding free energies are in good agreement with the experimental data, and the MM/GBSA free energy decomposition analysis sheds light on the importance of hydrogen bonds with Gln817, π–π stacks against Phe820 and hydrophobic residues for the PDE5 binding of the studied inhibitors. The structural and energetic insights obtained here are useful for understanding the molecular mechanism of ligand binding and designing novel potent and selective PDE5 inhibitors with new scaffolds.

Cyclic nucleotide phosphodiesterase type 5 (PDE5) is exclusively specific for the cyclic guanosine monophosphate (cGMP), and PDE5 is an important drug target for the treatment of erectile dysfunction and pulmonary arterial hypertension (PAH).

Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β2-Adrenergic Receptor

It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.

Discovery of Selective Cannabinoid CB2 Receptor Agonists by High-Throughput Screening

The endocannabinoid system (ECS) plays a diverse role in human physiology ranging from the regulation of mood and appetite to immune modulation and the response to pain. Drug development that targets the cannabinoid receptors (CB₁ and CB₂) has been explored; however, success in the clinic has been limited by the psychoactive side effects associated with modulation of the neuronally expressed CB₁ that are enriched in the CNS. CB₂, however, are expressed in peripheral tissues, primarily in immune cells, and thus development of CB₂-selective drugs holds the potential to modulate pain among other indications without eliciting anxiety and other undesirable side effects associated with CB₁ activation. As part of a collaborative effort among industry and academic laboratories, we performed a high-throughput screen designed to discover selective agonists or positive allosteric modulators (PAMs) of CB₂. Although no CB₂ PAMs were identified, 167 CB₂ agonists were discovered here, and further characterization of four select compounds revealed two with high selectivity for CB₂ versus CB₁. These results broaden drug discovery efforts aimed at the ECS and may lead to the development of novel therapies for immune modulation and pain management with improved side effect profiles.