Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Binding Selectivity Analysis from Alchemical Receptor Hopping and Swapping Free Energy Calculations

We present receptor hopping and receptor swapping free energy estimation protocols based on the Alchemical Transfer Method (ATM) to model the binding selectivity of a set of ligands to two arbitrary receptors. The receptor hopping protocol, where a ligand is alchemically transferred from one receptor to another in one simulation, directly yields the ligand's binding selectivity free energy (BSFE) for the two receptors, which is the difference between the two individual binding free energies. In the receptor swapping protocol, the first ligand of a pair is transferred from one receptor to another while the second ligand is simultaneously transferred in the opposite direction. The receptor swapping free energy yields the differences in binding selectivity free energies of a set of ligands, which, when combined using a generalized DiffNet algorithm, yield the binding selectivity free energies of the ligands. We test these algorithms on host-guest systems and show that they yield results that agree with experimental data and are consistent with differences in absolute and relative binding free energies obtained by conventional methods. Preliminary applications to the selectivity analysis of molecular fragments binding to the trypsin and thrombin serine protease confirm the potential of the receptor swapping technology in structure-based drug discovery. The novel methodologies presented in this work are a first step toward streamlined and computationally efficient protocols for ligand selectivity optimization between mutants and homologous proteins.

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins

Modulating the function of Intrinsically Disordered Proteins (IDPs) with small molecules is of considerable importance given the crucial roles of IDPs in the pathophysiology of numerous diseases. Reported binding affinities for ligands to diverse IDPs vary broadly, and little is known about the detailed molecular mechanisms that underpin ligand efficacy. Molecular simulations of IDP ligand binding mechanisms can help us understand the mode of action of small molecule inhibitors of IDP function, but it is still unclear how binding energies can be modeled rigorously for such a flexible class of proteins. Here, we compare alchemical absolute binding free energy calculations (ABFE) and Markov-State Modeling (MSM) protocols to model the binding of the small molecule 10058-F4 to a disordered peptide extracted from a segment of the oncoprotein c-Myc. The ABFE results produce binding energy estimates that are sensitive to the choice of reference structure. In contrast, the MSM results produce more reproducible binding energy estimates consistent with weak mM binding affinities and transient intermolecular contacts reported in the literature.

Analysis of Glycan Recognition by Concanavalin A Using Absolute Binding Free Energy Calculations

Carbohydrates are key biological mediators of molecular recognition and signaling processes. In this case study, we explore the ability of absolute binding free energy (ABFE) calculations to predict the affinities of a set of five related carbohydrate ligands for the lectin protein, concanavalin A, ranging from 27-atom monosaccharides to a 120-atom complex-type N-linked glycan core pentasaccharide. ABFE calculations quantitatively rank and estimate the affinity of the ligands in relation to microcalorimetry, with a mean signed error in the binding free energy of -0.63 ± 0.04 kcal/mol. Consequently, the diminished binding efficiencies of the larger carbohydrate ligands are closely reproduced: the ligand efficiency values from isothermal titration calorimetry for the glycan core pentasaccharide and its constituent trisaccharide and monosaccharide compounds are respectively -0.14, -0.22, and -0.41 kcal/mol per heavy atom. ABFE calculations predict these ligand efficiencies to be -0.14 ± 0.02, -0.24 ± 0.03, and -0.46 ± 0.06 kcal/mol per heavy atom, respectively. Consequently, the ABFE method correctly identifies the high affinity of the key anchoring mannose residue and the negligible contribution to binding of both β-GlcNAc arms of the pentasaccharide. While challenges remain in sampling the conformation and interactions of these polar, flexible, and weakly bound ligands, we nevertheless find that the ABFE method performs well for this lectin system. The approach shows promise as a quantitative tool for predicting and deconvoluting carbohydrate-protein interactions, with potential application to design of therapeutics, vaccines, and diagnostics.

Prevalence of the genetic variant rs61330082 and serum levels of the visfatin gene in Mexican individuals with metabolic syndrome: a clinical and bioinformatics approach

BACKGROUND

metabolic syndrome (MetS) is a group of clinical anomalies that share an inflammatory component of multifactorial etiology.

OBJECTIVES

the present study aims to relate the genetic variant (rs61330082 C/T) with dietary patterns in the presence of MetS and the application of molecular docking according to the genotype and associated transcription factors.

METHODS

197 individuals aged 18 to 65 were included, from whom anthropometric measurements were taken, and a blood sample from the forearm. DNA extraction and enzymatic digestion were performed to determine the genotype of each participant by PCR-RFLP. Dietary patterns were analyzed using a nutritional questionnaire validated for the Mexican population. Serum levels of the protein visfatin were assessed by ELISA. Finally, bioinformatics tools were used for molecular docking to infer the binding of transcriptional factors in the polymorphic region.

RESULTS

the TT genotype was present in only 10 % of the population. Women carrying the CT+TT genotype, according to the dominant genetic model, had higher serum levels of triglycerides and VDLD-C. Statistical analysis did not show a significant association between the presence of MetS and the dominant CT+TT model (OR = 1.41, 95 % CI = 0.61-3.44, p = 0.53). We identified PAX5 as a transcription factor binding to the polymorphic site of this genetic variant.

CONCLUSIONS

this study demonstrated a significant association between the genetic variant (rs61330082 C/T) and lipid parameters. Women carrying the T allele have a higher risk of high triglyceride levels, a criterion for metabolic syndrome.

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

Advances in molecular glues: exploring chemical space and design principles for targeted protein degradation

The discovery of the E3 ligase cereblon (CRBN) as the target of thalidomide and its analogs revolutionized the field of targeted protein degradation (TPD). This ubiquitin-mediated degradation pathway was first harnessed by bivalent degraders. Recently, the emergence of low-molecular-weight molecular glue degraders (MGDs) has expanded the TPD landscape, because MGDs operate via the same mechanism while offering attractive physicochemical properties that are consistent with small-molecule therapeutics. This review delves into the discovery and advancement of MGDs, with case studies on cyclin K and the zinc finger protein IKZF2, highlighting the design principles, biological assays and therapeutic applications. Additionally, it examines the chemical space of molecular glues and outlines the collaborative efforts that are fueling innovation in this field.

A Method for Treating Significant Conformational Changes in Alchemical Free Energy Simulations of Protein-Ligand Binding

Relative binding free energy (RBFE) simulation is a rigorous approach to the calculation of quantitatively accurate binding free energy values for protein-ligand binding in which a reference binder is gradually converted to a target binder through alchemical transformation during the simulation. The success of such simulations relies on being able to accurately sample the correct conformational phase space for each alchemical state; however, this becomes a challenge when a significant conformation change occurs between the reference and target binder-receptor complexes. Increasing the simulation time and using enhanced sampling methods can be helpful, but effects can be limited, especially when the free energy barrier between conformations is high or when the correct target complex conformation is difficult to find and maintain. Current RBFE protocols seed the reference complex structure into every alchemical window of the simulation. In our study, we describe an improved protocol in which the reference structure is seeded into the first half of the alchemical windows, and the target structure is seeded into the second half of the alchemical windows. By applying information about the relevant correct end point conformations to different simulation windows from the beginning, the need for large barrier crossings or simulation prediction of the correct structures during an alchemical simulation is in many cases obviated. In the diverse cases we examine below, the simulations yielded free energy predictions that are satisfactory as compared to experiment and superior to running the simulations utilizing the conventional protocol. The method is straightforward to implement for publicly available FEP workflows.

