<div dir="ltr">Thanks for your help! Unfortunately, in my hands, this as well produces a phosphate that blows up, both for your exact example and also as applied to my cases. I tried uninstalling SwissSidechain as well, but to no effect. <div><br></div><div>Generally, the phosphate becomes agitated with severe conformations and steric clashes (e.g. you see things like one O bonding to another O to form a three-membered OPO ring, but with a hydrogen superimposed on an O). The agitation increases across the molecule until it "explodes" into pieces.</div><div><br></div><div>I'm going to work with the well-behaved dibasic phosphate for now.... I'm just recording my experience here in the email record in case it helps someone else.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Mar 30, 2019 at 4:44 PM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Geoffrey,<br>
I don’t think anybody has tried the Chimera dynamics tool with phosphorylated sidechains.  The minimization/MD tools in Chimera are fairly limited and slow (no cutoffs or sophisticated handling of long-range electrostatics, no detailed control over residue parameters) compared to dedicated programs like GROMACS and AMBER, so most detailed research with MD, as opposed to brief explorations of flexibility or structure cleanup, would use one of those dedicated packages instead.<br>
<br>
I can answer the part about protonation and adding charges in Chimera, but whether it will attain stability in simulation may be a separate matter. <br>
<br>
Before you add hydrogens, just change the automatically assigned atom type of the oxygen you wish to protonate to O3.  I can see for example with the TPO in 1guf chain A that the three terminal oxygens are automatically given type O3-<br>
<br>
open 6guf<br>
delete ~protein<br>
delete ~:.a<br>
disp :tpo<br>
focus :tpo<br>
labelopt info idatmType<br>
label :tpo<br>
<br>
(… select the oxygen you want to change with Ctrl-click)<br>
<br>
setattr a idatmType O3 sel<br>
label sel<br>
<br>
addh<br>
addcharge   <br>
(...TPO is now estimated as charge -1, OK)<br>
<br>
I tried 100 steps steepest descent minimization, selecting a zone within 4A of TPO and allowing only those atoms to move, and it didn’t move much.<br>
<br>
I hope this helps,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D.                       <br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Mar 30, 2019, at 11:31 AM, Geoffrey Sametz <<a href="mailto:sametz@udel.edu" target="_blank">sametz@udel.edu</a>> wrote:<br>
> <br>
> I am having my students model dipeptides where one residue is a phosphorylated serine or threonine. Molecular dynamics with the dibasic phosphate works fine. However, the student NMR data for their dipeptides is acquired at low pH, so I would like them to be able to model the monobasic form. Unfortunately, all my attempts at creating a monobasic form have resulted in the molecule "blowing up" in MD, first by crazily shaking the phosphate as if it wants to dislodge the proton, and then the entire molecule until it fragments.<br>
> <br>
> I have tried building these from scratch by various methods (e.g. those listed in  <a href="http://plato.cgl.ucsf.edu/pipermail/chimera-users/2012-January/007068.html" rel="noreferrer" target="_blank">http://plato.cgl.ucsf.edu/pipermail/chimera-users/2012-January/007068.html</a> ) as well as using SwissSidechain and mutating residues. I have tried exporting as a .pdb, opening in PyMOL to make sure everything looks OK, resaving as .pdb and importing PyMOL's version. <br>
> <br>
> I have also used a protein from PDB that had a TPO residue assigned, and snipped it down to a dipeptide.For both this manner of construction, and for building the dipeptide from scratch with SwissSidechain, Chimera interprets SwissSidechain's TPO (pThr monobasic) and SEP (pSer monobasic) as their dibasic forms (TPO2 and SEP2 respectively). PyMOL interprets these structures with the correct level of protonation, however. <br>
> <br>
> Is there a known issue with Chimera and treatment of [-OPO3H]-1 groups, or have I managed to make the same error across multiple build methods?<br>
> <br>
> -- <br>
> Dr. Geoffrey Sametz<br>
> QDH 104<br>
> Department of Chemistry and Biochemistry<br>
> University of Delaware<br>
> <a href="mailto:sametz@udel.edu" target="_blank">sametz@udel.edu</a><br>
> (302) 831-3621<br>
> _______________________________________________<br>
> Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br>
> Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" rel="noreferrer" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Dr. Geoffrey Sametz</div><div>QDH 104<br>Department of Chemistry and Biochemistry<br>University of Delaware<br><a href="mailto:sametz@udel.edu" target="_blank">sametz@udel.edu</a></div><div>(302) 831-3621</div></div></div></div></div>