NaViA: a program for the visual analysis of complex mass spectra.

MOTIVATION: Native mass spectrometry is now a well-established method for the investigation of protein complexes, specifically their subunit stoichiometry and ligand binding properties. Recent advances allowing the analysis of complex mixtures lead to an increasing diversity and complexity in the spectra obtained. These spectra can be time-consuming to tackle through manual assignment and challenging for automated approaches. RESULTS: Native Mass Spectrometry Visual Analyser is a web-based tool to augment the manual process of peak assignment. In addition to matching masses to the stoichiometry of its component subunits, it allows raw data processing, assignment and annotation and permits mass spectra to be shared with their respective interpretation. AVAILABILITY AND IMPLEMENTATION: NaViA is open-source and can be accessed online under https://navia.ms. The source code and documentation can be accessed at https://github.com/d-que/navia, under the BSD 2-Clause licence. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PlayMolecule CrypticScout: Predicting Protein Cryptic Sites Using Mixed-Solvent Molecular Simulations.

Cryptic pockets are protein cavities that remain hidden in resolved apo structures and generally require the presence of a co-crystallized ligand to become visible. Finding new cryptic pockets is crucial for structure-based drug discovery to identify new ways of modulating protein activity and thus expand the druggable space. We present here a new method and associated web application leveraging mixed-solvent molecular dynamics (MD) simulations using benzene as a hydrophobic probe to detect cryptic pockets. Our all-atom MD-based workflow was systematically tested on 18 different systems and 5 additional kinases and represents the largest validation study of this kind. CrypticScout identifies benzene probe binding hotspots on a protein surface by mapping probe occupancy, residence time, and the benzene occupancy reweighed by the residence time. The method is presented to the scientific community in a web application available via www.playmolecule.org using a distributed computing infrastructure to perform the simulations.

The Aminotriazole Antagonist Cmpd-1 Stabilises a Distinct Inactive State of the Adenosine 2A Receptor.

The widely expressed G-protein coupled receptors (GPCRs) are versatile signal transducer proteins that are attractive drug targets but structurally challenging to study. GPCRs undergo a number of conformational rearrangements when transitioning from the inactive to the active state but have so far been believed to adopt a fairly conserved inactive conformation. Using 19 F NMR spectroscopy and advanced molecular dynamics simulations we describe a novel inactive state of the adenosine 2A receptor which is stabilised by the aminotriazole antagonist Cmpd-1. We demonstrate that the ligand stabilises a unique conformation of helix V and present data on the putative binding mode of the compound involving contacts to the transmembrane bundle as well as the extracellular loop 2.

Towards a Molecular Understanding of the Link between Imatinib Resistance and Kinase Conformational Dynamics.

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics.

Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators.

Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.

The different flexibility of c-Src and c-Abl kinases regulates the accessibility of a druggable inactive conformation.

c-Src and c-Abl are two closely related protein kinases that constitute important anticancer targets. Despite their high sequence identity, they show different sensitivities to the anticancer drug imatinib, which binds specifically to a particular inactive conformation in which the Asp of the conserved DFG motif points outward (DFG-out). We have analyzed the DFG conformational transition of the two kinases using massive molecular dynamics simulations, free energy calculations, and isothermal titration calorimetry. On the basis of the reconstruction of the free energy surfaces for the DFG-in to DFG-out conformational changes of c-Src and c-Abl, we propose that the different flexibility of the two kinases results in a different stability of the DFG-out conformation and might be the main determinant of imatinib selectivity.

Crystal structure of dopamine D1 receptor in complex with G protein and a non-catechol agonist.

