Ligand migration and steered molecular dynamics in drug discovery: Comment on "Ligand diffusion in proteins via enhanced sampling in molecular dynamics" by Jakub Rydzewski and Wieslaw Nowak

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction—when ligand unbinds along the RDM vector—results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function.

Probing the Binding Affinity by Jarzynski's Nonequilibrium Binding Free Energy and Rupture Time

Binding affinity of a small ligand to a receptor is the important quantity in drug design, and it might be characterized by different quantities. The most popular one is the binding free energy, which can be estimated by several methods in conventional molecular dynamics simulation. So far in steered molecular dynamics (SMD), one can use either the rupture force or nonequilibrium pulling work as a measure for binding affinity. In this paper, we have shown that the nonequilibrium binding free energy Δ G_neq^Jar, obtained by Jarzynski's equality at a finite pulling speed, has good correlation with experimental data on inhibition constants, implying that this quantity can be used as a good scoring function for binding affinity. A similar correlation has also been disclosed for binding and unbinding free energy barriers. Applying the SMD method to unbinding of 23 small compounds from the binding site of β-lactamase protein, a bacteria-produced enzyme, we have demonstrated that the rupture or unbinding time strongly correlates with experimental data with correlation level R ≈ 0.84. As follows from the Jarzynski's equality, the rupture time depends on the unbinding barrier exponentially. We show that Δ G_neq^Jar, the rupture time, and binding and unbinding free energy barriers are good descriptors for binding affinity. Our observation may be useful for fast screening of potential leads as the SMD simulation is not time-consuming. On the basis of nonequilibrium simulation, we disclosed that, in agreement with the experiment, the binding time is much longer than the unbinding one.