From owner-chemistry@ccl.net Sat Apr 16 10:52:00 2011 From: "Giovanni Finoto Caramori caramori(_)qmc.ufsc.br" To: CCL Subject: CCL: problem in EDA calculation Message-Id: <-44389-110415171336-22117-4FMTU538Fn6zfcXyHph0OA,,server.ccl.net> X-Original-From: Giovanni Finoto Caramori Content-Transfer-Encoding: 7bit Content-Type: multipart/alternative; boundary="=_5uy8r2lgm34y" Date: Fri, 15 Apr 2011 17:41:06 -0300 MIME-Version: 1.0 Sent to CCL by: Giovanni Finoto Caramori [caramori||qmc.ufsc.br] This message is in MIME format. --=_5uy8r2lgm34y Content-Type: text/plain; charset=windows-1256 Content-Description: Versao da Mensagem em Texto Puro Content-Disposition: inline Content-Transfer-Encoding: quoted-printable Dear Kaul, In your output check the SCF convergence. Probably it is unconverged. look f= or a sentence like this: SCF IS UNCONVERGED, TOO MANY ITERATIONS If the SCF is unconverged the wave functions is not reliable! regards, Giovanni Citando "indu kaul indu.kaul _ iiserpune.ac.in" : > > Sent to CCL by: "indu=A0 kaul" [indu.kaul,iiserpune.ac.in] > Dear All, > > I have newly started using GAMESS for EDA analysis of certain dimers. > 1)I dont know the representation of writing 6-311++g** basis,can=20 > somebody tell how to do it.?? > 2)I ran EDA for a dimer at b3lyp/aug-cc-pvdz and the result shows: > ALL BASIS SET=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 = =A0 =A0 =A0 =A0 =A0HARTREE=A0 =A0 =A0 =A0 =A0 KCAL/MOL >=A0 ------------- >=A0 ELECTROSTATIC ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 ES=3D=A0 =A0 = =A0 =A0 =A0-1.236592=A0 =A0 =A0 =A0 =A0=20 >=A0 -775.97 >=A0 EXCHANGE ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0EX=3D=A0 = =A0 =A0 =A0 =A0-1.002818=A0 =A0 =A0 =A0 =A0=20 >=A0 -629.28 >=A0 REPULSION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0REP=3D=A0 = =A0 =A0 =A0 =A0 5.291539=A0 =A0 =A0 =A0 =A0=20 >=A0 3320.49 >=A0 POLARIZATION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 POL=3D=A0 =A0 = =A0 =A0 =A0-3.962384=A0 =A0 =A0 =A0 =A0=20 > -2486.44 >=A0 DFT DISPERSION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0DISP=3D=A0 =A0 =A0 = =A0 =A0-0.095027=A0 =A0 =A0 =A0 =A0=20 >=A0 =A0-59.63 >=A0 TOTAL INTERACTION ENERGY HF OR DFT=A0 =A0 =A0E=3D=A0 =A0 =A0 =A0 =A0-1.= 005283=A0 =A0 =A0 =A0 =A0=20 >=A0 -630.83 > > How can the electrostatic energy be so huge..!! > AM i doing something wrong in giving input??..The input i gave is :- >=A0 $CONTRL SCFTYP=3DRHF RUNTYP=3DEDA DFTTYP=3DB3LYP MAXIT=3D200 >=A0 =A0 =A0MULT=3D1 ISPHER=3D1 COORD=3DUNIQUE NOSYM=3D1 $END >=A0 $SYSTEM TIMLIM=3D525600=A0 MWORDS=3D50 MEMDDI=3D500 $ENDI >=A0 $BASIS GBASIS=3DACCD $END >=A0 $SCF DIRSCF=3D.TRUE. CONV=3D1.0D-04=A0 $END >=A0 $LMOEDA MATOM(1)=3D15,11 MCHARG(1)=3D0,0 >=A0 =A0 =A0 MMULT(1)=3D1,1 $END >=A0 $DATA > > It will be a great pleasure if somebody can help me out of it.? > Thanks In Advance............. > Indu Kaul > > > > -=3D This is automatically added to each message by the mailing script =3D= ->=A0 =A0 =A0 =A0http://www.ccl.net/cgi-bin/ccl/send_ccl_message>=A0 =A0 =A0 =A0http://www.ccl.net/cgi-bin/ccl/send_ccl_message>=A0 =A0 =A0 =A0http://www.ccl.net/chemistry/sub_unsub.shtml>=A0 =A0 =A0 =A0http://www.ccl.net/spammers.txt> > > =A0=20 ______________________________________ Giovanni Finoto Caramori PhD Departamento de Quimica - CFM-UFSC Campus Universitario Trindade - C.P. 476 88040-900 Florianopolis/SC http://www.qmc.ufsc.br/~caramori Phone: ++55(48)3721-6844/6845/6846 Fax: ++55(48)3721-6852 Ramal 239 Mobile: ++55(48)9926-8342 ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. --=_5uy8r2lgm34y Content-Type: text/html; charset=windows-1256 Content-Description: Versao HTML da Mensagem Content-Disposition: inline Content-Transfer-Encoding: quoted-printable Dear Kaul,

