Topic Options
#37567 - 07/13/19 12:49 PM parameters for pyridoxal-5'-phosphate (PLP)-serine
Aravinda Offline
Forum Member

Registered: 08/23/16
Posts: 6
Hi,
I want to derive parameters for Ser-PLP (see attachment). For this first, I tried CGenFF web site and got a really high penalty. Also, it failed to identify the Imine bond (C=N). And gave N atom type of NG311. I have used PDB as well as mol2 with connectivity information.

With one of the advices I received from an earlier post (https://www.charmm.org/ubbthreads/ubbthreads.php?ubb=showflat&Number=37536#Post37536) I tried using parameters from analog compounds as much as I can.

For pyridinium, I was able to find parameters from two places. (1) 3-hydroxypyridine and ethylpyridine from toppar_all36_prot_pyridines.str, and also from all36_cgenff.rtf for protonated pyridine (CG2R62). I am planning to derive charges and other parameters for these compounds. However, what atom type would be most reasonable to take out of these two options. In 1st option, I dont have the protonated version. But in the second option I dont have other substituents. Any assistant is highly appreciated.


Attachments
ser_plp.png (18 downloads)


Top
#37568 - 07/13/19 02:33 PM Re: parameters for pyridoxal-5'-phosphate (PLP)-serine [Re: Aravinda]
rmv Online   content

Forum Member

Registered: 09/17/03
Posts: 8373
Loc: 39 03 48 N, 77 06 54 W
Note that the .mol2 file must have the correct bond order for the C=N bond, or the CGenFF server may be confused. The H atom placement in the molecule, esp. protonation states, must also be correct.

A high penalty doesn't necessarily mean the estimated parameters are wrong, just that the analogy to existing validated parameters is not very good. It means at least testing the estimated parameters, which should be done for whatever approach is used to get a first test set, CGenFF or chemical intuition. Multiple parameter sets are often tested in simulations, so it is not unreasonable to test both choices and the CGenFF estimates.

It should be noted that QM calculations are encouraged, if possible, when the penalties are very large.

The problem for testing is often the difficulty in finding good target data from experiments that help refine the conformational preferences in the appropriate environment, e.g. water and/or a protein binding site. Aggressive and persistent searching of the literature and other sources for the whole molecule or relevant fragments is often needed to find some useful data.

If the molecule will form an important part of the research project, the decisions made on these parameters are critical. Some additional, unplanned effort may be required in order to ensure producing a respectable result for the parameters.
_________________________
Rick Venable
computational chemist


Top
#37569 - 07/13/19 08:46 PM Re: parameters for pyridoxal-5'-phosphate (PLP)-serine [Re: Aravinda]
Aravinda Offline
Forum Member

Registered: 08/23/16
Posts: 6
Thank you very much sir!! I did discover why CGenFF server got confused based on your feedback. I manually adjusted the bond order earlier,and hence added the bond order of 2 for C=N and others accordingly. But I realized I should add "ar" in the bond order column since the pyridine system is a conjugated one. Updated mol2 file gave me CG2R62 for the C and others as I expected.

> Multiple parameter sets are often tested in simulations, so it is not unreasonable to test both choices and the CGenFF estimates.

Okay sir. I will test both types of parameters to evaluate them.

PLS-Ser does play an important role. Because of this, I did design a positive control and a negative control with our experimental collaborators to validate our simulation models.

Hopefully, through that I might be able to ensure the accuracy of the simulation setup (fingers crossed).

Thank you very much for your kind feedback (as always) on this.
Best,
Aravinda


Edited by Aravinda (07/13/19 08:47 PM)

Top

Moderator:  alex, chmgr, John Legato