|
|
(22 intermediate revisions not shown) |
Line 1: |
Line 1: |
- | {{EPF-Lausanne09}}
| |
- | <div CLASS="epfltrick">__TOC__
| |
- | </div><div CLASS="epfl09">
| |
| | | |
- |
| |
- | We start with the PDB file of the protein, obtained through the Protein Data Bank.
| |
- |
| |
- | 1. We open VMD, and lauch our protein
| |
- |
| |
- |
| |
- | 2. Creation of the pdb file
| |
- | In the Tk Console menu of VMD we open the VMD TkCon window, and type the following commands:
| |
- | : set our_protein [atomselect top not water and not GOL]
| |
- | : $our_protein writepdb our_protein.pdb
| |
- |
| |
- | This will select all the protein except the water and except the glycerol, with the cofactor.
| |
- | We have created the file our_protein.pdb, which contains the coordinates of the protein alone without hydrogens.
| |
- | Quit VMD.
| |
- |
| |
- |
| |
- | 3. Creation of the psf file
| |
- | There is 2 ways to create a psf:
| |
- | * in VCMD, Extensions → Modeling → Automatic PSF Builder
| |
- |
| |
- | * We first make a pgn file, which will be the target of psfgen.
| |
- | In a Terminal window, open a text editor and type:
| |
- | : package require psfgen
| |
- | : topology top all27 prot lipid.inp
| |
- | : pdbalias residue HIS HSE
| |
- | : pdbalias atom ILE CD1 CD
| |
- | : segment U
| |
- | : {pdb ubqp.pdb}
| |
- | : coordpdb ubqp.pdb U
| |
- | : guesscoord
| |
- | : writepdb ubq.pdb
| |
- | : writepsf ubq.psf
| |
- |
| |
- | In a Terminal window, type the following command:
| |
- | > vmd -dispdev text -e our_protein.pgn
| |
- | This will run the package psfgen on the file ubq.pgn and generate the psf and the pdb file of ubiquitin with hydrogens.
| |
- | A new pdb file with the complete coordinates of all atoms is written, including H; and a psf file with the complete structural information of the protein.
| |
- |
| |
- | 4. Solvating the Protein
| |
- | Now, the protein needs to be solvated, i.e., put inside water, to more closely resemble the cellular environment. This will be done by placing the protein in a water box, in preparation for minimization and equilibration with periodic boundary conditions.
| |
- |
| |
- | In the VMD Main window, open the Tk Console, and type:
| |
- | : package require solvate
| |
- | : solvate our_protein.psf our_protein.pdb -t 5 -o our_prot_in_a_water_box
| |
- |
| |
- |
| |
- |
| |
- |
| |
- | </div><div CLASS="epfl09bouchon"></div>
| |