Multiple Binding Sites Contribute to the Mechanism of Mixed Agonistic and Positive Allosteric Modulators of the Cannabinoid CB1 Receptor

Characterization of binding kinetics and intracellular signaling of new psychoactive substances targeting cannabinoid receptor using transition-based reweighting method

New psychoactive substances (NPS) targeting cannabinoid receptor 1 pose a significant threat to society as recreational abusive drugs that have pronounced physiological side effects. These greater adverse effects compared to classical cannabinoids have been linked to the higher downstream β-arrestin signaling. Thus, understanding the mechanism of differential signaling will reveal important structure-activity relationship essential for identifying and potentially regulating NPS molecules. In this study, we simulate the slow (un)binding process of NPS MDMB-Fubinaca and classical cannabinoid HU-210 from CB1 using multi-ensemble simulation to decipher the effects of ligand binding dynamics on downstream signaling. The transition-based reweighing method is used for the estimation of transition rates and underlying thermodynamics of (un)binding processes of ligands with nanomolar affinities. Our analyses reveal major interaction differences with transmembrane TM7 between NPS and classical cannabinoids. A variational autoencoder-based approach, neural relational inference (NRI), is applied to assess the allosteric effects on intracellular regions attributable to variations in binding pocket interactions. NRI analysis indicate a heightened level of allosteric control of NPxxY motif for NPS-bound receptors, which contributes to the higher probability of formation of a crucial triad interaction (Y7.53-Y5.58-T3.46) necessary for stronger β-arrestin signaling. Hence, in this work, MD simulation, data-driven statistical methods, and deep learning point out the structural basis for the heightened physiological side effects associated with NPS, contributing to efforts aimed at mitigating their public health impact.

Identifying Residues for Substrate Recognition in Human GPAT4 by Molecular Dynamics Simulations

Glycerol-3-phosphate acyltransferase (GPAT) catalyzes the first step in triacylglycerol synthesis. Understanding its substrate recognition mechanism may help to design drugs to regulate the production of glycerol lipids in cells. In this work, we investigate how the native substrate, glycerol-3-phosphate (G3P), and palmitoyl-coenzyme A (CoA) bind to the human GPAT isoform GPAT4 via molecular dynamics simulations (MD). As no experimentally resolved GPAT4 structure is available, the AlphaFold model is employed to construct the GPAT4-substrate complex model. Using another isoform, GPAT1, we demonstrate that once the ligand binding is properly addressed, the AlphaFold complex model can deliver similar results to the experimentally resolved structure in MD simulations. Following the validated protocol of complex construction, we perform MD simulations using the GPAT4-substrate complex. Our simulations reveal that R427 is an important residue in recognizing G3P via a stable salt bridge, but its motion can bring the ligand to different binding hotspots on GPAT4. Such high flexibility can be attributed to the flexible region that exists only on GPAT4 and not on GPAT1. Our study reveals the substrate recognition mechanism of GPAT4 and hence paves the way towards designing GPAT4 inhibitors.

Uracil/H+ Symport by FurE Refines Aspects of the Rocking-bundle Mechanism of APC-type Transporters

Transporters mediate the uptake of solutes, metabolites and drugs across the cell membrane. The eukaryotic FurE nucleobase/H+ symporter of Aspergillus nidulans has been used as a model protein to address structure-function relationships in the APC transporter superfamily, members of which are characterized by the LeuT-fold and seem to operate by the so-called 'rocking-bundle' mechanism. In this study, we reveal the binding mode, translocation and release pathway of uracil/H+ by FurE using path collective variable, funnel metadynamics and rational mutational analysis. Our study reveals a stepwise, induced-fit, mechanism of ordered sequential transport of proton and uracil, which in turn suggests that FurE, functions as a multi-step gated pore, rather than employing 'rocking' of compact domains, as often proposed for APC transporters. Finally, our work supports that specific residues of the cytoplasmic N-tail are involved in substrate translocation, in line with their essentiality for FurE function.

Novel Cannabinoid Receptor 2 (CB2) Low Lipophilicity Agonists Produce Distinct cAMP and Arrestin Signalling Kinetics without Bias

