From owner-chemistry@ccl.net Tue Jan 18 04:26:00 2011 From: "sobereva sobjubao|*|yahoo.com.cn" To: CCL Subject: CCL: Looking for a source code for calculating Rys Roots and Weights Message-Id: <-43676-110118042150-2711-h2pJNHhM45N0B78AA3Oj6g%server.ccl.net> X-Original-From: sobereva Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=iso-8859-1 Date: Tue, 18 Jan 2011 17:21:33 +0800 (CST) MIME-Version: 1.0 Sent to CCL by: sobereva [sobjubao|*|yahoo.com.cn] Hi, Lan You can make use of subroutines in ryspol.f of GAMESS(US) to generate Rys roots and weights. Tian Lu --- On Tue, 1/18/11, Youzhao Lan lyzhao(0)zjnu.cn wrote: > From: Youzhao Lan lyzhao(0)zjnu.cn > Subject: CCL: Looking for a source code for calculating Rys Roots and Weights > To: "Lu, Tian " > Date: Tuesday, January 18, 2011, 8:48 AM > > Sent to CCL by: "Youzhao  Lan" [lyzhao|a|zjnu.cn] > Dear All, >    Is there a source code for calculating > the roots and weights of Rys polynomial? > > Thanks for your help. > > Lan > > > > -= This is automatically added to each message by the > mailing script =- > To recover the email address of the author of the message, > please change > the strange characters on the top line to the |*| sign. You > can also> > E-mail to subscribers: CHEMISTRY|*|ccl.net > or use: >      > > E-mail to administrators: CHEMISTRY-REQUEST|*|ccl.net > or use >      >      >      > > > From owner-chemistry@ccl.net Tue Jan 18 06:50:00 2011 From: "Lukman Olawale Olasunkanmi waleolasunkanmi ~ gmail.com" To: CCL Subject: CCL: geometry optimization help Message-Id: <-43677-110118064731-20334-X9knLBMk2/AS8YmuSFjSVQ^server.ccl.net> X-Original-From: Lukman Olawale Olasunkanmi Content-Type: multipart/alternative; boundary=0016e64c1e765a10b9049a1d7654 Date: Tue, 18 Jan 2011 11:46:46 +0000 MIME-Version: 1.0 Sent to CCL by: Lukman Olawale Olasunkanmi [waleolasunkanmi:+:gmail.com] --0016e64c1e765a10b9049a1d7654 Content-Type: text/plain; charset=ISO-8859-1 {:} Victor Rosas & others, Thank you for your earlier reply to my mail on how to build complexes. I discover that optimization with Avogadro gives varying coordinates and energies of the molecule with continuous click on the optimization option. How do I know when to stop optimizing? Because the coordinates I use in MOPAC calculations is the determinant of the accuracy of my results. I will really appreciate your contributions on this. Thank you. -- Olasunkanmi Lukman Olawale +234-08052401564 --0016e64c1e765a10b9049a1d7654 Content-Type: text/html; charset=ISO-8859-1 {:} Victor Rosas & others,
Thank you for your earlier reply to my mail on how to build complexes.
I discover that optimization with Avogadro gives varying coordinates and energies of the molecule with continuous click on the optimization option.
How do I know when to stop optimizing? Because the coordinates I use in MOPAC calculations is the determinant of the accuracy of my results.
I will really appreciate your contributions on this.
Thank you.

--
Olasunkanmi Lukman Olawale
+234-08052401564
--0016e64c1e765a10b9049a1d7654-- From owner-chemistry@ccl.net Tue Jan 18 07:25:00 2011 From: "P.D.Jarowski^surrey.ac.uk" To: CCL Subject: CCL: TD Errors Message-Id: <-43678-110118071511-12632-S5jOet1eALF4GYQDvYxv2w,+,server.ccl.net> X-Original-From: Content-Language: en Content-Type: multipart/alternative; boundary="_000_C95B34BFC49pj0013surreyacuk_" Date: Tue, 18 Jan 2011 12:14:55 +0000 MIME-Version: 1.0 Sent to CCL by: [P.D.Jarowski:-:surrey.ac.uk] --_000_C95B34BFC49pj0013surreyacuk_ Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Hello, Anyone have a reference for the error bars in a TD calculation in the lamba= max and oscillator strengths? Thanks in Advance. Peter --_000_C95B34BFC49pj0013surreyacuk_ Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable TD Errors Hello,

Anyone have a reference for the error bars in a TD calculation in the lamba= max and oscillator strengths?

Thanks in Advance.

Peter
--_000_C95B34BFC49pj0013surreyacuk_-- From owner-chemistry@ccl.net Tue Jan 18 08:20:00 2011 From: "Andreas Klamt klamt(!)cosmologic.de" To: CCL Subject: CCL: Pioglitazone Tautomers Message-Id: <-43679-110118014940-12219-0RTv1Dgtf/qJeL5CIAe0hA|,|server.ccl.net> X-Original-From: Andreas Klamt Content-Transfer-Encoding: 7bit Content-Type: text/html; charset=ISO-8859-1 Date: Tue, 18 Jan 2011 07:49:24 +0100 MIME-Version: 1.0 Sent to CCL by: Andreas Klamt [klamt]![cosmologic.de] Nancy,

although I find Gabriele's reply a bit to harsh, I also was a bit frustrated when I read your last entry, which essentially repeats the initial question after you got a number of replies. You now should draw your own conclusions on this.

What you can learns from the diverse answers is that tautomerism is still a complicated and quite unsolved question. This was also shown in the special issue of J Comput Aided Mol Des. 2010 Jun,24 on tautomerism edited by Yvonne Martin and by the outcome of the SAMPLE2 blindtest (see Journal of Computer-Aided Molecular Design, Volume 24, Number 4).

For your case of Pioglitazone I would conclude that the quantum chemical results for the tautomer energy differences are so large, that a rather unique answer can be given: The diketonic form (neutral and deprotonated) are absolutely predominant. This is a lucky case: There are many cases of tatuomerism where the DFT energy differences are in the range of a few kcal/mol, i.e. within the error bars of DFT,  and then no such clear answer can be given.

I would go so far to conclude that the huge energy differences of ~ 12 kcal/mol to the other forms also makes the other forms rather unlikely in any protein docked form. I cannot see how the local environment of the protein can be so much more favoring the enolic forms to compensate such penalty.

If the pKa of 6.3 reported by Gabriele is an experimental one, I have no problem to admit that most likely his 85% deprotonated is the better answer than my answer based on DFT/COSMO-RS prediction of pKa=6.55. The difference of 0.25 pKa units is far within the error bars (~ 0.5) of our pKa prediction.

Good luck in going on from here. Perhaps you can send me your final paper, or thesis or what ever kind of document you will produce on this case.

Andreas



Am 17.01.2011 22:14, schrieb Gabriele Cruciani gabri- -chemiome.chm.unipg.it:

Sent to CCL by: Gabriele Cruciani [gabri|a|chemiome.chm.unipg.it]
Nancy,
you get the same answer from different people, with different background, using different methods, and you keep asking the same question.

I think CCL list should be used 'with caution' especially after you get all the answers.

What you summarized is not correct. The two forms are not tautomeric forms. They refer to ONE tautomeric form that is in equilibrium with two species, one protonated and one neutral. The percentages you reported refer to the solvent pH and to the specie pKa. Since the pKa of the pioglitazone is 6.3, the relative abundance at pH 7.4 is
85% for the charged and 15% for the neutral one.

As told in previous mails, this has nothing to do with the docking into a protein, where the local environment, dielectric, 3D interactions may change the 'form' of the interacting ligand. In this specific case, you cannot say that the charge form may interact better because is more abundant, since the equilibrium between the two forms can be modified by the protein.

Precise pKa calculations are very relevant for pharmacokinetic and metabolism studies, but for docking to receptors is better to produce the 'relevant' structures. Do not use too simple methods (like those you reported) since the problem is very complex and require good approaches. Often using the brain is better than using free available software.

Gabriele Cruciani



Thus, the two major tautomers of pioglitazone are as shown in the attached
figure?

Also, is there a published article, or other evidence that would confirm
that this is indeed the correct tautomeric distribution at pH 7?

Thanks you,
Nancy



On Mon, Jan 17, 2011 at 6:17 AM, Andreas Klamt klamt#cosmologic.de<
owner-chemistry/a\ccl.net>  wrote:

Sent to CCL by: Andreas Klamt [klamt-x-cosmologic.de]
I was just pointed by Leif Nørskov-Lauritsen to a stupid error which I made
in my first entry on this topic. I wrote that at pH=7 25% would be
deprotonated, 75% neutral. That was obviously wrong. It is the other way
round: 72% is ionic, 25% neutral. The total amount of zwitterions is then
~3%.

Sorry for my mistake. But the basic message about the most stable neutral
form stays valid.
The most stable anionic form is the di-keto anion.

Andreas

Am 16.01.2011 19:29, schrieb Nancy nancy5villa**gmail.com:

Using quantum chemical calculations, you found that the diketone form is
dominant; does this mean that the enol form predicted by MarvinSketch v3.5.4
in the published article is wrong?

"Structure-based design of a thiazolidinedione which targets the
mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 Feb
1;20(3):819-23

Thanks,
Nancy


On Thu, Jan 13, 2011 at 10:31 PM, Andreas Klamt klamt(a)cosmologic.de<
owner-chemistry(~)ccl.net<owner-chemistry%28%7E%29ccl.net>>  wrote:

  Sent to CCL by: Andreas Klamt [klamt,cosmologic.de]
  Hi Nancy,

as I wrote in the previous message, we find by quantum chemical
calculations within the COSMO-RS solvation that the diketone form is more
stable by abot 12 kcal/mol, that means it absolutely dominates!

Andreas

Am 13.01.2011 22:25, schrieb Nancy nancy5villa~!~gmail.com:

  Hi All,

I happened to find a published article where molecular docking simulations
of TZDs against a novel protein is detailed:

"Structure-based design of a thiazolidinedione which targets the
mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 Feb
1;20(3):819-23

In this paper, the authors, using MarvinSketch v3.5.4, determined that the
TZDs would exist predominantly in the enol form at pH 7.4, as opposed to the
diketone form (see attached figure).  It appears that the article's tautomer
prediction was based only on the pH prior to docking, and had nothing to do
with ligand-protein interactions.

Does anyone know if the diketone or enol, or perhaps a different tautomer,
would be predominant at pH 7.4?

Thanks in advance.
Nancy



   On Wed, Jan 12, 2011 at 9:53 AM, Andreas Klamt klamt-*-cosmologic.de<
owner-chemistry(_)ccl.net<owner-chemistry%28_%29ccl.net>>  wrote:


Sent to CCL by: Andreas Klamt [klamt[*]cosmologic.de]
Hi All,

let me put my 5 cents into the game:
We have done quantum calculations for the heterocyclic ring, with just a
methyl group, since electronically the remainder of the compound should have
almost no influence on the heterocycle.
DFT/COSMO-RS calculations prove the neutral species
SMILES:CC1C(=O)[NH]C(=O)S1 of the Pioglitazone heterocycle to be favored
by 11.3 kcal/mol compared to all other neutral species of the heterocycle,
which means that there is no other relvant neutral form.

We find with COSMOtherm that the pka of the heterocycle is 6.55, i.e. at
pH=7 it should be deprotonated with about 25%.
Since the pka of the conj. acid of the pyridine ring according to
COSMOtherm (we calculated diethylpyridine) is ~5.6 there should be ~10% in
the zwitterionic form.

It will be hard to decide this finally, unless someone measures it, but
we are quite confident to be qualitatively right here, because the energy
differences for the tautomers should not have more than 5 kcal/mol erro, as
we learned in in the SAMPL2 tautomer contest:
see Andreas Klamt and Michael Diedenhofen, Some conclusions regarding the
predictions of tautomeric equilibria in solution based on the SAMPL2
challenge, http://www.springerlink.com/content/l577667th758n6h1/

Andreas







Hi All,

Using MarvinSketch v5.3.3 (ChemAxon software), the predicted pKa value
of the acidic hydrogen on the thiazolidinedione ring of pioglitazone is
4.57 (see attached "Figure_1.gif"). Therefore, the predominant species
at pH 7.0 predicted by MarvinSketch is the one depicted in
"Figure_2.gif". The different tautomeric forms predicted by MarvinSketch
are shown in "Figure_3.gif".

Can anyone explain where these predictions come from, and if they are
correct?

Thanks,
Nancy


On Mon, Jan 10, 2011 at 4:16 AM, Gabriele Cruciani
  gabri(0)chemiome.chm.unipg.it<http://chemiome.chm.unipg.it>

