Answers.Body: <blockquote>
  <p>NB. This is a long comment as opposed to an answer which requires a colossal working chunk of code.</p>
  </blockquote>
  <h2>Control</h2>
  <p>Science is all about controls failing. A great control for docking is the barnase-barnstar complex (PDB:1BRS): both small, but form a strong interactions. If this is done, the RMSD against the original structure is very meaningful and useful.</p>
  <h2>RMSF</h2>
  <p>A RMSD is a metric of how much two structures differ. In an ensemble, the metric of how much a structure and the ensemble average is called a RMSF (root-mean square fluctuation). Here <em>the RMSF is sought</em> as the RMSD from first/top pose or worse the manually thrown together pre-docking pose does not mean much.</p>
  <h2>Clusters</h2>
  <p>The decoys (PBDs trails) need to be clustered into different ensembles first, before a RMSF can be calculated.</p>
  <h2>∆∆G and water</h2>
  <p>For the calculation of difference in Gibbs free energy between the docked minus the docked poses, the structures need to be energy minimised with the forcefield scorefunction that is to be used for scoring.
  Water is very important for protein-protein interactions, but adds a lot of complication. Say Gromacs is used, <code>gmx energy</code> will return the single state (conformer) Gibbs free energy, but it requires a water box. For the difference, the water box of undocked needs to be the same size and preferably with the same amount of molecules. A much simpler solution is using an implicit solvent model, e.g. Rosetta. As a bonus, the energy minimisation (<code>relax</code>) is relative fast. And specific explicit waters can be added. However, calculating ∆∆G for every decoy will take hundreds of CPU-hours and is a bit nonsensical.</p>
  <p>A better approach is to look at residue conservation (<a href="https://consurf.tau.ac.il/" rel="nofollow noreferrer">consurf</a>), cross-linking mass spec data if available, and if human post-translational modifications (<a href="https://www.phosphosite.org/homeAction.action" rel="nofollow noreferrer">phosphosite</a>) and healthy variants (<a href="https://gnomad.broadinstitute.org/" rel="nofollow noreferrer">gnomAD</a>). But ultimately wet lab validation by ITC of different mutants/labelling based upon the docking results <em>etc.</em> is required.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2021-07-07T15:34:13.817Z
Answers.LastActivityDate: 2021-07-07T15:34:13.817Z
Answers.UserDisplayName: Matteo Ferla
Body: <p>I have the following PDB file that looks like this.
  It is downloadable <a href="https://bit.ly/3xAlF5X" rel="nofollow noreferrer"><strong>here</strong></a>. This PDB already contains ligand and receptor together.</p>
  <p>The red molecule is the ligand and the yellow is the receptor.
  They're already docked with <a href="https://cluspro.bu.edu" rel="nofollow noreferrer"><strong>ClusPro</strong></a>.
  My question is how to calculate the RMSD and binding energy from it?
  Any existing tools (with or without command line) that can help?</p>
  <p>I don't have any data to try it out at hand but have a look at the documentation I linked. There you can find ways of constructing an Ensemble which can then be superposed and RMSDs calculated, even pairwise.</p>
  <p>Reference is the red one (receptor). I have hundreds of other complexed ligands. I want to calculate how good are those complexes. What's the reasonable measure to do that?</p>
  <p>Example code using Prody? please :-)</p>
  <p><a href="https://i.sstatic.net/YuHMK.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/YuHMK.jpg" alt="enter image description here" /></a></p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2023-04-25T20:59:08.053Z
Comments.Text: Could you figure it out how to calculate the binding energy? which software did you use? Thanks,
Comments.UserDisplayName: Giselle SG
Tags: <proteins><protein-structure><pdb><docking>
Title: How to calculate RMSD and binding energy from the already docked ligand-receptor
ThreadText: # How to calculate RMSD and binding energy from the already docked ligand-receptor
  
  (Asked by: username_0 on 2021-06-21)
  
  I have the following PDB file that looks like this.
  It is downloadable [**here**](https://bit.ly/3xAlF5X). This PDB already contains ligand and receptor together.
  
  The red molecule is the ligand and the yellow is the receptor.
  They're already docked with [**ClusPro**](https://cluspro.bu.edu).
  My question is how to calculate the RMSD and binding energy from it?
  Any existing tools (with or without command line) that can help?
  
  I don't have any data to try it out at hand but have a look at the documentation I linked. There you can find ways of constructing an Ensemble which can then be superposed and RMSDs calculated, even pairwise.
  
  Reference is the red one (receptor). I have hundreds of other complexed ligands. I want to calculate how good are those complexes. What's the reasonable measure to do that?
  
  Example code using Prody? please :-)
  
  [![enter image description here](https://i.sstatic.net/YuHMK.jpg)](https://i.sstatic.net/YuHMK.jpg)
  
  **Tags:** proteins, protein-structure, pdb, docking
  
  > **Comments:**
  > **username_1**: Could you figure it out how to calculate the binding energy? which software did you use? Thanks,  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_2 on 2021-07-07. Score: 1
  
  > NB. This is a long comment as opposed to an answer which requires a colossal working chunk of code.
  
  ## Control
  
  Science is all about controls failing. A great control for docking is the barnase-barnstar complex (PDB:1BRS): both small, but form a strong interactions. If this is done, the RMSD against the original structure is very meaningful and useful.
  
  ## RMSF
  
  A RMSD is a metric of how much two structures differ. In an ensemble, the metric of how much a structure and the ensemble average is called a RMSF (root-mean square fluctuation). Here *the RMSF is sought* as the RMSD from first/top pose or worse the manually thrown together pre-docking pose does not mean much.
  
  ## Clusters
  
  The decoys (PBDs trails) need to be clustered into different ensembles first, before a RMSF can be calculated.
  
  ## ∆∆G and water
  
  For the calculation of difference in Gibbs free energy between the docked minus the docked poses, the structures need to be energy minimised with the forcefield scorefunction that is to be used for scoring.
  Water is very important for protein-protein interactions, but adds a lot of complication. Say Gromacs is used, `gmx energy` will return the single state (conformer) Gibbs free energy, but it requires a water box. For the difference, the water box of undocked needs to be the same size and preferably with the same amount of molecules. A much simpler solution is using an implicit solvent model, e.g. Rosetta. As a bonus, the energy minimisation (`relax`) is relative fast. And specific explicit waters can be added. However, calculating ∆∆G for every decoy will take hundreds of CPU-hours and is a bit nonsensical.
  
  A better approach is to look at residue conservation ([consurf](https://consurf.tau.ac.il/)), cross-linking mass spec data if available, and if human post-translational modifications ([phosphosite](https://www.phosphosite.org/homeAction.action)) and healthy variants ([gnomAD](https://gnomad.broadinstitute.org/)). But ultimately wet lab validation by ITC of different mutants/labelling based upon the docking results *etc.* is required.
  
  
  ---
  
  