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

Halogen Bonding Hot Spots as a Constraint in Virtual Screening: A Case Study of 5-HT7R

The recently developed and used molecular modeling approach to search for privileged amino acids for halogen bonding (XB hot spots) through XSAR sets has been applied to 5-HT7R. Herein, among all identified 5-HT7R XB hot spots, the S5x42 was employed in a virtual screening protocol as a constraint. Through a designed virtual screening protocol, 63 XSAR sets (156 compounds) were selected from more than 8 million commercially available compounds and examined using in vitro assay toward 5-HT7R. A 68% accuracy was found in predicting halogenated derivatives with higher affinity for 5-HT7R than their unsubstituted analogs. Moreover, it was observed that a halogen bond formed between S5x42 and a chlorine atom at the 3-position of the arylpiperazine fragment caused the most remarkable, 35.4-fold increase in binding affinity for 5-HT7R when compared to the nonhalogenated analog. Interestingly, molecular dynamics simulations showed the formation of a bifurcated halogen bond with S5x42.

Unlocking the Reverse Targeting Mechanisms of Cannabidiol: Unveiling New Therapeutic Avenues

Cannabidiol (CBD) and Δ9-tetrahydrocannabinol (THC), the main components of Cannabis sativa plants, have attracted a significant amount of attention due to their biological activities. This study identified GPR18 as the target of partial agonist CBD activating the p42/p44 MAPK pathway leading to migration of endometrial epithelial cells. Induced fit docking (IFD) showed that the affinity of THC for GPR18 is higher than that of CBD, and molecular dynamics (MD) simulations showed that CBD-GPR18 complexes at 130/200 ns might have stable conformations, potentially activating GPR18 by changing the distances of key residues in its active pocket. In contrast, THC maintains "metastable" conformations, generating a "shrinking space" leading to full agonism of THC by adding mechanical constraints in GPR18's active pocket. Steered molecular dynamics (SMD) revealed GPR18's active pocket was influenced more by CBD's partial agonism compared with THC. This combined IFD-MD-SMD method may be used to explain the mechanism of activation of partial or full agonists of GPR18.

Synthesis and Biological Evaluation of Pyrazole-Pyrimidones as a New Class of Correctors of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR)

Cystic fibrosis (CF) is caused by the functional expression defect of the cystic fibrosis transmembrane conductance regulator (CFTR) protein. Despite the recent success in CFTR modulator development, the available correctors only partially restore the F508del-CFTR channel function, and several rare CF mutations show resistance to available drugs. We previously identified compound 4172 that synergistically rescued the F508del-CFTR folding defect in combination with the existing corrector drugs VX-809 and VX-661. Here, novel CFTR correctors were designed by applying a classical medicinal chemistry approach on the 4172 scaffold. Molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted to propose a plausible binding site and design more potent and effective analogs. We identified three optimized compounds, which, in combination with VX-809 and the investigational corrector 3151, increased the plasma membrane density and function of F508del-CFTR and other rare CFTR mutants resistant to the currently approved therapies.

Discovery of Novel, Selective, and Nonbasic Agonists for the Kappa-Opioid Receptor Determined by Salvinorin A-Based Virtual Screening

Modulating the kappa-opioid receptor (KOR) is a promising strategy for treating various human diseases. KOR agonists show potential for treating pain, pruritus, and epilepsy, while KOR antagonists show potential for treating depression, anxiety, and addiction. The diterpenoid Salvinorin A (SalA), a secondary metabolite of Salvia divinorum, is a potent and selective KOR agonist. Unlike typical opioids, SalA lacks a basic nitrogen, which encouraged us to search for nonbasic KOR ligands. Through structure-based virtual screening using 3D pharmacophore models based on the binding mode of SalA, we identified novel, nonbasic, potent, and selective KOR agonists. In vitro studies confirmed two virtual hits, SalA-VS-07 and SalA-VS-08, as highly selective for the KOR and showing G protein-biased KOR agonist activity. Both KOR ligands share a novel spiro-moiety and a nonbasic scaffold. Our findings provide novel starting points for developing therapeutics aimed at treating pain and other conditions in which KOR is a central player.

Elucidating the Activation Mechanism of the Proton-sensing GPR68 Receptor

GPR68 is a proton-sensing G-protein Coupled Receptor (GPCR) involved in a variety of physiological processes and disorders including neoplastic pathologies. While GPR68 and few other GPCRs have been shown to be activated by a decrease in the extracellular pH, the molecular mechanism of their activation remains largely unknown. In this work, we used a combined computational and in vitro approach to provide new insight into the activation mechanism of the receptor. Molecular Dynamics simulations of GPR68 were used to model the changes in residue interactions and motions triggered by pH. Global and local rearrangements consistent with partial activation were observed upon protonation of the inactive state. Selected extracellular histidine and transmembrane acidic residues were found to have significantly upshifted pKa values during the simulations, consistently with their previously hypothesised role in activation through changes in protonation state. Moreover, a novel pairing between histidine and acidic residues in the extracellular region was highlighted by both sequence analyses and simulation data and tested through site-directed mutagenesis. At last, we identified a previously unknown hydrophobic lock in the extracellular region that might stabilise the inactive conformation and regulate the transition to the active state.

Structural insights into the mechanism and dynamics of proteorhodopsin biogenesis and retinal scavenging

Microbial ion-pumping rhodopsins (MRs) are extensively studied retinal-binding membrane proteins. However, their biogenesis, including oligomerisation and retinal incorporation, remains poorly understood. The bacterial green-light absorbing proton pump proteorhodopsin (GPR) has emerged as a model protein for MRs and is used here to address these open questions using cryo-electron microscopy (cryo-EM) and molecular dynamics (MD) simulations. Specifically, conflicting studies regarding GPR stoichiometry reported pentamer and hexamer mixtures without providing possible assembly mechanisms. We report the pentameric and hexameric cryo-EM structures of a GPR mutant, uncovering the role of the unprocessed N-terminal signal peptide in the assembly of hexameric GPR. Furthermore, certain proteorhodopsin-expressing bacteria lack retinal biosynthesis pathways, suggesting that they scavenge the cofactor from their environment. We shed light on this hypothesis by solving the cryo-EM structure of retinal-free proteoopsin, which together with mass spectrometry and MD simulations suggests that decanoate serves as a temporary placeholder for retinal in the chromophore binding pocket. Further MD simulations elucidate possible pathways for the exchange of decanoate and retinal, offering a mechanism for retinal scavenging. Collectively, our findings provide insights into the biogenesis of MRs, including their oligomeric assembly, variations in protomer stoichiometry and retinal incorporation through a potential cofactor scavenging mechanism.

Functional Characterization of Splice Variants in the Diagnosis of Albinism

