Synthesis and Evaluation of 5-(Heteroarylmethylene)hydantoins as Glycogen Synthase Kinase-3β Inhibitors

Glycogen synthase kinase-3 (GSK-3) is a serine/threonine kinase which plays a center role in the phosphorylation of a wide variety of proteins, generally leading to their inactivation. As such, GSK-3 is viewed as a therapeutic target. An ever-increasing number of small organic molecule inhibitors of GSK-3 have been reported. Phenylmethylene hydantoins are known to exhibit a wide range of inhibitory activities including for GSK-3β. A family of fourteen 2-heterocycle substituted methylene hydantoins (14, 17-29) were prepared and evaluated for the inhibition of GSK-3β at 25 μM. The IC50 values of five of these compounds was determined; the two best inhibitors are 5-[(4'-chloro-2-pyridinyl)methylene]hydantoin (IC50 = 2.14 ± 0.18 μM) and 5-[(6'-bromo-2-pyridinyl)methylene]hydantoin (IC50 = 3.39 ± 0.16 μM). The computational docking of the compounds with GSK-3β (pdb 1q41) revealed poses with hydrogen bonding to the backbone at Val135. The 5-[(heteroaryl)methylene]hydantoins did not strongly inhibit other metalloenzymes, demonstrating poor inhibitory activity against matrix metalloproteinase-12 at 25 μM and against human carbonic anhydrase at 200 μM, and were not inhibitors for Staphylococcus aureus pyruvate carboxylase at concentrations >1000 μM.

Assessment of fragment docking and scoring with the endothiapepsin model system

Fragment-based screening has become indispensable in drug discovery. Yet, the weak binding affinities of these small molecules still represent a challenge for the reliable detection of fragment hits. The extent of this issue was illustrated in the literature for the aspartic protease endothiapepsin: When seven biochemical and biophysical in vitro screening methods were applied to screen a library of 361 fragments, very poor overlap was observed between the hit fragments identified by the individual approaches, resulting in high levels of false positive and/or false negative results depending on the mutually compared methods. Here, the reported in vitro findings are juxtaposed with the results from in silico docking and scoring approaches. The docking programs GOLD and Glide were considered with the scoring functions ASP, ChemScore, ChemPLP, GoldScore, DSXCSD, and GlideScore. First, the ranking power and scoring power were assessed for the named scoring functions. Second, the capability of reproducing the crystallized fragment binding modes was tested in a structure-based redocking approach. The redocking success notably depended on the ligand efficiency of the considered fragments. Third, a blinded virtual screening approach was employed to evaluate whether in silico screening can compete with in vitro methods in the enrichment of fragment databases.

Integration of Computational Docking into Anti-Cancer Drug Response Prediction Models

Cancer is a heterogeneous disease in that tumors of the same histology type can respond differently to a treatment. Anti-cancer drug response prediction is of paramount importance for both drug development and patient treatment design. Although various computational methods and data have been used to develop drug response prediction models, it remains a challenging problem due to the complexities of cancer mechanisms and cancer-drug interactions. To better characterize the interaction between cancer and drugs, we investigate the feasibility of integrating computationally derived features of molecular mechanisms of action into prediction models. Specifically, we add docking scores of drug molecules and target proteins in combination with cancer gene expressions and molecular drug descriptors for building response models. The results demonstrate a marginal improvement in drug response prediction performance when adding docking scores as additional features, through tests on large drug screening data. We discuss the limitations of the current approach and provide the research community with a baseline dataset of the large-scale computational docking for anti-cancer drugs.

Structure-Based Discovery of Small-Molecule Inhibitors of the Autocatalytic Proliferation of α-Synuclein Aggregates

The presence of amyloid fibrils of α-synuclein is closely associated with Parkinson's disease and related synucleinopathies. It is still very challenging, however, to systematically discover small molecules that prevent the formation of these aberrant aggregates. Here, we describe a structure-based approach to identify small molecules that specifically inhibit the surface-catalyzed secondary nucleation step in the aggregation of α-synuclein by binding to the surface of the amyloid fibrils. The resulting small molecules are screened using a range of kinetic and thermodynamic assays for their ability to bind α-synuclein fibrils and prevent the further generation of α-synuclein oligomers. This study demonstrates that the combination of structure-based and kinetic-based drug discovery methods can lead to the identification of small molecules that selectively inhibit the autocatalytic proliferation of α-synuclein aggregates.

Evaluation of AlphaFold2 structures as docking targets

AlphaFold2 is a promising new tool for researchers to predict protein structures and generate high-quality models, with low backbone and global root-mean-square deviation (RMSD) when compared with experimental structures. However, it is unclear if the structures predicted by AlphaFold2 will be valuable targets of docking. To address this question, we redocked ligands in the PDBbind datasets against the experimental co-crystallized receptor structures and against the AlphaFold2 structures using AutoDock-GPU. We find that the quality measure provided during structure prediction is not a good predictor of docking performance, despite accurately reflecting the quality of the alpha carbon alignment with experimental structures. Removing low-confidence regions of the predicted structure and making side chains flexible improves the docking outcomes. Overall, despite high-quality prediction of backbone conformation, fine structural details limit the naive application of AlphaFold2 models as docking targets.

X-ray structure and mechanism of ZgHAD, a l-2-haloacid dehalogenase from the marine Flavobacterium Zobellia galactanivorans

Haloacid dehalogenases are potentially involved in bioremediation of contaminated environments and few have been biochemically characterized from marine organisms. The l-2-haloacid dehalogenase (l-2-HAD) from the marine Bacteroidetes Zobellia galactanivorans Dsij^T (ZgHAD) has been shown to catalyze the dehalogenation of C2 and C3 short-chain l-2-haloalkanoic acids. To better understand its catalytic properties, its enzymatic stability, active site, and 3D structure were analyzed. ZgHAD demonstrates high stability to solvents and a conserved catalytic activity when heated up to 60°C, its melting temperature being at 65°C. The X-ray structure of the recombinant enzyme was solved by molecular replacement. The enzyme folds as a homodimer and its active site is very similar to DehRhb, the other known l-2-HAD from a marine Rhodobacteraceae. Marked differences are present in the putative substrate entrance sites of the two enzymes. The H179 amino acid potentially involved in the activation of a catalytic water molecule was confirmed as catalytic amino acid through the production of two inactive site-directed mutants. The crystal packing of 13 dimers in the asymmetric unit of an active-site mutant, ZgHAD-H179N, reveals domain movements of the monomeric subunits relative to each other. The involvement of a catalytic His/Glu dyad and substrate binding amino acids was further confirmed by computational docking. All together our results give new insights into the catalytic mechanism of the group of marine l-2-HAD.

Exploring the Inhibition of Quercetin on Acetylcholinesterase by Multispectroscopic and In Silico Approaches and Evaluation of Its Neuroprotective Effects on PC12 Cells

This study investigated the inhibitory mechanism of quercetin in acetylcholinesterase (AChE) and its neuroprotective effects on β-amyloid_25-35-induced oxidative stress injury in PC12 cells. Quercetin inhibited AChE in a reversible mixed manner with an IC₅₀ of 4.59 ± 0.27 µM. The binding constant of quercetin with AChE at 25 °C was (5.52 ± 0.05) × 10⁴ L mol^-1. Hydrogen bonding and van der Waals forces were the main interactions in forming the stable quercetin-AChE complex. Computational docking revealed that quercetin was dominant at the peripheral aromatic site in AChE and induced enzymatic allosterism; meanwhile, it extended deep into the active center of AChE and destabilized the hydrogen bond network, which caused the constriction of the gorge entrance and prevented the substrate from entering the enzyme, thus resulting in the inhibition of AChE. Molecular dynamics (MD) simulation emphasized the stability of the quercetin-AChE complex and corroborated the previous findings. Interestingly, a combination of galantamine hydrobromide and quercetin exhibited the synergistic inhibition effect by binding to different active sites of AChE. In a β-amyloid_25-35-induced oxidative stress injury model in PC12 cells, quercetin exerted neuroprotective effects by increasing the glutathione level and reducing the malondialdehyde content and reactive oxygen species levels. These findings may provide novel insights into the development and application of quercetin in the dietary treatment of Alzheimer's disease.

Computational Docking Reveals Co-Evolution of C4 Carbon Delivery Enzymes in Diverse Plants

Proteins are modular functionalities regulating multiple cellular activities in prokaryotes and eukaryotes. As a consequence of higher plants adapting to arid and thermal conditions, C4 photosynthesis is the carbon fixation process involving multi-enzymes working in a coordinated fashion. However, how these enzymes interact with each other and whether they co-evolve in parallel to maintain interactions in different plants remain elusive to date. Here, we report our findings on the global protein co-evolution relationship and local dynamics of co-varying site shifts in key C4 photosynthetic enzymes. We found that in most of the selected key C4 photosynthetic enzymes, global pairwise co-evolution events exist to form functional couplings. Besides, protein-protein interactions between these enzymes may suggest their unknown functionalities in the carbon delivery process. For PEPC and PPCK regulation pairs, pocket formation at the interactive interface are not necessary for their function. This feature is distinct from another well-known regulation pair in C4 photosynthesis, namely, PPDK and PPDK-RP, where the pockets are necessary. Our findings facilitate the discovery of novel protein regulation types and contribute to expanding our knowledge about C4 photosynthesis.

Antamanide Analogs as Potential Inhibitors of Tyrosinase

The tyrosinase enzyme, which catalyzes the hydroxylation of monophenols and the oxidation of o-diphenols, is typically involved in the synthesis of the dark product melanin starting from the amino acid tyrosine. Contributing to the browning of plant and fruit tissues and to the hyperpigmentation of the skin, leading to melasma or age spots, the research of possible tyrosinase inhibitors has attracted much interest in agri-food, cosmetic, and medicinal industries. In this study, we analyzed the capability of antamanide, a mushroom bioactive cyclic decapeptide, and some of its glycine derivatives, compared to that of pseudostellarin A, a known tyrosinase inhibitor, to hinder tyrosinase activity by using a spectrophotometric method. Additionally, computational docking studies were performed in order to elucidate the interactions occurring with the tyrosinase catalytic site. Our results show that antamanide did not exert any inhibitory activity. On the contrary, the three glycine derivatives AG9, AG6, and AOG9, which differ from each other by the position of a glycine that substitutes phenylalanine in the parent molecule, improving water solubility and flexibility, showed tyrosinase inhibition by spectrophotometric assays. Analytical data were confirmed by computational studies.

GC/MS characterization and computational kinome-wide screening of pomegranate fruit extract identifies key phytochemicals interacting to CDK kinases implicated in acute myeloid leukemia cells

Punica granatum (Pomegranate fruit) and its constituents are proven effective against various cancer types. However, a kinome-wide screening for the active phytochemicals against kinases is not reported. This study aims in validating pomegranate fruit extract (PFE) against acute myeloid leukemia (AML) cells, and computationally identifying the phytochemicals interacting with active kinases. PFE was made with Soxhlet extractor using absolute ethanol. Gas-chromatography-mass spectroscopy (GC-MS) for phytochemical identification and MTT assay for cytotoxicity in AML (THP-1, TF-1 and HL-60) cells were performed. Apoptosis, CDK5 and CDK8 were assessed with flow cytometry. Kinase profiling was performed using In silico kinome screening. GC-MS analysis revealed 38 bioactive phytochemicals in PFE including pyrazoles, aldehydes, phenols, esters, pyranosides, and octadecadienoic acids. The extract inhibited the AML cell proliferations with GI₅₀ values of 195.5 μg/ml, 289.1 μg/ml, and 353.5 μg/ml in THP-1, THP-1, and HL-60 cells, respectively. PFE also exhibited a dose-responsive increase in apoptotic cell populations when treated to the AML cells. Computational screening and modeling predicted three critical constituents, viz., Deoxyartemisinin, 3-Methyl-3-phenyl-3H-indazole, and 8-fluoro-5,6-dimethoxy-3,4-dihydro-2H-naphthalen-1-one of pomegranate extract to interact mainly with cyclin-dependent kinases, including CDK5 and CDK8. Proteinand ligand docking predicted binding energies, and binding pose for top candidate lead molecules. In vitro assay exhibited the anticancer properties of PFE in AML cells. Computational kinome screening predicted top three PFE constituents targeting CDKs which may be responsible for the demonstrated anticancer efficacy of the extract against AML. This hypothesis further aligned with observed efficacy of PFE to inhibit CDK5 and CDK8 in all AML cells tested. PRACTICAL APPLICATIONS: Though Punica granatum (Pomegranate fruit) and its constituents are proven effective against various cancer types, a kinome-wide screening for the active phytochemicals against kinases is not reported. In this study, we have conducted GC/MS characterization of the active phytochemicals of PFE and have performed a kinome-wide screening for all the 38 identified compounds toward 310 active kinases commonly expressed in cancers. These observations warrant isolation and further evaluation of these phytochemicals or their analogues as effective CDK inhibitors against AML proliferation. Further, the computational methods used in this study will throw light on literature for new options of kinome panel screening of active phytochemicals or small molecules.

Computational Prediction of Phosphoinositide Binding to Hyperpolarization-Activated Cyclic-Nucleotide Gated Channels

Protein-lipid interactions are key regulators of ion channel function. Numerous ion channels, including hyperpolarization-activated cyclic-nucleotide gated (HCN) channels have been shown to be regulated by phosphoinositides (PIPs), with important implications in cardiac and neuronal function. Specifically, PIPs have been shown to enhance HCN activation. Using computational approaches, we aim to identify potential binding sites for HCN1-PIP interactions. Computational docking and coarse-grained simulations indicate that PIP binding to HCN1 channels is not well coordinated, but rather occurs over a broad surface of charged residues primarily in the HCN-domain, S2 and S3 helices that can be loosely organized in 2 or 3 overlapping clusters. Thus, PIP-HCN1 interactions are more resembling of electrostatic interactions that occur in myristoylated alanine-rich C kinase substrate (MARCKS) proteins, than the specifically coordinated interactions that occur in pleckstrin homology domains (PH domains) or ion channels such as inward rectifier potassium (Kir) channels. Our results also indicate that phosphatidylinositol (PI) interactions with HCN1 are even lower affinity, explaining why unphosphorylated PI have no effect on HCN1 activation unlike phosphorylated PIPs.

The Non-phosphorylating Glyceraldehyde-3-Phosphate Dehydrogenase GapN Is a Potential New Drug Target in Streptococcus pyogenes

The strict human pathogen Streptococcus pyogenes causes infections of varying severity, ranging from self-limiting suppurative infections to life-threatening diseases like necrotizing fasciitis or streptococcal toxic shock syndrome. Here, we show that the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase GapN is an essential enzyme for S. pyogenes. GapN converts glyceraldehyde 3-phosphate into 3-phosphoglycerate coupled to the reduction of NADP to NADPH. The knock-down of gapN by antisense peptide nucleic acids (asPNA) significantly reduces viable bacterial counts of S. pyogenes laboratory and macrolide-resistant clinical strains in vitro. As S. pyogenes lacks the oxidative part of the pentose phosphate pathway, GapN appears to be the major NADPH source for the bacterium. Accordingly, other streptococci that carry a complete pentose phosphate pathway are not prone to asPNA-based gapN knock-down. Determination of the crystal structure of the S. pyogenes GapN apo-enzyme revealed an unusual cis-peptide in proximity to the catalytic binding site. Furthermore, using a structural modeling approach, we correctly predicted competitive inhibition of S. pyogenes GapN by erythrose 4-phosphate, indicating that our structural model can be used for in silico screening of specific GapN inhibitors. In conclusion, the data provided here reveal that GapN is a potential target for antimicrobial substances that selectively kill S. pyogenes and other streptococci that lack the oxidative part of the pentose phosphate pathway.

High throughput virtual screening and molecular dynamics simulation analysis of phytomolecules against BfmR of Acinetobacter baumannii: anti-virulent drug development campaign

Acinetobacter baumannii is a notorious multidrug resistant bacterium responsible for several hospital acquired infections assisted by its capacity to develop biofilms. A. baumannii BfmR (RstA), a response regulator from the BfmR/S two-component signal transduction system, is the major controller of A. baumannii biofilm development and formation. As a result, BfmR represents a novel target for anti-biofilm treatment against A. baumannii. The discovery of the high-resolution crystal structure of BfmR provides a good chance for computational screening of its probable inhibitors. Therefore, in this study we aim to search new, less toxic, and natural BfmR inhibitors from 8450 phytomolecules available in the Indian Medicinal Plants, Phytochemistry and Therapeutic (IMPPAT) database by analyzing molecular docking against BfmR (PDB ID: 6BR7). Out of these 8450 phytomolecules 6742 molecules were successfully docked with BfmR with the docking score range -6.305 kcal/mol to +5.120 kcal/mol. Structure based-molecular docking (SB-MD) and ADMET (absorption, distribution, metabolism, excretion, & toxicity) profile examination revealed that Norepinephrine, Australine, Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline phytocompounds strongly binds to the active site residues of BfmR. Furthermore, molecular dynamics simulation (MDS) studies for 100 ns and the binding free energy (MM/GBSA) analysis elucidated the binding mechanism of Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline to BfmR. In summary, these phytocompounds seems to have the promising molecules against BfmR, and thus necessitates further verification by both in vitro and in vivo experiments. HighlightsBfmR plays a key role in biofilm development and exopolysaccharide (EPS) synthesis in A. baumannii.Computational approach to search for promising BfmR inhibitors from IMPAAT database.The lead phytomolecules such as Calystegine B3, 7,7 A-Diepialexine, and Alpha-Methylnoradrenaline displayed significant binding with BfmR active site.The outcome of BfmR binding phytomolecules has broadened the scope of hit molecules validation.Communicated by Ramaswamy H. Sarma.

Interaction between miR4749 and Human Serum Albumin as Revealed by Fluorescence, FRET, Atomic Force Spectroscopy and Computational Modelling

The interaction of Human Serum Albumin (HSA) with the microRNA, miR4749, was investigated by Atomic Force Spectrscopy (AFS), static and time-resolved fluorescence spectroscopy and by computational methods. The formation of a HSA/miR4749 complex with an affinity of about 10⁴ M⁻¹ has been assessed through a Stern–Volmer analysis of steady-state fluorescence quenching of the lone Trp residue (Trp214) emission of HSA. Förster Resonance Energy Transfer (FRET) measurements of fluorescence lifetime of the HSA/miR4749 complex were carried out in the absence and in the presence of an acceptor chromophore linked to miR4749. This allowed us to determine a distance of 4.3 ± 0.5 nm between the lone Trp of HSA and the dye bound to miR4749 5p-end. Such a distance was exploited for a screening of the possible binding sites between HSA and miR4749, as predicted by computational docking. Such an approach, further refined by binding free energy calculations, led us to the identification of a consistent model for the structure of the HSA/miR4749 complex in which a positively charged HSA pocket accommodates the negatively charged miRNA molecule. These results designate native HSA as a suitable miRNA carrier under physiological conditions for delivering to appropriate targets.

A Comparative Study between Spanish and British SARS-CoV-2 Variants

The study of the interaction between the SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor is key to understanding binding affinity and stability. In the present report, we sought to investigate the differences between two already sequenced genome variants (Spanish and British) of SARS-CoV-2. Methods: In silico model evaluating the homology, identity and similarity in the genome sequence and the structure and alignment of the predictive spike by computational docking methods. Results: The identity results between the Spanish and British variants of the Spike protein were 28.67%. This close correspondence in the results between the Spanish and British SARS-CoV-2 variants shows that they are very similar (99.99%). The alignment obtained results in four deletions. There were 23 nucleotide substitutions also predicted which could affect the functionality of the proteins produced from this sequence. The interaction between the binding receptor domain from the spike protein and the ACE2 receptor produces some of the mutations found and, therefore, the energy of this ligand varies. However, the estimated antigenicity of the British variant is higher than its Spanish counterpart. Conclusions: Our results indicate that minimal mutations could interfere in the infectivity of the virus due to changes in the fitness between host cell recognition and interaction proteins. In particular, the N501Y substitution, situated in the RBD of the spike of the British variant, might be the reason for its extraordinary infective potential.