Hit me with your best shot: Integrated hit discovery for the next generation of drug targets

Identification of high-quality hit chemical matter is of vital importance to the success of drug discovery campaigns. However, this goal is becoming ever harder to achieve as the targets entering the portfolios of pharmaceutical and biotechnology companies are increasingly trending towards novel and traditionally challenging to drug. This demand has fuelled the development and adoption of numerous new screening approaches, whereby the contemporary hit identification toolbox comprises a growing number of orthogonal and complementary technologies including high-throughput screening, fragment-based ligand design, affinity screening (affinity-selection mass spectrometry, differential scanning fluorimetry, DNA-encoded library screening), as well as increasingly sophisticated computational predictive approaches. Herein we describe how an integrated strategy for hit discovery, whereby multiple hit identification techniques are tactically applied, selected in the context of target suitability and resource priority, represents an optimal and often essential approach to maximise the likelihood of identifying quality starting points from which to develop the next generation of medicines.

Calculated hydration free energies become less accurate with increases in molecular weight

In order for computer-aided drug design to fulfil its long held promise of delivering new medicines faster and cheaper, extensive development and validation work must be done first. This pertains particularly to molecular dynamics force fields where one important aspect-the hydration free energy (HFE) of small molecules-is often insufficiently analyzed. While most benchmarking studies report excellent accuracies of calculated hydration free energies-usually within 2 kcal/mol of experimental values-we find that deeper analysis reveals significant shortcomings. Herein, we report a dependence of HFE prediction errors on ligand molecular weight-the higher the weight, the bigger the prediction error and the higher the probability the calculated result is erroneous by a large amount. We show that in the drug-like molecular weight region, HFE predictions can easily be off by 5 kcal/mol or more. This is likely to be highly problematic in a drug discovery and development setting. We make our HFE results and molecular descriptors freely and fully available in order to encourage deeper analysis of future molecular dynamics results and facilitate development of the next generation of force fields.

Assessment of the in vitro anti-diabetic activity with molecular dynamic simulations of limonoids isolated from Adalia lemon peels

Limonoids are important constituents of citrus that have a significant impact on promoting human health. Therefore, the primary focus of this research was to assess the overall limonoid content and isolate limonoids from Adalia lemon (Citrus limon L.) peels for their potential use as antioxidants and anti-diabetic agents. The levels of limonoid aglycones in the C. limon peel extract were quantified through a colorimetric assay, revealing a concentration of 16.53 ± 0.93 mg/L limonin equivalent. Furthermore, the total concentration of limonoid glucosides was determined to be 54.38 ± 1.02 mg/L. The study successfully identified five isolated limonoids, namely limonin, deacetylnomilin, nomilin, obacunone 17-O-β-D-glucopyranoside, and limonin 17-O-β-D-glucopyranoside, along with their respective yields. The efficacy of the limonoids-rich extract and the five isolated compounds was evaluated at three different concentrations (50, 100, and 200 µg/mL). It was found that both obacunone 17-O-β-D-glucopyranoside and limonin 17-O-β-D-glucopyranoside possessed the highest antioxidant, free radical scavenging, and anti-diabetic activities, followed by deacetylnomilin, and then the limonoids-rich extract. The molecular dynamic simulations were conducted to predict the behavior of the isolated compounds upon binding to the protein's active site, as well as their interaction and stability. The results revealed that limonin 17-O-β-D-glucopyranoside bound to the protein complex system exhibited a relatively more stable conformation than the Apo system. The analysis of Solvent Accessible Surface Area (SASA), in conjunction with the data obtained from Root-Mean-Square Deviation (RMSD), Root-Mean-Square Fluctuation (RMSF), and Radius of Gyration (ROG) computations, provided further evidence that the limonin 17-O-β-D-glucopyranoside complex system remained stable within the catalytic domain binding site of the human pancreatic alpha-amylase (HPA)-receptor. The research findings suggest that the limonoids found in Adalia lemon peels have the potential to be used as effective natural substances in creating innovative therapeutic treatments for conditions related to oxidative stress and disorders in carbohydrate metabolism.

Effect of Protein-Polarized Ligand Charges on Relative Protein Ligand Binding Affinities

A major challenge in computer-aided drug design is predicting relative binding energies of different molecules to a target protein using fast and accurate free-energy calculation methods. Free-energy calculations are primarily computed by utilizing classical molecular dynamics simulations based on all-atom force fields (FF) to model the interactions in the system. The present standard classical all-atom FFs contain fixed partial charges on the atoms, and hence electrostatic interactions are modeled between them. The parametrization process to determine these partial charges usually relies on quantum mechanics or semiempirical calculations of the molecule in the gas phase or homogeneous water surrounding. These present standard parametrization schemes of the partial charges neglect, therefore, polarization effects from the protein surrounding. The absence of protein polarization effects can lead to significant errors in free-energy calculations in proteins. We present a parametrization scheme for the partial charges of ligands, named protein-induced polarization (PIP) charges, which account for the electrostatic polarization due to the protein surrounding. The scheme involves single-point quantum mechanics/molecular mechanics calculations of the ligand charges in the protein/water surrounding. Using PIP ligand partial charges, we have calculated the relative binding free energies (RBFEs) of well-studied protein-ligand systems. We show here that RBFEs computed with PIP charges are either significantly improved or at least comparable to those computed with nonpolarized standard GAFF charges. Overall, we present a simple-to-use parametrization scheme to include protein polarization in any type of binding free-energy calculations. The parametrization scheme increases the accuracy in RBFE calculations, while it does not add significant computation time to standard parametrization procedures.

Relative genotoxicity of polycyclic aromatic hydrocarbons inferred from free energy perturbation approaches

Utilizing molecular dynamics and free energy perturbation, we examine the relative binding affinity of several covalent polycyclic aromatic hydrocarbon - DNA (PAH-DNA) adducts at the central adenine of NRAS codon-61, a mutational hotspot implicated in cancer risk. Several PAHs classified by the International Agency for Research on Cancer as probable, possible, or unclassifiable as to carcinogenicity are found to have greater binding affinity than the known carcinogen, benzo[a]pyrene (B[a]P). van der Waals interactions between the intercalated PAH and neighboring nucleobases, and minimal disruption of the DNA duplex drive increases in binding affinity. PAH-DNA adducts may be repaired by global genomic nucleotide excision repair (GG-NER), hence we also compute relative free energies of complexation of PAH-DNA adducts with RAD4-RAD23 (the yeast ortholog of human XPC-RAD23) which constitutes the recognition step in GG-NER. PAH-DNA adducts exhibiting the greatest DNA binding affinity also exhibit the least RAD4-RAD23 complexation affinity and are thus predicted to resist the GG-NER machinery, contributing to their genotoxic potential. In particular, the fjord region PAHs dibenzo[a,l]pyrene, benzo[g]chrysene, and benzo[c]phenanthrene are found to have greater binding affinity while having weaker RAD4-RAD23 complexation affinity than their respective bay region analogs B[a]P, chrysene, and phenanthrene. We also find that the bay region PAHs dibenzo[a,j]anthracene, dibenzo[a,c]anthracene, and dibenzo[a,h]anthracene exhibit greater binding affinity and weaker RAD4-RAD23 complexation affinity than B[a]P. Thus, the study of PAH genotoxicity likely needs to be substantially broadened, with implications for public policy and the health sciences. This approach can be broadly applied to assess factors contributing to the genotoxicity of other unclassified compounds.