Cannabinoid Receptor 2 (CB2) is a promising target for treating inflammatory diseases. We designed derivatives of 3-carbamoyl-2-pyridone and 1,8-naphthyridin-2(1H)-one-3-carboxamide CB2-selective agonists with reduced lipophilicity. The new compounds were measured for their affinity (radioligand binding) and ability to elicit cyclic adenosine monophosphate (cAMP) signalling and β-arrestin-2 translocation with temporal resolution (BRET-based biosensors). For the 3-carbamoyl-2-pyridone derivatives, we found that modifying the previously reported compound UOSS77 (also known as S-777469) by appending a PEG2-alcohol via a 3-carbomylcyclohexyl carboxamide (UOSS75) lowered lipophilicity, and preserved binding affinity and signalling profile. The 1,8-naphthyridin-2(1H)-one-3-carboxamide UOMM18, containing a cis configuration at the 3-carboxamide cyclohexyl and with an alcohol on the 4-position of the cyclohexyl, had lower lipophilicity but similar CB2 affinity and biological activity to previously reported compounds of this class. Relative to CP55,940, the new compounds acted as partial agonists and did not exhibit signalling bias. Interestingly, while all compounds shared similar temporal trajectories for maximal efficacy, differing temporal trajectories for potency were observed. Consequently, when applied at sub-maximal concentrations, CP55,940 tended to elicit sustained (cAMP) or increasing (arrestin) responses, whereas responses to the new compounds tended to be transient (cAMP) or sustained (arrestin). In future studies, the compounds characterised here may be useful in elucidating the consequences of differential temporal signalling profiles on CB2-mediated physiological responses.

Study of tyramine-binding mechanism and insecticidal activity of oil extracted from Eucalyptus against Sitophilus oryzae

The rice weevil, Sitophilus oryzae (L.), is a major pest of stored grains throughout the world, which causes quantitative and qualitative losses of food commodities. Eucalyptus essential oils (EOs) possess insecticidal and repellent properties, which make them a potential option for insect control in stored grains with environmentally friendly properties. In the current study, the binding mechanism of tyramine (TA) as a control compound has been investigated by funnel metadynamics (FM) simulation toward the homology model of tyramine1 receptor (TyrR) to explore its binding mode and key residues involved in the binding mechanism. EO compounds have been extracted from the leaf and flower part of Eucalyptus camaldulensis and characterized by GC/MS, and their effectiveness has been evaluated by molecular docking and conventional molecular dynamic (CMD) simulation toward the TyrR model. The FM results suggested that Asp114 followed by Asp80, Asn91, and Asn427 are crucial residues in the binding and the functioning of TA toward TyrR in Sitophilus Oryzae. The GC/MS analysis confirmed a total of 54 and 31 constituents in leaf and flower, respectively, where most of the components (29) are common in both groups. This analysis also revealed the significant concentration of Eucalyptus and α-pinene in leaves and flower EOs. The docking followed by CMD was performed to find the most effective compound in Eucalyptus EOs. In this regard, butanoic acid, 3-methyl-, 3-methyl butyl ester (B12) and 2-Octen-1-ol, 3,7-dimethyl- (B23) from leaf and trans- β-Ocimene (G04) from flower showed the maximum dock score and binding free energy, making them the leading candidates to replace tyramine in TyrR. The MM-PB/GBSA and MD analysis proved that the B12 structure is the most effective compound in inhibition of TyrR.

Mechanistic Origin of Partial Agonism of Tetrahydrocannabinol for Cannabinoid Receptors

