Docking to large allosteric binding sites on protein surfaces

Re-Exploring the Ability of Common Docking Programs to Correctly Reproduce the Binding Modes of Non-Covalent Inhibitors of SARS-CoV-2 Protease Mpro

In the latest few decades, molecular docking has imposed itself as one of the most used approaches for computational drug discovery. Several docking benchmarks have been published, comparing the performance of different algorithms in respect to a molecular target of interest, usually evaluating their ability in reproducing the experimental data, which, in most cases, comes from X-ray structures. In this study, we elucidated the variation of the performance of three docking algorithms, namely GOLD, Glide, and PLANTS, in replicating the coordinates of the crystallographic ligands of SARS-CoV-2 main protease (Mpro). Through the comparison of the data coming from docking experiments and the values derived from the calculation of the solvent exposure of the crystallographic ligands, we highlighted the importance of this last variable for docking performance. Indeed, we underlined how an increase in the percentage of the ligand surface exposed to the solvent in a crystallographic complex makes it harder for the docking algorithms to reproduce its conformation. We further validated our hypothesis through molecular dynamics simulations, showing that the less stable protein–ligand complexes (in terms of root-mean-square deviation and root-mean-square fluctuation) tend to be derived from the cases in which the solvent exposure of the ligand in the starting system is higher.

Repurposing of FDA-approved drugs against active site and potential allosteric drug-binding sites of COVID-19 main protease

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still has serious negative effects on health, social life, and economics. Recently, vaccines from various companies have been urgently approved to control SARS-CoV-2 infections. However, any specific antiviral drug has not been confirmed so far for regular treatment. An important target is the main protease (Mpro ), which plays a major role in replication of the virus. In this study, Gaussian and residue network models are employed to reveal two distinct potential allosteric sites on Mpro that can be evaluated as drug targets besides the active site. Then, Food and Drug Administration (FDA)-approved drugs are docked to three distinct sites with flexible docking using AutoDock Vina to identify potential drug candidates. Fourteen best molecule hits for the active site of Mpro are determined. Six of these also exhibit high docking scores for the potential allosteric regions. Full-atom molecular dynamics simulations with MM-GBSA method indicate that compounds docked to active and potential allosteric sites form stable interactions with high binding free energy (∆Gbind ) values. ∆Gbind values reach -52.06 kcal/mol for the active site, -51.08 kcal/mol for the potential allosteric site 1, and - 42.93 kcal/mol for the potential allosteric site 2. Energy decomposition calculations per residue elucidate key binding residues stabilizing the ligands that can further serve to design pharmacophores. This systematic and efficient computational analysis successfully determines ivermectine, diosmin, and selinexor currently subjected to clinical trials, and further proposes bromocriptine, elbasvir as Mpro inhibitor candidates to be evaluated against SARS-CoV-2 infections.

Targeting the C-Terminal Domain Small Phosphatase 1

The human C-terminal domain small phosphatase 1 (CTDSP1/SCP1) is a protein phosphatase with a conserved catalytic site of DXDXT/V. CTDSP1’s major activity has been identified as dephosphorylation of the 5th Ser residue of the tandem heptad repeat of the RNA polymerase II C-terminal domain (RNAP II CTD). It is also implicated in various pivotal biological activities, such as acting as a driving factor in repressor element 1 (RE-1)-silencing transcription factor (REST) complex, which silences the neuronal genes in non-neuronal cells, G1/S phase transition, and osteoblast differentiation. Recent findings have denoted that negative regulation of CTDSP1 results in suppression of cancer invasion in neuroglioma cells. Several researchers have focused on the development of regulating materials of CTDSP1, due to the significant roles it has in various biological activities. In this review, we focused on this emerging target and explored the biological significance, challenges, and opportunities in targeting CTDSP1 from a drug designing perspective.

Porphobilinogen synthase: An equilibrium of different assemblies in human health

