On Thu, Dec 10, 2020, Seibold, Steve Allan wrote:
>I am attempting to build a peptide using tleap with a disulfide bond. I
>have no problem making the peptide; however, when I add the CYX to CYX
>disulfide, I get the error:
>> bond mol.4.SG mol.133.SG
>bond: Argument #1 is type String must be of type: [atom]
>usage: bond <atom1> <atom2> [order]
Use the "desc" command to see what you have, e.g. 'desc mol.4' will
show you what is in residue 4. See if there is an SG atom there.
>
>Yet the usage is correct here: mol = name.pdb and the atom numbers are
>4 and 133 for SG.
The second part of 'mol.4.SG' is the *residue* number, not the atom
number.
The problem maybe that the distance between the CYX
>atom 4 and the CYX atom 133 is long.
No: tleap will create a bond without looking at the distance. You
problem is the "mol.4.SG" does not refer to any atom that tleap knows
about.
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Dec 10 2020 - 10:30:03 PST
