Fine-tuning activation specificity of G-protein-coupled receptors via automated path searching.

Physics-based simulation methods can grant atomistic insights into the molecular origin of the function of biomolecules. However, the potential of such approaches has been hindered by their low efficiency, including in the design of selective agonists where simulations of myriad protein-ligand combinations are necessary. Here, we describe an automated input-free path searching protocol that offers (within 14 d using Graphics Processing Unit servers) a minimum free energy path (MFEP) defined in high-dimension configurational space for activating sphingosine-1-phosphate receptors (S1PRs) by arbitrary ligands. The free energy distributions along the MFEP for four distinct ligands and three S1PRs reached a remarkable agreement with Bioluminescence Resonance Energy Transfer (BRET) measurements of G-protein dissociation. In particular, the revealed transition state structures pointed out toward two S1PR3 residues F263/I284, that dictate the preference of existing agonists CBP307 and BAF312 on S1PR1/5. Swapping these residues between S1PR1 and S1PR3 reversed their response to the two agonists in BRET assays. These results inspired us to design improved agonists with both strong polar head and bulky hydrophobic tail for higher selectivity on S1PR1. Through merely three in silico iterations, our tool predicted a unique compound scaffold. BRET assays confirmed that both chiral forms activate S1PR1 at nanomolar concentration, 1 to 2 orders of magnitude less than those for S1PR3/5. Collectively, these results signify the promise of our approach in fine agonist design for G-protein-coupled receptors.

Structural insights into sphingosine-1-phosphate receptor activation.

As a critical sphingolipid metabolite, sphingosine-1-phosphate (S1P) plays an essential role in immune and vascular systems. There are five S1P receptors, designated as S1PR1 to S1PR5, encoded in the human genome, and their activities are governed by endogenous S1P, lipid-like S1P mimics, or nonlipid-like therapeutic molecules. Among S1PRs, S1PR1 stands out due to its nonredundant functions, such as the egress of T and B cells from the thymus and secondary lymphoid tissues, making it a potential therapeutic target. However, the structural basis of S1PR1 activation and regulation by various agonists remains unclear. Here, we report four atomic resolution cryo-electron microscopy (cryo-EM) structures of Gi-coupled human S1PR1 complexes: bound to endogenous agonist d18:1 S1P, benchmark lipid-like S1P mimic phosphorylated Fingolimod [(S)-FTY720-P], or nonlipid-like therapeutic molecule CBP-307 in two binding modes. Our results revealed the similarities and differences of activation of S1PR1 through distinct ligands binding to the amphiphilic orthosteric pocket. We also proposed a two-step "shallow to deep" transition process of CBP-307 for S1PR1 activation. Both binding modes of CBP-307 could activate S1PR1, but from shallow to deep transition may trigger the rotation of the N-terminal helix of Gαi and further stabilize the complex by increasing the Gαi interaction with the cell membrane. We combine with extensive biochemical analysis and molecular dynamic simulations to suggest key steps of S1P binding and receptor activation. The above results decipher the common feature of the S1PR1 agonist recognition and activation mechanism and will firmly promote the development of therapeutics targeting S1PRs.

Differential Regulation of Microglial Activation in Response to Different Degree of Ischemia.

Microglia are primary immune cells within the brain and are rapidly activated after cerebral ischemia. The degree of microglial activation is closely associated with the severity of ischemia. However, it remains largely unclear how microglial activation is differentially regulated in response to a different degree of ischemia. In this study, we used a bilateral common carotid artery ligation (BCAL) model and induced different degrees of ischemia by varying the duration of ligation to investigate the microglial response in CX3CR1GFP/+ mice. Confocal microscopy, immunofluorescence staining, RNA sequencing, and qRT-PCR were used to evaluate the de-ramification, proliferation, and differential gene expression associated with microglial activation. Our results showed that 30 min of ischemia induced rapid de-ramification of microglia but did not have significant influence on the microglial density. In contrast, 60 min of ischemia led to a significant decrease in microglial density and more pronounced de-ramification of microglial processes. Importantly, 30 min of ischemia did not induce proliferation of microglia, but 60 min of ischemia led to a marked increase in the density of proliferative microglia. Further analysis utilized transcriptome sequencing showed that microglial activation is differentially regulated in response to different degrees of ischemia. A total of 1,097 genes were differentially regulated after 60 min of ischemia, but only 68 genes were differentially regulated after 30 min of ischemia. Pathway enrichment analysis showed that apoptosis, cell mitosis, immune receptor activity and inflammatory-related pathways were highly regulated after 60 min of ischemia compared to 30 min of ischemia. Multiple microglia-related genes such as Cxcl10, Tlr7, Cd86, Tnfrsf1a, Nfkbia, Tgfb1, Ccl2 and Il-6, were upregulated with prolonged ischemia. Pharmacological inhibition of CSF1 receptor demonstrated that CSF1R signaling pathway contributed to microglial proliferation. Together, these results suggest that the proliferation of microglia is gated by the duration of ischemia and microglia were differentially activated in responding to different degrees of ischemia.

The complete mitochondrial genome of Great tit Parus major (Aves, Passeriformes, Paridae).

The complete mitochondrial genome sequence of Great tit Parus major was sequenced used polymerase chain reaction (PCR), long-and-accurate PCR and directly sequencing by primer walking. The Genbank accession was KP137624. The entire mitochondrial genome of P. major is a circular molecule of 16,776 bp in length and the content of A, T, C and G were 29.68%, 22.63%, 33.56% and 14.13%, respectively. The complete mitochondrial genome of P. major contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, plus 1 control regions and was similar to most of the other Aves birds in gene arrangement and composition. The complete mitochondrial genome of P. major could provide a useful data for resolving phylogenetic relationship problems related to Parus and P. major subspecies complex.