[gmx-users] Problem with simulation of Protein-DNA complex
Dear Gromacs users, I have been trying to simulate a Protein-DNA complex using gromacs, but each time I have been facing problems. I would like to point out that both the Protein and DNA have been modeled and after that docked in order to obtain a complex. Following are the parameters I am using: Force-Field: amber99sb-ildn water model: TIP3P 40 NA ions were added in order to neutralise the complex-solvent system. Gromacs doesn't show any error up until I proceed to do nvt equilibriation. Following is nvt.mdp: title = PFV_DNA_NVT define = -DPOSRES ; position restrain the protein ; Run parameters integrator = md; leap-frog integrator nsteps = 100 ; 2 * 5 = 100 ps dt = 0.001 ; 2 fs ; Output control nstxout = 100 ; save coordinates every 0.2 ps nstvout = 100 ; save velocities every 0.2 ps nstenergy = 100 ; save energies every 0.2 ps nstlog = 100 ; update log file every 0.2 ps ; Bond parameters continuation= no; first dynamics run constraint_algorithm = lincs; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracyi ; Neighborsearching ns_type = grid ; search neighboring grid cels nstlist = 5 ; 10 fs rlist = 1.0 ; short-range neighborlist cutoff (in nm) rcoulomb= 1.0 ; short-range electrostatic cutoff (in nm) rvdw= 1.0 ; short-range van der Waals cutoff (in nm) ; Electrostatics coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling is on tcoupl = V-rescale ; modified Berendsen thermostat tc-grps = Protein Non-Protein tau_t = 0.1 0.1 ; time constant, in ps ref_t = 300 300 ; reference temperature, one for each group, in K ; Pressure coupling is off pcoupl = no; no pressure coupling in NVT ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr= EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = yes ; assign velocities from Maxwell distribution gen_temp= 300 ; temperature for Maxwell distribution gen_seed= -1; generate a random seed It generates a tpr, but upon mpirun gives the following error: starting mdrun 'Protein in water' 100 steps, 1000.0 ps. Step 0, time 0 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.027514, max 0.979682 (between atoms 3460 and 3461) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3455 3457 65.20.1542 0.2535 0.1526 3455 3456 33.30.1093 0.1322 0.1090 3457 3459 72.70.1105 0.1739 0.1090 3457 3458 65.70.1105 0.1588 0.1090 3460 3469 78.60.1409 0.2693 0.1400 3460 3461 78.80.1409 0.2772 0.1400 3461 3463 75.10.1406 0.1489 0.1400 3463 3465 37.20.1403 0.1756 0.1400 3463 3464 40.70.1083 0.1408 0.1080 3465 3467 35.50.1403 0.1733 0.1400 3467 3469 74.00.1405 0.1373 0.1400 3467 3468 42.90.1083 0.1433 0.1080 Wrote pdb files with previous and current coordinates Step 1, time 0.001 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 364940.065742, max 11977066.00 (between atoms 3465 and 3466) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3475 3495 110.80.1526 20544.7148 0.1522 3475 3477 111.60.1529 20544.4336 0.1526 3475 3476 108.10.1094 20543.1309 0.1090 3477 3480 62.80.1527 4.6441 0.1526 3477 3479 57.50.1091 4.5683 0.1090 3477 3478 93.60.1091 4.6160 0.1090 3495 3497 138.10.1336 4.9133 0.1335 3495 3496 69.60.1230 4.9085 0.1229 3497 3499 40.20.1449 0.2036 0.1449 3497 3498 43.80.1010 0.1527 0.1010 3480 3483 35.60.1526 0.1902 0.1526 3480 3482 42.20.1090 0.1501 0.1090 3480 3481 41.70.1090 0.1485 0.1090 3453 3455 88.70.1701 117322.3672 0.1449 3453 3454 88.30.1035 65066.2227 0.1010 3455 3471 48.90.1763 171383.6562 0.1522 3455 3457 134.60.2535 572521.2500 0.1526 3455 3456 66.70.1322 143042. 0.1090 3457 3460 67.5
Re: [gmx-users] Problem with simulation of Protein-DNA complex
On Thu, Feb 9, 2012 at 8:44 PM, wrote: > Dear Gromacs users, > > I have been trying to simulate a Protein-DNA complex using gromacs, but each > time I have been facing problems. I would like to point out that both the > Protein and DNA have been modeled and after that docked in order to obtain a > complex. Following are the parameters I am using: > > Force-Field: amber99sb-ildn > water model: TIP3P > > 40 NA ions were added in order to neutralise the complex-solvent system. > Gromacs doesn't show any error up until I proceed to do nvt equilibriation. Before the NVT, you may do some energy minimization. > Following is nvt.mdp: > > title = PFV_DNA_NVT > define = -DPOSRES ; position restrain the protein > ; Run parameters > integrator = md ; leap-frog integrator > nsteps = 100 ; 2 * 5 = 100 ps > dt = 0.001 ; 2 fs > ; Output control > nstxout = 100 ; save coordinates every 0.2 ps > nstvout = 100 ; save velocities every 0.2 ps > nstenergy = 100 ; save energies every 0.2 ps > nstlog = 100 ; update log file every 0.2 ps > ; Bond parameters > continuation = no ; first dynamics run > constraint_algorithm = lincs ; holonomic constraints > constraints = all-bonds ; all bonds (even heavy atom-H bonds) > constrained > lincs_iter = 1 ; accuracy of LINCS > lincs_order = 4 ; also related to accuracyi > ; Neighborsearching > ns_type = grid ; search neighboring grid cels > nstlist = 5 ; 10 fs > rlist = 1.0 ; short-range neighborlist cutoff (in nm) > rcoulomb = 1.0 ; short-range electrostatic cutoff (in nm) > rvdw = 1.0 ; short-range van der Waals cutoff (in nm) > ; Electrostatics > coulombtype = PME ; Particle Mesh Ewald for long-range > electrostatics > pme_order = 4 ; cubic interpolation > fourierspacing = 0.16 ; grid spacing for FFT > ; Temperature coupling is on > tcoupl = V-rescale ; modified Berendsen thermostat > tc-grps = Protein Non-Protein > tau_t = 0.1 0.1 ; time constant, in ps > ref_t = 300 300 ; reference temperature, one for each group, > in K > ; Pressure coupling is off > pcoupl = no ; no pressure coupling in NVT > ; Periodic boundary conditions > pbc = xyz ; 3-D PBC > ; Dispersion correction > DispCorr = EnerPres ; account for cut-off vdW scheme > ; Velocity generation > gen_vel = yes ; assign velocities from Maxwell > distribution > gen_temp = 300 ; temperature for Maxwell distribution > gen_seed = -1 ; generate a random seed > > > It generates a tpr, but upon mpirun gives the following error: > > starting mdrun 'Protein in water' > 100 steps, 1000.0 ps. > > Step 0, time 0 (ps) LINCS WARNING > relative constraint deviation after LINCS: > rms 0.027514, max 0.979682 (between atoms 3460 and 3461) > bonds that rotated more than 30 degrees: > atom 1 atom 2 angle previous, current, constraint length > 3455 3457 65.2 0.1542 0.2535 0.1526 > 3455 3456 33.3 0.1093 0.1322 0.1090 > 3457 3459 72.7 0.1105 0.1739 0.1090 > 3457 3458 65.7 0.1105 0.1588 0.1090 > 3460 3469 78.6 0.1409 0.2693 0.1400 > 3460 3461 78.8 0.1409 0.2772 0.1400 > 3461 3463 75.1 0.1406 0.1489 0.1400 > 3463 3465 37.2 0.1403 0.1756 0.1400 > 3463 3464 40.7 0.1083 0.1408 0.1080 > 3465 3467 35.5 0.1403 0.1733 0.1400 > 3467 3469 74.0 0.1405 0.1373 0.1400 > 3467 3468 42.9 0.1083 0.1433 0.1080 > Wrote pdb files with previous and current coordinates > > Step 1, time 0.001 (ps) LINCS WARNING > relative constraint deviation after LINCS: > rms 364940.065742, max 11977066.00 (between atoms 3465 and 3466) > bonds that rotated more than 30 degrees: > atom 1 atom 2 angle previous, current, constraint length > 3475 3495 110.8 0.1526 20544.7148 0.1522 > 3475 3477 111.6 0.1529 20544.4336 0.1526 > 3475 3476 108.1 0.1094 20543.1309 0.1090 > 3477 3480 62.8 0.1527 4.6441 0.1526 > 3477 3479 57.5 0.1091 4.5683 0.1090 > 3477 3478 93.6 0.1091 4.6160 0.1090 > 3495 3497 138.1 0.1336 4.9133 0.1335 > 3495 3496 69.6 0.1230 4.9085 0.1229 > 3497 3499 40.2 0.1449 0.2036 0.1449 > 3497 3498 43.8 0.1010 0.1527 0.1010 > 3480 3483 35.6 0.1526 0.1902 0.1526 > 3480 3482 42.2 0.1090 0.1501 0.1090 > 3480 3481 41.7 0.1090 0.1485 0.1090 > 3453 3455 88.7 0.1701 117322.3672 0.14
Re: [gmx-users] Problem with simulation of Protein-DNA complex
Hi Lina, I am sorry, I think I forgot to mention that I did perform energy minimisation using Steep-Descent for 5000 steps, before NVT. I was so engrossed in other details that I forgot to mention it! Quoting lina : On Thu, Feb 9, 2012 at 8:44 PM, wrote: Dear Gromacs users, I have been trying to simulate a Protein-DNA complex using gromacs, but each time I have been facing problems. I would like to point out that both the Protein and DNA have been modeled and after that docked in order to obtain a complex. Following are the parameters I am using: Force-Field: amber99sb-ildn water model: TIP3P 40 NA ions were added in order to neutralise the complex-solvent system. Gromacs doesn't show any error up until I proceed to do nvt equilibriation. Before the NVT, you may do some energy minimization. Following is nvt.mdp: title = PFV_DNA_NVT define = -DPOSRES ; position restrain the protein ; Run parameters integrator = md ; leap-frog integrator nsteps = 100 ; 2 * 5 = 100 ps dt = 0.001 ; 2 fs ; Output control nstxout = 100 ; save coordinates every 0.2 ps nstvout = 100 ; save velocities every 0.2 ps nstenergy = 100 ; save energies every 0.2 ps nstlog = 100 ; update log file every 0.2 ps ; Bond parameters continuation = no ; first dynamics run constraint_algorithm = lincs ; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracyi ; Neighborsearching ns_type = grid ; search neighboring grid cels nstlist = 5 ; 10 fs rlist = 1.0 ; short-range neighborlist cutoff (in nm) rcoulomb = 1.0 ; short-range electrostatic cutoff (in nm) rvdw = 1.0 ; short-range van der Waals cutoff (in nm) ; Electrostatics coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling is on tcoupl = V-rescale ; modified Berendsen thermostat tc-grps = Protein Non-Protein tau_t = 0.1 0.1 ; time constant, in ps ref_t = 300 300 ; reference temperature, one for each group, in K ; Pressure coupling is off pcoupl = no ; no pressure coupling in NVT ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr = EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = yes ; assign velocities from Maxwell distribution gen_temp = 300 ; temperature for Maxwell distribution gen_seed = -1 ; generate a random seed It generates a tpr, but upon mpirun gives the following error: starting mdrun 'Protein in water' 100 steps, 1000.0 ps. Step 0, time 0 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.027514, max 0.979682 (between atoms 3460 and 3461) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3455 3457 65.2 0.1542 0.2535 0.1526 3455 3456 33.3 0.1093 0.1322 0.1090 3457 3459 72.7 0.1105 0.1739 0.1090 3457 3458 65.7 0.1105 0.1588 0.1090 3460 3469 78.6 0.1409 0.2693 0.1400 3460 3461 78.8 0.1409 0.2772 0.1400 3461 3463 75.1 0.1406 0.1489 0.1400 3463 3465 37.2 0.1403 0.1756 0.1400 3463 3464 40.7 0.1083 0.1408 0.1080 3465 3467 35.5 0.1403 0.1733 0.1400 3467 3469 74.0 0.1405 0.1373 0.1400 3467 3468 42.9 0.1083 0.1433 0.1080 Wrote pdb files with previous and current coordinates Step 1, time 0.001 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 364940.065742, max 11977066.00 (between atoms 3465 and 3466) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3475 3495 110.8 0.1526 20544.7148 0.1522 3475 3477 111.6 0.1529 20544.4336 0.1526 3475 3476 108.1 0.1094 20543.1309 0.1090 3477 3480 62.8 0.1527 4.6441 0.1526 3477 3479 57.5 0.1091 4.5683 0.1090 3477 3478 93.6 0.1091 4.6160 0.1090 3495 3497 138.1 0.1336 4.9133 0.1335 3495 3496 69.6 0.1230 4.9085 0.1229 3497 3499 40.2 0.1449 0.2036 0.1449 3497 3498 43.8 0.1010 0.1527 0.1010 3480 3483 35.6 0.1526 0.1902 0.1526 3480 3482 42.2 0.1090 0.1501 0.1090 3480 3481 41.7 0.1090 0.1485 0.1090 3453 3455 88.7
Re: [gmx-users] Problem with simulation of Protein-DNA complex
Hi Rohit, Have you checked the atoms around 3460? It seems there's something wrong with your starting structure. Cheers, Tsjerk On Feb 9, 2012 2:09 PM, wrote: Hi Lina, I am sorry, I think I forgot to mention that I did perform energy minimisation using Steep-Descent for 5000 steps, before NVT. I was so engrossed in other details that I forgot to mention it! Quoting lina : > On Thu, Feb 9, 2012 at 8:44 PM, http://lists.gromacs.org/mailman/listinfo/gmx-... -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] Problem with simulation of Protein-DNA complex
Hi Tsjerk, I have checked and re-checked the structure and it seems fine. I also made sure to minimise the starting structure really well after pdb2gmx. Although it did not achive required precision, but machine precision instead. Also, the problem starts with the NVT, before that everythings seems to work fine. Could there be any other source of error? I am not able to figure it out. Many Thanks. Quoting Tsjerk Wassenaar : Hi Rohit, Have you checked the atoms around 3460? It seems there's something wrong with your starting structure. Cheers, Tsjerk On Feb 9, 2012 2:09 PM, wrote: Hi Lina, I am sorry, I think I forgot to mention that I did perform energy minimisation using Steep-Descent for 5000 steps, before NVT. I was so engrossed in other details that I forgot to mention it! Quoting lina : > On Thu, Feb 9, 2012 at 8:44 PM, http://lists.gromacs.org/mailman/listinfo/gmx-... Rohit Arora Laboratoire de Biologie et de Pharmacologie Génétique Appliquée (CNRS UMR 8113) Ecole Normale Supérieure, Cachan France Bureau: +33 (0) 1 47 40 77 49 Portable: +33 (0) 6 23 85 12 46 -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists