A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics

Drug-protein binding kinetic parameters are key indicators of drug efficacy, but their experimental determination is often costly and time-consuming. Computational approaches require optimization of the in silico approximations to achieve sufficient accuracy while remaining computationally feasible....

Full description

Saved in:
Bibliographic Details
Main Authors: Muñiz Chicharro, Abraham (Author) , Ganotra, Gaurav K. (Author) , Wade, Rebecca C. (Author)
Format: Article (Journal)
Language:English
Published: October 4, 2025
In: Journal of chemical information and modeling
Year: 2025, Volume: 65, Issue: 20, Pages: 11215-11231
ISSN:1549-960X
DOI:10.1021/acs.jcim.5c01488
Online Access:Verlag, kostenfrei, Volltext: https://doi.org/10.1021/acs.jcim.5c01488
Get full text
Author Notes:Abraham Muñiz-Chicharro, Gaurav K. Ganotra, and Rebecca C. Wade

MARC

LEADER 00000naa a2200000 c 4500
001 1950433889
003 DE-627
005 20260128104105.0
007 cr uuu---uuuuu
008 260128s2025 xx |||||o 00| ||eng c
024 7 |a 10.1021/acs.jcim.5c01488  |2 doi 
035 |a (DE-627)1950433889 
035 |a (DE-599)KXP1950433889 
040 |a DE-627  |b ger  |c DE-627  |e rda 
041 |a eng 
084 |a 32  |2 sdnb 
100 1 |a Muñiz Chicharro, Abraham  |d 1990-  |e VerfasserIn  |0 (DE-588)1254724451  |0 (DE-627)179744140X  |4 aut 
245 1 2 |a A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics  |c Abraham Muñiz-Chicharro, Gaurav K. Ganotra, and Rebecca C. Wade 
264 1 |c October 4, 2025 
300 |b Diagramme, Illustrationen 
300 |a 17 
336 |a Text  |b txt  |2 rdacontent 
337 |a Computermedien  |b c  |2 rdamedia 
338 |a Online-Ressource  |b cr  |2 rdacarrier 
500 |a Gesehen am 28.01.2026 
520 |a Drug-protein binding kinetic parameters are key indicators of drug efficacy, but their experimental determination is often costly and time-consuming. Computational approaches require optimization of the in silico approximations to achieve sufficient accuracy while remaining computationally feasible. Here, the combination of Brownian dynamics (BD) and molecular dynamics (MD) has been investigated for this purpose. BD is used for simulating long-range diffusion and diffusional encounter complex formation, whereas MD captures the subsequent formation of the bound complex, providing a detailed treatment of short-range interactions and molecular flexibility. While existing methods that employ this approach have successfully yielded estimated association rate constants (kon), they often require extensive computational resources. In this work, we developed a multiscale pipeline that achieves improved computational efficiency by optimizing the sampling by BD simulation to generate an ensemble of diffusional encounter complexes in which the ligand comes very close to its protein binding site and then uses these as starting structures for MD simulation. Due to the much lower computational cost of BD simulation and the reduced MD simulation time, the approach is computationally efficient while preserving accuracy. The pipeline has been validated for a diverse set of protein-ligand complexes, varying in size, flexibility, and binding properties, yielding kon values that align well with experimental measurements, as well as insights into the physical determinants of the association rate. 
700 1 |a Ganotra, Gaurav K.  |d 1989-  |e VerfasserIn  |0 (DE-588)1195542747  |0 (DE-627)1677559217  |4 aut 
700 1 |a Wade, Rebecca C.  |e VerfasserIn  |0 (DE-588)102801774X  |0 (DE-627)730136000  |0 (DE-576)276591402  |4 aut 
773 0 8 |i Enthalten in  |t Journal of chemical information and modeling  |d Washington, DC : American Chemical Society, 2005  |g 65(2025), 20, Seite 11215-11231  |h Online-Ressource  |w (DE-627)302467327  |w (DE-600)1491237-5  |w (DE-576)090855124  |x 1549-960X  |7 nnas  |a A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics 
773 1 8 |g volume:65  |g year:2025  |g number:20  |g pages:11215-11231  |g extent:17  |a A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics 
856 4 0 |u https://doi.org/10.1021/acs.jcim.5c01488  |x Verlag  |x Resolving-System  |z kostenfrei  |3 Volltext  |7 0 
951 |a AR 
992 |a 20260128 
993 |a Article 
994 |a 2025 
998 |g 102801774X  |a Wade, Rebecca C.  |m 102801774X:Wade, Rebecca C.  |d 700000  |d 706000  |e 700000PW102801774X  |e 706000PW102801774X  |k 0/700000/  |k 1/700000/706000/  |p 3  |y j 
999 |a KXP-PPN1950433889  |e 4863032846 
BIB |a Y 
SER |a journal 
JSO |a {"note":["Gesehen am 28.01.2026"],"relHost":[{"title":[{"title_sort":"Journal of chemical information and modeling","title":"Journal of chemical information and modeling"}],"language":["eng"],"pubHistory":["Volume 45, issue 1 (January 2005)-"],"physDesc":[{"extent":"Online-Ressource"}],"note":["Gesehen am 03.05.2023"],"type":{"media":"Online-Ressource","bibl":"periodical"},"part":{"volume":"65","pages":"11215-11231","extent":"17","text":"65(2025), 20, Seite 11215-11231","year":"2025","issue":"20"},"origin":[{"publisherPlace":"Washington, DC","dateIssuedDisp":"[2005]-","publisher":"American Chemical Society"}],"id":{"zdb":["1491237-5"],"issn":["1549-960X"],"eki":["302467327"]},"disp":"A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamicsJournal of chemical information and modeling","recId":"302467327"}],"person":[{"display":"Muñiz Chicharro, Abraham","family":"Muñiz Chicharro","given":"Abraham","role":"aut"},{"display":"Ganotra, Gaurav K.","family":"Ganotra","given":"Gaurav K.","role":"aut"},{"family":"Wade","role":"aut","given":"Rebecca C.","display":"Wade, Rebecca C."}],"title":[{"title":"A multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics","title_sort":"multiscale simulation approach to compute protein-ligand association rate constants by combining Brownian dynamics and molecular dynamics"}],"physDesc":[{"extent":"17 S.","noteIll":"Diagramme, Illustrationen"}],"language":["eng"],"type":{"bibl":"article-journal","media":"Online-Ressource"},"id":{"eki":["1950433889"],"doi":["10.1021/acs.jcim.5c01488"]},"origin":[{"dateIssuedKey":"2025","dateIssuedDisp":"October 4, 2025"}],"name":{"displayForm":["Abraham Muñiz-Chicharro, Gaurav K. Ganotra, and Rebecca C. Wade"]},"recId":"1950433889"} 
SRT |a MUNIZCHICHMULTISCALE4202