Albinism is a genetically heterogeneous disease in which 21 genes are known so far. Its inheritance mode is autosomal recessive except for one X-linked form. The molecular analysis of exonic sequences of these genes allows for about a 70% diagnostic rate. About half (15%) of the unsolved cases are heterozygous for one pathogenic or probably pathogenic variant. Assuming that the missing variant may be located in non-coding regions, we performed sequencing for 122 such heterozygous patients of either the whole genome (27 patients) or our NGS panel (95 patients) that includes, in addition to all exons of the 21 genes, the introns and flanking sequences of five genes, TYR, OCA2, SLC45A2, GPR143 and HPS1. Rare variants (MAF < 0.01) in trans to the first variant were tested by RT-PCR and/or minigene assay. Of the 14 variants tested, nine caused either exon skipping or the inclusion of a pseudoexon, allowing for the diagnosis of 11 patients. This represents 9.8% (12/122) supplementary diagnosis for formerly unsolved patients and 75% (12/16) of those in whom the candidate variant was in trans to the first variant. Of note, one missense variant was demonstrated to cause skipping of the exon in which it is located, thus shedding new light on its pathogenic mechanism. Searching for non-coding variants and testing them for an effect on RNA splicing is warranted in order to increase the diagnostic rate.

Hit-to-Lead Optimization of Heterocyclic Carbonyloxycarboximidamides as Selective Antagonists at Human Adenosine A3 Receptor

Antagonism of the human adenosine A3 receptor (hA3R) has potential therapeutic application. Alchemical relative binding free energy calculations of K18 and K32 suggested that the combination of a 3-(2,6-dichlorophenyl)-isoxazolyl group with 2-pyridinyl at the ends of a carbonyloxycarboximidamide group should improve hA3R affinity. Of the 25 new analogues synthesized, 37 and 74 showed improved hA3R affinity compared to K18 (and K32). This was further improved through the addition of a bromine group to the 2-pyridinyl at the 5-position, generating compound 39. Alchemical relative binding free energy calculations, mutagenesis studies and MD simulations supported the compounds' binding pattern while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, so highlighting the importance of rigidification of the carbonyloxycarboximidamide moiety. MD simulations highlighted the importance of rigidification of the carbonyloxycarboximidamide, while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, which was supported through mutagenesis studies 39 also selectively antagonized endogenously expressed hA3R in nonsmall cell lung carcinoma cells, while pharmacokinetic studies indicated low toxicity enabling in vivo evaluation. We therefore suggest that 39 has potential for further development as a high-affinity hA3R antagonist.

Bilayer lipids modulate ligand binding to atypical chemokine receptor 3

Chemokine receptors belong to the large class of G protein-coupled receptors (GPCRs) and are involved in a number of (patho)physiological processes. Previous studies highlighted the importance of membrane lipids for modulating GPCR structure and function. However, the underlying mechanisms of how lipids regulate GPCRs are often poorly understood. Here, we report that anionic lipid bilayers increase the binding affinity of the chemokine CXCL12 for the atypical chemokine receptor 3 (ACKR3) by modulating the CXCL12 binding kinetics. Notably, the anionic bilayer favors CXCL12 over the more positively charged chemokine CXCL11, which we explained by bilayer interactions orienting CXCL12 but not CXCL11 for productive ACKR3 binding. Furthermore, our data suggest a stabilization of active ACKR3 conformations in anionic bilayers. Taken together, the described regulation of chemokine selectivity of ACKR3 by the lipid bilayer proposes an extended version of the classical model of chemokine binding including the lipid environment of the receptor.

Molecular Dynamics Simulations of the Eye Lens Water Channel Aquaporin 0 from Fish

Aquaporin 0 (AQP0) plays a key role in water circulation in the eye lens through a variety of functions. In contrast to mammalian genomes, zebrafish contains two aqp0 genes leading to a separation of AQP0 multiple functions between the two gene products, Aqp0a and Aqp0b. A notable feature of the zebrafish AQP0 paralogs is the increased water permeability of Aqp0b relative to Aqp0a as well as a severa lfold increase relative to mammalian AQP0. Here, we report equilibrium molecular dynamics (MD) simulations on the microsecond timescale to identify the structural basis underlying the differences in water permeability between zebrafish AQP0 paralogs and between AQP0 mammalian and fish orthologs. Our simulations are able to reproduce the experimental trends in water permeability. Our results suggest that a substitution of a key Y23 residue in mammalian AQP0 for F23 in fish AQP0 orthologs introduces significant changes in the conformational dynamics of the CS-I structural motif, which, in conjunction with different levels of hydration of the channel vestibule, can account for the differences in permeabilities between fish and mammalian AQP0 orthologs and between zebrafish AQP0 paralogs.

Piplartine attenuates aminoglycoside-induced TRPV1 activity and protects from hearing loss in mice

Hearing loss is a major health concern in our society, affecting more than 400 million people worldwide. Among the causes, aminoglycoside therapy can result in permanent hearing loss in 40% to 60% of patients receiving treatment, and despite these high numbers, no drug for preventing or treating this type of hearing loss has yet been approved by the US Food and Drug Administration. We have previously conducted high-throughput screenings of bioactive compounds, using zebrafish as our discovery platform, and identified piplartine as a potential therapeutic molecule. In the present study, we expanded this work and characterized piplartine's physicochemical and therapeutic properties. We showed that piplartine had a wide therapeutic window and neither induced nephrotoxicity in vivo in zebrafish nor interfered with aminoglycoside antibacterial activity. In addition, a fluorescence-based assay demonstrated that piplartine did not inhibit cytochrome C activity in microsomes. Coadministration of piplartine protected from kanamycin-induced hair cell loss in zebrafish and protected hearing function, outer hair cells, and presynaptic ribbons in a mouse model of kanamycin ototoxicity. Last, we investigated piplartine's mechanism of action by phospho-omics, immunoblotting, immunohistochemistry, and molecular dynamics experiments. We found an up-regulation of AKT1 signaling in the cochleas of mice cotreated with piplartine. Piplartine treatment normalized kanamycin-induced up-regulation of TRPV1 expression and modulated the gating properties of this receptor. Because aminoglycoside entrance to the inner ear is, in part, mediated by TRPV1, these results suggested that by regulating TRPV1 expression, piplartine blocked aminoglycoside's entrance, thereby preventing the long-term deleterious effects of aminoglycoside accumulation in the inner ear compartment.

Structure of the human dopamine transporter in complex with cocaine

The dopamine transporter (DAT) is crucial for regulating dopamine signalling and is the prime mediator for the rewarding and addictive effects of cocaine1. As part of the neurotransmitter sodium symporter family, DAT uses the Na+ gradient across cell membranes to transport dopamine against its chemical gradient2. The transport mechanism involves both intra- and extracellular gates that control substrate access to a central site. However, the molecular intricacies of this process and the inhibitory mechanism of cocaine have remained unclear. Here, we present the molecular structure of human DAT in complex with cocaine at a resolution of 2.66 Å. Our findings reveal that DAT adopts the expected LeuT-fold, posing in an outward-open conformation with cocaine bound at the central (S1) site. Notably, while an Na+ occupies the second Na+ site (Na2), the Na1 site seems to be vacant, with the side chain of Asn82 occupying the presumed Na+ space. This structural insight elucidates the mechanism for the cocaine inhibition of human DAT and deepens our understanding of neurotransmitter transport. By shedding light on the molecular underpinnings of how cocaine acts, our study lays a foundation for the development of targeted medications to combat addiction.

Structure of the human dopamine transporter and mechanisms of inhibition

