Structural diversity of leukotriene G-protein coupled receptors.

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

Regression-Based Active Learning for Accessible Acceleration of Ultra-Large Library Docking.

Structure-based drug discovery is a process for both hit finding and optimization that relies on a validated three-dimensional model of a target biomolecule, used to rationalize the structure-function relationship for this particular target. An ultralarge virtual screening approach has emerged recently for rapid discovery of high-affinity hit compounds, but it requires substantial computational resources. This study shows that active learning with simple linear regression models can accelerate virtual screening, retrieving up to 90% of the top-1% of the docking hit list after docking just 10% of the ligands. The results demonstrate that it is unnecessary to use complex models, such as deep learning approaches, to predict the imprecise results of ligand docking with a low sampling depth. Furthermore, we explore active learning meta-parameters and find that constant batch size models with a simple ensembling method provide the best ligand retrieval rate. Finally, our approach is validated on the ultralarge size virtual screening data set, retrieving 70% of the top-0.05% of ligands after screening only 2% of the library. Altogether, this work provides a computationally accessible approach for accelerated virtual screening that can serve as a blueprint for the future design of low-compute agents for exploration of the chemical space via large-scale accelerated docking. With recent breakthroughs in protein structure prediction, this method can significantly increase accessibility for the academic community and aid in the rapid discovery of high-affinity hit compounds for various targets.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: A Structural Biology Perspective.

G protein-coupled receptors (GPCRs) play a key role in the transduction of extracellular signals to cells and regulation of many biological processes, which makes these membrane proteins one of the most important targets for pharmacological agents. A significant increase in the number of resolved atomic structures of GPCRs has opened the possibility of developing pharmaceuticals targeting these receptors via structure-based drug design (SBDD). SBDD employs information on the structure of receptor-ligand complexes to search for selective ligands without the need for an extensive high-throughput experimental ligand screening and can significantly expand the chemical space for ligand search. In this review, we describe the process of deciphering GPCR structures using X-ray diffraction analysis and cryoelectron microscopy as an important stage in the rational design of drugs targeting this receptor class. Our main goal was to present modern developments and key features of experimental methods used in SBDD of GPCR-targeting agents to a wide range of specialists.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: Ligand Search and Screening.

G protein-coupled receptors (GPCRs) are transmembrane proteins that participate in many physiological processes and represent major pharmacological targets. Recent advances in structural biology of GPCRs have enabled the development of drugs based on the receptor structure (structure-based drug design, SBDD). SBDD utilizes information about the receptor-ligand complex to search for suitable compounds, thus expanding the chemical space of possible receptor ligands without the need for experimental screening. The review describes the use of structure-based virtual screening (SBVS) for GPCR ligands and approaches for the functional testing of potential drug compounds, as well as discusses recent advances and successful examples in the application of SBDD for the identification of GPCR ligands.

Sub-millisecond conformational dynamics of the A2A adenosine receptor revealed by single-molecule FRET.

The complex pharmacology of G-protein-coupled receptors (GPCRs) is defined by their multi-state conformational dynamics. Single-molecule Förster Resonance Energy Transfer (smFRET) is well suited to quantify dynamics for individual protein molecules; however, its application to GPCRs is challenging. Therefore, smFRET has been limited to studies of inter-receptor interactions in cellular membranes and receptors in detergent environments. Here, we performed smFRET experiments on functionally active human A2A adenosine receptor (A2AAR) molecules embedded in freely diffusing lipid nanodiscs to study their intramolecular conformational dynamics. We propose a dynamic model of A2AAR activation that involves a slow (>2 ms) exchange between the active-like and inactive-like conformations in both apo and antagonist-bound A2AAR, explaining the receptor's constitutive activity. For the agonist-bound A2AAR, we detected faster (390 ± 80 µs) ligand efficacy-dependent dynamics. Our work establishes a general smFRET platform for GPCR investigations that can potentially be used for drug screening and/or mechanism-of-action studies.

Functional GPCR Expression in Eukaryotic LEXSY System.

G protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins in the human genome, and represent one of the most important classes of drug targets. Their structural studies facilitate rational drug discovery. However, atomic structures of only about 20% of human GPCRs have been solved to date. Recombinant production of GPCRs for structural studies at a large scale is challenging due to their low expression levels and stability. Therefore, in this study, we explored the efficacy of the eukaryotic system LEXSY (Leishmania tarentolae) for GPCR production. We selected the human A2A adenosine receptor (A2AAR), as a model protein, expressed it in LEXSY, purified it, and compared with the same receptor produced in insect cells, which is the most popular expression system for structural studies of GPCRs. The A2AAR purified from both expression systems showed similar purity, stability, ligand-induced conformational changes and structural dynamics, with a remarkably higher protein yield in the case of LEXSY expression. Overall, our results suggest that LEXSY is a promising platform for large-scale production of GPCRs for structural studies.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors.

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P1-5. S1P5 is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P5 receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P5-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

The Highly Conservative Cysteine of Oncomodulin as a Feasible Redox Sensor.

Oncomodulin (Ocm), or parvalbumin ß, is an 11-12 kDa Ca2+-binding protein found inside and outside of vertebrate cells, which regulates numerous processes via poorly understood mechanisms. Ocm consists of two active Ca2+-specific domains of the EF-hand type ("helix-loop-helix" motif), covered by an EF-hand domain with inactive EF-hand loop, which contains a highly conservative cysteine with unknown function. In this study, we have explored peculiarities of the microenvironment of the conservative Cys18 of recombinant rat Ocm (rWT Ocm), redox properties of this residue, and structural/functional sensitivity of rWT Ocm to the homologous C18S substitution. We have found that pKa of the Cys18 thiol lays beyond the physiological pH range. The measurement of redox dependence of rWT Ocm thiol-disulfide equilibrium (glutathione redox pair) showed that redox potential of Cys18 for the metal-free and Ca2+-loaded protein is of -168 mV and -176 mV, respectively. Therefore, the conservative thiol of rWT Ocm is prone to disulfide dimerization under physiological redox conditions. The C18S substitution drastically reduces α-helices content of the metal-free and Mg2+-bound Ocm, increases solvent accessibility of its hydrophobic residues, eliminates the cooperative thermal transition in the apo-protein, suppresses Ca2+/Mg2+ affinity of the EF site, and accelerates Ca2+ dissociation from Ocm. The distinct structural and functional consequences of the minor structural modification of Cys18 indicate its possible redox sensory function. Since some other EF-hand proteins also contain a conservative redox-sensitive cysteine located in an inactive EF-hand loop, it is reasonable to suggest that in the course of evolution, some of the EF-hands attained redox sensitivity at the expense of the loss of their Ca2+ affinity.

Small-wedge synchrotron and serial XFEL datasets for Cysteinyl leukotriene GPCRs.

Structural studies of challenging targets such as G protein-coupled receptors (GPCRs) have accelerated during the last several years due to the development of new approaches, including small-wedge and serial crystallography. Here, we describe the deposition of seven datasets consisting of X-ray diffraction images acquired from lipidic cubic phase (LCP) grown microcrystals of two human GPCRs, Cysteinyl leukotriene receptors 1 and 2 (CysLT1R and CysLT2R), in complex with various antagonists. Five datasets were collected using small-wedge synchrotron crystallography (SWSX) at the European Synchrotron Radiation Facility with multiple crystals under cryo-conditions. Two datasets were collected using X-ray free electron laser (XFEL) serial femtosecond crystallography (SFX) at the Linac Coherent Light Source, with microcrystals delivered at room temperature into the beam within LCP matrix by a viscous media microextrusion injector. All seven datasets have been deposited in the open-access databases Zenodo and CXIDB. Here, we describe sample preparation and annotate crystallization conditions for each partial and full datasets. We also document full processing pipelines and provide wrapper scripts for SWSX and SFX data processing.

Author Correction: Small-wedge synchrotron and serial XFEL datasets for Cysteinyl leukotriene GPCRs.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-020-00759-w.

Effects of his-tags on physical properties of parvalbumins.

A comparative study of His-tagged and non-tagged rat ß-parvalbumin (rWT ß-PA), calcium binding protein with the EF-hand calcium binding domains, has been carried out. The attachment of His-tag increases α-helical content and decreases ß-sheets and ß-turns content of the metal free form (apo-state) of ß-PA. In contrast to this, the attachment of His-tag decreases α-helical content by more than 10% and increases contents of ß-sheets and ß-turns of the Ca2+-loaded state. According to the dynamic light scattering analysis, apo-state of His-tagged rat ß-PA seems to be less compact compared with the apo-state of non-tagged rat ß-PA. Surprisingly, the attachment of His-tag practically does not change mean hydrodynamic radius of Ca2+-loaded rat ß-PA. The attachment of His-tag shifts thermal denaturation peaks of both apo- and Ca2+-loaded states of rat ß-PA towards higher temperatures by 3-4 °C and slightly decreases its Ca2+ affinity. These results should be taken into consideration in the use of His-tagged parvalbumins.

Comprehensive analysis of the roles of 'black' and 'gray' clusters in structure and function of rat ß-parvalbumin.