Automated Adaptive Absolute Binding Free Energy Calculations

Alchemical absolute binding free energy (ABFE) calculations have substantial potential in drug discovery, but are often prohibitively computationally expensive. To unlock their potential, efficient automated ABFE workflows are required to reduce both computational cost and human intervention. We present a fully automated ABFE workflow based on the automated selection of λ windows, the ensemble-based detection of equilibration, and the adaptive allocation of sampling time based on inter-replicate statistics. We find that the automated selection of intermediate states with consistent overlap is rapid, robust, and simple to implement. Robust detection of equilibration is achieved with a paired t-test between the free energy estimates at initial and final portions of a an ensemble of runs. We determine reasonable default parameters for all algorithms and show that the full workflow produces equivalent results to a nonadaptive scheme over a variety of test systems, while often accelerating equilibration. Our complete workflow is implemented in the open-source package A3FE (https://github.com/michellab/a3fe).

NPEX: Never give up protein exploration with deep reinforcement learning

Elucidating unknown structures of proteins, such as metastable states, is critical in designing therapeutic agents. Protein structure exploration has been performed using advanced computational methods, especially molecular dynamics and Markov chain Monte Carlo simulations, which require untenably long calculation times and prior structural knowledge. Here, we developed an innovative method for protein structure determination called never give up protein exploration (NPEX) with deep reinforcement learning. The NPEX method leverages the soft actor-critic algorithm and the intrinsic reward system, effectively adding a bias potential without the need for prior knowledge. To demonstrate the method's effectiveness, we applied it to four models: a double well, a triple well, the alanine dipeptide, and the tryptophan cage. Compared with Markov chain Monte Carlo simulations, NPEX had markedly greater sampling efficiency. The significantly enhanced computational efficiency and lack of prior domain knowledge requirements of the NPEX method will revolutionize protein structure exploration.

Caught between a ROCK and a hard place: current challenges in structure-based drug design

The discipline of structure-based drug design (SBDD) is several decades old and it is tempting to think that the proliferation of experimental structures for many drug targets might make computer-aided drug design (CADD) straightforward. However, this is far from true. In this review, we illustrate some of the challenges that CADD scientists face every day in their work, even now. We use Rho-associated protein kinase (ROCK), and public domain structures and data, as an example to illustrate some of the challenges we have experienced during our project targeting this protein. We hope that this will help to prevent unrealistic expectations of what CADD can accomplish and to educate non-CADD scientists regarding the challenges still facing their CADD colleagues.

Time-asymmetric fluctuation theorem and efficient free-energy estimation

The free-energy difference ΔF between two high-dimensional systems is notoriously difficult to compute but very important for many applications such as drug discovery. We demonstrate that an unconventional definition of work introduced by Vaikuntanathan and Jarzynski (2008) satisfies a microscopic fluctuation theorem that relates path ensembles that are driven by protocols unequal under time reversal. It has been shown before that counterdiabatic protocols-those having additional forcing that enforces the system to remain in instantaneous equilibrium, also known as escorted dynamics or engineered swift equilibration-yield zero-variance work measurements for this definition. We show that this time-asymmetric microscopic fluctuation theorem can be exploited for efficient free-energy estimation by developing a simple (i.e., neural-network free) and efficient adaptive time-asymmetric protocol optimization algorithm that yields ΔF estimates that are orders of magnitude lower in mean squared error than the generic linear interpolation protocol with which it is initialized.

Breaking the Gender Imbalance: Female Presence in the Computational Chemistry Group at Viva Biotech

The growing field of computational chemistry and artificial intelligence is drawing unprecedented attention in drug discovery. However, the underrepresentation of female scientists persists in the field. While the gender disparity is acknowledged in academic settings, the issue is less addressed and intensifies in the pharmaceutical and biotech industry both at the entry-level and leadership positions. At Viva Biotech, we challenge this norm: over 60% of the full-time employees in our computational chemistry group are women, a testament of our commitment to the pursuit of gender equality. We share our vision of the evolving role of computational chemistry in drug discovery and where women stand and rise in the field. In this article, we discuss how we engage female empowerment with tactical approaches and from personal experiences. The intention is to offer actionable guidance for female computational chemists at early career stages to bring visibility to their impacts, and ascent to senior positions.

Machine Learning Guided AQFEP: A Fast and Efficient Absolute Free Energy Perturbation Solution for Virtual Screening

Structure-based methods in drug discovery have become an integral part of the modern drug discovery process. The power of virtual screening lies in its ability to rapidly and cost-effectively explore enormous chemical spaces to select promising ligands for further experimental investigation. Relative free energy perturbation (RFEP) and similar methods are the gold standard for binding affinity prediction in drug discovery hit-to-lead and lead optimization phases, but have high computational cost and the requirement of a structural analog with a known activity. Without a reference molecule requirement, absolute FEP (AFEP) has, in theory, better accuracy for hit ID, but in practice, the slow throughput is not compatible with VS, where fast docking and unreliable scoring functions are still the standard. Here, we present an integrated workflow to virtually screen large and diverse chemical libraries efficiently, combining active learning with a physics-based scoring function based on a fast absolute free energy perturbation method. We validated the performance of the approach in the ranking of structurally related ligands, virtual screening hit rate enrichment, and active learning chemical space exploration; disclosing the largest reported collection of free energy simulations to date.

BAT2: an Open-Source Tool for Flexible, Automated, and Low Cost Absolute Binding Free Energy Calculations

Absolute binding free energy (ABFE) calculations with all-atom molecular dynamics (MD) have the potential to greatly reduce costs in the first stages of drug discovery. Here, we introduce BAT2, the new version of the Binding Affinity Tool (BAT.py), designed to combine full automation of ABFE calculations with high-performance MD simulations, making it a potential tool for virtual screening. We describe and test several changes and new features that were incorporated into the code, such as relative restraints between the protein and the ligand instead of using fixed dummy atoms, support for the OpenMM simulation engine, a merged approach to the application/release of restraints, support for cobinders and proteins with multiple chains, and many others. We also reduced the simulation times for each ABFE calculation, assessing the effect on the expected robustness and accuracy of the calculations.

Development of human lactate dehydrogenase a inhibitors: high-throughput screening, molecular dynamics simulation and enzyme activity assay

Lactate dehydrogenase A (LDHA) is highly expressed in many tumor cells and promotes the conversion of pyruvate to lactic acid in the glucose pathway, providing energy and synthetic precursors for rapid proliferation of tumor cells. Therefore, inhibition of LDHA has become a widely concerned tumor treatment strategy. However, the research and development of highly efficient and low toxic LDHA small molecule inhibitors still faces challenges. To discover potential inhibitors against LDHA, virtual screening based on molecular docking techniques was performed from Specs database of more than 260,000 compounds and Chemdiv-smart database of more than 1,000 compounds. Through molecular dynamics (MD) simulation studies, we identified 12 potential LDHA inhibitors, all of which can stably bind to human LDHA protein and form multiple interactions with its active central residues. In order to verify the inhibitory activities of these compounds, we established an enzyme activity assay system and measured their inhibitory effects on recombinant human LDHA. The results showed that Compound 6 could inhibit the catalytic effect of LDHA on pyruvate in a dose-dependent manner with an EC50 value of 14.54 ± 0.83 µM. Further in vitro experiments showed that Compound 6 could significantly inhibit the proliferation of various tumor cell lines such as pancreatic cancer cells and lung cancer cells, reduce intracellular lactic acid content and increase intracellular reactive oxygen species (ROS) level. In summary, through virtual screening and in vitro validation, we found that Compound 6 is a small molecule inhibitor for LDHA, providing a good lead compound for the research and development of LDHA related targeted anti-tumor drugs.

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Calculating the Gibbs energies of binding of ligand-receptor systems with a thermochemical accuracy of ± 1 kcal mol-1 is a challenge to computational approaches. After exploration of the conformational space of the host, ligand and their resulting complexes upon coordination by semi-empirical GFN2 MD and meta-MD simulations, the systematic refinement through a multi-level improvement of binding modes in terms of electronic energies and solvation is able to give Gibbs energies of binding of drug molecules to CB[8] and β-CD macrocyclic receptors with such an accuracy. The accurate treatment of a small number of structures outperforms system-specific force-matching and alchemical transfer model approaches without an extensive sampling and integration.

Identification of novel LRRK2 inhibitors by structure-based virtual screening and alchemical free energy calculation

The Leucine-rich repeat kinase 2 (LRRK2) target has been identified as a promising drug target for Parkinson's disease (PD) treatment. This study focuses on optimizing the activity of LRRK2 inhibitors using alchemical relative binding free energy (RBFE) calculations. Initially, we assessed various free energy calculation methods across different LRRK2 kinase inhibitor scaffolds. The results indicate that alchemical free energy calculations are promising for prospective predictions on LRRK2 inhibitors, especially for the aminopyrimidine scaffold with an RMSE of 1.15 kcal mol-1 and Rp of 0.83. Following this, we optimized a potent LRRK2 kinase inhibitor identified from previous virtual screenings, featuring a novel scaffold. Guided by RBFE predictions using alchemical methods, this optimization led to the discovery of compound LY2023-001. This compound, with a [1,2,4]triazolo[5,6-b]indole scaffold, exhibited enhanced inhibitory activity against G2019S LRRK2 (IC50 = 12.9 nM). Molecular dynamics (MD) simulations revealed that LY2023-001 formed stable hydrogen bonds with Glu1948, and Ala1950 in the G2019S LRRK2 protein. Additionally, its phenyl substituents engage in strong electrostatic interactions with Lys1906 and van der Waals interactions with Leu1885, Phe1890, Val1893, Ile1933, Met1947, Leu1949, Leu2001, Ala2016, and Asp2017. Our findings underscore the potential of computational methods in the successful optimization of small molecules, offering important insights for the development of novel LRRK2 inhibitors.

Rapid, Accurate, Ranking of Protein-Ligand Binding Affinities with VM2, the Second-Generation Mining Minima Method

The structure-based technologies most widely used to rank the affinities of candidate small molecule drugs for proteins range from faster but less reliable docking methods to slower but more accurate explicit solvent free energy methods. In recent years, we have advanced another technology, which is called mining minima because it "mines" out the main contributions to the chemical potentials of the free and bound molecular species by identifying and characterizing their main local energy minima. The present study provides systematic benchmarks of the accuracy and computational speed of mining minima, as implemented in the VeraChem Mining Minima Generation 2 (VM2) code, across two well-regarded protein-ligand benchmark data sets, for which there are already benchmark data for docking, free energy, and other computational methods. A core result is that VM2's accuracy approaches that of explicit solvent free energy methods at a far lower computational cost. In finer-grained analyses, we also examine the influence of various run settings, such as the treatment of crystallographic water molecules, on the accuracy, and define the costs in time and dollars of representative runs on Amazon Web Services (AWS) compute instances with various CPU and GPU combinations. We also use the benchmark data to determine the importance of VM2's correction from generalized Born to finite-difference Poisson-Boltzmann results for each energy well and find that this correction affords a remarkably consistent improvement in accuracy at a modest computational cost. The present results establish VM2 as a distinctive technology for early-stage drug discovery, which provides a strong combination of efficiency and predictivity.

Do All Paths Lead to Rome? How Reliable is Umbrella Sampling Along a Single Path?

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

STORMM: Structure and topology replica molecular mechanics for chemical simulations

The Structure and TOpology Replica Molecular Mechanics (STORMM) code is a next-generation molecular simulation engine and associated libraries optimized for performance on fast, vectorized central processor units and graphics processing units (GPUs) with independent memory and tens of thousands of threads. STORMM is built to run thousands of independent molecular mechanical calculations on a single GPU with novel implementations that tune numerical precision, mathematical operations, and scarce on-chip memory resources to optimize throughput. The libraries are built around accessible classes with detailed documentation, supporting fine-grained parallelism and algorithm development as well as copying or swapping groups of systems on and off of the GPU. A primary intention of the STORMM libraries is to provide developers of atomic simulation methods with access to a high-performance molecular mechanics engine with extensive facilities to prototype and develop bespoke tools aimed toward drug discovery applications. In its present state, STORMM delivers molecular dynamics simulations of small molecules and small proteins in implicit solvent with tens to hundreds of times the throughput of conventional codes. The engineering paradigm transforms two of the most memory bandwidth-intensive aspects of condensed-phase dynamics, particle-mesh mapping, and valence interactions, into compute-bound problems for several times the scalability of existing programs. Numerical methods for compressing and streamlining the information present in stored coordinates and lookup tables are also presented, delivering improved accuracy over methods implemented in other molecular dynamics engines. The open-source code is released under the MIT license.

Insights and Challenges in Correcting Force Field Based Solvation Free Energies Using a Neural Network Potential

We present a comprehensive study investigating the potential gain in accuracy for calculating absolute solvation free energies (ASFE) using a neural network potential to describe the intramolecular energy of the solute. We calculated the ASFE for most compounds from the FreeSolv database using the Open Force Field (OpenFF) and compared them to earlier results obtained with the CHARMM General Force Field (CGenFF). By applying a nonequilibrium (NEQ) switching approach between the molecular mechanics (MM) description (either OpenFF or CGenFF) and the neural net potential (NNP)/MM level of theory (using ANI-2x as the NNP potential), we attempted to improve the accuracy of the calculated ASFEs. The predictive performance of the results did not change when this approach was applied to all 589 small molecules in the FreeSolv database that ANI-2x can describe. When selecting a subset of 156 molecules, focusing on compounds where the force fields performed poorly, we saw a slight improvement in the root-mean-square error (RMSE) and mean absolute error (MAE). The majority of our calculations utilized unidirectional NEQ protocols based on Jarzynski's equation. Additionally, we conducted bidirectional NEQ switching for a subset of 156 solutes. Notably, only a small fraction (10 out of 156) exhibited statistically significant discrepancies between unidirectional and bidirectional NEQ switching free energy estimates.

MDFit: automated molecular simulations workflow enables high throughput assessment of ligands-protein dynamics

Molecular dynamics (MD) simulation is a powerful tool for characterizing ligand-protein conformational dynamics and offers significant advantages over docking and other rigid structure-based computational methods. However, setting up, running, and analyzing MD simulations continues to be a multi-step process making it cumbersome to assess a library of ligands in a protein binding pocket using MD. We present an automated workflow that streamlines setting up, running, and analyzing Desmond MD simulations for protein-ligand complexes using machine learning (ML) models. The workflow takes a library of pre-docked ligands and a prepared protein structure as input, sets up and runs MD with each protein-ligand complex, and generates simulation fingerprints for each ligand. Simulation fingerprints (SimFP) capture protein-ligand compatibility, including stability of different ligand-pocket interactions and other useful metrics that enable easy rank-ordering of the ligand library for pocket optimization. SimFPs from a ligand library are used to build & deploy ML models that predict binding assay outcomes and automatically infer important interactions. Unlike relative free-energy methods that are constrained to assess ligands with high chemical similarity, ML models based on SimFPs can accommodate diverse ligand sets. We present two case studies on how SimFP helps delineate structure-activity relationship (SAR) trends and explain potency differences across matched-molecular pairs of (1) cyclic peptides targeting PD-L1 and (2) small molecule inhibitors targeting CDK9.

Adsorption of Flavonoids in a Transcriptional Regulator TtgR: Relative Binding Free Energies and Intermolecular Interactions

Antimicrobial resistance in bacteria often arises from their ability to actively identify and expel toxic compounds. The bacterium strain Pseudomonas putida DOT-T1E utilizes its TtgABC efflux pump to confer robust resistance against antibiotics, flavonoids, and organic solvents. This resistance mechanism is intricately regulated at the transcriptional level by the TtgR protein. Through molecular dynamics and alchemical free energy simulations, we systematically examine the binding of seven flavonoids and their derivatives with the TtgR transcriptional regulator. Our simulations reveal distinct binding geometries and free energies for the flavonoids in the active site of the protein, which are driven by a range of noncovalent forces encompassing van der Waals, electrostatic, and hydrogen bonding interactions. The interplay of molecular structures, substituent patterns, and intermolecular interactions effectively stabilizes the bound flavonoids, confining their movements within the TtgR binding pocket. These findings yield valuable insights into the molecular determinants that govern ligand recognition in TtgR and shed light on the mechanism of antimicrobial resistance in P. putida DOT-T1E.

Computational Investigation of BMAA and Its Carbamate Adducts as Potential GluR2 Modulators

Beta-N-methylamino-l-alanine (BMAA) is a potential neurotoxic nonprotein amino acid, which can reach the human body through the food chain. When BMAA interacts with bicarbonate in the human body, carbamate adducts are produced, which share a high structural similarity with the neurotransmitter glutamate. It is believed that BMAA and its l-carbamate adducts bind in the glutamate binding site of ionotropic glutamate receptor 2 (GluR2). Chronic exposure to BMAA and its adducts could cause neurological illness such as neurodegenerative diseases. However, the mechanism of BMAA action and its carbamate adducts bound to GluR2 has not yet been elucidated. Here, we investigate the binding modes and the affinity of BMAA and its carbamate adducts to GluR2 in comparison to the natural agonist, glutamate, to understand whether these can act as GluR2 modulators. Initially, we perform molecular dynamics simulations of BMAA and its carbamate adducts bound to GluR2 to examine the stability of the ligands in the S1/S2 ligand-binding core of the receptor. In addition, we utilize alchemical free energy calculations to compute the difference in the free energy of binding of the beta-carbamate adduct of BMAA to GluR2 compared to that of glutamate. Our findings indicate that carbamate adducts of BMAA and glutamate remain stable in the binding site of the GluR2 compared to BMAA. Additionally, alchemical free energy results reveal that glutamate and the beta-carbamate adduct of BMAA have comparable binding affinity to the GluR2. These results provide a rationale that BMAA carbamate adducts may be, in fact, the modulators of GluR2 and not BMAA itself.

Fluorinated Protein-Ligand Complexes: A Computational Perspective

Fluorine is an element renowned for its unique properties. Its powerful capability to modulate molecular properties makes it an attractive substituent for protein binding ligands; however, the rational design of fluorination can be challenging with effects on interactions and binding energies being difficult to predict. In this Perspective, we highlight how computational methods help us to understand the role of fluorine in protein-ligand binding with a focus on molecular simulation. We underline the importance of an accurate force field, present fluoride channels as a showcase for biomolecular interactions with fluorine, and discuss fluorine specific interactions like the ability to form hydrogen bonds and interactions with aryl groups. We put special emphasis on the disruption of water networks and entropic effects.

Targeting the conserved active site of splicing machines with specific and selective small molecule modulators

The self-splicing group II introns are bacterial and organellar ancestors of the nuclear spliceosome and retro-transposable elements of pharmacological and biotechnological importance. Integrating enzymatic, crystallographic, and simulation studies, we demonstrate how these introns recognize small molecules through their conserved active site. These RNA-binding small molecules selectively inhibit the two steps of splicing by adopting distinctive poses at different stages of catalysis, and by preventing crucial active site conformational changes that are essential for splicing progression. Our data exemplify the enormous power of RNA binders to mechanistically probe vital cellular pathways. Most importantly, by proving that the evolutionarily-conserved RNA core of splicing machines can recognize small molecules specifically, our work provides a solid basis for the rational design of splicing modulators not only against bacterial and organellar introns, but also against the human spliceosome, which is a validated drug target for the treatment of congenital diseases and cancers.

Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation

We introduce the lambda-Adaptive Biasing Force (lambda-ABF) method for the computation of alchemical free-energy differences. We propose a software implementation and showcase it on biomolecular systems. The method arises from coupling multiple-walker adaptive biasing force with λ-dynamics. The sampling of the alchemical variable is continuous and converges toward a uniform distribution, making manual optimization of the λ schedule unnecessary. Contrary to most other approaches, alchemical free-energy estimates are obtained immediately without any postprocessing. Free diffusion of λ improves orthogonal relaxation compared to fixed-λ thermodynamic integration or free-energy perturbation. Furthermore, multiple walkers provide generic orthogonal space coverage with minimal user input and negligible computational overhead. We show that our high-performance implementations coupling the Colvars library with NAMD and Tinker-HP can address real-world cases including ligand-receptor binding with both fixed-charge and polarizable models, with a demonstrably richer sampling than fixed-λ methods. The implementation is fully open-source, publicly available, and readily usable by practitioners of current alchemical methods. Thanks to the portable Colvars library, lambda-ABF presents a unified user interface regardless of the back-end (NAMD, Tinker-HP, or any software to be interfaced in the future), sparing users the effort of learning multiple interfaces. Finally, the Colvars Dashboard extension of the visual molecular dynamics (VMD) software provides an interactive monitoring and diagnostic tool for lambda-ABF simulations.

Understanding the impact of binding free energy and kinetics calculations in modern drug discovery

INTRODUCTION

For rational drug design, it is crucial to understand the receptor-drug binding processes and mechanisms. A new era for the use of computer simulations in predicting drug-receptor interactions at an atomic level has begun with remarkable advances in supercomputing and methodological breakthroughs.

AREAS COVERED

End-point free energy calculation methods such as Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) or Molecular-Mechanics/Generalized Born Surface Area (MM/GBSA), free energy perturbation (FEP), and thermodynamic integration (TI) are commonly used for binding free energy calculations in drug discovery. In addition, kinetic dissociation and association rate constants (k off and k on ) play critical roles in the function of drugs. Nowadays, Molecular Dynamics (MD) and enhanced sampling simulations are increasingly being used in drug discovery. Here, the authors provide a review of the computational techniques used in drug binding free energy and kinetics calculations.

EXPERT OPINION

The applications of computational methods in drug discovery and design are expanding, thanks to improved predictions of the binding free energy and kinetic rates of drug molecules. Recent microsecond-timescale enhanced sampling simulations have made it possible to accurately capture repetitive ligand binding and dissociation, facilitating more efficient and accurate calculations of ligand binding free energy and kinetics.

Subtractive proteomics-guided vaccine targets identification and designing of multi-epitopes vaccine for immune response instigation against Burkholderia pseudomallei

In this study, a methodical workflow using subtractive proteomics, vaccine designing, molecular simulation, and agent-based modeling approaches were used to annotate the whole proteome of Burkholderia pseudomallei (strain K96243) for vaccine designing. Among the total 5717 proteins in the whole proteome, 505 were observed to be essential for the pathogen's survival and pathogenesis predicted by the Database of Essential Genes. Among these, 23 vaccine targets were identified, of which fimbrial assembly chaperone (Q63UH5), Outer membrane protein (Q63UH1), and Hemolysin-like protein (Q63UE4) were selected for the subsequent analysis based on the systematic approaches. Using immunoinformatic approaches CTL (cytotoxic T lymphocytes), HTL (helper T lymphocytes), IFN-positive, and B cell epitopes were predicted for these targets. A total of 9 CTL epitopes were added using the GSS linker, 6 HTL epitopes using the GPGPG linker, and 6 B cell epitopes using the KK linker. An adjuvant was added for enhanced antigenicity, an HIV-TAT peptide for improved delivery, and a PADRE sequence was added to form a 466 amino acids long vaccine construct. The construct was classified as non-allergenic, highly antigenic, and experimentally feasible. Molecular docking results validated the robust interaction of MEVC with immune receptors such as TLR2/4. Furthermore, molecular simulation revealed stable dynamics and compact nature of the complexes. The binding free energy results further validated the robust binding. In silico cloning, results revealed GC contents of 50.73 % and a CIA value of 0.978 which shows proper downstream processing. Immune simulation results reported that after the three injections of the vaccine a robust secondary immune response, improved antigen clearance, and effective immune memory generation were observed highlighting its potential for effective and sustained immunity. Future directions should encompass experimental validations, animal model studies, and clinical trials to substantiate the vaccine's efficacy, safety, and immunogenicity.

Calculating Binding Free Energies in Model Host-Guest Systems with Unrestrained Advanced Sampling

Host-guest interactions are important to the design of pharmaceuticals and, more broadly, to soft materials as they can enable targeted, strong, and specific interactions between molecules. The binding process between the host and guest may be classified as a "rare event" when viewing the system at atomic scales, such as those explored in molecular dynamics simulations. To obtain equilibrium binding conformations and dissociation constants from these simulations, it is essential to resolve these rare events. Advanced sampling methods such as the adaptive biasing force (ABF) promote the occurrence of less probable configurations in a system, therefore facilitating the sampling of essential collective variables that characterize the host-guest interactions. Here, we present the application of ABF to a rod-cavitand coarse-grained model of host-guest systems to acquire the potential of mean force. We show that the employment of ABF enables the computation of the configurational and thermodynamic properties of bound and unbound states, including the free energy landscape. Moreover, we identify important dynamic bottlenecks that limit sampling and discuss how these may be addressed in more general systems.

Alchemical Enhanced Sampling with Optimized Phase Space Overlap

An alchemical enhanced sampling (ACES) method has recently been introduced to facilitate importance sampling in free energy simulations. The method achieves enhanced sampling from Hamiltonian replica exchange within a dual topology framework while utilizing new smoothstep softcore potentials. A common sampling problem encountered in lead optimization is the functionalization of aromatic rings that exhibit distinct conformational preferences when interacting with the protein. It is difficult to converge the distribution of ring conformations due to the long time scale of ring flipping events; however, the ACES method addresses this issue by modeling the syn and anti ring conformations within a dual topology. ACES thereby samples the conformer distributions by alchemically tunneling between states, as opposed to traversing a physical pathway with a high rotational barrier. We demonstrate the use of ACES to overcome conformational sampling issues involving ring flipping in ML300-derived noncovalent inhibitors of SARS-CoV-2 Main Protease (Mpro). The demonstrations explore how the use of replica exchange and the choice of softcore selection affects the convergence of the ring conformation distributions. Furthermore, we examine how the accuracy of the calculated free energies is affected by the degree of phase space overlap (PSO) between adjacent states (i.e., between neighboring λ-windows) and the Hamiltonian replica exchange acceptance ratios. Both of these factors are sensitive to the spacing between the intermediate states. We introduce a new method for choosing a schedule of λ values. The method analyzes short "burn-in" simulations to construct a 2D map of the nonlocal PSO. The schedule is obtained by optimizing an alchemical pathway on the 2D map that equalizes the PSO between the λ intervals. The optimized phase space overlap λ-spacing method (Opt-PSO) leads to more numerous end-to-end single passes and round trips due to the correlation between PSO and Hamiltonian replica exchange acceptance ratios. The improved exchange statistics enhance the efficiency of ACES method. The method has been implemented into the FE-ToolKit software package, which is freely available.

