Project ID: plumID:22.033
Name: Reciprocal barrier restraint. Application to PROTAC passive permeability prediction
Keywords: PROTAC, membrane permeability, PMF, restraint, meta-eABF, metadynamics, DRR
PLUMED version: 2.7
Contributor: Istvan Kolossvary
Submitted on: 29 Aug 2022
Publication: I. Kolossváry, W. Sherman, A Different Kind of Restraint Suitable for Molecular Dynamics Simulations (2022)
PLUMED input files
Last tested: 01 Jun 2023, 22:52:21
Project description and instructions
To run this simulation you need OpenMM 7.6 or newer and the OpenMM-PLUMED plugin installed. To start the job run the following command line.
python run_protac-permeab-meta-eabf.py n_steps do_relax do_plumed. For example,
python run_protac-permeab-meta-eabf.py 50000 1 0 will minimize the system and run 50,000 steps of MD relaxation while
python run_protac-permeab-meta-eabf.py 5000000 0 1 will read the checkpoint file and continue 5 million steps of meta-eABF simulation. The collective variable is the normalized and periodic Z coordinate of the center of mass (COM) of the PROTAC2 molecule within the domain -0.5, 0.5 where the floating origin of the coordinate system corresponds to the COM of the membrane.NOTES. (1) The COM of the membrane should be computed via the PHASES keyword, the traditional calculation gives highly inaccurate results. (2) We found it beneficial to delay the meta-eABF bias significantly by setting FULLSAMPLES=200000 in the DRR action.
[v1] 29 Aug 2022: original submission
Click on the image below and get the code to add the badge to your website!