The neurotransmitter dopamine has central roles in mood, appetite, arousal and movement1. Despite its importance in brain physiology and function, and as a target for illicit and therapeutic drugs, the human dopamine transporter (hDAT) and mechanisms by which it is inhibited by small molecules and Zn2+ are without a high-resolution structural context. Here we determine the structure of hDAT in a tripartite complex with the competitive inhibitor and cocaine analogue, (-)-2-β-carbomethoxy-3-β-(4-fluorophenyl)tropane2 (β-CFT), the non-competitive inhibitor MRS72923 and Zn2+ (ref. 4). We show how β-CFT occupies the central site, approximately halfway across the membrane, stabilizing the transporter in an outward-open conformation. MRS7292 binds to a structurally uncharacterized allosteric site, adjacent to the extracellular vestibule, sequestered underneath the extracellular loop 4 (EL4) and adjacent to transmembrane helix 1b (TM1b), acting as a wedge, precluding movement of TM1b and closure of the extracellular gate. A Zn2+ ion further stabilizes the outward-facing conformation by coupling EL4 to EL2, TM7 and TM8, thus providing specific insights into how Zn2+ restrains the movement of EL4 relative to EL2 and inhibits transport activity.

Olanzapine manipulates neuroactive signals and may onset metabolic disturbances

Olanzapine is one of the most prescribed atypical antipsychotics to treat psychiatric illness and is associated with weight gain and metabolic disturbance. The present study investigated the olanzapine-regulated metabolic pathways using functional enrichent analysis including binding affinity with G-protein-coupled receptors (GPCRs). Proteins modulated by olanzapine were retrieved from SwissTargetPrediction, DIGEP-Pred, and BindingDB and then enriched in Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) to assess molecular function, biological process, and cellular components including Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. We used homology modeling to improve the 3D structure for GPCR synapse proteins including dopamine, serotonin, muscarinic, and histamine receptors which were then optimized using molecular dynamics (MD) simulations. The protein-olanzapine binding mechanisms for different GPCR binders were evaluated using molecular docking; later refined by MD simulations. Binding mechanism of olanzapine with D2, 5HT1A, 5HT2A, 5HT2B, 5HT2C, M1, and M2 receptors were created using homology modeling and optimized using MD simulations. In target identification, it was observed that olanzapine majority targeted G-protein coupled receptors. Further, enrichment analysis identified around 76% of the total genes regulated in molecular function, biological process, and cellular components were common including KEGG pathways. Moreover, it was observed that olanzapine had a major potency over the neurotransmitter synapse including neuroactive signals . Olanzapine-induced weight gain and metabolic alterations could be due to the deregulation of multiple synapses like dopamine, serotonin, muscarinic, and histamine at the feeding center followed by cGMP-PKG, cAMP, and PI3K-Akt signaling pathways. HIGHLIGHTSOlanzapine is used in the management of psychiatric illnesses.Olanzapine causes disturbance in lipids and glucosehomeostasis and manipulates energy expenditure.Olanzapine-induced weight gain may occur due to the deregulation of the multiple synapse and cGMP-PKG, cAMP, and PI3K-Akt signaling pathwayCommunicated by Ramaswamy H. Sarma.

Deep learning structural insights into heterotrimeric alternatively spliced P2X7 receptors

P2X7 receptors (P2X7Rs) are membrane-bound ATP-gated ion channels that are composed of three subunits. Different subunit structures may be expressed due to alternative splicing of the P2RX7 gene, altering the receptor's function when combined with the wild-type P2X7A subunits. In this study, the application of the deep-learning method, AlphaFold2-Multimer (AF2M), for the generation of trimeric P2X7Rs was validated by comparing an AF2M-generated rat wild-type P2X7A receptor with a structure determined by cryogenic electron microscopy (cryo-EM) (Protein Data Bank Identification: 6U9V). The results suggested AF2M could firstly, accurately predict the structures of P2X7Rs and secondly, accurately identify the highest quality model through the ranking system. Subsequently, AF2M was used to generate models of heterotrimeric alternatively spliced P2X7Rs consisting of one or two wild-type P2X7A subunits in combination with one or two P2X7B, P2X7E, P2X7J, and P2X7L splice variant subunits. The top-ranking models were deemed valid based on AF2M's confidence measures, stability in molecular dynamics simulations, and consistent flexibility of the conserved regions between the models. The structure of the heterotrimeric receptors, which were missing key residues in the ATP binding sites and carboxyl terminal domains (CTDs) compared to the wild-type receptor, help to explain their observed functions. Overall, the models produced in this study (available as supplementary material) unlock the possibility of structure-based studies into the heterotrimeric P2X7Rs.

Effect of substituents at the C3´, C3´N, C10 and C2-meta-benzoate positions of taxane derivatives on their activity against resistant cancer cells

We tested the effect of substituents at the (1) C3´, C3´N, (2) C10, and (3) C2-meta-benzoate positions of taxane derivatives on their activity against sensitive versus counterpart paclitaxel-resistant breast (MCF-7) and ovarian (SK-OV-3) cancer cells. We found that (1) non-aromatic groups at both C3´ and C3´N positions, when compared with phenyl groups at the same positions of a taxane derivative, significantly reduced the resistance of ABCB1 expressing MCF-7/PacR and SK-OV-3/PacR cancer cells. This is, at least in the case of the SB-T-1216 series, accompanied by an ineffective decrease of intracellular levels in MCF-7/PacR cells. The low binding affinity of SB-T-1216 in the ABCB1 binding cavity can elucidate these effects. (2) Cyclopropanecarbonyl group at the C10 position, when compared with the H atom, seems to increase the potency and capability of the derivative in overcoming paclitaxel resistance in both models. (3) Derivatives with fluorine and methyl substituents at the C2-meta-benzoate position were variously potent against sensitive and resistant cancer cells. All C2 derivatives were less capable of overcoming acquired resistance to paclitaxel in vitro than non-substituted analogs. Notably, fluorine derivatives SB-T-121205 and 121,206 were more potent against sensitive and resistant SK-OV-3 cells, and derivatives SB-T-121405 and 121,406 were more potent against sensitive and resistant MCF-7 cells. (4) The various structure-activity relationships of SB-T derivatives observed in two cell line models known to express ABCB1 favor their complex interaction not based solely on ABCB1.

PI3 kinase inhibitor PI828 uncouples aminergic GPCRs and Ca2+ mobilization irrespectively of its primary target

The phosphoinositide 3-kinase (PI3K) is involved in regulation of multiple intracellular processes. Although the inhibitory analysis is generally employed for validating a physiological role of PI3K, increasing body of evidence suggests that PI3K inhibitors can exhibit PI3K-unrelated activity as well. Here we studied Ca2+ signaling initiated by aminergic agonists in a variety of different cells and analyzed effects of the PI3K inhibitor PI828 on cell responsiveness. It turned out that PI828 inhibited Ca2+ transients elicited by acetylcholine (ACh), histamine, and serotonin, but did not affect Ca2+ responses to norepinephrine and ATP. Another PI3K inhibitor wortmannin negligibly affected Ca2+ signaling initiated by any one of the tested agonists. Using the genetically encoded PIP3 sensor PH(Akt)-Venus, we confirmed that both PI828 and wortmannin effectively inhibited PI3K and ascertained that this kinase negligibly contributed to ACh transduction. These findings suggested that PI828 inhibited Ca2+ responses to aminergic agonists tested, involving an unknown cellular mechanism unrelated to the PI3K inhibition. Complementary physiological experiments provided evidence that PI828 could inhibit Ca2+ signals induced by certain agonists, by acting extracellularly, presumably, through their surface receptors. For the muscarinic M3 receptor, this possibility was verified with molecular docking and molecular dynamics. As demonstrated with these tools, wortmannin could be bound in the extracellular vestibule at the muscarinic M3 receptor but this did not preclude binding of ACh to the M3 receptor followed by its activation. In contrast, PI828 could sterically block the passage of ACh into the allosteric site, preventing activation of the muscarinic M3 receptor.

