Re: Unable to reproduce same results in minimization with GBIS

From: Chris Harrison (charris5_at_gmail.com)
Date: Thu Sep 13 2012 - 11:50:01 CDT

Vinod,

"Kasam, Vinod" <vkasam_at_coh.org> writes:
> I am running a minimizationin NAMD using the parameters shown below. When I ran minimization using the same parameters in two different folders with the same input (protein), I am ending up in different results.
> The co-ordinates of the proteins and the energy values after minimization are different.

How different? Please attempt to quantitate it? ie, rmsd plots,
energy difference plot, etc.

Chris

> However, if I do not use GBIS (comment out GBIS input parameters), the results are same and just FINE.
>
> I am new to NAMD and I am using NAMD2.8 version, can some explain me why are the results different for the same input if I use GBIS.
>
> Thank you,
> Vinod
>
>
>
>
> # MOLECULAR SYSTEM
> structure min_WT-chain_D.psf
> coordinates min_WT-chain_D.pdb
>
>
> # FORCE FIELD
> paraTypeCharmm on
> parameters par_all22_prot_cmap.inp
>
>
> # APPROXIMATIONS
> switching on
> switchdist 15
> cutoff 16
> pairlistdist 17.5
> #margin 3
> exclude scaled1-4
> 1-4scaling 1.0
> dielectric 3
> seed 20111012
>
> #INTEGRATOR
> timestep 1.0
> nonbondedFreq 1
> fullElectFrequency 2
> stepspercycle 20
>
>
> #GBIS - Implicit solvent
> GBIS on
> solventDielectric 78.5
> ionConcentration 0.2
> alphaCutoff 15
>
>
> # PERIODIC BOUNDARY CONDITIONS
> cellBasisVector1 51.76400065422058 0 0
> cellBasisVector2 0 45.07099914550781 0
> cellBasisVector3 0 0 43.320000648498535
> cellOrigin 28.250900268554688 39.77300262451172 7.405788421630859
>
>
> # WRAPPING
> wrapAll on
>
>
>
> # OUTPUT
> outputname monomerA-crystal-min
> dcdfreq 500
> xstfreq 10000
>
> outputenergies 10
> outputtiming 100
> binaryoutput no
> binaryrestart yes
>
>
> # RESTART FILES
> restartname res
> restartfreq 1000
> restartsave yes
>
>
> # TEMPERATURE CONTROL
> temperature 0
>
>
> minimize 1000
>
>
>
> ---------------------------------------------------------------------
> *SECURITY/CONFIDENTIALITY WARNING:
> This message and any attachments are intended solely for the individual or entity to which they are addressed. This communication may contain information that is privileged, confidential, or exempt from disclosure under applicable law (e.g., personal health information, research data, financial information). Because this e-mail has been sent without encryption, individuals other than the intended recipient may be able to view the information, forward it to others or tamper with the information without the knowledge or consent of the sender. If you are not the intended recipient, or the employee or person responsible for delivering the message to the intended recipient, any dissemination, distribution or copying of the communication is strictly prohibited. If you received the communication in error, please notify the sender immediately by replying to this message and deleting the message and any accompanying files from your system. If, due to the security risks, you do not wi!
> sh to receive further communications via e-mail, please reply to this message and inform the sender that you do not wish to receive further e-mail from the sender. (fpc5p)
> ---------------------------------------------------------------------
>
>

Best,
Chris

--
Chris Harrison, Ph.D.
NIH Center for Macromolecular Modeling and Bioinformatics
Theoretical and Computational Biophysics Group
Beckman Institute for Advanced Science and Technology
University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801
http://www.ks.uiuc.edu/Research/namd       Voice: 773-570-6078
http://www.ks.uiuc.edu/~char               Fax:   217-244-6078

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:22:33 CST