Porphobilinogen synthase (PBGS) is an essential enzyme that catalyzes an early step in heme biosynthesis. An unexpected human PBGS quaternary structure dynamic drove the definition of morpheeins, which are protein multimers that dissociate, change shape, and re-assemble differently with functional consequences. Each PBGS monomer has two domains that can reposition through a hinge motion. Human PBGS exists in an equilibrium among high activity octamer, low activity hexamer, and low mole-fraction dimer in which the hinge motion occurs. The dimer conformation dictates the multimer architecture. An octamer-specific inter-subunit interaction responds to pH, resulting in a pH-dependence to the octamer-hexamer equilibrium. An inborn error of metabolism, ALAD porphyria, is caused by single amino acid substitutions that stabilize the hexamer relative to octamer. Drugs that stabilize the PBGS hexamer result in a drug side effect that can exacerbate porphyria. PBGS is essential for all organisms that require respiration, photosynthesis, or methanogenesis. Consequently, phylogenetic variation in PBGS multimerization equilibria provides insight into how Nature has harnessed oligomeric variation in the control of protein function. The dynamic multimerization of PBGS revealed the morpheein mechanism for allostery, a structural basis for inborn errors of metabolism, a quaternary structure focus for drug discovery and/or drug side effects, and a pathway toward new antibiotics or herbicides. The fortuitous discovery of PBGS quaternary structure dynamics arose from characterization of a low-activity single amino acid variant that dramatically stabilized the hexamer, whose existence had previously gone unnoticed.

Systematic discovery of molecular probes targeting multiple non-orthosteric and spatially distinct sites in the botulinum neurotoxin subtype A (BoNT/A)

The development of molecular probes targeting proteins has traditionally relied on labeling compounds already known to bind to the protein of interest. These known ligands bind to orthosteric or allosteric sites in their target protein as a way to control their activity. Binding pockets other than known orthosteric or allosteric sites may exist that are large enough to accommodate a ligand without significantly disrupting protein activity. Such sites may provide opportunities to discriminate between subtypes or other closely related proteins, since they are under less evolutionary pressure to be conserved. The Protein Scanning with Virtual Ligand Screening (PSVLS) approach was previously used to identify a novel inhibitor and a fluorescent probe against the catalytic site of the botulinum neurotoxin subtype A (BoNT/A). PSVLS screens compound databases against multiple sites within a target protein, and the results for all the sites probed against BoNT/A, not only the catalytic site, are available online. Here, we analyze the PSVLS data for multiple sites in order to identify molecular probes with affinity for binding pockets other than the catalytic site of BoNT/A. BoNT/A is a large protein with a light (LC) and a heavy (HC) chain that can be assayed separately. We used scintillation proximity assay (SPA) to test experimentally 5 probe candidates predicted computationally to have affinity for different non-orthosteric binding regions within the HC and LC, and one compound predicted not to have affinity for either domain. The binding profiles obtained experimentally confirmed the targeting of multiple and spatially distinct pockets within BoNT/A. Moreover, inhibition assay results indicate that some of these probes do not significantly interfere with the catalytic activity of BoNT/A.

Enzyme immobilization and molecular modeling studies on an organic–inorganic polypyrrole–titanium(iv)phosphate nanocomposite

A polypyrrole–titanium(iv)phosphate nanocomposite was synthesized by using facile chemical oxidative polymerization of pyrrole in the presence of titanium(iv)phosphate for YADH immobilization.

Impact of quaternary structure dynamics on allosteric drug discovery

