Re: [AMBER] Parameterize SO4

From: Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu>
Date: Fri, 21 Jan 2022 00:57:57 -0500

Zizhang,

It was a while ago that I worked on my project with ammonium sulfate. Looking back on the files, what I apparently did is follow the tutorial for parameterizing a non-standard residue  <https://ambermd.org/tutorials/basic/tutorial4b/index.php>(RESP charge fitting and assignment of GAFF atom types). You could improve on this by editing the frcmod or lib files to reflect literature values. I’ll send you all the files I used to construct the frcmod and lib files off-list just so you have the full story, but please give some thought to whether the GAFF parameterization is good enough or if you need to refine what I am going to share.

Best,
Matthew

> On Jan 20, 2022, at 10:24 PM, zz sheng <shengzizhang.gmail.com> wrote:
>
> *Message sent from a system outside of UConn.*
>
>
> Dear Anselm, Matthew, and David,
>
> Thanks very much for the suggestions. I attached the pdb file for the structure. I am not very sure if the SO4 should be there for MD but it appears in the structure and is very close to protein-protein interface.
>
> The literature is very helpful.
> Matthew, do you mind sending me the parameter files from the literature (.lib and .frcmod) for tleap? I truly have no idea how to generate them.
>
>
> Best,
>
> Zizhang
>
>
>
>> On Jan 20, 2022, at 7:50 AM, Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu> wrote:
>>
>> Seconding Dr. Case’s points, you can check out Phys. Chem. Chem. Phys., 2017, 19, 20593.
>> I think I used the parameters reported there when modeling the crystal structure of a protein with sulfates at the protein-protein interface.
>>
>> Best,
>> Matthew
>>
>>
>>
>>> On Jan 20, 2022, at 7:05 AM, David A Case <david.case.rutgers.edu <mailto:david.case.rutgers.edu>> wrote:
>>>
>>> *Message sent from a system outside of UConn.*
>>>
>>>
>>> On Thu, Jan 20, 2022, zz sheng wrote:
>>>>
>>>> My protein-protein interface has a SO4, I wonder how to parameterize it for
>>>> MD? Could I use antechamber to generate the topology file and treat it as a
>>>> 'ligand'? I saw ions have different ways of parameterization in the manual.
>>>
>>> First, are you sure the sulphate is present under the conditions you are
>>> interested in? Unless you are simulating the crystal itself, that seems
>>> unlikely.
>>>
>>> Second, yes, antechamber will give you a good starting point. One could try
>>> to optimize those values based on experimental ion-water data. That is
>>> unlikely to make a big difference, and the choice depends on how important
>>> the sulphate is to your simulation.
>>>
>>> ....dac
>>>
>>>
>>> _______________________________________________
>>> AMBER mailing list
>>> AMBER.ambermd.org
>>> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C397468ca0d0d477e2c1e08d9dc8d9d4e%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637783323045596582%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=wWJFDmUE3FNVGNerwzaQN%2Ff%2FLMvsHGrdOsDzMJ1ezsc%3D&amp;reserved=0 <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C397468ca0d0d477e2c1e08d9dc8d9d4e%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637783323045596582%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=wWJFDmUE3FNVGNerwzaQN%2Ff%2FLMvsHGrdOsDzMJ1ezsc%3D&amp;reserved=0>
>>
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org <mailto:AMBER.ambermd.org>
>> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C397468ca0d0d477e2c1e08d9dc8d9d4e%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637783323045596582%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=wWJFDmUE3FNVGNerwzaQN%2Ff%2FLMvsHGrdOsDzMJ1ezsc%3D&amp;reserved=0 <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C397468ca0d0d477e2c1e08d9dc8d9d4e%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637783323045596582%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=wWJFDmUE3FNVGNerwzaQN%2Ff%2FLMvsHGrdOsDzMJ1ezsc%3D&amp;reserved=0>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C397468ca0d0d477e2c1e08d9dc8d9d4e%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637783323045596582%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=wWJFDmUE3FNVGNerwzaQN%2Ff%2FLMvsHGrdOsDzMJ1ezsc%3D&amp;reserved=0

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jan 20 2022 - 22:00:02 PST
Custom Search