Identification of new fisetin analogs as kinase inhibitors: Data on synthesis and anti-skin cancer activities evaluation

This article contains supplemental datasets of the recently published related research article “Synthesis, Inverse Docking-Assisted Identification and in vitro Biological Characterization of Flavonol-based Analogs of Fisetin as c-Kit, CDK2 and mTOR Inhibitors against Melanoma and Non-melanoma Skin Cancers” by Roy et al., [1]. It provides in-depth data not included in the original co-submission on the biophysical, molecular docking, and biological characterization of newly synthesized flavonol-based analogs of fisetin, a natural dietary small molecule with anticancer and anti-inflammatory properties. These synthetic small molecules were investigated as new, potential single and/or multi-kinase inhibitors of the cyclin-dependent kinase-2 (CDK2), receptor tyrosine kinases (c-KITs), and mammalian targets of rapamycin (mTOR) targets, potentially active against melanoma or non-melanoma skin cancers. Furthermore, this data-in-brief article comprises additional sets of results on several aspects of the properties of the dual and multiple kinase inhibitor compounds’ effects that were not presented in the associated article, including the activated targets that are dysregulated in skin cancers; the effects on markers of apoptosis; on colony formation; and in scratch wound healing assays. The study has identified a panel of novel fisetin analogs that are either single- or multi-kinase inhibitors, which may be further developed as active for the treatment of melanoma and non-melanoma skin cancers. The dataset presented herein will be utilized for additional studies aiming to establish a biological platform to steer for predictive and experimental screening of novel flavonoids and analogs in relevant organoids, humanized animal models and in vivo disease models. The present results should also serve as a key stepping-stone towards enabling target-structure-based design, synthesis and initial testing of novel analogs or derivatives of fisetin. The current study may eventually lead to the development of safe, promising and preclinical candidate entities for treatment of skin and other forms of cancers as well as various other human diseases, which can possibly add to the general armamentarium of promising and safe drugs for health promotion.

Known Drugs Identified by Structure-Based Virtual Screening Are Able to Bind Sigma-1 Receptor and Increase Growth of Huntington Disease Patient-Derived Cells

Huntington disease (HD) is a devastating and presently untreatable neurodegenerative disease characterized by progressively disabling motor and mental manifestations. The sigma-1 receptor (σ1R) is a protein expressed in the central nervous system, whose 3D structure has been recently determined by X-ray crystallography and whose agonists have been shown to have neuroprotective activity in neurodegenerative diseases. To identify therapeutic agents against HD, we have implemented a drug repositioning strategy consisting of: (i) Prediction of the ability of the FDA-approved drugs publicly available through the ZINC database to interact with σ1R by virtual screening, followed by computational docking and visual examination of the 20 highest scoring drugs; and (ii) Assessment of the ability of the six drugs selected by computational analyses to directly bind purified σ1R in vitro by Surface Plasmon Resonance and improve the growth of fibroblasts obtained from HD patients, which is significantly impaired with respect to control cells. All six of the selected drugs proved able to directly bind purified σ1R in vitro and improve the growth of HD cells from both or one HD patient. These results support the validity of the drug repositioning procedure implemented herein for the identification of new therapeutic tools against HD.

Structure-Guided Computational Methods Predict Multiple Distinct Binding Modes for Pyrazoloquinolinones in GABAA Receptors

Pyrazoloquinolinones (PQs) are a versatile class of GABAA receptor ligands. It has been demonstrated that high functional selectivity for certain receptor subtypes can be obtained by specific substitution patterns, but so far, no clear SAR rules emerge from the studies. As is the case for many GABAA receptor targeting chemotypes, PQs can interact with distinct binding sites on a given receptor pentamer. In pentamers of αβγ composition, such as the most abundant α1β2γ2 subtype, many PQs are high affinity binders of the benzodiazepine binding site at the extracellular α+/γ2- interfaces. There they display a functionally near silent, flumazenil-like allosteric activity. More recently, interactions with extracellular α+/β- interfaces have been investigated, where strong positive modulation can be steered toward interesting subtype preferences. The most prominent examples are functionally α6-selective PQs. Similar to benzodiazepines, PQs also seem to interact with sites in the transmembrane domain, mainly the sites used by etomidate and barbiturates. This promiscuity leads to potential contributions from multiple sites to net modulation. Developing ligands that interact exclusively with the extracellular α+/β- interfaces would be desired. Correlating functional profiles with binding sites usage is hampered by scarce and heterogeneous experimental data, as shown in our meta-analysis of aggregated published data. In the absence of experimental structures, bound states can be predicted with pharmacophore matching methods and with computational docking. We thus performed pharmacophore matching studies for the unwanted sites, and computational docking for the extracellular α1,6+/β3- interfaces. The results suggest that PQs interact with their binding sites with diverse binding modes. As such, rational design of improved ligands needs to take a complex structure-activity landscape with branches between sub-series of derivatives into account. We present a workflow, which is suitable to identify and explore potential branching points on the structure-activity landscape of any small molecule chemotype.

The AutoDock suite at 30

The AutoDock suite provides a comprehensive toolset for computational ligand docking and drug design and development. The suite builds on 30 years of methods development, including empirical free energy force fields, docking engines, methods for site prediction, and interactive tools for visualization and analysis. Specialized tools are available for challenging systems, including covalent inhibitors, peptides, compounds with macrocycles, systems where ordered hydration plays a key role, and systems with substantial receptor flexibility. All methods in the AutoDock suite are freely available for use and reuse, which has engendered the continued growth of a diverse community of primary users and third-party developers.

Curcumin Conjugates of Non-steroidal Anti-Inflammatory Drugs: Synthesis, Structures, Anti-proliferative Assays, Computational Docking, and Inflammatory Response

