MCL1 perturbation 23~66
See original GitHub issuePer our conversation during the meeting, I think we should also use MCL1 as a test system for ‘validation’ because it almost 3x smaller than CDK2. The perturbation I am working with is simple and is located at indices 19 > 13 in MCL1_ligands.sdf.zip. MCL1_protein_fixed.pdb.zip is the protein that I prepared using pdbfixer.
However after preparing the topology proposal with NonequilibriumFEPSetup, I get an error when trying to build the HybridTopologyFactory:
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-17-9565cf892879> in <module>()
1 snefep = NonequilibriumSwitchingFEP(sfesetup[0], sfesetup[1], sfesetup[2], ncmc_nsteps=1000000,
2 nsteps_per_iteration=1000, trajectory_directory='mcl1_solvated_011118',
----> 3 trajectory_prefix='traj')
~/PycharmProjects/perses/perses/dispersed/relative_setup.py in __init__(self, topology_proposal, pos_old, new_positions, use_dispersion_correction, forward_functions, n_equil_steps, ncmc_nsteps, nsteps_per_iteration, temperature, trajectory_directory, trajectory_prefix, atom_selection, scheduler_address)
362
363 #construct the hybrid topology factory object
--> 364 self._factory = HybridTopologyFactory(topology_proposal, pos_old, new_positions, use_dispersion_correction=use_dispersion_correction)
365
366 #use default functions if none specified
~/PycharmProjects/perses/perses/annihilation/new_relative.py in __init__(self, topology_proposal, current_positions, new_positions, use_dispersion_correction, functions)
146
147 #verify that no constraints are changing over the course of the switching.
--> 148 self._constraint_check_fast()
149
150 #construct dictionary of exceptions in old and new systems
~/PycharmProjects/perses/perses/annihilation/new_relative.py in _constraint_check_fast(self)
455 for constraint_pair in overlapping_constraints:
456 if old_system_constraints[constraint_pair] != new_system_constraints[constraint_pair]:
--> 457 raise ValueError("There is a changing constraint length in this system.")
458
459 def _determine_interaction_group(self, atoms_in_interaction):
ValueError: There is a changing constraint length in this system.
I am currently trying to debug this. @pgrinaway can you give it a try on your system and see if you get the same error?
Issue Analytics
- State:
- Created 6 years ago
- Comments:6 (4 by maintainers)
Top Results From Across the Web
Akt is activated in chronic lymphocytic leukemia cells and ...
As MCL1 can be cleaved by caspases during apoptosis in CLL cells, we incubated cells with the Akt inhibitor in the presence or...
Read more >Targeting OXPHOS de novo purine synthesis as the nexus of ...
It has been reported to trigger apoptosis and arrest cell growth. On-target CDK9 inhibition by dinaciclib has been shown to play major roles...
Read more >MicroRNA-Regulated Signaling Pathways: Potential ... - MDPI
P16 is a well-described cyclin-dependent kinase inhibitor which binds to CDK4 and CDK6 [126]; it is the product of the CDKN2A gene and...
Read more >BAX/BAK-Induced Apoptosis Results in Caspase-8
Mechanistically, vioprolides inhibit MCL-1 and BCL2, which in turn triggers. BAX/BAK-dependent mitochondrial outer mem- brane permeabilization ( ...
Read more >Systemic mastocytosis in adults: 2013 update on diagnosis ...
Systemic mastocytosis in adults: 2013 update on diagnosis, ... ISM = 19, SSM = 3, ASM = 12, SM-AHNMD = 6, and MCL...
Read more >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
Ah, right–I had written that code before you changed the
TopologyProposal
to contain theold_alchemical_atoms
argument. You’re right, we should update the construction code.Sounds good