Protein structure prediction provides comparable performance to crystallographic structures in docking-based virtual screening

Fragment-Based Drug Discovery against Mycobacteria: The Success and Challenges

The emergence of drug-resistant mycobacteria, including Mycobacterium tuberculosis (Mtb) and non-tuberculous mycobacteria (NTM), poses an increasing global threat that urgently demands the development of new potent anti-mycobacterial drugs. One of the approaches toward the identification of new drugs is fragment-based drug discovery (FBDD), which is the most ingenious among other drug discovery models, such as structure-based drug design (SBDD) and high-throughput screening. Specialized techniques, such as X-ray crystallography, nuclear magnetic resonance spectroscopy, and many others, are part of the drug discovery approach to combat the Mtb and NTM global menaces. Moreover, the primary drawbacks of traditional methods, such as the limited measurement of biomolecular toxicity and uncertain bioavailability evaluation, are successfully overcome by the FBDD approach. The current review focuses on the recognition of fragment-based drug discovery as a popular approach using virtual, computational, and biophysical methods to identify potent fragment molecules. FBDD focuses on designing optimal inhibitors against potential therapeutic targets of NTM and Mtb (PurC, ArgB, MmpL3, and TrmD). Additionally, we have elaborated on the challenges associated with the FBDD approach in the identification and development of novel compounds. Insights into the applications and overcoming the challenges of FBDD approaches will aid in the identification of potential therapeutic compounds to treat drug-sensitive and drug-resistant NTMs and Mtb infections.

Computational modeling and druggability assessment of Aggregatibacter actinomycetemcomitans leukotoxin

The leukotoxin (LtxA) of Aggregatibacter actinomycetemcomitans (A. actinomycetemcomitans) is a protein exotoxin belonging to the repeat-in-toxin family (RTX). Numerous studies have demonstrated that LtxA may play a critical role in the pathogenicity of A. actinomycetemcomitans since hyper-leukotoxic strains have been associated with severe disease. Accordingly, considerable effort has been made to elucidate the mechanisms by which LtxA interacts with host cells and induce their death. However, these attempts have been hampered by the unavailability of a tertiary structure of the toxin, which limits the understanding of its molecular properties and mechanisms. In this paper, we used homology and template free modeling algorithms to build the complete tertiary model of LtxA at atomic level in its calcium-bound Holo-state. The resulting model was refined by energy minimization, validated by Molprobity and ProSA tools, and subsequently subjected to a cumulative 600ns of all-atom classical molecular dynamics simulation to evaluate its structural aspects. The druggability of the proposed model was assessed using Fpocket and FTMap tools, resulting in the identification of four putative cavities and fifteen binding hotspots that could be targeted by rational drug design tools to find new ligands to inhibit LtxA activity.

Structure-Based Design of First-Generation Small Molecule Inhibitors Targeting the Catalytic Pockets of AID, APOBEC3A, and APOBEC3B

Activation-induced cytidine deaminase (AID) initiates antibody diversification by mutating immunoglobulin loci in B lymphocytes. AID and related APOBEC3 (A3) enzymes also induce genome-wide mutations and lesions implicated in tumorigenesis and tumor progression. The most prevalent mutation signatures across diverse tumor genomes are attributable to the mistargeted mutagenic activities of AID/A3s. Thus, inhibiting AID/A3s has been suggested to be of therapeutic benefit. We previously used a computational-biochemical approach to gain insight into the structure of AID’s catalytic pocket, which resulted in the discovery of a novel type of regulatory catalytic pocket closure that regulates AID/A3s that we termed the “Schrodinger’s CATalytic pocket”. Our findings were subsequently confirmed by direct structural studies. Here, we describe our search for small molecules that target the catalytic pocket of AID. We identified small molecules that inhibit purified AID, AID in cell extracts, and endogenous AID of lymphoma cells. Analogue expansion yielded derivatives with improved potencies. These were found to also inhibit A3A and A3B, the two most tumorigenic siblings of AID. Two compounds exhibit low micromolar IC₅₀ inhibition of AID and A3A, exhibiting the strongest potency for A3A. Docking suggests key interactions between their warheads and residues lining the catalytic pockets of AID, A3A, and A3B and between the tails and DNA-interacting residues on the surface proximal to the catalytic pocket opening. Accordingly, mutants of these residues decreased inhibition potency. The chemistry and abundance of key stabilizing interactions between the small molecules and residues within and immediately outside the catalytic pockets are promising for therapeutic development.

