From: Josh Vermaas (joshua.vermaas_at_gmail.com)
Date: Wed Mar 04 2020 - 09:22:20 CST
Ok, fix your psf so that they are bonded. How are you currently making your
psf? Using extrabonds or colvars for this isn't the right approach, but
instead you should be focusing on why there isn't a bond where you think
there should be one.
-Josh
On Wed, Mar 4, 2020 at 10:16 AM soodabeh ghaffari <
soodabehghaffari_at_gmail.com> wrote:
> Hello Dr Fiorin,
> These two atoms ( backbone N of one residue and carbonyl carbon of the
> other) should be covalently bound but in my pdb file they are not. Is
> extraBonds fine(the first link you provided)?Do you have any example of
> usage of extraBonds?
>
> 1. If so, how I can get the parameters k, ref, n. And how can I
> select atom?
>
> 2.Is this right to add these line to MD.conf?
> *extraBonds*
> *extraBondsCosAngles*
> *extraBondsFile myfile.txt*
>
> 3. In myfile.txt,is <atom> atom serial number in pdb file? and <ref> is
> the bond length?
> bond <atom> <atom> <k> <ref>
> For instance, is it right to write below(if atom serial number is 100 and
> and 101 and k=250 and bond length should be 1.3 (but currently, it is more
> than 1.3))?
> bond 100 101 250 1.3
>
> 4.and one more question, the location of myfile.txt should be the same as
> conf file?
>
>
> I appreciate your help
> Thank you so much
> Best Regards,
> Soodabeh
>
>
> On Wed, Mar 4, 2020 at 9:04 AM Giacomo Fiorin <giacomo.fiorin_at_gmail.com>
> wrote:
>
>> Hello Soodabeh, if these are not meant to be covalently bound (are
>> they??), the best option to create a new "peptide bond" sounds like
>> extraBonds:
>>
>> http://www.ks.uiuc.edu/Research/namd/2.13/ug/node27.html#SECTION00086400000000000000
>>
>> Other more flexible but also more complex options are Colvars or
>> tclForces:
>> http://www.ks.uiuc.edu/Research/namd/2.13/ug/node51.html
>> http://www.ks.uiuc.edu/Research/namd/2.13/ug/node48.html
>>
>> On Wed, Mar 4, 2020 at 12:31 AM soodabeh ghaffari <
>> soodabehghaffari_at_gmail.com> wrote:
>>
>>> Dear Sir/Madam,
>>> I have a pdb file in which two residues are not connected and their
>>> distances are greater than peptide bond. I want to do MD simulation and add
>>> constraint to MD.conf so that their distance should be with peptide bond
>>> and I have a connected pdb file at the end.
>>> Can you help me? which constraint should I have use and how?
>>> Do you have any example?
>>> Thanks
>>> Best Regards,
>>> Soodabeh
>>>
>>
>>
>> --
>> Giacomo Fiorin
>> Associate Professor of Research, Temple University, Philadelphia, PA
>> Research collaborator, National Institutes of Health, Bethesda, MD
>> http://goo.gl/Q3TBQU
>> https://github.com/giacomofiorin
>>
>
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:08 CST