Previous Thread
Next Thread
Print Thread
Joined: Aug 2016
Posts: 6
A
Forum Member
OP Offline
Forum Member
A
Joined: Aug 2016
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.

Attached Images
ser_plp.png (59.55 KB, 297 downloads)
Joined: Sep 2003
Posts: 8,660
Likes: 26
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,660
Likes: 26
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

Joined: Aug 2016
Posts: 6
A
Forum Member
OP Offline
Forum Member
A
Joined: Aug 2016
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

Last edited by Aravinda; 07/14/19 12:47 AM.

Moderated by  alex, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u5 Page Time: 0.009s Queries: 21 (0.005s) Memory: 0.7407 MB (Peak: 0.8039 MB) Data Comp: Off Server Time: 2023-12-10 13:28:02 UTC
Valid HTML 5 and Valid CSS