Protein Flexibility in Drug Discovery: From Theory to Computation

Molecular Modelling in Bioactive Peptide Discovery and Characterisation

Molecular modelling is a vital tool in the discovery and characterisation of bioactive peptides, providing insights into their structural properties and interactions with biological targets. Many models predicting bioactive peptide function or structure rely on their intrinsic properties, including the influence of amino acid composition, sequence, and chain length, which impact stability, folding, aggregation, and target interaction. Homology modelling predicts peptide structures based on known templates. Peptide-protein interactions can be explored using molecular docking techniques, but there are challenges related to the inherent flexibility of peptides, which can be addressed by more computationally intensive approaches that consider their movement over time, called molecular dynamics (MD). Virtual screening of many peptides, usually against a single target, enables rapid identification of potential bioactive peptides from large libraries, typically using docking approaches. The integration of artificial intelligence (AI) has transformed peptide discovery by leveraging large amounts of data. AlphaFold is a general protein structure prediction tool based on deep learning that has greatly improved the predictions of peptide conformations and interactions, in addition to providing estimates of model accuracy at each residue which greatly guide interpretation. Peptide function and structure prediction are being further enhanced using Protein Language Models (PLMs), which are large deep-learning-derived statistical models that learn computer representations useful to identify fundamental patterns of proteins. Recent methodological developments are discussed in the context of canonical peptides, as well as those with modifications and cyclisations. In designing potential peptide therapeutics, the main outstanding challenge for these methods is the incorporation of diverse non-canonical amino acids and cyclisations.

Toward the Prediction of Binding Events in Very Flexible, Allosteric, Multidomain Proteins

Knowledge of the structures formed by proteins and small molecules is key to understand the molecular principles of chemotherapy and for designing new and more effective drugs. During the early stage of a drug discovery program, it is customary to predict ligand-protein complexes in silico, particularly when screening large compound databases. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology to generate bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites by exploiting only information on the unbound structure and the putative binding sites. The protocol is validated on the paradigm enzyme adenylate kinase, for which we generated a significant fraction of bound-like structures. A fraction of these conformations, employed in ensemble-docking calculations, allowed to find native-like poses of substrates and inhibitors (binding to the active form of the enzyme), as well as catalytically incompetent analogs (binding the inactive form). Our protocol provides a general framework for the generation of bound-like conformations of challenging drug targets that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening, in the prediction of the impact of amino acid mutations on structure and dynamics, and in protein engineering.

Biological Implications of the Intrinsic Deformability of Human Acetylcholinesterase Induced by Diverse Compounds: A Computational Study

The enzyme acetylcholinesterase (AChE) plays a crucial role in the termination of nerve impulses by hydrolyzing the neurotransmitter acetylcholine (ACh). The inhibition of AChE has emerged as a promising therapeutic approach for the management of neurological disorders such as Lewy body dementia and Alzheimer's disease. The potential of various compounds as AChE inhibitors was investigated. In this study, we evaluated the impact of natural compounds of interest on the intrinsic deformability of human AChE using computational biophysical analysis. Our approach incorporates classical dynamics, elastic networks (ENM and NMA), statistical potentials (CUPSAT and SWOTein), energy frustration (Frustratometer), and volumetric cavity analyses (MOLE and PockDrug). The results revealed that cyanidin induced significant changes in the flexibility and rigidity of AChE, especially in the distribution and volume of internal cavities, compared to model inhibitors such as TZ2PA6, and through a distinct biophysical-molecular mechanism from the other inhibitors considered. These findings suggest that cyanidin could offer potential mechanistic pathways for future research and applications in the development of new treatments for neurodegenerative diseases.

Simulation- and AI-directed optimization of 4,6-substituted 1,3,5-triazin-2(1H)-ones as inhibitors of human DNA topoisomerase IIα

The 4,6-substituted-1,3,5-triazin-2(1H)-ones are promising inhibitors of human DNA topoisomerase IIα. To further develop this chemical class targeting the enzyme´s ATP binding site, the triazin-2(1H)-one substitution position 6 was optimized. Inspired by binding of preclinical substituted 9H-purine derivative, bicyclic substituents were incorporated at position 6 and the utility of this modification was validated by a combination of molecular simulations, dynamic pharmacophores, and free energy calculations. Considering also predictions of Deepfrag, a software developed for structure-based lead optimization based on deep learning, compounds with both bicyclic and monocyclic substitutions were synthesized and investigated for their inhibitory activity. The SAR data showed that the bicyclic substituted compounds exhibited good inhibition of topo IIα, comparable to their mono-substituted counterparts. Further evaluation on a panel of human protein kinases showed selectivity for the inhibition of topo IIα. Mechanistic studies indicated that the compounds acted predominantly as catalytic inhibitors, with some exhibiting topo IIα poison effects at higher concentrations. Integration of STD NMR experiments and molecular simulations, provided insights into the binding model and highlighted the importance of the Asn120 interaction and hydrophobic interactions with substituents at positions 4 and 6. In addition, NCI-60 screening demonstrated cytotoxicity of the compounds with bicyclic substituents and identified sensitive human cancer cell lines, underlining the translational relevance of our findings for further preclinical development of this class of compounds. The study highlights the synergy between simulation and AI-based approaches in efficiently guiding molecular design for drug optimization, which has implications for further preclinical development of this class of compounds.

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small molecules is key to understand the molecular principles of chemotherapy and for designing new and more effective drugs. During the early stage of a drug discovery program, it is customary to predict ligand-protein complexes in silico, particularly when screening large compound databases. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology to generate bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites by exploiting only information on the unbound structure and the putative binding sites. The protocol is validated on the paradigm enzyme adenylate kinase, for which we generated a significant fraction of bound-like structures. A fraction of these conformations, employed in ensemble-docking calculations, allowed to find native-like poses of substrates and inhibitors (binding to the active form of the enzyme), as well as catalytically incompetent analogs (binding the inactive form). Our protocol provides a general framework for the generation of bound-like conformations of challenging drug targets that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening, in the prediction of the impact of amino acid mutations on structure and dynamics, and in protein engineering.

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.

