Re: [AMBER] Combining CpHMD and aMD

From: Cruzeiro,Vinicius Wilian D <>
Date: Mon, 28 Jan 2019 12:36:28 +0000

Hello Cecilia,

I (highly) suspect that your issue is coming from the solvent relaxation stage. In the constant pH code, once a protonation state change is accepted some MD steps are executed where the solute degrees of freedom are frozen. This stage could me messing with the aMD actions. If my hypothesis is correct, running without solvent relaxation (ntrelax=0) should be fine and it would explain why CpHMD+aMD works in Implicit solvent (as you mentioned in your first email). As you shouldn’t run CpHMD without solvent relaxation, some coding adjustments are necessary.

All the best,

Vinícius Wilian D Cruzeiro

PhD Candidate
Department of Chemistry, Physical Chemistry Division
University of Florida, United States

Voice: +1(352)846-1633<tel:+1(352)846-1633>

On Jan 28, 2019, at 2:27 AM, Cecilia Lindgren <<>> wrote:

Hello Vinícius,

Thank you for your reply. I have tried to run it with only aMD as well, and then it looks ok. So it is only when I combine the two that I experience this odd behavior.

Best regards,

Cecilia Lindgren, Ph.D.
Umeå University
Department of Chemistry
Linnaeus väg
SE-901 87 Umeå

-----Original Message-----
From: Cruzeiro,Vinicius Wilian D <<>>
Sent: den 25 januari 2019 23:00
To: '<>' <<>>
Subject: Re: [AMBER] Combining CpHMD and aMD

Hello Cecilia,

Combining different techniques in Amber is something not necessarily trivial. I am aware of some cases in which additional coding would be required for the combination of techniques to work properly. However, I have never attempted performing CpHMD with aMD, thus I do not know before hand of any eventual issues that would only raise when both are executed together.

With that in mind, I would suggest you to run your simulation that fail again but with no CpHMD (icnstph=0). This would tell you to if your issue happens with aMD or only when CpHMD and aMD are executed together.

I hope this helps,

All the best,

Vinícius Wilian D. Cruzeiro

PhD Candidate
Department of Chemistry, Physical Chemistry Division University of Florida, United States

Voice: +1(352)846-1633

From: Cecilia Lindgren <<>>
Sent: Thursday, January 24, 2019 10:37:26 AM
To: '<>'
Subject: [AMBER] Combining CpHMD and aMD

Dear Amber Users,

I am trying to run constant pH MD (CpHMD) with AMBER, in combination with accelerated MD (aMD). I am using explicit water.
However, when I use these settings the protein unfolds, as if it is being "pulled apart". I have attached a figure from the beginning of the MD simulation and from when the protein starts to unfold. This simulation included 10 ns of CpHMD, followed by 10 ns of CpHMD in combination with aMD. The unfolding starts almost immediately after aMD is started.
I have restarted the simulation twice but the problem remains.

If I use CpHMD without trying to combine it with aMD, it looks ok during the 70 ns I have run it for. (However, it is clearly more dynamic than when running normal MD.)

Has anyone else tried this? I can only find examples of implicit MD where CpHMD and aMD are combined. Is it possible, or do I need to use implicit water instead of explicit in case I want to combine CpHMD with aMD?

I have also attached the settings I use for the production run.

Best regards,

AMBER mailing list<>

AMBER mailing list<>
AMBER mailing list
Received on Mon Jan 28 2019 - 05:00:02 PST
Custom Search