Comparison of structure- and ligand-based scoring functions for deep generative models: a GPCR case study

Deep generative models have shown the ability to devise both valid and novel chemistry, which could significantly accelerate the identification of bioactive compounds. Many current models, however, use molecular descriptors or ligand-based predictive methods to guide molecule generation towards a desirable property space. This restricts their application to relatively data-rich targets, neglecting those where little data is available to sufficiently train a predictor. Moreover, ligand-based approaches often bias molecule generation towards previously established chemical space, thereby limiting their ability to identify truly novel chemotypes. In this work, we assess the ability of using molecular docking via Glide-a structure-based approach-as a scoring function to guide the deep generative model REINVENT and compare model performance and behaviour to a ligand-based scoring function. Additionally, we modify the previously published MOSES benchmarking dataset to remove any induced bias towards non-protonatable groups. We also propose a new metric to measure dataset diversity, which is less confounded by the distribution of heavy atom count than the commonly used internal diversity metric. With respect to the main findings, we found that when optimizing the docking score against DRD2, the model improves predicted ligand affinity beyond that of known DRD2 active molecules. In addition, generated molecules occupy complementary chemical and physicochemical space compared to the ligand-based approach, and novel physicochemical space compared to known DRD2 active molecules. Furthermore, the structure-based approach learns to generate molecules that satisfy crucial residue interactions, which is information only available when taking protein structure into account. Overall, this work demonstrates the advantage of using molecular docking to guide de novo molecule generation over ligand-based predictors with respect to predicted affinity, novelty, and the ability to identify key interactions between ligand and protein target. Practically, this approach has applications in early hit generation campaigns to enrich a virtual library towards a particular target, and also in novelty-focused projects, where de novo molecule generation either has no prior ligand knowledge available or should not be biased by it.

Prediction of Functional Consequences of Missense Mutations in ANO4 Gene

The anoctamin (TMEM16) family of transmembrane protein consists of ten members in vertebrates, which act as Ca²⁺-dependent ion channels and/or Ca²⁺-dependent scramblases. ANO4 which is primarily expressed in the CNS and certain endocrine glands, has been associated with various neuronal disorders. Therefore, we focused our study on prioritizing missense mutations that are assumed to alter the structure and stability of ANO4 protein. We employed a wide array of evolution and structure based in silico prediction methods to identify potentially deleterious missense mutations in the ANO4 gene. Identified pathogenic mutations were then mapped to the modeled human ANO4 structure and the effects of missense mutations were studied on the atomic level using molecular dynamics simulations. Our data show that the G80A and A500T mutations significantly alter the stability of the mutant proteins, thus providing new perspective on the role of missense mutations in ANO4 gene. Results obtained in this study may help to identify disease associated mutations which affect ANO4 protein structure and function and might facilitate future functional characterization of ANO4.

Protein X-ray Crystallography and Drug Discovery

With the advent of structural biology in the drug discovery process, medicinal chemists gained the opportunity to use detailed structural information in order to progress screening hits into leads or drug candidates. X-ray crystallography has proven to be an invaluable tool in this respect, as it is able to provide exquisitely comprehensive structural information about the interaction of a ligand with a pharmacological target. As fragment-based drug discovery emerged in the recent years, X-ray crystallography has also become a powerful screening technology, able to provide structural information on complexes involving low-molecular weight compounds, despite weak binding affinities. Given the low numbers of compounds needed in a fragment library, compared to the hundreds of thousand usually present in drug-like compound libraries, it now becomes feasible to screen a whole fragment library using X-ray crystallography, providing a wealth of structural details that will fuel the fragment to drug process. Here, we review theoretical and practical aspects as well as the pros and cons of using X-ray crystallography in the drug discovery process.

