Answers.Body: <p>Yes you can modify the reference PDB file and look for the changes and for this purpose you need <a href="https://www.rcsb.org/pdb/static.do?p=software/software_links/molecular_graphics.html" rel="noreferrer">visualizers.</a> One among these is <a href="https://www.cgl.ucsf.edu/chimera/" rel="noreferrer">Chimera</a>. You can easily carry out the energy minimization steps using Chimera, first mutate your PDB file manually and then provide it as an input for Chimera and the steps for doing energy minimization are fairly simple as given in Chimera docs <a href="https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/minimize/minimize.html" rel="noreferrer">here.</a></p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-06-12T11:04:18.347Z
Answers.LastActivityDate: 2018-06-12T11:04:18.347Z
Answers.UserDisplayName: Abdul Ahad
Answers.Body: <p>The previous answers suggest Chimera and FoldX, which require a license for commercial use.</p>
  <p>Commercial users can use <a href="https://github.com/openmm/pdbfixer/" rel="nofollow noreferrer">pdbfixer</a> or a similar tool to perform mutations (using e.g. <a href="https://github.com/openmm/pdbfixer/blob/master/pdbfixer/pdbfixer.py#L729" rel="nofollow noreferrer">applyMutation</a> in the API).</p>
  <p>You can then use <a href="https://openmm.org/" rel="nofollow noreferrer">OpenMM</a> to <a href="http://docs.openmm.org/latest/userguide/application/02_running_sims.html#a-first-example" rel="nofollow noreferrer">refine the mutated structure</a>.</p>
  <p>Granted, this requires some comfort writing Python scripts to use the APIs.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2023-11-13T20:18:23.347Z
Answers.LastActivityDate: 2023-11-13T20:18:23.347Z
Answers.UserDisplayName: Maximilian Press
Answers.Body: <p>I would recommend you try the FoldX mutation engine.  It is very easy to use and state of the art (it is the main competitor to expensive commercial software).  It can do both of what you want - the mutation and energy optimisation.  It may not work as well as molecular dynamics software, but for your simple task, it is a reasonable solution.</p>
  
  <p>A simple command line can do it:</p>
  
  <p>First, repair your structure (since, usually PDBs are not energetically optimised, so you could get better energy anyway without doing this):</p>
  
  <pre><code>FoldX --command=repairpdb --pdb=1AJY.pdb
  </code></pre>
  
  <p>Then, either use AlaScan or PositionScan to mutate.  AlaScan only mutates to Alanine.  PositionScan can mutate to any amino acid. More advanced mutations can be done with the BuildModel command.</p>
  
  <pre><code>FoldX --command=PositionScan --pdb=1AJY_Repair.pdb --output-file=1AJY_Repair_MA30d --out-pdb=false --positions=MA30d
  </code></pre>
  
  <p>Note: 'd' means try mutate the residue to 24 different amino acids; or, you can specify any of those 24, including the standard 20.</p>
  
  <p>Which produces two files:</p>
  
  <ol>
  <li>PS_1AJY_Repair_MA30d_scanning_output.txt</li>
  <li>energies_30_1AJY_Repair</li>
  </ol>
  
  <p>(It can also output the mutated pdb file by setting --out-pdb=true.)</p>
  
  <p>Either of those files can be used to determine change in free energy upon mutation (ddg ~ delta delta g).  The second file has additional information (the individual values used to calculate total change in free energy).</p>
  
  <p>Note: The FoldX website has very good documentation for each FoldX command.  You can find any additional information you require there.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-06-12T21:23:42.003Z
Answers.LastActivityDate: 2018-06-12T21:31:51.003Z
Answers.LastEditDate: 2018-06-12T21:31:51.003Z
Answers.UserDisplayName: Aalawlx
Body: <p>I have a series of proteins that are all phylogenetically related. Only one of this proteins is currently X-ray crystallized. </p>
  
  <p>Is it valid to modify the reference PDB file (I mean by hand) to match the mutated sequence? </p>
  
  <p>And if so, is there any way to minimize the protein energy to find the most likely structure? (Trying to match as closely as possible the reference structure)</p>
  
Tags: <proteins><protein-structure><pdb>
Title: If I modify a PDB file with a specific mutation, how to minimize energy?
ThreadText: # If I modify a PDB file with a specific mutation, how to minimize energy?
  
  (Asked by: username_0 on 2018-05-02)
  
  I have a series of proteins that are all phylogenetically related. Only one of this proteins is currently X-ray crystallized.
  
  Is it valid to modify the reference PDB file (I mean by hand) to match the mutated sequence?
  
  And if so, is there any way to minimize the protein energy to find the most likely structure? (Trying to match as closely as possible the reference structure)
  
  **Tags:** proteins, protein-structure, pdb
  
  
  ---
  
  ## 3 answers
  
  ### Answer by username_3 on 2018-06-12. Score: 4 | Accepted answer
  
  I would recommend you try the FoldX mutation engine. It is very easy to use and state of the art (it is the main competitor to expensive commercial software). It can do both of what you want - the mutation and energy optimisation. It may not work as well as molecular dynamics software, but for your simple task, it is a reasonable solution.
  
  A simple command line can do it:
  
  First, repair your structure (since, usually PDBs are not energetically optimised, so you could get better energy anyway without doing this):
  
  ```
  FoldX --command=repairpdb --pdb=1AJY.pdb
  
  ```
  
  Then, either use AlaScan or PositionScan to mutate. AlaScan only mutates to Alanine. PositionScan can mutate to any amino acid. More advanced mutations can be done with the BuildModel command.
  
  ```
  FoldX --command=PositionScan --pdb=1AJY_Repair.pdb --output-file=1AJY_Repair_MA30d --out-pdb=false --positions=MA30d
  
  ```
  
  Note: 'd' means try mutate the residue to 24 different amino acids; or, you can specify any of those 24, including the standard 20.
  
  Which produces two files:
  
  1. PS\_1AJY\_Repair\_MA30d\_scanning\_output.txt
  2. energies\_30\_1AJY\_Repair
  
  (It can also output the mutated pdb file by setting --out-pdb=true.)
  
  Either of those files can be used to determine change in free energy upon mutation (ddg ~ delta delta g). The second file has additional information (the individual values used to calculate total change in free energy).
  
  Note: The FoldX website has very good documentation for each FoldX command. You can find any additional information you require there.
  
  
  ---
  
  ### Answer by username_1 on 2018-06-12. Score: 5
  
  Yes you can modify the reference PDB file and look for the changes and for this purpose you need [visualizers.](https://www.rcsb.org/pdb/static.do?p=software/software_links/molecular_graphics.html) One among these is [Chimera](https://www.cgl.ucsf.edu/chimera/). You can easily carry out the energy minimization steps using Chimera, first mutate your PDB file manually and then provide it as an input for Chimera and the steps for doing energy minimization are fairly simple as given in Chimera docs [here.](https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/minimize/minimize.html)
  
  
  ---
  
  ### Answer by username_2 on 2023-11-13. Score: 2
  
  The previous answers suggest Chimera and FoldX, which require a license for commercial use.
  
  Commercial users can use [pdbfixer](https://github.com/openmm/pdbfixer/) or a similar tool to perform mutations (using e.g. [applyMutation](https://github.com/openmm/pdbfixer/blob/master/pdbfixer/pdbfixer.py#L729) in the API).
  
  You can then use [OpenMM](https://openmm.org/) to [refine the mutated structure](http://docs.openmm.org/latest/userguide/application/02_running_sims.html#a-first-example).
  
  Granted, this requires some comfort writing Python scripts to use the APIs.
  
  
  ---
  
  