<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
The problem with prodrg is that it removes all my hydrogen atoms and
that it will not accept mocules of more than 2 parts. A rotaxane is a
molecule which consists of a linear part, mostly saturated hydorcarbon
(the thread) and a ring which goes around the thread. This ring (aka
macorcycle) is held in place at one position by hydrogen bonds the tow
ends of the thread are then blocked by bulky groups&nbsp; and the ring is
mechanically locked onto the thread. (see Brouwer et al SCIENCE 291
(5511): 2124-2128 MAR 16 2001). By applying a photochemical stimulus
the macrocycle moves from one position to another (since the hydrogen
bond affinity is changed) and back again.<br>
<br>
So i need all the polar hydrogen atoms and prodrg deletes them. Writing
a topology by hand is I think a rather time consuming job for a system
of more then 250 atoms. This besides the fact that i would'nt really
know how to write a topology from scratch.&nbsp; Or could i take the
topology from prodrg and easily add some hydrogen atoms?<br>
<br>
Peter<br>
<br>
David van der Spoel wrote:<br>
<blockquote cite="mid1096631179.26920.1.camel@vangogh.bmc.uu.se"
 type="cite">
  <pre wrap="">On Fri, 2004-10-01 at 12:07, Peter Zoon wrote:
  </pre>
  <blockquote type="cite">
    <pre wrap="">Hi all,

Here's the story. We wrote a program which converts trajectories from 
Macromodel to xtc. We obtain a initial gro file and the trajectory in 
xtc format. We would like to use the gromacs analysis routines to 
analyze the trajectory. We a looking at a rotaxane in particular we want 
to study the hydrogen bonds between the macrocycle and the thread and 
the dihedrals within the thread.
    </pre>
  </blockquote>
  <pre wrap=""><!---->the problem I presume it the topology. How about generating one using
prodrg? Many analysis tools can substitute the tpr file for a gro/pdb
file but g_hbond needs the atom names and the connectivity.
g_dih actually also needs the topology, although some things can be done
with just an index file (g_angle, which also analyses dihedrals)
  </pre>
  <blockquote type="cite"><br>
  </blockquote>
</blockquote>
</body>
</html>