Dataset Information

MixMD Probeview: Robust Binding Site Prediction from Cosolvent Simulations.

ABSTRACT: Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique for identifying binding hotspots and specific favorable interactions on a protein's surface. MixMD studies have the ability to identify these biologically relevant sites by examining the occupancy of the cosolvent over the course of the simulation. However, previous MixMD analysis required a great deal of manual inspection to identify relevant sites. To address this limitation, we have developed MixMD Probeview as a plugin for the freely available, open-source version of the molecular visualization program PyMOL. MixMD Probeview incorporates two analysis procedures: (1) identifying and ranking whole binding sites and (2) identifying and ranking local maxima for each probe type. These functionalities were validated using four common benchmark proteins, including two with both active and allosteric sites. In addition, three different cosolvent procedures were compared to examine the impact of including more than one cosolvent in the simulations. For all systems tested, MixMD Probeview successfully identified known active and allosteric sites based on the total occupancy of neutral probe molecules. As an easy-to-use PyMOL plugin, we expect that MixMD Probeview will facilitate identification and analysis of binding sites from cosolvent simulations performed on a wide range of systems.

SUBMITTER: Graham SE

PROVIDER: S-EPMC6186161 | biostudies-literature | 2018 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

MixMD Probeview: Robust Binding Site Prediction from Cosolvent Simulations.

Graham Sarah E SE Leja Noah N Carlson Heather A HA

Journal of chemical information and modeling 20180626 7

Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique for identifying binding hotspots and specific favorable interactions on a protein's surface. MixMD studies have the ability to identify these biologically relevant sites by examining the occupancy of the cosolvent over the course of the simulation. However, previous MixMD analysis required a great deal of manual inspection to identify relevant sites. To address this limitation, we have developed MixMD Probeview as a plu ...[more]

PMID: 29905479

Similar Datasets

Project description:Mixed-solvent molecular dynamics (MixMD) is a hotspot-mapping technique that relies on molecular dynamics simulations of proteins in binary solvent mixtures. Previous work on MixMD has established the technique's effectiveness in capturing binding sites of small organic compounds. In this work, we show that MixMD can identify both competitive and allosteric sites on proteins. The MixMD approach embraces full protein flexibility and allows competition between solvent probes and water. Sites preferentially mapped by probe molecules are more likely to be binding hotspots. There are two important requirements for the identification of ligand-binding hotspots: (1) hotspots must be mapped at very high signal-to-noise ratio and (2) the hotspots must be mapped by multiple probe types. We have developed our mapping protocol around acetonitrile, isopropanol, and pyrimidine as probe solvents because they allowed us to capture hydrophilic, hydrophobic, hydrogen-bonding, and aromatic interactions. Charged probes were needed for mapping one target, and we introduce them in this work. In order to demonstrate the robust nature and wide applicability of the technique, a combined total of 5 μs of MixMD was applied across several protein targets known to exhibit allosteric modulation. Most notably, all the protein crystal structures used to initiate our simulations had no allosteric ligands bound, so there was no preorganization of the sites to predispose the simulations to find the allosteric hotspots. The protein test cases were ABL Kinase, Androgen Receptor, CHK1 Kinase, Glucokinase, PDK1 Kinase, Farnesyl Pyrophosphate Synthase, and Protein-Tyrosine Phosphatase 1B. The success of the technique is demonstrated by the fact that the top-four sites solely map the competitive and allosteric sites. Lower-ranked sites consistently map other biologically relevant sites, multimerization interfaces, or crystal-packing interfaces. Lastly, we highlight the importance of including protein flexibility by demonstrating that MixMD can map allosteric sites that are not detected in half the systems using FTMap applied to the same crystal structures.

Project description:In drug design, chemical groups are sequentially added to improve a weak-binding fragment into a tight-binding lead molecule. Often, the direction to make these additions is unclear, and there are numerous chemical modifications to choose. Lead development can be guided by crystal structures of the fragment-bound protein, but this alone is unable to capture structural changes like closing or opening of the binding site and any side-chain movements. Accounting for adaptation of the site requires a dynamic approach. Here, we use molecular dynamics calculations of small organic solvents with protein-fragment pairs to reveal the nearest "hot spots". These close hot spots show the direction to make appropriate additions and suggest types of chemical modifications that could improve binding affinity. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that is well established for finding binding "hot spots" in active sites and allosteric sites of proteins. We simulated 20 fragment-bound and apo forms of key pharmaceutical targets to map out hot spots for potential lead space. Furthermore, we analyzed whether the presence of a fragment facilitates the probes' binding in the lead space, a type of binding cooperativity. To the best of our knowledge, this is the first use of cosolvent MD conducted with bound inhibitors in the simulation. Our work provides a general framework to extract molecular features of binding sites to choose chemical groups for growing lead molecules. Of the 20 systems, 17 systems were well mapped by MixMD. For the three not-mapped systems, two had lead growth out into solution away from the protein, and the third had very small modifications which indicated no nearby hot spots. Therefore, our lack of mapping in three systems was appropriate given the experimental data (true-negative cases). The simulations are run for very short time scales, making this method tractable for use in the pharmaceutical industry.

Project description:Protein dynamics play an important role in small molecule binding and can pose a significant challenge in the identification of potential binding sites. Cryptic binding sites have been defined as sites which require significant rearrangement of the protein structure to become physically accessible to a ligand. Mixed-solvent MD (MixMD) is a computational protocol which maps the surface of the protein using molecular dynamics (MD) of the unbound protein solvated in a 5% box of probe molecules with explicit water. This method has successfully identified known active and allosteric sites which did not require reorganization. In this study, we apply the MixMD protocol to identify known cryptic sites of 12 proteins characterized by a wide range of conformational changes. Of these 12 proteins, three require reorganization of side chains, five require loop movements, and four require movement of more significant structures such as whole helices. In five cases, we find that standard MixMD simulations are able to map the cryptic binding sites with at least one probe type. In two cases (guanylate kinase and TIE-2), accelerated MD, which increases sampling of torsional angles, was necessary to achieve mapping of portions of the cryptic binding site missed by standard MixMD. For more complex systems where movement of a helix or domain is necessary, MixMD was unable to map the binding site even with accelerated dynamics, possibly due to the limited timescale (100 ns for individual simulations). In general, similar conformational dynamics are observed in water-only simulations and those with probe molecules. This could imply that the probes are not driving opening events but rather take advantage of mapping sites that spontaneously open as part of their inherent conformational behavior. Finally, we show that docking to an ensemble of conformations from the standard MixMD simulations performs better than docking the apo crystal structure in nine cases and even better than half of the bound crystal structures. Poorer performance was seen in docking to ensembles of conformations from the accelerated MixMD simulations.

Dataset Information

MixMD Probeview: Robust Binding Site Prediction from Cosolvent Simulations.

Publications

MixMD Probeview: Robust Binding Site Prediction from Cosolvent Simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets