Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design

Exploring an Intracellular Allosteric Site of CC-Chemokine Receptor 4 from 3D Models, Probe Simulations, and Mutagenesis

We applied our previously developed probe confined dynamic mapping protocol, which combines enhanced sampling molecular dynamics (MD) simulations and fragment-based approaches, to identify the binding site of GSK2239633A (N-[[3-[[3-[(5-chlorothiophen-2-yl)sulfonylamino]-4-methoxyindazol-1-yl]methyl]phenyl]methyl]-2-hydroxy-2-methylpropanamide), a selective CC-chemokine receptor type 4 (CCR4) negative allosteric modulator, using CCR4 homology and AlphaFold models. By comparing the performance across five computational models, we identified conserved (K3108.49 and Y3047.53) and non-conserved (M2436.36) residue hotspots for GSK2239633A binding, which were validated by mutagenesis and bioluminescence resonance energy transfer assay. Further analysis of 3D models and MD simulations highlighted the pair of residues 6.36 and 7.56 that might account for antagonist selectivity among chemokine receptors. Our in silico protocol provides a promising approach for characterizing ligand binding sites in membrane proteins, considering receptor dynamics and adaptability and guiding protein template selection for ligand design.

Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms

In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.

Identification of 3, 4-dihydroxy complexes as potential antiviral via DFT, molecular docking, molecular dynamics and MM/PBSA against rabies and dengue receptors

The quest to identify antiviral drug candidates for dengue and rabies viral diseases is a great challenge for the researchers. While different research is being conducted on the repurposed drugs against these two viruses, no drug compound has gained success in treating them. Therefore, in this study, 3, 4-dihydroxy complexes have been virtually designed to investigate their antiviral properties and analyze their efficiency in interaction with the concerned viral diseases. DFT calculations are carried out to study the electronic and thermodynamic properties to understand the stability and reactivity of the reported compounds. These compounds were subjected to molecular docking studies to understand the binding interactions with NS5 Dengue virus mRNA 2'-O-methyltransferase and phosphoprotein C-terminal domain of Rabies virus. MD simulation, hydrogen bond analysis, and MM/PBSA were performed at 100 ns to support the obtained docking results.Communicated by Ramaswamy H. Sarma.

Expanded ensemble predictions of absolute binding free energies in the SAMPL9 host-guest challenge

As part of the SAMPL9 community-wide blind host-guest challenge, we implemented an expanded ensemble workflow to predict absolute binding free energies for 13 small molecules against pillar[6]arene. Notable features of our protocol include consideration of a variety of protonation and enantiomeric states for both host and guests, optimization of alchemical intermediates, and analysis of free energy estimates and their uncertainty using large numbers of simulation replicates performed using distributed computing. Our predictions of absolute binding free energies resulted in a mean absolute error of 2.29 kcal mol-1 and an R2 of 0.54. Overall, results show that expanded ensemble calculations using all-atom molecular dynamics simulations are a valuable and efficient computational tool in predicting absolute binding free energies.

Integrated Covalent Drug Design Workflow Using Site Identification by Ligand Competitive Saturation

Covalent drug design is an important component in drug discovery. Traditional drugs interact with their target in a reversible equilibrium, while irreversible covalent drugs increase the drug-target interaction duration by forming a covalent bond with targeted residues and thus may offer a more effective therapeutic approach. To facilitate the design of this class of ligands, computational methods can be used to help identify reactive nucleophilic residues, frequently cysteines, on a target protein for covalent binding, to test various warhead groups for their potential reactivities, and to predict noncovalent contributions to binding that can facilitate drug-target interactions that are important for binding specificity. To further aid covalent drug design, we extended a functional group mapping approach based on explicit solvent all-atom molecular simulations (SILCS: site identification by ligand competitive saturation) that intrinsically considers protein flexibility, functional group, and protein desolvation along with functional group-protein interactions. Through docking of a library of representative warhead fragments using SILCS-Monte Carlo (SILCS-MC), reactive cysteines can be correctly identified for proteins being tested. Furthermore, a machine learning model was trained to quantify the effectiveness of various warhead groups for proteins using metrics from SILCS-MC as well as experimental model compound warhead reactivity data. The ability to rank covalent molecular binders with similar warheads using SILCS ligand grid free energy (LGFE) ranking was also tested for several proteins. Based on these tools, an integrated SILCS-based workflow was developed, named SILCS-Covalent, which can both qualitatively and quantitatively inform covalent drug discovery.

Computer-Aided Drug Design: An Update

Computer-aided drug design (CADD) approaches are playing an increasingly important role in understanding the fundamentals of ligand-receptor interactions and helping medicinal chemists design therapeutics. About 5 years ago, we presented a chapter devoted to an overview of CADD methods and covered typical CADD protocols including structure-based drug design (SBDD) and ligand-based drug design (LBDD) approaches that were frequently used in the antibiotic drug design process. Advances in computational hardware and algorithms and emerging CADD methods are enhancing the accuracy and ability of CADD in drug design and development. In this chapter, an update to our previous chapter is provided with a focus on new CADD approaches from our laboratory and other peers that can be employed to facilitate the development of antibiotic therapeutics.

