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.This work has been funded by grants BIO2016-79930-R and SEV-2015-0493 from the Spanish Ministry of Economy, Industry and Competitiveness, and grant EFA086/15 from EU Interreg V POCTEFA. M Rosell is supported by an FPI fellowship from the Severo Ochoa program. The authors are grateful for the support of the the Joint BSC-CRG-IRB Programme in Computational Biology.Peer ReviewedPostprint (author's final draft

Surface probing by fragment-based screening and computational methods identifies ligandable pockets on the von Hippel-Lindau (VHL) E3 ubiquitin ligase

Beyond the targeting of E3 ubiquitin ligases to inhibit protein homeostasis, E3 ligase binders can be repurposed as targeted protein degraders (PROTACs or molecular glues). We sought to identify new binders of the VHL E3 ligase by biophysical fragment-based screening followed by X-ray crystallographic soaking. We identified fragments binding at the ElonginC:Cullin2 interface and a new cryptic pocket in VHL, along with other potential ligandable sites predicted computationally and found to bind solvent molecules in crystal structures. The elucidated interactions provide starting points for future ligand development

Pyrone-based inhibitors of metalloproteinase types 2 and 3 may work as conformation-selective inhibitors.

Matrix metalloproteinases are zinc-containing enzymes capable of degrading all components of the extracellular matrix. Owing to their role in human disease, matrix metalloproteinase have been the subject of extensive study. A bioinorganic approach was recently used to identify novel inhibitors based on a maltol zinc-binding group, but accompanying molecular-docking studies failed to explain why one of these inhibitors, AM-6, had approximately 2500-fold selectivity for MMP-3 over MMP-2. A number of studies have suggested that the matrix-metalloproteinase active site is highly flexible, leading some to speculate that differences in active-site flexibility may explain inhibitor selectivity. To extend the bioinorganic approach in a way that accounts for MMP-2 and MMP-3 dynamics, we here investigate the predicted binding modes and energies of AM-6 docked into multiple structures extracted from matrix-metalloproteinase molecular dynamics simulations. Our findings suggest that accounting for protein dynamics is essential for the accurate prediction of binding affinity and selectivity. Additionally, AM-6 and other similar inhibitors likely select for and stabilize only a subpopulation of all matrix-metalloproteinase conformations sampled by the apo protein. Consequently, when attempting to predict ligand affinity and selectivity using an ensemble of protein structures, it may be wise to disregard protein conformations that cannot accommodate the ligand

Genes Suggest Ancestral Colour Polymorphisms Are Shared across Morphologically Cryptic Species in Arctic Bumblebees

email Suzanne orcd idCopyright: © 2015 Williams et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Unveiling a Novel Transient Druggable Pocket in BACE-1 through Molecular Simulations: Conformational Analysis and Binding Mode of Multisite Inhibitors

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer's disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challenging. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments - huprine and rhein - that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8-14, 154-169 and 307-318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable 'floppy' pocket would enable the binding of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligomethylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications introduced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors

Insights from quantitative and mathematical modelling on the proposed 2030 goal for gambiense human African trypanosomiasis (gHAT)

Gambiense human African trypanosomiasis (gHAT) is a parasitic, vector-borne neglected tropical disease that has historically affected populations across West and Central Africa and can result in death if untreated. Following from the success of recent intervention programmes against gHAT, the World Health Organization (WHO) has defined a 2030 goal of global elimination of transmission (EOT). The key proposed indicator to measure achievement of the goal is to have zero reported cases. Results of previous mathematical modelling and quantitative analyses are brought together to explore both the implications of the proposed indicator and the feasibility of achieving the WHO goal. Whilst the indicator of zero case reporting is clear and measurable, it is an imperfect proxy for EOT and could arise either before or after EOT is achieved. Lagging reporting of infection and imperfect diagnostic specificity could result in case reporting after EOT, whereas the converse could be true due to underreporting, lack of coverage, and cryptic human and animal reservoirs. At the village-scale, the WHO recommendation of continuing active screening until there are three years of zero cases yields a high probability of local EOT, but extrapolating this result to larger spatial scales is complex. Predictive modelling of gHAT has consistently found that EOT by 2030 is unlikely across key endemic regions if current medical-only strategies are not bolstered by improved coverage, reduced time to detection and/or complementary vector control. Unfortunately, projected costs for strategies expected to meet EOT are high in the short term and strategies that are cost-effective in reducing burden are unlikely to result in EOT by 2030. Future modelling work should aim to provide predictions while taking into account uncertainties in stochastic dynamics and infection reservoirs, as well as assessment of multiple spatial scales, reactive strategies, and measurable proxies of EOT

Mafic alkaline metasomatism in the lithosphere underneath East Serbia: evidence from the study of xenoliths and the host alkali basalts

Effects of mafic alkaline metasomatism have been investigated by a combined study of the East Serbian mantle xenoliths and their host alkaline rocks. Fertile xenoliths and tiny mineral assemblages found in depleted xenoliths have been investigated. Fertile lithologies are represented by clinopyroxene (cpx)-rich lherzolite and spinel (sp)-rich olivine websterite containing Ti–Al-rich Cr-augite, Fe-rich olivine, Fe–Al-rich orthopyroxene and Al-rich spinel. Depleted xenoliths, which are the predominant lithology in the suite of East Serbian xenoliths, are harzburgite, cpx-poor lherzolite and rare Mg-rich dunite. They contain small-scale assemblages occurring as pocket-like, symplectitic or irregular, deformation-assisted accumulations of metasomatic phases, generally composed of Ti–Al- and incompatible element-rich Cr-diopside, Cr–Fe–Ti-rich spinel, altered glass, olivine, apatite, ilmenite, carbonate, feldspar, and a high-TiO2 (c. 11 wt%) phlogopite. The fertile xenoliths are too rich in Al, Ca and Fe to simply represent undepleted mantle. By contrast, their composition can be reproduced by the addition of 5–20 wt% of a basanitic melt to refractory mantle. However, textural relationships found in tiny mineral assemblages inside depleted xenoliths imply the following reaction: opx+sp1 (primary mantle Cr-spinel) ±phlogopite+Si-poor alkaline melt=Ti–Al-cpx+sp2 (metasomatic Ti-rich spinel)±ol±other minor phases. Inversion modelling, performed on the least contaminated and most isotopically uniform host basanites (87Sr/86Sr=c. 0.7031; 143Nd/144Nd=c. 0.5129), implies a source that was enriched in highly and moderately incompatible elements (c. 35–40× chondrite for U–Th–Nb–Ta, 2× chondrite for heavy rare earth elements (HREE), made up of clinopyroxene, carbonate (c. 5%), and traces of ilmenite (c. 1%) and apatite (c. 0.05%). A schematic model involves: first, percolation of CO2- and H2O-rich fluids and precipitation of metasomatic hydrous minerals; and, second, the subsequent breakdown of these hydrous minerals due to the further uplift of hot asthenospheric mantle. This model links intraplate alkaline magmatism to lithospheric mantle sources enriched by sublithospheric melts at some time in the past

Predicting locations of cryptic pockets from single protein structures using the PocketMiner graph neural network

Cryptic pockets expand the scope of drug discovery by enabling targeting of proteins currently considered undruggable because they lack pockets in their ground state structures. However, identifying cryptic pockets is labor-intensive and slow. The ability to accurately and rapidly predict if and where cryptic pockets are likely to form from a structure would greatly accelerate the search for druggable pockets. Here, we present PocketMiner, a graph neural network trained to predict where pockets are likely to open in molecular dynamics simulations. Applying PocketMiner to single structures from a newly curated dataset of 39 experimentally confirmed cryptic pockets demonstrates that it accurately identifies cryptic pockets (ROC-AUC: 0.87) \u3e1,000-fold faster than existing methods. We apply PocketMiner across the human proteome and show that predicted pockets open in simulations, suggesting that over half of proteins thought to lack pockets based on available structures likely contain cryptic pockets, vastly expanding the potentially druggable proteome

Virtual Screening of Plant Volatile Compounds Reveals a High Affinity of Hylamorpha elegans (Coleoptera: Scarabaeidae) Odorant-Binding Proteins for Sesquiterpenes From Its Native Host

Indexación: Web of ScienceHylamorpha elegans (Burmeister) is a native Chilean scarab beetle considered to be a relevant agricultural pest to pasture and cereal and small fruit crops. Because of their cryptic habits, control with conventional methods is difficult; therefore, alternative and environmentally friendly control strategies are highly desirable. The study of proteins that participate in the recognition of odorants, such as odorant-binding proteins (OBPs), offers interesting opportunities to identify new compounds with the potential to modify pest behavior and computational screening of compounds, which is commonly used in drug discovery, may help to accelerate the discovery of new semiochemicals. Here, we report the discovery of four OBPs in H. elegans as well as six new volatiles released by its native host Nothofagus obliqua (Mirbel). Molecular docking performed between OBPs and new and previously reported volatiles from N. obliqua revealed the best binding energy values for sesquiterpenic compounds. Despite remarkable divergence at the amino acid level, three of the four OBPs evaluated exhibited the best interaction energy for the same ligands. Molecular dynamics investigation reinforced the importance of sesquiterpenes, showing that hydrophobic residues of the OBPs interacted most frequently with the tested ligands, and binding free energy calculations demonstrated van der Waals and hydrophobic interactions to be the most important. Altogether, the results suggest that sesquiterpenes are interesting candidates for in vitro and in vivo assays to assess their potential application in pest management strategies.http://jinsectscience.oxfordjournals.org/content/16/1/3

Inclusion of Enclosed Hydration Effects in the Binding Free Energy Estimation of Dopamine D3 Receptor Complexes

Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.Comment: This is the first report of using enclosed hydration in estimating binding free energies of protein-ligand complexes using implicit solvatio