Opened 3 months ago

Closed 6 days ago

Last modified 5 days ago

#20109 closed enhancement (fixed)

Protonate parts of a model

Reported by: a.rohou@… Owned by: Eric Pettersen
Priority: moderate Milestone:
Component: Structure Editing Version:
Keywords: Cc: Elaine Meng
Blocked By: Blocking:
Notify when closed: Platform: all
Project: ChimeraX

Description

From chimerax-users thread "addh to only parts of a model?":

Alexis (partial): However I am facing a situation where I'd really, really like to add hydrogens only to waters (which were added by Servalcat). For reasons that I haven't got to the bottom of, addh adds hydrogens in the wrong places in a ligand I'm working on (a peptide with a mixture of natural and non-canonical amino acids, though I'm almost certain I've observed this behavior with other, less exotic, ligands).

Eric (partial): Probably the only feasible method for partially protonating a model would be to remember the initial protonation state of the model, add hydrogens, then remove any hydrogens that were added outside the target area.

Change History (4)

comment:1 by Eric Pettersen, 6 days ago

Resolution: fixed
Status: assignedclosed

addh now accepts an atom spec instead of a model spec and only protonates atoms in the atom spec

comment:2 by a.rohou@…, 6 days ago

Thanks Eric!

On Tue, Jun 16, 2026 at 11:41 AM ChimeraX <ChimeraX-bugs-admin@cgl.ucsf.edu>
wrote:

>
>
>
>

comment:3 by a.rohou@…, 6 days ago

Hi Eric,

If Claude summarized your fix correctly to me, I think your implementation
will have undesirable side effects. Suppose for example that I have an
already-protonated ligand, which I place into a protein that isn't yet
protonated. I then want to use "addh #1 & ~ligand" to protonate the protein
but leave the ligand alone. Claude reckons that addh would actually strip
hydrogens off of my ligand. I think this is undesirable behavior. The
command is called "addh" - it should add zero or more hydrogens, not remove
any...

#20109: Protonate parts of a model
>
>
>
>

comment:4 by Eric Pettersen, 5 days ago

That diabolical Claude is correct. I have committed a fix.

Note: See TracTickets for help on using tickets.