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.

Two distinct mechanisms of flavoprotein spectral tuning revealed by low-temperature and time-dependent spectroscopy.

Flavins such as flavin mononucleotide or flavin adenine dinucleotide are bound by diverse proteins, yet have very similar spectra when in the oxidized state. Recently, we developed new variants of flavin-binding protein CagFbFP exhibiting notable blue (Q148V) or red (I52V A85Q) shifts of fluorescence emission maxima. Here, we use time-resolved and low-temperature spectroscopy to show that whereas the chromophore environment is static in Q148V, an additional protein-flavin hydrogen bond is formed upon photoexcitation in the I52V A85Q variant. Consequently, in Q148V, excitation, emission, and phosphorescence spectra are shifted, whereas in I52V A85Q, excitation and low-temperature phosphorescence spectra are relatively unchanged, while emission spectrum is altered. We also determine the x-ray structures of the two variants to reveal the flavin environment and complement the spectroscopy data. Our findings illustrate two distinct color-tuning mechanisms of flavin-binding proteins and could be helpful for the engineering of new variants with improved optical properties.

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.

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 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.

The Role of Tyr-His-Trp Triad and Water Molecule Near the N1-Atom of 2-Hydroperoxycoelenterazine in Bioluminescence of Hydromedusan Photoproteins: Structural and Mutagenesis Study.

Hydromedusan photoproteins responsible for the bioluminescence of a variety of marine jellyfish and hydroids are a unique biochemical system recognized as a stable enzyme-substrate complex consisting of apoprotein and preoxygenated coelenterazine, which is tightly bound in the protein inner cavity. The binding of calcium ions to the photoprotein molecule is only required to initiate the light emission reaction. Although numerous experimental and theoretical studies on the bioluminescence of these photoproteins were performed, many features of their functioning are yet unclear. In particular, which ionic state of dioxetanone intermediate decomposes to yield a coelenteramide in an excited state and the role of the water molecule residing in a proximity to the N1 atom of 2-hydroperoxycoelenterazine in the bioluminescence reaction are still under discussion. With the aim to elucidate the function of this water molecule as well as to pinpoint the amino acid residues presumably involved in the protonation of the primarily formed dioxetanone anion, we constructed a set of single and double obelin and aequorin mutants with substitutions of His, Trp, Tyr, and Ser to residues with different properties of side chains and investigated their bioluminescence properties (specific activity, bioluminescence spectra, stopped-flow kinetics, and fluorescence spectra of Ca2+-discharged photoproteins). Moreover, we determined the spatial structure of the obelin mutant with a substitution of His64, the key residue of the presumable proton transfer, to Phe. On the ground of the bioluminescence properties of the obelin and aequorin mutants as well as the spatial structures of the obelin mutants with the replacements of His64 and Tyr138, the conclusion was made that, in fact, His residue of the Tyr-His-Trp triad and the water molecule perform the "catalytic function" by transferring the proton from solvent to the dioxetanone anion to generate its neutral ionic state in complex with water, as only the decomposition of this form of dioxetanone can provide the highest light output in the light-emitting reaction of the hydromedusan photoproteins.

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.

Structural insights into thrombolytic activity of destabilase from medicinal leech.

Destabilase from the medical leech Hirudo medicinalis belongs to the family of i-type lysozymes. It has two different enzymatic activities: microbial cell walls destruction (muramidase activity), and dissolution of the stabilized fibrin (isopeptidase activity). Both activities are known to be inhibited by sodium chloride at near physiological concentrations, but the structural basis remains unknown. Here we present two crystal structures of destabilase, including a 1.1 Å-resolution structure in complex with sodium ion. Our structures reveal the location of sodium ion between Glu34/Asp46 residues, which were previously recognized as a glycosidase active site. While sodium coordination with these amino acids may explain inhibition of the muramidase activity, its influence on previously suggested Ser49/Lys58 isopeptidase activity dyad is unclear. We revise the Ser49/Lys58 hypothesis and compare sequences of i-type lysozymes with confirmed destabilase activity. We suggest that the general base for the isopeptidase activity is His112 rather than Lys58. pKa calculations of these amino acids, assessed through the 1 µs molecular dynamics simulation, confirm the hypothesis. Our findings highlight the ambiguity of destabilase catalytic residues identification and build foundations for further research of structure-activity relationship of isopeptidase activity as well as structure-based protein design for potential anticoagulant drug development.

Protein Design Strategies for the Structural-Functional Studies of G Protein-Coupled Receptors.

G protein-coupled receptors (GPCRs) are an important family of membrane proteins responsible for many physiological functions in human body. High resolution GPCR structures are required to understand their molecular mechanisms and perform rational drug design, as GPCRs play a crucial role in a variety of diseases. That is difficult to obtain for the wild-type proteins because of their low stability. In this review, we discuss how this problem can be solved by using protein design strategies developed to obtain homogeneous stabilized GPCR samples for crystallization and cryoelectron microscopy.

Structural insights into the effects of glycerol on ligand binding to cytochrome P450.

New antitubercular drugs are vital due to the spread of resistant strains. Carbethoxyhexyl imidazole (CHImi) inhibits cytochrome P450 CYP124, which is a steroid-metabolizing enzyme that is important for the survival of Mycobacterium tuberculosis in macrophages. The available crystal structure of the CYP124-CHImi complex reveals two glycerol molecules in the active site. A 1.15 Šresolution crystal structure of the glycerol-free CYP124-CHimi complex reported here shows multiple conformations of CHImi and the CYP124 active site which were previously restricted by glycerol. Complementary molecular dynamics simulations show coherence of the ligand and enzyme conformations. Spectrophotometric titration confirmed the influence of glycerol on CHImi binding: the affinity decreases more than tenfold in glycerol-containing buffer. In addition, it also showed that glycerol has a similar effect on other azole and triazole CYP124 ligands. Together, these data show that glycerol may compromise structural-functional studies and impede rational drug-design campaigns.

Crystal structure of semi-synthetic obelin-v after calcium induced bioluminescence implies coelenteramine as the main reaction product.

Coelenterazine-v (CTZ-v), a synthetic vinylene-bridged π-extended derivative, is able to significantly alter bioluminescence spectra of different CTZ-dependent luciferases and photoproteins by shifting them towards longer wavelengths. However, Ca2+-regulated photoproteins activated with CTZ-v display very low bioluminescence activities that hampers its usage as a substrate of photoprotein bioluminescence. Here, we report the crystal structure of semi-synthetic Ca2+-discharged obelin-v bound with the reaction product determined at 2.1 Å resolution. Comparison of the crystal structure of Ca2+-discharged obelin-v with those of other obelins before and after bioluminescence reaction reveals no considerable changes in the overall structure. However, the drastic changes in CTZ-binding cavity are observed owing to the completely different reaction product, coelenteramine-v (CTM-v). Since CTM-v is certainly the main product of obelin-v bioluminescence and is considered to be a product of the "dark" pathway of dioxetanone intermediate decomposition, it explains the low bioluminescence activity of obelin and apparently of other photoproteins with CTZ-v.

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.

Structure-based rational design of an enhanced fluorogen-activating protein for fluorogens based on GFP chromophore.

"Fluorescence-Activating and absorption-Shifting Tag" (FAST) is a well-studied fluorogen-activating protein with high brightness and low size, able to activate a wide range of fluorogens. This makes FAST a promising target for both protein and fluorogen optimization. Here, we describe the structure-based rational design of the enhanced FAST mutants, optimized for the N871b fluorogen. Using the spatial structure of the FAST/N871b complex, NMR relaxation analysis, and computer simulations, we identify the mobile regions in the complex and suggest mutations that could stabilize both the protein and the ligand. Two of our mutants appear brighter than the wild-type FAST, and these mutants provide up to 35% enhancement for several other fluorogens of similar structure, both in vitro and in vivo. Analysis of the mutants by NMR reveals that brighter mutants demonstrate the highest stability and lowest length of intermolecular H-bonds. Computer simulations provide the structural basis for such stabilization.

Structural insights into 3Fe-4S ferredoxins diversity in M. tuberculosis highlighted by a first redox complex with P450.

Ferredoxins are small iron-sulfur proteins and key players in essential metabolic pathways. Among all types, 3Fe-4S ferredoxins are less studied mostly due to anaerobic requirements. Their complexes with cytochrome P450 redox partners have not been structurally characterized. In the present work, we solved the structures of both 3Fe-4S ferredoxins from M. tuberculosis-Fdx alone and the fusion FdxE-CYP143. Our SPR analysis demonstrated a high-affinity binding of FdxE to CYP143. According to SAXS data, the same complex is present in solution. The structure reveals extended multipoint interactions and the shape/charge complementarity of redox partners. Furthermore, FdxE binding induced conformational changes in CYP143 as evident from the solved CYP143 structure alone. The comparison of FdxE-CYP143 and modeled Fdx-CYP51 complexes further revealed the specificity of ferredoxins. Our results illuminate the diversity of electron transfer complexes for the production of different secondary metabolites.

A new species of the grenadier genus Coelorinchus from the seamounts in the southeastern Atlantic Ocean, with redefinition of C. vityazae (Teleostei: Gadiformes: Macrouridae).

Re-investigation of the grenadier Coelorinchus vityazae endemic to the West Wind Drift Islands Province reveals species-level differences between the populations from the southeastern Atlantic and southwestern Indian oceans. The southeastern Atlantic populations (from Discovery and Gough seamounts) are described as a new species, C. inventionis sp. nov., characterized by a moderately short snout (27.7-33.9 % HL, vs. 32.0-38.7 % in C. vityazae) tipped with short, weakly tripartite terminal scute (vs. triangular and sharply pointed); uniformly thick, unpigmented lips (vs. fleshy, partly blackish upper lip with lateral portions expanded at middle in C. vityazae); modally i+17 or i+18 pectoral-fin rays (vs. i+15 or i+16), and anal-fin rays conspicuously darkened distally (vs. uniformly and finely peppered). Statistically significant differences between these two species were found for 28 of 39 morphometric characters. The Discovery and Gough specimens show a drastic difference in size of the light organ, which may reflect an initial stage of speciation within C. inventionis sp. nov. Iwamoto & Graham's (2008) key to the species of the C. fasciatus group is modified for inclusion of the new species.