Rational design of agonists for bitter taste receptor TAS2R14: from modeling to bench and back

Human bitter taste receptors (TAS2Rs) are a subfamily of 25 G protein-coupled receptors that mediate bitter taste perception. TAS2R14 is the most broadly tuned bitter taste receptor, recognizing a range of chemically diverse agonists with micromolar-range potency. The receptor is expressed in several extra-oral tissues and is suggested to have physiological roles related to innate immune responses, male fertility, and cancer. Higher potency ligands are needed to investigate TAS2R14 function and to modulate it for future clinical applications. Here, a structure-based modeling approach is described for the design of TAS2R14 agonists beginning from flufenamic acid, an approved non-steroidal anti-inflammatory analgesic that activates TAS2R14 at sub-micromolar concentrations. Structure-based molecular modeling was integrated with experimental data to design new TAS2R14 agonists. Subsequent chemical synthesis and in vitro profiling resulted in new TAS2R14 agonists with improved potency compared to the lead. The integrated approach provides a validated and refined structural model of ligand-TAS2R14 interactions and a general framework for structure-based discovery in the absence of closely related experimental structures.

A Chemosensory GPCR as a Potential Target to Control the Root-Knot Nematode Meloidogyne incognita Parasitism in Plants

Root-knot nematodes (RKN), from the Meloidogyne genus, have a worldwide distribution and cause severe economic damage to many life-sustaining crops. Because of their lack of specificity and danger to the environment, most chemical nematicides have been banned from use. Thus, there is a great need for new and safe compounds to control RKN. Such research involves identifying beforehand the nematode proteins essential to the invasion. Since G protein-coupled receptors GPCRs are the target of a large number of drugs, we have focused our research on the identification of putative nematode GPCRs such as those capable of controlling the movement of the parasite towards (or within) its host. A datamining procedure applied to the genome of Meloidogyne incognita allowed us to identify a GPCR, belonging to the neuropeptide GPCR family that can serve as a target to carry out a virtual screening campaign. We reconstructed a 3D model of this receptor by homology modeling and validated it through extensive molecular dynamics simulations. This model was used for large scale molecular dockings which produced a filtered limited set of putative antagonists for this GPCR. Preliminary experiments using these selected molecules allowed the identification of an active compound, namely C260-2124, from the ChemDiv provider, which can serve as a starting point for further investigations.

In silico prediction of prolactin molecules as a tool for equine genomics reproduction

The prolactin hormone is involved in several biological functions, although its main role resides on reproduction. As it interferes on fertility changes, studies focused on human health have established a linkage of this hormone to fertility losses. Regarding animal research, there is still a lack of information about the structure of prolactin. In case of horse breeding, prolactin has a particular influence; once there is an individualization of these animals and equines are known for presenting several reproductive disorders. As there is no molecular structure available for the prolactin hormone and receptor, we performed several bioinformatics analyses through prediction and refinement softwares, as well as manual modifications. Aiming to elucidate the first computational structure of both molecules and analyse structural and functional aspects related to these proteins, here we provide the first known equine model for prolactin and prolactin receptor, which obtained high global quality scores in diverse software's for quality assessment. QMEAN overall score obtained for ePrl was (- 4.09) and QMEANbrane for ePrlr was (- 8.45), which proves the structures' reliability. This study will implement another tool in equine genomics in order to give light to interactions of these molecules, structural and functional alterations and therefore help diagnosing fertility problems, contributing in the selection of a high genetic herd.

Dynamics fingerprints of active conformers of epidermal growth factor receptor kinase

