Skip to content

[docs] we should point out some example(s) on how to read out the trajectories from our output files  #916

Open

Description

@mikemhenry

@RiesBen
I checked again. State A is Valine and State B is Alanine. So, the ddg should be positive.

Also, I am wondering why the intermediates are needed in OpenFE if they are not used in perses and FEP+. I think FEP+ uses restraints on the backbone of the protein/peptide which may help with convergence. Is there a way to do the same with openFE to test if this is the source of convergence problems?

I think it would also be useful to visualize the trajectories. Are there any scripts available to extract trajectories from the saved simulation.nc file?

Originally posted by @amin-sagar in #889 (reply in thread)

Activity

  1. RiesBen commented on Aug 6, 2024

    @RiesBen
    Contributor

    @mikemhenry

    regarding the traj extraction have a look here: :)

    import MDAnalysis as mda
    from openfe_analysis import FEReader
    
    with mda.Writer('out.xtc', n_atoms=len(u.atoms)) as w:
        for ts in u.trajectory:
            w.write(u.atoms)

    Source: https://github.com/OpenFreeEnergy/openfe_analysis

  2. self-assigned this
    on Nov 9, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions