A highly accurate metadynamics-based Dissociation Free Energy method to calculate protein-protein and protein-ligand binding potencies

A Mode Evolution Metric to Extract Reaction Coordinates for Biomolecular Conformational Transitions

The complex, multidimensional energy landscape of biomolecules makes the extraction of suitable, nonintuitive collective variables (CVs) that describe their conformational transitions challenging. At present, dimensionality reduction approaches and machine learning (ML) schemes are employed to obtain CVs from molecular dynamics (MD)/Monte Carlo (MC) trajectories or structural databanks for biomolecules. However, minimum sampling conditions to generate reliable CVs that accurately describe the underlying energy landscape remain unclear. Here, we address this issue by developing a Mode evolution Metric (MeM) to extract CVs that can pinpoint new states and describe local transitions in the vicinity of a reference minimum from nonequilibrated MD/MC trajectories. We present a general mathematical formulation of MeM for both statistical dimensionality reduction and machine learning approaches. Application of MeM to MC trajectories of model potential energy landscapes and MD trajectories of solvated alanine dipeptide reveals that the principal components which locate new states in the vicinity of a reference minimum emerge well before the trajectories locally equilibrate between the associated states. Finally, we demonstrate a possible application of MeM in designing efficient biased sampling schemes to construct accurate energy landscape slices that link transitions between states. MeM can help speed up the search for new minima around a biomolecular conformational state and enable the accurate estimation of thermodynamics for states lying on the energy landscape and the description of associated transitions.

Computational Insights into Papaveroline as an In Silico Drug Candidate for Alzheimer's Disease via Fyn Tyrosine Kinase Inhibition

Alzheimer's disease (AD) poses a significant global health challenge, necessitating the exploration of novel therapeutic strategies. Fyn Tyrosine Kinase has emerged as a key player in AD pathogenesis, making it an attractive target for drug development. This study focuses on investigating the potential of Papaveroline as a drug candidate for AD by targeting Fyn Tyrosine Kinase. The research employed high-throughput virtual screening and QSAR analysis were conducted to identify compounds with optimal drug-like properties, emphasizing adherence to ADMET parameters for further evaluation. Molecular dynamics simulations to analyze the binding interactions between Papaveroline and Staurosporine with Fyn Tyrosine Kinase over a 200-ns period. The study revealed detailed insights into the binding mechanisms and stability of the Papaveroline-Fyn complex, showcasing the compound's potential as an inhibitor of Fyn Tyrosine Kinase. Comparative analysis with natural compounds and a reference compound highlighted Papaveroline's unique characteristics and promising therapeutic implications for AD treatment. Overall, the findings underscore Papaveroline's potential as a valuable drug candidate for targeting Fyn Tyrosine Kinase in AD therapy, offering new avenues for drug discovery in neurodegenerative diseases. This study contributes to advancing our understanding of molecular interactions in AD pathogenesis and paves the way for further research and development in this critical area.

Targeting fusion oncoproteins in childhood cancers: challenges and future opportunities for developing therapeutics

Fusion oncoproteins are associated with childhood cancers and have proven challenging to target, aside from those that include kinases. As part of its efforts for targeting childhood cancers, the National Cancer Institute recently conducted a series on Novel Chemical Approaches for Targeting Fusion Oncoproteins. Key learnings on leading platforms and technologies that can be used to advance the development of molecular therapeutics that target fusion oncoproteins in childhood cancers are described. Recent breakthroughs in medicinal chemistry and chemical biology provide new ground and creative strategies to exploit for the development of targeted agents for improving outcomes against these recalcitrant cancers.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media

Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.

Molecular Dynamics Simulations in Protein-Protein Docking

Concerted interactions between all the cell components form the basis of biological processes. Protein-protein interactions (PPIs) constitute a tremendous part of this interaction network. Deeper insight into PPIs can help us better understand numerous diseases and lead to the development of new diagnostic and therapeutic strategies. PPI interfaces, until recently, were considered undruggable. However, it is now believed that the interfaces contain "hot spots," which could be targeted by small molecules. Such a strategy would require high-quality structural data of PPIs, which are difficult to obtain experimentally. Therefore, in silico modeling can complement or be an alternative to in vitro approaches. There are several computational methods for analyzing the structural data of the binding partners and modeling of the protein-protein dimer/oligomer structure. The major problem with in silico structure prediction of protein assemblies is obtaining sufficient sampling of protein dynamics. One of the methods that can take protein flexibility and the effects of the environment into account is Molecular Dynamics (MD). While sampling of the whole protein-protein association process with plain MD would be computationally expensive, there are several strategies to harness the method to PPI studies while maintaining reasonable use of resources. This chapter reviews known applications of MD in the PPI investigation workflows.

