Systematic exploration of multiple drug binding sites

Synthesis and estrogenic activity of BODIPY-labeled estradiol conjugates

Novel BODIPY-estradiol conjugates have been synthesized by selecting position C-3-O for labeling. The conjugation strategy was based on Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) or etherification. Estradiol derivatives used as azide partners bearing an ω-azidoalkyl function through C4-C8-long linkers have been prepared. CuAAC reactions of estradiol azides with BODIPY alkyne furnished fluorescent 3-O-labeled conjugates bearing the triazole ring as a coupling moiety. Williamson etherifications of 3-O-(ω-bromoalkyl)-17β-estradiol derivatives with BODIPY-OH resulted in labeled conjugates connected with an ether moiety. Interactions of the conjugates with estrogen receptor (ER) were investigated using molecular docking calculations in comparison with estradiol. The conjugates occupied both the classical and alternative binding sites on human ERα, with slightly lower binding affinity to references estradiol and diethystilbestrol. All compounds have displayed reasonable estrogenic activity. They increased the proliferation of ER-positive breast cancer cell line MCF7 contrary to ER-negative SKBR-3 cell line. The most potent compound 13a induced the transcriptional activity of ER in dose-dependent manner in dual luciferase recombinant reporter model and increased progesterone receptor's expression, proving the retained estrogenic activity. The fluorescence of candidate compound 13a co-localised with the ERα. The newly synthesized labeled compounds might serve as good starting point for further development of fluorescent probes for modern biological applications. In addition to studying steroid uptake and transport in cells, e.g. in the processes of biodegradation of estrogen-hormones micropollutants, they could also be utilized in examination of estrogen-binding proteins.

Efficient Refinement of Complex Structures of Flexible Histone Peptides Using Post-Docking Molecular Dynamics Protocols

Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To address this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In this present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32% improvement over the docked structures in terms of the change in root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements are also discussed to help with their implementation in future methods and applications.

Brassica oleracea L. var. italica Aquaporin Reconstituted Proteoliposomes as Nanosystems for Resveratrol Encapsulation

Aquaporins (AQPs), membrane proteins responsible for facilitating water transport, found in plant membrane vesicles (MV), have been related to the functionality and stability of MV. We focused on AQPs obtained from broccoli, as they show potential for biotechnological applications. To gain further insight into the role of AQPs in MV, we describe the heterologous overexpression of two broccoli AQPs (BoPIP1;2 and BoPIP2;2) in Pichia pastoris, resulting in their purification with high yield (0.14 and 0.99 mg per gram cells for BoPIP1;2 and BoPIP2;2). We reconstituted AQPs in liposomes to study their functionality, and the size of proteoliposomes did not change concerning liposomes. BoPIP2;2 facilitated water transport, which was preserved for seven days at 4 °C and at room temperature but not at 37 °C. BoPIP2;2 was incorporated into liposomes to encapsulate a resveratrol extract, resulting in increased entrapment efficiency (EE) compared to conventional liposomes. Molecular docking was utilized to identify binding sites in PIP2s for resveratrol, highlighting the role of aquaporins in the improved EE. Moreover, interactions between plant AQP and human integrin were shown, which may increase internalization by the human target cells. Our results suggest AQP-based alternative encapsulation systems can be used in specifically targeted biotechnological applications.

Construction of Histone-Protein Complex Structures by Peptide Growing

The structures of histone complexes are master keys to epigenetics. Linear histone peptide tails often bind to shallow pockets of reader proteins via weak interactions, rendering their structure determination challenging. In the present study, a new protocol, PepGrow, is introduced. PepGrow uses docked histone fragments as seeds and grows the full peptide tails in the reader-binding pocket, producing atomic-resolution structures of histone-reader complexes. PepGrow is able to handle the flexibility of histone peptides, and it is demonstrated to be more efficient than linking pre-docked peptide fragments. The new protocol combines the advantages of popular program packages and allows fast generation of solution structures. AutoDock, a force-field-based program, is used to supply the docked peptide fragments used as structural seeds, and the building algorithm of Modeller is adopted and tested as a peptide growing engine. The performance of PepGrow is compared to ten other docking methods, and it is concluded that in situ growing of a ligand from a seed is a viable strategy for the production of complex structures of histone peptides at atomic resolution.