Alternative splicing generates a novel ferroportin isoform with a shorter C-terminal and an intact iron- and hepcidin-binding property

Ferroportin (FPN) is a transmembrane protein and is the only known iron exporter that helps in maintaining iron homeostasis in vertebrates. To maintain stable iron equilibrium in the body, ferroportin works in conjunction with a peptide called hepcidin. In this study, we have identified an alternatively spliced novel isoform of the human SLC40A1 gene, which encodes for the FPN protein and is found to be expressed in different tissues. The novel transcript has an alternate last exon and encodes 31-amino acid long peptide sequence that replaces 104 amino acids at C-terminal in the novel transcript. Molecular modelling and molecular dynamics (MD) simulation studies revealed key structural features of the novel isoform (FPN-N). FPN-N was predicted to have 12 transmembrane domains similar to the reported isoform (FPN), despite being much smaller in size. FPN-N was found to interact with hepcidin, a key regulator of ferroportin activity. Also, the iron-binding sites were retained in the novel isoform as revealed by the MD simulation of FPN-N in bilipid membrane. The novel isoform identified in this study may play important role in iron homeostasis. However, further studies are required to characterize the FPN-N isoform and decipher its role inside the cell.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Activation of the influenza B M2 proton channel (BM2)

Influenza B viruses have co-circulated during most seasonal flu epidemics and can cause significant human morbidity and mortality due to their rapid mutation, emerging drug resistance, and severe impact on vulnerable populations. The influenza B M2 proton channel (BM2) plays an essential role in viral replication, but the mechanisms behind its symmetric proton conductance and the involvement of a second histidine (His27) cluster remain unclear. Here we perform the membrane-enabled continuous constant-pH molecular dynamics simulations on wildtype BM2 and a key H27A mutant to explore its pH-dependent conformational switch. Simulations capture the activation as the first histidine (His19) protonates and reveal the transition at lower pH values compared to AM2 is a result of electrostatic repulsions between His19 and pre-protonated His27. Crucially, we provide an atomic-level understanding of the symmetric proton conduction by identifying pre-activating channel hydration in the C-terminal portion. This research advances our understanding of the function of BM2 function and lays the groundwork for further chemically reactive modeling of the explicit proton transport process as well as possible anti-flu drug design efforts.

Atomic insights into the signaling landscape of E. coli PhoQ histidine kinase from molecular dynamics simulations

Bacteria rely on two-component systems to sense environmental cues and regulate gene expression for adaptation. The PhoQ/PhoP system exemplifies this crucial role, playing a key part in sensing magnesium (Mg2+) levels, antimicrobial peptides, mild acidic pH, osmotic upshift, and long-chain unsaturated fatty acids, promoting virulence in certain bacterial species. However, the precise details of PhoQ activation remain elusive. To elucidate PhoQ's signaling mechanism at atomic resolution, we combined AlphaFold2 predictions with molecular modeling and carried out extensive Molecular Dynamics (MD) simulations. Our MD simulations revealed three distinct PhoQ conformations that were validated by experimental data. Notably, one conformation was characterized by Mg2+ bridging the acidic patch in the sensor domain to the membrane, potentially representing a repressed state. Furthermore, the high hydration observed in a putative intermediate state lends support to the hypothesis of water-mediated conformational changes during PhoQ signaling. Our findings not only revealed specific conformations within the PhoQ signaling pathway, but also hold significant promise for understanding the broader histidine kinase family due to their shared structural features. Our approach paves the way for a more comprehensive understanding of histidine kinase signaling mechanisms across various bacterial species and opens the door for developing novel therapeutics that target PhoQ modulation.

Interior pH-sensing residue of human voltage-gated proton channel Hv1 is histidine 168

The molecular mechanisms governing the human voltage-gated proton channel hHv1 remain elusive. Here, we used membrane-enabled hybrid-solvent continuous constant pH molecular dynamics (CpHMD) simulations with pH replica exchange to further evaluate the structural models of hHv1 in the closed (hyperpolarized) and open (depolarized) states recently obtained with MD simulations and explore potential pH-sensing residues. The CpHMD titration at a set of symmetric pH conditions revealed three residues that can gain or lose protons upon channel depolarization. Among them, residue H168 at the intracellular end of the S3 helix switches from the deprotonated to the protonated state and its protonation is correlated with the increased tilting of the S3 helix during the transition from the closed to the open state. Thus, the simulation data suggest H168 as an interior pH sensor, in support of a recent finding based on electrophysiological experiments of Hv1 mutants. We propose that protonation of H168 acts as a key that unlocks the closed channel configuration by increasing the flexibility of the S2-S3 linker, which increases the tilt angle of S3 and enhances the mobility of the S4 helix, thus promoting channel opening. Our work represents an important step toward deciphering the pH-dependent gating mechanism of hHv1.

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

Allosterism in the adenosine A2A and cannabinoid CB2 heteromer

BACKGROUND AND PURPOSE

Allosterism is a regulatory mechanism for GPCRs that can be attained by ligand-binding or protein-protein interactions with another GPCR. We have studied the influence of the dimer interface on the allosteric properties of the A2A receptor and CB2 receptor heteromer.

EXPERIMENTAL APPROACH

We have evaluated cAMP production, phosphorylation of signal-regulated kinases (pERK1/2), label-free dynamic mass redistribution, β-arrestin 2 recruitment and bimolecular fluorescence complementation assays in the absence and presence of synthetic peptides that disrupt the formation of the heteromer. Molecular dynamic simulations provided converging evidence that the heteromeric interface influences the allosteric properties of the A2AR-CB2R heteromer.

KEY RESULTS

Apo A2AR blocks agonist-induced signalling of CB2R. The disruptive peptides, with the amino acid sequence of transmembrane (TM) 6 of A2AR or CB2R, facilitate CB2R activation, suggesting that A2AR allosterically prevents the outward movement of TM 6 of CB2R for G protein binding. Significantly, binding of the selective antagonist SCH 58261 to A2AR also facilitated agonist-induced activation of CB2R.

CONCLUSIONS AND IMPLICATIONS

It is proposed that the A2AR-CB2R heteromer contains distinct dimerization interfaces that govern its functional properties. The molecular interface between protomers of the A2AR-CB2R heteromer interconverted from TM 6 for apo or agonist-bound A2AR, blocking CB2R activation, to mainly the TM 1/7 interface for antagonist-bound A2AR, facilitating the independent opening of intracellular cavities for G protein binding. These novel results shed light on a different type of allosteric mechanism and extend the repertoire of GPCR heteromer signalling.

Structural insights into ligand recognition and activation of the succinate receptor SUCNR1