FEP Augmentation as a Means to Solve Data Paucity Problems for Machine Learning in Chemical Biology

In the realm of medicinal chemistry, the primary objective is to swiftly optimize a multitude of chemical properties of a set of compounds to yield a clinical candidate poised for clinical trials. In recent years, two computational techniques, machine learning (ML) and physics-based methods, have evolved substantially and are now frequently incorporated into the medicinal chemist's toolbox to enhance the efficiency of both hit optimization and candidate design. Both computational methods come with their own set of limitations, and they are often used independently of each other. ML's capability to screen extensive compound libraries expediently is tempered by its reliance on quality data, which can be scarce especially during early-stage optimization. Contrarily, physics-based approaches like free energy perturbation (FEP) are frequently constrained by low throughput and high cost by comparison; however, physics-based methods are capable of making highly accurate binding affinity predictions. In this study, we harnessed the strength of FEP to overcome data paucity in ML by generating virtual activity data sets which then inform the training of algorithms. Here, we show that ML algorithms trained with an FEP-augmented data set could achieve comparable predictive accuracy to data sets trained on experimental data from biological assays. Throughout the paper, we emphasize key mechanistic considerations that must be taken into account when aiming to augment data sets and lay the groundwork for successful implementation. Ultimately, the study advocates for the synergy of physics-based methods and ML to expedite the lead optimization process. We believe that the physics-based augmentation of ML will significantly benefit drug discovery, as these techniques continue to evolve.

