CrypticScout: Mixed Solvent Molecular Dynamics For Drug Discovery
This protocol simulates the target of interest (GPCR, kinase, ion channel, etc.) in a solution of water and a co-solvent, like benzene or any other fragment-like molecule.
During the simulation, the co-solvent molecules interact with the surface of the protein, revealing binding hotspots. These hotspots have been proven to correlate very well with actual pockets.
Furthermore, binding modes for the probes can be extracted, which can prove very valuable, as the binding mode of bigger, drug-like molecules could mirror that of the probe.
We can run multiple different fragments in parallel to conduct a virtual screening campaign.
- Binding site prediction: Identify the most likely binding sites of your target of interest, including orthosteric, allosteric and cryptic sites.
- Binding mode prediction: Obtain binding mode predictions and pharmacophoric insights from your co-solvent molecule, which can be used to precisely guide a docking campaign.
- Pocket ensemble: For each pocket, we will select the snapshots from the simulation where the pocket is in a “binding-ready” state. This ensemble can be used to model protein flexibility while docking.
- Binding pathway: For each pocket, we will provide a trajectory linking the unbound state to the bound state. CrypticScout job typically takes only one week to finish (~10 replicas of 80 nanoseconds each).
Results of a CrypticScout job for a Ras protein (PDB code:4Q21). At high isovalues, the orange wireframe only highligths the pockets which the probe visited the most (or the longest). One of the identified pockets is the actual binding pocket of the endogenous ligand. The extracted pose for the chosen co-solvent (imidazole) overlaps very well with the actual ligand.
- A PlayMolecule scene with binding hotspots and binding mode predictions.
- The full simulations (.xtc and .pdb files).
- The cloud of hotspots as a .cube file (orange wireframe in the image).
- An extensive report summarizing the structural insights obtained.
For each identified pocket:
- Binding mode prediction can be used to do pharmacophoric docking or scaffold docking. (.pdb + .sdf)
- Ensemble of druggable conformations, to model protein flexibility in docking (.pdb files)
- Binding pathway of the co-solvent molecule from the bulk into the pocket (.xtc file)