29 September 2023 0 4K Report

I hope this message finds you well. I am currently working on a research project focused on MD simulation study of drug desorption behavior on a boron nitride nanosheet with a non-covalent polymer in gromacs . As part of my research, I am facing challenges with the protonation procedure,

Here is a brief overview of the procedure I have been following:

  • I generated an itp file for the protonated drug and polymer using a web server.
  • I converted the gro file from the last molecular dynamics (MD) run and created an index file containing only the drug polymer and sheets. I then converted this to a pdb file.
  • To protonate the system, I added hydrogen atoms to the drug and polymer using GaussView and saved the structure as a pdb file. Subsequently, I converted the pdb file back to a gro file.
  • I used an 'include' statement for the protonated drug and polymer itp files in the previous topology.
  • However, I have encountered issues with this procedure. The conversions seem to alter atom types, leading to discrepancies between the atom types in the topology (itp) and gro files, resulting in errors during the simulation.

    I would greatly appreciate it if anyone with expertise in this area could share their knowledge and suggest the correct procedure for protonating the system and running MD simulations. Additionally, if you have recommendations for software tools that you have successfully used for drug protonation in similar research, I would be eager to hear your suggestions.

    More Aqsa Kamran's questions See All
    Similar questions and discussions