<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div><blockquote type="cite" class=""><div class="">On Feb 8, 2017, at 9:07 AM, Stellamaris Sotomayor Burneo <<a href="mailto:stellamaris_sb@yahoo.com" class="">stellamaris_sb@yahoo.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class=""><div style="background-color: rgb(255, 255, 255); font-family: verdana, helvetica, sans-serif; font-size: 13px;" class=""><div id="yiv0449843958" class=""><div id="yui_3_16_0_ym19_1_1486571803181_22378" class=""><div id="yui_3_16_0_ym19_1_1486571803181_22377" style="background-color: rgb(255, 255, 255); font-family: verdana, helvetica, sans-serif; font-size: 13px;" class=""><div id="yiv0449843958" class=""><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_7640" class=""><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_7639" style="background-color: rgb(255, 255, 255); font-family: verdana, helvetica, sans-serif; font-size: 13px;" class=""><div id="yiv0449843958yui_3_16_0_ym19_1_1486571020879_2956" class=""><br class=""></div><div dir="ltr" id="yiv0449843958yui_3_16_0_ym19_1_1486571020879_6450" class="">Dear Dr. Eric,</div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571020879_6487" class=""><br class=""></div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571020879_6532" class=""><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_7880" class="">My name is Stellamaris. I am an undergraduate student from Ecuador. There is almost no one working in molecular docking in my country so far. So, it is very difficult to find somebody to ask for help. I would like to kindly ask you a question about the problem I am having with Chimera and I would really appreciate your help.</div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_7881" class=""><br class=""></div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_8034" dir="ltr" class=""> I am trying to convert a set of receptors (named receptor.pdb) to receptor.mol2. using chimera. However i am having problems with the charges of FAD, NAD, NDP, and FMN. Receptors that have problems in those groups are not being transformed.</div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_8035" dir="ltr" class=""><br class=""></div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_8036" dir="ltr" class=""> Is there a way to transform those files without having to change to charges? I have tried changing the charges and it seems to work. However, I am very afraid of changing the conditions of my simulations if I change the charges.</div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571803181_8080" dir="ltr" class=""><br class=""></div><div id="yiv0449843958yui_3_16_0_ym19_1_1486571020879_6533" class=""><div id="yui_3_16_0_ym19_1_1486571803181_22392" class="">I have gotten the following message: <br class=""></div><div dir="ltr" id="yui_3_16_0_ym19_1_1486571803181_22449" class="">NAD) Error: cannot run "/opt/UCSF/Chimera64-1.10.2/bin/amber14/bin/sqm -O -i sqm.in -o sqm.out" of bcc() in charge.c properly, exit<br class=""></div><div id="yui_3_16_0_ym19_1_1486571803181_22450" class=""> Failure running ANTECHAMBER for residue NAD</div></div></div></div></div></div></div></div></div></div></div></div></blockquote><br class=""></div><div>Hi Stellamaris,</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>As you have discovered, the underlying tools that Chimera relies on to derive charges (AmberTools/antechamber/sqm) sometimes fail for entities with concentrated areas of negative charge. It is for this reason that Chimera already includes pre-calculated charges for ADP, ATP, GDP, and GTP. These charges were taken from the <a href="http://research.bmh.manchester.ac.uk/bryce/amber" class="">AMBER parameter database</a>. Charges for the moieties you mentioned are also in that database, and I’ve gone ahead and added NAD and NDP to the ones that Chimera knows about without needing to do a calculation. Precalculated NAD/NDP charges will be available in daily builds dated Feb. 9 or later.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>Unfortunately, FMN and FAD are more problematic in that Chimera doesn’t protonate them in the standard way used in the database. It adds a proton (and positive charge) to the isoalloxazine ring system so that the central ring becomes aromatic. Obviously, the simplest way to deal with this if you don’t need FAD/FMN as a cofactor is to simply delete them. If you do need them then things get hairy. I guess what I would do is:</div><div><br class=""></div><div>1) Delete FMN/FAD.</div><div>2) Run Dock Prep which will get you a Mol2 file of the remainder of the system with charges.</div><div>3) Close your system and reopen it.</div><div>4) Delete everything but the FMN/FAD.</div><div>5) Set the IDATM types of the FMN/FAD so that they correspond to the standard protonation pattern.</div><div>6) Add hydrogens.</div><div>7) Write another Mol2 file with just FMN/FAD.</div><div>8) Edit this second Mol2 file and transcribe the charges found in the corresponding Amber parameter database file to the atoms.</div><div>9) Combine your two Mol2 files.</div><div><br class=""></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>Step 6 is a little tricky in that the atom names in the “OFF” parameter files aren’t all the same as the PDB standard names, but there is connectivity information in the OFF file, so you can figure out which atom is which.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>You can do step 5 by running the “defattr” command with the file I have attached, which will set the types of FMN/FAD heavy atoms to the standard protonation pattern.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>Also, instead of step 8, you could work up another “defattr” file to assign charges to the FMN/FAD before you write out the Mol2 file. Probably easier.</div><div><br class=""></div><div>—Eric</div><div><br class=""></div><div><div class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><span class="Apple-tab-span" style="white-space: pre;"> </span>Eric Pettersen</div><div class=""><span class="Apple-tab-span" style="white-space: pre;"> </span>UCSF Computer Graphics Lab</div><div class=""><br class=""></div></div><br class="Apple-interchange-newline"></div></div></div></body></html>