Ligand binding free energy evaluation by Monte Carlo Recursion

Autor(es)

de Souza, Joao Victor ; Nogueira, Victor H. R. ; Nascimento, Alessandro S.

• Tipo de documento: Artigo Científico
• Fonte: COMPUTATIONAL BIOLOGY AND CHEMISTRY; v. 103, p. 12-pg., 2023-02-20.

Resumo

The correct evaluation of ligand binding free energies by computational methods is still a very challenging active area of research. The most employed methods for these calculations can be roughly classified into four groups: (i) the fastest and less accurate methods, such as molecular docking, designed to sample a large number of mole-cules and rapidly rank them according to the potential binding energy; (ii) the second class of methods use a thermodynamic ensemble, typically generated by molecular dynamics, to analyze the endpoints of the ther-modynamic cycle for binding and extract differences, in the so-called 'end-point' methods; (iii) the third class of methods is based on the Zwanzig relationship and computes the free energy difference after a chemical change of the system (alchemical methods); and (iv) methods based on biased simulations, such as metadynamics, for example. These methods require increased computational power and as expected, result in increased accuracy for the determination of the strength of binding. Here, we describe an intermediate approach, based on the Monte Carlo Recursion (MCR) method first developed by Harold Scheraga. In this method, the system is sampled at increasing effective temperatures, and the free energy of the system is assessed from a series of terms W(b, T), computed from Monte Carlo (MC) averages at each iteration. We show the application of the MCR for ligand binding with datasets of guest-hosts systems (N = 75) and we observed that a good correlation is obtained be-tween experimental data and the binding energies computed with MCR. We also compared the experimental data with an end-point calculation from equilibrium Monte Carlo calculations that allowed us to conclude that the lower-energy (lower-temperature) terms in the calculation are the most relevant to the estimation of the binding energies, resulting in similar correlations between MCR and MC data and the experimental values. On the other hand, the MCR method provides a reasonable view of the binding energy funnel, with possible connections with the ligand binding kinetics, as well. The codes developed for this analysis are publicly available on GitHub as a part of the LiBELa/MCLiBELa project (https://github.com/alessandronascimento/LiBELa). (AU)

• Processo FAPESP: 17/18173-0 - Mecanismos Envolvidos em Resistência a Antibióticos: Parede de Ácidos Teicóicos e Biofilmes como Alvos Moleculares
• Beneficiário: Alessandro Silva Nascimento
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 20/03983-9 - Investigação estrutural da síntese de ramnose e implicações na resistência a antibióticos
• Beneficiário: Alessandro Silva Nascimento
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 14/06565-2 - Estendendo as fronteiras em interações biomoleculares: docking e avaliação da energia livre
• Beneficiário: Alessandro Silva Nascimento
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 10/15376-8 - Estudo prospectivo de ligantes do receptor PPAR gama
• Beneficiário: Alessandro Silva Nascimento
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 14/01751-2 - Aplicação de Monte Carlo em um algoritmo de docking para a geração de ensembles e cálculos de energia livre
• Beneficiário: João Victor de Souza Cunha
• Modalidade de apoio: Bolsas no Brasil - Mestrado
• Processo FAPESP: 15/01709-9 - Cálculos de energia livre para medidas computacionais da entalpia e entropia em interações proteína-ligante
• Beneficiário: João Victor de Souza Cunha
• Modalidade de apoio: Bolsas no Exterior - Estágio de Pesquisa - Mestrado
• Processo FAPESP: 15/26722-8 - Drug discovery contra doenças infecciosas humanos
• Beneficiário: Carsten Wrenger
• Modalidade de apoio: Auxílio à Pesquisa - Temático