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

A shared mechanism for TNP-ATP recognition by members of the P2X receptor family

P2X receptors (P2X1-7) are non-selective cation channels involved in many physiological activities such as synaptic transmission, immunological modulation, and cardiovascular function. These receptors share a conserved mechanism to sense extracellular ATP. TNP-ATP is an ATP derivative acting as a nonselective competitive P2X antagonist. Understanding how it occupies the orthosteric site in the absence of agonism may help reveal the key allostery during P2X gating. However, TNP-ATP/P2X complexes (TNP-ATP/human P2X3 (hP2X3) and TNP-ATP/chicken P2X7 (ckP2X7)) with distinct conformations and different mechanisms of action have been proposed. Whether these represent species and subtype variations or experimental differences remains unclear. Here, we show that a common mechanism of TNP-ATP recognition exists for the P2X family members by combining enhanced conformation sampling, engineered disulfide bond analysis, and covalent occupancy. In this model, the polar triphosphate moiety of TNP-ATP interacts with the orthosteric site, while its TNP-moiety is deeply embedded in the head and dorsal fin (DF) interface, creating a restrictive allostery in these two domains that results in a partly enlarged yet ion-impermeable pore. Similar results were obtained from multiple P2X subtypes of different species, including ckP2X7, hP2X3, rat P2X2 (rP2X2), and human P2X1 (hP2X1). Thus, TNP-ATP uses a common mechanism for P2X recognition and modulation by restricting the movements of the head and DF domains which are essential for P2X activation. This knowledge is applicable to the development of new P2X inhibitors.

Mechanistic insights into P-glycoprotein ligand transport and inhibition revealed by enhanced molecular dynamics simulations

P-glycoprotein (P-gp) plays a crucial role in cellular detoxification and drug efflux processes, transitioning between inward-facing (IF) open, occluded, and outward-facing (OF) states to facilitate substrate transport. Its role is critical in cancer therapy, where P-gp contributes to the multidrug resistance phenotype. In our study, classical and enhanced molecular dynamics (MD) simulations were conducted to dissect the structural and functional features of the P-gp conformational states. Our advanced MD simulations, including kinetically excited targeted MD (ketMD) and adiabatic biasing MD (ABMD), provided deeper insights into state transition and translocation mechanisms. Our findings suggest that the unkinking of TM4 and TM10 helices is a prerequisite for correctly achieving the outward conformation. Simulations of the IF-occluded conformations, characterized by kinked TM4 and TM10 helices, consistently demonstrated altered communication between the transmembrane domains (TMDs) and nucleotide binding domain 2 (NBD2), suggesting the implication of this interface in inhibiting P-gp's efflux function. A particular emphasis was placed on the unstructured linker segment connecting the NBD1 to TMD2 and its role in the transporter's dynamics. With the linker present, we specifically noticed a potential entrance of cholesterol (CHOL) through the TM4-TM6 portal, shedding light on crucial residues involved in accommodating CHOL. We therefore suggest that this entry mechanism could be employed for some P-gp substrates or inhibitors. Our results provide critical data for understanding P-gp functioning and developing new P-gp inhibitors for establishing more effective strategies against multidrug resistance.

Structural insights into the allosteric effects of the antiepileptic drug topiramate on the CaV2.3 channel

The R-type voltage-gated calcium channel CaV2.3 is predominantly located in the presynapse and is implicated in distinct types of epileptic seizures. It has consequently emerged as a molecular target in seizure treatment. Here, we determined the cryo-EM structure of the CaV2.3-α2δ1-β1 complex in the topiramate-bound state at a 3.0 Å resolution. We provide a snapshot of the binding site of topiramate, a widely prescribed antiepileptic drug, on a voltage-gated ion channel. The binding site is located at an intracellular juxtamembrane hydrophilic cavity. Further structural analysis revealed that topiramate may allosterically facilitate channel inactivation. These findings provide fundamental insights into the mechanism underlying the inhibitory effect of topiramate on CaV and NaV channels, elucidating a previously unseen modulator binding site and thus pointing toward a route for the development of new drugs.

5-HT6 receptor neutral antagonists protect astrocytes: A lesson from 2-phenylpyrrole derivatives

The serotonin type 6 receptor (5-HT6R) displays a strong constitutive activity, suggesting it participates largely in the physiological and pathological processes controlled by the receptor. The active states of 5-HT6R engage particular signal transduction pathways that lead to different biological responses. In this study, we present the development of 5-HT6R neutral antagonists at Gs signaling built upon the 2-phenylpyrrole scaffold. Using molecular dynamics simulations, we outline the relationship between the exposure of the basic center of the molecules and their ability to target the agonist-activated state of the receptor. Our study identifies compound 30 as a potent and selective neutral antagonist at 5-HT6R-operated Gs signaling. Furthermore, we demonstrate the cytoprotective effects of 30 and structurally diverse 5-HT6R neutral antagonists at Gs signaling in C8-D1A cells and human astrocytes exposed to rotenone. This effect is not observed for 5-HT6R agonists or inverse agonists. In light of these findings, we propose compound 30 as a valuable molecular probe to study the biological effects associated with the agonist-activated state of 5-HT6R and provide insight into the glioprotective properties of 5-HT6R neutral antagonists at Gs signaling.

The C-terminus of Rad is required for membrane localization and L-type calcium channel regulation