Adaptive lambda schemes for efficient relative binding free energy calculation

The relative free energy perturbation (RFEP) calculation is one of the most theoretically sound computational chemistry approaches for the binding affinity prediction. However, its application is often hindered by the complexity of the calculation choices and the requirement of expertise in the field. Improper lambda scheme of RFEP may result in deviations from an accurate description of the perturbation process and is prone to erroneous affinity predictions. To address such challenges, an automated adaptive lambda method is proposed where the adaptive lambda schemes are obtained through a split-and-merge algorithm based on the pilot runs. The newly established workflow along with a series of improvements to the perturbation settings increases the consistency of the RFEP calculation results. Comparing the pilot and adaptive lambda schemes, the latter demonstrated improvements in convergence and reproducibility and lowered the mean unsigned error and the root-mean-square error. Overall, the adaptive lambda method is a reliable and robust choice to predict small molecule relative binding free energy and can be capitalized to benefit routine RFEP calculations for drug discovery projects.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

BioProtIS: Streamlining protein-ligand interaction pipeline for analysis in genomic and transcriptomic exploration

The identification of protein-ligand interactions plays a pivotal role in elucidating biological processes and discovering potential bioproducts. Harnessing the capabilities of computational methods in drug discovery, we introduce an innovative Inverted Virtual Screening (IVS) pipeline. This pipeline Integrated molecular dynamics and docking analyses to ensure that protein structures are not only energetically favorable but also representative of stable conformations. The primary objective of this pipeline is to automate and streamline the analysis of protein-ligand interactions at both genomic and transcriptomic scales. In the contemporary post-genomic era, high-throughput computational screening for bioproducts, biological systems, and therapeutic drugs has become a cornerstone practice. This approach offers the promise of cost-effectiveness, time efficiency, and optimization of laboratory work. Nevertheless, a notable deficiency persists in the availability of efficient pipelines capable of automating the virtual screening process, seamlessly integrating input and output, and leveraging the full potential of open-source tools. To bridge this critical gap, we have developed a versatile pipeline known as BioProtIS. This tool seamlessly integrates a suite of state-of-the-art tools, including Modeller, AlphaFold, Gromacs, FPOCKET, and AutoDock Vina, thus facilitating the streamlined docking of ligands with an expansive repertoire of proteins sourced from genomes and transcriptomes, and substrates. To assess the pipeline's performance, we employed the transcriptomes of Cereus jamacaru (a cactus species) and Aspisoma lineatum (firefly), along with the genome of Homo sapiens. This integration not only improves the accuracy of ligand-protein interactions by minimizing replicability deviations but also optimizes the discovery process by enabling the simultaneous evaluation of multiple substrates. Furthermore, our pipeline accommodates distinct testing scenarios, such as blind docking or site-specific targeting, which are invaluable in applications ranging from drug repositioning to the exploration of new allosteric binding sites and toxicity assessments. BioProtIS has been designed with modularity at its core. This inherent flexibility empowers users to make custom modifications directly within the source code, tailoring the pipeline to their specific research needs. Moreover, it lays the foundation for seamless integration of diverse docking algorithms in future iterations, promising ongoing advancements in the field of computational biology. This pipeline is available for free distribution and can be download at: https://github.com/BBMDO/BioProtIS.

Elucidating the Mechanism of Metabolism of Cannabichromene by Human Cytochrome P450s

Cannabichromene (CBC) is a nonpsychoactive phytocannabinoid well-known for its wide-ranging health advantages. However, there is limited knowledge regarding its human metabolism following CBC consumption. This research aimed to explore the metabolic pathways of CBC by various human liver cytochrome P450 (CYP) enzymes and support the outcomes using in vivo data from mice. The results unveiled two principal CBC metabolites generated by CYPs: 8'-hydroxy-CBC and 6',7'-epoxy-CBC, along with a minor quantity of 1″-hydroxy-CBC. Notably, among the examined CYPs, CYP2C9 demonstrated the highest efficiency in producing these metabolites. Moreover, through a molecular dynamics simulation spanning 1 μs, it was observed that CBC attains stability at the active site of CYP2J2 by forming hydrogen bonds with I487 and N379, facilitated by water molecules, which specifically promotes the hydroxy metabolite's formation. Additionally, the presence of cytochrome P450 reductase (CPR) amplified CBC's binding affinity to CYPs, particularly with CYP2C8 and CYP3A4. Furthermore, the metabolites derived from CBC reduced cytokine levels, such as IL6 and NO, by approximately 50% in microglia cells. This investigation offers valuable insights into the biotransformation of CBC, underscoring the physiological importance and the potential significance of these metabolites.

LC-MS/MS profiling of Tipuana tipu flower, HPLC-DAD quantification of its bioactive components, and interrelationships with antioxidant, and anti-inflammatory activity: in vitro and in silico approaches

BACKGROUND

Fabaceae plays a crucial role in African traditional medicine as a source of large number of important folk medication, agriculture and food plants. In a search of potential antioxidant and anti-inflammatory candidates derived from locally cultivated plants, the flowers of Tipuana tipu (Benth.) Lillo growing in Egypt were subjected to extensive biological and phytochemical studies. The impact of the extraction technique on the estimated biological activities was investigated.

METHODS

The flowers were extracted using different solvents (aqueous, methanol, water/methanol (1:1), methanol/methylene chloride (1:1), and methylene chloride). The different extracts were subjected to antioxidant (DPPH, ABTS, and FRAP) and anti-inflammatory (COX-2 and 5-LOX) assays. The methanol extract was assessed for its inhibitory activity against iNOS, NO production, and pro-inflammatory cytokines (NF-KB, TNF-R2, TNF-α, IL-1β, and IL-6) in LPS-activated RAW 264.7 macrophages. The composition-activity relationship of the active methanol extract was further investigated using a comprehensive LC-QTOF-MS/MS analysis. The major identified phenolic compounds were further quantified using HPLC-DAD technique. The affinity of representative compounds to iNOS, COX-2, and 5-LOX target active sites was investigated using molecular docking and molecular dynamics simulations.

RESULTS

The methanol extract exhibited the highest radical scavenging capacity and enzyme inhibitory activities against COX-2 and 5-LOX enzymes with IC50 values of 10.6 ± 0.4 and 14.4 ± 1.0 µg/mL, respectively. It also inhibited iNOS enzyme activity, suppressed NO production, and decreased the secretion of pro-inflammatory cytokines. In total, 62 compounds were identified in the extract including flavonoids, coumarins, organic, phenolic, and fatty acids. Among them 18 phenolic compounds were quantified by HPLC-DAD. The highest docking scores were achieved by kaempferol-3-glucoside and orientin. Additionally, molecular dynamics simulations supported the docking findings.

CONCLUSION

The flower could be considered a potentially valuable component in herbal medicines owing to its unique composition and promising bioactivities. These findings encourage increased propagation of T. tipu or even tissue culturing of its flowers for bioprospecting of novel anti-inflammatory drugs. Such applications could be adopted as future approaches that benefit the biomedical field.

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the λ-exchange of free energy perturbation (FEP). Two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/λ-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

Binding Sites of Bicarbonate in Phosphoenolpyruvate Carboxylase

