Normally yes there is a TER between the both chains in the initial pdb
as well as a TER in the end.. :-)
here is the tleap log
/teLeap: Warning!
Converting N-terminal residue name to PDB format: NGLY -> GLY
/teLeap: Warning!
Converting C-terminal residue name to PDB format: CGLY -> GLY
/teLeap: Warning!
Converting N-terminal residue name to PDB format: NASP -> ASP
/teLeap: Warning!
Converting C-terminal residue name to PDB format: CARG -> ARG
Checking Unit.
Building topology.
Building atom parameters.
Building bond parameters.
Building angle parameters.
Building proper torsion parameters.
Building improper torsion parameters.
total 672 improper torsions applied
Building H-Bond parameters.
Incorporating Non-Bonded adjustments.
Not Marking per-residue atom chain types.
Marking per-residue atom chain types.
(Residues lacking connect0/connect1 -
these don't have chain types marked:
res total affected
CARG 1
CGLY 1
NASP 1
NGLY 1
WAT 37077
mb this is just a visual effect ?
Il giorno mer 5 lug 2023 alle ore 11:41 Steinbrecher, Thomas
<thomas.steinbrecher.roche.com> ha scritto:
>
> Hi Enrico,
>
> does your input PDB contain the right "TER" lines between separate protein chains? Otherwise, leap may consider this as one continuous protein chain but should warn you about any "long bonds" that are created.
>
> Kind Regards,
>
> Thomas
>
> On Wed, Jul 5, 2023 at 11:34 AM Enrico Martinez via AMBER <amber.ambermd.org> wrote:
>>
>> Dear Amber users !
>>
>> I am preparing system of water soluble protein using tleap with the
>> following script :
>>
>> source leaprc.protein.ff19SB
>> source leaprc.water.opc
>> prot = loadpdb input.pdb
>> addIons prot Cl- 0
>> addIons prot Na+ 0
>> setbox prot centers
>> solvateoct prot TIP3PBOX 20.0
>> savepdb output.pdb
>> saveamberparm output.prmtop output.prmcrd
>> quit
>>
>>
>> Briefly, everything works OK but in the output.pdb I see that the N
>> and C terminus of the proteins are connected by the artificial
>> covalent bond despite the large physical distance between them... I
>> performed a short MD simulation for the system and extracted snapshots
>> after it and in pymol I see the same thing.
>>
>> Should the C and N parts be capped to prevent this bug ? If yes, could
>> you please suggest the modification of my tleap script in order that
>> the first and the last residues would be automatically capped without
>> explicit definition of the amino acid ?
>>
>> Many thanks in advance
>>
>> Cheers
>>
>> Enrico
>>
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>
>
>
> --
> Dr. Thomas Steinbrecher
> Principal Scientist CADD
>
> Roche Pharma Research and Early Development
> Roche Innovation Center Basel
> F. Hoffmann-La Roche Ltd
> Bldg. 092/3.92
> Grenzacherstrasse 124
> 4070 Basel
> Switzerland
>
> Phone +41 61 682 1319
> mailto: thomas.steinbrecher.roche.com
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jul 05 2023 - 03:30:02 PDT