Re: [AMBER] Suggestions on crossing energy barrier to break strong HBs

From: Amit Sharma (Asstt. Prof., MCARS) <"Amit>
Date: Sat, 30 Jan 2021 11:57:43 +0530

I am trying to replicate what is already done in one of the published works
where it is shown through long amber simulation that these bond breaking
happens.

Although following exactly what is already done, I do not get to see the
bond breaking in my case.

Wrote to the authors, but no reply yet since longtime.....I thought why not
I try forcing the system to break the bond and then see if further steps
replicate or not.



On Sat, Jan 30, 2021 at 8:39 AM David A Case <david.case.rutgers.edu> wrote:

> On Sat, Jan 30, 2021, Amit Sharma (Asstt. Prof., MCARS) wrote:
> >
> >I am simulating a protein system in which a heme molecule is strongly
> >hydrogen bonded (via 2 HB's) to one of the protein residues. These HB's
> are
> >expected to break during the course (around ~700 ns) of NPT run and shall
> >lead to the further cascade of bond breaking/making in the protein.
>
> Why do you think the hydrogen bonds should break at some particular time?
> Even if you have a good reason for thinking you know the rate constant,
> in any particular simulation the actual time may be much shorter of much
> longer than the average time. Brute force MD is not a very efficient
> method
> for studying kinetic processes on such long (by MD simulation standards)
> time scales.
>
> And, it could well be that your estimate of the strength of the hydrogen
> bonds in question is too low. Further, if you have to break both hydrogen
> bonds in order to proceed further that could slow things down even more.
>
> >
> > imin = 0, irest = 1, ntx = 7, ntb = 2,
> > pres0 = 1.0, ntp = 1, taup = 0.5, gamma_ln = 1.0,
> > cut = 8.0, ntr = 0, ntc = 2, ntf = 2,
> > tempi = 298.0, temp0 = 298.0, ntt = 3
>
> These look pretty standard, and generally OK. Others on the list might be
> able to suggest some "accelerated" technique that would allow particular
> hydrogen bonds to break more rapidly.
>
> ....good luck....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>

-- 
_
_
*Please consider the environment before printing this email.*
The 
information contained in this electronic message and any *attachments to 
this message are intended for the exclusive use of addressee(s) and may 
contain proprietary, confidential or privileged information. If you are not 
the intended recipient, you should not disseminate, distribute **or copy 
this e-mail. Please notify the
sender immediately and destroy all copies of 
this message and any attachments.***
**
**WARNING***: The recipient should 
check this email and any attachments for
the presence of viruses. JMI 
accepts no liability for any damage caused by any
virus transmitted by this 
email.*
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jan 29 2021 - 22:30:02 PST
Custom Search