In your output check the SCF convergence. Probably it is unconverged. look f= or a sentence like this:

SCF IS UNCONVERGED, TOO MANY ITERATIONS

If the SCF is unconverged the wave functions is not reliable!

regards,

Giovanni

Citando "indu kaul indu.kaul _ iiserpune.ac.in" <owner-chemistr= y _ ccl.net>:

>
> Sent to CCL by: "indu=A0 kaul" [indu.kaul,iiserpune.ac.in] > Dear All,
>
> I have newly started using GAMESS for EDA analysis of certain dimers. > 1)I dont know the representation of writing 6-311++g** basis,can
> somebody tell how to do it.??
> 2)I ran EDA for a dimer at b3lyp/aug-cc-pvdz and the result shows:
> ALL BASIS SET=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 = =A0 =A0 =A0 =A0 =A0 =A0HARTREE=A0 =A0 =A0 =A0 =A0 KCAL/MOL
>=A0 -------------
>=A0 ELECTROSTATIC ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 ES=3D=A0 =A0= =A0 =A0 =A0-1.236592=A0 =A0 =A0 =A0 =A0
>=A0 -775.97
>=A0 EXCHANGE ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0EX=3D= =A0 =A0 =A0 =A0 =A0-1.002818=A0 =A0 =A0 =A0 =A0
>=A0 -629.28
>=A0 REPULSION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0REP=3D=A0= =A0 =A0 =A0 =A0 5.291539=A0 =A0 =A0 =A0 =A0
>=A0 3320.49
>=A0 POLARIZATION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0 POL=3D=A0 =A0= =A0 =A0 =A0-3.962384=A0 =A0 =A0 =A0 =A0
> -2486.44
>=A0 DFT DISPERSION ENERGY=A0 =A0 =A0 =A0 =A0 =A0 =A0 =A0DISP=3D=A0 =A0 = =A0 =A0 =A0-0.095027=A0 =A0 =A0 =A0 =A0
>=A0 =A0-59.63
>=A0 TOTAL INTERACTION ENERGY HF OR DFT=A0 =A0 =A0E=3D=A0 =A0 =A0 =A0 =A0= -1.005283=A0 =A0 =A0 =A0 =A0
>=A0 -630.83
>
> How can the electrostatic energy be so huge..!!
> AM i doing something wrong in giving input??..The input i gave = is :-
>=A0 $CONTRL SCFTYP=3DRHF RUNTYP=3DEDA DFTTYP=3DB3LYP MAXIT=3D200
>=A0 =A0 =A0MULT=3D1 ISPHER=3D1 COORD=3DUNIQUE NOSYM=3D1 $END
>=A0 $SYSTEM TIMLIM=3D525600=A0 MWORDS=3D50 MEMDDI=3D500 $ENDI
>=A0 $BASIS GBASIS=3DACCD $END
>=A0 $SCF DIRSCF=3D.TRUE. CONV=3D1.0D-04=A0 $END
>=A0 $LMOEDA MATOM(1)=3D15,11 MCHARG(1)=3D0,0
>=A0 =A0 =A0 MMULT(1)=3D1,1 $END
>=A0 $DATA
>
> It will be a great pleasure if somebody can help me out of it.?
> Thanks In Advance.............
> Indu Kaul
>
>
>
> -=3D This is automatically added to each message by the mailing script = =3D-
> To recover the email address of the author of the message, please chang= e
>
>
>
>
>=A0 =A0 =A0 =A0http://www.ccl.net/cgi-bin/ccl/send_ccl_message
>
>
>=A0 =A0 =A0 =A0http://www.ccl.net/cgi-bin/ccl/send_ccl_message
>
> Subscribe/Unsubscribe:
>=A0 =A0 =A0 =A0http://www.ccl.net/chemistry/sub_unsub.shtml
>
> Before posting, check wait time at: http://www.ccl.net
>
> Job: http://www.c= cl.net/jobs
> Conferences: http://server.ccl.net/chemistry/announcements= /conferences/
>
> Search Messages: http://www.ccl.net/chemistry/searchccl/index.shtm= l
>
>
>=A0 =A0 =A0 =A0http://www.ccl.net/spammers.txt
>
> RTFI: http://www.ccl.net/chemistry/aboutccl/instructions/
>
>
>

=A0

___________________________= ___________
Giovanni Finoto Caramori PhD
Departamento de Quimica - CFM-UFSC
Campus Universitario Trindade - C.P. 476
88040-900 Florianopolis/SC
http://www.qm= c.ufsc.br/~caramori
Phone: ++55(48)3721-6844/6845/6846
Fax: ++55(48)3721-6852 Ramal 239
Mobile: ++55(48)9926-8342



----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
--=_5uy8r2lgm34y-- From owner-chemistry@ccl.net Sat Apr 16 11:27:00 2011 From: "Jeff Hammond jeff.science\a/gmail.com" To: CCL Subject: CCL: computational chemistry on accelerators session at SAAHPC 2011 Message-Id: <-44390-110416105835-31984-2wRwpJ7st8wW+FLlbl6Vtg^server.ccl.net> X-Original-From: "Jeff Hammond" Date: Sat, 16 Apr 2011 10:58:32 -0400 Sent to CCL by: "Jeff Hammond" [jeff.science\a/gmail.com] Members of CCL working with GPUs may be interested in submitting to the computational chemistry session of SAAHPC 2011. Please note the paper submission deadline of May 6, 2011. Regards, Jeff Hammond Argonne Leadership Computing Facility jhammond%x%alcf.anl.gov http://www.linkedin.com/in/jeffhammond =================================================================== SAAHPC 2011 call for papers The 2011 Symposium on Application Accelerators in High-Performance Computing will be held July 19-21, 2011, in Knoxville, Tennessee, bringing together developers of hybrid computing architectures and computing accelerators and experts involved with the development of applications on such systems to exchange ideas and discuss issues and emerging challenges. Presentations from technology developers and the academic user community are invited on the following topics: - novel accelerator processors, systems, and architectures - integration of accelerators with high-performance computing systems - programming models for accelerator-based computing - languages and compilers for accelerator-based computing - run-time environments, profiling and debugging tools for accelerator-based computing - scientific and engineering applications that use application accelerators In addition to the general session, submissions are invited for the following domain-specific topics: - Computational chemistry on accelerators (Chair: Jeff Hammond, Argonne Leadership Computing Facility) - Lattice QCD (Chair: Steven Gottlieb, Indiana University, Bloomington) - Weather and climate modeling (Chair: John Michalakes, National Renewable Energy Laboratory) - Bioinformatics (Chair: Jeremy Buhler, Washington University in St. Louis) Visit the EasyChair conference submission system to submit either a short paper (up to 4 pages, for a poster presentation) or a long paper (up to 10 pages, for an oral presentation), by May 6, 2011. All submissions will be reviewed by the Technical Program Committee and the accepted submissions will be presented as either oral presentations or posters. Papers presented at the symposium will be included in the proceedings published by Conference Publishing Services (pending approval). Important dates: Submissions due: May 6, 2011 Presentation acceptance notification: June 6, 2011 Final papers due: June 30, 2011 From owner-chemistry@ccl.net Sat Apr 16 12:02:00 2011 From: "Henry Martinez hmartine~~gmail.com" To: CCL Subject: CCL:G: Do Calculation without nitrogen lone pair Message-Id: <-44391-110416113332-18278-C85t6km/arj/+S/oavk9/g[A]server.ccl.net> X-Original-From: "Henry Martinez" Date: Sat, 16 Apr 2011 11:33:29 -0400 Sent to CCL by: "Henry Martinez" [hmartine=gmail.com] Hi everyone, I do have access to Gaussian 03, and I would like to know if there is any way to perform a calculation "telling" the calculation to not "use" the nitrogen lone pair during the optimization. I want to do this in order to compare the ground state of a molecule, with and without the nitrogen Lone pair. Thanks a lot for the help. From owner-chemistry@ccl.net Sat Apr 16 12:48:01 2011 From: "Esteban Gabriel Vega Hissi egvega/./gmail.com" To: CCL Subject: CCL:G: Do Calculation without nitrogen lone pair Message-Id: <-44392-110416124610-19669-CZNnVtlkZD1Nqvt/rGGAqQ],[server.ccl.net> X-Original-From: Esteban Gabriel Vega Hissi Content-Type: multipart/alternative; boundary=90e6ba3fcdcbbc718e04a10be44a Date: Sat, 16 Apr 2011 18:45:41 +0200 MIME-Version: 1.0 Sent to CCL by: Esteban Gabriel Vega Hissi [egvega.:.gmail.com] --90e6ba3fcdcbbc718e04a10be44a Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Hi, Try NBO implemented in Gaussian. I haven't perform anything like this but I've read it's possible. http://www.chem.wisc.edu/~nbo5/TUT_DEL.HTM Hope this help --=20 Esteban Gabriel Vega Hissi =C1rea de Qu=EDmica F=EDsica Departamento de Qu=EDmica Universidad Nacional de San Luis Argentina --90e6ba3fcdcbbc718e04a10be44a Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
Hi,

