<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-2">
<META content="MSHTML 6.00.2900.2627" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV>&nbsp;</DIV>
<DIV><FONT face=Arial size=2>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US 
style="FONT-SIZE: 10pt; FONT-FAMILY: Arial">Hi STEVE </SPAN><SPAN lang=EN-US 
style="FONT-SIZE: 10pt; FONT-FAMILY: Wingdings; mso-ascii-font-family: Arial; mso-hansi-font-family: Arial; mso-bidi-font-family: Arial; mso-char-type: symbol; mso-symbol-font-family: Wingdings"><SPAN 
style="mso-char-type: symbol; mso-symbol-font-family: Wingdings">J</SPAN></SPAN><SPAN 
lang=EN-US><?xml:namespace prefix = o ns = 
"urn:schemas-microsoft-com:office:office" /><o:p></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US 
style="FONT-SIZE: 10pt; FONT-FAMILY: Arial">I found it in mdrun.log.&nbsp; Look 
for&nbsp; "ewald correction".&nbsp;&nbsp; <o:p></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US 
style="FONT-SIZE: 10pt; FONT-FAMILY: Arial">Do you know how to evaluate this LR 
coulomb <o:p></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US 
style="FONT-SIZE: 10pt; FONT-FAMILY: Arial">interaction energy between protein 
and ligand?</SPAN><SPAN lang=EN-US><o:p></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US><FONT 
face="Times New Roman" size=3>David wrote that I missed some terms 
</FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US><FONT 
size=3><FONT face="Times New Roman">but I don&#8217;t know which ones!.<SPAN 
style="mso-spacerun: yes">&nbsp; </SPAN></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US><FONT 
face="Times New Roman" size=3>My evaluated energy <SPAN 
style="mso-spacerun: yes">&nbsp;</SPAN>seems wrong.</FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US><FONT 
face="Times New Roman" size=3>Please let me know if you find a right way to 
</FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US><FONT 
face="Times New Roman" size=3>do it.</FONT></SPAN></P><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: PL; mso-bidi-language: AR-SA">Michal&nbsp;<SPAN 
style="mso-spacerun: yes">&nbsp;&nbsp;</SPAN></SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV>Hi Michal,<BR><BR>It sounds like we are doing similar things.<BR><BR>Can 
you explain to me how you were able to use mdrun -rerun -debug to get<BR>the 
correction terms?&nbsp; I tried running this yesterday, and it generated<BR>some 
files (ctab.xvg, dtab.xvg, rtab.xvg, mdrun.log), but I have no idea<BR>where to 
look for the correction.<BR><BR>Thanks!<BR>Steve<BR><BR><BR>On Wed, 4 May 2005, 
Michal Kolinski wrote:<BR><BR>&gt; I tried to evaluate LR Coulomb interaction 
energy for ligand and protein.&nbsp; My system contains<BR>&gt; membrane, 
ligand, TM protein and ions.&nbsp; I did three reruns of my trajectory:<BR>&gt; 
1) leaving charges on ligand, zeroing rest<BR>&gt; 2) leaving charges on 
protein, zeroing rest<BR>&gt; 3) leaving charges on protein and ligand, zeroing 
rest.<BR>&gt; Then I calculated LR Coulomb interaction energy for prot-lig 
(Elr_c(ligand-protein)):<BR>&gt;<BR>&gt;<BR>&gt; Elr_c(protein-ligand) = 
Elr_c(ligand-ligand) + Elr_c(ligand-protein) 
+<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Elr_c(protein-protein)<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; in my 
case:<BR>&gt;<BR>&gt; 3 = 1 +&nbsp; Elr_c(ligand-protein) + 
2<BR>&gt;<BR>&gt;&nbsp; I also included long range corrections for each of the 
reruns that I found using -rerun -debug option.<BR>&gt; Is this a right way to 
estimate LR coulomb interaction energy between two molecules, my estimated 
energy seems wrong&nbsp; 13.567.<BR>&gt;&nbsp; Please give me some comment on 
this.<BR>&gt; Thank you in 
advance.<BR>&gt;<BR>&gt;<BR><BR><BR></DIV></BODY></HTML>