From chemistry-request -x- at -x- server.ccl.net Sat Dec 15 17:29:52 2001
Received: from bute.st-andrews.ac.uk (im17 # - at - # [138.251.12.1])
	by server.ccl.net (8.11.6/8.11.0) with ESMTP id fBFMTqi31935
	for <chemistry -8 at 8- ccl.net>; Sat, 15 Dec 2001 17:29:52 -0500
Received: from localhost (im17 |-at-| localhost)
	by bute.st-andrews.ac.uk (8.9.1a/8.9.1) with SMTP id WAA02981
	for <chemistry:~at~:ccl.net>; Sat, 15 Dec 2001 22:29:32 GMT
Date: Sat, 15 Dec 2001 22:29:32 +0000 (GMT)
From: Ibrahim Moustafa <im17 ":at:" st-andrews.ac.uk>
X-Sender: im17 ^at^ st-andrews.ac.uk
To: chemistry: at :ccl.net
Subject: Molecular dynamics with implicit solvent
Message-ID: <Pine.SOL.3.96.1011215215908.2325B-100000 \\at// bute>
MIME-Version: 1.0
Content-Type: TEXT/PLAIN; charset=US-ASCII

 Dear ccl users,

   I'm doing MD on my protein structure without explicit solvents because
of the limiting computing power. I use Tripos FF, in SYBYL6.6; cutoff
distance 15 A, distance dependent dielectric function with dielectric
constant 2.

   First, I heat the system to 300 K in steps ( 10, 50, 100, 200, 300 K),
2 ps for each step.
   
  Comparing the resultant configuration of the heating step with the
starting configuration gives ( rmsd 3.79 A ). Of course the absence of
solvents is responsible for this nonsense movements. Now I'm thinking to
do the simulation at low temperature ( 50 K ) to supress these large
movements. I want to discuss this with ccl users before doing so.



  What I want to ask about is: suppose i conduct this simulation at 50 K
instead of 300 k, non physiological condition, ( I have seen some people
doing that " Journal of molecular structure (theochem) 368;1996: 205-212")
Can we rely on the output results of that kind of simulation? can we realy
trust the answer for the question under investigation?

  -Another question: in absence of solvents do we need to equilebrate the
system before the production phase in the MD? if so, WHY?

  - Is it better to start the simulation from the average structure of the
heating step ( or the equilebration step ) or start from the last
configuration of each step??

   I'd be really appreciated if someone answered these questions. I
promise I'll sumarize the answer to the CCL users.

  Many thanks,
   Ibrahim
 

Name         :Ibrahim M.Moustafa
Mail address :Center for Biomolecular Science,
              BIOMOLECULAR SCIENCE BUILDING,
              North Haugh,St-Andrews,
              Fife,KY16 9ST
              Scotland,U.K.
         Tel :+44(0)1334-467257
      E-mail :im17 -8 at 8- st-andrews.ac.uk