Solvent Mapping Research Articles

The use of fake ligands from computational solvent mapping in ligand and structure-based virtual screening.

Virtual screening selects compounds that resemble a known modulator or compounds that fit into the binding site of a target protein. Computational solvent mapping defines important chemical features for binding to a target protein. Results/methodology: We have tested the ability to use solvent mapping for generating a 'fake' ligand that is a negative image of the binding site. We used this fake ligand as a query for the program ROCS and to define the search space of the docking programs FRED and HYBRID. The fake ligands perform comparably to or better than the ligands from crystal structures across a set of ten targets. Thus, the approach is suitable for guiding virtual screening and hit-to-lead optimization.

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket.

Ligand-regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug-drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely "expanded" and "contracted", in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single "ligand-adaptable" conformational state. Ensemble-based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named "α8 pocket", whose opening-closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure-based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a "tightly-contracted" conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.

The Multifaceted Roles of Molecular Dynamics Simulations in Drug Discovery.

Discovery of new therapeutics is a very challenging, expensive and time-consuming process. With the number of approved drugs declining steadily, combined with increasing costs, a rational approach is needed to facilitate, expedite and streamline the drug discovery process. In silico methods are playing key roles in the discovery of a growing number of marketed drugs. The use of computational approaches, particularly molecular dynamics, in drug design is rapidly gaining momentum and acceptance as an essential part of the toolkit for modern drug discovery. From analysing atomistic details for explaining experimentally observed phenomena, to designing drugs with increased efficacy and specificity, the insight that such simulations can provide is generating new ideas and applications that have previously been unexplored. Here we discuss physics-based simulation methodologies and applications in drug design: from locating pockets to designing novel lead compounds, from small molecules to peptides. With developments in hardware, software and theory, the improved predictive abilities of in silico efforts are becoming an essential part of efficient, economic and accurate drug development strategies.

Discovery of the glycogen phosphorylase-modulating activity of a resveratrol glucoside by using a virtual screening protocol optimized for solvation effects.

The identification of natural products that can modulate blood glucose levels is of great interest as it can possibly facilitate the utilization of mild interventions such as herbal medicine or functional foods in the treatment of chronic diseases like diabetes. One of the established drug targets for antihyperglycemic therapy is glycogen phosphorylase. To evaluate the glycogen phosphorylase inhibitory properties of an in-house compound collection consisting to a large extent of natural products, a stepwise virtual and experimental screening protocol was devised and implemented. The fact that the active site of glycogen phosphorylase is highly hydrated emphasized that a methodological aspect needed to be efficiently addressed prior to an in silico evaluation of the compound collection. The effect of water molecules on docking calculations was regarded as a key parameter in terms of virtual screening protocol optimization. Statistical analysis of 125 structures of glycogen phosphorylase and solvent mapping focusing on the active site hydration motif in combination with a retrospective screening revealed the importance of a set of 29 crystallographic water molecules for achieving high enrichment as to the discrimination between active compounds and inactive decoys. The scaling of Van der Waals radii of system atoms had an additional effect on screening performance. Having optimized the in silico protocol, a prospective evaluation of the in-house compound collection derived a set of 18 top-ranked natural products that were subsequently evaluated in vitro for their activity as glycogen phosphorylase inhibitors. Two phenolic glucosides with glycogen phosphorylase-modulating activity were identified, whereas the most potent compound affording mid-micromolar inhibition was a glucosidic derivative of resveratrol, a stilbene well-known for its wide range of biological activities. Results show the possible phytotherapeutic and nutraceutical potential of products common in the Mediterranean countries, such as red wine and Vitis products in general or green raw salads and herbal preparations, where such compounds are abundant.

Computational solvent mapping in structure-based drug design.

Over the past two decades, solvent mapping has emerged as a useful tool for identifying hot spots within binding sites on proteins for drug-like molecules and suggesting properties of potential binders. While the experimental technique requires solving multiple crystal structures of a protein in different solvents, computational solvent mapping allows for fast analysis of a protein for potential binding sites and their druggability. Recent advances in genomics, systems biology and interactomics provide a multitude of potential targets for drug development and solvent mapping can provide useful information to help prioritize targets for drug discovery projects. Here, we review various approaches to computational solvent mapping, highlight some key advances and provide our opinion on future directions in the field.

Identification of inhibitors against the potential ligandable sites in the active cholera toxin