Inhibition and structure-activity relationship of dietary flavones against three Loop 1-type human gut microbial β-glucuronidases

Gut microbial β-glucuronidases (GUSs) inhibition is a new approach for managing some diseases and medication therapy. However, the structural and functional complexity of GUSs have posed tremendous challenges to discover specific or broad-spectrum GUSs inhibitors using Escherichia coli GUS (EcoGUS) alone. This study first assessed the effects of twenty-one dietary flavones employing three Loop 1-type GUSs of different taxonomic origins, which were considered to be the main GUSs involved in deglucuronidation of small molecules, on p-nitrophenyl-β-D-glucuronide hydrolysis and a structure-activity relationship is preliminarily proposed based on both in vitro assays and a docking study with representative compounds. EcoGUS and Staphylococcus pasteuri GUS showed largely similar inhibition propensities with potencies positively correlating with the total hydroxyl groups and those at ring B of flavones, while docking results revealed strong interactions developed via ring A and/or C. Streptococcus agalactiae GUS (SagaGUS) exhibited distinct inhibition propensities, displaying late-onset inhibition and steep dose-response profiles with most tested compounds. The α-helix in loop 1 region of SagaGUS which causes spatial hindrance but offers a hydrophobic surface for contacting with the carbonyl group on ring C of flavones is believed to be essential for the allosteric inhibition of SagaGUS. Taken together, the study with a series of flavones revealed varied preferences for GUSs belonging to the same Loop 1-type, highlighting the necessity of adopting multi-GUSs instead of EcoGUS alone for screening broad-spectrum GUSs inhibitors or tailoring the inhibition based on specific GUS structure.

The Cholesterol-Modulating Effect of the New Herbal Medicinal Recipe from Yellow Vine (Coscinium fenestratum (Goetgh.)), Ginger (Zingiber officinale Roscoe.), and Safflower (Carthamus tinctorius L.) on Suppressing PCSK9 Expression to Upregulate LDLR Expression in HepG2 Cells

PCSK9 is a promising target for developing novel cholesterol-lowering drugs. We developed a recipe that combined molecular docking, GC-MS/MS, and real-time PCR to identify potential PCSK9 inhibitors for herb ratio determination. Three herbs, Carthamus tinctorius, Coscinium fenestratum, and Zingiber officinale, were used in this study. This work aimed to evaluate cholesterol-lowering through a PCSK9 inhibitory mechanism of these three herbs for defining a suitable ratio. Chemical constituents were identified using GC-MS/MS. The PCSK9 inhibitory potential of the compounds was determined using molecular docking, real-time PCR, and Oil red O staining. It has been shown that most of the active compounds of C. fenestratum and Z. officinale inhibit PCSK9 when extracted with water, and C. fenestratum has been shown to yield tetraacetyl-d-xylonic nitrile (27.92%) and inositol, 1-deoxy-(24.89%). These compounds could inhibit PCSK9 through the binding of 6 and 5 hydrogen bonds, respectively, while the active compound in Z. officinale is 2-Formyl-9-[.beta.-d-ribofuranosyl] hypoxanthine (4.37%) inhibits PCSK9 by forming 8 hydrogen bonds. These results suggest that a recipe comprising three parts C. fenestratum, two parts Z. officinale, and one part C. tinctorius is a suitable herbal ratio for reducing lipid levels in the bloodstream through a PCSK9 inhibitory mechanism.

Binding Networks Identify Targetable Protein Pockets for Mechanism-Based Drug Design

The human genome codes only a few thousand druggable proteins, mainly receptors and enzymes. While this pool of available drug targets is limited, there is an untapped potential for discovering new drug-binding mechanisms and modes. For example, enzymes with long binding cavities offer numerous prerequisite binding sites that may be visited by an inhibitor during migration from a bulk solution to the destination site. Drug design can use these prerequisite sites as new structural targets. However, identifying these ephemeral sites is challenging. Here, we introduce a new method called NetBinder for the systematic identification and classification of prerequisite binding sites at atomic resolution. NetBinder is based on atomistic simulations of the full inhibitor binding process and provides a networking framework on which to select the most important binding modes and uncover the entire binding mechanism, including previously undiscovered events. NetBinder was validated by a study of the binding mechanism of blebbistatin (a potent inhibitor) to myosin 2 (a promising target for cancer chemotherapy). Myosin 2 is a good test enzyme because, like other potential targets, it has a long internal binding cavity that provides blebbistatin with numerous potential prerequisite binding sites. The mechanism proposed by NetBinder of myosin 2 structural changes during blebbistatin binding shows excellent agreement with experimentally determined binding sites and structural changes. While NetBinder was tested on myosin 2, it may easily be adopted to other proteins with long internal cavities, such as G-protein-coupled receptors or ion channels, the most popular current drug targets. NetBinder provides a new paradigm for drug design by a network-based elucidation of binding mechanisms at an atomic resolution.

Exploration of Somatostatin Binding Mechanism to Somatostatin Receptor Subtype 4

Somatostatin (also named as growth hormone-inhibiting hormone or somatotropin release-inhibiting factor) is a regulatory peptide important for the proper functioning of the endocrine system, local inflammatory reactions, mood and motor coordination, and behavioral responses to stress. Somatostatin exerts its effects via binding to G-protein-coupled somatostatin receptors of which the fourth subtype (SSTR4) is a particularly important receptor mediating analgesic, anti-inflammatory, and anti-depressant effects without endocrine actions. Thus, SSTR4 agonists are promising drug candidates. Although the knowledge of the atomic resolution-binding modes of SST would be essential for drug development, experimental elucidation of the structures of SSTR4 and its complexes is still awaiting. In the present study, structures of the somatostatin–SSTR4 complex were produced using an unbiased, blind docking approach. Beyond the static structures, the binding mechanism of SST was also elucidated in the explicit water molecular dynamics (MD) calculations, and key binding modes (external, intermediate, and internal) were distinguished. The most important residues on both receptor and SST sides were identified. An energetic comparison of SST binding to SSTR4 and 2 offered a residue-level explanation of receptor subtype selectivity. The calculated structures show good agreement with available experimental results and indicate that somatostatin binding is realized via prerequisite binding modes and an induced fit mechanism. The identified binding modes and the corresponding key residues provide useful information for future drug design targeting SSTR4.

Unveiling the Possible Oryzalin-Binding Site in the α-Tubulin of Toxoplasma gondii

Dinitroaniline derivatives have been widely used as herbicidal agents to control weeds and grass. Previous studies demonstrated that these compounds also exhibit good antiparasitic activity against some protozoan parasites. Oryzalin (ORY), a representative dinitroaniline derivative, exerts its antiprotozoal activity against Toxoplasma gondii by inhibiting the microtubule polymerization process. Moreover, the identification of ORY-resistant T. gondii lines obtained by chemical mutagenesis confirmed that this compound binds selectively to α-tubulin. Based on experimental information reported so far and a multiple sequence analysis carried out in this work, we propose that the pironetin (PIR) site is the potential ORY-binding site. Therefore, we employed state-of-the-art computational approaches to characterize the interaction profile of ORY at the proposed site in the α-tubulin of T. gondii. An exhaustive search for other possible binding sites was performed using the Wrap "N" Shake method, which showed that ORY exhibits highest stability and affinity for the PIR site. Moreover, our molecular dynamics simulations revealed that the dipropylamine substituent of ORY interacts with a hydrophobic pocket, while the sulfonamide group formed hydrogen bonds with water molecules at the site entrance. Overall, our results suggest that ORY binds to the PIR site on the α-tubulin of the protozoan parasite T. gondii. This information will be very useful for designing less toxic and more potent antiprotozoal agents.

Insights into the binding of morin to human γD-crystallin