<owner-chemistry[a]ccl.net<mailto:owner-chemistry[a]ccl.net>>  wrote:


    Sent to CCL by: Gabriele Cruciani [gabri##chemiome.chm.unipg.it
     <http://chemiome.chm.unipg.it>]

    Nancy,
    the first form you reported is the most stable in water at pH 7.0.
    However, the fact that one form is more stable than another in water
    does not help you to understand which form will be more 'relevant'
    for docking. In protein the tautomeric equilibria may produce and
    stabilize different forms according to the complementary site. There
    are examples of tautomeric form in protein not stable in water,
    where the energy difference is more than 5 Kcal/mol.

     MoKa software (www.moldiscovery.com<http://www.moldiscovery.com>)

    is fast and accurate to produce tautomer stable in water, but it can
    produce also all the (plausible) tautomeric forms.

    Then, a possibility is to dock them into your protein, to see if
    docking methods may differentiate their binding.

    Gabriele Cruciani




        Hi All,

        I am performing molecular docking and molecular dynamics
        simulations of the
        thiazolidinedione pioglitazone binding to the PPAR-gamma
        receptor protein
        (PDB ID: 1ZGY). The thiazolidinedione ring can exist in numerous
        different
        tautomeric states; I have attached a figure depicting several of
        them.
        Which tautomer would be dominant at the physiological pH of ~7.0?

        Also, are there any software programs that can predict which
        tautomer would
        be correct?

        Thanks in advance,
        Nancy>


     E-mail to subscribers: CHEMISTRY[a]ccl.net
    <mailto:CHEMISTRY[a]ccl.net>  or use:>
    E-mail to administrators: CHEMISTRY-REQUEST[a]ccl.net
    <mailto:CHEMISTRY-REQUEST[a]ccl.net>  or useConferences:

    http://server.ccl.net/chemistry/announcements/conferences/>




--
Dr. Jens Reinisch
Chemist / Customer Support
COSMOlogic GmbH&  Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone   +49-2171-363664
mobile  +49-163-7337310
fax     +49-2171-731689
e-mail  reinisch[*]cosmologic.de
web     www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt






  E-mail to subscribers: CHEMISTRY(_)ccl.net<CHEMISTRY%28_%29ccl.net>  or
use:>>>
E-mail to administrators: CHEMISTRY-REQUEST(_)ccl.net<CHEMISTRY-REQUEST%28_%29ccl.net>or useConferences:
http://server.ccl.net/chemistry/announcements/conferences/>>>




--
PD. Dr. Andreas Klamt
CEO / Geschäftsführer
COSMOlogic GmbH&  Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone      +49-2171-731681
fax        +49-2171-731689
e-mail     klamt(_)cosmologic.de<klamt%28_%29cosmologic.de>

web        www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt

CHEMISTRY(~)ccl.net<CHEMISTRY%28%7E%29ccl.net>  or use:E-mail to administrators:
CHEMISTRY-REQUEST(~)ccl.net<CHEMISTRY-REQUEST%28%7E%29ccl.net>  or usehttp://www.ccl.net/chemistry/sub_unsub.shtml<http://www.ccl.net/cgi-bin/ccl/send_ccl_message>Before posting, check wait time at:
http://www.ccl.netConferences:
http://server.ccl.net/chemistry/announcements/conferences/ Search
Messages: http://www.ccl.net/chemistry/searchccl/index.shtmlRTFI:
http://www.ccl.net/chemistry/aboutccl/instructions/




--
PD. Dr. Andreas Klamt
CEO / Geschäftsführer
COSMOlogic GmbH&  Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone      +49-2171-731681
fax        +49-2171-731689
e-mail     klamt:cosmologic.de
web        www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas KlamtTo
recover the email address of the author of the message, please change the
strange characters on the top line to the /a\ sign. You can also look up the
X-Original-From: line in the mail header. E-mail to subscribers:
CHEMISTRY/a\ccl.net or use:E-mail to administrators:
CHEMISTRY-REQUEST/a\ccl.net or useBefore posting, check wait
time at: http://www.ccl.netConferences:
http://server.ccl.net/chemistry/announcements/conferences/ Search
Messages: http://www.ccl.net/chemistry/searchccl/index.shtml If your mail
bounces from CCL with 5.7.1 error, check:RTFI:
http://www.ccl.net/chemistry/aboutccl/instructions/





E-mail to subscribers: CHEMISTRY]*[ccl.net or use:
     http://www.ccl.net/cgi-bin/ccl/send_ccl_message

E-mail to administrators: CHEMISTRY-REQUEST]*[ccl.net or use
     http://www.ccl.net/cgi-bin/ccl/send_ccl_message
http://www.ccl.net/chemistry/sub_unsub.shtml

Before posting, check wait time at: http://www.ccl.net

Job: http://www.ccl.net/jobs Conferences: http://server.ccl.net/chemistry/announcements/conferences/

Search Messages: http://www.ccl.net/chemistry/searchccl/index.shtml
     http://www.ccl.net/spammers.txt

RTFI: http://www.ccl.net/chemistry/aboutccl/instructions/





-- 
PD. Dr. Andreas Klamt
CEO / Geschäftsführer
COSMOlogic GmbH & Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone  	+49-2171-731681
fax    	+49-2171-731689
e-mail 	klamt]*[cosmologic.de
web    	www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt
From owner-chemistry@ccl.net Tue Jan 18 08:55:00 2011 From: "sumit kumar sumitkumar!^!iiserpune.ac.in" To: CCL Subject: CCL:G: About negative frequency for stable conformers in NWCHEM Message-Id: <-43680-110118021527-18801-w0OTmlK8tzb5w6a91DEmag[#]server.ccl.net> X-Original-From: "sumit kumar" Date: Tue, 18 Jan 2011 02:14:46 -0500 Sent to CCL by: "sumit kumar" [sumitkumar---iiserpune.ac.in] Dear Sir, We have performed geometry optimization of mixed dimer of indole and pyridine at m05-2x/6-31+g* level using Gaussian 03 and NWCHEM softwares. In case of gaussian calculation we have obtained three sable conformers,where all of them have positive frequency only. But when same three conformers were run using NWCHEM, we have found that two of them are giving few negative frequency. Also NH stretch frequency of one of the conformers is very much different in NWCHEM (3730cm-1) than that obtained in gaussian (3636cm-1) We don't understand how same theory and same basis set using two softwares are giving this different result. I have attached here both gaussian and NWCHEM files. These three conformers when were run also at higher level of theory ( MP2 ) as well as higher basis set ( cc-pVDZ,aug-cc-pVDZ ), using gaussian, don't give any negative frequency. But when they are run in NWCHEM gives negative frequency. Using both Gaussian and NWCHEM, DFT level calculations are done using Grid=Ultrafine command. I will appreciate help us for solving this problem in NWCHEM . sincerely Yours -- -- Sumit Kumar Research Scholar (laser spectroscopy) Indian Institutes of Science Education and Research, (iiser,pune) mobile no:-9579450484 From owner-chemistry@ccl.net Tue Jan 18 09:30:00 2011 From: "Olasunkanmi Lukman Olawale walecomuk ~ yahoo.co.uk" To: CCL Subject: CCL: complexes Message-Id: <-43681-110118045051-8397-rIOMkiC5z4nv1eQeEmh0GQ ~ server.ccl.net> X-Original-From: Olasunkanmi Lukman Olawale Content-Type: multipart/alternative; boundary="0-747588781-1295344236=:94901" Date: Tue, 18 Jan 2011 09:50:36 +0000 (GMT) MIME-Version: 1.0 Sent to CCL by: Olasunkanmi Lukman Olawale [walecomuk-.-yahoo.co.uk] --0-747588781-1295344236=:94901 Content-Type: text/plain; charset=utf-8 Content-Transfer-Encoding: quoted-printable ,+, Victor Rosas & others,=0AThank you for your earlier reply to my mail on h= ow to build complexes.=0AI discover that optimization with Avogadro gives = varying coordinates and =0Aenergies of the molecule with continuous click = on the optimization option.=0AHow do I know when to stop optimizing? Becau= se the coordinates I use in MOPAC =0Acalculations is the determinant of th= e accuracy of my results.=0AI will really appreciate your contributions on = this.=0AThank you.=0A=0A Olasunkanmi Lukman Olawale=0A_____________________= ___________=0ACurrent Address:=0ADepartment of Chemistry,=0AObafemi Awolowo= University,=0AIle-Ife, Osun State.=0ANigeria.=0A=0A+234-0-80-52401564 Or += 234-0-80-67161091=0A________________________________=0A =0A=0A=0A=0A=0A____= ____________________________=0AFrom: Victor Rosas Garcia rosas.victor##gmai= l.com =0ATo: "Olasunkanmi, Olawale Lukman = " =0ASent: Thu, 13 January, 2011 15:38:39=0ASubject:= CCL: complexes=0A=0AWell, if you are able to sketch the structure in a pie= ce of paper, you should be =0Aable to build the molecule in Avogadro. Buil= ding mode is the default when =0AAvogadro starts. Personally I tick off th= e option of adjusting hydrogens, as I =0Afind it annoying most of the time.= Of course whatever you sketch will have lots =0Aof strained bonds and angl= es so you will have to "clean" your structure by =0Aenergy minimization (un= der the "Extensions" menu). You may have to resort to a =0Adifferent force= field than the default, because the Fe(II) bonding may not be =0Awell repr= esented by the MMFF94 forcefield. Probably UFF will work, or some =0Asemie= mpirical method (if you have mopac installed). That will give you a =0Area= sonable starting structure for higher level treatment.=0A=0AAlternatively, = if you can get the XYZ coordinates from the crystal structure, =0AAvogadro = will read those and you will get a reasonable structure in one step.=0A=0AH= ope this helps=0A=0AVictor=0A=0A=0A2011/1/13 Olawale Lukman Olasunkanmi wal= eolasunkanmi:_:gmail.com =0A=0A=0A=0A>Sent to CCL = by: "Olawale Lukman Olasunkanmi" [waleolasunkanmi^_^gmail.com]=0A>Thank you= all for your earlier contributions. I am getting there.=0A>I have got Avog= adro installed and have started using it.=0A>Please, can anyone put me thro= ugh how to build complexes like =0A>iron(II)phenanthroline complex with the= package? There is no procedure for such =0A>molecules in the manual.=0A>Th= ank you.=0A>=0A>=0A>=0A>-=3D This is automatically added to each message by= the mailing script =3D-=0A>=0A>E-mail to subscribers: CHEMISTRY*ccl.net or= use:=0A=0A>=0A>E-mail= to administrators: CHEMISTRY-REQUEST*ccl.net or use=0A> http://www.ccl= .net/cgi-bin/ccl/send_ccl_message=0A>=0A>Subscribe/Unsubscribe:=0A> htt= p://www.ccl.net/chemistry/sub_unsub.shtml=0A>=0A>Before posting, check wait= time at: http://www.ccl.net=0A>=0A>Job: http://www.ccl.net/jobs=0A>Confere= nces: http://server.ccl.net/chemistry/announcements/conferences/=0A>=0A>Sea= rch Messages: http://www.ccl.net/chemistry/searchccl/index.shtml=0A>=0A> ==0A>=0A>RTFI: http://www.ccl.net/chemistry= /aboutccl/instructions/=0A>=0A>=0A>=0A=0A=0A=0A --0-747588781-1295344236=:94901 Content-Type: text/html; charset=utf-8 Content-Transfer-Encoding: quoted-printable
,+, Victor Rosas & others,
Thank you for your earlier rep= ly to my mail on how to build complexes.
I=0A discover that optimization= with Avogadro gives varying coordinates and =0Aenergies of the molecule wi= th continuous click on the optimization =0Aoption.
=0AHow do I know when= to stop optimizing? Because the coordinates I use in =0AMOPAC calculations= is the determinant of the accuracy of my results.
I will really appreci= ate your contributions on this.
Thank you.
 
Olasunkanmi Lukman Olawale=

Current Address:
Department of Chemistry,
Obafemi Awolowo University,
Ile-Ife, Osun State.
Nigeria.
+234-0-80-52401564= Or +234-0-80-67161091

<= /div>


From: Victor = Rosas Garcia rosas.victor##gmail.com <owner-chemistry,+,ccl.net>
= To: "Olasunkanmi, Olawale Luk= man " <walecomuk,+,yahoo.co.uk>
Sent: Thu, 13 January, 2011 15:38:39
Subject: CCL: complexes

Well, if you are = able to sketch the structure in a piece of paper, you should be able to build the molecule in= Avogadro.  Building mode is the default when Avogadro starts.  P= ersonally I tick off the option of adjusting hydrogens, as I find it annoyi= ng most of the time. Of course whatever you sketch will have lots of strain= ed bonds and angles so you will have to "clean" your structure by energy mi= nimization (under the "Extensions" menu).  You may have to resort to a= different force field than the default, because the Fe(II) bonding may not= be well represented by the MMFF94 forcefield.  Probably UFF will work= , or some semiempirical method (if you have mopac installed).  That wi= ll give you a reasonable starting structure for higher level treatment.
= =0A
Alternatively, if you can get the XYZ coordinates from the crystal s= tructure, Avogadro will read those and you will get a reasonable structure = in one step.

Hope this helps

Victor

=0A2011/1/13 Olawale Lukman Olasunkanmi waleolasunkanmi:_:gmail.com <owner-chem= istry*ccl.net>
=0A
=0ASent to CCL by: "Olawale Lukman Olasunkanmi" [wale= olasunkanmi^_^gmail.com]
=0AThank you all for your earlier contributions. I am = getting there.
=0AI have got Avogadro installed and have started using i= t.
=0APlease, can anyone put me through how to build complexes like iron= (II)phenanthroline complex with the package? There is no procedure for such= molecules in the manual.
=0AThank you.
=0A
=0A
=0A
=0A-=3D = This is automatically added to each message by the mailing script =3D-
= =0A
=0AE-mail to subscribers: CHEMISTR= Y*ccl.net or use:
=0A      http:/= /www.ccl.net/cgi-bin/ccl/send_ccl_message
=0A
=0AE-mail to admini= strators: CHEMISTRY-REQUE= ST*ccl.net or use
=0A      http:/= /www.ccl.net/cgi-bin/ccl/send_ccl_message
=0A
=0ASubscribe/Unsubs= cribe:
=0A      http://www.ccl.net/chemi= stry/sub_unsub.shtml
=0A
=0ABefore posting, check wait time at: <= a rel=3D"nofollow" target=3D"_blank" href=3D"http://www.ccl.net">http://www= .ccl.net
=0A
=0AJob: http://www.ccl.net/jobs
=0AConferences:= http://server.ccl.net/chemistry/announcemen= ts/conferences/
=0A
=0ASearch Messages: ht= tp://www.ccl.net/chemistry/searchccl/index.shtml
=0A
=0A   &= nbsp;  http://www.ccl.net/spammers.txt
=0A
=0ARTFI: http://www.ccl.net/chemistry/aboutccl/instructions/=
=0A
=0A
=0A

=0A
=0A

=0A=0A=0A=0A <= /body> --0-747588781-1295344236=:94901-- From owner-chemistry@ccl.net Tue Jan 18 10:04:00 2011 From: "Eduard Matito ematito()gmail.com" To: CCL Subject: CCL:G: Energy decomposition programs, such as SAPT Message-Id: <-43682-110118083717-8716-Ot1ZYHIXyWKeIkAhrqneBg!A!server.ccl.net> X-Original-From: Eduard Matito Content-Transfer-Encoding: 7bit Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes Date: Tue, 18 Jan 2011 14:36:58 +0100 Mime-Version: 1.0 (Apple Message framework v936) Sent to CCL by: Eduard Matito [ematito.:.gmail.com] Istvan Mayer and Pedro Salvador have some programs to perform an energy decomposition analysis. They are for free. Most of them use gaussian files (fchk or similar). http://occam.chemres.hu/programs/ Eduard On Jan 17, 2011, at 11:15 PM, Victor Rosas Garcia rosas.victor(!)gmail.com wrote: > There is Natural Energy Decomposition Analysis, in NBO 5.0. Not > free but very cheap. > > Victor > > 2011/1/17 Tapas Kar tapas.kar|*|usu.edu > > Sent to CCL by: "Tapas Kar" [tapas.kar[#]usu.edu] > Hello, > > Is MOLPRO the standard QM program that runs SAPT and SAPT-DFT, or > are there others? > > Also any other Energy decomposition method, besides SAPT (Symmetry- > adapted perturbation theory) and Morokuma Energy decomposition > method (implemented in GAMESS-US), are developed? > > Thanks > Tapas > > > > -= This is automatically added to each message by the mailing script > =- > > > > E-mail to subscribers: CHEMISTRY,ccl.net or use:> > E-mail to administrators: CHEMISTRY-REQUEST,ccl.net or use> Conferences: http://server.ccl.net/chemistry/announcements/ > conferences/> > > -- Eduard Matito Postdoctoral researcher Institute of Physics University of Szczecin Wielkopolska 15 70451 Szczecin (Poland) Phone: +48914441436 http://ematito.webs.com From owner-chemistry@ccl.net Tue Jan 18 10:40:00 2011 From: "Victor Rosas Garcia rosas.victor**gmail.com" To: CCL Subject: CCL: geometry optimization help Message-Id: <-43683-110118090436-8982-ATFRtM6U9SIKoNCZnEShYQ .. server.ccl.net> X-Original-From: Victor Rosas Garcia Content-Type: multipart/alternative; boundary=0015174c43e215695c049a1f604b Date: Tue, 18 Jan 2011 08:04:13 -0600 MIME-Version: 1.0 Sent to CCL by: Victor Rosas Garcia [rosas.victor/a\gmail.com] --0015174c43e215695c049a1f604b Content-Type: text/plain; charset=ISO-8859-1 Yes. Avogadro does minimization mostly as a "clean up" of the structure. It yields chemically reasonable structures but it is not really a thorough optimization. The good thing is that, once you feed those coordinates to MOPAC, the optimizer in MOPAC will look for an energetic stationary point. That is, if you start from similar but not identical geometries, MOPAC will still find the same local stationary point. Hope this helps Victor 2011/1/18 Lukman Olawale Olasunkanmi waleolasunkanmi ~ gmail.com < owner-chemistry-.-ccl.net> > (0) Victor Rosas & others, > Thank you for your earlier reply to my mail on how to build complexes. > I discover that optimization with Avogadro gives varying coordinates and > energies of the molecule with continuous click on the optimization option. > How do I know when to stop optimizing? Because the coordinates I use in > MOPAC calculations is the determinant of the accuracy of my results. > I will really appreciate your contributions on this. > Thank you. > > -- > Olasunkanmi Lukman Olawale > +234-08052401564 > --0015174c43e215695c049a1f604b Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Yes.=A0 Avogadro does minimization mostly as a "clean up" of the = structure.=A0 It yields chemically reasonable structures but it is not real= ly a thorough optimization.=A0 The good thing is that, once you feed those = coordinates to MOPAC, the optimizer in MOPAC will look for an energetic sta= tionary point.=A0 That is, if you start from similar but not identical geom= etries, MOPAC will still find the same local stationary point.

Hope this helps

Victor

2011/1/= 18 Lukman Olawale Olasunkanmi waleolasunkanmi ~ gmail.com <owner-chemistry-.-ccl.net>
(0) Victor Rosas = & others,
Thank you for your earlier reply to my mail on how to buil= d complexes.
I discover that optimization with Avogadro gives varying coordinates and=20 energies of the molecule with continuous click on the optimization=20 option.
How do I know when to stop optimizing? Because the coordinates I use in=20 MOPAC calculations is the determinant of the accuracy of my results.
I w= ill really appreciate your contributions on this.
Thank you.

--
Olasunkanmi Lukman Olawale
+234= -08052401564

--0015174c43e215695c049a1f604b-- From owner-chemistry@ccl.net Tue Jan 18 11:26:00 2011 From: "Ol Ga eurisco1 _ pochta.ru" To: CCL Subject: CCL:G: About negative frequency for stable conformers in NWCHEM Message-Id: <-43684-110118112532-18268-EjjkVOMaJJhYl8vZB4KiOw_-_server.ccl.net> X-Original-From: "Ol Ga" Date: Tue, 18 Jan 2011 11:25:30 -0500 Sent to CCL by: "Ol Ga" [eurisco1_._pochta.ru] Dear Sumit Kumar, I would suggest to optimize conformers by means of NWChem with convergence criteria of stationary point search similar to Tight in Gaussian and calculate frequencies again. I hope it will solve the problem. How large is difference of geometries for conformers optimized by means of Gaussian and NWChem? Sincerely, Ol Ga ----- ----- > From: sumit kumar sumitkumar!^!iiserpune.ac.in Sent: Tuesday, January 18, 2011 10:14 AM To: Ga, Ol Subject: CCL:G: About negative frequency for stable conformers in NWCHEM Sent to CCL by: "sumit kumar" [sumitkumar---iiserpune.ac.in] Dear Sir, We have performed geometry optimization of mixed dimer of indole and pyridine at m05-2x/6-31+g* level using Gaussian 03 and NWCHEM softwares. In case of gaussian calculation we have obtained three sable conformers,where all of them have positive frequency only. But when same three conformers were run using NWCHEM, we have found that two of them are giving few negative frequency. Also NH stretch frequency of one of the conformers is very much different in NWCHEM (3730cm-1) than that obtained in gaussian (3636cm-1) We don't understand how same theory and same basis set using two softwares are giving this different result. I have attached here both gaussian and NWCHEM files. These three conformers when were run also at higher level of theory ( MP2 ) as well as higher basis set ( cc-pVDZ,aug-cc-pVDZ ), using gaussian, don't give any negative frequency. But when they are run in NWCHEM gives negative frequency. Using both Gaussian and NWCHEM, DFT level calculations are done using Grid=Ultrafine command. I will appreciate help us for solving this problem in NWCHEM . sincerely Yours -- -- Sumit Kumar Research Scholar (laser spectroscopy) Indian Institutes of Science Education and Research, (iiser,pune) mobile no:-9579450484 -= From owner-chemistry@ccl.net Tue Jan 18 12:57:00 2011 From: "Pietro Ottati pietro.ottati-.-gmail.com" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43685-110118095736-29056-7qYS4o1cQqV8XrgTJXNjBA===server.ccl.net> X-Original-From: Pietro Ottati Content-Type: multipart/alternative; boundary=0016e64dda762e59af049a201ebe Date: Tue, 18 Jan 2011 15:57:22 +0100 MIME-Version: 1.0 Sent to CCL by: Pietro Ottati [pietro.ottati|a|gmail.com] --0016e64dda762e59af049a201ebe Content-Type: text/plain; charset=ISO-8859-1 Hi all, I'm a Junior Researcher in these days I'm starting my first time work in Computational Chemistry. My job will base on QM/MM the software that I have to use is Qsite Schrodinger Suite. Unfortunately there aren't a lot of guide, tutorial etc... about this software... Can you help me to learn? Do you have some links, books, articles, wiky etc? Anybody use it? Please take in account that I'm newbie and this is my first time. Thanks in advance Pietro --0016e64dda762e59af049a201ebe Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable
Hi all, I'm a Junior Researcher in these=
 days I'm starting my first time work=20
in Computational Chemistry. My job will base on QM/MM the software that I h=
ave=20
to use is Qsite Schrodinger Suite. Unfortunately there aren't a lot of =
guide,=20
tutorial etc... about this software... Can you help me to learn? Do you hav=
e=20
some links, books, articles, wiky etc? Anybody use it?

Please take in account that I'm newbie and this is my first time.

Thanks in advance=20

Pietro
--0016e64dda762e59af049a201ebe-- From owner-chemistry@ccl.net Tue Jan 18 13:32:00 2011 From: "Markus Reiher stc2011(0)ethz.ch" To: CCL Subject: CCL: Conference Announcement: Designing Molecular Functionality: Challenges Message-Id: <-43686-110118075718-29479-rzScLnB16JK1ng4xaeSNwg!A!server.ccl.net> X-Original-From: "Markus Reiher" Date: Tue, 18 Jan 2011 07:57:16 -0500 Sent to CCL by: "Markus Reiher" [stc2011|,|ethz.ch] Dear colleague, we are happy to announce the international conference "47th Symposium on Theoretical Chemistry - Designing Molecular Functionality: Challenges for Theoretical Approaches" which will be held in Sursee, Switzerland from August 21st to August 25th, 2011. The following invited speakers will cover the conference topic: A. Aspuru-Guzik ( Cambridge, MA) J. Autschbach (Buffalo, NY) P.W. Ayers (Hamilton) D.N. Beratan (Durham, NC) D. Crawford (Blacksburg, VA) A. Dreuw (Frankfurt) G. Frenking (Marburg) P. Geerlings (Brussels) L. De Gioia (Milano-Bicocca) S. Grimme (Mnster) W.F. van Gunsteren (Zurich) P. Hnenberger (Zurich) B. Kuhlman (Chapel Hill, NC) O.A. von Lilienfeld (Albuquerque, NM) R. Lindh (Uppsala) M. Parrinello (Lugano) A. Rizzo (Pisa) K. Ruud (Troms) T. Saue (Toulouse) F. Schoenebeck (Zurich) S. Shaik (Jerusalem) G.C. Solomon (Copenhagen) H. Su (Singapore) W. Thiel (Mlheim) L. Visscher (Amsterdam) H.J. Werner (Stuttgart) Experimental keynote lecture: F. Diederich (Zurich) Contributed Talks and an expected number of about 250 posters will cover all topics of Theoretical Chemistry. On August 23rd, the Hellmann Prize will be awarded to a young theoretical chemist. Young researchers are especially invited to present their contributions. Sursee is located in the heart of Switzerland close to the lake of Sempach and near Lucerne. The location can be reached via Zurich Airport and EuroAirport Basel Mulhouse Freiburg, by train or by car. The conference venue is at Campus Sursee which is about a 15 minutes walk from Sursee. A bus service is available from Sursee train station to Campus Sursee. For further information please visit the conference website: http://www.stc2011.ethz.ch/ Early bird registration is open until March 15th, 2011. http://www.stc2011.ethz.ch/registration Please note that the number of participants is restricted to 250. Early registration is therefore recommended. We are looking forward to welcoming you among many scientists working in the field of theoretical and computer-assisted chemistry, presenting and discussing their latest research results. Best regards, Markus Reiher ETH Zurich Laboratory for Physical Chemistry Wolfgang-Pauli-Strasse 10 CH-8093 Zurich/Switzerland email: stc2011*ethz.ch http://www.stc2011.ethz.ch P.S.: The 2nd circular will only be distributed to those who have registered for the conference. From owner-chemistry@ccl.net Tue Jan 18 14:43:00 2011 From: "Jans Alzate-Morales jalzate()utalca.cl" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43687-110118140005-22248-hIpLNh51mGC+FRYbN9Gk/g-x-server.ccl.net> X-Original-From: Jans Alzate-Morales Content-Type: multipart/alternative; boundary=Apple-Mail-23--757596058 Date: Tue, 18 Jan 2011 15:59:43 -0300 Mime-Version: 1.0 (Apple Message framework v1082) Sent to CCL by: Jans Alzate-Morales [jalzate++utalca.cl] --Apple-Mail-23--757596058 Content-Transfer-Encoding: quoted-printable Content-Type: text/plain; charset=iso-8859-1 Dear Pietro, If you are planning to use QSite, there is a lot of information in the = Schrodinger web page (https://www.schrodinger.com/supportfaq/18/19/). Also, you can search the User manual in the documentation folder which = is available when you make the software installation. Best regards, Jans El Jan 18, 2011, a las 11:57 AM, Pietro Ottati pietro.ottati-.-gmail.com = escribi=F3: > Hi all, I'm a Junior Researcher in these days I'm starting my first = time work=20 > in Computational Chemistry. My job will base on QM/MM the software = that I have=20 > to use is Qsite Schrodinger Suite. Unfortunately there aren't a lot of = guide,=20 > tutorial etc... about this software... Can you help me to learn? Do = you have=20 > some links, books, articles, wiky etc? Anybody use it? >=20 > Please take in account that I'm newbie and this is my first time. >=20 > Thanks in advance=20 >=20 > Pietro ----------------------------------------------------------- Jans Alzate-Morales Pharmaceutical Chemist, PhD. in Chemistry Centre for Bioinformatics and Molecular Simulations (CBSM) Universidad de Talca 2 Norte 685, Casilla 721, Talca - Chile =20 Tel=E9fono: 56-71-201 798,=20 Fax: 56-71-201 561 jalzate,utalca.cl --Apple-Mail-23--757596058 Content-Transfer-Encoding: quoted-printable Content-Type: text/html; charset=iso-8859-1 Dear = Pietro,

If you are planning to use QSite, there is a = lot of information in the Schrodinger web page (https://www.schrodi= nger.com/supportfaq/18/19/).
Also, you can search the User = manual in the documentation folder which is available when you make the = software installation.

Best = regards,

Jans


<= div>El Jan 18, 2011, a las 11:57 AM, Pietro Ottati pietro.ottati-.-gmail.com = escribi=F3:

Hi all, I'm a Junior =
Researcher in these days I'm starting my first time work=20
in Computational Chemistry. My job will base on QM/MM the software that =
I have=20
to use is Qsite Schrodinger Suite. Unfortunately there aren't a lot of =
guide,=20
tutorial etc... about this software... Can you help me to learn? Do you =
have=20
some links, books, articles, wiky etc? Anybody use it?

Please take in account that I'm newbie and this is my first time.

Thanks in advance=20

Pietro

Jans = Alzate-Morales
Pharmaceutical Chemist, PhD. in = Chemistry
Centre for Bioinformatics and Molecular Simulations = (CBSM)
Universidad de Talca
2 Norte 685, Casilla 721, Talca = - Chile  
Tel=E9fono: 56-71-201 798, 
Fax: 56-71-201 = 561




= --Apple-Mail-23--757596058-- From owner-chemistry@ccl.net Tue Jan 18 15:21:00 2011 From: "Nancy nancy5villa++gmail.com" To: CCL Subject: CCL: Pioglitazone Tautomers Message-Id: <-43688-110118135433-16548-gnj+5L9ixXE0dGg0NMHGlQ . server.ccl.net> X-Original-From: Nancy Content-Type: multipart/alternative; boundary=00248c11dcc9d1e3e5049a236c84 Date: Tue, 18 Jan 2011 13:54:07 -0500 MIME-Version: 1.0 Sent to CCL by: Nancy [nancy5villa.^-^.gmail.com] --00248c11dcc9d1e3e5049a236c84 Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Thank you very much; that answers all my questions. Nancy On Tue, Jan 18, 2011 at 1:49 AM, Andreas Klamt klamt(!)cosmologic.de < owner-chemistry^-^ccl.net> wrote: > Sent to CCL by: Andreas Klamt [klamt]![cosmologic.de] > Nancy, > > although I find Gabriele's reply a bit to harsh, I also was a bit > frustrated when I read your last entry, which essentially repeats the > initial question after you got a number of replies. You now should draw y= our > own conclusions on this. > > What you can learns from the diverse answers is that tautomerism is still= a > complicated and quite unsolved question. This was also shown in the speci= al > issue of J Comput Aided Mol Des. 2010 Jun,24 on tautomerism edited by > Yvonne Martin and by the outcome of the SAMPLE2 blindtest (see Journal of > Computer-Aided Molecular Design, Volume 24, Number 4). > > For your case of Pioglitazone I would conclude that the quantum chemical > results for the tautomer energy differences are so large, that a rather > unique answer can be given: The diketonic form (neutral and deprotonated) > are absolutely predominant. This is a lucky case: There are many cases of > tatuomerism where the DFT energy differences are in the range of a few > kcal/mol, i.e. within the error bars of DFT, and then no such clear answ= er > can be given. > > I would go so far to conclude that the huge energy differences of ~ 12 > kcal/mol to the other forms also makes the other forms rather unlikely in > any protein docked form. I cannot see how the local environment of the > protein can be so much more favoring the enolic forms to compensate such > penalty. > > If the pKa of 6.3 reported by Gabriele is an experimental one, I have no > problem to admit that most likely his 85% deprotonated is the better answ= er > than my answer based on DFT/COSMO-RS prediction of pKa=3D6.55. The differ= ence > of 0.25 pKa units is far within the error bars (~ 0.5) of our pKa > prediction. > > Good luck in going on from here. Perhaps you can send me your final paper= , > or thesis or what ever kind of document you will produce on this case. > > Andreas > > > > Am 17.01.2011 22:14, schrieb Gabriele Cruciani gabri- - > chemiome.chm.unipg.it: > > > Sent to CCL by: Gabriele Cruciani [gabri|a|chemiome.chm.unipg.it] > Nancy, > you get the same answer from different people, with different background, > using different methods, and you keep asking the same question. > > I think CCL list should be used 'with caution' especially after you get a= ll > the answers. > > What you summarized is not correct. The two forms are not tautomeric form= s. > They refer to ONE tautomeric form that is in equilibrium with two species= , > one protonated and one neutral. The percentages you reported refer to the > solvent pH and to the specie pKa. Since the pKa of the pioglitazone is 6.= 3, > the relative abundance at pH 7.4 is > 85% for the charged and 15% for the neutral one. > > As told in previous mails, this has nothing to do with the docking into a > protein, where the local environment, dielectric, 3D interactions may cha= nge > the 'form' of the interacting ligand. In this specific case, you cannot s= ay > that the charge form may interact better because is more abundant, since = the > equilibrium between the two forms can be modified by the protein. > > Precise pKa calculations are very relevant for pharmacokinetic and > metabolism studies, but for docking to receptors is better to produce the > 'relevant' structures. Do not use too simple methods (like those you > reported) since the problem is very complex and require good approaches. > Often using the brain is better than using free available software. > > Gabriele Cruciani > > > > Thus, the two major tautomers of pioglitazone are as shown in the attache= d > figure? > > Also, is there a published article, or other evidence that would confirm > that this is indeed the correct tautomeric distribution at pH 7? > > Thanks you, > Nancy > > > > On Mon, Jan 17, 2011 at 6:17 AM, Andreas Klamt klamt#cosmologic.de< > owner-chemistry/a\ccl.net> wrote: > > Sent to CCL by: Andreas Klamt [klamt-x-cosmologic.de] > I was just pointed by Leif N=F8rskov-Lauritsen to a stupid error which I = made > > in my first entry on this topic. I wrote that at pH=3D7 25% would be > deprotonated, 75% neutral. That was obviously wrong. It is the other way > round: 72% is ionic, 25% neutral. The total amount of zwitterions is then > ~3%. > > Sorry for my mistake. But the basic message about the most stable neutral > form stays valid. > The most stable anionic form is the di-keto anion. > > Andreas > > Am 16.01.2011 19:29, schrieb Nancy nancy5villa**gmail.com: > > Using quantum chemical calculations, you found that the diketone form is > dominant; does this mean that the enol form predicted by MarvinSketch > v3.5.4 > in the published article is wrong? > > "Structure-based design of a thiazolidinedione which targets the > mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 Feb > 1;20(3):819-23 > > Thanks, > Nancy > > > On Thu, Jan 13, 2011 at 10:31 PM, Andreas Klamt klamt(a)cosmologic.de< > owner-chemistry(~)ccl.net> wrote: > > Sent to CCL by: Andreas Klamt [klamt,cosmologic.de] > Hi Nancy, > > as I wrote in the previous message, we find by quantum chemical > calculations within the COSMO-RS solvation that the diketone form is more > stable by abot 12 kcal/mol, that means it absolutely dominates! > > Andreas > > Am 13.01.2011 22:25, schrieb Nancy nancy5villa~!~gmail.com: > > Hi All, > > I happened to find a published article where molecular docking simulation= s > of TZDs against a novel protein is detailed: > > "Structure-based design of a thiazolidinedione which targets the > mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 Feb > 1;20(3):819-23 > > In this paper, the authors, using MarvinSketch v3.5.4, determined that th= e > TZDs would exist predominantly in the enol form at pH 7.4, as opposed to > the > diketone form (see attached figure). It appears that the article's > tautomer > prediction was based only on the pH prior to docking, and had nothing to = do > > with ligand-protein interactions. > > Does anyone know if the diketone or enol, or perhaps a different tautomer= , > would be predominant at pH 7.4? > > Thanks in advance. > Nancy > > > > On Wed, Jan 12, 2011 at 9:53 AM, Andreas Klamt klamt-*-cosmologic.de< > owner-chemistry(_)ccl.net> wrote: > > > Sent to CCL by: Andreas Klamt [klamt[*]cosmologic.de] > Hi All, > > let me put my 5 cents into the game: > We have done quantum calculations for the heterocyclic ring, with just a > methyl group, since electronically the remainder of the compound should > have > almost no influence on the heterocycle. > DFT/COSMO-RS calculations prove the neutral species > SMILES:CC1C(=3DO)[NH]C(=3DO)S1 of the Pioglitazone heterocycle to be favo= red > by 11.3 kcal/mol compared to all other neutral species of the heterocycle= , > which means that there is no other relvant neutral form. > > We find with COSMOtherm that the pka of the heterocycle is 6.55, i.e. at > pH=3D7 it should be deprotonated with about 25%. > Since the pka of the conj. acid of the pyridine ring according to > COSMOtherm (we calculated diethylpyridine) is ~5.6 there should be ~10% i= n > the zwitterionic form. > > It will be hard to decide this finally, unless someone measures it, but > we are quite confident to be qualitatively right here, because the energy > differences for the tautomers should not have more than 5 kcal/mol erro, = as > > we learned in in the SAMPL2 tautomer contest: > see Andreas Klamt and Michael Diedenhofen, Some conclusions regarding the > predictions of tautomeric equilibria in solution based on the SAMPL2 > challenge, http://www.springerlink.com/content/l577667th758n6h1/ > > Andreas > > > > > > > > Hi All, > > Using MarvinSketch v5.3.3 (ChemAxon software), the predicted pKa value > of the acidic hydrogen on the thiazolidinedione ring of pioglitazone is > 4.57 (see attached "Figure_1.gif"). Therefore, the predominant species > at pH 7.0 predicted by MarvinSketch is the one depicted in > "Figure_2.gif". The different tautomeric forms predicted by MarvinSketch > are shown in "Figure_3.gif". > > Can anyone explain where these predictions come from, and if they are > correct? > > Thanks, > Nancy > > > On Mon, Jan 10, 2011 at 4:16 AM, Gabriele Cruciani > gabri(0)chemiome.chm.unipg.it > > > > wrote: > > > Sent to CCL by: Gabriele Cruciani [gabri##chemiome.chm.unipg.it > ] > > Nancy, > the first form you reported is the most stable in water at pH 7.0. > However, the fact that one form is more stable than another in water > does not help you to understand which form will be more 'relevant' > for docking. In protein the tautomeric equilibria may produce and > stabilize different forms according to the complementary site. There > are examples of tautomeric form in protein not stable in water, > where the energy difference is more than 5 Kcal/mol. > > MoKa software (www.moldiscovery.com) > > > is fast and accurate to produce tautomer stable in water, but it can > produce also all the (plausible) tautomeric forms. > > Then, a possibility is to dock them into your protein, to see if > docking methods may differentiate their binding. > > Gabriele Cruciani > > > > > Hi All, > > I am performing molecular docking and molecular dynamics > simulations of the > thiazolidinedione pioglitazone binding to the PPAR-gamma > receptor protein > (PDB ID: 1ZGY). The thiazolidinedione ring can exist in numerous > different > tautomeric states; I have attached a figure depicting several of > them. > Which tautomer would be dominant at the physiological pH of ~7.0? > > Also, are there any software programs that can predict which > tautomer would > be correct? > > Thanks in advance, > Nancy> > > > E-mail to subscribers: CHEMISTRY[a]ccl.net > or use:> > E-mail to administrators: CHEMISTRY-REQUEST[a]ccl.net > o= r > useConferences: > > http://server.ccl.net/chemistry/announcements/conferences/> > > > > > -- > Dr. Jens Reinisch > Chemist / Customer Support > COSMOlogic GmbH& Co. KG > Burscheider Strasse 515 > D-51381 Leverkusen, Germany > > phone +49-2171-363664 > mobile +49-163-7337310 > fax +49-2171-731689 > e-mail reinisch[*]cosmologic.de > web www.cosmologic.de > > HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt > Komplementaer: COSMOlogic Verwaltungs GmbH > HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt > > > > > > > E-mail to subscribers: CHEMISTRY(_)ccl.net or > use:>>> > E-mail to administrators: CHEMISTRY-REQUEST(_)ccl.net > or useConferences: > http://server.ccl.net/chemistry/announcements/conferences/>>> > > > > > -- > PD. Dr. Andreas Klamt > CEO / Gesch=E4ftsf=FChrer > COSMOlogic GmbH& Co. KG > Burscheider Strasse 515 > D-51381 Leverkusen, Germany > > phone +49-2171-731681 > fax +49-2171-731689 > e-mail klamt(_)cosmologic.de > > web www.cosmologic.de > > HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt > Komplementaer: COSMOlogic Verwaltungs GmbH > HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt > > CHEMISTRY(~)ccl.net or use:E-mail to > administrators: > CHEMISTRY-REQUEST(~)ccl.net or > usehttp://www.ccl.net/chemistry/sub_unsub.shtml > Before > posting, check wait time at: > http://www.ccl.netConferences: > http://server.ccl.net/chemistry/announcements/conferences/ Search > Messages: http://www.ccl.net/chemistry/searchccl/index.shtmlRTFI: > http://www.ccl.net/chemistry/aboutccl/instructions/ > > > > > -- > PD. Dr. Andreas Klamt > CEO / Gesch=E4ftsf=FChrer > COSMOlogic GmbH& Co. KG > Burscheider Strasse 515 > D-51381 Leverkusen, Germany > > phone +49-2171-731681 > fax +49-2171-731689 > e-mail klamt:cosmologic.de > web www.cosmologic.de > > HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt > Komplementaer: COSMOlogic Verwaltungs GmbH > HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas KlamtTo > recover the email address of the author of the message, please change the > strange characters on the top line to the /a\ sign. You can also look up > the > X-Original-From: line in the mail header. E-mail to subscribers: > CHEMISTRY/a\ccl.net or use:E-mail to administrators: > CHEMISTRY-REQUEST/a\ccl.net or useBefore posting, check wait > time at: http://www.ccl.netConferences: > http://server.ccl.net/chemistry/announcements/conferences/ Search > Messages: http://www.ccl.net/chemistry/searchccl/index.shtml If your mail > bounces from CCL with 5.7.1 error, check:RTFI: > http://www.ccl.net/chemistry/aboutccl/instructions/ > > > > > > E-mail to subscribers: CHEMISTRY%ccl.net or use:> > E-mail to administrators: CHEMISTRY-REQUEST%ccl.net or useConferences: > http://server.ccl.net/chemistry/announcements/conferences/> > > > > > -- > PD. Dr. Andreas Klamt > CEO / Gesch=E4ftsf=FChrer > COSMOlogic GmbH & Co. KG > Burscheider Strasse 515 > D-51381 Leverkusen, Germany > > phone +49-2171-731681 > fax +49-2171-731689 > e-mail klamt%cosmologic.de > > web www.cosmologic.de > > HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt > Komplementaer: COSMOlogic Verwaltungs GmbH > HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt > > -=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 change the > strange characters on the top line to the ^-^ sign. You can also look up th= e > X-Original-From: line in the mail header. E-mail to subscribers: > CHEMISTRY^-^ccl.net or use:= E-mail to administrators: > CHEMISTRY-REQUEST^-^ccl.net or useBefore posting, check wait > time at: http://www.ccl.netConferences: > http://server.ccl.net/chemistry/announcements/conferences/ Search > Messages: http://www.ccl.net/chemistry/searchccl/index.shtml If your mail > bounces from CCL with 5.7.1 error, check:= RTFI: > http://www.ccl.net/chemistry/aboutccl/instructions/ > --00248c11dcc9d1e3e5049a236c84 Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Thank you very much; that answers all my questions.

Nancy


On Tue, Jan 18, 2011 at 1:49 AM, Andreas Kl= amt klamt(!)cosmologic.de <owner-chemistry^-^ccl.ne= t> wrote:
Sent to CCL by: Andreas Klamt [klamt]![cosmologic.de] =20 =20 =20 =20
Nancy,

although I find Gabriele's reply a bit to harsh, I also was a bit frustrated when I read your last entry, which essentially repeats the initial question after you got a number of replies. You now should draw your own conclusions on this.

What you can learns from the diverse answers is that tautomerism is still a complicated and quite unsolved question. This was also shown in the special issue of J Comput Aided Mol Des. 2010 Jun,24 on tautomerism edited by Yvonne Martin and by the outcome of the SAMPLE2 blindtest (see Journal of Computer-Aided Molecular Design, Volume 24, Number 4).

For your case of Pioglitazone I would conclude that the quantum chemical results for the tautomer energy differences are so large, that a rather unique answer can be given: The diketonic form (neutral and deprotonated) are absolutely predominant. This is a lucky case: There are many cases of tatuomerism where the DFT energy differences are in the range of a few kcal/mol, i.e. within the error bars of DFT,=A0 and then no such clear answer can be given.

I would go so far to conclude that the huge energy differences of ~ 12 kcal/mol to the other forms also makes the other forms rather unlikely in any protein docked form. I cannot see how the local environment of the protein can be so much more favoring the enolic forms to compensate such penalty.

If the pKa of 6.3 reported by Gabriele is an experimental one, I have no problem to admit that most likely his 85% deprotonated is the better answer than my answer based on DFT/COSMO-RS prediction of pKa=3D6.55. The difference of 0.25 pKa units is far within the error bars (~ 0.5) of our pKa prediction.

Good luck in going on from here. Perhaps you can send me your final paper, or thesis or what ever kind of document you will produce on this case.

Andreas



Am 17.01.2011 22:14, schrieb Gabriele Cruciani gabri- -chemiome.ch= m.unipg.it:

Sent to CCL by: Gabriele Cruciani [gabri|a|chemiome.chm.unipg.it]
Nancy,
you get the same answer from different people, with different background, using different methods, and you keep asking the same question.

I think CCL list should be used 'with caution' especially aft= er you get all the answers.

What you summarized is not correct. The two forms are not tautomeric forms. They refer to ONE tautomeric form that is in equilibrium with two species, one protonated and one neutral. The percentages you reported refer to the solvent pH and to the specie pKa. Since the pKa of the pioglitazone is 6.3, the relative abundance at pH 7.4 is
85% for the charged and 15% for the neutral one.

As told in previous mails, this has nothing to do with the docking into a protein, where the local environment, dielectric, 3D interactions may change the 'form' of the interacting ligand.= In this specific case, you cannot say that the charge form may interact better because is more abundant, since the equilibrium between the two forms can be modified by the protein.

Precise pKa calculations are very relevant for pharmacokinetic and metabolism studies, but for docking to receptors is better to produce the 'relevant' structures. Do not use too simple meth= ods (like those you reported) since the problem is very complex and require good approaches. Often using the brain is better than using free available software.

Gabriele Cruciani



Thus, the two major tautomers of pioglitazone are as shown in the attached
figure?

Also, is there a published article, or other evidence that would confirm
that this is indeed the correct tautomeric distribution at pH 7?

Thanks you,
Nancy



On Mon, Jan 17, 2011 at 6:17 AM, Andreas Klamt klamt#cosmologic= .de<
owner-chemistry/a\ccl.= net>=A0 wrote:

Sent to CCL by: Andreas Klamt [klamt= -x-cosmologic.de]
I was just pointed by Leif N=F8rskov-Lauritsen to a stupid error which I made
in my first entry on this topic. I wrote that at pH=3D7 25% would be
deprotonated, 75% neutral. That was obviously wrong. It is the other way
round: 72% is ionic, 25% neutral. The total amount of zwitterions is then
~3%.

Sorry for my mistake. But the basic message about the most stable neutral
form stays valid.
The most stable anionic form is the di-keto anion.

Andreas

Am 16.01.2011 19:29, schrieb Nancy nancy5villa**gmail.com:

Using quantum chemical calculations, you found that the diketone form is
dominant; does this mean that the enol form predicted by MarvinSketch v3.5.4
in the published article is wrong?

"Structure-based design of a thiazolidinedione which targets the
mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 F= eb
1;20(3):819-23

Thanks,
Nancy


On Thu, Jan 13, 2011 at 10:31 PM, Andreas Klamt klamt(a)cosmol= ogic.de<
owner-chemistry(~)cc= l.net<owner-chemistry%28%7E%29ccl.net>>=A0 wrote:

=A0 Sent to CCL by: Andreas Klamt [klamt,cosmo= logic.de]
=A0 Hi Nancy,

as I wrote in the previous message, we find by quantum chemical
calculations within the COSMO-RS solvation that the diketone form is more
stable by abot 12 kcal/mol, that means it absolutely dominates!

Andreas

Am 13.01.2011 22:25, schrieb Nancy nancy5villa~!~gmail.com:

=A0 Hi All,

I happened to find a published article where molecular docking simulations
of TZDs against a novel protein is detailed:

"Structure-based design of a thiazolidinedione which targe= ts the
mitochondrial protein mitoNEET" Bioorg Med Chem Lett. 2010 Feb
1;20(3):819-23

In this paper, the authors, using MarvinSketch v3.5.4, determined that the
TZDs would exist predominantly in the enol form at pH 7.4, as opposed to the
diketone form (see attached figure).=A0 It appears that the article's tautomer
prediction was based only on the pH prior to docking, and had nothing to do
with ligand-protein interactions.

Does anyone know if the diketone or enol, or perhaps a different tautomer,
would be predominant at pH 7.4?

Thanks in advance.
Nancy



=A0=A0 On Wed, Jan 12, 2011 at 9:53 AM, Andreas Klamt klamt-*-cosm= ologic.de<
owner-chemistry(_)= ccl.net<owner-chemistry%28_%29ccl.net>>=A0 wrote:


Sent to CCL by: Andreas Klamt [klamt[*]cosmologic.de]
Hi All,

let me put my 5 cents into the game:
We have done quantum calculations for the heterocyclic ring, with just a
methyl group, since electronically the remainder of the compound should have
almost no influence on the heterocycle.
DFT/COSMO-RS calculations prove the neutral species
SMILES:CC1C(=3DO)[NH]C(=3DO)S1 of the Pioglitazone heterocycl= e to be favored
by 11.3 kcal/mol compared to all other neutral species of the heterocycle,
which means that there is no other relvant neutral form.

We find with COSMOtherm that the pka of the heterocycle is 6.55, i.e. at
pH=3D7 it should be deprotonated with about 25%.
Since the pka of the conj. acid of the pyridine ring according to
COSMOtherm (we calculated diethylpyridine) is ~5.6 there should be ~10% in
the zwitterionic form.

It will be hard to decide this finally, unless someone measures it, but
we are quite confident to be qualitatively right here, because the energy
differences for the tautomers should not have more than 5 kcal/mol erro, as
we learned in in the SAMPL2 tautomer contest:
see Andreas Klamt and Michael Diedenhofen, Some conclusions regarding the
predictions of tautomeric equilibria in solution based on the SAMPL2
challenge, http://www.springerlink.com/content/l577667th758n6h1= /

Andreas







Hi All,

Using MarvinSketch v5.3.3 (ChemAxon software), the predicted pKa value
of the acidic hydrogen on the thiazolidinedione ring of pioglitazone is
4.57 (see attached "Figure_1.gif"). Therefore, th= e predominant species
at pH 7.0 predicted by MarvinSketch is the one depicted in
"Figure_2.gif". The different tautomeric forms pr= edicted by MarvinSketch
are shown in "Figure_3.gif".

Can anyone explain where these predictions come from, and if they are
correct?

Thanks,
Nancy


On Mon, Jan 10, 2011 at 4:16 AM, Gabriele Cruciani
=A0 gabri(0)chemiome.chm.unipg.it<http://chemiome.chm.unipg.it>

<owner-chemistry[a]ccl.net<mailto:owner-chemistry[a]ccl.net>>=A0 wrote:


=A0=A0=A0 Sent to CCL by: Gabriele Cruciani [gabri##chemiome.chm.unipg.it
=A0=A0=A0=A0 <http://chemiome.chm.unipg.it>]

=A0=A0=A0 Nancy,
=A0=A0=A0 the first form you reported is the most stable in water at pH 7.0.
=A0=A0=A0 However, the fact that one form is more stable th= an another in water
=A0=A0=A0 does not help you to understand which form will b= e more 'relevant'
=A0=A0=A0 for docking. In protein the tautomeric equilibria may produce and
=A0=A0=A0 stabilize different forms according to the complementary site. There
=A0=A0=A0 are examples of tautomeric form in protein not stable in water,
=A0=A0=A0 where the energy difference is more than 5 Kcal/m= ol.

=A0=A0=A0=A0 MoKa software (= www.moldiscovery.com<http://www.moldiscovery.com>)

=A0=A0=A0 is fast and accurate to produce tautomer stable i= n water, but it can
=A0=A0=A0 produce also all the (plausible) tautomeric forms= .

=A0=A0=A0 Then, a possibility is to dock them into your protein, to see if
=A0=A0=A0 docking methods may differentiate their binding.

=A0=A0=A0 Gabriele Cruciani




=A0=A0=A0=A0=A0=A0=A0 Hi All,

=A0=A0=A0=A0=A0=A0=A0 I am performing molecular docking and= molecular dynamics
=A0=A0=A0=A0=A0=A0=A0 simulations of the
=A0=A0=A0=A0=A0=A0=A0 thiazolidinedione pioglitazone bindin= g to the PPAR-gamma
=A0=A0=A0=A0=A0=A0=A0 receptor protein
=A0=A0=A0=A0=A0=A0=A0 (PDB ID: 1ZGY). The thiazolidinedione= ring can exist in numerous
=A0=A0=A0=A0=A0=A0=A0 different
=A0=A0=A0=A0=A0=A0=A0 tautomeric states; I have attached a = figure depicting several of
=A0=A0=A0=A0=A0=A0=A0 them.
=A0=A0=A0=A0=A0=A0=A0 Which tautomer would be dominant at t= he physiological pH of ~7.0?

=A0=A0=A0=A0=A0=A0=A0 Also, are there any software programs= that can predict which
=A0=A0=A0=A0=A0=A0=A0 tautomer would
=A0=A0=A0=A0=A0=A0=A0 be correct?

=A0=A0=A0=A0=A0=A0=A0 Thanks in advance,
=A0=A0=A0=A0=A0=A0=A0 Nancy>


=A0=A0=A0=A0 E-mail to subscribers: CHEMISTRY[a]ccl.net
=A0=A0=A0 <mailto:CHEMISTRY[a]ccl.net>=A0 or use:>
=A0=A0=A0 E-mail to administrators: CHEMISTRY-REQUEST[a]ccl.net
=A0=A0=A0 <mailto:CHEMISTRY-REQUEST[a]ccl.net>=A0 or useConferences:

=A0=A0=A0 http://server.ccl.net/chemistry/announcements= /conferences/>




--
Dr. Jens Reinisch
Chemist / Customer Support
COSMOlogic GmbH&=A0 Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone=A0=A0 +49-2171-363664
mobile=A0 +49-163-7337310
fax=A0=A0=A0=A0 +49-2171-731689
e-mail=A0 reinisch[*]cosmologic.de
web=A0=A0=A0=A0 www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt






=A0 E-mail to subscribers: CHEMISTRY(_)ccl.= net<CHEMISTRY%28_%29c= cl.net>=A0 or
use:>>>
E-mail to administrators: CHEMISTRY-REQUEST(_)ccl.net<CHEMISTRY-REQUEST%28_%29ccl.net>or useConferences:
http://server.ccl.net/chemistry/announcements/conferences/>>>




--
PD. Dr. Andreas Klamt
CEO / Gesch=E4ftsf=FChrer
COSMOlogic GmbH&=A0 Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone=A0=A0=A0=A0=A0 +49-2171-731681
fax=A0=A0=A0=A0=A0=A0=A0 +49-2171-731689
e-mail=A0=A0=A0=A0 klamt(_)
cosm= ologic.de<klamt%28_%29cosmologic.de>

web=A0=A0=A0=A0=A0=A0=A0 www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt

CHEMISTRY(~)ccl.ne= t<CHEMISTRY%28%7E%29c= cl.net>=A0 or use:E-mail to administrators:
CHEMISTRY-REQUEST(~)ccl.net<CHEMISTRY-REQUEST%28%7E%29ccl.net>=A0 or usehttp://www.ccl.net/chemistry/sub_unsub.shtml<http:= //www.ccl.net/cgi-bin/ccl/send_ccl_message>Before posting, check wait time at:
ht= tp://www.ccl.netConferences:
http://server.ccl.net/chemistry/announcements/con= ferences/ Search
Messages: http://www.ccl.net/chemistry/searchccl/index.shtmlRT= FI:
http://www.ccl.net/chemistry/aboutccl/instructions/




--
PD. Dr. Andreas Klamt
CEO / Gesch=E4ftsf=FChrer
COSMOlogic GmbH&=A0 Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone=A0=A0=A0=A0=A0 +49-2171-731681
fax=A0=A0=A0=A0=A0=A0=A0 +49-2171-731689
e-mail=A0=A0=A0=A0 klamt:cosmologic.de
web=A0=A0=A0=A0=A0=A0=A0 www.cosmologic.de

HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas KlamtTo
recover the email address of the author of the message, please change the
strange characters on the top line to the /a\ sign. You can also look up the
X-Original-From: line in the mail header. E-mail to subscribers:
CHEMISTRY/a\ccl.net<= /a> or use:E-mail to administrators:
CHEMISTRY-REQUEST/a\
= ccl.net or useBefore posting, check wait
time at: http://www.ccl.netConferences:
http://server.ccl.net/chemistry/announcements/confe= rences/ Search
Messages: http://www.ccl.net/chemistry/searchccl/index.shtm= l If your mail
bounces from CCL with 5.7.1 error, check:RTFI:
http://www.ccl.net/chemistry/aboutccl/instructions/




=20
E-mail to subscribers: CHEMISTRY%ccl.net or use:
=A0=A0=A0=A0 http://www.ccl.net/cgi-bin/ccl/send_ccl_message

E-mail to administrators: CHEMISTRY-REQUEST%ccl.net or use
=A0=A0=A0=A0 http://www.ccl.net/cgi-bin/ccl/send_ccl_message


--=20
PD. Dr. Andreas Klamt
CEO / Gesch=E4ftsf=FChrer
COSMOlogic GmbH & Co. KG
Burscheider Strasse 515
D-51381 Leverkusen, Germany

phone  	+49-2171-731681
fax    	+49-2171-731689
e-mail 	klamt%co=
smologic.de
web www.cosmolog= ic.de HRA 20653 Amtsgericht Koeln, GF: Dr. Andreas Klamt Komplementaer: COSMOlogic Verwaltungs GmbH HRB 49501 Amtsgericht Koeln, GF: Dr. Andreas Klamt

--00248c11dcc9d1e3e5049a236c84-- From owner-chemistry@ccl.net Tue Jan 18 15:56:00 2011 From: "Vincent Leroux vincent.leroux _ loria.fr" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43689-110118150104-30541-ukWNFQhoqDNVhaM7hzydHQ|*|server.ccl.net> X-Original-From: Vincent Leroux Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=ISO-8859-1 Date: Tue, 18 Jan 2011 21:00:24 +0100 MIME-Version: 1.0 Sent to CCL by: Vincent Leroux [vincent.leroux[*]loria.fr] Hi, "schrodinger qsite manual" on Google gives you the pdf manual in the top answers... "schrodinger qsite" on Google Scholar gives you plenty of published papers that made use of this software... Ok, you are a newbie, but come on... VL Le 18/01/11 15:57, Pietro Ottati pietro.ottati-.-gmail.com a écrit : > Hi all, I'm a Junior Researcher in these days I'm starting my first time work > in Computational Chemistry. My job will base on QM/MM the software that I have > to use is Qsite Schrodinger Suite. Unfortunately there aren't a lot of guide, > tutorial etc... about this software... Can you help me to learn? Do you have > some links, books, articles, wiky etc? Anybody use it? > > Please take in account that I'm newbie and this is my first time. > > Thanks in advance > > Pietro > From owner-chemistry@ccl.net Tue Jan 18 17:23:00 2011 From: "Pietro Ottati pietro.ottati,,gmail.com" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43690-110118171953-4436-9286R8MfinECTAe7UrbDWQ:-:server.ccl.net> X-Original-From: Pietro Ottati Content-Type: multipart/alternative; boundary=0015175ca7bca9e4c0049a264b9e Date: Tue, 18 Jan 2011 23:19:35 +0100 MIME-Version: 1.0 Sent to CCL by: Pietro Ottati [pietro.ottati(-)gmail.com] --0015175ca7bca9e4c0049a264b9e Content-Type: text/plain; charset=ISO-8859-1 Oh come on people do you think that I post here and I didn't read the manual? I've just read it... But Qsite manual is really uncomplete.... I've just studied Qsite manual, jaguar quick guide, Maestro quick guide, maestro installation guide and some more articles... But this is really not enough I need a complete manual to understand my new work. This is my actual position: I graduated in Medicinal chemistry in Italy (5 years degree) my thesis work was about organic synthesis but now I'm moving to computational chemistry because I love Computers, I'm good at operating system, network administration and a little bit at programming. Now I need some basis in quantum mechanics chemistry applied in computational work. For me now QM/MM is a parallel world and I need to discover it... I'm sorry if my discussion is not so clear... Thanks in advance byebye --0015175ca7bca9e4c0049a264b9e Content-Type: text/html; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Oh come on people do you think that I post here and I didn't read the m= anual? I've just read it... But Qsite manual is really uncomplete.... I= 've just studied Qsite manual, jaguar quick guide, Maestro quick guide,= maestro installation guide and some more articles... But this is really no= t enough I need a complete manual to understand my new work. This is my act= ual position: I graduated in Medicinal chemistry in Italy (5 years degree) = my thesis work was about organic synthesis but now I'm moving to comput= ational chemistry because I love Computers, I'm good at operating syste= m, network administration and a little bit at programming. Now I need some = basis in quantum mechanics chemistry applied in computational work.

For me now QM/MM is a parallel world and I need to discover it...
<= br>I'm sorry if my discussion is not so clear...

Thanks in advan= ce

byebye
--0015175ca7bca9e4c0049a264b9e-- From owner-chemistry@ccl.net Tue Jan 18 20:07:00 2011 From: "Jun Zhang coolrainbow{:}yahoo.cn" To: CCL Subject: CCL: Looking for a source code for calculating Rys Roots and Weights Message-Id: <-43691-110118200621-16038-R8K9N5rHGN6PotM3Y4j8wA|-|server.ccl.net> X-Original-From: Jun Zhang Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=utf-8 Date: Wed, 19 Jan 2011 09:06:04 +0800 (CST) MIME-Version: 1.0 Sent to CCL by: Jun Zhang [coolrainbow a yahoo.cn] Hello: It seems that GAMESS can only compute the roots of rys polynomial up to 13 degree. If you need a general one(such as computing ERIs up to i angular), you can refer to Dalton's Hermite module. Cheers up! Jun Zhang Nankai University coolrainbow++yahoo.cn From owner-chemistry@ccl.net Tue Jan 18 20:42:00 2011 From: "Vincent Leroux vincent.leroux(~)loria.fr" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43692-110118185407-28184-QmQDXrXaqrJcdlE9K+OPLw:_:server.ccl.net> X-Original-From: Vincent Leroux Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=ISO-8859-1 Date: Wed, 19 Jan 2011 00:53:50 +0100 MIME-Version: 1.0 Sent to CCL by: Vincent Leroux [vincent.leroux_+_loria.fr] Hi Pietro, Sorry if I was a little bit rude. If you have some experience with computer science you might be used to large, extensive documentation. Most of the times, this simply does not exist in the computational chemistry field. Welcome a parallel world! If you are into programming, you probably dislike working with undocumented code from others. Better get used to it. And more often than you think, this will be old-school Fortran code, to add insult to injury ;-) QSite, with its 80-page incomplete manual and plenty of published papers that made use of it, is indeed, IMHO, above average regarding the level of documentation you usually get. So if you encounter some obscure option you have no idea of the purpose or effects, this is all normal. You can ask Schrödinger support for explanations but do not expect much. You have nothing to lose here if you have specific questions. Another option is to contact the corresponding authors of some QSite-using papers you find the most interesting. This might pay the most in the long term. I can certify you that your current skills will not be useless if you go forward into computer-aided chemistry. You can also try to get used to strategies closer to your past expertise, like QSAR and other ligand-based drug design techniques, before jumping into more difficult, low-level theoretical playing grounds. Good luck VL Le 18/01/11 23:19, Pietro Ottati pietro.ottati,,gmail.com a écrit : > Oh come on people do you think that I post here and I didn't read the > manual? I've just read it... But Qsite manual is really uncomplete.... > I've just studied Qsite manual, jaguar quick guide, Maestro quick guide, > maestro installation guide and some more articles... But this is really > not enough I need a complete manual to understand my new work. This is > my actual position: I graduated in Medicinal chemistry in Italy (5 years > degree) my thesis work was about organic synthesis but now I'm moving to > computational chemistry because I love Computers, I'm good at operating > system, network administration and a little bit at programming. Now I > need some basis in quantum mechanics chemistry applied in computational > work. > > For me now QM/MM is a parallel world and I need to discover it... > > I'm sorry if my discussion is not so clear... > > Thanks in advance > > byebye From owner-chemistry@ccl.net Tue Jan 18 21:17:00 2011 From: "Brian Austin brianmaustin*gmail.com" To: CCL Subject: CCL: Hello list... Help me first time. Message-Id: <-43693-110118182714-6488-9286R8MfinECTAe7UrbDWQ{=}server.ccl.net> X-Original-From: Brian Austin Content-Type: text/plain; charset=ISO-8859-1 Date: Tue, 18 Jan 2011 15:26:40 -0800 MIME-Version: 1.0 Sent to CCL by: Brian Austin [brianmaustin#gmail.com] Hello Pietro, I haven't used QSite for about ten years, but I do remember that Schrodinger's customer support is friendly and generous. Consider contacting them directly http://www.schrodinger.com/supportcenter/ They also have training sessions fairly regularly. http://www.schrodinger.com/events/ -Brian On Tue, Jan 18, 2011 at 2:19 PM, Pietro Ottati pietro.ottati,,gmail.com wrote: > Oh come on people do you think that I post here and I didn't read the > manual? I've just read it... But Qsite manual is really uncomplete.... I've > just studied Qsite manual, jaguar quick guide, Maestro quick guide, maestro > installation guide and some more articles... But this is really not enough I > need a complete manual to understand my new work. This is my actual > position: I graduated in Medicinal chemistry in Italy (5 years degree) my > thesis work was about organic synthesis but now I'm moving to computational > chemistry because I love Computers, I'm good at operating system, network > administration and a little bit at programming. Now I need some basis in > quantum mechanics chemistry applied in computational work. > > For me now QM/MM is a parallel world and I need to discover it... > > I'm sorry if my discussion is not so clear... > > Thanks in advance > > byebye >