OPM database and PPM web server: resources for positioning of proteins in membranes

Ligand recognition and G-protein coupling of trace amine receptor TAAR1

Trace-amine-associated receptors (TAARs), a group of biogenic amine receptors, have essential roles in neurological and metabolic homeostasis1. They recognize diverse endogenous trace amines and subsequently activate a range of G-protein-subtype signalling pathways2,3. Notably, TAAR1 has emerged as a promising therapeutic target for treating psychiatric disorders4,5. However, the molecular mechanisms underlying its ability to recognize different ligands remain largely unclear. Here we present nine cryo-electron microscopy structures, with eight showing human and mouse TAAR1 in a complex with an array of ligands, including the endogenous 3-iodothyronamine, two antipsychotic agents, the psychoactive drug amphetamine and two identified catecholamine agonists, and one showing 5-HT1AR in a complex with an antipsychotic agent. These structures reveal a rigid consensus binding motif in TAAR1 that binds to endogenous trace amine stimuli and two extended binding pockets that accommodate diverse chemotypes. Combined with mutational analysis, functional assays and molecular dynamic simulations, we elucidate the structural basis of drug polypharmacology and identify the species-specific differences between human and mouse TAAR1. Our study provides insights into the mechanism of ligand recognition and G-protein selectivity by TAAR1, which may help in the discovery of ligands or therapeutic strategies for neurological and metabolic disorders.

Structural basis of antibody inhibition and chemokine activation of the human CC chemokine receptor 8

The C-C motif chemokine receptor 8 (CCR8) is a class A G-protein coupled receptor that has emerged as a promising therapeutic target in cancer. Targeting CCR8 with an antibody has appeared to be an attractive therapeutic approach, but the molecular basis for chemokine-mediated activation and antibody-mediated inhibition of CCR8 are not fully elucidated. Here, we obtain an antagonist antibody against human CCR8 and determine structures of CCR8 in complex with either the antibody or the endogenous agonist ligand CCL1. Our studies reveal characteristic antibody features allowing recognition of the CCR8 extracellular loops and CCL1-CCR8 interaction modes that are distinct from other chemokine receptor - ligand pairs. Informed by these structural insights, we demonstrate that CCL1 follows a two-step, two-site binding sequence to CCR8 and that antibody-mediated inhibition of CCL1 signaling can occur by preventing the second binding event. Together, our results provide a detailed structural and mechanistic framework of CCR8 activation and inhibition that expands our molecular understanding of chemokine - receptor interactions and offers insight into the development of therapeutic antibodies targeting chemokine GPCRs.

An in silico investigation of Kv2.1 potassium channel: Model building and inhibitors binding sites analysis

Kv2.1 is widely expressed in brain, and inhibiting Kv2.1 is a potential strategy to prevent cell death and achieve neuroprotection in ischemic stroke. Herein, an in silico model of Kv2.1 tetramer structure was constructed by employing the AlphaFold-Multimer deep learning method to facilitate the rational discovery of Kv2.1 inhibitors. GaMD was utilized to create an ion transporting trajectory, which was analyzed with HMM to generate multiple representative receptor conformations. The binding site of RY785 and RY796(S) under the P-loop was defined with Fpocket program together with the competitive binding electrophysiology assay. The docking poses of the two inhibitors were predicted with the aid of the semi-empirical quantum mechanical calculation, and the IGMH results suggested that Met375, Thr376, and Thr377 of the P-helix and Ile405 of the S6 segment made significant contributions to the binding affinity. These results provided insights for rational molecular design to develop novel Kv2.1 inhibitors.

Deciphering specificity and cross-reactivity in tachykinin NK1 and NK2 receptors

The tachykinin receptors neurokinin 1 (NK1R) and neurokinin 2 (NK2R) are G protein-coupled receptors that bind preferentially to the natural peptide ligands substance P and neurokinin A, respectively, and have been targets for drug development. Despite sharing a common C-terminal sequence of Phe-X-Gly-Leu-Met-NH2 that helps direct biological function, the peptide ligands exhibit some degree of cross-reactivity toward each other's non-natural receptor. Here, we investigate the detailed structure-activity relationships of the ligand-bound receptor complexes that underlie both potent activation by the natural ligand and cross-reactivity. We find that the specificity and cross-reactivity of the peptide ligands can be explained by the interactions between the amino acids preceding the FxGLM consensus motif of the bound peptide ligand and two regions of the receptor: the β-hairpin of the extracellular loop 2 (ECL2) and a N-terminal segment leading into transmembrane helix 1. Positively charged sidechains of the ECL2 (R177 of NK1R and K180 of NK2R) are seen to play a vital role in the interaction. The N-terminal positions 1 to 3 of the peptide ligand are entirely dispensable. Mutated and chimeric receptor and ligand constructs neatly swap around ligand specificity as expected, validating the structure-activity hypotheses presented. These findings will help in developing improved agonists or antagonists for NK1R and NK2R.

Structural basis of G protein-Coupled receptor CMKLR1 activation and signaling induced by a chemerin-derived agonist

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Local Disordered Region Sampling (LDRS) for ensemble modeling of proteins with experimentally undetermined or low confidence prediction segments

SUMMARY

The Local Disordered Region Sampling (LDRS, pronounced loaders) tool is a new module developed for IDPConformerGenerator, a previously validated approach to model intrinsically disordered proteins (IDPs). The IDPConformerGenerator LDRS module provides a method for generating all-atom conformations of intrinsically disordered protein regions at N- and C-termini of and in loops or linkers between folded regions of an existing protein structure. These disordered elements often lead to missing coordinates in experimental structures or low confidence in predicted structures. Requiring only a pre-existing PDB or mmCIF formatted structural template of the protein with missing coordinates or with predicted confidence scores and its full-length primary sequence, LDRS will automatically generate physically meaningful conformational ensembles of the missing flexible regions to complete the full-length protein. The capabilities of the LDRS tool of IDPConformerGenerator include modeling phosphorylation sites using enhanced Monte Carlo-Side Chain Entropy, transmembrane proteins within an all-atom bilayer, and multi-chain complexes. The modeling capacity of LDRS capitalizes on the modularity, the ability to be used as a library and via command-line, and the computational speed of the IDPConformerGenerator platform.

AVAILABILITY AND IMPLEMENTATION

The LDRS module is part of the IDPConformerGenerator modeling suite, which can be downloaded from GitHub at https://github.com/julie-forman-kay-lab/IDPConformerGenerator. IDPConformerGenerator is written in Python3 and works on Linux, Microsoft Windows, and Mac OS versions that support DSSP. Users can utilize LDRS's Python API for scripting the same way they can use any part of IDPConformerGenerator's API, by importing functions from the "idpconfgen.ldrs_helper" library. Otherwise, LDRS can be used as a command line interface application within IDPConformerGenerator. Full documentation is available within the command-line interface as well as on IDPConformerGenerator's official documentation pages (https://idpconformergenerator.readthedocs.io/en/latest/).

Structural insights into ligand recognition and selectivity of the human hydroxycarboxylic acid receptor HCAR2

Hydroxycarboxylic acid receptor 2 (HCAR2) belongs to the family of class A G protein-coupled receptors with key roles in regulating lipolysis and free fatty acid formation in humans. It is deeply involved in many pathophysiological processes and serves as an attractive target for the treatment of cardiovascular, neoplastic, autoimmune, neurodegenerative, inflammatory, and metabolic diseases. Here, we report four cryo-EM structures of human HCAR2-Gi1 complexes with or without agonists, including the drugs niacin (2.69 Å) and acipimox (3.23 Å), the highly subtype-specific agonist MK-6892 (3.25 Å), and apo form (3.28 Å). Combined with molecular dynamics simulation and functional analysis, we have revealed the recognition mechanism of HCAR2 for different agonists and summarized the general pharmacophore features of HCAR2 agonists, which are based on three key residues R1113.36, S17945.52, and Y2847.43. Notably, the MK-6892-HCAR2 structure shows an extended binding pocket relative to other agonist-bound HCAR2 complexes. In addition, the key residues that determine the ligand selectivity between the HCAR2 and HCAR3 are also illuminated. Our findings provide structural insights into the ligand recognition, selectivity, activation, and G protein coupling mechanism of HCAR2, which shed light on the design of new HCAR2-targeting drugs for greater efficacy, higher selectivity, and fewer or no side effects.

A Combined Ligand- and Structure-Based Virtual Screening To Identify Novel NaV1.2 Blockers: In Vitro Patch Clamp Validation and In Vivo Anticonvulsant Activity

Epilepsy is a neurological disorder characterized by recurrent seizures that arise from abnormal electrical activity in the brain. Voltage-gated sodium channels (NaVs), responsible for the initiation and propagation of action potentials in neurons, play a critical role in the pathogenesis of epilepsy. This study sought to discover potential anticonvulsant compounds that interact with NaVs, specifically, the brain subtype hNaV1.2. A ligand-based QSAR model and a docking model were constructed, validated, and applied in a parallel virtual screening over the DrugBank database. Montelukast, Novobiocin, and Cinnarizine were selected for in vitro testing, using the patch-clamp technique, and all of them proved to inhibit hNaV1.2 channels heterologously expressed in HEK293 cells. Two hits were evaluated in the GASH/Sal model of audiogenic seizures and demonstrated promising activity, reducing the severity of sound-induced seizures at the doses tested. The combination of ligand- and structure-based models presents a valuable approach for identifying potential NaV inhibitors. These findings may provide a basis for further research into the development of new antiseizure drugs for the treatment of epilepsy.

Molecular Photoswitches Regulating the Activity of the Human Serotonin Transporter

Serotonin is an essential mediator regulating diverse neural processes, and its deregulation is related to the development of debilitating neurological diseases. In particular, the human serotonin transporter (hSERT) is fundamental in completing the synaptic neural cycle by allowing reuptake of serotonin. Its inhibition is particularly attractive, especially as a pharmacological target against depressive syndrome. Here, we analyze, by using long-range molecular dynamic simulations, the behavior of a molecular photoswitch whose cis- and trans-isomers inhibit the hSERT differently. In particular, we evidence the structural and molecular basis behind the higher inhibiting capacity of the cis-isomer, which blocks more efficiently the hSERT conformational cycle, leading to serotonin uptake.

Engineered odorant receptors illuminate structural principles of odor discrimination

A central challenge in olfaction is understanding how the olfactory system detects and distinguishes odorants with diverse physicochemical properties and molecular configurations. Vertebrate animals perceive odors via G protein-coupled odorant receptors (ORs). In humans, ~400 ORs enable the sense of smell. The OR family is composed of two major classes: Class I ORs are tuned to carboxylic acids while Class II ORs, representing the vast majority of the human repertoire, respond to a wide variety of odorants. How ORs recognize chemically diverse odorants remains poorly understood. A fundamental bottleneck is the inability to visualize odorant binding to ORs. Here, we uncover fundamental molecular properties of odorant-OR interactions by employing engineered ORs crafted using a consensus protein design strategy. Because such consensus ORs (consORs) are derived from the 17 major subfamilies of human ORs, they provide a template for modeling individual native ORs with high sequence and structural homology. The biochemical tractability of consORs enabled four cryoEM structures of distinct consORs with unique ligand recognition properties. The structure of a Class I consOR, consOR51, showed high structural similarity to the native human receptor OR51E2 and yielded a homology model of a related member of the human OR51 family with high predictive power. Structures of three Class II consORs revealed distinct modes of odorant-binding and activation mechanisms between Class I and Class II ORs. Thus, the structures of consORs lay the groundwork for understanding molecular recognition of odorants by the OR superfamily.

Analysis of structure-activity relationship of indol-3-yl-N-phenylcarbamic amides as potent STING inhibitors

A structure-activity relationship (SAR) study of stimulator of interferon gene (STING) inhibition was performed using a series of indol-3-yl-N-phenylcarbamic amides and indol-2-yl-N-phenylcarbamic amides. Among these analogs, compounds 10, 13, 15, 19, and 21 inhibited the phosphorylation of STING and interferon regulatory factor 3 (IRF3) to a greater extent than the reference compound, H-151. All five analogs showed stronger STING inhibition than H-151 on the 2',3'-cyclic GMP-AMP-induced expression of interferon regulatory factors (IRFs) in a STINGR232 knock-in THP-1 reporter cell line. The half-maximal inhibitory concentration of the most potent compound, 21, was 11.5 nM. The molecular docking analysis of compound 21 and STING combined with the SAR study suggested that the meta- and para-positions of the benzene ring of the phenylcarbamic amide moiety could be structurally modified by introducing halides or alkyl substituents.

Zika virus prM protein contains cholesterol binding motifs required for virus entry and assembly

For successful infection of host cells and virion production, enveloped viruses, including Zika virus (ZIKV), extensively rely on cellular lipids. However, how virus protein-lipid interactions contribute to the viral life cycle remains unclear. Here, we employ a chemo-proteomics approach with a bifunctional cholesterol probe and show that cholesterol is closely associated with the ZIKV structural protein prM. Bioinformatic analyses, reverse genetics alongside with photoaffinity labeling assays, and atomistic molecular dynamics simulations identified two functional cholesterol binding motifs within the prM transmembrane domain. Loss of prM-cholesterol association has a bipartite effect reducing ZIKV entry and leading to assembly defects. We propose a model in which membrane-resident M facilitates cholesterol-supported lipid exchange during endosomal entry and, together with cholesterol, creates a platform promoting virion assembly. In summary, we identify a bifunctional role of prM in the ZIKV life cycle by mediating viral entry and virus assembly in a cholesterol-dependent manner.

Automatic Rhodopsin Modeling with Multiple Protonation Microstates

Automatic Rhodopsin Modeling (ARM) is a simulation protocol providing QM/MM models of rhodopsins capable of reproducing experimental electronic absorption and emission trends. Currently, ARM is restricted to a single protonation microstate for each rhodopsin model. Herein, we incorporate an extension of the minimal electrostatic model (MEM) into the ARM protocol to account for all relevant protonation microstates at a given pH. The new ARM+MEM protocol determines the most important microstates contributing to the description of the absorption spectrum. As a test case, we have applied this methodology to simulate the pH-dependent absorption spectrum of a toy model, showing that the single-microstate picture breaks down at certain pH values. Subsequently, we applied ARM+MEM toAnabaenasensory rhodopsin, confirming an improved description of its absorption spectrum when the titration of several key residues is considered.

Molecular Mechanisms behind Conformational Transitions of the Influenza Virus Hemagglutinin Membrane Anchor

Membrane fusion is a fundamental process that is exploited by enveloped viruses to enter host cells. In the case of the influenza virus, fusion is facilitated by the trimeric viral hemagglutinin protein (HA). So far, major focus has been put on its N-terminal fusion peptides, which are directly responsible for fusion initiation. A growing body of evidence points also to a significant functional role of the HA C-terminal domain, which however remains incompletely understood. Our computational study aimed to elucidate the structural and functional interdependencies within the HA C-terminal region encompassing the transmembrane domain (TMD) and the cytoplasmic tail (CT). In particular, we were interested in the conformational shift of the TMD in response to varying cholesterol concentration in the viral membrane and in its modulation by the presence of CT. Using free-energy calculations based on atomistic molecular dynamics simulations, we characterized transitions between straight and tilted metastable TMD configurations under varying conditions. We found that the presence of CT is essential for achieving a stable, highly tilted TMD configuration. As we demonstrate, such a configuration of HA membrane anchor likely supports the tilting motion of its ectodomain, which needs to be executed during membrane fusion. This finding highlights the functional role of, so far, the relatively overlooked CT region.

Assembly of Protein Complexes In and On the Membrane with Predicted Spatial Arrangement Constraints

Membrane proteins play crucial roles in various cellular processes, and their interactions with other proteins in and on the membrane are essential for their proper functioning. While an increasing number of structures of more membrane proteins are being determined, the available structure data is still sparse. To gain insights into the mechanisms of membrane protein complexes, computational docking methods are necessary due to the challenge of experimental determination. Here, we introduce Mem-LZerD, a rigid-body membrane docking algorithm designed to take advantage of modern membrane modeling and protein docking techniques to facilitate the docking of membrane protein complexes. Mem-LZerD is based on the LZerD protein docking algorithm, which has been constantly among the top servers in many rounds of CAPRI protein docking assessment. By employing a combination of geometric hashing, newly constrained by the predicted membrane height and tilt angle, and model scoring accounting for the energy of membrane insertion, we demonstrate the capability of Mem-LZerD to model diverse membrane protein-protein complexes. Mem-LZerD successfully performed unbound docking on 13 of 21 (61.9%) transmembrane complexes in an established benchmark, more than shown by previous approaches. It was additionally tested on new datasets of 44 transmembrane complexes and 92 peripheral membrane protein complexes, of which it successfully modeled 35 (79.5%) and 15 (16.3%) complexes respectively. When non-blind orientations of peripheral targets were included, the number of successes increased to 54 (58.7%). We further demonstrate that Mem-LZerD produces complex models which are suitable for molecular dynamics simulation. Mem-LZerD is made available at https://lzerd.kiharalab.org.

Superiority of the Triple-Acting 5-HT6R/5-HT3R Antagonist and MAO-B Reversible Inhibitor over 5-HT6R Antagonist Intepirdine in Alleviation of Cognitive Deficits in Rats

The multifactorial origin and neurochemistry of Alzheimer's disease (AD) call for the development of multitarget treatment strategies. We report a first-in-class triple acting compound that targets serotonin type 6 and 3 receptors (5-HT-Rs) and monoamine oxidase type B (MAO-B) as an approach for treating AD. The key structural features required for MAO-B inhibition and 5-HT6R antagonism and interaction with 5-HT3R were determined using molecular dynamic simulations and cryo-electron microscopy, respectively. Bioavailable PZ-1922 reversed scopolamine-induced cognitive deficits in the novel object recognition test. Furthermore, it displayed superior pro-cognitive properties compared to intepirdine (a 5-HT6R antagonist) in the AD model, which involved intracerebroventricular injection of an oligomeric solution of amyloid-β peptide (oAβ) in the T-maze test in rats. PZ-1922, but not intepirdine, restored levels of biomarkers characteristic of the debilitating effects of oAβ. These data support the potential of a multitarget approach involving the joint modulation of 5-HT6R/5-HT3R/MAO-B in AD.

Niacinamide enhances cathelicidin mediated SARS-CoV-2 membrane disruption

The continual emergence of SARS-CoV-2 variants threatens to compromise the effectiveness of worldwide vaccination programs, and highlights the need for complementary strategies for a sustainable containment plan. An effective approach is to mobilize the body's own antimicrobial peptides (AMPs), to combat SARS-CoV-2 infection and propagation. We have found that human cathelicidin (LL37), an AMP found at epithelial barriers as well as in various bodily fluids, has the capacity to neutralise multiple strains of SARS-CoV-2. Biophysical and computational studies indicate that LL37's mechanism of action is through the disruption of the viral membrane. This antiviral activity of LL37 is enhanced by the hydrotropic action of niacinamide, which may increase the bioavailability of the AMP. Interestingly, we observed an inverse correlation between LL37 levels and disease severity of COVID-19 positive patients, suggesting enhancement of AMP response as a potential therapeutic avenue to mitigate disease severity. The combination of niacinamide and LL37 is a potent antiviral formulation that targets viral membranes of various variants and can be an effective strategy to overcome vaccine escape.

Synthesis, computational and experimental pharmacological studies for (thio)ether-triazine 5-HT6R ligands with noticeable action on AChE/BChE and chalcogen-dependent intrinsic activity in search for new class of drugs against Alzheimer's disease

Alzheimer's disease is becoming a growing problem increasing at a tremendous rate. Serotonin 5-HT6 receptors appear to be a particularly attractive target from a therapeutic perspective, due to their involvement not only in cognitive processes, but also in depression and psychosis. In this work, we present the synthesis and broad biological characterization of a new series of 18 compounds with a unique 1,3,5-triazine backbone, as potent 5-HT6 receptor ligands. The main aim of this research is to compare the biological activity of the newly synthesized sulfur derivatives with their oxygen analogues and their N-demethylated O- and S-metabolites obtained for the first time. Most of the new triazines displayed high affinity (Ki < 200 nM) and selectivity towards 5-HT6R, with respect to 5-HT2AR, 5-HT7R, and D2R, in the radioligand binding assays. For selected, active compounds crystallographic studies, functional bioassays, and ADME-Tox profile in vitro were performed. The exciting novelty is that the sulfur derivatives exhibit an agonistic mode of action contrary to all other compounds obtained to date in this chemical class herein and previously reported. Advanced computational studies indicated that this intriguing functional shift might be caused by presence of chalcogen bonds formed only by the sulfur atom. In addition, the N-demethylated derivatives have emerged highly potent antioxidants and, moreover, show a significant improvement in metabolic stability compared to the parent structures. The cholinesterase study present micromolar inhibitory AChE and BChE activity for both 5-HT6 agonist 19 and potent antagonist 5. Finally, the behavioral experiments of compound 19 demonstrated its antidepressant-like properties and slight ability to improve cognitive deficits, without inducing memory impairments by itself. Described pharmacological properties of both compounds (5 and 19) allow to give a design clue for the development of multitarget compounds with 5-HT6 (both agonist and antagonist)/AChE and/or BChE mechanism in the group of 1,3,5-triazine derivatives.

Structural dynamics of moonlighting intrinsically disordered proteins - A black box in multiple sclerosis

Multiple Sclerosis (MS) is a demyelinating disease of the central nervous system that disturbs the flow of brain signals to other parts of the body. The actual cause of the disease is still not apparent. The intrinsically disordered proteins (IDP) play a crucial role in neurodegenerative diseases like Alzheimer's, Lewy bodies, Parkinson's, Amyotrophic Lateral Sclerosis, Multiple Sclerosis, etc. In MS, it was known that the immune system attacks the proteins like Myelin Basic Protein (MBP), Myelin-associated Oligodendrocyte Basic protein (MOBP), Myelin-Associated Protein (MAG), and Myelin Proteolipid Protein (PLP) and this leads to demyelination causing MS. Here the proteins MBP and MOBP are both moonlighting intrinsically disordered proteins and exist between the myelin sheath, unlike MAG which is a transmembrane protein. The main focus of the article was to examine the significant role of proteins intrinsically disordered regions (IDR) in maintaining their function. Molecular dynamics simulation studies were performed to study the conformational dynamics of these protein IDRs both in water and near the myelin sheath. The results suggest that the IDR dominates the structural dynamics of these proteins and IDR in both proteins was responsible for their interaction with the myelin sheath. Interestingly, it was noted that the known epitopes MBP83-96 and MOBP65-87 in the IDR have no interaction with the myelin sheath. Thus when the protein remains intrinsically disordered it maintains the proper function and myelin integrity and if it adopts folds the region was identified as an epitope by the immune system leading to demyelination causing MS.

Membrane regulation of 15LOX-1/PEBP1 complex prompts the generation of ferroptotic signals, oxygenated PEs

Ferroptosis is a regulated form of cell death, the mechanism of which is still to be understood. 15-lipoxygenase (15LOX) complex with phosphatidylethanolamine (PE)-binding protein 1 (PEBP1) catalyzes the generation of pro-ferroptotic cell death signals, hydroperoxy-polyunsaturated PE. We focused on gaining new insights into the molecular basis of these pro-ferroptotic interactions using computational modeling and liquid chromatography-mass spectrometry experiments. Simulations of 15LOX-1/PEBP1 complex dynamics and interactions with lipids revealed that association with the membrane triggers a conformational change in the complex. This conformational change facilitates the access of stearoyl/arachidonoyl-PE (SAPE) substrates to the catalytic site. Furthermore, the binding of SAPE promotes tight interactions within the complex and induces further conformational changes that facilitate the oxidation reaction. The reaction yields two hydroperoxides as products, 15-HpETE-PE and 12-HpETE-PE, at a ratio of 5:1. A significant effect of PEBP1 is observed only on the predominant product. Moreover, combined experiments and simulations consistently demonstrate the significance of PEBP1 P112E mutation in generating ferroptotic cell death signals.

Computational and structural insights into the pre- and post-hydrolysis states of bovine multidrug resistance-associated protein 1

ATP-binding cassette C-family drug membrane transporters play an important role in local pharmacokinetics, that is, drug concentration in cellular compartments. From the structural point of view, only the bovine ortholog of the multidrug resistance-associated protein 1 (bMRP1) has been resolved. We here used μs-scaled molecular dynamics simulations to investigate the structure and dynamics of the bovine multidrug resistance-associated protein 1 in pre- and post-hydrolysis functional states. The present work aims to examine the slight but likely relevant structural differences between pre- and post-hydrolysis states of outward-facing conformations as well as the interactions between the multidrug resistance-associated protein 1 and the surrounding lipid bilayer. Global conformational dynamics show unfavourable extracellular opening associated with nucleotide-binding domain dimerization indicating that the post-hydrolysis state adopts a close-cleft conformation rather than an outward-open conformation. Our present simulations also highlight persistent interactions with annular cholesterol molecules and the expected active role of lipid bilayer in the allosteric communication between distant domains of the transporter.

Inhibition of altered Orai1 channels in Müller cells protects photoreceptors in retinal degeneration

The expressions of ion channels by Müller glial cells (MGCs) may change in response to various retinal pathophysiological conditions. There remains a gap in our understanding of MGCs' responses to photoreceptor degeneration towards finding therapies. The study explores how an inhibition of store-operated Ca2+ entry (SOCE) and its major component, Orai1 channel, in MGCs protects photoreceptors from degeneration. The study revealed increased Orai1 expression in the MGCs of retinal degeneration 10 (rd10) mice. Enhanced expression of oxidative stress markers was confirmed as a crucial pathological mechanism in rd10 retina. Inducing oxidative stress in rat MGCs resulted in increasing SOCE and Ca2+ release-activated Ca2+ (CRAC) currents. SOCE inhibition by 2-Aminoethoxydiphenyl borate (2-APB) protected photoreceptors in degenerated retinas. Finally, molecular simulations proved the structural and dynamical features of 2-APB to the target structure Orai1. Our results provide new insights into the physiology of MGCs regarding retinal degeneration and shed a light on SOCE and Orai1 as new therapeutic targets.

Utilizing Heteroatom Types and Numbers from Extensive Ligand Libraries to Develop Novel hERG Blocker QSAR Models Using Machine Learning-Based Classifiers

The human ether-à-go-go-related gene (hERG) channel plays a crucial role in membrane repolarization. Any disruptions in its function can lead to severe cardiovascular disorders such as long QT syndrome (LQTS), which increases the risk of serious cardiovascular problems such as tachyarrhythmia and sudden cardiac death. Drug-induced LQTS is a significant concern and has resulted in drug withdrawals from the market in the past. The main objective of this study is to pinpoint crucial heteroatoms present in ligands that initiate interactions leading to the effective blocking of the hERG channel. To achieve this aim, ligand-based quantitative structure-activity relationships (QSAR) models were constructed using extensive ligand libraries, considering the heteroatom types and numbers, and their associated hERG channel blockage pIC50 values. Machine learning-assisted QSAR models were developed to analyze the key structural components influencing compound activity. Among the various methods, the KPLS method proved to be the most efficient, allowing the construction of models based on eight distinct fingerprints. The study delved into investigating the influence of heteroatoms on the activity of hERG blockers, revealing their significant role. Furthermore, by quantifying the effect of heteroatom types and numbers on ligand activity at the hERG channel, six compound pairs were selected for molecular docking. Subsequent molecular dynamics simulations and per residue MM/GBSA calculations were performed to comprehensively analyze the interactions of the selected pair compounds.

Intracellular Plasmodium aquaporin 2 is important for sporozoite production in the mosquito vector and malaria transmission

Malaria remains a devastating disease and, with current measures failing to control its transmission, there is a need for novel interventions. A family of proteins that have long been pursued as potential intervention targets are aquaporins, which are channels facilitating the movement of water and other solutes across membranes. We identify an aquaporin in malaria parasites and demonstrate that it is important for completion of Plasmodium development in the mosquito vector. Disruption of AQP2 in the human parasite Plasmodium falciparum and the rodent parasite Plasmodium berghei blocks sporozoite production inside oocysts established on mosquito midguts, greatly limiting parasite infection of salivary glands and transmission to a new host. In vivo epitope tagging of AQP2 in P. berghei, combined with immunofluorescence assays, reveals that the protein is localized in vesicle-like organelles found in the cytoplasm of gametocytes, ookinetes, and sporozoites. The number of these organelles varies between individual parasites and lifecycle stages suggesting that they are likely part of a dynamic endomembrane system. Phylogenetic analysis confirms that AQP2 is unique to malaria and closely related parasites and most closely resembles intracellular aquaporins. Structure prediction analyses identify several unusual features, including a large accessory extracellular loop and an arginine-to-phenylalanine substitution in the selectivity filter principally determining pore function, a unique feature among known aquaporins. This in conjunction with the importance of AQP2 for malaria transmission suggests that AQP2 may be a fruitful target of antimalarial interventions.

Structure Elucidation and Interaction Dynamics of MefA-MsrD Efflux Proteins in Streptococcus pneumoniae: Impact on Macrolide Susceptibility

Macrolides are empirically used to treat bacterial community-acquired pneumonia (CAP). Streptococcus pneumoniae, being the major pathogen responsible for bacterial CAP with high mortality rates, express MefA-MsrD efflux pumps to hinder macrolide susceptibility. Despite its importance, the structural features of the efflux-protein complex and its impact on macrolide susceptibility have not yet been elucidated explicitly. Therefore, in the present study, combining homology, threading, and dynamics approaches, MefA and MsrD proteins in pathogenic S. pneumoniae were modeled. Both membrane (lipid-bilayer) and cytoplasmic (aqueous) environments were considered to simulate the MefA and MsrD proteins in their ideal cellular conditions followed by dynamics analyses. The simulated MefA structure represented a typical major facilitator superfamily protein structure with 13 transmembrane helices. MefA-MsrD interaction via clustering-based docking revealed low-energy conformers with stable intermolecular interactions. The higher clinical MIC value of azithromycin over erythromycin was reflected upon erythromycin eliciting stronger interactions (dissociation constant or ki = ∼52 μM) with the cytoplasmic ATP-binding MsrD than azithromycin (ki = ∼112 μM). The strong (binding energy = -132.1 ± 9.5 kcal/mol) and highly stable (root-mean-square fluctuation < 1.0 Å) physical association between MefA with MsrD was validated and was found to be unaffected by the antibiotic binding. Higher propensity of the macrolides to interact with MsrD than MefA established the importance of the former in macrolide susceptibility. Ours is probably the first report on the structural arrangements in the MefA-MsrD efflux complex and the macrolide susceptibility in S. pneumoniae. This study provides a novel lead for experimental explorations and efflux-pump inhibitor designs.