Crystallin aggregation in the eye lens is one of the leading causes of cataract formation. The increase in the human γD-crystallin (HγDC) aggregation propensity has been associated with the oligomerization of its partially folded and fully unfolded structure. A recent study demonstrated that the binding of flavonoid morin (MOR) to HγDC inhibits the fibrillation of this protein. In this work, we carry out an exhaustive search for the possible binding site of MOR on HγDC by combining an ensemble docking approach with the Wrap 'N' Shake protocol. In agreement with previous results, we found a potential MOR-binding site in the cleft formed between the N-terminal and C-terminal domains of HγDC. MOR preference for the cleft residues was observed even with the interface-opened intermediate conformers of HγDC. In addition, metadynamics simulations were carried out to corroborate the stabilizing activity of MOR on HγDC structure and to identify the structural regions implicated during the unfolding inhibition. Overall, this study provides relevant insights into the identification of new HγDC aggregation inhibitors.

Prerequisite Binding Modes Determine the Dynamics of Action of Covalent Agonists of Ion Channel TRPA1

Transient receptor potential ankyrin 1 (TRPA1) is a transmembrane protein channeling the influx of calcium ions. As a polymodal nocisensor, TRPA1 can be activated by thermal, mechanical stimuli and a wide range of chemically damaging molecules including small volatile environmental toxicants and endogenous algogenic lipids. After activation by such compounds, the ion channel opens up, its central pore widens allowing calcium influx into the cytosol inducing signal transduction pathways. Afterwards, the calcium influx desensitizes irritant evoked responses and results in an inactive state of the ion channel. Recent experimental determination of structures of apo and holo forms of TRPA1 opened the way towards the design of new agonists, which can activate the ion channel. The present study is aimed at the elucidation of binding dynamics of agonists using experimental structures of TRPA1-agonist complexes at the atomic level applying molecular docking and dynamics methods accounting for covalent and non-covalent interactions. Following a test of docking methods focused on the final, holo structures, prerequisite binding modes were detected involving the apo forms. It was shown how reversible interactions with prerequisite binding sites contribute to structural changes of TRPA1 leading to covalent bonding of agonists. The proposed dynamics of action allowed a mechanism-based forecast of new, druggable binding sites of potent agonists.

Determination of Ligand Binding Modes in Hydrated Viral Ion Channels to Foster Drug Design and Repositioning

Target-based design and repositioning are mainstream strategies of drug discovery. Numerous drug design and repositioning projects have been launched to fight the ongoing COVID-19 pandemic. The resulting drug candidates have often failed due to the misprediction of their target-bound structures. The determination of water positions of such structures is particularly challenging due to the large number of possible drugs and the diversity of their hydration patterns. To answer this challenge and help correct predictions, we introduce a new protocol HydroDock, which can build hydrated drug-target complexes from scratch. HydroDock requires only the dry target and drug structures and produces their complexes with appropriately positioned water molecules. As a test application of the protocol, we built the structures of amantadine derivatives in complex with the influenza M2 transmembrane ion channel. The repositioning of amantadine derivatives from this influenza target to the SARS-CoV-2 envelope protein was also investigated. Excellent agreement was observed between experiments and the structures determined by HydroDock. The atomic resolution complex structures showed that water plays a similar role in the binding of amphipathic amantadine derivatives to transmembrane ion channels of both influenza A and SARS-CoV-2. While the hydrophobic regions of the channels capture the bulky hydrocarbon group of the ligand, the surrounding waters direct its orientation parallel with the axes of the channels via bridging interactions with the ionic ligand head. As HydroDock supplied otherwise undetermined structural details, it can be recommended to improve the reliability of future design and repositioning of antiviral drug candidates and many other ligands with an influence of water structure on their mechanism of action.

Finding Druggable Sites in Proteins Using TACTICS

Structure-based drug discovery efforts require knowledge of where drug-binding sites are located on target proteins. To address the challenge of finding druggable sites, we developed a machine-learning algorithm called TACTICS (trajectory-based analysis of conformations to identify cryptic sites), which uses an ensemble of molecular structures (such as molecular dynamics simulation data) as input. First, TACTICS uses k-means clustering to select a small number of conformations that represent the overall conformational heterogeneity of the data. Then, TACTICS uses a random forest model to identify potentially bindable residues in each selected conformation, based on protein motion and geometry. Lastly, residues in possible binding pockets are scored using fragment docking. As proof-of-principle, TACTICS was applied to the analysis of simulations of the SARS-CoV-2 main protease and methyltransferase and the Yersinia pestis aryl carrier protein. Our approach recapitulates known small-molecule binding sites and predicts the locations of sites not previously observed in experimentally determined structures. The TACTICS code is available at https://github.com/Albert-Lau-Lab/tactics_protein_analysis.

