Re: [AMBER] Fw: parameterization of ligand+protein complex

From: Hira Jabeen <hira_bioinfo.yahoo.com>
Date: Thu, 20 Jun 2019 11:04:27 +0000 (UTC)

 Hi,
Like, this is pdb file: 
HETATM 4017 HO72 SFR A 269      22.818  12.980  10.549  0.84 36.34           H
HETATM 4018 HO62 SFR A 269      21.824  15.969   6.233  0.84 37.27           H
HETATM 4019 HO31 SFR A 269      18.325  17.622  12.473  0.84 26.32           H
HETATM 4020 H62B SFR A 269      24.127  13.353   8.201  0.84 37.22           H
HETATM 4021 H62A SFR A 269      22.801  12.946   7.085  0.84 37.22           H
HETATM 4022 H17A SFR A 269      19.459  22.308   7.884  0.84 46.88           H
HETATM 4023 H16A SFR A 269      20.971  22.079   9.766  0.84 46.47           H
HETATM 4024 H15A SFR A 269      20.450  20.229  10.526  0.84 45.48           H
And this is mol2:
     17 C17         20.2460    21.5770     7.6980 c3       301 SFR       0.123600
     18 O18         19.7130    20.2900     7.3990 os       301 SFR      -0.407600
     19 O32         17.6250    17.4610    10.3710 o        301 SFR      -0.482300
     20 O31         19.1320    17.4760    11.9360 o        301 SFR      -0.482300
     21 H6          21.1830    14.7760     8.2140 hc       301 SFR       0.142700
     22 H4          21.3740    16.5880    10.8730 hn       301 SFR       0.376700
     23 H5          22.8750    17.1310     8.9520 h2       301 SFR       0.184700
     24 H2          18.9100    17.7500     7.8870 h1       301 SFR       0.171700

how should I change the highlighted atoms? Also in mol2 file few atoms are missing-how should I recover that?
Thanks,Hira

    On Wednesday, June 19, 2019, 03:37:49 PM GMT+1, Matias Machado <mmachado.pasteur.edu.uy> wrote:
 
 Dear Hira Jabeen,

Please dig into the mailing list, most answers to your problems are there...

Indeed, I've just answered the issue of properly loading mol2 files two weeks ago!... [http://archive.ambermd.org/201906/0032.html]

Other very frequent mistakes are:

1) Atom and residue names in the PDB file HAVE TO match those in the residue library (e.g. mol2).

  That's causing this error:
  "Created a new atom named: O72 within residue: .R<SFR 269>" ... bla bla bla

  Which later generates the following issues because no atom types were set to the newly created atoms:
  "FATAL:  Atom .R<SFR 269>.A<O72 1> does not have a type."... bla bla bla

  Solution: rename atom in the PDB or mol2 to be the same.

2) Atom types in the residue library (e.g. mol2) HAVE TO match those at the parameter file (e.g. frcmod).
 
  You haven't seen this error yet, but you would:
  "For atom: .R<GOL 291>.A<C1 1> Could not find vdW (or other) parameters for type: C.3"

  Solution: It seems gol.mol2 takes parameters from gaff as gol.frcmod is empty, then the former's are the atom types to match.
            In case of sfr.mol2, rename atom types to match sfr.frcmod.

Best,

Matias

PS: There is no magic nor artificial intelligence behind AMBER, yet... so user are expected to provide that sense...

------------------------------------
PhD.
Researcher at Biomolecular Simulations Lab.
Institut Pasteur de Montevideo | Uruguay
[http://pasteur.uy/en/labs/biomolecular-simulations-laboratory]
[http://www.sirahff.com]

----- Mensaje original -----
De: "Hira Jabeen" <hira_bioinfo.yahoo.com>
Para: "AMBER Mailing List" <amber.ambermd.org>
Enviados: MiƩrcoles, 19 de Junio 2019 8:20:19
Asunto: Re: [AMBER] Fw: parameterization of ligand+protein complex

 Hi Dan,
I tried changing the case sensitive thing...changed all GOL and SFR as uppercase now but still facing the same problem. Looking forward
Thanks,Hira
    On Tuesday, June 18, 2019, 05:32:20 PM GMT+1, Daniel Roe <daniel.r.roe.gmail.com> wrote: 
 
 Hi,

I think your problem might be that your ligands GOL and SFR are not
being recognized due to case sensitivity. If you look at the leap
output messages you see:

> > loadoff gol.lib
> Loading library: ./gol.lib
> Loading: gol

But in the PDB you're trying to load the residue name is GOL (all
uppercase). I think it might be enough to add PDB resmap entries prior
to trying to load the PDB, e.g.

addpdbresmap { { "GOL" "gol"} {"SFR" "sfr"} }

However, you may want to correct the names in your library files as a
more long term solution.

Hope this helps,

-Dan

On Tue, Jun 18, 2019 at 12:14 PM Hira Jabeen <hira_bioinfo.yahoo.com> wrote:
>
>  Hi Dan,I am using Amber18 with Ambertool19. The leap and frcmod files are attached here with. Thanks
>    On Tuesday, June 18, 2019, 05:10:52 PM GMT+1, Daniel Roe <daniel.r.roe.gmail.com> wrote:
>
>  Hi,
>
> You should state the version of AmberTools you are using, as well as
> the exact input you are giving to leap.
>
> -Dan
>
> On Tue, Jun 18, 2019 at 6:28 AM Hira Jabeen <hira_bioinfo.yahoo.com> wrote:
> >
> >
> >
> >    ----- Forwarded Message ----- From: Hira Jabeen <hira_bioinfo.yahoo.com>To: amber.ambermd.org <amber.ambermd.org>Sent: Friday, June 14, 2019, 05:09:46 PM GMT+1Subject: parameterization of ligand+protein complex
> >  Hi Amber Team,
> > I am using Amber-18 for the parameterization of ligand-enzyme complex through tleap. I have used pdb files to generate mol2 and then libraries in tleap.I have edited mol2 file according to pdb (replaced UNK with ligand names)Still I get this error:Failed to generate parameters
> > Parameter file was not saved.Protein and ligand files are attached here with.Please advice!
> > Thanks,Hira
> >  _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>  _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
 
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
  
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jun 20 2019 - 04:30:02 PDT
Custom Search