Recently we found two highly conserved structural motifs in the proteins of the EF-hand calcium binding protein family. These motifs provide a supporting scaffold for the Ca2+ binding loops and contribute to the hydrophobic core of the EF-hand domain. Each structural motif forms a cluster of three amino acids called cluster I ('black' cluster) and cluster II ('grey' cluster). Cluster I is much more conserved and mostly incorporates aromatic amino acids. In contrast, cluster II includes a mix of aromatic, hydrophobic, and polar amino acids. The 'black' and 'gray' clusters in rat ß-parvalbumin consist of F48, A100, F103 and G61, L64, M87, respectively. In the present work, we sequentially substituted these amino acids residues by Ala, except Ala100, which was substituted by Val. Physical properties of the mutants were studied by circular dichroism, scanning calorimetry, dynamic light scattering, chemical crosslinking, and fluorescent probe methods. The Ca2+ and Mg2+ binding affinities of these mutants were evaluated by intrinsic fluorescence and equilibrium dialysis methods. In spite of a rather complicated pattern of contributions of separate amino acid residues of the 'black' and 'gray' clusters into maintenance of rat ß-parvalbumin structural and functional status, the alanine substitutions in the cluster I cause noticeably more pronounced changes in various structural parameters of proteins, such as hydrodynamic radius of apo-form, thermal stability of Ca2+/Mg2+-loaded forms, and total energy of Ca2+ binding in comparison with the changes caused by amino acid substitutions in the cluster II. These findings were further supported by the outputs of computational analysis of the effects of these mutations on the intrinsic disorder predisposition of rat ß-parvalbumin, which also indicated that local intrinsic disorder propensities and the overall levels of predicted disorder were strongly affected by mutations in the cluster I, whereas mutations in cluster II had less pronounced effects. These results demonstrate that amino acids of the cluster I provide more essential contribution to the maintenance of structuraland functional properties of the protein in comparison with the residues of the cluster II.

In search for globally disordered apo-parvalbumins: Case of parvalbumin ß-1 from coho salmon.

Parvalbumin (PA) is a classical EF-hand calcium-binding protein of muscle, neuronal, and other tissues, and a major fish allergen. Although certain apo-PAs lack tertiary structure, functional implications of that feature and its structural prerequisites remain unclear. In a search for unstable PAs, we probed conformational stability of parvalbumin ß-1 from coho salmon (csPA), a cold water fish species, using circular dichroism, scanning calorimetry, hydrophobic probe fluorescence, limited proteolysis, chemical crosslinking and dynamic light scattering techniques. Apo-csPA is shown to be mainly monomeric protein with markedly disorganized secondary structure and lack of rigid tertiary structure. Examination of per-residue propensity for intrinsic disorder in the PA groups with either folded or unfolded apo-form using the average PONDR® VSL2P profiles revealed that the N-terminal region that includes α-helix A, AB-loop and N-terminal half of α-helix B is predicted to be less ordered in PAs with disordered apo-state. Application of the structural criteria developed for discrimination of disordered PAs indicate that the latter comprise about 16-19% of all PAs. We show that structural instability of apo-ß-PA serves as a hallmark of elevated calcium affinity of the protein. Therefore, the successful predictions of unstable apo-PAs might facilitate search for PAs with maximal calcium affinity and possibly serving as calcium sensors.

Interleukin-11 binds specific EF-hand proteins via their conserved structural motifs.

Interleukin-11 (IL-11) is a hematopoietic cytokine engaged in numerous biological processes and validated as a target for treatment of various cancers. IL-11 contains intrinsically disordered regions that might recognize multiple targets. Recently we found that aside from IL-11RA and gp130 receptors, IL-11 interacts with calcium sensor protein S100P. Strict calcium dependence of this interaction suggests a possibility of IL-11 interaction with other calcium sensor proteins. Here we probed specificity of IL-11 to calcium-binding proteins of various types: calcium sensors of the EF-hand family (calmodulin, S100B and neuronal calcium sensors: recoverin, NCS-1, GCAP-1, GCAP-2), calcium buffers of the EF-hand family (S100G, oncomodulin), and a non-EF-hand calcium buffer (α-lactalbumin). A specific subset of the calcium sensor proteins (calmodulin, S100B, NCS-1, GCAP-1/2) exhibits metal-dependent binding of IL-11 with dissociation constants of 1-19 µM. These proteins share several amino acid residues belonging to conservative structural motifs of the EF-hand proteins, 'black' and 'gray' clusters. Replacements of the respective S100P residues by alanine drastically decrease its affinity to IL-11, suggesting their involvement into the association process. Secondary structure and accessibility of the hinge region of the EF-hand proteins studied are predicted to control specificity and selectivity of their binding to IL-11. The IL-11 interaction with the EF-hand proteins is expected to occur under numerous pathological conditions, accompanied by disintegration of plasma membrane and efflux of cellular components into the extracellular milieu.