Capturing Protein-Ligand Recognition Pathways in Coarse-Grained Simulation

Protein-ligand recognition is dynamic and complex. A key approach in deciphering the mechanism underlying the recognition process is to capture the kinetic process of the ligand in its act of binding to its designated protein cavity. Toward this end, ultralong all-atom molecular dynamics simulation has recently emerged as a popular method of choice because of its ability to record these events at high spatial and temporal resolution. However, success via this route comes at an exorbitant computational cost. Herein, we demonstrate that coarse-grained models of the protein, when systematically optimized to maintain its tertiary fold, can capture the complete process of spontaneous protein-ligand binding from bulk media to the cavity at crystallographic precision and within wall clock time that is orders of magnitude shorter than that of all-atom simulations. The exhaustive sampling of ligand exploration in protein and solvent, harnessed by coarse-grained simulation, leads to elucidation of new ligand recognition pathways and discovery of non-native binding poses.

Molecular Structure, Binding Affinity, and Biological Activity in the Epigenome

Development of valid structure–activity relationships (SARs) is a key to the elucidation of pathomechanisms of epigenetic diseases and the development of efficient, new drugs. The present review is based on selected methodologies and applications supplying molecular structure, binding affinity and biological activity data for the development of new SARs. An emphasis is placed on emerging trends and permanent challenges of new discoveries of SARs in the context of proteins as epigenetic drug targets. The review gives a brief overview and classification of the molecular background of epigenetic changes, and surveys both experimental and theoretical approaches in the field. Besides the results of sophisticated, cutting edge techniques such as cryo-electron microscopy, protein crystallography, and isothermal titration calorimetry, examples of frequently used assays and fast screening techniques are also selected. The review features how different experimental methods and theoretical approaches complement each other and result in valid SARs of the epigenome.

A Comprehensive Computational Platform to Guide Drug Development Using Graph-Based Signature Methods

High-throughput computational techniques have become invaluable tools to help increase the overall success, process efficiency, and associated costs of drug development. By designing ligands tailored to specific protein structures in a disease of interest, an understanding of molecular interactions and ways to optimize them can be achieved prior to chemical synthesis. This understanding can help direct crucial chemical and biological experiments by maximizing available resources on higher quality leads. Moreover, predicting molecular binding affinity within specific biological contexts, as well as ligand pharmacokinetics and toxicities, can aid in filtering out redundant leads early on within the process. We describe a set of computational tools which can aid in drug discovery at different stages, from hit identification (EasyVS) to lead optimization and candidate selection (CSM-lig, mCSM-lig, Arpeggio, pkCSM). Incorporating these tools along the drug development process can help ensure that candidate leads are chemically and biologically feasible to become successful and tractable drugs.

Towards Unraveling the Histone Code by Fragment Blind Docking

Histones serve as protein spools for winding the DNA in the nucleosome. High variability of their post-translational modifications result in a unique code system often responsible for the pathomechanisms of epigenetics-based diseases. Decoding is performed by reader proteins via complex formation with the N-terminal peptide tails of histones. Determination of structures of histone-reader complexes would be a key to unravel the histone code and the design of new drugs. However, the large number of possible histone complex variations imposes a true challenge for experimental structure determination techniques. Calculation of such complexes is difficult due to considerable size and flexibility of peptides and the shallow binding surfaces of the readers. Moreover, location of the binding sites is often unknown, which requires a blind docking search over the entire surface of the target protein. To accelerate the work in this field, a new approach is presented for prediction of the structure of histone H3 peptide tails docked to their targets. Using a fragmenting protocol and a systematic blind docking method, a collection of well-positioned fragments of the H3 peptide is produced. After linking the fragments, reconstitution of anchoring regions of the target-bound H3 peptide conformations was possible. As a first attempt of combination of blind and fragment docking approaches, our new method is named fragment blind docking (FBD).