Re: [AMBER] Peptide's charge is not 0 when getting parm7 and rst7 files

From: Carlos Simmerling <carlos.simmerling.gmail.com>
Date: Mon, 13 Dec 2021 17:15:34 -0500

>
>
> res total affected
>
> ACE 1
> NHE 1
> )
> (no restraints)
>

this is fine, there is no error and it is just saying they aren't connected
at both sides.

>
> Should I fix the charge of the system so it is 0 before doing anything
> else? Like try adding Cl- ion to neutralize it?
>

this is up to you, and not really about this specific system. People often
add ions to neutralize the system, but it is not needed. You may also
consider if you want extra ions for salt, but this is up to you and may
depend on the conditions you are trying to model. Ions can be added in
tleap. More info can be found in this tutorial:
https://ambermd.org/tutorials/basic/tutorial7/index.php



> I also don't know if the ACE 1 and NHE1 not being recognized is going to
> cause a big problem in the simulation.
>
see above - it was recognized fine and there is no problem with this part.

On Mon, Dec 13, 2021 at 5:14 PM Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:

> I replied to this earlier today, search your email and check your spam
> folder
>
> On Mon, Dec 13, 2021 at 5:13 PM Yen Nguyen <ynguyen3.gmu.edu> wrote:
>
>> Hello,
>>
>> I tried building the structure from scratch in tleap using the following
>> command
>>
>> AbSOD-15aa-modified = sequence { ACE TRP GLU HID ALA TYR TYR ILE ASP PHE
>> ARG ASN LEU ARG PRO LYS NHE}
>>
>> and the output structure is correct when I viewed it in chimera.
>>
>> But when I tried to save the parm7 and rst7 files after solvating the
>> peptide, I got these messages.
>>
>> saveamberparm AbSOD-15aa-modified parm7 rst7
>> Checking Unit.
>>
>> Warning: The unperturbed charge of the unit (1.000000) is not zero.
>>
>> Note: Ignoring the warnings from Unit Checking.
>>
>> Building topology.
>> Building atom parameters.
>> Building bond parameters.
>> Building angle parameters.
>> Building proper torsion parameters.
>> Building improper torsion parameters.
>> total 73 improper torsions applied
>> Building H-Bond parameters.
>> Incorporating Non-Bonded adjustments.
>> Not Marking per-residue atom chain types.
>> Marking per-residue atom chain types.
>> (Residues lacking connect0/connect1 -
>> these don't have chain types marked:
>>
>> res total affected
>>
>> ACE 1
>> NHE 1
>> )
>> (no restraints)
>>
>> Should I fix the charge of the system so it is 0 before doing anything
>> else? Like try adding Cl- ion to neutralize it?
>> I also don't know if the ACE 1 and NHE1 not being recognized is going to
>> cause a big problem in the simulation.
>>
>> Thank you.
>>
>> Yen Nguyen
>> Graduate Research Assistant
>> George Mason University
>> Department of Chemistry and Biochemistry
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>>
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Dec 13 2021 - 14:30:04 PST
Custom Search