L-type CaV1.2 current (ICa,L) links electrical excitation to contraction in cardiac myocytes. ICa,L is tightly regulated to control cardiac output. Rad is a Ras-related, monomeric protein that binds to L-type calcium channel β subunits (CaVβ) to promote inhibition of ICa,L. In addition to CaVβ interaction conferred by the Rad core motif, the highly conserved Rad C-terminus can direct membrane association in vitro and inhibition of ICa,L in immortalized cell lines. In this work, we test the hypothesis that in cardiomyocytes the polybasic C-terminus of Rad confers t-tubular localization, and that membrane targeting is required for Rad-dependent ICa,L regulation. We introduced a 3xFlag epitope to the N-terminus of the endogenous mouse Rrad gene to facilitate analysis of subcellular localization. Full-length 3xFlag-Rad (Flag-Rad) mice were compared with a second transgenic mouse model, in which the extended polybasic C-termini of 3xFlag-Rad was truncated at alanine 277 (Flag-RadΔCT). Ventricular cardiomyocytes were isolated for anti-Flag-Rad immunocytochemistry and ex vivo electrophysiology. Full-length Flag-Rad showed a repeating t-tubular pattern whereas Flag-RadΔCT failed to display membrane association. ICa,L in Flag-RadΔCT cardiomyocytes showed a hyperpolarized activation midpoint and an increase in maximal conductance. Additionally, current decay was faster in Flag-RadΔCT cells. Myocardial ICa,L in a Rad C-terminal deletion model phenocopies ICa,L modulated in response to β-AR stimulation. Mechanistically, the polybasic Rad C-terminus confers CaV1.2 regulation via membrane association. Interfering with Rad membrane association constitutes a specific target for boosting heart function as a treatment for heart failure with reduced ejection fraction.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Isolation, Structure Elucidation, and Biological Activity of the Selective TACR2 Antagonist Tumonolide and its Aldehyde from a Marine Cyanobacterium

The macrocyclic tumonolide (1) with enamide functionality and the linear tumonolide aldehyde (2) are new interconverting natural products from a marine cyanobacterium with a peptide-polyketide skeleton, representing a hybrid of apratoxins and palmyrolides or laingolides. The planar structures were established by NMR and mass spectrometry. The relative configuration of the stereogenically-rich apratoxin-like polyketide portion was determined using J-based configuration analysis. The absolute configuration of tumonolide (1) was determined by chiral analysis of the amino acid units and computational methods, followed by NMR chemical shift and ECD spectrum prediction, indicating all-R configuration for the polyketide portion, as in palmyrolide A and contrary to the all-S configuration in apratoxins. Functional screening against a panel of 168 GPCR targets revealed tumonolide (1) as a selective antagonist of TACR2 with an IC50 of 7.0 μM, closely correlating with binding affinity. Molecular docking studies established the binding mode and rationalized the selectivity for TACR2 over TACR1 and TACR3. RNA sequencing upon treatment of HCT116 colorectal cancer cells demonstrated activation of the pulmonary fibrosis idiopathic signaling pathway and the insulin secretion signaling pathway at 20 μM, indicating its potential to modulate these pathways.

Enhancement of high-density lipoprotein-associated protease inhibitor activity prevents atherosclerosis progression

BACKGROUND AND AIMS

Inflammatory cells within atherosclerotic lesions secrete proteolytic enzymes that contribute to lesion progression and destabilization, increasing the risk for an acute cardiovascular event. Elastase is a serine protease, secreted by macrophages and neutrophils, that may contribute to the development of unstable plaque. We previously reported interaction of endogenous protease-inhibitor proteins with high-density lipoprotein (HDL), including alpha-1-antitrypsin, an inhibitor of elastase. These findings support a potential role for HDL as a modulator of protease activity. In this study, we test the hypothesis that enhancement of HDL-associated elastase inhibitor activity is protective against atherosclerotic lesion progression.

METHODS

We designed an HDL-targeting protease inhibitor (HTPI) that binds to HDL and confers elastase inhibitor activity. Lipoprotein binding and the impact of HTPI on atherosclerosis were examined using mouse models. Histology and immunofluorescence staining of aortic root sections were used to examine the impact of HTPI on lesion morphology and inflammatory features.

RESULTS

HTPI is a small (1.6 kDa) peptide with an elastase inhibitor domain, a soluble linker, and an HDL-targeting domain. When incubated with human plasma ex vivo, HTPI predominantly binds to HDL. Intravenous administration of HTPI to mice resulted in its binding to plasma HDL and increased elastase inhibitor activity on isolated HDL. Accumulation of HTPI within plaque was observed after administration to Apoe-/- mice. To examine the effect of HTPI treatment on atherosclerosis, prevention and progression studies were performed using Ldlr-/- mice fed Western diet. In both study designs, HTPI-treated mice had reduced lipid deposition in plaque.

CONCLUSIONS

These data support the hypothesis that HDL-associated anti-elastase activity can improve the atheroprotective potential of HDL and highlight the potential utility of HDL enrichment with anti-protease activity as an approach for stabilization of atherosclerotic lesions.

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.

They all rock: A systematic comparison of conformational movements in LeuT-fold transporters

Many membrane transporters share the LeuT fold-two five-helix repeats inverted across the membrane plane. Despite hundreds of structures, whether distinct conformational mechanisms are supported by the LeuT fold has not been systematically determined. After annotating published LeuT-fold structures, we analyzed distance difference matrices (DDMs) for nine proteins with multiple available conformations. We identified rigid bodies and relative movements of transmembrane helices (TMs) during distinct steps of the transport cycle. In all transporters, the bundle (first two TMs of each repeat) rotates relative to the hash (third and fourth TMs). Motions of the arms (fifth TM) to close or open the intracellular and outer vestibules are common, as is a TM1a swing, with notable variations in the opening-closing motions of the outer vestibule. Our analyses suggest that LeuT-fold transporters layer distinct motions on a common bundle-hash rock and demonstrate that systematic analyses can provide new insights into large structural datasets.

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.

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.

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.

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.

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.

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.

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.

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.

Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work, we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi, and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.