Search

Atomistic modeling and simulation methods facilitate biomedical research from many respects, including structure-based drug design. The ability of these methods to address biologically relevant problems is largely determined by the accuracy of the treatment of complex solvation effects in target biomolecules surrounded by water. The implicit solvent model – which treats solvent as a continuum with the dielectric and non-polar properties of water – offers a good balance between accuracy and speed. Simple and efficient, generalized Born (GB) model has become a widely used implicit solvent responsible for the estimation of key electrostatic interactions. The main goal of this research is to improve the accuracy of protein-ligand binding calculations in the implicit solvent framework. To address the problem (1) GBNSR6, an accurate yet efficient flavor of GB, has been thoroughly explored in the context of protein-ligand binding, (2) a global multidimensional optimization pipeline is developed to find the optimal dielectric boundary made of atomic and water probe radii specifically for protein-ligand binding calculations using GBNSR6. The pipeline includes (3) two novel post-processing steps for optimum robustness analysis and optimization landscape visualization. In the final step of this research, (4) accuracy gain the optimal dielectric boundary can bring in practice is explored on binding benchmarks, including the SARS-CoV-2 spike receptor-binding domain and the human ACE2 receptor. Doctor of Philosophy Drug discovery is one of the most challenging tasks in biological sciences as it takes about 10-15 years and $1.5-2 billion on average to discover a new drug. Therefore, efforts to speed up this process or lower its costs are highly valuable. Computer-aided drug design (CADD) plays a crucial role in the early stage of drug discovery. In CADD, computational approaches are used in order to discover, develop, and analyze drugs and similar biologically active molecules, such as proteins. Proteins are an important class of biological macromolecules that perform their functionality mainly through interactions with other molecules, for example, binding to small molecules so-called ligands. Thorough understanding of protein-ligand interactions is central to comprehending biology at the molecular level. In this study, we introduce and analyze a computational model used for protein-ligand binding free energy calculations. A global multidimensional optimization pipeline is developed to find the optimal parameters of the model,˘aparticularly˘athose parameters involved in the dielectric boundary. In order to examine the robustness of the optimal model to unavoidable perturbations and uncertainties, virtually inevitable in any complex system being optimized, a novel robustness metric is introduced. Finally, the robust optimal model is tested on protein-ligand benchmarks, including a complex related to the novel coronavirus. Results demonstrate relatively higher accuracy in terms of binding free energy calculations compared to reference models.

The ability to estimate protein-protein binding free energy in a computationally efficient via a physics-based approach is beneficial to research focused on the mechanism of viruses binding to their target proteins. Implicit solvation methodology may be particularly useful in the early stages of such research, as it can offer valuable insights into the binding process, quickly. Here we evaluate the potential of the related molecular mechanics generalized Born surface area (MMGB/SA) approach to estimate the binding free energy Δ G bind between the SARS-CoV-2 spike receptor-binding domain and the human ACE2 receptor. The calculations are based on a recent flavor of the generalized Born model, GBNSR6. Two estimates of Δ G bind are performed: one based on standard bondi radii, and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein-ligand binding. We take the average of the resulting two Δ G bind values as the consensus estimate. For the well-studied Ras-Raf protein-protein complex, which has similar binding free energy to that of the SARS-CoV-2/ACE2 complex, the consensus Δ G bind = -11.8 ± 1 kcal/mol, vs. experimental -9.7 ± 0.2 kcal/mol. The consensus estimates for the SARS-CoV-2/ACE2 complex is Δ G bind = -9.4 ± 1.5 kcal/mol, which is in near quantitative agreement with experiment (-10.6 kcal/mol). The availability of a conceptually simple MMGB/SA-based protocol for analysis of the SARS-CoV-2 /ACE2 binding may be beneficial in light of the need to move forward fast.