In an effort to combine the anti-proliferative effect of CUR-BF2 and CUR compounds with anti-inflammatory benefits of non-steroidal anti-inflammatory drugs (NSAIDs), a library of the bis- and mono-NSAID/CUR-BF2 and NSAID/CUR conjugates were synthesized by coupling flufenamic acid, flurbiprofen, naproxen, indomethacin, and ibuprofen to diversely substituted hydroxy-benzaldehydes via an ester linkage, and by subsequent reaction with acetylacetone-BF2 to form the bis- and the mono-NSAID/CUR-BF2 adducts. Since conversion to NSAID/CUR by the previously developed decomplexation protocol showed limited success, a set of NSAID/CUR conjugates were independently prepared by directly coupling the NSAIDs with parent curcumin. The bis-NSAID/CUR-BF2 and bis-NSAID-CUR hybrids exhibited low cytotoxicity in NCI-60 assay, and in independent cell viability assay on colorectal cancer (CRC) cells (HCT116, HT29, DLD-1, RKO, SW837, CaCo2) and in normal CR cells (CCD841CoN). By contrast, the mono-naproxin and mono-flurbiprofen CUR-BF2 adducts exhibited remarkable anti-proliferative and apoptopic activity in NCI-60 assay most notably against HCT-116 (colon), OVCAR-3 (ovarian), and ACHN (renal) cells. Computational molecular docking calculations showed favorable binding energies to HER2, VEGFR2, BRAF, and Bcl-2 as well as to COX-1 and COX-2, which in several cases exceeded known inhibitors. The main interactions between the ligands and the proteins were hydrophobic, although several hydrogen bonds were also observed. A sub-set of six compounds that had exhibited little or no cytotoxicity were tested for their anti-inflammatory response with THP-1 human macrophages in comparison to parent NSAIDs or parent curcumin.

The small molecule antibody mimic SH7139 targets a family of HLA-DRs expressed by B-cell lymphomas and other solid cancers

Selective high-affinity ligands (SHALs) belong to a novel class of small-molecule cancer therapeutics that function as targeted prodrugs. SH7139, the most advanced of the SHAL drugs designed to bind to a unique β-subunit structural epitope located on HLA-DR10, has exhibited exceptional preclinical efficacy and safety profiles. A comparison of SH7139 and SH7129, a biotin derivative of the drug developed for use as a diagnostic, showed the incorporation of a biotin tag did not alter the SHALs ability to target or kill HLA-DR10 expressing Raji cells. The use of SH7129 in an immuno-histochemical type assay to stain peripheral blood mononuclear cells (PBMCs) obtained from individuals expressing specific HLA-DRB1 alleles has also revealed that in addition to HLA-DR10, seven other more commonly expressed HLA-DRs are targeted by the drug. Computational dockings of the SHAL's recognition ligands to a number of HLA-DR structures explain, in part, why the targeting domains of SH7129 and SH7139 bind to some HLA-DRs but not others. The results also substantiate the selectivity of SH7129 and suggest it may prove useful as a companion diagnostic for pre-screening biopsy samples to identify those patients whose tumours should respond to SH7139 therapy.

Investigation of a Direct Interaction between miR4749 and the Tumor Suppressor p53 by Fluorescence, FRET and Molecular Modeling

The interactions between the DNA binding domain (DBD) of the tumor suppressor p53 and miR4749, characterized by a high sequence similarity with the DNA Response Element (RE) of p53, was investigated by fluorescence spectroscopy combined with computational modeling and docking. Fluorescence quenching experiments witnessed the formation of a specific complex between DBD and miR4749 with an affinity of about 10⁵ M. Förster Resonance Energy Transfer (FRET) allowed us to measure a distance of 3.9 ± 0.3 nm, between the lone tryptophan of DBD and an acceptor dye suitably bound to miR4749. Such information, combined with a computational modeling approach, allowed us to predict possible structures for the DBD-miR4749 complex. A successive docking refinement, complemented with binding free energy calculations, led us to single out a best model for the DBD-miR4749 complex. We found that the interaction of miR4749 involves the DBD L₃ loop and the H₁ helix, close to the Zn-finger motif; with this suggesting that miR4749 could directly inhibit the p53 interaction with DNA. These results might inspire new therapeutic strategies finalized to restore the p53 functional activity.

Short peptides secreted by Bacillus subtilis inhibit the growth of mold on fresh-cut pumpkin (Cucurbita pepo)

BACKGROUND

This study investigates the efficacy of short peptides secreted by Bacillus subtilis for fungal inhibition in fresh-cut pumpkin and for maintaining its shelf life.

RESULTS

Low-molecular-weight filtrate (LC < 1000 Da) of B. subtilis culture (BC) significantly lowered the total number of molds on fresh-cut pumpkin compared with the untreated control and a BC group after storage. Low-molecular-weight filtrate prevented the deterioration of sensory quality in a pumpkin incision, and reduced pectinase activity. It also inhibited the growth of Phytophthora capsici and Penicillium chrysogenum, and the activity of β-1,3-glucan synthase (GS) secreted by both molds. Fifty-seven GS-inhibiting peptides were screened from 95 LC peptides with two to five amino acid residues. The two most potent peptides, AWYW and HWWY, had strongly suppressive effects on the growth of P. capsici and P. chrysogenum.

CONCLUSION

Our study demonstrated that short peptides present in B. subtilis culture can play an important role in the maintenance of fresh-cut pumpkin by suppressing fungal growth. © 2019 Society of Chemical Industry.

Identification of Aethina tumida Kir Channels as Putative Targets of the Bee Venom Peptide Tertiapin Using Structure-Based Virtual Screening Methods

Venoms are comprised of diverse mixtures of proteins, peptides, and small molecules. Identifying individual venom components and their target(s) with mechanism of action is now attainable to understand comprehensively the effectiveness of venom cocktails and how they collectively function in the defense and predation of an organism. Here, structure-based computational methods were used with bioinformatics tools to screen and identify potential biological targets of tertiapin (TPN), a venom peptide from Apis mellifera (European honey bee). The small hive beetle (Aethina tumida (A. tumida)) is a natural predator of the honey bee colony and was found to possess multiple inwardly rectifying K⁺ (Kir) channel subunit genes from a genomic BLAST search analysis. Structure-based virtual screening of homology modelled A. tumida Kir (atKir) channels found TPN to interact with a docking profile and interface “footprint” equivalent to known TPN-sensitive mammalian Kir channels. The results support the hypothesis that atKir channels, and perhaps other insect Kir channels, are natural biological targets of TPN that help defend the bee colony from infestations by blocking K⁺ transport via atKir channels. From these in silico findings, this hypothesis can now be subsequently tested in vitro by validating atKir channel block as well as in vivo TPN toxicity towards A. tumida. This study highlights the utility and potential benefits of screening in virtual space for venom peptide interactions and their biological targets, which otherwise would not be feasible.

Probing Multi-Target Action of Phlorotannins as New Monoamine Oxidase Inhibitors and Dopaminergic Receptor Modulators with the Potential for Treatment of Neuronal Disorders