Standard Binding Free-Energy Calculations: How Far Are We from Automation?

Recent success stories suggest that in silico protein-ligand binding free-energy calculations are approaching chemical accuracy. However, their widespread application remains limited by the extensive human intervention required, posing challenges for the neophyte. As such, it is critical to develop automated workflows for estimating protein-ligand binding affinities with minimum personal involvement. Key human efforts include setting up and tuning enhanced-sampling or alchemical-transformation algorithms as a preamble to computational binding free-energy estimations. Additionally, preparing input files, bookkeeping, and postprocessing represent nontrivial tasks. In this Perspective, we discuss recent progress in automating standard binding free-energy calculations, featuring the development of adaptive or parameter-free algorithms, standardization of binding free-energy calculation workflows, and the implementation of user-friendly software. We also assess the current state of automated standard binding free-energy calculations and evaluate the limitations of existing methods. Last, we outline the requirements for future algorithms and workflows to facilitate automated free-energy calculations for diverse protein-ligand complexes.

Computational generation and characterization of IsdA-binding aptamers with single-molecule FRET analysis

Staphylococcus aureus is a major foodborne bacterial pathogen. Early detection of S. aureus is crucial to prevent infections and ensure food quality. The iron-regulated surface determinant protein A (IsdA) of S. aureus is a unique surface protein necessary for sourcing vital iron from host cells for the survival and colonization of the bacteria. The function, structure, and location of the IsdA protein make it an important protein for biosensing applications relating to the pathogen. Here, we report an in-silico approach to develop and validate high-affinity binding aptamers for the IsdA protein detection using custom-designed in-silico tools and single-molecule Fluorescence Resonance Energy Transfer (smFRET) measurements. We utilized in-silico oligonucleotide screening methods and metadynamics-based methods to generate 10 aptamer candidates and characterized them based on the Dissociation Free Energy (DFE) of the IsdA-aptamer complexes. Three of the aptamer candidates were shortlisted for smFRET experimental analysis of binding properties. Limits of detection in the low picomolar range were observed for the aptamers, and the results correlated well with the DFE calculations, indicating the potential of the in-silico approach to support aptamer discovery. This study showcases a computational SELEX method in combination with single-molecule binding studies deciphering effective aptamers against S. aureus IsdA, protein. The established approach demonstrates the ability to expedite aptamer discovery that has the potential to cut costs and predict binding efficacy. The application can be extended to designing aptamers for various protein targets, enhancing molecular recognition, and facilitating the development of high-affinity aptamers for multiple uses.

Assessing Metadynamics and Docking for Absolute Binding Free Energy Calculations Using Severe Acute Respiratory Syndrome Coronavirus 2 Main Protease Inhibitors

Absolute binding free energy (ABFE) calculations can be an important part of the drug discovery process by identifying molecules that have the potential to be strong binders for a biomolecular target. Recent work has used free energy perturbation (FEP) theory for these calculations, focusing on a set of 16 inhibitors of the severe acute respiratory syndrome coronavirus 2 main protease (Mpro). Herein, the same data set is evaluated by metadynamics (MetaD), four different docking programs, and molecular mechanics with generalized Born and surface area solvation. MetaD yields a Kendall τ distance of 0.28 and Pearson r2 of 0.49, which reflect somewhat less accuracy than that from the ABFE FEP results. Notably, it is demonstrated that an ensemble docking protocol by which each ligand is docked into the 13 crystal structures in this data set provides improved performance, particularly when docking is carried out with Glide XP (Kendall τ distance = 0.20, Pearson r2 = 0.71), Glide SP (Kendall τ distance = 0.19, Pearson r2 = 0.66), or AutoDock 4 (Kendall τ distance = 0.21, Pearson r2 = 0.55). The best results are obtained with "superconsensus" docking by averaging the 52 results for each compound using the 4 docking protocols and all 13 crystal structures (Kendall τ distance = 0.18, Pearson r2 = 0.73).

Theoretical investigation of hydroxylated analogues of valinomycin as potassium transporter