The morpheein model of allosteric regulation draws attention to proteins that can exist as an equilibrium of functionally distinct assemblies where: one subunit conformation assembles into one multimer; a different subunit conformation assembles into a different multimer; and the various multimers are in a dynamic equilibrium whose position can be modulated by ligands that bind to a multimer-specific ligand binding site. The case study of porphobilinogen synthase (PBGS) illustrates how such an equilibrium holds lessons for disease mechanisms, drug discovery, understanding drug side effects, and identifying proteins wherein drug discovery efforts might focus on quaternary structure dynamics. The morpheein model of allostery has been proposed as applicable for a wide assortment of disease-associated proteins (Selwood, T., Jaffe, E., (2012) Arch. Bioch. Biophys, 519:131-143). Herein we discuss quaternary structure dynamics aspects to drug discovery for the disease-associated putative morpheeins phenylalanine hydroxylase, HIV integrase, pyruvate kinase, and tumor necrosis factor α. Also highlighted is the quaternary structure equilibrium of transthyretin and successful drug discovery efforts focused on controlling its quaternary structure dynamics.

A new model for allosteric regulation of phenylalanine hydroxylase: implications for disease and therapeutics

The structural basis for allosteric regulation of phenylalanine hydroxylase (PAH), whose dysfunction causes phenylketonuria (PKU), is poorly understood. A new morpheein model for PAH allostery is proposed to consist of a dissociative equilibrium between two architecturally different tetramers whose interconversion requires a ∼90° rotation between the PAH catalytic and regulatory domains, the latter of which contains an ACT domain. This unprecedented model is supported by in vitro data on purified full length rat and human PAH. The conformational change is both predicted to and shown to render the tetramers chromatographically separable using ion exchange methods. One novel aspect of the activated tetramer model is an allosteric phenylalanine binding site at the intersubunit interface of ACT domains. Amino acid ligand-stabilized ACT domain dimerization follows the multimerization and ligand binding behavior of ACT domains present in other proteins in the PDB. Spectroscopic, chromatographic, and electrophoretic methods demonstrate a PAH equilibrium consisting of two architecturally distinct tetramers as well as dimers. We postulate that PKU-associated mutations may shift the PAH quaternary structure equilibrium in favor of the low activity assemblies. Pharmacological chaperones that stabilize the ACT:ACT interface can potentially provide PKU patients with a novel small molecule therapeutic.

The morpheein model of allostery: evaluating proteins as potential morpheeins

An equilibrium mixture of alternate quaternary structure assemblies can form a basis for allostery. The morpheein model of allostery is a concerted dissociative model that describes an equilibrium of alternate quaternary structure assemblies whose architectures are dictated by alternate conformations in the dissociated state. Kinetic and biophysical anomalies that suggest that the morpheein model of allostery applies for a given protein of interest are briefly described. Two methods are presented for evaluating proteins as potential morpheeins. One is a subunit interchange method that uses chromatography, dialysis, and mass spectroscopy to monitor changes in multimer composition. The other is a two-dimensional native gel electrophoresis method to monitor ligand-induced changes in an equilibrium of alternate multimeric assemblies.

Allostery and the dynamic oligomerization of porphobilinogen synthase

The structural basis for allosteric regulation of porphobilinogen synthase (PBGS) is modulation of a quaternary structure equilibrium between octamer and hexamer (via dimers), which is represented schematically as 8mer ⇔ 2mer ⇔ 2mer∗⇔ 6mer∗. The "∗" represents a reorientation between two domains of each subunit that occurs in the dissociated state because it is sterically forbidden in the larger multimers. Allosteric effectors of PBGS are both intrinsic and extrinsic and are phylogenetically variable. In some species this equilibrium is modulated intrinsically by magnesium which binds at a site specific to the 8mer. In other species this equilibrium is modulated intrinsically by pH with the guanidinium group of an arginine being spatially equivalent to the allosteric magnesium ion. In humans, disease associated variants all shift the equilibrium toward the 6mer∗ relative to wild type. The 6mer∗ has a surface cavity that is not present in the 8mer and is proposed as a small molecule allosteric binding site. In silico and in vitro approaches have revealed species-specific allosteric PBGS inhibitors that stabilize the 6mer∗. Some of these inhibitors are drugs in clinical use leading to the hypothesis that extrinsic allosteric inhibition of human PBGS could be a mechanism for drug side effects.