Structure and Function of ArnD. A Deformylase Essential for Lipid A Modification with 4-Amino-4-deoxy-l-arabinose and Polymyxin Resistance

Covalent modification of lipid A with 4-deoxy-4-amino-l-arabinose (Ara4N) mediates resistance to cationic antimicrobial peptides and polymyxin antibiotics in Gram-negative bacteria. The proteins required for Ara4N biosynthesis are encoded in the pmrE and arnBCADTEF loci, with ArnT ultimately transferring the amino sugar from undecaprenyl-phospho-4-deoxy-4-amino-l-arabinose (C55P-Ara4N) to lipid A. However, Ara4N is N-formylated prior to its transfer to undecaprenyl-phosphate by ArnC, requiring a deformylase activity downstream in the pathway to generate the final C55P-Ara4N donor. Here, we show that deletion of the arnD gene in an Escherichia coli mutant that constitutively expresses the arnBCADTEF operon leads to accumulation of the formylated ArnC product undecaprenyl-phospho-4-deoxy-4-formamido-l-arabinose (C55P-Ara4FN), suggesting that ArnD is the downstream deformylase. Purification of Salmonella typhimurium ArnD (stArnD) shows that it is membrane-associated. We present the crystal structure of stArnD revealing a NodB homology domain structure characteristic of the metal-dependent carbohydrate esterase family 4 (CE4). However, ArnD displays several distinct features: a 44 amino acid insertion, a C-terminal extension in the NodB fold, and sequence divergence in the five motifs that define the CE4 family, suggesting that ArnD represents a new family of carbohydrate esterases. The insertion is responsible for membrane association as its deletion results in a soluble ArnD variant. The active site retains a metal coordination H-H-D triad, and in the presence of Co2+ or Mn2+, purified stArnD efficiently deformylates C55P-Ara4FN confirming its role in Ara4N biosynthesis. Mutations D9N and H233Y completely inactivate stArnD implicating these two residues in a metal-assisted acid-base catalytic mechanism.

Structural basis for triacylglyceride extraction from mycobacterial inner membrane by MFS transporter Rv1410

Mycobacterium tuberculosis is protected from antibiotic therapy by a multi-layered hydrophobic cell envelope. Major facilitator superfamily (MFS) transporter Rv1410 and the periplasmic lipoprotein LprG are involved in transport of triacylglycerides (TAGs) that seal the mycomembrane. Here, we report a 2.7 Å structure of a mycobacterial Rv1410 homologue, which adopts an outward-facing conformation and exhibits unusual transmembrane helix 11 and 12 extensions that protrude ~20 Å into the periplasm. A small, very hydrophobic cavity suitable for lipid transport is constricted by a functionally important ion-lock likely involved in proton coupling. Combining mutational analyses and MD simulations, we propose that TAGs are extracted from the core of the inner membrane into the central cavity via lateral clefts present in the inward-facing conformation. The functional role of the periplasmic helix extensions is to channel the extracted TAG into the lipid binding pocket of LprG.

Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes

YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.

Phospholipid transport to the bacterial outer membrane through an envelope-spanning bridge

The outer membrane of Gram-negative bacteria provides a formidable barrier, essential for both pathogenesis and antimicrobial resistance. Biogenesis of the outer membrane requires the transport of phospholipids across the cell envelope. Recently, YhdP was implicated as a major protagonist in the transport of phospholipids from the inner membrane to the outer membrane however the molecular mechanism of YhdP mediated transport remains elusive. Here, utilising AlphaFold, we observe YhdP to form an elongated assembly of 60 β strands that curve to form a continuous hydrophobic groove. This architecture is consistent with our negative stain electron microscopy data which reveals YhdP to be approximately 250 Å in length and thus sufficient to span the bacterial cell envelope. Furthermore, molecular dynamics simulations and in vivo bacterial growth assays indicate essential helical regions at the N- and C-termini of YhdP, that may embed into the inner and outer membranes respectively, reinforcing its envelope spanning nature. Our in vivo crosslinking data reveal phosphate-containing substrates captured along the length of the YhdP groove, providing direct evidence that YhdP transports phospholipids. This finding is congruent with our molecular dynamics simulations which demonstrate the propensity for inner membrane lipids to spontaneously enter the groove of YhdP. Collectively, our results support a model in which YhdP bridges the cell envelope, providing a hydrophobic environment for the transport of phospholipids to the outer membrane.

Interaction of Synthetic Cannabinoid Receptor Agonists with Cannabinoid Receptor I: Insights into Activation Molecular Mechanism

Synthetic cannabinoid receptor agonists (SCRAs) have become a wide group of new psychoactive substances since the 2010s. For the last few years, the X-ray structures of the complexes of cannabinoid receptor I (CB1) with SCRAs as well as the complexes of CB1 with its antagonist have been published. Based on those data, SCRA-CB1 interactions are analyzed in detail, using molecular modeling and molecular dynamics simulations. The molecular mechanism of the conformational transformation of the transmembrane domain of CB1 caused by its interaction with SCRA is studied. These conformational changes allosterically modulate the CB1-Gi complex, providing activation of the Gi protein. Based on the X-ray-determined structures of the CB1-ligand complexes, a stable apo conformation of inactive CB1 with a relatively low potential barrier of receptor activation was modeled. For that model, molecular dynamic simulations of SCRA binding to CB1 led to the active state of CB1, which allowed us to explore the key features of this activation and the molecular mechanism of the receptor's structural transformation. The simulated CB1 activation is in accordance with the previously published experimental data for the activation at protein mutations or structural changes of ligands. The key feature of the suggested activation mechanism is the determination of the stiff core of the CB1 transmembrane domain and the statement that the entire conformational transformation of the receptor to the active state is caused by a shift of alpha helix TM7 relative to this core. The shift itself is caused by protein-ligand interactions. It was verified via steered molecular dynamics simulations of the X-ray-determined structures of the inactive receptor, which resulted in the active conformation of CB1 irrespective of the placement of agonist ligand in the receptor's active site.

TMH Stab-pred: Predicting the stability of α-helical membrane proteins using sequence and structural features

