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.

Structural Aspects of Photopharmacology: Insight into the Binding of Photoswitchable and Photocaged Inhibitors to the Glutamate Transporter Homologue.

Photopharmacology addresses the challenge of drug selectivity and side effects through creation of photoresponsive molecules activated with light with high spatiotemporal precision. This is achieved through incorporation of molecular photoswitches and photocages into the pharmacophore. However, the structural basis for the light-induced modulation of inhibitory potency in general is still missing, which poses a major design challenge for this emerging field of research. Here we solved crystal structures of the glutamate transporter homologue GltTk in complex with photoresponsive transport inhibitors-azobenzene derivative of TBOA (both in trans and cis configuration) and with the photocaged compound ONB-hydroxyaspartate. The essential role of glutamate transporters in the functioning of the central nervous system renders them potential therapeutic targets in the treatment of neurodegenerative diseases. The obtained structures provide a clear structural insight into the origins of photocontrol in photopharmacology and lay the foundation for application of photocontrolled ligands to study the transporter dynamics by using time-resolved X-ray crystallography.

A new twist of rubredoxin function in M. tuberculosis.

Electron transfer mediated by metalloproteins drives many biological processes. Rubredoxins are a ubiquitous [1Fe-0S] class of electron carriers that play an important role in bacterial adaptation to changing environmental conditions. In Mycobacterium tuberculosis, oxidative and acidic stresses as well as iron starvation induce rubredoxins expression. However, their functions during M. tuberculosis infection are unknown. In the present work, we show that rubredoxin B (RubB) is able to efficiently shuttle electrons from cognate reductases, FprA and FdR to support catalytic activity of cytochrome P450s, CYP124, CYP125, and CYP142, which are important for bacterial viability and pathogenicity. We solved the crystal structure of RubB and characterized the interaction between RubB and CYPs using site-directed mutagenesis. Mutations that not only neutralize single charge but also change the specific residues on the surface of RubB did not dramatically decrease activity of studied CYPs. Together with isothermal calorimetry (ITC) experiments, the obtained results suggest that interactions are transient and not highly specific. The redox potential of RubB is -264 mV vs. Ag/AgCl and the measured extinction coefficients are 9931 M-1cm-1 and 8371 M-1cm-1 at 380 nm and 490 nm, respectively. Characteristic parameters of RubB along with the discovered function might be useful for biotechnological applications. Our findings suggest that a switch from ferredoxins to rubredoxins might be crucial for M. tuberculosis to support CYPs activity during the infection.

Hydroxylation of Antitubercular Drug Candidate, SQ109, by Mycobacterial Cytochrome P450.

Spreading of the multidrug-resistant (MDR) strains of the one of the most harmful pathogen Mycobacterium tuberculosis (Mtb) generates the need for new effective drugs. SQ109 showed activity against resistant Mtb and already advanced to Phase II/III clinical trials. Fast SQ109 degradation is attributed to the human liver Cytochrome P450s (CYPs). However, no information is available about interactions of the drug with Mtb CYPs. Here, we show that Mtb CYP124, previously assigned as a methyl-branched lipid monooxygenase, binds and hydroxylates SQ109 in vitro. A 1.25 Å-resolution crystal structure of the CYP124-SQ109 complex unambiguously shows two conformations of the drug, both positioned for hydroxylation of the ω-methyl group in the trans position. The hydroxylated SQ109 presumably forms stabilizing H-bonds with its target, Mycobacterial membrane protein Large 3 (MmpL3). We anticipate that Mtb CYPs could function as analogs of drug-metabolizing human CYPs affecting pharmacokinetics and pharmacodynamics of antitubercular (anti-TB) drugs.

Custom Design of a Humidifier Chamber for InMeso Crystallization.

Membrane proteins are indispensable for every living organism, yet their structural organization remains underexplored. Despite the recent advancements in single-particle cryogenic electron microscopy and cryogenic electron tomography, which have significantly increased the structural coverage of membrane proteins across various kingdoms, certain scientific methods, such as time-resolved crystallography, still mostly rely on crystallization techniques, such as lipidic cubic phase (LCP) or in meso crystallization. In this study, we present an open-access blueprint for a humidity control chamber designed for LCP/in meso crystallization experiments using a Gryphon crystallization robot. Using this chamber, we have obtained crystals of a transmembrane aspartate transporter GltTk from Thermococcus kodakarensis in a lipidic environment using in meso crystallization. The data collected from these crystals allowed us to perform an analysis of lipids bound to the transporter. With this publication of our open-access design of a humidity chamber, we aim to improve the accessibility of in meso protein crystallization for the scientific community.

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