Fully Flexible Molecular Alignment Enables Accurate Ligand Structure Modeling

Accurate protein-ligand binding poses are the prerequisites of structure-based binding affinity prediction and provide the structural basis for in-depth lead optimization in small molecule drug design. However, it is challenging to provide reasonable predictions of binding poses for different molecules due to the complexity and diversity of the chemical space of small molecules. Similarity-based molecular alignment techniques can effectively narrow the search range, as structurally similar molecules are likely to have similar binding modes, with higher similarity usually correlated to higher success rates. However, molecular similarity is not consistently high because molecules often require changes to achieve specific purposes, leading to reduced alignment precision. To address this issue, we propose a new alignment method─Z-align. This method uses topological structural information as a criterion for evaluating similarity, reducing the reliance on molecular fingerprint similarity. Our method has achieved success rates significantly higher than those of other methods at moderate levels of similarity. Additionally, our approach can comprehensively and flexibly optimize bond lengths and angles of molecules, maintaining a high accuracy even when dealing with larger molecules. Consequently, our proposed solution helps in achieving more accurate binding poses in protein-ligand docking problems, facilitating the development of small molecule drugs. Z-align is freely available as a web server at https://cloud.zelixir.com/zalign/home.

Virtual high throughput screening of natural peptides against ErbB1 and ErbB2 to identify potential inhibitors for cancer chemotherapy

Human epidermal growth factor receptors (EGFR), namely ErbB1/HER1, ErbB2/HER2/neu, ErbB3/HER3, and ErbB4/HER4, the trans-membrane family of tyrosine kinase receptors, are overexpressed in many types of cancers. These receptors play an important role in cell proliferation, differentiation, invasion, metastasis and angiogenesis including unregulated activation of cancer cells. Overexpression of ErbB1 and ErbB2 that occurs in several types of cancers is associated with poor prognosis leading to resistance to ErbB1-directed therapies. In this connection, promising strategy to overcome the disadvantages of the existing chemotherapeutic drugs is the use of short peptides as anticancer agents. In the present study, we have performed virtual high throughput screening of natural peptides against ErbB1 and ErbB2 to identify potential dual inhibitors and identified five inhibitors based on their binding affinities, ADMET analysis, MD simulation studies and calculation of free energy of binding. These natural peptides could be further exploited for developing drugs for treating cancer.Communicated by Ramaswamy H. Sarma.

Nature-inspired substituted 3-(imidazol-2-yl) morpholines targeting human topoisomerase IIα: Dynophore-derived discovery

The molecular nanomachine, human DNA topoisomerase IIα, plays a crucial role in replication, transcription, and recombination by catalyzing topological changes in the DNA, rendering it an optimal target for cancer chemotherapy. Current clinical topoisomerase II poisons often cause secondary tumors as side effects due to the accumulation of double-strand breaks in the DNA, spurring the development of catalytic inhibitors. Here, we used a dynamic pharmacophore approach to develop catalytic inhibitors targeting the ATP binding site of human DNA topoisomerase IIα. Our screening of a library of nature-inspired compounds led to the discovery of a class of 3-(imidazol-2-yl) morpholines as potent catalytic inhibitors that bind to the ATPase domain. Further experimental and computational studies identified hit compound 17, which exhibited selectivity against the human DNA topoisomerase IIα versus human protein kinases, cytotoxicity against several human cancer cells, and did not induce DNA double-strand breaks, making it distinct from clinical topoisomerase II poisons. This study integrates an innovative natural product-inspired chemistry and successful implementation of a molecular design strategy that incorporates a dynamic component of ligand-target molecular recognition, with comprehensive experimental characterization leading to hit compounds with potential impact on the development of more efficient chemotherapies.

Sequentially Controlled Recognition of Different Proteins Using Programmable Protein Imprinted Nanospheres

Although protein imprinted materials with multiple templates are developed to selectively separate different proteins, it is difficult to achieve the programmed adsorption and separation of different proteins using one material, because the available protein imprinted materials are constructed through irreversible crosslinking and their structures are unprogrammable and non-reconstructive. Herein, a novel nanosphere (MS@PTL-g-PNIPAM) is designed, which not only is temperature and pH responsive but also can dynamically reversibly crosslink/de-crosslink under ultraviolet light of different wavelengths. With the help of the dynamically reversible photo-crosslinking, the nanospheres can be repeatedly programmed into protein imprinted nanospheres toward different target proteins. Moreover, the prepared imprinted nanospheres can easily achieve the controlled rebinding and release of target proteins, benefiting from the introduced temperature- and pH-responsive moieties. As a consequence, this study realizes the specific separation of different target proteins from protein mixture and the real bovine blood sequentially by programming one material. It is resource saving, time saving, recyclable, and it will provide convenience for protein imprinted materials to use in the blood purification, drug delivery, and virus detection.