The folding and stability of transmembrane proteins (TMPs) are governed by the insertion of secondary structural elements into the cell membrane followed by their assembly. Understanding the important features that dictate the stability of TMPs is important for elucidating their functions. In this work, we related sequence and structure-based parameters with free energy (ΔG0) of α-helical membrane proteins. Our results showed that the free energy transfer of hydrophobic peptides, relative contact order, total interaction energy, number of hydrogen bonds and lipid accessibility of transmembrane regions are important for stability. Further, we have developed multiple-regression models to predict the stability of α-helical membrane proteins using these features and our method can predict the stability with a correlation and mean absolute error (MAE) of 0.89 and 1.21 kcal/mol, respectively, on jack-knife test. The method was validated with a blind test set of three recently reported experimental ΔG0, which could predict the stability within an average MAE of 0.51 kcal/mol. Further, we developed a webserver for predicting the stability and it is freely available at (https://web.iitm.ac.in/bioinfo2/TMHS/). The importance of selected parameters and limitations are discussed.

Structural diversity of leukotriene G-protein coupled receptors

Dihydroxy acid leukotriene (LTB4) and cysteinyl leukotrienes (LTC4, LTD4, and LTE4) are inflammatory mediators derived from arachidonic acid via the 5-lipoxygenase pathway. While structurally similar, these two types of leukotrienes (LTs) exert their functions through interactions with two distinct G protein-coupled receptor (GPCR) families, BLT and CysLT receptors, which share low sequence similarity and belong to phylogenetically divergent GPCR groups. Selective antagonism of LT receptors has been proposed as a promising strategy for the treatment of many inflammation-related diseases including asthma and chronic obstructive pulmonary disease, rheumatoid arthritis, cystic fibrosis, diabetes, and several types of cancer. Selective CysLT1R antagonists are currently used as antiasthmatic drugs, however, there are no approved drugs targeting CysLT2 and BLT receptors. In this review, we highlight recently published structures of BLT1R and CysLTRs revealing unique structural features of the two receptor families. X-ray and cryo-EM data shed light on their overall conformations, differences in functional motifs involved in receptor activation, and details of the ligand-binding pockets. An unexpected binding mode of the selective antagonist BIIL260 in the BLT1R structure makes it the first example of a compound targeting the sodium-binding site of GPCRs and suggests a novel strategy for the receptor activity modulation. Taken together, these recent structural data reveal dramatic differences in the molecular architecture of the two LT receptor families and pave the way to new therapeutic strategies of selective targeting individual receptors with novel tool compounds obtained by the structure-based drug design approach.

Hesperidin potentially interacts with the catalytic site of gamma-secretase and modifies notch sensitive genes and cancer stemness marker expression in colon cancer cells and colonosphere

Gamma secretase (GS) produces Notch Intracellular Domain (NICD) by trans-membrane cleavage of notch receptor. The NICD enters the nucleus and activates the notch signaling pathway (NSP) by activating notch-responsive gene transcription. Hyperactivation of NSP is related to cancer aggressiveness, therapy resistance, and poor therapy outcome, and decreased overall disease-free survival in patients. Till date, none of the GS inhibitors (GSI) has been clinically approved due to their toxicity in patients. Thus in the present study, we explored the GS catalytic site binding potential of hesperidin (natural flavone glycoside) and its effect on notch responsive gene expression in HCT-116 cells. Molecular docking, MM-GBSA binding energy calculations, and molecular dynamics (MD) simulation experiments were performed to study the GS catalytic site binding potential of hesperidin. The compound showed better GS catalytic site binding potential at the active site compared to experimentally validated GSI, N-N-(3, 5-Difluorophenacetyl)-L-alanyl-S-phenylglycine t-butyl ester (DAPT) in molecular docking and MM-GBSA experiments. MD simulation results showed that hesperidin forms stable and energetically favorable complex with gamma secretase in comparison to standard inhibitor (DAPT)-GS complex. Further, in vitro experiments showed that hesperidin inhibited cell growth and sphere formation potential in HCT-116 cells. Further, hesperidin treatment altered notch responsive genes (Hes1, Hey1, and E-cad) and cancer stemness/self-renewal markers expression at transcription levels. In conclusion, hesperidin produces toxicity in HCT-116 cells and decreases colonosphere formation by inhibiting transcription of notch signaling pathway target genes and stemness markers.Communicated by Ramaswamy H. Sarma.

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.

Sigma Receptors: Novel Regulators of Iron/Heme Homeostasis and Ferroptosis

Sigma receptors are non-opiate/non-phencyclidine receptors that bind progesterone and/or heme and also several unrelated xenobiotics/chemicals. They reside in the plasma membrane and in the membranes of the endoplasmic reticulum, mitochondria, and nucleus. Until recently, the biology/pharmacology of these proteins focused primarily on their role in neuronal functions in the brain/retina. However, there have been recent developments in the field with the discovery of unexpected roles for these proteins in iron/heme homeostasis. Sigma receptor 1 (S1R) regulates the oxidative stress-related transcription factor NRF2 and protects against ferroptosis, an iron-induced cell death process. Sigma receptor 2 (S2R), which is structurally unrelated to S1R, complexes with progesterone receptor membrane components PGRMC1 and PGRMC2. S2R, PGRMC1, and PGRMC2, either independently or as protein-protein complexes, elicit a multitude of effects with a profound influence on iron/heme homeostasis. This includes the regulation of the secretion of the iron-regulatory hormone hepcidin, the modulation of the activity of mitochondrial ferrochelatase, which catalyzes iron incorporation into protoporphyrin IX to form heme, chaperoning heme to specific hemoproteins thereby influencing their biological activity and stability, and protection against ferroptosis. Consequently, S1R, S2R, PGRMC1, and PGRMC2 potentiate disease progression in hemochromatosis and cancer. These new discoveries usher this intriguing group of non-traditional progesterone receptors into an unchartered territory in biology and medicine.

Rational design and characterization of cell-selective antimicrobial peptides based on a bioactive peptide from Crocodylus siamensis hemoglobin

Antimicrobial resistance is a growing health concern. Antimicrobial peptides are a potential solution because they bypass conventional drug resistance mechanisms. Previously, we isolated a peptide from Crocodylus siamensis hemoglobin hydrolysate, which has antimicrobial activity and identified the main peptide from this mixture (QL17). The objective of this work was to evaluate and rationally modify QL17 in order to: (1) control its mechanism of action through bacterial membrane disruption; (2) improve its antimicrobial activity; and (3) ensure it has low cytotoxicity against normal eukaryotic cells. QL17 was rationally designed using physicochemical and template-based methods. These new peptide variants were assessed for: (1) their in vitro inhibition of microbial growth, (2) their cytotoxicity against normal cells, (3) their selectivity for microbes, and (4) the mode of action against bacteria using scanning electron microscopy (SEM), transmission electron microscopy (TEM) and confocal microscopy. The results indicate that all designed peptides have more potent antimicrobial efficacy than QL17 and IL15 peptides. However, only the most rationally modified peptides showed strong antimicrobial activity and minimal toxicity against normal cells. In particular, IL15.3 (hydrophobicity of 47% and net charge of + 6) was a potent antimicrobial agent (MIC = 4-12 μg/mL; MBC = 6-25 μg/mL) and displayed excellent selectivity for microbes (cf. human cells) via FACS assays. Microscopy confirmed that IL15.3 acts against bacteria by disrupting the cell membrane integrity and penetrating into the membrane. This causes the release of intracellular content into the outer environment leading to the death of bacteria. Moreover, IL15.3 can also interact with DNA suggesting it could have dual mode of action. Overall, a novel variant of QL17 is described that increases antimicrobial activity by over 1000-fold (~ 5 μg/mL MIC) and has minimal cytotoxicity. It may have applications in clinical use to treat and safeguard against bacteria.

Orthosteric ligand selectivity and allosteric probe dependence at Hydroxycarboxylic acid receptor HCAR2

Hydroxycarboxylic acid receptor 2 (HCAR2), a member of Class A G-protein-coupled receptor (GPCR) family, plays a pivotal role in anti-lipolytic and anti-inflammatory effects, establishing it as a significant therapeutic target for treating dyslipidemia and inflammatory diseases. However, the mechanism underlying the signaling of HCAR2 induced by various types of ligands remains elusive. In this study, we elucidate the cryo-electron microscopy (cryo-EM) structure of Gi-coupled HCAR2 in complex with a selective agonist, MK-6892, resolved to a resolution of 2.60 Å. Our structural analysis reveals that MK-6892 occupies not only the orthosteric binding pocket (OBP) but also an extended binding pocket (EBP) within HCAR2. Pharmacological assays conducted in this study demonstrate that the OBP is a critical determinant for ligand selectivity among the HCARs subfamily. Moreover, we investigate the pharmacological properties of the allosteric modulator compound 9n, revealing its probe-dependent behavior on HCAR2 in response to varying orthosteric agonists. Collectively, our findings provide invaluable structural insights that contribute to a deeper understanding of the regulatory mechanisms governing HCAR2 signaling transduction mediated by both orthosteric and allosteric ligands.

The Leu/Val6.51 Side Chain of Cannabinoid Receptors Regulates the Binding Mode of the Alkyl Chain of Δ9-Tetrahydrocannabinol

(-)-Δ9-trans-tetrahydrocannabinol (THC), which is the principal psychoactive constituent of Cannabis, mediates its action by binding to two members of the G-protein-coupled receptor (GPCR) family: the cannabinoid CB1 (CB1R) and CB2 (CB2R) receptors. Molecular dynamics simulations showed that the pentyl chain of THC could adopts an I-shape conformation, filling an intracellular cavity between Phe3.36 and Trp6.48 for initial agonist-induced receptor activation, in CB1R but not in CB2R. This cavity opens to the five-carbon chain of THC by the conformational change of the γ-branched, flexible, Leu6.51 side chain of CB1R, which is not feasible by the β-branched, mode rigid, Val6.51 side chain of CB2R. In agreement with our computational results, THC could not decrease the forskolin-induced cAMP levels in cells expressing mutant CB1RL6.51V receptor but could activate the mutant CB2RV6.51L receptor as efficiently as wild-type CB1R. Additionally, JWH-133, a full CB2R agonist, contains a branched dimethyl moiety in the ligand chain that bridges Phe3.36 and Val6.51 for receptor activation. In this case, the substitution of Val6.51 to Leu in CB2R makes JWH-133 unable to activate CB2RV6.51L. In conclusion, our combined computational and experimental results have shown that the amino acid at position 6.51 is a key additional player in the initial mechanism of activation of GPCRs that recognize signaling molecules derived from lipid species.

Mapping Structural Distribution and Gating-Property Impacts of Disease-Associated Missense Mutations in Voltage-Gated Sodium Channels

Thousands of voltage-gated sodium (Nav) channel variants contribute to a variety of disorders, including epilepsy, autism, cardiac arrhythmia, and pain disorders. Yet variant effects of more mutations remain unclear. The conventional gain-of-function (GoF) or loss-of-function (LoF) classifications is frequently employed to interpret of variant effects on function and guide precision therapy for sodium channelopathies. Our study challenges this binary classification by analyzing 525 mutations associated with 34 diseases across 366 electrophysiology studies, revealing that diseases with similar phenotypic effects can stem from unique molecular mechanisms. Our results show a high biophysical agreement (86%) between homologous disease-associated variants in different Nav genes, significantly surpassing the 60% phenotype (GoFo/LoFo) agreement among homologous mutants, suggesting the need for more nuanced disease categorization and treatment based on specific gating-property changes. Using UniProt data, we mapped over 2,400 disease-associated missense variants across nine human Nav channels and identified three clusters of mutation hotspots. Our findings indicate that mutations near the selectivity filter generally diminish the maximal current amplitude, while those in the fast inactivation region lean towards a depolarizing shift in half-inactivation voltage in steady-state activation, and mutations in the activation gate commonly enhance persistent current. In contrast to mutations in the PD, those within the VSD exhibit diverse impacts and subtle preferences on channel activity. This study shows great potential to enhance prediction accuracy for variant effects based on the structural context, laying the groundwork for targeted drug design in precision medicine.

Structural basis of hydroxycarboxylic acid receptor signaling mechanisms through ligand binding

Hydroxycarboxylic acid receptors (HCA) are expressed in various tissues and immune cells. HCA2 and its agonist are thus important targets for treating inflammatory and metabolic disorders. Only limited information is available, however, on the active-state binding of HCAs with agonists. Here, we present cryo-EM structures of human HCA2-Gi and HCA3-Gi signaling complexes binding with multiple compounds bound. Agonists were revealed to form a salt bridge with arginine, which is conserved in the HCA family, to activate these receptors. Extracellular regions of the receptors form a lid-like structure that covers the ligand-binding pocket. Although transmembrane (TM) 6 in HCAs undergoes dynamic conformational changes, ligands do not directly interact with amino acids in TM6, suggesting that indirect signaling induces a slight shift in TM6 to activate Gi proteins. Structural analyses of agonist-bound HCA2 and HCA3 together with mutagenesis and molecular dynamics simulation provide molecular insights into HCA ligand recognition and activation mechanisms.

Structural insights into FSP1 catalysis and ferroptosis inhibition

Ferroptosis suppressor protein 1 (FSP1, also known as AIMF2, AMID or PRG3) is a recently identified glutathione-independent ferroptosis suppressor1-3, but its underlying structural mechanism remains unknown. Here we report the crystal structures of Gallus gallus FSP1 in its substrate-free and ubiquinone-bound forms. The structures reveal a FAD-binding domain and a NAD(P)H-binding domain, both of which are shared with AIF and NADH oxidoreductases4-9, and a characteristic carboxy-terminal domain as well. We demonstrate that the carboxy-terminal domain is crucial for the catalytic activity and ferroptosis inhibition of FSP1 by mediating the functional dimerization of FSP1, and the formation of two active sites located on two sides of FAD, which are responsible for ubiquinone reduction and a unique FAD hydroxylation respectively. We also identify that FSP1 can catalyze the production of H2O2 and the conversion of FAD to 6-hydroxy-FAD in the presence of oxygen and NAD(P)H in vitro, and 6-hydroxy-FAD directly inhibits ferroptosis in cells. Together, these findings further our understanding on the catalytic and ferroptosis suppression mechanisms of FSP1 and establish 6-hydroxy-FAD as an active cofactor in FSP1 and a potent radical-trapping antioxidant in ferroptosis inhibition.

Chronic cough relief by allosteric modulation of P2X3 without taste disturbance

P2X receptors are cation channels that sense extracellular ATP. Many therapeutic candidates targeting P2X receptors have begun clinical trials or acquired approval for the treatment of refractory chronic cough (RCC) and other disorders. However, the present negative allosteric modulation of P2X receptors is primarily limited to the central pocket or the site below the left flipper domain. Here, we uncover a mechanism of allosteric regulation of P2X3 in the inner pocket of the head domain (IP-HD), and show that the antitussive effects of quercetin and PSFL2915 (our nM-affinity P2X3 inhibitor optimized based on quercetin) on male mice and guinea pigs were achieved by preventing allosteric changes of IP-HD in P2X3. While being therapeutically comparable to the newly licensed P2X3 RCC drug gefapixant, quercetin and PSFL2915 do not have an adverse effect on taste as gefapixant does. Thus, allosteric modulation of P2X3 via IP-HD may be a druggable strategy to alleviate RCC.

Discovery of WLB-89462, a New Drug-like and Highly Selective σ2 Receptor Ligand with Neuroprotective Properties

The synthesis and pharmacological activity of a new series of isoxazolylpyrimidines as sigma-2 receptor (σ2R) ligands are reported. Modification of a new hit retrieved in an HTS campaign allowed the identification of the compound WLB-89462 (20c) with good σ2R affinity (Ki = 13 nM) and high selectivity vs both the σ1R (Ki = 1777 nM) and a general panel of 180 targets. It represents one of the first σ2R ligands with drug-like properties, linked to a good physicochemical and ADMET profile (good solubility, no CYP inhibition, good metabolic stability, high permeability, brain penetration, and high oral exposure in rodents). Compound 20c shows neuroprotective activity in vitro and improves short-term memory impairment induced by hippocampal injection of amyloid β peptide in rats. Together with the promising effects in the chronic models where 20c is currently being evaluated, these results pave the way toward its clinical development as a neuroprotective agent.

First Potent Macrocyclic A3 Adenosine Receptor Agonists Reveal G-Protein and β-Arrestin2 Signaling Preferences

(N)-Methanocarba adenosine derivatives (A3 adenosine receptor (AR) agonists containing bicyclo[3.1.0]hexane replacing furanose) were chain-extended at N6 and C2 positions with terminal alkenes for ring closure. The resulting macrocycles of 17-20 atoms retained affinity, indicating a spatially proximal orientation of these receptor-bound chains, consistent with molecular modeling of 12. C2-Arylethynyl-linked macrocycle 19 was more A3AR-selective than 2-ether-linked macrocycle 12 (both 5'-methylamides, human (h) A3AR affinities (Ki): 22.1 and 25.8 nM, respectively), with lower mouse A3AR affinities. Functional hA3AR comparison of two sets of open/closed analogues in β-arrestin2 and Gi/o protein assays showed certain signaling preferences divergent from reference agonist Cl-IB-MECA 1. The potencies of 1 at all three Gαi isoforms were slightly less than its hA3AR binding affinity (Ki: 1.4 nM), while the Gαi1 and Gαi2 potencies of macrocycle 12 were roughly an order of magnitude higher than its radioligand binding affinity. Gαi2-coupling was enhanced in macrocycle 12 (EC50 2.56 nM, ∼40% greater maximal efficacy than 1). Di-O-allyl precursor 18 cyclized to form 19, increasing the Gαi1 potency by 7.5-fold. The macrocycles 12 and 19 and their open precursors 11 and 18 potently stimulated β-arrestin2 recruitment, with EC50 values (nM) of 5.17, 4.36, 1.30, and 4.35, respectively, and with nearly 50% greater efficacy compared to 1. This example of macrocyclization altering the coupling pathways of small-molecule (nonpeptide) GPCR agonists is the first for potent and selective macrocyclic AR agonists. These initial macrocyclic derivatives can serve as a guide for the future design of macrocyclic AR agonists displaying unanticipated pharmacology.

Biased allosteric activation of ketone body receptor HCAR2 suppresses inflammation

Hydroxycarboxylic acid receptor 2 (HCAR2), modulated by endogenous ketone body β-hydroxybutyrate and exogenous niacin, is a promising therapeutic target for inflammation-related diseases. HCAR2 mediates distinct pathophysiological events by activating Gi/o protein or β-arrestin effectors. Here, we characterize compound 9n as a Gi-biased allosteric modulator (BAM) of HCAR2 and exhibit anti-inflammatory efficacy in RAW264.7 macrophages via a specific HCAR2-Gi pathway. Furthermore, four structures of HCAR2-Gi complex bound to orthosteric agonists (niacin or monomethyl fumarate), compound 9n, and niacin together with compound 9n simultaneously reveal a common orthosteric site and a unique allosteric site. Combined with functional studies, we decipher the action framework of biased allosteric modulation of compound 9n on the orthosteric site. Moreover, co-administration of compound 9n with orthosteric agonists could enhance anti-inflammatory effects in the mouse model of colitis. Together, our study provides insight to understand the molecular pharmacology of the BAM and facilitates exploring the therapeutic potential of the BAM with orthosteric drugs.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the β-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or β-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

Interactions of the male contraceptive target EPPIN with semenogelin-1 and small organic ligands

Novel male contraceptives will promote gender equality in sharing contraceptive responsibility. The sperm-associated protein epididymal protease inhibitor (EPPIN) is a promising target for non-hormonal male contraception. EPPIN interacts with the semen coagulum protein semenogelin-1 (SEMG1) on the sperm surface, leading to transient inhibition of sperm motility after ejaculation. Small organic molecules targeting EPPIN's SEMG1-binding are under development as male contraceptives. Here, we combined computational approaches to uncover key aspects underlying EPPIN binding to SEMG1 and small organic ligands. We generated a human EPPIN model showing a typical arrangement of the WFDC (Whey-acid four disulfide core)-type and Kunitz-type domains, connected by a hinge region. Determining the EPPIN model's intrinsic motion by molecular dynamics simulations and normal mode analysis revealed a conformation, presenting a binding pocket that accommodates SEMG1Glu229-Gln247, EP055, and EP012. EPPIN's residues Phe63 and Lys68 (WFDC domain), Asp71 (hinge region), and Asn113, Asn114, and Asn115 (Kunitz domain) were identified as hot spots for SEMG1, EP055, and EP012 binding. Moreover, hydrophobic and hydrophilic residues in the WFDC and Kunitz domains allow plasma membrane anchoring, orienting the EPPIN binding pocket to the solvent. Targeting EPPIN's essential residues for its biomolecular interactions may improve the rational design of EPPIN ligands as spermiostatic compounds.

Signal peptide mimicry primes Sec61 for client-selective inhibition

Preventing the biogenesis of disease-relevant proteins is an attractive therapeutic strategy, but attempts to target essential protein biogenesis factors have been hampered by excessive toxicity. Here we describe KZR-8445, a cyclic depsipeptide that targets the Sec61 translocon and selectively disrupts secretory and membrane protein biogenesis in a signal peptide-dependent manner. KZR-8445 potently inhibits the secretion of pro-inflammatory cytokines in primary immune cells and is highly efficacious in a mouse model of rheumatoid arthritis. A cryogenic electron microscopy structure reveals that KZR-8445 occupies the fully opened Se61 lateral gate and blocks access to the lumenal plug domain. KZR-8445 binding stabilizes the lateral gate helices in a manner that traps select signal peptides in the Sec61 channel and prevents their movement into the lipid bilayer. Our results establish a framework for the structure-guided discovery of novel therapeutics that selectively modulate Sec61-mediated protein biogenesis.

Investigating membrane-binding properties of lipoxygenases using surface plasmon resonance

Lipoxygenases (LOXs) catalyze the oxidation of polyunsaturated fatty acids and synthesize oxylipin products that drive important cellular signaling processes in plants and animals. While there has been indirect evidence presented for the interaction of mammalian LOXs with membranes, a quantitative study of the molecular details of LOX-membrane interactions is lacking. Here, we mimicked biological membranes using surface plasmon resonance (SPR) sensor chips derivatized with 2-D planar lipophilic anchors (2D LP) to capture liposomes of varying phospholipid compositions that self-assemble into lipid bilayers on the SPR chip. The sensor chip surfaces were then used to investigate the membrane-binding properties of model LOX enzymes. SPR binding assays displayed reproducible and stable liposome capture to the sensor chip surface that allowed for the detailed characterization of LOX-membrane interactions. Our studies demonstrate a calcium-dependence for the membrane binding activities of coral 8R-LOX and human 15-LOX-2. Furthermore, our data confirm the importance of key membrane insertion loop residues in each of these LOX enzymes for membrane binding activity. Experiments utilizing model plant and human LOXs reveal differences in membrane-binding specificities. Our study establishes and validates a robust SPR-based platform using 2D LP sensor chips that allows for the detailed study of LOX-membrane interactions under different experimental conditions, including altered membrane compositions. Collectively, this investigation improves our overall understanding of LOX-membrane interaction properties, and our SPR-based approach holds potential for future use in the development of LOX-based therapeutics.

Dynamics of the Apo µ-Opioid Receptor in Complex with Gi Protein

Opioid receptors, particularly the µ-opioid receptor (μOR), play a pivotal role in mediating the analgesic and addictive effects of opioid drugs. G protein signaling is an important pathway of μOR function, usually associated with painkilling effects. However, the molecular mechanisms underlying the interaction between the μOR and G protein remain poorly understood. In this study, we employed classical all-atom molecular dynamics simulations to investigate the structural changes occurring with the μOR-G protein complex under two different conditions: with the G protein in the apo form (open) and with the GDP bound G protein (closed, holo form). The receptor was in the apo form and active conformation in both cases, and the simulation time comprised 1µs for each system. In order to assess the effect of the G protein coupling on the receptor activation state, three parameters were monitored: the correlation of the distance between TM3 and TM6 and the RMSD of the NPxxYA motif; the universal activation index (A100); and the χ2 dihedral distribution of residue W2936.48. When complexed with the open G protein, receptor conformations with intermediate activation state prevailed throughout the molecular dynamics, whereas in the condition with the closed G protein, mostly inactive conformations of the receptor were observed. The major effect of the G protein in the receptor conformation comes from a steric hindrance involving an intracellular loop of the receptor and a β-sheet region of the G protein. This suggests that G-protein precoupling is essential for receptor activation, but this fact is not sufficient for complete receptor activation.