Mirror proteorhodopsins.

Proteorhodopsins (PRs), bacterial light-driven outward proton pumps comprise the first discovered and largest family of rhodopsins, they play a significant role in life on the Earth. A big remaining mystery was that up-to-date there was no described bacterial rhodopsins pumping protons at acidic pH despite the fact that bacteria live in different pH environment. Here we describe conceptually new bacterial rhodopsins which are operating as outward proton pumps at acidic pH. A comprehensive function-structure study of a representative of a new clade of proton pumping rhodopsins which we name "mirror proteorhodopsins", from Sphingomonas paucimobilis (SpaR) shows cavity/gate architecture of the proton translocation pathway rather resembling channelrhodopsins than the known rhodopsin proton pumps. Another unique property of mirror proteorhodopsins is that proton pumping is inhibited by a millimolar concentration of zinc. We also show that mirror proteorhodopsins are extensively represented in opportunistic multidrug resistant human pathogens, plant growth-promoting and zinc solubilizing bacteria. They may be of optogenetic interest.

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

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.

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.

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.

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.

Molecular mechanism of light-driven sodium pumping.

The light-driven sodium-pumping rhodopsin KR2 from Krokinobacter eikastus is the only non-proton cation active transporter with demonstrated potential for optogenetics. However, the existing structural data on KR2 correspond exclusively to its ground state, and show no sodium inside the protein, which hampers the understanding of sodium-pumping mechanism. Here we present crystal structure of the O-intermediate of the physiologically relevant pentameric form of KR2 at the resolution of 2.1 Å, revealing a sodium ion near the retinal Schiff base, coordinated by N112 and D116 of the characteristic NDQ triad. We also obtained crystal structures of D116N and H30A variants, conducted metadynamics simulations and measured pumping activities of putative pathway mutants to demonstrate that sodium release likely proceeds alongside Q78 towards the structural sodium ion bound between KR2 protomers. Our findings highlight the importance of pentameric assembly for sodium pump function, and may be used for rational engineering of enhanced optogenetic tools.

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.

Structural basis of ligand selectivity and disease mutations in cysteinyl leukotriene receptors.

Cysteinyl leukotriene G protein-coupled receptors CysLT1 and CysLT2 regulate pro-inflammatory responses associated with allergic disorders. While selective inhibition of CysLT1R has been used for treating asthma and associated diseases for over two decades, CysLT2R has recently started to emerge as a potential drug target against atopic asthma, brain injury and central nervous system disorders, as well as several types of cancer. Here, we describe four crystal structures of CysLT2R in complex with three dual CysLT1R/CysLT2R antagonists. The reported structures together with the results of comprehensive mutagenesis and computer modeling studies shed light on molecular determinants of CysLTR ligand selectivity and specific effects of disease-related single nucleotide variants.

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

The G protein-coupled cysteinyl leukotriene receptor CysLT1R mediates inflammatory processes and plays a major role in numerous disorders, including asthma, allergic rhinitis, cardiovascular disease, and cancer. Selective CysLT1R 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 CysLT1R 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.

Structural insights into ion conduction by channelrhodopsin 2.

The light-gated ion channel channelrhodopsin 2 (ChR2) from Chlamydomonas reinhardtii is a major optogenetic tool. Photon absorption starts a well-characterized photocycle, but the structural basis for the regulation of channel opening remains unclear. We present high-resolution structures of ChR2 and the C128T mutant, which has a markedly increased open-state lifetime. The structure reveals two cavities on the intracellular side and two cavities on the extracellular side. They are connected by extended hydrogen-bonding networks involving water molecules and side-chain residues. Central is the retinal Schiff base that controls and synchronizes three gates that separate the cavities. Separate from this network is the DC gate that comprises a water-mediated bond between C128 and D156 and interacts directly with the retinal Schiff base. Comparison with the C128T structure reveals a direct connection of the DC gate to the central gate and suggests how the gating mechanism is affected by subtle tuning of the Schiff base's interactions.