Valinomycin is a potent ionophore known for its ability to transport potassium ions across biological membranes. The study focuses on the hydroxylated analogues of valinomycin (HyVLMs) and compares their energy profiles and capabilities for transporting potassium ions across phospholipid membranes. Using metadynamics, we investigated the energy profiles of wildtype valinomycin (VLM_1) and its three hydroxylated analogues (VLM_2, VLM_3, and VLM_4). We observed that all analogues exhibited energy maxima in the centre of the membrane and preferred positions below the phospholipid heads. Furthermore, the entry barriers for membrane penetration were similar among the analogues, suggesting that the hydroxyl group did not significantly affect their passage through the membrane. Transition state calculations provided insights into the ability of valinomycin analogues to capture potassium ions, with VLM_4 showing the lowest activation energy and VLM_2 displaying the highest. Our findings contribute to understanding the mechanisms of potassium transport by valinomycin analogues and highlight their potential as ionophores. The presence of the hydroxyl group is of particular importance because it paves the way for subsequent chemical modifications and the synthesis of new antiviral agents with reduced intrinsic toxicity.

De Novo Design of Anti-COVID Drugs Using Machine Learning-Based Equivariant Diffusion Model Targeting the Spike Protein

The drug discovery and research for an anti-COVID-19 drug has been ongoing despite repurposed drugs in the market. Over time, these drugs were discontinued due to side effects. The search for effective drugs is still under process. The role of Machine Learning (ML) is critical in the search for novel drug compounds. In the current work, using the equivariant diffusion model, we built novel compounds targeting the spike protein of SARS-CoV-2. Using the ML models, 196 de novo compounds were generated which had no hits on any major chemical databases. These novel compounds fulfilled all the criteria of ADMET properties to be lead-like and drug-like compounds. Of the 196 compounds, 15 were docked with high confidence in the target. These compounds were further subjected to molecular docking, the best compound having an IUPAC name of (4aS,4bR,8aS,8bS)-4a,8a-dimethylbiphenylene-1,4,5,8(4aH,4bH,8aH,8bH)-tetraone and a binding score of -6.930 kcal/mol. The principal compound is labeled as CoECG-M1. Density Function Theory (DFT) and Quantum optimization was carried out along with the study of ADMET properties. This suggests that the compound has potential drug-like properties. The docked complex was further subjected to MD simulations, GBSA, and metadynamics simulations to gain insights into the stability of binding. The model can be in the future modified to improve the positive docking rate.

Protein-protein interfaces in molecular glue-induced ternary complexes: classification, characterization, and prediction

Molecular glues are a class of small molecules that stabilize the interactions between proteins. Naturally occurring molecular glues are present in many areas of biology where they serve as central regulators of signaling pathways. Importantly, several clinical compounds act as molecular glue degraders that stabilize interactions between E3 ubiquitin ligases and target proteins, leading to their degradation. Molecular glues hold promise as a new generation of therapeutic agents, including those molecular glue degraders that can redirect the protein degradation machinery in a precise way. However, rational discovery of molecular glues is difficult in part due to the lack of understanding of the protein-protein interactions they stabilize. In this review, we summarize the structures of known molecular glue-induced ternary complexes and the interface properties. Detailed analysis shows different mechanisms of ternary structure formation. Additionally, we also review computational approaches for predicting protein-protein interfaces and highlight the promises and challenges. This information will ultimately help inform future approaches for rational molecular glue discovery.

Mycobacterium Time-Series Genome Analysis Identifies AAC2′ as a Potential Drug Target with Naloxone Showing Potential Bait Drug Synergism

The World Health Organization has put drug resistance in tuberculosis on its list of significant threats, with a critical emphasis on resolving the genetic differences in Mycobacterium tuberculosis. This provides an opportunity for a better understanding of the evolutionary progression leading to anti-microbial resistance. Anti-microbial resistance has a great impact on the economic stability of the global healthcare sector. We performed a timeline genomic analysis from 2003 to 2021 of 578 mycobacterium genomes to understand the pattern underlying genomic variations. Potential drug targets based on functional annotation was subjected to pharmacophore-based screening of FDA-approved phyto-actives. Reaction search, MD simulations, and metadynamics studies were performed. A total of 4,76,063 mutations with a transition/transversion ratio of 0.448 was observed. The top 10 proteins with the least number of mutations were high-confidence drug targets. Aminoglycoside 2′-N-acetyltransferase protein (AAC2′), conferring resistance to aminoglycosides, was shortlisted as a potential drug target based on its function and role in bait drug synergism. Gentamicin-AAC2′ binding pose was used as a pharmacophore template to screen 10,570 phyto-actives. A total of 66 potential hits were docked to obtain naloxone as a lead—active with a docking score of −6.317. Naloxone is an FDA-approved drug that rapidly reverses opioid overdose. This is a classic case of a repurposed phyto-active. Naloxone consists of an amine group, but the addition of the acetyl group is unfavorable, with a reaction energy of 612.248 kcal/mol. With gentamicin as a positive control, molecular dynamic simulation studies were performed for 200 ns to check the stability of binding. Metadynamics-based studies were carried out to compare unbinding energy with gentamicin. The unbinding energies were found to be −68 and −74 kcal/mol for naloxone and gentamycin, respectively. This study identifies naloxone as a potential drug candidate for a bait drug synergistic approach against Mycobacterium tuberculosis.