Epidermal growth factor receptor (EGFR) is a prototypical cell-surface receptor that plays a key role in the regulation of cellular signaling, proliferation and differentiation. Mutations of its kinase domain have been associated with the development of a variety of cancers and, therefore, it has been the target of drug design. Single amino acid substitutions (SASs) in this domain have been proven to alter the equilibrium of pre-existing conformer populations. Despite the advances in structural descriptions of its so-called active and inactive conformations, the associated dynamics aspects that characterize them have not been thoroughly studied yet. As the dynamic behaviors and molecular motions of proteins are important for a complete understanding of their structure-function relationships we present a novel procedure, using (or based on) normal mode analysis, to identify the collective dynamics shared among different conformers in EGFR kinase. The method allows the comparison of patterns of low-frequency vibrational modes defining representative directions of motions. Our procedure is able to emphasize the main similarities and differences between the collective dynamics of different conformers. In the case of EGFR kinase, two representative directions of motions have been found as dynamics fingerprints of the active conformers. Protein motion along both directions reveals to have a significant impact on the cavity volume of the main pocket of the active site. Otherwise, the inactive conformers exhibit a more heterogeneous distribution of collective motions. © 2018 Wiley Periodicals, Inc.

Homology Model and Docking-Based Virtual Screening for Ligands of Human Dyskerin as New Inhibitors of Telomerase for Cancer Treatment

Immortality is one of the main features of cancer cells. Tumor cells have an unlimited replicative potential, principally due to the holoenzyme telomerase. Telomerase is composed mainly by dyskerin (DKC1), a catalytic retrotranscriptase (hTERT) and an RNA template (hTR). The aim of this work is to develop new inhibitors of telomerase, selecting the interaction between hTR⁻DKC1 as a target. We designed two models of the human protein DKC1: homology and ab initio. These models were evaluated by different procedures, revealing that the homology model parameters were the most accurate. We selected two hydrophobic pockets contained in the PUA (pseudouridine synthase and archaeosine transglycosylase) domain, using structural and stability analysis. We carried out a docking-based virtual screen on these pockets, using the reported mutation K314 as the center of the docking. The hDKC1 model was tested against a library of 450,000 drug-like molecules. We selected the first 10 molecules that showed the highest affinity values to test their inhibitory activity on the cell line MDA MB 231 (Monroe Dunaway Anderson Metastasis Breast cancer 231), obtaining three compounds that showed inhibitory effect. These results allowed us to validate our design and set the basis to continue with the study of telomerase inhibitors for cancer treatment.

Homology Modeling of 5-alpha-Reductase 2 Using Available Experimental Data

5-Alpha-reductase 2 is an interesting pharmaceutical target for the treatment of several diseases, including prostate cancer, benign prostatic hyperplasia, male pattern baldness, acne, and hirsutism. One of the main approaches in computer aided drug design is structure-based drug discovery. However, the experimental 3D structure of 5-alpha-reductase 2 is not available at present. Therefore, a homology modeling method and molecular dynamics simulation were used to develop a reliable model of 5-alpha-reductase 2 for inhibitor pose prediction and virtual screening. Despite the low sequence identity between the target and template sequences, a useful 3D model of 5-alpha-reductase 2 was generated by the inclusion of experimental data.

Computational protein structure refinement: Almost there, yet still so far to go

Protein structures are essential in modern biology yet experimental methods are far from being able to catch up with the rapid increase in available genomic data. Computational protein structure prediction methods aim to fill the gap while the role of protein structure refinement is to take approximate initial template-based models and bring them closer to the true native structure. Current methods for computational structure refinement rely on molecular dynamics simulations, related sampling methods, or iterative structure optimization protocols. The best methods are able to achieve moderate degrees of refinement but consistent refinement that can reach near-experimental accuracy remains elusive. Key issues revolve around the accuracy of the energy function, the inability to reliably rank multiple models, and the use of restraints that keep sampling close to the native state but also limit the degree of possible refinement. A different aspect is the question of what exactly the target of high-resolution refinement should be as experimental structures are affected by experimental conditions and different biological questions require varying levels of accuracy. While improvement of the global protein structure is a difficult problem, high-resolution refinement methods that improves local structural quality such as favorable stereochemistry and the avoidance of atomic clashes are much more successful.

