[gmx-users] MD at different pH
jalemkul at vt.edu
Thu May 4 14:29:28 CEST 2017
On 5/4/17 1:02 AM, Anu George wrote:
> Through the simulation, I am trying to find out if the binding free energy of the drug-protein complex is strong enough for the drug to be a good inhibitor of the protein.
> so, does the pH of the simulation matter as in reality the active monomer exists in acidic pH conditions?
You have to model whatever the realistic conditions are, whether that's in vivo
or in vitro. This is part of good experimental design. If your target exists
in a mildly acidic microenvironment, that's what you should model. Binding
thermodynamics do depend on protonation states (in some cases, quite strongly)
so if you model one state and try to relate it to data that exist in another,
you're comparing apples to oranges and you've wasted a lot of time and resources.
Also, dispel with the notion of "pH of the simulation" because such a thing does
not exist. You'll be modeling a constant protonation state and there are no
dissociable protons floating around. You're going to do a simulation using the
dominant protonation state of a given set of molecules at a corresponding
> Anu George
> ----- Original Message -----
> From: "Justin Lemkul" <jalemkul at vt.edu>
> To: gmx-users at gromacs.org
> Sent: Wednesday, May 3, 2017 8:15:33 PM
> Subject: Re: [gmx-users] MD at different pH
> On 5/3/17 10:04 AM, Anu George wrote:
>> Thanks Justin for the information.
>> But I had a further question regarding the simulation-if the active monomer exists only in the acidic pH, then will carrying out a simulation of the monomer with a drug at neutral pH help?
> That depends on what your scientific goals/questions are and what you are
> actually trying to model.
>> Anu George
>> ----- Original Message -----
>> From: "Justin Lemkul" <jalemkul at vt.edu>
>> To: gmx-users at gromacs.org
>> Sent: Wednesday, May 3, 2017 5:37:33 PM
>> Subject: Re: [gmx-users] MD at different pH
>> On 5/2/17 11:31 PM, Anu George wrote:
>>> Dear Gromacs users,
>>> I am working on a tetramer(human beta-2 tryptase) protein which exist as an active monomer at an acidic pH of about 6..If I have to do a molecular dynamic simulation of the monomer with a small molecule, should the simulation be carried out at this pH ? what are the steps I should follow?
>> Calculate the pKa values for titratable residues and select their dominant
>> protonation states at that pH when running pdb2gmx.
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
More information about the gromacs.org_gmx-users