FEP Protocol Builder: Optimization of Free Energy Perturbation Protocols Using Active Learning

被引:7
作者
de Oliveira, Cesar [1 ]
Leswing, Karl [2 ]
Feng, Shulu [2 ]
Kanters, Reneï [2 ]
Abel, Robert [2 ]
Bhat, Sathesh [2 ]
机构
[1] Schrodinger Inc, San Diego, CA 92121 USA
[2] Schrodinger Inc, New York, NY 10036 USA
关键词
DRUG DISCOVERY; LEAD OPTIMIZATION; REPLICA EXCHANGE; PROTEIN; INHIBITORS; DESIGN; WATER; SIMULATIONS; PREDICTIONS; ACCURACY;
D O I
10.1021/acs.jcim.3c00681
中图分类号
R914 [药物化学];
学科分类号
100701 ;
摘要
Significant improvements have been made in the past decadeto methodsthat rapidly and accurately predict binding affinity through freeenergy perturbation (FEP) calculations. This has been driven by recentadvances in small-molecule force fields and sampling algorithms combinedwith the availability of low-cost parallel computing. Predictive accuraciesof & SIM;1 kcal mol(-1) have been regularly achieved,which are sufficient to drive potency optimization in modern drugdiscovery campaigns. Despite the robustness of these FEP approachesacross multiple target classes, there are invariably target systemsthat do not display expected performance with default FEP settings.Traditionally, these systems required labor-intensive manual protocoldevelopment to arrive at parameter settings that produce a predictiveFEP model. Due to the (a) relatively large parameter space to be explored,(b) significant compute requirements, and (c) limited understandingof how combinations of parameters can affect FEP performance, manualFEP protocol optimization can take weeks to months to complete, andoften does not involve rigorous train-test set splits, resulting inpotential overfitting. These manual FEP protocol development timelinesdo not coincide with tight drug discovery project timelines, essentiallypreventing the use of FEP calculations for these target systems. Here,we describe an automated workflow termed FEP Protocol Builder (FEP-PB)to rapidly generate accurate FEP protocols for systems that do notperform well with default settings. FEP-PB uses an active-learningworkflow to iteratively search the protocol parameter space to developaccurate FEP protocols. To validate this approach, we applied it topharmaceutically relevant systems where default FEP settings couldnot produce predictive models. We demonstrate that FEP-PB can rapidlygenerate accurate FEP protocols for the previously challenging MCL1system with limited human intervention. We also apply FEP-PB in areal-world drug discovery setting to generate an accurate FEP protocolfor the p97 system. FEP-PB is able to generate a more accurate protocolthan the expert user, rapidly validating p97 as amenable to free energycalculations. Additionally, through the active-learning workflow,we are able to gain insight into which parameters are most importantfor a given system. These results suggest that FEP-PB is a robusttool that can aid in rapidly developing accurate FEP protocols andincreasing the number of targets that are amenable to the technology.
引用
收藏
页码:5592 / 5603
页数:12
相关论文
共 69 条
[1]   Advancing Drug Discovery through Enhanced Free Energy Calculations [J].
Abel, Robert ;
Wang, Lingle ;
Harder, Edward D. ;
Berne, B. J. ;
Friesner, Richard A. .
ACCOUNTS OF CHEMICAL RESEARCH, 2017, 50 (07) :1625-1632
[2]   Is Structure-Based Drug Design Ready for Selectivity Optimization? [J].
Albanese, Steven K. ;
Chodera, John D. ;
Volkamer, Andrea ;
Keng, Simon ;
Abel, Robert ;
Wang, Lingle .
JOURNAL OF CHEMICAL INFORMATION AND MODELING, 2020, 60 (12) :6211-6227
[3]   Accurate calculation of the absolute free energy of binding for drug molecules [J].
Aldeghi, Matteo ;
Heifetz, Alexander ;
Bodkin, Michael J. ;
Knappcd, Stefan ;
Biggin, Philip C. .
CHEMICAL SCIENCE, 2016, 7 (01) :207-218
[4]   EFFICIENT ESTIMATION OF FREE-ENERGY DIFFERENCES FROM MONTE-CARLO DATA [J].
BENNETT, CH .
JOURNAL OF COMPUTATIONAL PHYSICS, 1976, 22 (02) :245-268
[5]   MOLECULAR-DYNAMICS WITH COUPLING TO AN EXTERNAL BATH [J].
BERENDSEN, HJC ;
POSTMA, JPM ;
VANGUNSTEREN, WF ;
DINOLA, A ;
HAAK, JR .
JOURNAL OF CHEMICAL PHYSICS, 1984, 81 (08) :3684-3690
[6]  
Berendsen HJC., 1981, Interaction Model for Water in Relation to Protein Hydration, DOI [DOI 10.1007/978-94-015-7658-121, 10.1007/978-94-015-7658-121, 10.1007/978-94-015-7658-1_21]
[7]   MCL-1 inhibitors, fast-lane development of a new class of anti-cancer agents [J].
Bolomsky, Arnold ;
Vogler, Meike ;
Kose, Murat Cem ;
Heckman, Caroline A. ;
Ehx, Gregory ;
Ludwig, Heinz ;
Caers, Jo .
JOURNAL OF HEMATOLOGY & ONCOLOGY, 2020, 13 (01)
[8]   Water Network Perturbation in Ligand Binding: Adenosine A2A Antagonists as a Case Study [J].
Bortolato, Andrea ;
Tehan, Ben G. ;
Bodnarchuk, Michael S. ;
Essex, Jonathan W. ;
Mason, Jonathan S. .
JOURNAL OF CHEMICAL INFORMATION AND MODELING, 2013, 53 (07) :1700-1713
[9]   AutoDesigner, a De Novo Design Algorithm for Rapidly Exploring Large Chemical Space for Lead Optimization: Application to the Design and Synthesis of D-Amino Acid Oxidase Inhibitors [J].
Bos, Pieter H. ;
Houang, Evelyne M. ;
Ranalli, Fabio ;
Leffler, Abba E. ;
Boyles, Nicholas A. ;
Eyrich, Volker A. ;
Luria, Yuval ;
Katz, Dana ;
Tang, Haifeng ;
Abel, Robert ;
Bhat, Sathesh .
JOURNAL OF CHEMICAL INFORMATION AND MODELING, 2022, 62 (08) :1905-1915
[10]   Walking Outcome After Traumatic Paraplegic Spinal Cord Injury: The Function of Which Myotomes Makes a Difference? [J].
Cathomen, Adrian ;
Maier, Doris ;
Kriz, Jiri ;
Abel, Rainer ;
Roehrich, Frank ;
Baumberger, Michael ;
Scivoletto, Giorgio ;
Weidner, Norbert ;
Rupp, Ruediger ;
Jutzeler, Catherine R. ;
Steeves, John D. ;
Curt, Armin ;
Bolliger, Marc .
NEUROREHABILITATION AND NEURAL REPAIR, 2023, 37 (05) :316-327