mutating_a_lipid
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
mutating_a_lipid [2019/09/23 20:32] – edit | mutating_a_lipid [2024/02/20 16:53] (current) – admin | ||
---|---|---|---|
Line 4: | Line 4: | ||
======Mutating a lipid====== | ======Mutating a lipid====== | ||
- | Author: Eric Wang | + | Authors: Eric Wang & Jeff Klauda |
- | Please email me at //ericzwang@gmail.com// if something isn't clear | + | Please email me at //jbklauda@umd.edu// if something isn't clear |
Mutating a lipid means changing its chemical structure. We do this when our membrane contains a lipid which is not available by default in the CHARMM force field, but the structure is similar to an existing lipid. To mutate a lipid, you'll need to create your membrane with the pre-mutated lipids using CHARMM-GUI first. | Mutating a lipid means changing its chemical structure. We do this when our membrane contains a lipid which is not available by default in the CHARMM force field, but the structure is similar to an existing lipid. To mutate a lipid, you'll need to create your membrane with the pre-mutated lipids using CHARMM-GUI first. | ||
- | Once that is done, modify the relevant .str file to include your mutated lipid. See the following path on DT2 for an example for mutating CER240 to CER6. | + | Once that is done, modify the relevant .str file to include your mutated lipid. See the following path on ZT1 for an example for mutating CER240 to CER6. For small mutations, the IC table isn't necessary. \\ |
- | ''/ | + | ''/ |
+ | |||
+ | This mutation involved adding a few hydroxyls. When you add new atoms, you need to decide what the partial charges will be. You also may need to add new parameters for dihedrals or angles at the end of the .str file. CHARMM will give you an error if these are missing. | ||
+ | |||
+ | Mutate the lipids using **mutate.inp** An example is in the path\\ | ||
+ | ''/ | ||
+ | The lipid names and how the atoms are deleted/ | ||
+ | Execute mutate.inp in CHARMM \\ | ||
+ | '' | ||
+ | '' | ||
+ | |||
+ | Note that you will need to use the *inp files in this directory to build the files needed to run the MD simulation. The current CHARMM-GUI files will not work. | ||
+ | |||
+ | Run **step5_assembly.inp** to create the step5 input files. | ||
+ | |||
+ | You need to edit the *restraint.namd.str file to add your new lipid/s. Copy and paste an if/endif block of a similar lipid and rename it as a new lipid. If you have a more complex lipid you might also need to edit the *restraint2.namd.str file, so discuss this with Dr. Klauda. | ||
+ | |||
+ | Then run **step5_input.inp** to create the namd input files. If you look in **step5_assembly.pdb**, | ||
+ | |||
+ | Finally, you will need to go into the namd subdirectory and validate that the restraints subdirectory has files with your new lipids listed. Moreover, you will need to make edits to the '' | ||
+ | |||
mutating_a_lipid.1569285122.txt.gz · Last modified: 2019/09/23 20:32 by edit