Thanks. I'll start tomorrow learning about the procedures you suggest. As I
said, I am new to MD. I am not bound to find pathways, I'll be happy initially
to sample the conformational space in an explicit solvent.
90:10 population of conformers, ie 1.3 kcal/mol at rt, is an ideal situation
for conformational studies (i.e NMR techniques can deal with so that one is not
computing blind). Conformational energy differences are always in this range
with natural products (or organic synthetic drugs). If MD can't deal with such
energy differences, it will be of little use in this area and the interaction
of with receptors. However, I vaguely remember of published MD dealing with
such small energy differences. Perhaps I can find out the paper spending some
time.
At any event, hope someone can suggest more explicitly how to deal with the
procedures below. I would be surprised that no one is on such problems. Surely
there are people with Monte Carlo Metropolis.
Thanks
francesco pietra
--- Carlos Simmerling <carlos.simmerling.gmail.com> wrote:
> I don't have much experience with using NEB for this type of change, perhaps
>
> if you just want relative free energies you might be better off with free
> energy
> methods such as thermodynamic integration. NEB does not provide free
> energies and is not yet implemented in Amber for use with explicit solvent.
> If you can define a specific reaction coordinate you may use umbrella
> sampling.
> If you want to try to reproduce a 10:90 ratio, you may want to spend some
> time generating accurate parameters since this is a relatively small free
> energy difference (a little over 1 kcal).
>
> On 7/5/07, Francesco Pietra <chiendarret.yahoo.com> wrote:
> >
> > I fear that my answer to Carlos Simmerling about NEB was not dispatched
> > correctly, or not dispatched at all.
> >
> > I said to have carried out tutorials antechamcer-sustiva and NBE with
> > amber9,
> > aimed at setting up simulated annealing for natural products.
> > Antechamber-gaff,
> > with little additional computing, proved able to deal with such large
> > molecules
> > without repeating units. Partial charges that I had determined by formal
> > cutting of the molecule and applying HF-6-21G* methodology for RESP,
> > differ
> > substantially from those calculated by divcon-AM1.
> >
> > I was asked what I aimed to: querying about all coordinates, with only
> > possible
> > restraint for benzenoid rings. Inversion about chirality axes (which
> > determines
> > interconversion of diastereomers because of other chiral center present)
> > faces
> > a barrier of about 24 kcal/mol, values being solvent-dependent. Energy
> > difference between disteromers involved in the conformational change is
> > also
> > solvent dependent. A 90:10 population was detected by NMR in chloroform
> > solution.
> >
> > Simulation is intended for explicit solvents, also water. The final
> > (future)
> > aim is the study of interaction with receptors.
> >
> > I pointed out that I am new at simulated annealing, so that a modification
> > of
> > scripts for the NEB tutorial would be the best, if possible at all.
> >
> > Finally, I also said to have already carried out a global space
> > conformational
> > search by MM (MMFF94 ff) in vacuum. I also have computed the energy
> > difference
> > between disteromers in vacuum by ab initio methods, carrying out a
> > vibrational
> > analysis, correcting for ZPE and temperature, finding fair agreement with
> > experimental observation for chloroform solution only.
> >
> > Tanks
> >
> > francesco pietra
> >
> >
> >
> >
>
____________________________________________________________________________________
> > Park yourself in front of a world of choices in alternative vehicles.
> > Visit the Yahoo! Auto Green Center.
> > http://autos.yahoo.com/green_center/
> > -----------------------------------------------------------------------
> > The AMBER Mail Reflector
> > To post, send mail to amber.scripps.edu
> > To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
> >
>
____________________________________________________________________________________
Fussy? Opinionated? Impossible to please? Perfect. Join Yahoo!'s user panel and lay it on us.
http://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jul 08 2007 - 06:07:28 PDT