Succinate, a citric acid cycle intermediate, serves important functions in energy homeostasis and metabolic regulation. Extracellular succinate acts as a stress signal through succinate receptor (SUCNR1), a class A G protein-coupled receptor. Research on succinate signaling is hampered by the lack of high-resolution structures of the agonist-bound receptor. We present cryoelectron microscopy (cryo-EM) structures of SUCNR1-Gi complexes bound to succinate and its non-metabolite derivative cis-epoxysuccinate. Key determinants for the recognition of succinate in cis conformation include R2817.39 and Y832.64, while Y301.39 and R993.29 participate in the binding of both succinate and cis-epoxysuccinate. Extracellular loop 2, through F175ECL2 in its β-hairpin, forms a hydrogen bond with succinate and caps the binding pocket. At the receptor-Gi interface, agonist binding induces the rearrangement of a hydrophobic network on transmembrane (TM)5 and TM6, leading to TM signaling through TM3 and TM7. These findings extend our understanding of succinate recognition by SUCNR1, aiding the development of therapeutics for the succinate receptor.

Lipid Organization by the Caveolin-1 Complex

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S-complex. The CAV1 8S-complex comprises a disk-like structure, about 15nm in diameter, with a central beta barrel. Further oligomerization of 8S-complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration. However, the molecular mechanisms behind membrane remodeling and cholesterol filtering are still not understood. Performing atomistic Molecular Dynamics simulations in combination with advanced sampling techniques, we describe how the CAV1-8S complex bends the membrane and accumulates cholesterol. Here, our simulations show an enhancing effect by the palmitoylations of CAV1, and we predict that the CAV1-8S complex can extract cholesterol molecules from the lipid bilayer and accommodate them in its beta barrel. Through backmapping to the all-atom level we also conclude that the Martini v2 coarse-grained forcefield overestimates membrane bending, as the atomistic simulations exhibit only very localized bending.

Exploring bat-inspired cyclic tryptophan diketopiperazines as ABCB1 Inhibitors

Chemotherapy-induced drug resistance remains a major cause of cancer recurrence and patient mortality. ATP binding cassette subfamily B member 1 (ABCB1) transporter overexpression in tumors contributes to resistance, yet current ABCB1 inhibitors have been unsuccessful in clinical trials. To address this challenge, we propose a new strategy using tryptophan as a lead molecule for developing ABCB1 inhibitors. Our idea stems from our studies on bat cells, as bats have low cancer incidences and high ABCB1 expression. We hypothesized that potential ABCB1 substrates in bats could act as competitive inhibitors in humans. By molecular simulations of ABCB1-substrate interactions, we generated a benzylated Cyclo-tryptophan (C3N-Dbn-Trp2) that inhibits ABCB1 activity with efficacy comparable to or better than the classical inhibitor, verapamil. C3N-Dbn-Trp2 restored chemotherapy sensitivity in drug-resistant human cancer cells with no adverse effect on cell proliferation. Our unique approach presents a promising lead toward developing effective ABCB1 inhibitors to treat drug-resistant cancers.

Vaccine design and development: Exploring the interface with computational biology and AI

Computational biology involves applying computer science and informatics techniques in biology to understand complex biological data. It allows us to collect, connect, and analyze biological data at a large scale and build predictive models. In the twenty first century, computational resources along with Artificial Intelligence (AI) have been widely used in various fields of biological sciences such as biochemistry, structural biology, immunology, microbiology, and genomics to handle massive data for decision-making, including in applications such as drug design and vaccine development, one of the major areas of focus for human and animal welfare. The knowledge of available computational resources and AI-enabled tools in vaccine design and development can improve our ability to conduct cutting-edge research. Therefore, this review article aims to summarize important computational resources and AI-based tools. Further, the article discusses the various applications and limitations of AI tools in vaccine development.

Rejuvenation of Meropenem by Conjugation with Tilapia Piscidin-4 Peptide Targeting NDM-1 Escherichia coli

Gram-negative pathogens that produce β-lactamases pose a serious public health threat as they can render β-lactam antibiotics inactive via hydrolysis. This action contributes to the waning effectiveness of clinical antibiotics and creates an urgent need for new antimicrobials. Antimicrobial peptides (AMPs) exhibiting multimodal functions serve as a potential source in spite of a few limitations. Thus, the conjugation of conventional antibiotics with AMPs may be an effective strategy to leverage the advantages of each component. In this study, we conjugated meropenem to the AMP Tilapia piscidin 4 (TP4) using a typical coupling reaction. The conjugate was characterized by using HPLC-MS, HR-MS, and MS-MS fragmentation analysis. It was then evaluated in terms of antibacterial potency, hemolysis, and cytotoxicity toward RAW264.7 and CCD-966SK cell lines. The conjugation of meropenem with TP4 significantly reduced the cytotoxicity compared to TP4. Conjugation of unprotected TP4 with meropenem resulted in cross-linking at the N-terminal and lysine sites. The structural activity relationship of the two isomers of the TP4-meropenem conjugate was investigated. Both the isomers showed notable antibacterial activities against NDM-1 Escherichia coli and reduced red blood cell hemolysis as compared to TP4. Lysine conjugate (TP4-K-Mero) showed lesser hemolysis than the N-terminal conjugate (TP4-N-Mero). Molecular modeling further revealed that the conjugates can bind to lipopolysaccharides and inhibit NDM-1 β-lactamase. Together, these data show that conjugation of antibiotics with AMP can be a feasible approach to increase the therapeutic profile and effectively target multidrug-resistant pathogens. Furthermore, antibiotic conjugation at different AMP sites tends to show unique biological properties.

POTRA domains of the TamA insertase interact with the outer membrane and modulate membrane properties

The outer membrane (OM) of gram-negative bacteria serves as a vital organelle that is densely populated with OM proteins (OMPs) and plays pivotal roles in cellular functions and virulence. The assembly and insertion of these OMPs into the OM represent a fundamental process requiring specialized molecular chaperones. One example is the translocation and assembly module (TAM), which functions as a transenvelope chaperone promoting the folding of specific autotransporters, adhesins, and secretion systems. The catalytic unit of TAM, TamA, comprises a catalytic β-barrel domain anchored within the OM and three periplasmic polypeptide-transport-associated (POTRA) domains that recruit the TamB subunit. The latter acts as a periplasmic ladder that facilitates the transport of unfolded OMPs across the periplasm. In addition to their role in recruiting the auxiliary protein TamB, our data demonstrate that the POTRA domains mediate interactions with the inner surface of the OM, ultimately modulating the membrane properties. Through the integration of X-ray crystallography, molecular dynamic simulations, and biomolecular interaction methodologies, we located the membrane-binding site on the first and second POTRA domains. Our data highlight a binding preference for phosphatidylglycerol, a minor lipid constituent present in the OM, which has been previously reported to facilitate OMP assembly. In the context of the densely OMP-populated membrane, this association may serve as a mechanism to secure lipid accessibility for nascent OMPs through steric interactions with existing OMPs, in addition to creating favorable conditions for OMP biogenesis.

Structural Studies of the Taurine Transporter: A Potential Biological Target from the GABA Transporter Subfamily in Cancer Therapy