CSAlign and CSAlign-Dock: Structure alignment of ligands considering full flexibility and application to protein-ligand docking

Structure prediction of protein-ligand complexes, called protein-ligand docking, is a critical computational technique that can be used to understand the underlying principle behind the protein functions at the atomic level and to design new molecules regulating the functions. Protein-ligand docking methods have been employed in structure-based drug discovery for hit discovery and lead optimization. One of the important technical challenges in protein-ligand docking is to account for protein conformational changes induced by ligand binding. A small change such as a single side-chain rotation upon ligand binding can hinder accurate docking. Here we report an increase in docking performance achieved by structure alignment to known complex structures. First, a fully flexible compound-to-compound alignment method CSAlign is developed by global optimization of a shape score. Next, the alignment method is combined with a docking algorithm to dock a new ligand to a target protein when a reference protein-ligand complex structure is available. This alignment-based docking method, called CSAlign-Dock, showed superior performance to ab initio docking methods in cross-docking benchmark tests. Both CSAlign and CSAlign-Dock are freely available as a web server at https://galaxy.seoklab.org/csalign.

Modeling receptor flexibility in the structure-based design of KRASG12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS^G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS^G12C—sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS^G12C. In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS^G12C, we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS^G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

Machine learning assessment of the binding region as a tool for more efficient computational receptor-ligand docking

We present a combined computational approach to protein-ligand binding, which consists of two steps: (1) a deep neural network is used to locate a binding region on a target protein, and (2) molecular docking of a ligand is performed within the specified region to obtain the best pose using Autodock Vina. Our in-house designed neural network was trained using the PepBDB dataset. Although the training dataset consisted of protein-peptide complexes, we show that the approach is not limited to peptides, but also works remarkably well for a large class of non-peptide ligands. The results are compared with those in which the binding region (first step) was provided by Accluster. In cases where no prior experimental data on the binding region are available, our deep neural network provides a fast and effective alternative to classical software for its localization. Our code is available at https://github.com/mksmd/NNforDocking.

Design, Synthesis and Evaluation of Fused Bicyclo[2.2.2]octene as a Potential Core Scaffold for the Non-Covalent Inhibitors of SARS-CoV-2 3CLpro Main Protease

The emergence of SARS-CoV-2, responsible for the global COVID-19 pandemic, requires the rapid development of novel antiviral drugs that would contribute to an effective treatment alongside vaccines. Drug repurposing and development of new molecules targeting numerous viral targets have already led to promising drug candidates. To this end, versatile molecular scaffolds with high functionalization capabilities play a key role. Starting with the clinically used conformationally flexible HIV-1 protease inhibitors that inhibit replication of SARS-CoV-2 and bind major protease 3CL^pro, we designed and synthesized a series of rigid bicyclo[2.2.2]octenes fused to N-substituted succinimides to test whether this core scaffold could support the development of non-covalent 3CL^pro inhibitors. Inhibition assays confirmed that some compounds can inhibit the SARS-CoV-2 main protease; the most promising compound 11a inhibited 3CL^pro in micromolar range (IC₅₀ = 102.2 μM). Molecular simulations of the target-ligand complex in conjunction with dynophore analyses and endpoint free energy calculations provide additional insight and first recommendations for future optimization. The fused bicyclo[2.2.2]octenes can be used as a new potential starting point in the development of non-covalent SARS-CoV-2 3CL^pro protease inhibitors and the study also substantiates the potential of this versatile scaffold for the development of biologically active molecules.

Expanding the Accessible Chemical Space of SIRT2 Inhibitors through Exploration of Binding Pocket Dynamics

Considerations of binding pocket dynamics are one of the crucial aspects of the rational design of binders. Identification of alternative conformational states or cryptic subpockets could lead to the discovery of completely novel groups of the ligands. However, experimental characterization of pocket dynamics, besides being expensive, may not be able to elucidate all of the conformational states relevant for drug discovery projects. In this study, we propose the protocol for computational simulations of sirtuin 2 (SIRT2) binding pocket dynamics and its integration into the structure-based virtual screening (SBVS) pipeline. Initially, unbiased molecular dynamics simulations of SIRT2:inhibitor complexes were performed using optimized force field parameters of SIRT2 inhibitors. Time-lagged independent component analysis (tICA) was used to design pocket-related collective variables (CVs) for enhanced sampling of SIRT2 pocket dynamics. Metadynamics simulations in the tICA eigenvector space revealed alternative conformational states of the SIRT2 binding pocket and the existence of a cryptic subpocket. Newly identified SIRT2 conformational states outperformed experimentally resolved states in retrospective SBVS validation. After performing prospective SBVS, compounds from the under-represented portions of the SIRT2 inhibitor chemical space were selected for in vitro evaluation. Two compounds, NDJ18 and NDJ85, were identified as potent and selective SIRT2 inhibitors, which validated the in silico protocol and opened up the possibility for generalization and broadening of its application. The anticancer effects of the most potent compound NDJ18 were examined on the triple-negative breast cancer cell line. Results indicated that NDJ18 represents a promising structure suitable for further evaluation.

