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.

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.

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.

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 A_2A adenosine receptor (A_2AAR) molecules embedded in freely diffusing lipid nanodiscs to study their intramolecular conformational dynamics. We propose a dynamic model of A_2AAR activation that involves a slow (>2 ms) exchange between the active-like and inactive-like conformations in both apo and antagonist-bound A_2AAR, explaining the receptor's constitutive activity. For the agonist-bound A_2AAR, 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 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.

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.

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) - S1P_1-5. S1P₅ 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 S1P₅ 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 S1P₅-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

S1P₅ is a sphingosine-1-phosphate (S1P) receptor implicated in immune and neurodegenerative disorders. Here, authors report a crystal structure of the S1P₅ receptor in complex with a selective inverse agonist, revealing an allosteric subpocket and shedding light on inverse agonism in S1P receptors.

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.

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.

Structure-based mechanism of cysteinyl leukotriene receptor inhibition by antiasthmatic drugs

The G protein-coupled cysteinyl leukotriene receptor CysLT₁R mediates inflammatory processes and plays a major role in numerous disorders, including asthma, allergic rhinitis, cardiovascular disease, and cancer. Selective CysLT₁R antagonists are widely prescribed as antiasthmatic drugs; however, these drugs demonstrate low effectiveness in some patients and exhibit a variety of side effects. To gain deeper understanding into the functional mechanisms of CysLTRs, we determined the crystal structures of CysLT₁R bound to two chemically distinct antagonists, zafirlukast and pranlukast. The structures reveal unique ligand-binding modes and signaling mechanisms, including lateral ligand access to the orthosteric pocket between transmembrane helices TM4 and TM5, an atypical pattern of microswitches, and a distinct four-residue-coordinated sodium site. These results provide important insights and structural templates for rational discovery of safer and more effective drugs.

An outlook on using serial femtosecond crystallography in drug discovery

Introduction: X-ray crystallography has made important contributions to modern drug development but its application to many important drug targets has been extremely challenging. The recent emergence of X-ray free electron lasers (XFELs) and advancements in serial femtosecond crystallography (SFX) have offered new opportunities to overcome limitations of traditional crystallography to accelerate the structure-based drug discovery (SBDD) process. Areas covered: In this review, the authors describe the general principles of X-ray generation and the main properties of XFEL beams, outline details of SFX data collection and processing, and summarize the progress in the development of associated instrumentation for sample delivery and X-ray detection. An overview of the SFX applications to various important drug targets such as membrane proteins is also provided. Expert opinion: While SFX has already made clear advancements toward the understanding of the structure and dynamics of several major drug targets, its robust application in SBDD still needs further developments of new high-throughput techniques for sample production, automation of crystal delivery and data collection, as well as for processing and storage of large amounts of data. The expansion of the available XFEL beamtime is a key to the success of SFX in SBDD.

A thermostable flavin-based fluorescent protein from Chloroflexus aggregans: a framework for ultra-high resolution structural studies

A new thermostable fluorescent protein is shown to be a promising model for ultra-high resolution structural studies of LOV domains and for application as a fluorescent reporter.

Efficient non-cytotoxic fluorescent staining of halophiles

Research on halophilic microorganisms is important due to their relation to fundamental questions of survival of living organisms in a hostile environment. Here we introduce a novel method to stain halophiles with MitoTracker fluorescent dyes in their growth medium. The method is based on membrane-potential sensitive dyes, which were originally used to label mitochondria in eukaryotic cells. We demonstrate that these fluorescent dyes provide high staining efficiency and are beneficial for multi-staining purposes due to the spectral range covered (from orange to deep red). In contrast with other fluorescent dyes used so far, MitoTracker does not affect growth rate, and remains in cells after several washing steps and several generations in cell culture. The suggested dyes were tested on three archaeal (Hbt. salinarum, Haloferax sp., Halorubrum sp.) and two bacterial (Salicola sp., Halomonas sp.) strains of halophilic microorganisms. The new staining approach provides new insights into biology of Hbt. salinarum. We demonstrated the interconversion of rod-shaped cells of Hbt. salinarium to spheroplasts and submicron-sized spheres, as well as the cytoplasmic integrity of giant rod Hbt. salinarum species. By expanding the variety of tools available for halophile detection, MitoTracker dyes overcome long-standing limitations in fluorescence microscopy studies of halophiles.

Chemically Stable Lipids for Membrane Protein Crystallization

Lipidic cubic phase (LCP) has been widely recognized as a promising membrane-mimicking matrix for biophysical studies of membrane proteins and their crystallization in a lipidic environment. Application of this material to a wide variety of membrane proteins, however, is hindered due to a limited number of available host lipids, mostly monoacylglycerols (MAGs). Here, we designed, synthesized and characterized a series of chemically stable lipids resistant to hydrolysis, with properties complementary to the widely used MAGs. In order to assess their potential to serve as host lipids for crystallization, we characterized the phase properties and lattice parameters of mesophases made of two most promising lipids at a variety of different conditions by polarized light microscopy and small-angle X-ray scattering. Both lipids showed remarkable chemical stability and an extended LCP region in the phase diagram covering a wide range of temperatures down to 4 °C. One of these lipids has been used for crystallization and structure determination of a prototypical membrane protein bacteriorhodopsin at 4 °C and 20 °C.

Fast iodide-SAD phasing for high-throughput membrane protein structure determination