Accessing Mitochondrial Protein Import in Living Cells by Protein Microinjection.

Mitochondrial protein biogenesis relies almost exclusively on the expression of nuclear-encoded polypeptides. The current model postulates that most of these proteins have to be delivered to their final mitochondrial destination after their synthesis in the cytoplasm. However, the knowledge of this process remains limited due to the absence of proper experimental real-time approaches to study mitochondria in their native cellular environment. We developed a gentle microinjection procedure for fluorescent reporter proteins allowing a direct non-invasive study of protein transport in living cells. As a proof of principle, we visualized potential-dependent protein import into mitochondria inside intact cells in real-time. We validated that our approach does not distort mitochondrial morphology and preserves the endogenous expression system as well as mitochondrial protein translocation machinery. We observed that a release of nascent polypeptides chains from actively translating cellular ribosomes by puromycin strongly increased the import rate of the microinjected pre-protein. This suggests that a substantial amount of mitochondrial translocase complexes was involved in co-translational protein import of endogenously expressed pre-proteins. Our protein microinjection method opens new possibilities to study the role of mitochondrial protein import in cell models of various pathological conditions as well as aging processes.

NanoFAST: structure-based design of a small fluorogen-activating protein with only 98 amino acids.

One of the essential characteristics of any tag used in bioscience and medical applications is its size. The larger the label, the more it may affect the studied object, and the more it may distort its behavior. In this paper, using NMR spectroscopy and X-ray crystallography, we have studied the structure of fluorogen-activating protein FAST both in the apo form and in complex with the fluorogen. We showed that significant change in the protein occurs upon interaction with the ligand. While the protein is completely ordered in the complex, its apo form is characterized by higher mobility and disordering of its N-terminus. We used structural information to design the shortened FAST (which we named nanoFAST) by truncating 26 N-terminal residues. Thus, we created the shortest genetically encoded tag among all known fluorescent and fluorogen-activating proteins, which is composed of only 98 amino acids.

Structural insights into ligand recognition and activation of angiotensin receptors.

G protein-coupled angiotensin II receptors, AT1R and AT2R, are integral components of the renin-angiotensin system (RAS) that regulates blood pressure and fluid balance in humans. While AT1R is a well-established target of angiotensin receptor blockers (ARBs) for managing hypertension and a prime system for studying biased signaling, AT2R has been recognized as a promising target against neuropathic pain and lung fibrosis. In this review, we discuss how recent structural advances illuminate ligand-binding modes and subtype selectivity, shared and distinct features of the receptors, their transducer-coupling patterns, and downstream signaling responses. We also underscore the key ATR aspects that require further studies to fully appreciate the mechanistic framework that fine-tunes their cellular and physiological functions, providing untapped potential for drug discovery.

Metabolic Fate of Human Immunoactive Sterols in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) infection is among top ten causes of death worldwide, and the number of drug-resistant strains is increasing. The direct interception of human immune signaling molecules by Mtb remains elusive, limiting drug discovery. Oxysterols and secosteroids regulate both innate and adaptive immune responses. Here we report a functional, structural, and bioinformatics study of Mtb enzymes initiating cholesterol catabolism and demonstrated their interrelation with human immunity. We show that these enzymes metabolize human immune oxysterol messengers. Rv2266 - the most potent among them - can also metabolize vitamin D3 (VD3) derivatives. High-resolution structures show common patterns of sterols binding and reveal a site for oxidative attack during catalysis. Finally, we designed a compound that binds and inhibits three studied proteins. The compound shows activity against Mtb H37Rv residing in macrophages. Our findings contribute to molecular understanding of suppression of immunity and suggest that Mtb has its own transformation system resembling the human phase I drug-metabolizing system.

Structure-Based Virtual Screening of Ultra-Large Library Yields Potent Antagonists for a Lipid GPCR.

Cysteinyl leukotriene G protein-coupled receptors, CysLT1R and CysLT2R, regulate bronchoconstrictive and pro-inflammatory effects and play a key role in allergic disorders, cardiovascular diseases, and cancer. CysLT1R antagonists have been widely used to treat asthma disorders, while CysLT2R is a potential target against uveal melanoma. However, very few selective antagonist chemotypes for CysLT receptors are available, and the design of such ligands has proved to be challenging. To overcome this obstacle, we took advantage of recently solved crystal structures of CysLT receptors and an ultra-large Enamine REAL library, representing a chemical space of 680 M readily available compounds. Virtual ligand screening employed 4D docking models comprising crystal structures of CysLT1R and CysLT2R and their corresponding ligand-optimized models. Functional assessment of the candidate hits yielded discovery of five novel antagonist chemotypes with sub-micromolar potencies and the best Ki = 220 nM at CysLT1R. One of the hits showed inverse agonism at the L129Q constitutively active mutant of CysLT2R, with potential utility against uveal melanoma.