Getting Insights into Structural and Energetic Properties of Reciprocal Peptide-Protein Interactions

Peptide-protein interactions play a key role for many cellular and metabolic processes involved in the onset of largely spread diseases such as cancer and neurodegenerative pathologies. Despite the progress in the structural characterization of peptide-protein interfaces, the in-depth knowledge of the molecular details behind their interactions is still a daunting task. Here, we present the first comprehensive in silico morphological and energetic study of peptide binding sites by focusing on both peptide and protein standpoints. Starting from the PixelDB database, a nonredundant benchmark collection of high-quality 3D crystallographic structures of peptide-protein complexes, a classification analysis of the most representative categories based on the nature of each cocrystallized peptide has been carried out. Several interpretable geometrical and energetic descriptors have been computed both from peptide and target protein sides in the attempt to unveil physicochemical and structural causative correlations. Finally, we investigated the most frequent peptide-protein residue pairs at the binding interface and made extensive energetic analyses, based on GRID MIFs, with the aim to study the peptide affinity-enhancing interactions to be further exploited in rational drug design strategies.

Dynophore-Based Approach in Virtual Screening: A Case of Human DNA Topoisomerase IIα

In this study, we utilized human DNA topoisomerase IIα as a model target to outline a dynophore-based approach to catalytic inhibitor design. Based on MD simulations of a known catalytic inhibitor and the native ATP ligand analog, AMP-PNP, we derived a joint dynophore model that supplements the static structure-based-pharmacophore information with a dynamic component. Subsequently, derived pharmacophore models were employed in a virtual screening campaign of a library of natural compounds. Experimental evaluation identified flavonoid compounds with promising topoisomerase IIα catalytic inhibition and binding studies confirmed interaction with the ATPase domain. We constructed a binding model through docking and extensively investigated it with molecular dynamics MD simulations, essential dynamics, and MM-GBSA free energy calculations, thus reconnecting the new results to the initial dynophore-based screening model. We not only demonstrate a new design strategy that incorporates a dynamic component of molecular recognition, but also highlight new derivates in the established flavonoid class of topoisomerase II inhibitors.

Analytical and functional aspects of protein-ligand interactions: Beyond induced fit and conformational selection

Ligand-dependent changes in protein conformation are foundational to biology. Historical mechanistic models for substrate-specific proteins are induced fit (IF) and conformational selection (CS), which invoke a change in protein conformation after ligand binds or before ligand binds, respectively. These mechanisms have important, but rarely discussed, functional relevance because IF vs. CS can differentially affect a protein's substrate specificity or promiscuity, and its regulatory properties. The modern view of proteins as conformational ensembles in both ligand free and bound states, together with the realization that most proteins exhibit some substrate promiscuity, demands a deeper interpretation of the historical models and provides an opportunity to improve mechanistic analyses. Here we describe alternative analytical strategies for distinguishing the historical models, including the more complex expanded versions of IF and CS. Functional implications of the different models are described. We provide an alternative perspective based on protein ensembles interacting with ligand ensembles that clarifies how a single protein can 'apparently' exploit different mechanisms for different ligands. Mechanistic information about protein ensembles can be optimized when they are probed with multiple ligands.

Decrypting a Cryptic Allosteric Pocket in H. pylori Glutamate Racemase

One of our greatest challenges in drug design is targeting cryptic allosteric pockets in enzyme targets. Drug leads that do bind to these cryptic pockets are often discovered during HTS campaigns, and the mechanisms of action are rarely understood. Nevertheless, it is often the case that the allosteric pocket provides the best option for drug development against a given target. In the current studies we present a successful way forward in rationally exploiting the cryptic allosteric pocket of H. pylori glutamate racemase, an essential enzyme in this pathogen's life cycle. A wide range of computational and experimental methods are employed in a workflow leading to the discovery of a series of natural product allosteric inhibitors which occupy the allosteric pocket of this essential racemase. The confluence of these studies reveals a fascinating source of the allosteric inhibition, which centers on the abolition of essential monomer-monomer coupled motion networks.

Efficiency of Homology Modeling Assisted Molecular Docking in G-protein Coupled Receptors

BACKGROUND

Molecular docking is in regular practice to assess ligand affinity on a target protein crystal structure. In the absence of protein crystal structure, the homology modeling or comparative modeling is the best alternative to elucidate the relationship details between a ligand and protein at the molecular level. The development of accurate homology modeling (HM) and its integration with molecular docking (MD) is essential for successful, rational drug discovery.

OBJECTIVE

The G-protein coupled receptors (GPCRs) are attractive therapeutic targets due to their immense role in human pharmacology. The GPCRs are membrane-bound proteins with the complex constitution, and the understanding of their activation and inactivation mechanisms is quite challenging. Over the past decade, there has been a rapid expansion in the number of solved G-protein-coupled receptor (GPCR) crystal structures; however, the majority of the GPCR structures remain unsolved. In this context, HM guided MD has been widely used for structure-based drug design (SBDD) of GPCRs.

METHODS

The focus of this review is on the recent (i) developments on HM supported GPCR drug discovery in the absence of GPCR crystal structures and (ii) application of HM in understanding the ligand interactions at the binding site, virtual screening, determining receptor subtype selectivity and receptor behaviour in comparison with GPCR crystal structures.