Modulation of multiple protein targets with a single compound is essential for the effective treatment of central nervous system disorders. In our previous G protein-coupled receptor (GPCR) cell-based study, a selective human monoamine oxidase (hMAO)-A inhibitor, eckol, stimulated activity of dopamine D3 and D4 receptors. This result led to our interest in marine phlorotannin-mediated modulation of hMAO enzymes and related GPCRs in neuronal disorders. Here, we evaluate the multi-target effects of phloroglucinol, phlorofucofuroeckol-A (PFF-A), and dieckol by screening their modulatory activity against hMAO-A and -B and various neuronal GPCRs. Among the tested phlorotannins, PFF-A showed the strongest inhibitory activity against both hMAO isoforms, with higher selectivity toward hMAO-B than hMAO-A. Enzyme kinetics and docking data revealed that PFF-A noncompetitively acts on hMAOs into the alternative binding pocket of enzymes with allosteric functions. In a functional assay for GPCR screening, dieckol and PFF-A exhibited a multi-target combination of D3R/D4R agonism and D1/5HT1A/NK1 antagonism. In particular, they effectively stimulated D3R and D4R, compared to other GPCRs. Docking analysis confirmed that dieckol and PFF-A successfully docked into the conserved active sites of D3R and D4R and interacted with aspartyl and serine residues in the orthosteric binding pockets of the respective receptors. Based on our experimental and computational data, we established the structure-activity relationship between tested phlorotannins and target proteins, including hMAOs and GPCRs. Our current findings suggest that hMAO inhibitors dieckol and PFF-A, major phlorotannins of edible brown algae with multi-action on GPCRs, are potential agents for treatment of psychological disorders and Parkinson's disease.

Synthesis, Biological Evaluation, and In Silico Studies of Novel Aminated Xanthones as Potential p53-Activating Agents

Xanthone scaffold has been regarded as an attractive chemical tool in the search for bioactive molecules with antitumor activity, and in particular two xanthone derivatives, 12-hydroxy-2,2-dimethyl-3,4-dihydro-2H,6H-pyrano [3,2-b]xanthen-6-one (4) and 3,4-dimethoxy-9-oxo-9H-xanthene-1-carbaldehyde (5), were described as a murine double minute 2 (MDM2)-p53 inhibitor and a TAp73 activator, respectively. The xanthone 5 was used as a starting point for the construction of a library of 3,4-dioxygenated xanthones bearing chemical moieties of described MDM2-p53 inhibitors. Eleven aminated xanthones were successfully synthesized and initially screened for their ability to disrupt the MDM2-p53 interaction using a yeast cell-based assay. With this approach, xanthone 37 was identified as a putative p53-activating agent through inhibition of interaction with MDM2. Xanthone 37 inhibited the growth of human colon adenocarcinoma HCT116 cell lines in a p53-dependent manner. The growth inhibitory effect of xanthone 37 was associated with the induction of G1-phase cell cycle arrest and increased protein expression levels of p53 transcriptional targets. These results demonstrated the potential usefulness of coupling amine-containing structural motifs of known MDM2-p53 disruptors into a 3,4-dioxygenated xanthone scaffold in the design of novel and potent p53 activators with antitumor activity and favorable drug-like properties. Moreover, in silico docking studies were performed in order to predict the binding poses and residues involved in the potential MDM2-p53 interaction.

Structural and Computational Characterization of Disease-Related Mutations Involved in Protein-Protein Interfaces

One of the known potential effects of disease-causing amino acid substitutions in proteins is to modulate protein-protein interactions (PPIs). To interpret such variants at the molecular level and to obtain useful information for prediction purposes, it is important to determine whether they are located at protein-protein interfaces, which are composed of two main regions, core and rim, with different evolutionary conservation and physicochemical properties. Here we have performed a structural, energetics and computational analysis of interactions between proteins hosting mutations related to diseases detected in newborn screening. Interface residues were classified as core or rim, showing that the core residues contribute the most to the binding free energy of the PPI. Disease-causing variants are more likely to occur at the interface core region rather than at the interface rim (p < 0.0001). In contrast, neutral variants are more often found at the interface rim or at the non-interacting surface rather than at the interface core region. We also found that arginine, tryptophan, and tyrosine are over-represented among mutated residues leading to disease. These results can enhance our understanding of disease at molecular level and thus contribute towards personalized medicine by helping clinicians to provide adequate diagnosis and treatments.

Epimagnolin targeting on an active pocket of mammalian target of rapamycin suppressed cell transformation and colony growth of lung cancer cells

Mammalian target of rapamycin (mTOR) has a pivotal role in carcinogenesis and cancer cell proliferation in diverse human cancers. In this study, we observed that epimagnolin, a natural compound abundantly found in Shin-Yi, suppressed cell proliferation by inhibition of epidermal growth factor (EGF)-induced G1/S cell-cycle phase transition in JB6 Cl41 cells. Interestingly, epimagnolin suppressed EGF-induced Akt phosphorylation strongly at Ser473 and weakly at Thr308 without alteration of phosphorylation of MAPK/ERK kinases (MEKs), extracellular signal-regulated kinase (ERKs), and RSK1, resulting in abrogation of the phosphorylation of GSK3β at Ser9 and p70S6K at Thr389. Moreover, we found that epimagnolin suppressed c-Jun phosphorylation at Ser63/73, resulting in the inhibition of activator protein 1 (AP-1) transactivation activity. Computational docking indicated that epimagnolin targeted an active pocket of the mTOR kinase domain by forming three hydrogen bonds and three hydrophobic interactions. The prediction was confirmed by using in vitro kinase and adenosine triphosphate-bead competition assays. The inhibition of mTOR kinase activity resulted in the suppression of anchorage-independent cell transformation. Importantly, epimagnolin efficiently suppressed cell proliferation and anchorage-independent colony growth of H1650 rather than H460 lung cancer cells with dependency of total and phosphorylated protein levels of mTOR and Akt. Inhibitory signaling of epimagnolin on cell proliferation of lung cancer cells was observed mainly in mTOR-Akt-p70S6K and mTOR-Akt-GSK3β-AP-1, which was similar to that shown in JB6 Cl41 cells. Taken together, our results indicate that epimagnolin potentiates as chemopreventive or therapeutic agents by direct active pocket targeting of mTOR kinase, resulting in sensitizing cancer cells harboring enhanced phosphorylation of the mTORC2-Akt-p70S6k signaling pathway.

Arctigenin inhibits STAT3 and exhibits anticancer potential in human triple-negative breast cancer therapy