The taurine transporter (TauT, SLC6A6) is a member of the solute carrier 6 (SLC6) family, which plays multiple physiological roles. The SLC6 family is divided into four subfamilies: GABA (γ-aminobutyric acid), monoamine, glycine and neutral amino acid transporters. Proteins from the GABA group, including the taurine transporter, are primarily considered therapeutic targets for treating central nervous system disorders. However, recent studies have suggested that inhibitors of SLC6A6 could also serve as anticancer agents. Overexpression of TauT has been associated with the progression of colon and gastric cancer. The pool of known ligands of this transporter is limited and the exact spatial structure of taurine transporter remains unsolved. Understanding its structure could aid in the development of novel inhibitors. Therefore, we utilized homology modelling techniques to create models of TauT. Docking studies and molecular dynamics simulations were conducted to describe protein-ligand interactions. We compared the obtained information for TauT with literature data on other members of the GABA transporter group. Our in silico analysis allowed us to characterize the transporter structure and point out amino acids crucial for ligand binding: Glu406, Gly62 and Tyr138. The significance of selected residues was confirmed through structural studies of mutants. These results will aid in the development of novel taurine transporter inhibitors, which can be explored as anticancer agents.

A Membrane-Assisted Mechanism for the Release of Ceramide from the CERT START Domain

Ceramide transfer protein CERT is the mediator of nonvesicular transfer of ceramide from the ER to Golgi. In CERT, START is the domain responsible for the binding and transport of ceramide. A wealth of structural data has revealed a helix-grip fold surrounding a large hydrophobic cavity holding the ceramide. Yet, little is known about the mechanisms by which START releases the ceramide through the polar region and into the packed environment of cellular membranes. As such events do not lend themselves easily to experimental investigations, we used multiple unbiased microsecond-long molecular simulations. We propose a membrane-assisted mechanism in which the membrane acts as an allosteric effector initiating the release of ceramide and where the passage of the ceramide acyl chains is facilitated by the intercalation of a single phosphatidylcholine lipid in the cavity, practically greasing the ceramide way out. We verify using free energy calculation and experimental lipidomics data that CERT forms stable complexes with phosphatidylcholine lipids, in addition to ceramide, thus providing validation for the proposed mechanism.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Functional Characterization of the Lysosomal Peptide/Histidine Transporter PHT1 (SLC15A4) by Solid Supported Membrane Electrophysiology (SSME)

The peptide/histidine transporter PHT1 (SLC15A4) is expressed in the lysosomal membranes of immune cells where it plays an important role in metabolic and inflammatory signaling. PHT1 is an H+-coupled/histidine symporter that can transport a wide range of oligopeptides, including a variety of bacterial-derived peptides. Moreover, it enables the scaffolding of various metabolic signaling molecules and interacts with key regulatory elements of the immune response. Not surprisingly, PHT1 has been implicated in the pathogenesis of autoimmune diseases such as systemic lupus erythematosus (SLE). Unfortunately, the pharmacological development of PHT1 modulators has been hampered by the lack of suitable transport assays. To address this shortcoming, a novel transport assay based on solid-supported membrane-based electrophysiology (SSME) is presented. Key findings of the present SSME studies include the first recordings of electrophysiological properties, a pH dependence analysis, an assessment of PHT1 substrate selectivity, as well as the transport kinetics of the identified substrates. In contrast to previous work, PHT1 is studied in its native lysosomal environment. Moreover, observed substrate selectivity is validated by molecular docking. Overall, this new SSME-based assay is expected to contribute to unlocking the pharmacological potential of PHT1 and to deepen the understanding of its functional properties.

Structural basis for odorant recognition of the insect odorant receptor OR-Orco heterocomplex

Insects detect and discriminate a diverse array of chemicals using odorant receptors (ORs), which are ligand-gated ion channels comprising a divergent odorant-sensing OR and a conserved odorant receptor co-receptor (Orco). In this work, we report structures of the ApOR5-Orco heterocomplex from the pea aphid Acyrthosiphon pisum alone and bound to its known activating ligand, geranyl acetate. In these structures, three ApOrco subunits serve as scaffold components that cannot bind the ligand and remain relatively unchanged. Upon ligand binding, the pore-forming helix S7b of ApOR5 shifts outward from the central pore axis, causing an asymmetrical pore opening for ion influx. Our study provides insights into odorant recognition and channel gating of the OR-Orco heterocomplex and offers structural resources to support development of innovative insecticides and repellents for pest control.

Variety of Plant Oils: Species-Specific Lipid Biosynthesis

Plant oils represent a large group of neutral lipids with important applications in food, feed and oleochemical industries. Most plants accumulate oils in the form of triacylglycerol within seeds and their surrounding tissues, which comprises three fatty acids attached to a glycerol backbone. Different plant species accumulate unique fatty acids in their oils, serving a range of applications in pharmaceuticals and oleochemicals. To enable the production of these distinctive oils, select plant species have adapted specialized oil metabolism pathways, involving differential gene co-expression networks and structurally divergent enzymes/proteins. Here, we summarize some of the recent advances in our understanding of oil biosynthesis in plants. We compare expression patterns of oil metabolism genes from representative species, including Arabidopsis thaliana, Ricinus communis (castor bean), Linum usitatissimum L. (flax) and Elaeis guineensis (oil palm) to showcase the co-expression networks of relevant genes for acyl metabolism. We also review several divergent enzymes/proteins associated with key catalytic steps of unique oil accumulation, including fatty acid desaturases, diacylglycerol acyltransferases and oleosins, highlighting their structural features and preference toward unique lipid substrates. Lastly, we briefly discuss protein interactomes and substrate channeling for oil biosynthesis and the complex regulation of these processes.

Deciphering the conformational dynamics of Myelin Oligodendrocyte glycoprotein in the myelin sheath

Myelin Oligodendrocyte Glycoprotein (MOG) is a transmembrane protein in the myelin sheath. It acts as an auto-antigen under certain unknown conditions causing demyelination, thus resulting in Myelin Oligodendrocyte Glycoprotein Antibody-associated Disease (MOGAD). The significance of glycosylation in the conformational dynamics of the extracellular region (EC1) of the MOG were evident from the previous computational studies. Here, in this study, we performed the molecular dynamics simulation of the entire human MOG in the myelin sheath for 100 ns using the NAMD program. The results indicated that the EC1 and cytoplasmic region (CP) dominate the conformational rigidity of the protein, and enhance its interaction with lipids. This in turn helps in maintaining the myelin integrity in the presence of glycan. The transmembrane regions have reduced interaction with lipids in the glycosylated system. Moreover, the C-terminal extracellular region 2 (EC2) behaves exactly opposite to that of EC1 in the glycan presence. This may be attributed to the glycosylation site in the EC1 region. Hence, not only the region EC1 (having 3 crucial epitopes) but even the CP region were important for understanding the proper function of MOG in the glycan presence.Communicated by Ramaswamy H. Sarma.

Structural basis of tolvaptan binding to the vasopressin V2 receptor

The vasopressin V2 receptor (V2R) is a validated therapeutic target for autosomal dominant polycystic kidney disease (ADPKD), with tolvaptan being the first FDA-approved antagonist. Herein, we used Gaussian accelerated molecular dynamics simulations to investigate the spontaneous binding of tolvaptan to both active and inactive V2R conformations at the atomic-level. Overall, the binding process consists of two stages. Tolvaptan binds initially to extracellular loops 2 and 3 (ECL2/3) before overcoming an energy barrier to enter the pocket. Our simulations result highlighted key residues (e.g., R181, Y205, F287, F178) involved in this process, which were experimentally confirmed by site-directed mutagenesis. This work provides structural insights into tolvaptan-V2R interactions, potentially aiding the design of novel antagonists for V2R and other G protein-coupled receptors.