Try NBO implemented in Gaussian. I haven't= perform anything like this but I've read it's possible.
http://www.chem.wisc.edu/~= nbo5/TUT_DEL.HTM

Hope this help

--
Esteban Gabriel Veg= a Hissi
=C1rea de Qu=EDmica F=EDsica
Departamento de Qu=EDmica
Uni= versidad Nacional de San Luis
Argentina

--90e6ba3fcdcbbc718e04a10be44a-- From owner-chemistry@ccl.net Sat Apr 16 13:31:00 2011 From: "Rajan Chaudhari rchaudhari-,-mail.usp.edu" To: CCL Subject: CCL: Docking small molecules in the binding site of metalloproteins Message-Id: <-44393-110416131027-14187-Kk3j8A9mNVWgouhxb/SUVg++server.ccl.net> X-Original-From: Rajan Chaudhari Content-Type: multipart/alternative; boundary="_a46d62f5-2e51-45b6-b4fd-911a42f696c1_" Date: Sat, 16 Apr 2011 13:10:16 -0400 MIME-Version: 1.0 Sent to CCL by: Rajan Chaudhari [rchaudhari]-[mail.usp.edu] --_a46d62f5-2e51-45b6-b4fd-911a42f696c1_ Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hi Poonam=2C Read the following article: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2744578/ PLOP is free for academic use. Manual: http://www.jacobsonlab.org/plop_manu= al/plop_manual.htm=20 Autodock: http://autodock.scripps.edu/faqs-help/faq/faqsection_view?section=3DScienti= fic%20Questions How can I set up the protonation state of my histidine sidechains?Histidine= s can be neutral or positively charged. When neutral=2C they can be protona= ted at the delta (HD1) or epsilon (HE2) positions. How can I set these up?U= sing ADT There is a command in ADT to help you decide on the protonation of the Hist= idines=2C but you have to load the commands before you can use it: go to "F= ile > Load Module" and then scroll down=2C click on "repairCommands"=2C and= then click "Load Module" followed by "Dismiss". Now=2C go to "Edit > Hydro= gens > Edit Histidine Hydrogens". If there are any histidines in your molecule=2C a panel will open up listin= g each histidine residue along with a row of radio buttons. You can use the= se to choose whether each histidine should be neutral=2C HD1=3B neutral=2C = HE2=3B or protonated. Using Reduce/Molprobity There is a very nice tool called Reduce (with a web-accessible front end ca= lled Molprobity) that can be used for adding hydrogens and optimising the h= ydrogen-bond network by flipping amido groups in Asn and Gln sidechains=2C = and His imidazole rings by 180=BA. It can also be used for evaluating the q= uality of your protein structure. See: Word=2C et al. (1999) "Asparagine and glutamine: using hydrogen atom contac= ts in the choice of sidechain amide orientation" J. Mol. Biol. 285=2C 1733-= 1745. You need to study your receptor's active site as well as mechanism of actio= n of binding ligand. you will get idea about protonation state of amino aci= ds before the ligand interaction and after ligand interaction. Hope this helps. Best Luck! Regards=2CRajan > From: owner-chemistry,ccl.net To: rchaudhari,mail.usp.edu Subject: CCL: Docking small molecules in the binding site of metalloprotein= s Date: Fri=2C 15 Apr 2011 16:08:01 -0600 Hello I am trying to dock some small molecules with a metalloprotein having manga= nese ions in its binding site.The binding site of my receptor protein comprises mostly His=2CAsp=2Cand Glu residues and it is known that the pres= ence of metal ions can alter the hydration and protonation state of these amino aci= d residues in the binding site.Please suggest me as how should I optimize the metal ion parameters and the protonation states of these amino acid residue= s in the binding site.I am using Surflex-Dock (Tripos Inc.) and Autodock for doc= king studies.Also=2C please suggest me some relevant articles in this regar= d. Thank you =20 Prija Ponnan=20 =20 Research Student Department of Chemistry University of Delhi Delhi-110007=2C India & Canadian Commonwealth Scholarship Program Scholar College of Pharmacy and Nutrition University of Saskatchewan Saskatoon=2CSKN5C9=2CCanada =20 = --_a46d62f5-2e51-45b6-b4fd-911a42f696c1_ Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hi Poonam=2C

Read the following a= rticle:

=
http= ://www.ncbi.nlm.nih.gov/pmc/articles/PMC2744578/

PLOP is free for = academic use. =3B
Manual: http://www.jacobsonlab.org/plop_manual/plop= _manual.htm =3B

Autodock:

http://autodock.scripps.ed= u/faqs-help/faq/faqsection_view?section=3DScientific%20Questions

<= meta charset=3D"utf-8">

How can I set up the protonation state of my histidine sidechains?<= /h3>

Histidines can be neutral or positively charged. When n= eutral=2C they can be protonated at the delta (HD1) or epsilon (HE2) positi= ons. How can I set these up?

Using ADT


There is a command in ADT to help= you decide on the protonation of the Histidines=2C but you have to load th= e commands before you can use it: go to "File >=3B Load Module" and then = scroll down=2C click on "repairCommands"=2C and then click "Load Module" fo= llowed by "Dismiss". Now=2C go to "Edit >=3B Hydrogens >=3B Edit Histid= ine Hydrogens".

If there are any histidines in your molecule=2C a pa= nel will open up listing each histidine residue along with a row of radio b= uttons. You can use these to choose whether each histidine should be neutra= l=2C HD1=3B neutral=2C HE2=3B or protonated.

Using Reduce/Molprobity


There is a very nice tool called&nb= sp=3BReduce&nbs= p=3B(with a web-accessible front end called =3BMolprobity) that can be used for adding hydrogens =3B= and =3Boptimising the hydrogen-bond network by flipping amido groups in= Asn and Gln sidechains=2C and His imidazole rings by 180=BA. It can also b= e used for evaluating the quality of your protein structure. See:

Wo= rd=2C =3Bet al. =3B(1= 999) "Asparagine and glutamine: using hydrogen atom contacts in the choice = of sidechain amide orientation" = =3BJ. Mol. Biol. =3B285=2C 1733-1745.


<= div>You need to= study your receptor's active site as well as mechanism of action of bindin= g ligand. you will get idea about protonation state of amino acids before t= he ligand interaction and after ligand interaction.


Hope th= is helps. =3BBest Luck!

Regards=2C
Rajan<= /font>





From: owner-chemistry,ccl.net
<= font class=3D"Apple-style-span" face=3D"Tahoma" size=3D"2">To: rchaudhari,m= ail.usp.edu

Subject: CCL: Docking small molecules in the binding site of metallo= proteins
Date: Fri=2C 15 Apr 2011 16:08:01 -0600

Hello

I am trying to dock some s= mall molecules with a metalloprotein having manganese ions in its binding s= ite.The binding site of my receptor protein comprises mostly His=2CAsp=2Cand Glu residues and it is known that the pres= ence of metal ions can alter the hydration and protonation state of these amino aci= d residues in the binding site.Please suggest me as how should I optimize the metal ion parameters and the protonation states of these amino acid residue= s in the binding site.I am using Surflex-Dock (Tripos Inc.) and Autodock for docking stu= dies.Also=2C please suggest me some relevant articles in this regard.

Thank = you
&= nbsp=3B





Prija Ponnan=20
 =3B
Research Student<= /div>
Department of Che= mistry
University of Del= hi
Delhi-110007=2C I= ndia
&=3B
Canadian Commonwe= alth Scholarship Program Scholar
College of Pharma= cy and Nutrition
University of Sas= katchewan
Saskatoon=2CSKN5C= 9=2CCanada
 =3B
= --_a46d62f5-2e51-45b6-b4fd-911a42f696c1_-- From owner-chemistry@ccl.net Sat Apr 16 22:55:00 2011 From: "Hyunbok Lee mutebeat^-^gmail.com" To: CCL Subject: CCL: treating unrestrict fragment on ADF Message-Id: <-44394-110416224948-11661-9286R8MfinECTAe7UrbDWQ[]server.ccl.net> X-Original-From: "Hyunbok Lee" Date: Sat, 16 Apr 2011 22:49:38 -0400 Sent to CCL by: "Hyunbok Lee" [mutebeat[-]gmail.com] Dear all, I'm trying to calculate transfer integrals between the identical molecules which have odd number electrons using ADF program. However, I faced the error message "ERROR DETECTED: UNRESTRICTED FRAGMENT". Can't unrestrict fragment be used on ADF? If so, how can I solve this problem? Thanks in advance.