Triple-negative breast cancers (TNBCs) are the most aggressive and hard-to-treat breast tumors with poor prognosis, and exploration for novel therapeutic drugs is impending. Arctigenin (Atn), a bioactive lignan isolated from seeds of Arctium lappa L, has been reported to inhibit many cancer types; however, the effect of Atn on TNBC remains unclear. In this study, we demonstrated that Atn decreased proliferation, and induced apoptosis in TNBC cells. Furthermore, we explored the underlying mechanism of Atn inhibition on TNBC cells. Computational docking and affinity assay showed that Atn bound to the SH2 domain of STAT3. Atn inhibited STAT3 binding to genomic DNA by disrupting hydrogen bond linking between DNA and STAT3. In addition, Atn augmented Taxotere^®-induced TNBC cell cytotoxicity. TNBC xenograft tests also confirmed the antitumor effect of Atn in vivo. These characteristics render Atn as a promising candidate drug for further development and for designing new effective STAT3 inhibitors.

Hot-spot analysis for drug discovery targeting protein-protein interactions

INTRODUCTION

Protein-protein interactions are important for biological processes and pathological situations, and are attractive targets for drug discovery. However, rational drug design targeting protein-protein interactions is still highly challenging. Hot-spot residues are seen as the best option to target such interactions, but their identification requires detailed structural and energetic characterization, which is only available for a tiny fraction of protein interactions. Areas covered: In this review, the authors cover a variety of computational methods that have been reported for the energetic analysis of protein-protein interfaces in search of hot-spots, and the structural modeling of protein-protein complexes by docking. This can help to rationalize the discovery of small-molecule inhibitors of protein-protein interfaces of therapeutic interest. Computational analysis and docking can help to locate the interface, molecular dynamics can be used to find suitable cavities, and hot-spot predictions can focus the search for inhibitors of protein-protein interactions. Expert opinion: A major difficulty for applying rational drug design methods to protein-protein interactions is that in the majority of cases the complex structure is not available. Fortunately, computational docking can complement experimental data. An interesting aspect to explore in the future is the integration of these strategies for targeting PPIs with large-scale mutational analysis.

Exploring Heteroaryl-pyrazole Carboxylic Acids as Human Carbonic Anhydrase XII Inhibitors

We report the synthesis, biological evaluation, and structural study of a series of substituted heteroaryl-pyrazole carboxylic acid derivatives. These compounds have been developed as inhibitors of specific isoforms of carbonic anhydrase (CA), with potential as prototypes of a new class of chemotherapeutics. Both X-ray crystallography and computational modeling provide insights into the CA inhibition mechanism. Results indicate that this chemotype produces an indirect interference with the zinc ion, thus behaving differently from other related nonclassical inhibitors. Among the tested compounds, 2c with K_i = 0.21 μM toward hCA XII demonstrated significant antiproliferative activity against hypoxic tumor cell lines. Taken together, the results thus provide the basis of structural determinants for the development of novel anticancer agents.

Buried Hydrogen Bond Interactions Contribute to the High Potency of Complement Factor D Inhibitors

Aberrant activation of the complement system is associated with diseases, including paroxysmal nocturnal hemoglobinuria and age-related macular degeneration. Complement factor D is the rate-limiting enzyme for activating the alternative pathway in the complement system. Recent development led to a class of potent amide containing pyrrolidine derived factor D inhibitors. Here, we used biochemical enzymatic and biolayer interferometry assays to demonstrate that the amide group improves the inhibitor potency by more than 80-fold. Our crystal structures revealed buried hydrogen bond interactions are important. Molecular orbital analysis from quantum chemistry calculations dissects the chemical groups participating in these interactions. Free energy calculation supports the differential contributions of the amide group to the binding affinities of these inhibitors. Cell-based hemolysis assay confirmed these compounds inhibit factor D mediated complement activation via the alternative pathway. Our study highlights the important interactions contributing to the high potency of factor D inhibitors reported recently.

Covalent docking using autodock: Two-point attractor and flexible side chain methods