RESULTS

The HM in GPCRs has been extremely challenging due to the scarcity in template structures. In such a scenario, it is difficult to get accurate HM that can facilitate understanding of the ligand-receptor interactions. This problem has been alleviated to some extent by developing refined HM based on incorporating active /inactive ligand information and inducing protein flexibility. In some cases, HM proteins were found to outscore crystal structures.

CONCLUSION

The developments in HM have been highly operative to gain insights about the ligand interaction at the binding site and receptor functioning at the molecular level. Thus, HM guided molecular docking may be useful for rational drug discovery for the GPCRs mediated diseases.

Metadynamics-Based Approaches for Modeling the Hypoxia-Inducible Factor 2α Ligand Binding Process

Several methods based on enhanced-sampling molecular dynamics have been proposed for studying ligand binding processes. Here, we developed a protocol that combines the advantages of steered molecular dynamics (SMD) and metadynamics. While SMD is proposed for investigating possible unbinding pathways of the ligand and identifying the preferred one, metadynamics, with the path collective variable (PCV) formalism, is suggested to explore the binding processes along the pathway defined on the basis of SMD, by using only two CVs. We applied our approach to the study of binding of two known ligands to the hypoxia-inducible factor 2α, where the buried binding cavity makes simulation of the process a challenging task. Our approach allowed identification of the preferred entrance pathway for each ligand, highlighted the features of the bound and intermediate states in the free-energy surface, and provided a binding affinity scale in agreement with experimental data. Therefore, it seems to be a suitable tool for elucidating ligand binding processes of similar complex systems.

Current and Future Challenges in Modern Drug Discovery