GPCRs from fusarium graminearum detection, modeling and virtual screening - the search for new routes to control head blight disease

BACKGOUND

Fusarium graminearum (FG) is one of the major cereal infecting pathogens causing high economic losses worldwide and resulting in adverse effects on human and animal health. Therefore, the development of new fungicides against FG is an important issue to reduce cereal infection and economic impact. In the strategy for developing new fungicides, a critical step is the identification of new targets against which innovative chemicals weapons can be designed. As several G-protein coupled receptors (GPCRs) are implicated in signaling pathways critical for the fungi development and survival, such proteins could be valuable efficient targets to reduce Fusarium growth and therefore to prevent food contamination.

RESULTS

In this study, GPCRs were predicted in the FG proteome using a manually curated pipeline dedicated to the identification of GPCRs. Based on several successive filters, the most appropriate GPCR candidate target for developing new fungicides was selected. Searching for new compounds blocking this particular target requires the knowledge of its 3D-structure. As no experimental X-Ray structure of the selected protein was available, a 3D model was built by homology modeling. The model quality and stability was checked by 100 ns of molecular dynamics simulations. Two stable conformations representative of the conformational families of the protein were extracted from the 100 ns simulation and were used for an ensemble docking campaign. The model quality and stability was checked by 100 ns of molecular dynamics simulations previously to the virtual screening step. The virtual screening step comprised the exploration of a chemical library with 11,000 compounds that were docked to the GPCR model. Among these compounds, we selected the ten top-ranked nontoxic molecules proposed to be experimentally tested to validate the in silico simulation.

CONCLUSIONS

This study provides an integrated process merging genomics, structural bioinformatics and drug design for proposing innovative solutions to a world wide threat to grain producers and consumers.

Illustrating and homology modeling the proteins of the Zika virus

The Zika virus (ZIKV) is a flavivirus of the family Flaviviridae, which is similar to dengue virus, yellow fever and West Nile virus. Recent outbreaks in South America, Latin America, the Caribbean and in particular Brazil have led to concern for the spread of the disease and potential to cause Guillain-Barré syndrome and microcephaly. Although ZIKV has been known of for over 60 years there is very little in the way of knowledge of the virus with few publications and no crystal structures. No antivirals have been tested against it either in vitro or in vivo. ZIKV therefore epitomizes a neglected disease. Several suggested steps have been proposed which could be taken to initiate ZIKV antiviral drug discovery using both high throughput screens as well as structure-based design based on homology models for the key proteins. We now describe preliminary homology models created for NS5, FtsJ, NS4B, NS4A, HELICc, DEXDc, peptidase S7, NS2B, NS2A, NS1, E stem, glycoprotein M, propeptide, capsid and glycoprotein E using SWISS-MODEL. Eleven out of 15 models pass our model quality criteria for their further use. While a ZIKV glycoprotein E homology model was initially described in the immature conformation as a trimer, we now describe the mature dimer conformer which allowed the construction of an illustration of the complete virion. By comparing illustrations of ZIKV based on this new homology model and the dengue virus crystal structure we propose potential differences that could be exploited for antiviral and vaccine design. The prediction of sites for glycosylation on this protein may also be useful in this regard. While we await a cryo-EM structure of ZIKV and eventual crystal structures of the individual proteins, these homology models provide the community with a starting point for structure-based design of drugs and vaccines as well as a for computational virtual screening.

The Hevea brasiliensis XIP aquaporin subfamily: genomic, structural and functional characterizations with relevance to intensive latex harvesting

X-Intrinsic Proteins (XIP) were recently identified in a narrow range of plants as a full clade within the aquaporins. These channels reportedly facilitate the transport of a wide range of hydrophobic solutes. The functional roles of XIP in planta remain poorly identified. In this study, we found three XIP genes (HbXIP1;1, HbXIP2;1 and HbXIP3;1) in the Hevea brasiliensis genome. Comprehensive bioinformatics, biochemical and structural analyses were used to acquire a better understanding of this AQP subfamily. Phylogenetic analysis revealed that HbXIPs clustered into two major groups, each distributed in a specific lineage of the order Malpighiales. Tissue-specific expression profiles showed that only HbXIP2;1 was expressed in all the vegetative tissues tested (leaves, stem, bark, xylem and latex), suggesting that HbXIP2;1 could take part in a wide range of cellular processes. This is particularly relevant to the rubber-producing laticiferous system, where this isoform was found to be up-regulated during tapping and ethylene treatments. Furthermore, the XIP transcriptional pattern is significantly correlated to latex production level. Structural comparison with SoPIP2;1 from Spinacia oleracea species provides new insights into the possible role of structural checkpoints by which HbXIP2;1 ensures glycerol transfer across the membrane. From these results, we discuss the physiological involvement of glycerol and HbXIP2;1 in water homeostasis and carbon stream of challenged laticifers. The characterization of HbXIP2;1 during rubber tree tapping lends new insights into molecular and physiological response processes of laticifer metabolism in the context of latex exploitation.

Illustrating and homology modeling the proteins of the Zika virus

The Zika virus (ZIKV) is a flavivirus of the family Flaviviridae, which is similar to dengue virus, yellow fever and West Nile virus. Recent outbreaks in South America, Latin America, the Caribbean and in particular Brazil have led to concern for the spread of the disease and potential to cause Guillain-Barré syndrome and microcephaly. Although ZIKV has been known of for over 60 years there is very little in the way of knowledge of the virus with few publications and no crystal structures. No antivirals have been tested against it either in vitro or in vivo. ZIKV therefore epitomizes a neglected disease. Several suggested steps have been proposed which could be taken to initiate ZIKV antiviral drug discovery using both high throughput screens as well as structure-based design based on homology models for the key proteins. We now describe preliminary homology models created for NS5, FtsJ, NS4B, NS4A, HELICc, DEXDc, peptidase S7, NS2B, NS2A, NS1, E stem, glycoprotein M, propeptide, capsid and glycoprotein E using SWISS-MODEL. Eleven out of 15 models pass our model quality criteria for their further use. While a ZIKV glycoprotein E homology model was initially described in the immature conformation as a trimer, we now describe the mature dimer conformer which allowed the construction of an illustration of the complete virion. By comparing illustrations of ZIKV based on this new homology model and the dengue virus crystal structure we propose potential differences that could be exploited for antiviral and vaccine design. The prediction of sites for glycosylation on this protein may also be useful in this regard. While we await a cryo-EM structure of ZIKV and eventual crystal structures of the individual proteins, these homology models provide the community with a starting point for structure-based design of drugs and vaccines as well as a for computational virtual screening.

Can template-based protein models guide the design of sequence fitness for enhanced thermal stability of single domain antibodies?

We investigate the practical use of comparative (template-based) protein models in replica-exchange simulations of single-domain antibody (sdAb) chains to evaluate if the models can correctly predict in rank order the thermal susceptibility to unfold relative to experimental melting temperatures. The baseline model system is the recently determined crystallographic structure of a llama sdAb (denoted as A3), which exhibits an unusually high thermal stability. An evaluation of the simulation results for the A3 comparative model and crystal structure shows that, despite the overall low Cα root-mean-square deviation between the two structures, the model contains misfolded regions that yields a thermal profile of unraveling at a lower temperature. Yet comparison of the simulations of four different comparative models for sdAb A3, C8, A3C8 and E9, where A3C8 is a design of swapping the sequence of the complementarity determining regions of C8 onto the A3 framework, discriminated among the sequences to detect the highest and lowest experimental melting transition temperatures. Further structural analysis of A3 for selected alanine substitutions by a combined computational and experimental study found unexpectedly that the comparative model performed admirably in recognizing substitution 'hot spots' when using a support-vector machine algorithm.