We describe two methods of automated covalent docking using Autodock4: the two-point attractor method and the flexible side chain method. Both methods were applied to a training set of 20 diverse protein-ligand covalent complexes, evaluating their reliability in predicting the crystallographic pose of the ligands. The flexible side chain method performed best, recovering the pose in 75% of cases, with failures for the largest inhibitors tested. Both methods are freely available at the AutoDock website (http://autodock.scripps.edu).

Prediction of noncovalent Drug/DNA interaction using computational docking models: studies with over 1350 launched drugs

Noncovalent chemical/DNA interactions, for example, intercalation and groove-binding, may be more important to genomic integrity than previously appreciated, and there may very well be genotoxic consequences of that binding. It is of importance, then, to develop methods allowing a determination or prediction of such interactions. This would have particular utility in the pharmaceutical industry where genotoxicity is, for the most part, disallowed in new drug entities. We have previously used DNA docking simulations to assess if molecules had structure and charge characteristics which could accommodate noncovalent binding via, for example, electrostatic/hydrogen bonding. We here extend those earlier studies by examining a series of over 1,350 "launched" drugs for ability to noncovalently bind 10 different DNA sequences using two computational programs: Autodock and Surflex. These drugs were also evaluated for binding to the crystallographic ATP-binding site of human topoisomerase II. The results obtained clearly demonstrate multiple series of noncovalent DNA binding structure activity relationships which would not have been predicted based on cursory structural examination. Many drugs within these series are genotoxic although not via any commonly recognized structural covalent alerts. The present studies confirm previously implicated features such as N-dialkyl groups and specific N-aryl ketones as potential genotoxic chemical moieties acting through noncovalent mechanisms. These initial studies provide considerable evidence that DNA intercalation may be an important, largely overlooked, source of drug-induced genotoxicity and further suggest involvement of topoisomerase in that genotoxicity.

Designed Trpzip-3 β-Hairpin Inhibits Amyloid Formation in Two Different Amyloid Systems

The trpzip peptides are small, monomeric, and extremely stable β-hairpins that have become valuable tools for studying protein folding. Here, we show that trpzip-3 inhibits aggregation in two very different amyloid systems: transthyretin and Aβ(1-42). Interestingly, Trp → Leu mutations renders the peptide ineffective against transthyretin, but Aβ inhibition remains. Computational docking was used to predict the interactions between trpzip-3 and transthyretin, suggesting that inhibition occurs via binding to the outer region of the thyroxine-binding site, which is supported by dye displacement experiments.

DockTrina: docking triangular protein trimers

In spite of the abundance of oligomeric proteins within a cell, the structural characterization of protein-protein interactions is still a challenging task. In particular, many of these interactions involve heteromeric complexes, which are relatively difficult to determine experimentally. Hence there is growing interest in using computational techniques to model such complexes. However, assembling large heteromeric complexes computationally is a highly combinatorial problem. Nonetheless the problem can be simplified greatly by considering interactions between protein trimers. After dimers and monomers, triangular trimers (i.e. trimers with pair-wise contacts between all three pairs of proteins) are the most frequently observed quaternary structural motifs according to the three-dimensional (3D) complex database. This article presents DockTrina, a novel protein docking method for modeling the 3D structures of nonsymmetrical triangular trimers. The method takes as input pair-wise contact predictions from a rigid body docking program. It then scans and scores all possible combinations of pairs of monomers using a very fast root mean square deviation test. Finally, it ranks the predictions using a scoring function which combines triples of pair-wise contact terms and a geometric clash penalty term. The overall approach takes less than 2 min per complex on a modern desktop computer. The method is tested and validated using a benchmark set of 220 bound and seven unbound protein trimer structures. DockTrina will be made available at http://nano-d.inrialpes.fr/software/docktrina.

Crystal structure of malonyl CoA-Acyl carrier protein transacylase from Xanthomanous oryzae pv. oryzae and its proposed binding with ACP

Xanthomonas oryzae pv. oryzae (Xoo) is a plant bacterial pathogen that causes bacterial blight (BB) disease, resulting in serious production losses of rice. The crystal structure of malonyl CoA-acyl carrier protein transacylase (XoMCAT), encoded by the gene fabD (Xoo0880) from Xoo, was determined at 2.3 Å resolution in complex with N-cyclohexyl-2-aminoethansulfonic acid. Malonyl CoA-acyl carrier protein transacylase transfers malonyl group from malonyl CoA to acyl carrier protein (ACP). The transacylation step is essential in fatty acid synthesis. Based on the rationale, XoMCAT has been considered as a target for antibacterial agents against BB. Protein-protein interaction between XoMCAT and ACP was also extensively investigated using computational docking, and the proposed model revealed that ACP bound to the cleft between two XoMCAT subdomains.

Y95 and E444 interaction required for high-affinity S-citalopram binding in the human serotonin transporter

The human serotonin (5-hydroxytryptamine, 5-HT) transporter (hSERT) is responsible for the reuptake of 5-HT following synaptic release, as well as for import of the biogenic amine into several non-5-HT synthesizing cells including platelets. The antidepressant citalopram blocks SERT and thereby inhibits the transport of 5-HT. To identify key residues establishing high-affinity citalopram binding, we have built a comparative model of hSERT and Drosophila melanogaster SERT (dSERT) based on the Aquifex aeolicus leucine transporter (LeuT(Aa)) crystal structure. In this study, citalopram has been docked into the homology model of hSERT and dSERT using RosettaLigand. Our models reproduce the differential binding affinities for the R- and S-isomers of citalopram in hSERT and the impact of several hSERT mutants. Species-selective binding affinities for hSERT and dSERT also can be reproduced. Interestingly, the model predicts a hydrogen bond between E444 in transmembrane domain 8 (TM8) and Y95 in TM1 that places Y95 in a downward position, thereby removing Y95 from a direct interaction with S-citalopram. Mutation of E444D results in a 10-fold reduced binding affinity for S-citalopram, supporting the hypothesis that Y95 and E444 form a stabilizing interaction in the S-citalopram/hSERT complex.

Computational docking of antibody-antigen complexes, opportunities and pitfalls illustrated by influenza hemagglutinin

Antibodies play an increasingly important role in both basic research and the pharmaceutical industry. Since their efficiency depends, in ultimate analysis, on their atomic interactions with an antigen, studying such interactions is important to understand how they function and, in the long run, to design new molecules with desired properties. Computational docking, the process of predicting the conformation of a complex from its separated components, is emerging as a fast and affordable technique for the structural characterization of antibody-antigen complexes. In this manuscript, we first describe the different computational strategies for the modeling of antibodies and docking of their complexes, and then predict the binding of two antibodies to the stalk region of influenza hemagglutinin, an important pharmaceutical target. The purpose is two-fold: on a general note, we want to illustrate the advantages and pitfalls of computational docking with a practical example, using different approaches and comparing the results to known experimental structures. On a more specific note, we want to assess if docking can be successful in characterizing the binding to the same influenza epitope of other antibodies with unknown structure, which has practical relevance for pharmaceutical and biological research. The paper clearly shows that some of the computational docking predictions can be very accurate, but the algorithm often fails to discriminate them from inaccurate solutions. It is of paramount importance, therefore, to use rapidly obtained experimental data to validate the computational results.

Computational modeling of the potential interactions of the proteasome beta5 subunit and catechol-O-methyltransferase-resistant EGCG analogs

(-)-Epigallocatechin gallate [(-)-EGCG] has been implicated in cancer chemoprevention and has been shown as an inhibitor of tumor proteasomal chymotrypsin-like activity in vitro and in vivo. However, EGCG is subjected to rapid biotransforming modifications such as methylation by catechol-Omicron-methyltransferase (COMT) that limits its action. We recently reported that structure 7, an EGCG analog which should be resistant to COMT-mediated methylation and inactivation in cells, was able to inhibit the activity of purified 20S proteasome and cellular 26S proteasome. However, the involved molecular mechanism is unknown. Herein, we applied computational solution to understand the possible interaction between EGCG analogs including structure 7 and the proteasome beta5 subunit which is responsible for the chymotrypsin-like activity. We report that the ester carbonyls at C2 and C3 carbon atoms may be the active sites for nucleophilic attack in structure 7 and 5. Equally spaced carbon atoms in COMT-resistant structure 7 give more stable conformation and lower docked free energy than other EGCG analogs. The absence of a second gallate group in structure 16 and 21 significantly decreases the ability to inhibit the proteasome.

Identification of novel urease inhibitors by high-throughput virtual and in vitro screening

Ureases are important in both agriculture and human health. Bacterial ureases are directly involved in many farm-field problems and pathological conditions. Here, we report a structure-based virtual screening of an in-house compound bank of about 6000 molecular entities by computational docking and binding free energy calculations followed by in vitro screening. Applied protocol leads to the identification of novel urease inhibitors, which can serve as starting points for structural optimization.

Binding site on human immunoglobulin G for the affinity ligand HWRGWV

Affinity ligand HWRGWV has demonstrated the ability to isolate human immunoglobulin G (hIgG) from mammalian cell culture media. The ligand specifically binds hIgG through its Fc portion. This work shows that deglycosylation of hIgG has no influence on its binding to the HWRGWV ligand and the ligand does not compete with Protein A or Protein G in binding hIgG. It is suggested by the mass spectrometry (MS) data and docking simulation that HWRGWV binds to the pFc portion of hIgG and interacts with the amino acids in the loop Ser383-Asn389 (SNGQPEN) located in the C(H)3 domain. Subsequent modeling has suggested a possible three-dimensional minimized solution structure for the interaction of hIgG and the HWRGWV ligand. The results support the fact that a peptide as small as a hexamer can have specific interactions with large proteins such as hIgG.