Phosphoenolpyruvate carboxylase (PEPC) is used in plant metabolism for fruit maturation or seed development as well as in the C4 and crassulacean acid metabolism (CAM) mechanisms in photosynthesis, where it is used for the capture of hydrated CO2 (bicarbonate). To find the yet unknown binding site of bicarbonate in this enzyme, we have first identified putative binding sites with nonequilibrium molecular dynamics simulations and then ranked these sites with alchemical free energy calculations with corrections of computational artifacts. Fourteen pockets where bicarbonate could bind were identified, with three having realistic binding free energies with differences with the experimental value below 1 kcal/mol. One of these pockets is found far from the active site at 14 Å and predicted to be an allosteric binding site. In the two other binding sites, bicarbonate is in direct interaction with the magnesium ion; neither sequence alignment nor the study of mutant K606N allowed to discriminate between these two pockets, and both are good candidates as the binding site of bicarbonate in phosphoenolpyruvate carboxylase.

Water Network-Augmented Two-State Model for Protein-Ligand Binding Affinity Prediction

Water network rearrangement from the ligand-unbound state to the ligand-bound state is known to have significant effects on the protein-ligand binding interactions, but most of the current machine learning-based scoring functions overlook these effects. In this study, we endeavor to construct a comprehensive and realistic deep learning model by incorporating water network information into both ligand-unbound and -bound states. In particular, extended connectivity interaction features were integrated into graph representation, and graph transformer operator was employed to extract features of the ligand-unbound and -bound states. Through these efforts, we developed a water network-augmented two-state model called ECIFGraph::HM-Holo-Apo. Our new model exhibits satisfactory performance in terms of scoring, ranking, docking, screening, and reverse screening power tests on the CASF-2016 benchmark. In addition, it can achieve superior performance in large-scale docking-based virtual screening tests on the DEKOIS2.0 data set. Our study highlights that the use of a water network-augmented two-state model can be an effective strategy to bolster the robustness and applicability of machine learning-based scoring functions, particularly for targets with hydrophilic or solvent-exposed binding pockets.

Directional ΔG Neural Network (DrΔG-Net): A Modular Neural Network Approach to Binding Free Energy Prediction

The protein-ligand binding free energy is a central quantity in structure-based computational drug discovery efforts. Although popular alchemical methods provide sound statistical means of computing the binding free energy of a large breadth of systems, they are generally too costly to be applied at the same frequency as end point or ligand-based methods. By contrast, these data-driven approaches are typically fast enough to address thousands of systems but with reduced transferability to unseen systems. We introduce DrΔG-Net (or simply Dragnet), an equivariant graph neural network that can blend ligand-based and protein-ligand data-driven approaches. It is based on a 3D fingerprint representation of the ligand alone and in complex with the protein target. Dragnet is a global scoring function to predict the binding affinity of arbitrary protein-ligand complexes, but can be easily tuned via transfer learning to specific systems or end points, performing similarly to common 2D ligand-based approaches in these tasks. Dragnet is evaluated on a total of 28 validation proteins with a set of congeneric ligands derived from the Binding DB and one custom set extracted from the ChEMBL Database. In general, a handful of experimental binding affinities are sufficient to optimize the scoring function for a particular protein and ligand scaffold. When not available, predictions from physics-based methods such as absolute free energy perturbation can be used for the transfer learning tuning of Dragnet. Furthermore, we use our data to illustrate the present limitations of data-driven modeling of binding free energy predictions.

Enhancing Generalizability in Protein-Ligand Binding Affinity Prediction with Multimodal Contrastive Learning

Improving the generalization ability of scoring functions remains a major challenge in protein-ligand binding affinity prediction. Many machine learning methods are limited by their reliance on single-modal representations, hindering a comprehensive understanding of protein-ligand interactions. We introduce a graph-neural-network-based scoring function that utilizes a triplet contrastive learning loss to improve protein-ligand representations. In this model, three-dimensional complex representations and the fusion of two-dimensional ligand and coarse-grained pocket representations converge while distancing from decoy representations in latent space. After rigorous validation on multiple external data sets, our model exhibits commendable generalization capabilities compared to those of other deep learning-based scoring functions, marking it as a promising tool in the realm of drug discovery. In the future, our training framework can be extended to other biophysical- and biochemical-related problems such as protein-protein interaction and protein mutation prediction.

Kartograf: A Geometrically Accurate Atom Mapper for Hybrid-Topology Relative Free Energy Calculations

Relative binding free energy (RBFE) calculations have emerged as a powerful tool that supports ligand optimization in drug discovery. Despite many successes, the use of RBFEs can often be limited by automation problems, in particular, the setup of such calculations. Atom mapping algorithms are an essential component in setting up automatic large-scale hybrid-topology RBFE calculation campaigns. Traditional algorithms typically employ a 2D subgraph isomorphism solver (SIS) in order to estimate the maximum common substructure. SIS-based approaches can be limited by time-intensive operations and issues with capturing geometry-linked chemical properties, potentially leading to suboptimal solutions. To overcome these limitations, we have developed Kartograf, a geometric-graph-based algorithm that uses primarily the 3D coordinates of atoms to find a mapping between two ligands. In free energy approaches, the ligand conformations are usually derived from docking or other previous modeling approaches, giving the coordinates a certain importance. By considering the spatial relationships between atoms related to the molecule coordinates, our algorithm bypasses the computationally complex subgraph matching of SIS-based approaches and reduces the problem to a much simpler bipartite graph matching problem. Moreover, Kartograf effectively circumvents typical mapping issues induced by molecule symmetry and stereoisomerism, making it a more robust approach for atom mapping from a geometric perspective. To validate our method, we calculated mappings with our novel approach using a diverse set of small molecules and used the mappings in relative hydration and binding free energy calculations. The comparison with two SIS-based algorithms showed that Kartograf offers a fast alternative approach. The code for Kartograf is freely available on GitHub (https://github.com/OpenFreeEnergy/kartograf). While developed for the OpenFE ecosystem, Kartograf can also be utilized as a standalone Python package.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Enthalpies and entropies of hydration from Monte Carlo simulations

The changes in free energy, enthalpy, and entropy for transfer of a solute from the gas phase into solution are the fundamental thermodynamic quantities that characterize the solvation process. Owing to the development of methods based on free-energy perturbation theory, computation of free energies of solvation has become routine in conjunction with Monte Carlo (MC) statistical mechanics and molecular dynamics (MD) simulations. Computation of the enthalpy change and by inference the entropy change is more challenging. Two methods are considered in this work corresponding to direct averaging for the solvent and solution and to computing the temperature derivative of the free energy in the van't Hoff approach. The application is for neutral organic solutes in TIP4P water using long MC simulations to improve precision. Definitive results are also provided for pure TIP4P water. While the uncertainty in computed free energies of hydration is ca. 0.05 kcal mol-1, it is ca. 0.4 kcal mol-1 for the enthalpy changes from either van't Hoff plots or the direct method with sampling for 5 billion MC configurations. Partial molar volumes of hydration are also computed by the direct method; they agree well with experimental data with an average deviation of 3 cm3 mol-1. In addition, the results permit breakdown of the errors in the free energy changes from the OPLS-AA force field into their enthalpic and entropic components. The excess hydrophobicity of organic solutes is enthalpic in origin.