Using CPPTRAJ [59], one frame every 0.1 ns was retrieved from the MD trajectories and used as input for the MM-GBSA analysis. MM-GBSA analysis was carried out with the MMPBSA.py [60] module, determining the binding affinity between mTHPC and the different proteins. Electrostatic and van der Waals (vdW) interactions were calculated considering an infinite cut-off. The polar solvation term was calculated using the Generalized Born (GB) model, whereas the non-polar solvation term was determined using solvent-accessible, surface-area-dependent terms.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.