Dopamine D1 receptor (D1R) is an important drug target implicated in many psychiatric and neurological disorders. Selective agonism of D1R are sought to be the therapeutic strategy for these disorders. Most selective D1R agonists share a dopamine-like catechol moiety in their molecular structure, and their therapeutic potential is therefore limited by poor pharmacological properties in vivo. Recently, a class of non-catechol D1R selective agonists with a distinct scaffold and pharmacological properties were reported. Here, we report the crystal structure of D1R in complex with stimulatory G protein (Gs) and a non-catechol agonist Compound 1 at 3.8 Å resolution. The structure reveals the ligand bound to D1R in an extended conformation, spanning from the orthosteric site to extracellular loop 2 (ECL2). Structural analysis reveals that the unique features of D1R ligand binding pocket explains the remarkable selectivity of this scaffold for D1R over other aminergic receptors, and sheds light on the mechanism for D1R activation by the non-catechol agonist.

Reconstruction of apo A2A receptor activation pathways reveal ligand-competent intermediates and state-dependent cholesterol hotspots.

G-protein coupled receptors (GPCRs) play a pivotal role in transmitting signals at the cellular level. Structural insights can be exploited to support GPCR structure-based drug discovery endeavours. Despite advances in GPCR crystallography, active state structures are scarce. Molecular dynamics (MD) simulations have been used to explore the conformational landscape of GPCRs. Efforts have been made to retrieve active state conformations starting from inactive structures, however to date this has not been possible without using an energy bias. Here, we reconstruct the activation pathways of the apo adenosine receptor (A2A), starting from an inactive conformation, by applying adaptive sampling MD combined with a goal-oriented scoring function. The reconstructed pathways reconcile well with experiments and help deepen our understanding of A2A regulatory mechanisms. Exploration of the apo conformational landscape of A2A reveals the existence of ligand-competent states, active intermediates and state-dependent cholesterol hotspots of relevance for drug discovery. To the best of our knowledge this is the first time an activation process has been elucidated for a GPCR starting from an inactive structure only, using a non-biased MD approach, opening avenues for the study of ligand binding to elusive yet pharmacologically relevant GPCR states.

Defining the architecture of KPC-2 Carbapenemase: identifying allosteric networks to fight antibiotics resistance.

The rise of multi-drug resistance in bacterial pathogens is one of the grand challenges facing medical science. A major concern is the speed of development of ß-lactamase-mediated resistance in Gram-negative species, thus putting at risk the efficacy of the most recently approved antibiotics and inhibitors, including carbapenems and avibactam, respectively. New strategies to overcome resistance are urgently required, which will ultimately be facilitated by a deeper understanding of the mechanisms that regulate the function of ß-lactamases such as the Klebsiella Pneumoniae carbapenemases (KPCs). Using enhanced sampling computational methods together with site-directed mutagenesis, we report the identification of two "hydrophobic networks" in the KPC-2 enzyme, the integrity of which has been found to be essential for protein stability and corresponding resistance. Present throughout the structure, these networks are responsible for the structural integrity and allosteric signaling. Disruption of the networks leads to a loss of the KPC-2 mediated resistance phenotype, resulting in restored susceptibility to different classes of ß-lactam antibiotics including carbapenems and cephalosporins. The "hydrophobic networks" were found to be highly conserved among class-A ß-lactamases, which implies their suitability for exploitation as a potential target for therapeutic intervention.

Conformational Selection and Induced Fit Mechanisms in the Binding of an Anticancer Drug to the c-Src Kinase.

Understanding the conformational changes associated with the binding of small ligands to their biological targets is a fascinating and meaningful question in chemistry, biology and drug discovery. One of the most studied and important is the so-called "DFG-flip" of tyrosine kinases. The conserved three amino-acid DFG motif undergoes an "in to out" movement resulting in a particular inactive conformation to which "type II" kinase inhibitors, such as the anti-cancer drug Imatinib, bind. Despite many studies, the details of this prototypical conformational change are still debated. Here we combine various NMR experiments and surface plasmon resonance with enhanced sampling molecular dynamics simulations to shed light into the conformational dynamics associated with the binding of Imatinib to the proto-oncogene c-Src. We find that both conformational selection and induced fit play a role in the binding mechanism, reconciling opposing views held in the literature. Moreover, an external binding pose and local unfolding (cracking) of the aG helix are observed.