Browsing by Subject "Drug design"

Now showing 1 - 4 of 4
  • Thumbnail Image
    Item
    Editorial: Machine learning for peptide structure, function, and design
    (Frontiers Media, 2022-09-20) Ge, Ruiquan; Dong, Chuan; Wang, Juexin; Wei, Yanjie; Biomedical Engineering and Informatics, Luddy School of Informatics, Computing, and Engineering
  • Thumbnail Image
    Item
    Fast free energy estimates from λ-dynamics with bias-updated Gibbs sampling
    (Springer Nature, 2023-12-21) Robo, Michael T.; Hayes, Ryan L.; Ding, Xinqiang; Pulawski, Brian; Vilseck, Jonah Z.; Biochemistry and Molecular Biology, School of Medicine
    Relative binding free energy calculations have become an integral computational tool for lead optimization in structure-based drug design. Classical alchemical methods, including free energy perturbation or thermodynamic integration, compute relative free energy differences by transforming one molecule into another. However, these methods have high operational costs due to the need to perform many pairwise perturbations independently. To reduce costs and accelerate molecular design workflows, we present a method called λ-dynamics with bias-updated Gibbs sampling. This method uses dynamic biases to continuously sample between multiple ligand analogues collectively within a single simulation. We show that many relative binding free energies can be determined quickly with this approach without compromising accuracy. For five benchmark systems, agreement to experiment is high, with root mean square errors near or below 1.0 kcal mol-1. Free energy results are consistent with other computational approaches and within statistical noise of both methods (0.4 kcal mol-1 or less). Notably, large efficiency gains over thermodynamic integration of 18-66-fold for small perturbations and 100-200-fold for whole aromatic ring substitutions are observed. The rapid determination of relative binding free energies will enable larger chemical spaces to be more readily explored and structure-based drug design to be accelerated.
  • Thumbnail Image
    Item
    How to Sample Dozens of Substitutions per Site with λ Dynamics
    (American Chemical Society, 2024) Hayes, Ryan L.; Cervantes, Luis F.; Abad Santos, Justin Cruz; Samadi, Amirmasoud; Vilseck, Jonah Z.; Brooks, Charles L., III; Biochemistry and Molecular Biology, School of Medicine
    Alchemical free energy methods are useful in computer-aided drug design and computational protein design because they provide rigorous statistical mechanics-based estimates of free energy differences from molecular dynamics simulations. λ dynamics is a free energy method with the ability to characterize combinatorial chemical spaces spanning thousands of related systems within a single simulation, which gives it a distinct advantage over other alchemical free energy methods that are mostly limited to pairwise comparisons. Recently developed methods have improved the scalability of λ dynamics to perturbations at many sites; however, the size of chemical space that can be explored at each individual site has previously been limited to fewer than ten substituents. As the number of substituents increases, the volume of alchemical space corresponding to nonphysical alchemical intermediates grows exponentially relative to the size corresponding to the physical states of interest. Beyond nine substituents, λ dynamics simulations become lost in an alchemical morass of intermediate states. In this work, we introduce new biasing potentials that circumvent excessive sampling of intermediate states by favoring sampling of physical end points relative to alchemical intermediates. Additionally, we present a more scalable adaptive landscape flattening algorithm for these larger alchemical spaces. Finally, we show that this potential enables more efficient sampling in both protein and drug design test systems with up to 24 substituents per site, enabling, for the first time, simultaneous simulation of all 20 amino acids.
  • Thumbnail Image
    Item
    Optimizing Multisite λ-Dynamics Throughput with Charge Renormalization
    (American Chemical Society, 2022) Vilseck, Jonah Z.; Cervantes, Luis F.; Hayes, Ryan L.; Brooks, Charles L., III.; Biochemistry and Molecular Biology, School of Medicine
    With the ability to sample combinations of alchemical perturbations at multiple sites off a small molecule core, multisite λ-dynamics (MSλD) has become an attractive alternative to conventional alchemical free energy methods for exploring large combinatorial chemical spaces. However, current software implementations dictate that combinatorial sampling with MSλD must be performed with a multiple topology model (MTM), which is nontrivial to create by hand, especially for a series of ligand analogues which may have diverse functional groups attached. This work introduces an automated workflow, referred to as msld_py_prep, to assist in the creation of a MTM for use with MSλD. One approach for partitioning partial atomic charges between ligands to create a MTM, called charge renormalization, is also presented and rigorously evaluated. We find that msld_py_prep greatly accelerates the preparation of MSλD ready-to-use files and that charge renormalization can provide a successful approach for MTM generation, as long as bookending calculations are applied to correct small differences introduced by charge renormalization. Charge renormalization also facilitates the use of many different force field parameters with MSλD, broadening the applicability of MSλD for computer-aided drug design.