We describe a fast, easy, and potentially universal method for the de novo solution of the crystal structures of membrane proteins via iodide-single-wavelength anomalous diffraction (I-SAD). The potential universality of the method is based on a common feature of membrane proteins-the availability at the hydrophobic-hydrophilic interface of positively charged amino acid residues with which iodide strongly interacts. We demonstrate the solution using I-SAD of four crystal structures representing different classes of membrane proteins, including a human G protein-coupled receptor (GPCR), and we show that I-SAD can be applied using data collection strategies based on either standard or serial x-ray crystallography techniques.

Structural and functional investigation of flavin binding center of the NqrC subunit of sodium-translocating NADH:quinone oxidoreductase from Vibrio harveyi

Na⁺-translocating NADH:quinone oxidoreductase (NQR) is a redox-driven sodium pump operating in the respiratory chain of various bacteria, including pathogenic species. The enzyme has a unique set of redox active prosthetic groups, which includes two covalently bound flavin mononucleotide (FMN) residues attached to threonine residues in subunits NqrB and NqrC. The reason of FMN covalent bonding in the subunits has not been established yet. In the current work, binding of free FMN to the apo-form of NqrC from Vibrio harveyi was studied showing very low affinity of NqrC to FMN in the absence of its covalent bonding. To study structural aspects of flavin binding in NqrC, its holo-form was crystallized and its 3D structure was solved at 1.56 Å resolution. It was found that the isoalloxazine moiety of the FMN residue is buried in a hydrophobic cavity and that its pyrimidine ring is squeezed between hydrophobic amino acid residues while its benzene ring is extended from the protein surroundings. This structure of the flavin-binding pocket appears to provide flexibility of the benzene ring, which can help the FMN residue to take the bended conformation and thus to stabilize the one-electron reduced form of the prosthetic group. These properties may also lead to relatively weak noncovalent binding of the flavin. This fact along with periplasmic location of the FMN-binding domains in the vast majority of NqrC-like proteins may explain the necessity of the covalent bonding of this prosthetic group to prevent its loss to the external medium.

Low-dose X-ray radiation induces structural alterations in proteins

X-ray-radiation-induced alterations to protein structures are still a severe problem in macromolecular crystallography. One way to avoid the influence of radiation damage is to reduce the X-ray dose absorbed by the crystal during data collection. However, here it is demonstrated using the example of the membrane protein bacteriorhodopsin (bR) that even a low dose of less than 0.06 MGy may induce structural alterations in proteins. This dose is about 500 times smaller than the experimental dose limit which should ideally not be exceeded per data set (i.e. 30 MGy) and 20 times smaller than previously detected specific radiation damage at the bR active site. To date, it is the lowest dose at which radiation modification of a protein structure has been described. Complementary use was made of high-resolution X-ray crystallography and online microspectrophotometry to quantitatively study low-dose X-ray-induced changes. It is shown that structural changes of the protein correlate with the spectroscopically observed formation of the so-called bR orange species. Evidence is provided for structural modifications taking place at the protein active site that should be taken into account in crystallographic studies which aim to elucidate the molecular mechanisms of bR function.

Radioecological studies at the Kraton-3 underground nuclear explosion site in 1978-2007: a review

Within this paper, radioecological data concerning the "peaceful" underground nuclear explosion Kraton-3, conducted at a remote Arctic location (65.9 degrees N, 112.3 degrees E) within the former USSR in 1978, are reviewed. The data and estimates published in the available literature sources before September 2008 could be grouped as following: (a) characterisation of the current radioactive contamination (gamma-, beta- and alpha-emitters) of environmental compartments in terms of radionuclides composition, activity concentration, area contamination density; (b) determination of current gamma dose rates in air, including mapping using GPS; (c) evaluation of cumulative gamma doses in air (with calculations and thermoluminiscence measurements in ceramic objects); (d) description of the visually distinguishable changes in the terrestrial ecosystem; (e) description and quantitative evaluation of morphological abnormalities in the organs of adult plants as well as in seeds and seedlings of some herbs and shrubs, and in small mammals; (f) application of countermeasures. Knowledge gaps and possible further studies are indicated.

Time domain correlation analysis of heart rate variability in preterm neonates

BACKGROUND AND AIM

A fuller understanding of the neural control mechanisms of heart rate during the early stages of human development would be of great value to obstetric and neonatal management. In this paper, we investigate the correlation between heart rate variability (HRV) and other physiological parameters such as blood pressure and respiration in preterm neonates with the aim of developing a numerical model to explain and predict heart rate variability.

STUDY DESIGN AND SUBJECTS

All the required data are readily available for premature babies who are routinely monitored while being nursed in intensive care, and we have collected large data sets for a random group of such neonates. For the quantitative analysis of the data, we have developed a time domain correlation method, which has a number of advantages over the more commonly used power spectral analysis. We have been able to study the dynamics of the different frequency components of HRV by this method.

RESULTS

Highly correlated behaviour of the different HRV components, previously observed in our work on fetal HRV, is also present in the neonate, with similar characteristic time constants. Furthermore, the correlation of high-frequency (HF) oscillations of HRV with respiration and that of low-frequency (LF) oscillations of HRV with blood pressure are demonstrated on timescales of a single oscillation. In neonates receiving artificial ventilation, the correlation between HRV and respiration depends on the type of ventilation involved and assumes opposite polarities for the two main types of equipment currently in use.

CONCLUSION

We demonstrate that it is possible to analyse HRV quantitatively by calculating the relative gains and characteristic time constants for the correlated parameters and components.