The active cholera toxin responsible for the massive loss of water and ions in cholera patients via its ADP ribosylation activity is a heterodimer of the A1 subunit of the bacterial holotoxin and the human cytosolic ARF6 (ADP Ribosylation Factor 6). The active toxin is a potential target for the design of inhibitors against cholera. In this study we identified the potential ligandable sites of the active cholera toxin which can serve as binding sites for drug-like molecules. By employing an energy-based approach to identify ligand binding sites, and comparison with the results of computational solvent mapping, we identified two potential ligandable sites in the active toxin which can be targeted during structure-based drug design against cholera. Based on the probe affinities of the identified ligandable regions, docking-based virtual screening was employed to identify probable inhibitors against these sites. Several indole-based alkaloids and phosphates showed strong interactions to the important residues of the ligandable region at the A1 active site. On the other hand, 26 top scoring hits were identified against the ligandable region at the A1 ARF6 interface which showed strong hydrogen bonding interactions, including guanidines, phosphates, Leucopterin and Aristolochic acid VIa. This study has important implications in the application of hybrid structure-based and ligand-based methods against the identified ligandable sites using the identified inhibitors as reference ligands, for drug design against the active cholera toxin.

Conformational diversity of bacterial FabH: Implications for molecular recognition specificity

The molecular basis of variable substrate and inhibitor specificity of the highly conserved bacterial fatty acid synthase enzyme, FabH, across different bacterial species remains poorly understood. In the current work, we explored the conformational diversity of FabH enzymes to understand the determinants of diverse interaction specificity across Gram-positive and Gram-negative bacteria. Atomistic molecular dynamics simulations reveal that FabH from E. coli and E. faecalis exhibit distinct native state conformational ensembles and dynamic behaviors. Despite strikingly similar substrate binding pockets, hot spot assessment using computational solvent mapping identified quite different favorable binding interactions between the two homologs. Our data suggest that FabH utilizes protein dynamics and seemingly minor sequence and structural differences to modulate its molecular recognition and substrate specificity across bacterial species. These insights will potentially facilitate the rational design and development of antibacterial inhibitors against FabH enzymes.

Molecular insights into the interaction between Plasmodium falciparum apical membrane antigen 1 and an invasion-inhibitory peptide.

Apical membrane antigen 1 (AMA1) of the human malaria parasite Plasmodium falciparum has been implicated in invasion of the host erythrocyte. It interacts with malarial rhoptry neck (RON) proteins in the moving junction that forms between the host cell and the invading parasite. Agents that block this interaction inhibit invasion and may serve as promising leads for anti-malarial drug development. The invasion-inhibitory peptide R1 binds to a hydrophobic cleft on AMA1, which is an attractive target site for small molecules that block parasite invasion. In this work, truncation and mutational analyses show that Phe5-Phe9, Phe12 and Arg15 in R1 are the most important residues for high affinity binding to AMA1. These residues interact with two well-defined binding hot spots on AMA1. Computational solvent mapping reveals that one of these hot spots is suitable for small molecule targeting. We also confirm that R1 in solution binds to AMA1 with 1∶1 stoichiometry and adopts a secondary structure consistent with the major form of R1 observed in the crystal structure of the complex. Our results provide a basis for designing high affinity inhibitors of the AMA1-RON2 interaction.

Evidence of conformational selection driving the formation of ligand binding sites in protein-protein interfaces.

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small “probe” molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.

Ligand binding and activation of PPARγ by Firemaster® 550: effects on adipogenesis and osteogenesis in vitro.

Background: The use of alternative flame retardants has increased since the phase out of pentabromodiphenyl ethers (pentaBDEs). One alternative, Firemaster® 550 (FM550), induces obesity in rats. Triphenyl phosphate (TPP), a component of FM550, has a structure similar to that of organotins, which are obesogenic in rodents.Objectives: We tested the hypothesis that components of FM550 are biologically active peroxisome proliferator-activated receptor γ (PPARγ) ligands and estimated indoor exposure to TPP.Methods: FM550 and its components were assessed for ligand binding to and activation of human PPARγ. Solvent mapping was used to model TPP in the PPARγ binding site. Adipocyte and osteoblast differentiation were assessed in bone marrow multipotent mesenchymal stromal cell models. We estimated exposure of children to TPP using a screening-level indoor exposure model and house dust concentrations determined previously.Results: FM550 bound human PPARγ, and binding appeared to be driven primarily by TPP. Solvent mapping revealed that TPP interacted with binding hot spots within the PPARγ ligand binding domain. FM550 and its organophosphate components increased human PPARγ1 transcriptional activity in a Cos7 reporter assay and induced lipid accumulation and perilipin protein expression in BMS2 cells. FM550 and TPP diverted osteogenic differentiation toward adipogenesis in primary mouse bone marrow cultures. Our estimates suggest that dust ingestion is the major route of exposure of children to TPP.Conclusions: Our findings suggest that FM550 components bind and activate PPARγ. In addition, in vitro exposure initiated adipocyte differentiation and antagonized osteogenesis. TPP likely is a major contributor to these biological actions. Given that TPP is ubiquitous in house dust, further studies are warranted to investigate the health effects of FM550.Citation: Pillai HK, Fang M, Beglov D, Kozakov D, Vajda S, Stapleton HM, Webster TF, Schlezinger JJ. 2014. Ligand binding and activation of PPARγ by Firemaster® 550: effects on adipogenesis and osteogenesis in vitro. Environ Health Perspect 122:1225–1232; http://dx.doi.org/10.1289/ehp.1408111

Utilizing a Dynamical Description of IspH to Aid in the Development of Novel Antimicrobial Drugs

The nonmevalonate pathway is responsible for isoprenoid production in microbes, including H. pylori, M. tuberculosis and P. falciparum, but is nonexistent in humans, thus providing a desirable route for antibacterial and antimalarial drug discovery. We coordinate a structural study of IspH, a [4Fe-4S] protein responsible for converting HMBPP to IPP and DMAPP in the ultimate step in the nonmevalonate pathway. By performing accelerated molecular dynamics simulations on both substrate-free and HMBPP-bound [Fe4S4]2+ IspH, we elucidate how substrate binding alters the dynamics of the protein. Using principal component analysis, we note that while substrate-free IspH samples various open and closed conformations, the closed conformation observed experimentally for HMBPP-bound IspH is inaccessible in the absence of HMBPP. In contrast, simulations with HMBPP bound are restricted from accessing the open states sampled by the substrate-free simulations. Further investigation of the substrate-free simulations reveals large fluctuations in the HMBPP binding pocket, as well as allosteric pocket openings – both of which are achieved through the hinge motions of the individual domains in IspH. Coupling these findings with solvent mapping and various structural analyses reveals alternative druggable sites that may be exploited in future drug design efforts.

Molecular Simulations of Aromatase Reveal New Insights Into the Mechanism of Ligand Binding

CYP19A1, also known as aromatase or estrogen synthetase, is the rate-limiting enzyme in the biosynthesis of estrogens from their corresponding androgens. Several clinically used breast cancer therapies target aromatase. In this work, explicitly solvated all-atom molecular dynamics simulations of aromatase with a model of the lipid bilayer and the transmembrane helix are performed. The dynamics of aromatase and the role of titration of an important amino acid residue involved in aromatization of androgens are investigated via two 250-ns long simulations. One simulation treats the protonated form of the catalytic aspartate 309, which appears more consistent with crystallographic data for the active site, while the simulation of the deprotonated form shows some notable conformational shifts. Ensemble-based computational solvent mapping experiments indicate possible novel druggable binding sites that could be utilized by next-generation inhibitors. In addition, the effects of protonation on the ligand positioning and channel dynamics are investigated using geometrical models that estimate the opening width of critical channels. Significant differences in channel dynamics between the protonated and deprotonated trajectories are exhibited, suggesting that the mechanism for substrate and product entry and the aromatization process may be coupled to a "locking" mechanism and channel opening. Our results may be particularly relevant in the design of novel drugs, which may be useful therapeutic treatments of cancers such as those of the breast and prostate.

Synergistic Inhibitor Binding to the Papain‐Like Protease of Human SARS Coronavirus: Mechanistic and Inhibitor Design Implications

We previously developed two potent chemical classes that inhibit the essential papain‐like protease (PLpro) of severe acute respiratory syndrome coronavirus. In this study, we applied a novel approach to identify small fragments that act synergistically with these inhibitors. A fragment library was screened in combination with four previously developed lead inhibitors by fluorescence‐based enzymatic assays. Several fragment compounds synergistically enhanced the inhibitory activity of the lead inhibitors by approximately an order of magnitude. Surface plasmon resonance measurements showed that three fragments bind specifically to the PLpro enzyme. Mode of inhibition, computational solvent mapping, and molecular docking studies suggest that these fragments bind adjacent to the binding site of the lead inhibitors and further stabilize the inhibitor‐bound state. We propose potential next‐generation compounds based on a computational fragment‐merging approach. This approach provides an alternative strategy for lead optimization for cases in which direct co‐crystallization is difficult.

