[gmx-users] co-ordinate number mismatch
Ankita naithani
ankitanaithani at gmail.com
Fri Oct 12 21:40:26 CEST 2012
Hi,
I am trying to simulate a proein-ligand simulation. My protein is a
tetramer and has 3 ligands, ATP, OXL and FDP. As I am using
amber99sb-ildn forcefield, I generated the topology and co-ordiante
files fusing Acpype. I have appended the information to cinluded in
the topology file i.e. complex.top and also the atomtypes information
in tandem to the acpype instructions. Also, I appended the
co-ordinates of the ligands into my file generated after pdb2gmx step
where the topology of protein was generated. However, after defining
my box, when I preprocess my system before adding ions, I come across
the final error.
Fatal error:
number of coordinates in coordinate file (system_solv.gro, 315767)
does not match topology (complex.top, 315771)
I do understand that somewhere either the information is not updated,
which I have seen several times to be getting updated or else there is
a problem with SOL molecules. I have tried running a simulation of
another protein-ligand system in exactly the same manner and didn't
encounter any problem. It would be really helpful if someone possibly
give me the suggestion and possible solution. I have appended partial
contents of both the files below:
complex.top file contents:
[ molecules ]
; Compound #mols
Protein_chain_A 1
Ion_chain_A2 1
Other_chain_A3 1
Ion_chain_A4 1
Protein_chain_B 1
Ion_chain_B2 1
Other_chain_B3 1
Ion_chain_B4 1
Protein_chain_C 1
Ion_chain_C2 1
Other_chain_C3 1
Ion_chain_C4 1
Protein_chain_D 1
Ion_chain_D2 1
Other_chain_D3 1
Ion_chain_D4 1
SOL 1
SOL 1
SOL 165
SOL 1
SOL 1
SOL 377
SOL 1
SOL 1
SOL 303
SOL 1
SOL 1
SOL 398
FDP_A 1
FDP_B 1
FDP_C 1
FDP_D 1
ATP1 1
ATP2 1
ATP3 1
ATP4 1
OA_m 1
OB_m 1
OC_m 1
OD_m 1
SOL 92850
contents from system_solv.gro:
92849SOL HW215761 15.757 16.224 11.514
92850SOL OW15762 16.249 14.948 11.340
92850SOL HW115763 16.284 14.907 11.256
92850SOL HW215764 16.159 14.910 11.360
92851SOL OW15765 14.935 15.974 11.472
92851SOL HW115766 15.034 15.963 11.463
92851SOL HW215767 14.895 15.888 11.504
16.39712 16.39712 11.59451 0.00000 0.00000 0.00000 0.00000
8.19856 8.19856
--
Ankita Naithani
More information about the gromacs.org_gmx-users
mailing list