In silico identification of a β2-adrenoceptor allosteric site that selectively augments canonical β2AR-Gs signaling and function

Activation of β₂-adrenoceptors (β₂ARs) causes airway smooth muscle (ASM) relaxation and bronchodilation, and β₂AR agonists (β-agonists) are front-line treatments for asthma and other obstructive lung diseases. However, the therapeutic efficacy of β-agonists is limited by agonist-induced β₂AR desensitization and noncanonical β₂AR signaling involving β-arrestin that is shown to promote asthma pathophysiology. Accordingly, we undertook the identification of an allosteric site on β₂AR that could modulate the activity of β-agonists to overcome these limitations. We employed the site identification by ligand competitive saturation (SILCS) computational method to comprehensively map the entire 3D structure of in silico-generated β₂AR intermediate conformations and identified a putative allosteric binding site. Subsequent database screening using SILCS identified drug-like molecules with the potential to bind to the site. Experimental assays in HEK293 cells (expressing recombinant wild-type human β₂AR) and human ASM cells (expressing endogenous β₂AR) identified positive and negative allosteric modulators (PAMs and NAMs) of β₂AR as assessed by regulation of β-agonist-stimulation of cyclic AMP generation. PAMs/NAMs had no effect on β-agonist-induced recruitment of β-arrestin to β₂AR- or β-agonist-induced loss of cell surface expression in HEK293 cells expressing β₂AR. Mutagenesis analysis of β₂AR confirmed the SILCS identified site based on mutants of amino acids R131, Y219, and F282. Finally, functional studies revealed augmentation of β-agonist-induced relaxation of contracted human ASM cells and bronchodilation of contracted airways. These findings identify a allosteric binding site on the β₂AR, whose activation selectively augments β-agonist-induced Gs signaling, and increases relaxation of ASM cells, the principal therapeutic effect of β-agonists.

SILCS-RNA: Toward a Structure-Based Drug Design Approach for Targeting RNAs with Small Molecules

RNA molecules can act as potential drug targets in different diseases, as their dysregulated expression or misfolding can alter various cellular processes. Noncoding RNAs account for ∼70% of the human genome, and these molecules can have complex tertiary structures that present a great opportunity for targeting by small molecules. In the present study, the site identification by ligand competitive saturation (SILCS) computational approach is extended to target RNA, termed SILCS-RNA. Extensions to the method include an enhanced oscillating excess chemical potential protocol for the grand canonical Monte Carlo calculations and individual simulations of the neutral and charged solutes from which the SILCS functional group affinity maps (FragMaps) are calculated for subsequent binding site identification and docking calculations. The method is developed and evaluated against seven RNA targets and their reported small molecule ligands. SILCS-RNA provides a detailed characterization of the functional group affinity pattern in the small molecule binding sites, recapitulating the types of functional groups present in the ligands. The developed method is also shown to be useful for identification of new potential binding sites and identifying ligand moieties that contribute to binding, granular information that can facilitate ligand design. However, limitations in the method are evident including the ability to map the regions of binding sites occupied by ligand phosphate moieties and to fully account for the wide range of conformational heterogeneity in RNA associated with binding of different small molecules, emphasizing inherent challenges associated with applying computer-aided drug design methods to RNA. While limitations are present, the current study indicates how the SILCS-RNA approach may enhance drug discovery efforts targeting RNAs with small molecules.

hERG Blockade Prediction by Combining Site Identification by Ligand Competitive Saturation and Physicochemical Properties

Human ether-a-go-go-related gene (hERG) potassium channel is well-known contributor to drug-induced cardiotoxicity and therefore an extremely important target when performing safety assessments of drug candidates. Ligand-based approaches in connection with quantitative structure active relationships (QSAR) analyses have been developed to predict hERG toxicity. Availability of the recent published cryogenic electron microscopy (cryo-EM) structure for the hERG channel opened the prospect for using structure-based simulation and docking approaches for hERG drug liability predictions. In recent time, the idea of combining structure- and ligand-based approaches for modeling hERG drug liability has gained momentum offering improvements in predictability when compared to ligand-based QSAR practices alone. The present article demonstrates uniting the structure-based SILCS (site-identification by ligand competitive saturation) approach in conjunction with physicochemical properties to develop predictive models for hERG blockade. This combination leads to improved model predictability based on Pearson's R and percent correct (represents rank-ordering of ligands) metric for different validation sets of hERG blockers involving diverse chemical scaffold and wide range of pIC50 values. The inclusion of the SILCS structure-based approach allows determination of the hERG region to which compounds bind and the contribution of different chemical moieties in the compounds to blockade, thereby facilitating the rational ligand design to minimize hERG liability.