Structures of synaptic vesicle protein 2A and 2B bound to anticonvulsants

Epilepsy is a common neurological disorder characterized by abnormal activity of neuronal networks, leading to seizures. The racetam class of anti-seizure medications bind specifically to a membrane protein found in the synaptic vesicles of neurons called synaptic vesicle protein 2 (SV2) A (SV2A). SV2A belongs to an orphan subfamily of the solute carrier 22 organic ion transporter family that also includes SV2B and SV2C. The molecular basis for how anti-seizure medications act on SV2s remains unknown. Here we report cryo-electron microscopy structures of SV2A and SV2B captured in a luminal-occluded conformation complexed with anticonvulsant ligands. The conformation bound by anticonvulsants resembles an inhibited transporter with closed luminal and intracellular gates. Anticonvulsants bind to a highly conserved central site in SV2s. These structures provide blueprints for future drug design and will facilitate future investigations into the biological function of SV2s.

PIP2 modulates TRPC3 activity via TRP helix and S4-S5 linker

The transient receptor potential canonical type 3 (TRPC3) channel plays a pivotal role in regulating neuronal excitability in the brain via its constitutive activity. The channel is intricately regulated by lipids and has previously been demonstrated to be positively modulated by PIP2. Using molecular dynamics simulations and patch clamp techniques, we reveal that PIP2 predominantly interacts with TRPC3 at the L3 lipid binding site, located at the intersection of pre-S1 and S1 helices. We demonstrate that PIP2 sensing involves a multistep mechanism that propagates from L3 to the pore domain via a salt bridge between the TRP helix and S4-S5 linker. Notably, we find that both stimulated and constitutive TRPC3 activity require PIP2. These structural insights into the function of TRPC3 are invaluable for understanding the role of the TRPC subfamily in health and disease, in particular for cardiovascular diseases, in which TRPC3 channels play a major role.

TMC6 functions as a GPCR-like receptor to sense noxious heat via Gαq signaling

Thermosensation is vital for the survival, propagation, and adaption of all organisms, but its mechanism is not fully understood yet. Here, we find that TMC6, a membrane protein of unknown function, is highly expressed in dorsal root ganglion (DRG) neurons and functions as a Gαq-coupled G protein-coupled receptor (GPCR)-like receptor to sense noxious heat. TMC6-deficient mice display a substantial impairment in noxious heat sensation while maintaining normal perception of cold, warmth, touch, and mechanical pain. Further studies show that TMC6 interacts with Gαq via its intracellular C-terminal region spanning Ser780 to Pro810. Specifically disrupting such interaction using polypeptide in DRG neurons, genetically ablating Gαq, or pharmacologically blocking Gαq-coupled GPCR signaling can replicate the phenotype of TMC6 deficient mice regarding noxious heat sensation. Noxious heat stimulation triggers intracellular calcium release from the endoplasmic reticulum (ER) of TMC6- but not control vector-transfected HEK293T cell, which can be significantly inhibited by blocking PLC or IP3R. Consistently, noxious heat-induced intracellular Ca2+ release from ER and action potentials of DRG neurons largely reduced when ablating TMC6 or blocking Gαq/PLC/IP3R signaling pathway as well. In summary, our findings indicate that TMC6 can directly function as a Gαq-coupled GPCR-like receptor sensing noxious heat.

Structural and biophysical insights into targeting of claudin-4 by a synthetic antibody fragment

Claudins are a 27-member family of ~25 kDa membrane proteins that integrate into tight junctions to form molecular barriers at the paracellular spaces between endothelial and epithelial cells. As the backbone of tight junction structure and function, claudins are attractive targets for modulating tissue permeability to deliver drugs or treat disease. However, structures of claudins are limited due to their small sizes and physicochemical properties-these traits also make therapy development a challenge. Here we report the development of a synthetic antibody fragment (sFab) that binds human claudin-4 and the determination of a high-resolution structure of it bound to claudin-4/enterotoxin complexes using cryogenic electron microscopy. Structural and biophysical results reveal this sFabs mechanism of select binding to human claudin-4 over other homologous claudins and establish the ability of sFabs to bind hard-to-target claudins to probe tight junction structure and function. The findings provide a framework for tight junction modulation by sFabs for tissue-selective therapies.

The novel function of an orphan pheromone receptor reveals the sensory specializations of two potential distinct types of sex pheromones in noctuid moth

Sex pheromones play crucial role in mating behavior of moths, involving intricate recognition mechanisms. While insect chemical biology has extensively studied type I pheromones, type II pheromones remain largely unexplored. This study focused on Helicoverpa armigera, a representative species of noctuid moth, aiming to reassess its sex pheromone composition. Our research unveiled two previously unidentified candidate type II sex pheromones-3Z,6Z,9Z-21:H and 3Z,6Z,9Z-23:H-in H. armigera. Furthermore, we identified HarmOR11 as an orphan pheromone receptor of 3Z,6Z,9Z-21:H. Through AlphaFold2 structural prediction, molecular docking, and molecular dynamics simulations, we elucidated the structural basis and key residues governing the sensory nuances of both type I and type II pheromone receptors, particularly HarmOR11 and HarmOR13. This study not only reveals the presence and recognition of candidate type II pheromones in a noctuid moth, but also establishes a comprehensive structural framework for PRs, contributing to the understanding of connections between evolutionary adaptations and the emergence of new pheromone types.

In Silico Assisted Identification, Synthesis, and In Vitro Pharmacological Characterization of Potent and Selective Blockers of the Epilepsy-Associated KCNT1 Channel

Gain-of-function (GoF) variants in KCNT1 channels cause severe, drug-resistant forms of epilepsy. Quinidine is a known KCNT1 blocker, but its clinical use is limited due to severe drawbacks. To identify novel KCNT1 blockers, a homology model of human KCNT1 was built and used to screen an in-house library of compounds. Among the 20 molecules selected, five (CPK4, 13, 16, 18, and 20) showed strong KCNT1-blocking ability in an in vitro fluorescence-based assay. Patch-clamp experiments confirmed a higher KCNT1-blocking potency of these compounds when compared to quinidine, and their selectivity for KCNT1 over hERG and Kv7.2 channels. Among identified molecules, CPK20 displayed the highest metabolic stability; this compound also blocked KCNT2 currents, although with a lower potency, and counteracted GoF effects prompted by 2 recurrent epilepsy-causing KCNT1 variants (G288S and A934T). The present results provide solid rational basis for future design of novel compounds to counteract KCNT1-related neurological disorders.

Atomistic description of the OCTN1 recognition mechanism via in silico methods

