[gmx-users] During mdrun my protein moves towards edges

Dear gromacs users I have been running the real mdrun after equilibration. mdrun has been submitted for 500ns. After some time(150ns) the protein moves towards edges. I want the internal moments only, but whole protein moves all sides of box rather than at centre of the box. Kindly tell me how to

Re: [gmx-users] During mdrun my protein moves towards edges

You do not have a problem you only have the normal protein diffusion. Take a look in trjconv tool. http://manual.gromacs.org/programs/gmx-trjconv.html Option -center will be your best friend I guess. Diogo Vila-Viçosa Em seg, 11 de mai de 2015 às 10:14, Seera Suryanarayana

Re: [gmx-users] Subject: Re: how to avoid multiple 'aprun' on batch job script?

Hi, You should check out the other recent discussion in this list on performance variation. Getting your GROMACS nodes allocated close together is an important part of mitigating such problems. Rather than manually splitting your job into small pieces, you can have mdrun do that automatically

Re: [gmx-users] Why h_bond can not form

On 5/10/15 8:48 PM, 范聪 wrote: Hello everyone! I offered a picture called ask.png to illistrate my question. According to the referense papers, the H atom of the -OH group of TYR306 should form H_bond with the carbonyl O atom of the ligand. Depending on the methods used in those papers,

Re: [gmx-users] Sanjay (Re: Energy minimization problem)

On 5/10/15 3:28 PM, sanjay choubey wrote: HI,I am running membrane protein dynamics. After energy minimization i am getting the following error. I am not getting how to rectify it.Energy is minimized in fewer steps. Without proper minimization i cant proceed to next step of NVT run.

Re: [gmx-users] Why h_bond can not form

It's so kind of you to answer my question in such a short time. I'll try other forcefields. Thank you very much! 在 2015-05-11 19:46:27，Justin Lemkul jalem...@vt.edu 写道： On 5/10/15 8:48 PM, 范聪 wrote: Hello everyone! I offered a picture called ask.png to illistrate my question. According

Re: [gmx-users] HYL force field

On 5/11/15 8:19 AM, Ming Tang wrote: Dear Justin, The protein contains 15 constituent HYL residue modified from LYS. I realized that I need to modify the .rtp, but have no idea about how to derive the parameters correctly. Do you have any suggestions for this? As I said, it depends on

Re: [gmx-users] first residue in chains warning issue

On 5/10/15 11:36 PM, Ming Tang wrote: Dear all, I built a triple helix with 2 GLN and 1 SER being the first residue of the three chains. However, while using pdb2gmx -f 1.pdb -o complex.gro -ignh -ter -ff gromos54a7 -water SPC and choosing NH2 or NH3+ as the start terminus type for GLN-1

Re: [gmx-users] HYL force field

Dear Justin, The protein contains 15 constituent HYL residue modified from LYS. I realized that I need to modify the .rtp, but have no idea about how to derive the parameters correctly. Do you have any suggestions for this? thanks a lot. From:

Re: [gmx-users] Harmonic COM Restraints with the Pull Code

On 5/9/15 9:52 PM, Eric Dybeck wrote: Hello All, I am looking to run a set of simulations in which molecules are slowly restrained to their COM positions. To provide some context, I am going to melt a solid, and I would like to restrain the molecules first so that they remain in their lattice

[gmx-users] implicit solvent concerns

Hi All, I am thinking of running an implicit solvent simulation using Gromacs (.mdp file pasted at the end of the email). However, I am a bit concerned about its accuracy as I get different results for the Polarisation term. More specifically, the first step of an energy minimisation gives

Re: [gmx-users] Subject: Re: how to avoid multiple 'aprun' on batch job script?

Dear Dr. Mark Abraham, Thank you for all the suggestions. I have a query now about the 'mdrun' restart command line: First run: export OMP_NUM_THREADS=2 aprun -n 144 -N 24 -j1 mdrun_mpi -npme 48 -maxh 24 -v -deffnm simu-100ns (simu-100ns.tpr for whole 100ns) Restart-1,2,3...etc: export

Re: [gmx-users] Gromacs 4.6.7 with MPI and OpenMP

Szilárd, On Friday 08 May 2015 21:18:12 Szilárd Páll wrote: What is your goal with using CPUSETs? Node sharing? Correct. While it might be possible to see the cores that have been assigned to the job and do the correct 'pin setting' it would probably be ugly. Not sure what you

Re: [gmx-users] Gromacs 4.6.7 with MPI and OpenMP

Mark, On Friday 08 May 2015 15:15:31 Mark Abraham wrote: FWIW, I ran the same GROMACs run outside of the queuing system to verify that the CPUSETs were not causing the issue. MPI gets a chance to play with OMP_NUM_THREADS (and pinning!), too, so your tests suggest the issue lies there.

Re: [gmx-users] implicit solvent concerns

Hi, Stuff was probably broken here before 4.6. I would trust only single threaded in 4.6 and later, and only after comparing with 4.5. Mark On Mon, 11 May 2015 14:50 Albert Solernou a.soler...@leeds.ac.uk wrote: Hi All, I am thinking of running an implicit solvent simulation using Gromacs

Re: [gmx-users] Subject: Re: how to avoid multiple 'aprun' on batch job script?

Hi, Please see http://www.gromacs.org/Documentation/How-tos/Extending_Simulations and try things out. We can't bless everybody's scripts... Mark On Mon, 11 May 2015 15:30 Satyabrata Das satyabrata...@gmail.com wrote: Dear Dr. Mark Abraham, Thank you for all the suggestions. I have a query

Re: [gmx-users] implicit solvent concerns

Thanks for being so clear Mark. Best, Albert On 05/11/2015 03:18 PM, Mark Abraham wrote: Hi, Stuff was probably broken here before 4.6. I would trust only single threaded in 4.6 and later, and only after comparing with 4.5. Mark On Mon, 11 May 2015 14:50 Albert Solernou

Re: [gmx-users] Giving an initial velocity to a group of atoms

Thanks again for your reply Dr. Lemkul. I was also expecting a disaster because I gave an unrealistic starting velocity of the structure. But the simulation was running fine. I tried using the gmxdump and found that the velocities I am getting are not the velocities I've given to the atoms. In

[gmx-users] LINCS warnings and ligand topology

Dear Users, I prepared ligand topologies by two different methods, one with Acpype and the other with Prodrg. It is sad that independent simulations give the same (or similar) numerous Lincs warnings and the system blew up at the end. Errors look like LINCS WARNING relative constraint

[gmx-users] question

Dear gromacs users I have experience in MD simulation of protein and peptides by gromos force field. Now I am working on a new project in which I should connect two peptides (based on experimental data) to obtain one new peptide. This connection is as follows: CD atom of GLU residue from

Re: [gmx-users] question

In fact, I want to make an isopeptide bond. On Mon, May 11, 2015 at 11:04 PM, Andrew Bostick andrew.bosti...@gmail.com wrote: Dear gromacs users I have experience in MD simulation of protein and peptides by gromos force field. Now I am working on a new project in which I should connect two

Re: [gmx-users] [PyMOL] how to connect to atoms?

Hi Leila, This too is more a question for the Gromacs mailing list. Gromacs does not read bond information (CONECT records) from a PDB file. It only reads the ATOM/HETATM records and sets bonds based on the force field building block definitions and the special bonds (such as the isopeptide link)

[gmx-users] Dihedral copying from itp to ff

Dear User, I have done MD simulation successfully with ligand.itp file produced by Antechamber. Now I have to include the ff information of the ligand into AMBER99SB forcefield and in the process I am so for successful in creating rtp entry and bond and angles parameter in ffbonded.itp file but I

[gmx-users] question

Dear Justin Very thanks for your answer I should use my pdb file in which there is no bond between CD atom of GLU residue from peptide 1 and NZ atom of LYS residue from peptide 2. Creating new entry in specbond.dat file will consider this new bond. Is my vision true? Based on your answer, only

Re: [gmx-users] all-atom pca

Hi Jio, No, the 'extreme structures' are filtered along the eigenvector and are not real structures. But you wanted the corresponding all-atom ones. Cheers, Tsjerk On May 12, 2015 00:45, gromacs query gromacsqu...@gmail.com wrote: Hi Tsjerk In -proj plot of v1 vs Time I used extreme values