FTFlex: accounting for binding site flexibility to improve fragment-based identification of druggable hot spots

Computational solvent mapping finds binding hot spots, determines their druggability and provides information for drug design. While mapping of a ligand-bound structure yields more accurate results, usually the apo structure serves as the starting point in design. The FTFlex algorithm, implemented as a server, can modify an apo structure to yield mapping results that are similar to those of the respective bound structure. Thus, FTFlex is an extension of our FTMap server, which only considers rigid structures. FTFlex identifies flexible residues within the binding site and determines alternative conformations using a rotamer library. In cases where the mapping results of the apo structure were in poor agreement with those of the bound structure, FTFlex was able to yield a modified apo structure, which lead to improved FTMap results. In cases where the mapping results of the apo and bound structures were in good agreement, no new structure was predicted. FTFlex is freely available as a web-based server at http://ftflex.bu.edu/.

Variable ligand- and receptor-binding hot spots in key strains of influenza neuraminidase

Influenza A continues to be a major public health concern due to its ability to cause epidemic and pandemic disease outbreaks in humans. Computational investigations of structural dynamics of the major influenza glycoproteins, especially the neuraminidase (NA) enzyme, are able to provide key insights beyond what is currently accessible with standard experimental techniques. In particular, all-atom molecular dynamics simulations reveal the varying degrees of flexibility for such enzymes. Here we present an analysis of the relative flexibility of the ligand- and receptor-binding area of three key strains of influenza A: highly pathogenic H5N1, the 2009 pandemic H1N1, and a human N2 strain. Through computational solvent mapping, we investigate the various ligand- and receptor-binding "hot spots" that exist on the surface of NA which interacts with both sialic acid receptors on the host cells and antiviral drugs. This analysis suggests that the variable cavities found in the different strains and their corresponding capacities to bind ligand functional groups may play an important role in the ability of NA to form competent reaction encounter complexes with other species of interest, including antiviral drugs, sialic acid receptors on the host cell surface, and the hemagglutinin protein. Such considerations may be especially useful for the prediction of how such complexes form and with what binding capacity.

Analysis of Binding Site Hot Spots on the Surface of Ras GTPase

We have recently discovered an allosteric switch in Ras, bringing an additional level of complexity to this GTPase whose mutants are involved in nearly 30% of cancers. Upon activation of the allosteric switch, there is a shift in helix 3/loop 7 associated with a disorder to order transition in the active site. Here, we use a combination of multiple solvent crystal structures and computational solvent mapping (FTMap) to determine binding site hot spots in the “off” and “on” allosteric states of the GTP-bound form of H-Ras. Thirteen sites are revealed, expanding possible target sites for ligand binding well beyond the active site. Comparison of FTMaps for the H and K isoforms reveals essentially identical hot spots. Furthermore, using NMR measurements of spin relaxation, we determined that K-Ras exhibits global conformational dynamics very similar to those we previously reported for H-Ras. We thus hypothesize that the global conformational rearrangement serves as a mechanism for allosteric coupling between the effector interface and remote hot spots in all Ras isoforms. At least with respect to the binding sites involving the G domain, H-Ras is an excellent model for K-Ras and probably N-Ras as well. Ras has so far been elusive as a target for drug design. The present work identifies various unexplored hot spots throughout the entire surface of Ras, extending the focus from the disordered active site to well-ordered locations that should be easier to target.

Structural conservation of druggable hot spots in protein–protein interfaces

Despite the growing number of examples of small-molecule inhibitors that disrupt protein-protein interactions (PPIs), the origin of druggability of such targets is poorly understood. To identify druggable sites in protein-protein interfaces we combine computational solvent mapping, which explores the protein surface using a variety of small "probe" molecules, with a conformer generator to account for side-chain flexibility. Applications to unliganded structures of 15 PPI target proteins show that the druggable sites comprise a cluster of binding hot spots, distinguishable from other regions of the protein due to their concave topology combined with a pattern of hydrophobic and polar functionality. This combination of properties confers on the hot spots a tendency to bind organic species possessing some polar groups decorating largely hydrophobic scaffolds. Thus, druggable sites at PPI are not simply sites that are complementary to particular organic functionality, but rather possess a general tendency to bind organic compounds with a variety of structures, including key side chains of the partner protein. Results also highlight the importance of conformational adaptivity at the binding site to allow the hot spots to expand to accommodate a ligand of drug-like dimensions. The critical components of this adaptivity are largely local, involving primarily low energy side-chain motions within 6 Å of a hot spot. The structural and physicochemical signature of druggable sites at PPI interfaces is sufficiently robust to be detectable from the structure of the unliganded protein, even when substantial conformational adaptation is required for optimal ligand binding.

Domain motion and interdomain hot spots in a multidomain enzyme

The aim of this article is to analyze conformational changes by comparing 10 different structures of Pseudomonas aeruginosa phosphomannomutase/phosphoglucomutase (PMM/PGM), a four-domain enzyme in which both substrate binding and catalysis require substantial movement of the C-terminal domain. We focus on changes in interdomain and active site crevices using a method called computational solvent mapping rather than superimposing the structures. The method places molecular probes (i.e., small organic molecules containing various functional groups) around the protein to find hot spots. One of the most important hot spots is in the active site, consistent with the ability of the enzyme to bind both glucose and mannose phosphosugar substrates. The protein has eight additional hot spots at domain-domain interfaces and hinge regions. The locations and nature of six of these hot spots vary between the open, half-open, and closed conformers of the enzyme, in good agreement with the ligand-induced conformational changes. In the closed structures the number of probe clusters at the hinge region significantly depends on the position of the phosphorylated oxygen in the substrate (e.g., glucose 1-phosphate versus glucose 6-phosphate), but the protein remains almost unchanged in terms of the overall RMSD, indicating that computational solvent mapping is a more sensitive approach to detect changes in binding sites and interdomain crevices. Focusing on multidomain proteins we show that the subresolution conformational differences revealed by the mapping are in fact significant, and present a general statistical method of analysis to determine the significance of rigid body domain movements in X-ray structures.

The Structural Basis of Pregnane X Receptor Binding Promiscuity

The steroid and xenobiotic-responsive human pregnane X receptor (PXR) binds a broad range of structurally diverse compounds. The structures of the apo and ligand-bound forms of PXR are very similar, in contrast to most promiscuous proteins that generally adapt their shape to different ligands. We investigated the structural origins of PXR's recognition promiscuity using computational solvent mapping, a technique developed for the identification and characterization of hot spots, i.e., regions of the protein surface that are major contributors to the binding free energy. Results reveal that the smooth and nearly spherical binding site of PXR has a well-defined hot spot structure, with four hot spots located on four different sides of the pocket and a fifth close to its center. Three of these hot spots are already present in the ligand-free protein. The most important hot spot is defined by three structurally and sequentially conserved residues, W299, F288, and Y306. This largely hydrophobic site is not very specific and interacts with all known PXR ligands. Depending on their sizes and shapes, individual PXR ligands extend into two, three, or four more hot spot regions. The large number of potential arrangements within the binding site explains why PXR is able to accommodate a large variety of compounds. All five hot spots include at least one important residue, which is conserved in all mammalian PXRs, suggesting that the hot spot locations have remained largely invariant during mammalian evolution. The same side chains also show a high level of structural conservation across hPXR structures. However, each of the hPXR hot spots also includes residues with moveable side chains, further increasing the size variation in ligands that PXR can bind. Results also suggest a unique signal transduction mechanism between the PXR homodimerization interface and its coactivator binding site.

Binding Hot Spots and Amantadine Orientation in the Influenza A Virus M2 Proton Channel

Structures of truncated versions of the influenza A virus M2 proton channel have been determined recently by x-ray crystallography in the open conformation of the channel, and by NMR in the closed state. The structures differ in the position of the bound inhibitors. The x-ray structure shows a single amantadine molecule in the middle of the channel, whereas in the NMR structure four drug molecules bind at the channel's outer surface. To study this controversy we applied computational solvent mapping, a technique developed for the identification of the most druggable binding hot spots of proteins. The method moves molecular probes—small organic molecules containing various functional groups—around the protein surface, finds favorable positions using empirical free energy functions, clusters the conformations, and ranks the clusters on the basis of the average free energy. The results of the mapping show that in both structures the primary hot spot is an internal cavity overlapping the amantadine binding site seen in the x-ray structure. However, both structures also have weaker hot spots at the exterior locations that bind rimantadine in the NMR structure, although these sites are partially due to the favorable interactions with the interfacial region of the lipid bilayer. As confirmed by docking calculations, the open channel binds amantadine at the more favorable internal site, in good agreement with the x-ray structure. In contrast, the NMR structure is based on a peptide/micelle construct that is able to accommodate the small molecular probes used for the mapping, but has a too narrow pore for the rimantadine to access the internal hot spot, and hence the drug can bind only at the exterior sites.