The Organic Cation Transporter Novel 1 (OCTN1), also known as SLC22A4, is widely expressed in various human tissues, and involved in numerous physiological and pathological processes remains. It facilitates the transport of organic cations, zwitterions, with selectivity for positively charged solutes. Ergothioneine, an antioxidant compound, and acetylcholine (Ach) are among its substrates. Given the lack of experimentally solved structures of this protein, this study aimed at generating a reliable 3D model of OCTN1 to shed light on its substrate-binding preferences and the role of sodium in substrate recognition and transport. A chimeric model was built by grafting the large extracellular loop 1 (EL1) from an AlphaFold-generated model onto a homology model. Molecular dynamics simulations revealed domain-specific mobility, with EL1 exhibiting the highest impact on overall stability. Molecular docking simulations identified cytarabine and verapamil as highest affinity ligands, consistent with their known inhibitory effects on OCTN1. Furthermore, MM/GBSA analysis allowed the categorization of substrates into weak, good, and strong binders, with molecular weight strongly correlating with binding affinity to the recognition site. Key recognition residues, including Tyr211, Glu381, and Arg469, were identified through interaction analysis. Ach demonstrated a low interaction energy, supporting the hypothesis of its one-directional transport towards to outside of the membrane. Regarding the role of sodium, our model suggested the involvement of Glu381 in sodium binding. Molecular dynamics simulations of systems at increasing levels of Na+ concentrations revealed increased sodium occupancy around Glu381, supporting experimental data associating Na+ concentration to molecule transport. In conclusion, this study provides valuable insights into the 3D structure of OCTN1, its substrate-binding preferences, and the role of sodium in the recognition. These findings contribute to the understanding of OCTN1 involvement in various physiological and pathological processes and may have implications for drug development and disease management.

Molecular Determinant Underlying Selective Coupling of Primary G-Protein by Class A GPCRs

G-protein-coupled receptors (GPCRs) transmit downstream signals predominantly via G-protein pathways. However, the conformational basis of selective coupling of primary G-protein remains elusive. Histamine receptors H2R and H3R couple with Gs- or Gi-proteins respectively. Here, three cryo-EM structures of H2R-Gs and H3R-Gi complexes are presented at a global resolution of 2.6-2.7 Å. These structures reveal the unique binding pose for endogenous histamine in H3R, wherein the amino group interacts with E2065.46 of H3R instead of the conserved D1143.32 of other aminergic receptors. Furthermore, comparative analysis of the H2R-Gs and H3R-Gi complexes reveals that the structural geometry of TM5/TM6 determines the primary G-protein selectivity in histamine receptors. Machine learning (ML)-based structuromic profiling and functional analysis of class A GPCR-G-protein complexes illustrate that TM5 length, TM5 tilt, and TM6 outward movement are key determinants of the Gs and Gi/o selectivity among the whole Class A family. Collectively, the findings uncover the common structural geometry within class A GPCRs that determines the primary Gs- and Gi/o-coupling selectivity.

Cryo-EM structure of human sphingomyelin synthase and its mechanistic implications for sphingomyelin synthesis

Sphingomyelin (SM) has key roles in modulating mammalian membrane properties and serves as an important pool for bioactive molecules. SM biosynthesis is mediated by the sphingomyelin synthase (SMS) family, comprising SMS1, SMS2 and SMS-related (SMSr) members. Although SMS1 and SMS2 exhibit SMS activity, SMSr possesses ceramide phosphoethanolamine synthase activity. Here we determined the cryo-electron microscopic structures of human SMSr in complexes with ceramide, diacylglycerol/phosphoethanolamine and ceramide/phosphoethanolamine (CPE). The structures revealed a hexameric arrangement with a reaction chamber located between the transmembrane helices. Within this structure, a catalytic pentad E-H/D-H-D was identified, situated at the interface between the lipophilic and hydrophilic segments of the reaction chamber. Additionally, the study unveiled the two-step synthesis process catalyzed by SMSr, involving PE-PLC (phosphatidylethanolamine-phospholipase C) hydrolysis and the subsequent transfer of the phosphoethanolamine moiety to ceramide. This research provides insights into the catalytic mechanism of SMSr and expands our understanding of sphingolipid metabolism.

Orphan receptor GPR88 as a potential therapeutic target for CNS disorders - an in silico approach

The G-protein-coupled receptors are a part of the largest and most physiologically relevant family of membrane proteins. One-third of the medications, now on the market, target the GPCR receptor family, which is one of the most important therapeutic targets for many disorders. In the reported work, we have focussed on orphan GPR88 receptor which is a part of the GPCR protein family and a potential target for central nervous system disorders. GPR88 is known to show the highest expression in the striatum, which is a key region in motor control and cognitive functions. Recent studies have reported that GPR88 is activated by two agonists, 2-PCCA and RTI-13951-33. In this study, we have predicted the three-dimensional protein structure for the orphan GPR88 using the homology modeling approach. We then used shape-based screening techniques based on known agonists and structure-based virtual screening methods employing docking to uncover novel GPR88 ligands. The screened GPR88-ligand complexes were further subjected to molecular dynamics simulation studies. The selected ligands could fasten the development of novel treatments for the vast list of movement and central nervous system disorders.Communicated by Ramaswamy H. Sarma.

One-shot design elevates functional expression levels of a voltage-gated potassium channel

Membrane proteins play critical physiological roles as receptors, channels, pumps, and transporters. Despite their importance, however, low expression levels often hamper the experimental characterization of membrane proteins. We present an automated and web-accessible design algorithm called mPROSS (https://mPROSS.weizmann.ac.il), which uses phylogenetic analysis and an atomistic potential, including an empirical lipophilicity scale, to improve native-state energy. As a stringent test, we apply mPROSS to the Kv1.2-Kv2.1 paddle chimera voltage-gated potassium channel. Four designs, encoding 9-26 mutations relative to the parental channel, were functional and maintained potassium-selective permeation and voltage dependence in Xenopus oocytes with up to 14-fold increase in whole-cell current densities. Additionally, single-channel recordings reveal no significant change in the channel-opening probability nor in unitary conductance, indicating that functional expression levels increase without impacting the activity profile of individual channels. Our results suggest that the expression levels of other dynamic channels and receptors may be enhanced through one-shot design calculations.

In-silico screening and identification of glycomimetic as potential human sodium-glucose co-transporter 2 inhibitor

Sodium-glucose co-transporter 2 (SGLT2) is one of the important targets against type II diabetes mellitus. A typical SGLT2 inhibitor acts by inhibiting glucose reabsorption, thus lowering the blood glucose level. Unlike SGLT1, SGLT2 is responsible for almost 90% glucose reabsorption from glomerular filtrate. The current SGLT2 inhibitors include gliflozins, often prescribed as second or third-line agents in diabetes mellitus. The SGLT2 inhibitors also benefit patients with heart and kidney disease. Due to instability issues with the natural O-aryl glycoside analogues C-glycoside analogues were developed and showed improved stability. Despite enhanced bioavailability and selectivity of newer derivatives, some serious side effects are associated with gliflozin analogues. At the present study, we applied in-silico approaches to find new glycomimetic compounds as potent SGLT2 inhibitors that could show improvement in side effects associated with current analogues. This work applied both ligand-based and structure-based drug approaches to find potential compounds. We developed a 3D-QSAR method to screen potential inhibitors from a library of ten thousand compounds and performed docking studies. The compounds were ranked based on predicted pIC50 and docking score. An initial screening of five thousand compounds was conducted, and the subsequently selected top 12 compounds were based on binding free energy calculations. These selected compounds were subjected to molecular dynamics (MD) simulations. Remarkably, our simulations identified nine compounds that exhibited significant and sustained binding affinity compared to the co-crystallized Empagliflozin. Collectively, considering the anticipated pharmacokinetic profiles and toxicity assessments, several of these compounds emerged as promising candidates for further in-depth evaluation.