Cannabinoid receptor 1 (CB₁) is a therapeutically relevant drug target for controlling pain, obesity, and other central nervous system disorders. However, full agonists and antagonists of CB₁ have been reported to cause serious side effects in patients. Therefore, partial agonists have emerged as a viable alternative as they can mitigate overstimulation and side effects. One of the key bottlenecks in the design of partial agonists, however, is the lack of understanding of the molecular mechanism of partial agonism itself. In this study, we examine two mechanistic hypotheses for the origin of partial agonism in cannabinoid receptors and predict the mechanistic basis of partial agonism exhibited by Δ⁹-Tetrahydrocannabinol (THC) against CB₁. In particular, we inspect whether partial agonism emerges from the ability of THC to bind in both agonist and antagonist-binding poses or from its ability to only partially activate the receptor. We used extensive molecular dynamics simulations and Markov state modeling to capture the THC binding in both antagonist and agonist-binding poses in the CB₁ receptor. Furthermore, we predict that binding of THC in the agonist-binding pose leads to rotation of toggle switch residues and causes partial outward movement of intracellular transmembrane helix 6 (TM6). Our simulations also suggest that the alkyl side chain of THC plays a crucial role in determining partial agonism by stabilizing the ligand in the agonist and antagonist-like poses within the pocket. Taken together, this study provides important insights into the mechanistic origin of the partial agonism of THC.

Discovery of Homobivalent Bitopic Ligands of the Cannabinoid CB2 Receptor*

Single chemical entities with potential to simultaneously interact with two binding sites are emerging strategies in medicinal chemistry. We have designed, synthesized and functionally characterized the first bitopic ligands for the CB2 receptor. These compounds selectively target CB2 versus CB1 receptors. Their binding mode was studied by molecular dynamic simulations and site-directed mutagenesis.

Agonist Binding and G Protein Coupling in Histamine H2 Receptor: A Molecular Dynamics Study

The histamine H2 receptor (H₂R) plays an important role in the regulation of gastric acid secretion. Therefore, it is a main drug target for the treatment of gastroesophageal reflux or peptic ulcer disease. However, there is as of yet no 3D-structural information available hampering a mechanistic understanding of H₂R. Therefore, we created a model of the histamine-H₂R-G_s complex based on the structure of the ternary complex of the β2-adrenoceptor and investigated the conformational stability of this active GPCR conformation. Since the physiologically relevant motions with respect to ligand binding and conformational changes of GPCRs can only partly be assessed on the timescale of conventional MD (cMD) simulations, we also applied metadynamics and Gaussian accelerated molecular dynamics (GaMD) simulations. A multiple walker metadynamics simulation in combination with cMD was applied for the determination of the histamine binding mode. The preferential binding pose detected is in good agreement with previous data from site directed mutagenesis and provides a basis for rational ligand design. Inspection of the H₂R-G_s interface reveals a network of polar interactions that may contribute to H₂R coupling selectivity. The cMD and GaMD simulations demonstrate that the active conformation is retained on a μs-timescale in the ternary histamine-H₂R-G_s complex and in a truncated complex that contains only G_s helix α5 instead of the entire G protein. In contrast, histamine alone is unable to stabilize the active conformation, which is in line with previous studies of other GPCRs.

Metadynamics simulations of ligand binding to GPCRs

Recent developments in metadynamics simulation techniques for ligand binding to Class A GPCRs are described and the results obtained elucidated. The computational protocol makes good use of modern massively parallel hardware, making simulations of the binding/unbinding process routine. The simulations reveal unprecedented details of the ligand-binding pathways, including multiple binding sites in many cases. Free energies of binding are reproduced very well and the simulations allow prediction of the efficacy (agonist, antagonist etc.) of ligands.

Translational potential of allosteric modulators targeting the cannabinoid CB1 receptor

The cannabinoid type-1 (CB₁) receptor, a G-protein-coupled receptor, is an attractive target for drug discovery due to its involvement in many physiological processes. Historically, drug discovery efforts targeting the CB₁ receptor have focused on the development of orthosteric ligands that interact with the active site to which endogenous cannabinoids bind. Research performed over the last several decades has revealed substantial difficulties in translating CB₁ orthosteric ligands into druggable candidates. The difficulty is mainly due to the adverse effects associated with orthosteric CB₁ ligands. Recent discoveries of allosteric CB₁ modulators provide tremendous opportunities to develop CB₁ ligands with novel mechanisms of action; these ligands may potentially improve the pharmacological effects and enhance drug safety in treating the disorders by regulating the functions of the CB₁ receptor. In this paper, we review and summarize the complex pharmacological profiles of each class of CB₁ allosteric modulators, the development of new classes of CB₁ allosteric modulators and the results from in vivo assessments of their therapeutic value.