Drug discovery is an expensive, time-consuming, and risky business. To avoid late-stage failure, learnings from past projects and the development of new approaches are crucial. New modalities and emerging new target spaces allow the exploration of unprecedented indications or to address so far undrugable targets. Late-stage attrition is usually attributed to the lack of efficacy or to compound-related safety issues. Efficacy has been shown to be related to a strong genetic link to human disease, a better understanding of the target biology, and the availability of biomarkers to bridge from animals to humans. Compound safety can be improved by ligand optimization, which is becoming increasingly demanding for difficult targets. Therefore, new strategies include the design of allosteric ligands, covalent binders, and other modalities. Design methods currently heavily rely on artificial intelligence and advanced computational methods such as free energy calculations and quantum chemistry. Especially for quantum chemical methods, a more detailed overview is given in this chapter.

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a "sample and select" method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls (¹H, ¹⁵N, and ¹³C'). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.

Molecular Docking: Challenges, Advances and its Use in Drug Discovery Perspective

Molecular docking is a process through which small molecules are docked into the macromolecular structures for scoring its complementary values at the binding sites. It is a vibrant research area with dynamic utility in structure-based drug-designing, lead optimization, biochemical pathway and for drug designing being the most attractive tools. Two pillars for a successful docking experiment are correct pose and affinity prediction. Each program has its own advantages and drawbacks with respect to their docking accuracy, ranking accuracy and time consumption so a general conclusion cannot be drawn. Moreover, users don't always consider sufficient diversity in their test sets which results in certain programs to outperform others. In this review, the prime focus has been laid on the challenges of docking and troubleshooters in existing programs, underlying algorithmic background of docking, preferences regarding the use of docking programs for best results illustrated with examples, comparison of performance for existing tools and algorithms, state of art in docking, recent trends of diseases and current drug industries, evidence from clinical trials and post-marketing surveillance are discussed. These aspects of the molecular drug designing paradigm are quite controversial and challenging and this review would be an asset to the bioinformatics and drug designing communities.

Combined strategies in structure-based virtual screening

The identification and optimization of lead compounds are inalienable components in drug design and discovery pipelines. As a powerful computational approach for the identification of hits with novel structural scaffolds, structure-based virtual screening (SBVS) has exhibited a remarkably increasing influence in the early stages of drug discovery. During the past decade, a variety of techniques and algorithms have been proposed and tested with different purposes in the scope of SBVS. Although SBVS has been a common and proven technology, it still shows some challenges and problems that are needed to be addressed, where the negative influence regardless of protein flexibility and the inaccurate prediction of binding affinity are the two major challenges. Here, focusing on these difficulties, we summarize a series of combined strategies or workflows developed by our group and others. Furthermore, several representative successful applications from recent publications are also discussed to demonstrate the effectiveness of the combined SBVS strategies in drug discovery campaigns.

A Comprehensive Review on Current Advances in Peptide Drug Development and Design

Protein-protein interactions (PPIs) execute many fundamental cellular functions and have served as prime drug targets over the last two decades. Interfering intracellular PPIs with small molecules has been extremely difficult for larger or flat binding sites, as antibodies cannot cross the cell membrane to reach such target sites. In recent years, peptides smaller size and balance of conformational rigidity and flexibility have made them promising candidates for targeting challenging binding interfaces with satisfactory binding affinity and specificity. Deciphering and characterizing peptide-protein recognition mechanisms is thus central for the invention of peptide-based strategies to interfere with endogenous protein interactions, or improvement of the binding affinity and specificity of existing approaches. Importantly, a variety of computation-aided rational designs for peptide therapeutics have been developed, which aim to deliver comprehensive docking for peptide-protein interaction interfaces. Over 60 peptides have been approved and administrated globally in clinics. Despite this, advances in various docking models are only on the merge of making their contribution to peptide drug development. In this review, we provide (i) a holistic overview of peptide drug development and the fundamental technologies utilized to date, and (ii) an updated review on key developments of computational modeling of peptide-protein interactions (PepPIs) with an aim to assist experimental biologists exploit suitable docking methods to advance peptide interfering strategies against PPIs.

Probing the local conformational flexibility in receptor recognition: mechanistic insight from an atomic-scale investigation

Inherent protein conformational flexibility is important for biomolecular recognition, but this critical property is often neglected in several studies. This event can lead to large deviations in the research results. In the current contribution, we disclose the effects of the local conformational flexibility on receptor recognition by using an atomic-scale computational method. The results indicated that both static and dynamic reaction modes have noticeable differences, and these originated from the structural features of the protein molecules. Dynamic interaction results displayed that the structural stability and conformational flexibility of the proteins had a significant influence on the recognition processes. This point related closely to the characteristics of the flexible loop regions where bixin located within the protein structures. The energy decomposition analyses and circular dichroism results validated the rationality of the recognition studies. More importantly, the conformational and energy changes of some residues around the bixin binding domain were found to be vital to biological reactions. These microscopic findings clarified the nature of the phenomenon that the local conformational flexibility could intervene in receptor recognition. Obviously, this report may provide biophysical evidence for the exploration of the structure-function relationships of the biological receptors in the human body.

Holo-like and Druggable Protein Conformations from Enhanced Sampling of Binding Pocket Volume and Shape

Understanding molecular recognition of small molecules by proteins in atomistic detail is key for drug design. Molecular docking is a widely used computational method to mimic ligand-protein association in silico. However, predicting conformational changes occurring in proteins upon ligand binding is still a major challenge. Ensemble docking approaches address this issue by considering a set of different conformations of the protein obtained either experimentally or from computer simulations, e.g., molecular dynamics. However, holo structures prone to host (the correct) ligands are generally poorly sampled by standard molecular dynamics simulations of the apo protein. In order to address this limitation, we introduce a computational approach based on metadynamics simulations called ensemble docking with enhanced sampling of pocket shape (EDES) that allows holo-like conformations of proteins to be generated by exploiting only their apo structures. This is achieved by defining a set of collective variables that effectively sample different shapes of the binding site, ultimately mimicking the steric effect due to the ligand. We assessed the method on three challenging proteins undergoing different extents of conformational changes upon ligand binding. In all cases our protocol generates a significant fraction of structures featuring a low RMSD from the experimental holo geometry. Moreover, ensemble docking calculations using those conformations yielded in all cases native-like poses among the top-ranked ones.

In Silico Target Prediction for Small Molecules

Drugs modulate disease states through their actions on targets in the body. Determining these targets aids the focused development of new treatments, and helps to better characterize those already employed. One means of accomplishing this is through the deployment of in silico methodologies, harnessing computational analytical and predictive power to produce educated hypotheses for experimental verification. Here, we provide an overview of the current state of the art, describe some of the well-established methods in detail, and reflect on how they, and emerging technologies promoting the incorporation of complex and heterogeneous data-sets, can be employed to improve our understanding of (poly)pharmacology.

Empirical Scoring Functions for Structure-Based Virtual Screening: Applications, Critical Aspects, and Challenges

Structure-based virtual screening (VS) is a widely used approach that employs the knowledge of the three-dimensional structure of the target of interest in the design of new lead compounds from large-scale molecular docking experiments. Through the prediction of the binding mode and affinity of a small molecule within the binding site of the target of interest, it is possible to understand important properties related to the binding process. Empirical scoring functions are widely used for pose and affinity prediction. Although pose prediction is performed with satisfactory accuracy, the correct prediction of binding affinity is still a challenging task and crucial for the success of structure-based VS experiments. There are several efforts in distinct fronts to develop even more sophisticated and accurate models for filtering and ranking large libraries of compounds. This paper will cover some recent successful applications and methodological advances, including strategies to explore the ligand entropy and solvent effects, training with sophisticated machine-learning techniques, and the use of quantum mechanics. Particular emphasis will be given to the discussion of critical aspects and further directions for the development of more accurate empirical scoring functions.

Selection of protein conformations for structure-based polypharmacology studies

Several drugs exert their therapeutic effect through the modulation of multiple targets. Structure-based approaches hold great promise for identifying compounds with the desired polypharmacological profiles. These methods use knowledge of the protein binding sites to identify stereoelectronically complementary ligands. The selection of the most suitable protein conformations to be used in the design process is vital, especially for multitarget drug design in which the same ligand has to be accommodated in multiple binding pockets. Herein, we focus on currently available techniques for the selection of the most suitable protein conformations for multitarget drug design, compare the potential advantages and limitations of each method, and comment on how their combination could help in polypharmacology drug design.

A selective method for optimizing ensemble docking-based experiments on an InhA Fully-Flexible receptor model

BACKGROUND

In the rational drug design process, an ensemble of conformations obtained from a molecular dynamics simulation plays a crucial role in docking experiments. Some studies have found that Fully-Flexible Receptor (FFR) models predict realistic binding energy accurately and improve scoring to enhance selectiveness. At the same time, methods have been proposed to reduce the high computational costs involved in considering the explicit flexibility of proteins in receptor-ligand docking. This study introduces a novel method to optimize ensemble docking-based experiments by reducing the size of an InhA FFR model at docking runtime and scaling docking workflow invocations on cloud virtual machines.

RESULTS

First, in order to find the most affordable cost-benefit pool of virtual machines, we evaluated the performance of the docking workflow invocations in different configurations of Azure instances. Second, we validated the gains obtained by the proposed method based on the quality of the Reduced Fully-Flexible Receptor (RFFR) models produced using AutoDock4.2. The analyses show that the proposed method reduced the model size by approximately 50% while covering at least 86% of the best docking results from the 74 ligands tested. Third, we tested our novel method using AutoDock Vina, a different docking software, and showed the positive accuracy achieved in the resulting RFFR models. Finally, our results demonstrated that the method proposed optimized ensemble docking experiments and is applicable to different docking software. In addition, it detected new binding modes, which would be unreachable if employing only the rigid structure used to generate the InhA FFR model.

CONCLUSIONS

Our results showed that the selective method is a valuable strategy for optimizing ensemble docking-based experiments using different docking software. The RFFR models produced by discarding non-promising snapshots from the original model are accurately shaped for a larger number of ligands, and the elapsed time spent in the ensemble docking experiments are considerably reduced.

Protein-peptide docking: opportunities and challenges

Peptides have recently attracted much attention as promising drug candidates. Rational design of peptide-derived therapeutics usually requires structural characterization of the underlying protein-peptide interaction. Given that experimental characterization can be difficult, reliable computational tools are needed. In recent years, a variety of approaches have been developed for 'protein-peptide docking', that is, predicting the structure of the protein-peptide complex, starting from the protein structure and the peptide sequence, including variable degrees of information about the peptide binding site and/or conformation. In this review, we provide an overview of protein-peptide docking methods and outline their capabilities, limitations, and applications in structure-based drug design. Key challenges are also briefly discussed, such as modeling of large-scale conformational changes upon binding, scoring of predicted models, and optimal inclusion of varied types of experimental data and theoretical predictions into an integrative modeling process.

Fully Flexible Docking via Reaction-Coordinate-Independent Molecular Dynamics Simulations

Predicting the geometry of protein-ligand binding complexes is of primary importance for structure-based drug discovery. Molecular dynamics (MD) is emerging as a reliable computational tool for use in conjunction with, or an alternative to, docking methods. However, simulating the protein-ligand binding process often requires very expensive simulations. This drastically limits the practical application of MD-based approaches. Here, we propose a general framework to accelerate the generation of putative protein-ligand binding modes using potential-scaled MD simulations. The proposed dynamical protocol has been applied to two pharmaceutically relevant systems (GSK-3β and the N-terminal domain of HSP90α). Our approach is fully independent of any predefined reaction coordinate (or collective variable). It identified the correct binding mode of several ligands and can thus save valuable computational time in dynamic docking simulations.

Approaching Pharmacological Space: Events and Components

With a view to introducing the concept of pharmacological space and its potential applications in investigating and predicting the toxic mechanisms of xenobiotics, this opening chapter describes the logical relations between conformational behavior, physicochemical properties and binding spaces, which are seen as the three key elements composing the pharmacological space. While the concept of conformational space is routinely used to encode molecular flexibility, the concepts of property spaces and, particularly, of binding spaces are more innovative. Indeed, their descriptors can find fruitful applications (a) in describing the dynamic adaptability a given ligand experiences when inserted into a specific environment, and (b) in parameterizing the flexibility a ligand retains when bound to a biological target. Overall, these descriptors can conveniently account for the often disregarded entropic factors and as such they prove successful when inserted in ligand- or structure-based predictive models. Notably, and although binding space parameters can clearly be derived from MD simulations, the chapter will illustrate how docking calculations, despite their static nature, are able to evaluate ligand's flexibility by analyzing several poses for each ligand. Such an approach, which represents the founding core of the binding space concept, can find various applications in which the related descriptors show an impressive enhancing effect on the statistical performances of the resulting predictive models.

Designing interactions by control of protein-ligand complex conformation: tuning arginine-arene interaction geometry for enhanced electrostatic protein-ligand interactions

3-Benzamido-2-O-sulfo-galactosides can be designed to control protein conformation into forming entropically favourable galectin-3-arginine salt bridges with ligand sulfates.

We investigated galectin-3 binding to 3-benzamido-2-O-sulfo-galactoside and -thiodigalactoside ligands using a combination of site-specific mutagenesis, X-ray crystallography, computational approaches, and binding thermodynamics measurements. The results reveal a conformational variability in a surface-exposed arginine (R144) side chain in response to different aromatic C3-substituents of bound galactoside-based ligands. Fluorinated C3-benzamido substituents induced a shift in the side-chain conformation of R144 to allow for an entropically favored electrostatic interaction between its guanidine group and the 2-O-sulfate of the ligand. By contrast, binding of ligands with non-fluorinated substituents did not trigger a conformational change of R144. Hence, a sulfate–arginine electrostatic interaction can be tuned by the choice of ligand C3-benzamido structures to favor specific interaction modes and geometries. These results have important general implications for ligand design, as the proper choice of arginine–aromatic interacting partners opens up for ligand-controlled protein conformation that in turn may be systematically exploited in ligand design.

Dynamic Docking: A Paradigm Shift in Computational Drug Discovery

Molecular docking is the methodology of choice for studying in silico protein-ligand binding and for prioritizing compounds to discover new lead candidates. Traditional docking simulations suffer from major limitations, mostly related to the static or semi-flexible treatment of ligands and targets. They also neglect solvation and entropic effects, which strongly limits their predictive power. During the last decade, methods based on full atomistic molecular dynamics (MD) have emerged as a valid alternative for simulating macromolecular complexes. In principle, compared to traditional docking, MD allows the full exploration of drug-target recognition and binding from both the mechanistic and energetic points of view (dynamic docking). Binding and unbinding kinetic constants can also be determined. While dynamic docking is still too computationally expensive to be routinely used in fast-paced drug discovery programs, the advent of faster computing architectures and advanced simulation methodologies are changing this scenario. It is feasible that dynamic docking will replace static docking approaches in the near future, leading to a major paradigm shift in in silico drug discovery. Against this background, we review the key achievements that have paved the way for this progress.

Tailored Approaches in Drug Development and Diagnostics: From Molecular Design to Biological Model Systems

Approaches to increase the efficiency in developing drugs and diagnostics tools, including new drug delivery and diagnostic technologies, are needed for improved diagnosis and treatment of major diseases and health problems such as cancer, inflammatory diseases, chronic wounds, and antibiotic resistance. Development within several areas of research ranging from computational sciences, material sciences, bioengineering to biomedical sciences and bioimaging is needed to realize innovative drug development and diagnostic (DDD) approaches. Here, an overview of recent progresses within key areas that can provide customizable solutions to improve processes and the approaches taken within DDD is provided. Due to the broadness of the area, unfortunately all relevant aspects such as pharmacokinetics of bioactive molecules and delivery systems cannot be covered. Tailored approaches within (i) bioinformatics and computer-aided drug design, (ii) nanotechnology, (iii) novel materials and technologies for drug delivery and diagnostic systems, and (iv) disease models to predict safety and efficacy of medicines under development are focused on. Current developments and challenges ahead are discussed. The broad scope reflects the multidisciplinary nature of the field of DDD and aims to highlight the convergence of biological, pharmaceutical, and medical disciplines needed to meet the societal challenges of the 21st century.

Evaluation of the novel algorithm of flexible ligand docking with moveable target-protein atoms

We present the novel docking algorithm based on the Tensor Train decomposition and the TT-Cross global optimization. The algorithm is applied to the docking problem with flexible ligand and moveable protein atoms. The energy of the protein-ligand complex is calculated in the frame of the MMFF94 force field in vacuum. The grid of precalculated energy potentials of probe ligand atoms in the field of the target protein atoms is not used. The energy of the protein-ligand complex for any given configuration is computed directly with the MMFF94 force field without any fitting parameters. The conformation space of the system coordinates is formed by translations and rotations of the ligand as a whole, by the ligand torsions and also by Cartesian coordinates of the selected target protein atoms. Mobility of protein and ligand atoms is taken into account in the docking process simultaneously and equally. The algorithm is realized in the novel parallel docking SOL-P program and results of its performance for a set of 30 protein-ligand complexes are presented. Dependence of the docking positioning accuracy is investigated as a function of parameters of the docking algorithm and the number of protein moveable atoms. It is shown that mobility of the protein atoms improves docking positioning accuracy. The SOL-P program is able to perform docking of a flexible ligand into the active site of the target protein with several dozens of protein moveable atoms: the native crystallized ligand pose is correctly found as the global energy minimum in the search space with 157 dimensions using 4700 CPU ∗ h at the Lomonosov supercomputer.

Protein-ligand (un)binding kinetics as a new paradigm for drug discovery at the crossroad between experiments and modelling

In the last three decades, protein and nucleic acid structure determination and comprehension of the mechanisms, leading to their physiological and pathological functions, have become a cornerstone of biomedical sciences. A deep understanding of the principles governing the fates of cells and tissue at the molecular level has been gained over the years, offering a solid basis for the rational design of drugs aimed at the pharmacological treatment of numerous diseases. Historically, affinity indicators (i.e. Kd and IC50/EC50) have been assumed to be valid indicators of the in vivo efficacy of a drug. However, recent studies pointed out that the kinetics of the drug-receptor binding process could be as important or even more important than affinity in determining the drug efficacy. This eventually led to a growing interest in the characterisation and prediction of the rate constants of protein-ligand association and dissociation. For instance, a drug with a longer residence time can kinetically select a given receptor over another, even if the affinity for both receptors is comparable, thus increasing its therapeutic index. Therefore, understanding the molecular features underlying binding and unbinding processes is of central interest towards the rational control of drug binding kinetics. In this review, we report the theoretical framework behind protein-ligand association and highlight the latest advances in the experimental and computational approaches exploited to investigate the binding kinetics.

Enantioselective recognition of an isomeric ligand by a biomolecule: mechanistic insights into static and dynamic enantiomeric behavior and structural flexibility

Enantioselective biorecognition is a vital trigger that results in remarkable enantiomeric differences in the biochemical behavior of chiral substances.