Drug Discovery by Automated Adaptation of Chemical Structure and Identity

Computer-aided drug design offers the potential to dramatically reduce the cost and effort required for drug discovery. While screening-based methods are valuable in the early stages of hit identification, they are frequently succeeded by iterative, hypothesis-driven computations that require recurrent investment of human time and intuition. To increase automation, we introduce a computational method for lead refinement that combines concerted dynamics of the ligand/protein complex via molecular dynamics simulations with integrated Monte Carlo-based changes in the chemical formula of the ligand. This approach, which we refer to as ligand-exchange Monte Carlo molecular dynamics, accounts for solvent- and entropy-based contributions to competitive binding free energies by coupling the energetics of bound and unbound states during the ligand-exchange attempt. Quantitative comparison of relative binding free energies to reference values from free energy perturbation, conducted in vacuum, indicates that ligand-exchange Monte Carlo molecular dynamics simulations sample relevant conformational ensembles and are capable of identifying strongly binding compounds. Additional simulations demonstrate the use of an implicit solvent model. We speculate that the use of chemical graphs in which exchanges are only permitted between ligands with sufficient similarity may enable an automated search to capture some of the benefits provided by human intuition during hypothesis-guided lead refinement.

Hazardous Shortcuts in Standard Binding Free Energy Calculations

Calculating the standard binding free energies of protein-protein and protein-ligand complexes from atomistic molecular dynamics simulations in explicit solvent is a problem of central importance in computational biophysics. A rigorous strategy for carrying out such calculations is the so-called "geometrical route". In this method, two molecular objects are progressively separated from one another in the presence of orientational and conformational restraints serving to control the change in configurational entropy that accompanies the dissociation process, thereby allowing the computations to converge within simulations of affordable length. Although the geometrical route provides a rigorous theoretical framework, a tantalizing computational shortcut consists of simply leaving out such orientational and conformational degrees of freedom during the separation process. Here the accuracy and convergence of the two approaches are critically compared in the case of two protein-ligand complexes (Abl kinase-SH3:p41 and MDM2-p53:NVP-CGM097) and three protein-protein complexes (pig insulin dimer, SARS-CoV-2 spike RBD:ACE2, and CheA kinase-P2:CheY). The results of the simulations that strictly follow the geometrical route match the experimental standard binding free energies within chemical accuracy. In contrast, simulations bereft of geometrical restraints converge more poorly, yielding inconsistent results that are at variance with the experimental measurements. Furthermore, the orientational and positional time correlation functions of the protein in the unrestrained simulations decay over several microseconds, a time scale that is far longer than the typical simulation times of the geometrical route, which explains why those simulations fail to sample the relevant degrees of freedom during the separation process of the complexes.

Chasing molecular glue degraders: screening approaches

Protein-protein interactions (PPIs) govern all biological processes. Some small molecules modulate PPIs through induced protein proximity. In particular, molecular glue degraders are monovalent compounds that orchestrate interactions between a target protein and an E3 ubiquitin ligase, prompting the proteasomal degradation of the former. This and other pharmacological strategies of targeted protein degradation (e.g. proteolysis-targeting chimeras - PROTACs) overcome some limitations of traditional occupancy-based therapeutics. Here, we provide an overview of the "molecular glue" concept, with a special focus on natural and synthetic inducers of proximity to E3s. We then briefly highlight the serendipitous discoveries of some clinical and preclinical molecular glue degraders, and discuss the first examples of intentional discoveries. Specifically, we outline the different screening strategies reported in this rapidly evolving arena and our thoughts on future perspectives. By mastering the ability to influence PPIs, molecular glue degraders can induce the degradation of unligandable proteins, thus providing an exciting path forward to broaden the targetable proteome.

Challenges and frontiers of computational modelling of biomolecular recognition

Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modelling. Here, we review the challenges and computational approaches developed to characterise biomolecular binding, including molecular docking, molecular dynamics simulations (especially enhanced sampling) and machine learning. Further improvements are still needed in order to accurately and efficiently characterise binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.