From G Protein-coupled Receptor Structure Resolution to Rational Drug Design

Integrating biocomputational techniques for Breast cancer drug discovery via the HER-2, BCRA, VEGF and ER protein targets

Computational modelling remains an indispensable technique in drug discovery. With myriad of high computing resources, and improved modelling algorithms, there has been a high-speed in the drug development cycle with promising success rate compared to the traditional route. For example, lapatinib; a well-known anticancer drug with clinical applications was discovered with computational drug design techniques. Similarly, molecular modelling has been applied to various disease areas ranging from cancer to neurodegenerative diseases. The techniques ranges from high-throughput virtual screening, molecular mechanics with generalized Born and surface area solvation (MM/GBSA) to molecular dynamics simulation. This review focuses on the application of computational modelling tools in the identification of drug candidates for Breast cancer. First, we begin with a succinct overview of molecular modelling in the drug discovery process. Next, we take note of special efforts on the developments and applications of combining these techniques with particular emphasis on possible breast cancer therapeutic targets such as estrogen receptor (ER), human epidermal growth factor receptor 2 (HER2), vascular endothelial growth factor (VEGF), breast cancer gene 1 (BRCA1), and breast cancer gene 2 (BRCA2). Finally, we discussed the search for covalent inhibitors against these receptors using computational techniques, advances, pitfalls, possible solutions, and future perspectives.

In silico drug repurposing carvedilol and its metabolites against SARS-CoV-2 infection using molecular docking and molecular dynamic simulation approaches

The pandemic of coronavirus disease 2019 (COVID-19) caused by the infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a significant impact on the economy and public health worldwide. Therapeutic options such as drugs and vaccines for this newly emerged disease are eagerly desired due to the high mortality. Using the U.S. Food and Drug Administration (FDA) approved drugs to treat a new disease or entirely different diseases, in terms of drug repurposing, minimizes the time and cost of drug development compared to the de novo design of a new drug. Drug repurposing also has some other advantages such as reducing safety evaluation to accelerate drug application on time. Carvedilol, a non-selective beta-adrenergic blocker originally designed to treat high blood pressure and manage heart disease, has been shown to impact SARS-CoV-2 infection in clinical observation and basic studies. Here, we applied computer-aided approaches to investigate the possibility of repurposing carvedilol to combat SARS-CoV-2 infection. The molecular mechanisms and potential molecular targets of carvedilol were identified by evaluating the interactions of carvedilol with viral proteins. Additionally, the binding affinities of in vivo metabolites of carvedilol with selected targets were evaluated. The docking scores for carvedilol and its metabolites with RdRp were - 10.0 kcal/mol, - 9.8 kcal/mol (1-hydroxyl carvedilol), - 9.7 kcal/mol (3-hydroxyl carvedilol), - 9.8 kcal/mol (4-hydroxyl carvedilol), - 9.7 kcal/mol (5-hydroxyl carvedilol), - 10.0 kcal/mol (8-hydroxyl carvedilol), and - 10.1 kcal/mol (O-desmethyl carvedilol), respectively. Using the molecular dynamics simulation (100 ns) method, we further confirmed the stability of formed complexes of RNA-dependent RNA polymerase (RdRp) and carvedilol or its metabolites. Finally, the drug-target interaction mechanisms that contribute to the complex were investigated. Overall, this study provides the molecular targets and mechanisms of carvedilol and its metabolites as repurposed drugs to fight against SARS-CoV-2 infection.

G protein-coupled receptor-targeting antibody-drug conjugates: Current status and future directions

In recent years, antibody-drug conjugates (ADCs) have emerged as promising anti-cancer therapeutic agents with several having already received market approval for the treatment of solid tumor and hematological malignancies. As ADC technology continues to improve and the range of indications treatable by ADCs increases, the repertoire of target antigens has expanded and will undoubtedly continue to grow. G protein-coupled receptors (GPCRs) are well-characterized therapeutic targets implicated in many human pathologies, including cancer, and represent a promising emerging target of ADCs. In this review, we will discuss the past and present therapeutic targeting of GPCRs and describe ADCs as therapeutic modalities. Moreover, we will summarize the status of existing preclinical and clinical GPCR-targeted ADCs and address the potential of GPCRs as novel targets for future ADC development.

Updates on drug designing approach through computational strategies: a review

The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.

Hotspot Identification and Drug Design of Protein-Protein Interaction Modulators Using the Fragment Molecular Orbital Method

Protein-protein interactions (PPIs) are essential for the function of many proteins. Aberrant PPIs have the potential to lead to disease, making PPIs promising targets for drug discovery. There are over 64,000 PPIs in the human interactome reference database; however, to date, very few PPI modulators have been approved for clinical use. Further development of PPI-specific therapeutics is highly dependent on the availability of structural data and the existence of reliable computational tools to explore the interface between two interacting proteins. The fragment molecular orbital (FMO) quantum mechanics method offers comprehensive and computationally inexpensive means of identifying the strength (in kcal/mol) and the chemical nature (electrostatic or hydrophobic) of the molecular interactions taking place at the protein-protein interface. We have integrated FMO and PPI exploration (FMO-PPI) to identify the residues that are critical for protein-protein binding (hotspots). To validate this approach, we have applied FMO-PPI to a dataset of protein-protein complexes representing several different protein subfamilies and obtained FMO-PPI results that are in agreement with published mutagenesis data. We observed that critical PPIs can be divided into three major categories: interactions between residues of two proteins (intermolecular), interactions between residues within the same protein (intramolecular), and interactions between residues of two proteins that are mediated by water molecules (water bridges). We extended our findings by demonstrating how this information obtained by FMO-PPI can be utilized to support the structure-based drug design of PPI modulators (SBDD-PPI).

Lipid Nanodiscs for High-Resolution NMR Studies of Membrane Proteins

Membrane proteins (MPs) play essential roles in numerous cellular processes. Because around 70% of the currently marketed drugs target MPs, a detailed understanding of their structure, binding properties, and functional dynamics in a physiologically relevant environment is crucial for a more detailed understanding of this important protein class. We here summarize the benefits of using lipid nanodiscs for NMR structural investigations and provide a detailed overview of the currently used lipid nanodisc systems as well as their applications in solution-state NMR. Despite the increasing use of other structural methods for the structure determination of MPs in lipid nanodiscs, solution NMR turns out to be a versatile tool to probe a wide range of MP features, ranging from the structure determination of small to medium-sized MPs to probing ligand and partner protein binding as well as functionally relevant dynamical signatures in a lipid nanodisc setting. We will expand on these topics by discussing recent NMR studies with lipid nanodiscs and work out a key workflow for optimizing the nanodisc incorporation of an MP for subsequent NMR investigations. With this, we hope to provide a comprehensive background to enable an informed assessment of the applicability of lipid nanodiscs for NMR studies of a particular MP of interest.

XRRpred: Accurate Predictor of Crystal Structure Quality from Protein Sequence

MOTIVATION

X-ray crystallography was used to produce nearly 90% of protein structures. These efforts were supported by numerous sequence-based tools that accurately predict crystallizable proteins. However, protein structures vary widely in their quality, typically measured with resolution and R-free. This impacts the ability to use these structures for some applications including rational drug design and molecular docking and motivates development of methods that accurately predict structure quality.

RESULTS

We introduce XRRpred, the first predictor of the resolution and R-free values from protein sequences. XRRpred relies on original sequence profiles, hand-crafted features, empirically selected and parametrized regressors, and modern resampling techniques. Using an independent test dataset, we show that XRRpred provides accurate predictions of resolution and R-free. We demonstrate that XRRpred's predictions correctly model relationship between the resolution and R-free and reproduce structure quality relations between structural classes of proteins. We also show that XRRpred significantly outperforms indirect alternative ways to predict the structure quality that include predictors of crystallization propensity and an alignment-based approach. XRRpred is available as a convenient webserver that allows batch predictions and offers informative visualization of the results.

AVAILABILITY

http://biomine.cs.vcu.edu/servers/XRRPred/.

Guiding Medicinal Chemistry with Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between a protein and a small molecule plays a key role in the rationalization of potency and selectivity and in design of new ideas. However, even when a target of interest is structurally enabled, visual inspection and force field-based molecular mechanics calculations cannot always explain the full complexity of the molecular interactions that are critical in drug design. Quantum mechanical methods have the potential to address this shortcoming, but traditionally, computational expense has made the application of these calculations impractical. The fragment molecular orbital (FMO) method offers a solution that combines accuracy, speed, and the ability to characterize important interactions (i.e. its strength in kcal/mol and chemical nature: hydrophobic, electrostatic, etc) that would otherwise be hard to detect. In this chapter, we describe the FMO method and illustrate its application in the discovery of the benzothiazole (BZT) series as novel tyrosine kinase ITK inhibitors for treatment of allergic asthma.

Analyzing GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

G-protein-coupled receptors (GPCRs) have enormous physiological and biomedical importance, and therefore it is not surprising that they are the targets of many prescribed drugs. Further progress in GPCR drug discovery is highly dependent on the availability of protein structural information. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanics (QM) approaches are often too computationally expensive to be of practical use in time-sensitive situations, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed, and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule toward ligand binding, including an analysis of their chemical nature. Such information is essential for an efficient structure-based drug design (SBDD) process. In this chapter, we describe how to use FMO in the characterization of GPCR-ligand interactions.

Direct ligand screening against membrane proteins on live cells enabled by DNA-programmed affinity labelling

DNA-programmed affinity labelling (DPAL) enables the screening of chemical compounds against membrane proteins directly on live cells.

Pharmacokinetic and molecular docking studies to design antimalarial compounds targeting Actin I

OBJECTIVE

Malaria is an ancient disease that still causes more than 200 million of cases 7 with high mortality globally. Identification of new drug targets and development of novel antimalarial drugs with unique mode of action encounter the drug resistance and reduce the mortality by Plasmodium parasites. Actin protein is one of the key proteins in Plasmodium falciparum playing multifarious important roles including transport, cell motility, cell division, and shape determination. This study investigated Actin I as a drug target, in silico screening of diverse molecules through molecular docking was considered. Further, pharmacokinetic parameters of the selected molecules from the docking and interaction studies were planned to propose the lead molecules.b.

METHODS

Molecules were selected according to score and protein ligand interaction and selected molecules were subjected for pharmacokinetic studies to investigate important drug parameters.

RESULTS

The docked molecules were ranked according to the binding score and good interaction pattern was observed with Actin I within top 20 scoring molecules. The selected molecules also had optimum pharmacokinetic parameters.

CONCLUSION

The current study provides a set of hit molecules which can be further explored through in vitro and in vivo experiments for the development of potential drugs against malaria, there by encountering drug resistance and establishing Actin I as an important drug target.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Identifying GPCR-drug interaction based on wordbook learning from sequences

BACKGROUND

G protein-coupled receptors (GPCRs) mediate a variety of important physiological functions, are closely related to many diseases, and constitute the most important target family of modern drugs. Therefore, the research of GPCR analysis and GPCR ligand screening is the hotspot of new drug development. Accurately identifying the GPCR-drug interaction is one of the key steps for designing GPCR-targeted drugs. However, it is prohibitively expensive to experimentally ascertain the interaction of GPCR-drug pairs on a large scale. Therefore, it is of great significance to predict the interaction of GPCR-drug pairs directly from the molecular sequences. With the accumulation of known GPCR-drug interaction data, it is feasible to develop sequence-based machine learning models for query GPCR-drug pairs.

RESULTS

In this paper, a new sequence-based method is proposed to identify GPCR-drug interactions. For GPCRs, we use a novel bag-of-words (BoW) model to extract sequence features, which can extract more pattern information from low-order to high-order and limit the feature space dimension. For drug molecules, we use discrete Fourier transform (DFT) to extract higher-order pattern information from the original molecular fingerprints. The feature vectors of two kinds of molecules are concatenated and input into a simple prediction engine distance-weighted K-nearest-neighbor (DWKNN). This basic method is easy to be enhanced through ensemble learning. Through testing on recently constructed GPCR-drug interaction datasets, it is found that the proposed methods are better than the existing sequence-based machine learning methods in generalization ability, even an unconventional method in which the prediction performance was further improved by post-processing procedure (PPP).

CONCLUSIONS

The proposed methods are effective for GPCR-drug interaction prediction, and may also be potential methods for other target-drug interaction prediction, or protein-protein interaction prediction. In addition, the new proposed feature extraction method for GPCR sequences is the modified version of the traditional BoW model and may be useful to solve problems of protein classification or attribute prediction. The source code of the proposed methods is freely available for academic research at https://github.com/wp3751/GPCR-Drug-Interaction.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Toward G protein-coupled receptor structure-based drug design using X-ray lasers

Rational structure-based drug design (SBDD) relies on the availability of a large number of co-crystal structures to map the ligand-binding pocket of the target protein and use this information for lead-compound optimization via an iterative process. While SBDD has proven successful for many drug-discovery projects, its application to G protein-coupled receptors (GPCRs) has been limited owing to extreme difficulties with their crystallization. Here, a method is presented for the rapid determination of multiple co-crystal structures for a target GPCR in complex with various ligands, taking advantage of the serial femtosecond crystallography approach, which obviates the need for large crystals and requires only submilligram quantities of purified protein. The method was applied to the human β2-adrenergic receptor, resulting in eight room-temperature co-crystal structures with six different ligands, including previously unreported structures with carvedilol and propranolol. The generality of the proposed method was tested with three other receptors. This approach has the potential to enable SBDD for GPCRs and other difficult-to-crystallize membrane proteins.

The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization

BACKGROUND

Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug.

METHODS

In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan.

RESULTS

A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given.

CONCLUSION

The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.

Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.

Recent Insights from Molecular Dynamics Simulations for G Protein-Coupled Receptor Drug Discovery

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. Some ligands that bind to GPCRs activate different downstream signaling pathways. G protein activation, or -arrestin biased signaling, involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. -arrestin biased signaling has become a hot target for structure-based drug discovery. However, challenges include that there are few crystal structures available in the Protein Data Bank and that GPCRs are highly dynamic. Hence, molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions with receptors and ligands. Here, we highlight recent MD simulation studies and enhanced sampling methods used to study biased G protein-coupled receptor signaling and their conformational dynamics as well as applications to drug discovery.

Exploring the Potential of Spherical Harmonics and PCVM for Compounds Activity Prediction

Biologically active chemical compounds may provide remedies for several diseases. Meanwhile, Machine Learning techniques applied to Drug Discovery, which are cheaper and faster than wet-lab experiments, have the capability to more effectively identify molecules with the expected pharmacological activity. Therefore, it is urgent and essential to develop more representative descriptors and reliable classification methods to accurately predict molecular activity. In this paper, we investigate the potential of a novel representation based on Spherical Harmonics fed into Probabilistic Classification Vector Machines classifier, namely SHPCVM, to compound the activity prediction task. We make use of representation learning to acquire the features which describe the molecules as precise as possible. To verify the performance of SHPCVM ten-fold cross-validation tests are performed on twenty-one G protein-coupled receptors (GPCRs). Experimental outcomes (accuracy of 0.86) assessed by the classification accuracy, precision, recall, Matthews’ Correlation Coefficient and Cohen’s kappa reveal that using our Spherical Harmonics-based representation which is relatively short and Probabilistic Classification Vector Machines can achieve very satisfactory performance results for GPCRs.

Diverse GPCRs exhibit conserved water networks for stabilization and activation

G protein-coupled receptors (GPCRs) represent both the largest class of drug targets and the largest family of human membrane proteins. Recent three-dimensional structures have revealed the presence of many water molecules buried inside GPCRs, but the functional role of these water molecules has been unclear. Using extensive atomic-level computer simulations, we find that, although most of these waters are highly mobile, a few are stable. These stable water molecules mediate state-dependent and state-independent polar networks that are conserved across diverse GPCRs. These findings promise to help guide the rational design of GPCR-targeted drugs.

G protein-coupled receptors (GPCRs) have evolved to recognize incredibly diverse extracellular ligands while sharing a common architecture and structurally conserved intracellular signaling partners. It remains unclear how binding of diverse ligands brings about GPCR activation, the common structural change that enables intracellular signaling. Here, we identify highly conserved networks of water-mediated interactions that play a central role in activation. Using atomic-level simulations of diverse GPCRs, we show that most of the water molecules in GPCR crystal structures are highly mobile. Several water molecules near the G protein-coupling interface, however, are stable. These water molecules form two kinds of polar networks that are conserved across diverse GPCRs: (i) a network that is maintained across the inactive and the active states and (ii) a network that rearranges upon activation. Comparative analysis of GPCR crystal structures independently confirms the striking conservation of water-mediated interaction networks. These conserved water-mediated interactions near the G protein-coupling region, along with diverse water-mediated interactions with extracellular ligands, have direct implications for structure-based drug design and GPCR engineering.

A first attempt investigation on crystallization screening and crystal quality of lysozyme under different simulated gravities in a large-gradient magnetic field

A magnetic field has been proved useful in protein crystallization in that it can help to improve the crystal quality, which is essential for high-resolution diffraction using crystallography.

GPCR drug discovery: integrating solution NMR data with crystal and cryo-EM structures

The 826 G protein-coupled receptors (GPCRs) in the human proteome regulate key physiological processes and thus have long been attractive drug targets. With the crystal structures of more than 50 different human GPCRs determined over the past decade, an initial platform for structure-based rational design has been established for drugs that target GPCRs, which is currently being augmented with cryo-electron microscopy (cryo-EM) structures of higher-order GPCR complexes. Nuclear magnetic resonance (NMR) spectroscopy in solution is one of the key approaches for expanding this platform with dynamic features, which can be accessed at physiological temperature and with minimal modification of the wild-type GPCR covalent structures. Here, we review strategies for the use of advanced biochemistry and NMR techniques with GPCRs, survey projects in which crystal or cryo-EM structures have been complemented with NMR investigations and discuss the impact of this integrative approach on GPCR biology and drug discovery.

Membrane protein engineering to the rescue

The inherent hydrophobicity of membrane proteins is a major barrier to membrane protein research and understanding. Their low stability and solubility in aqueous environments coupled with poor expression levels make them a challenging area of research. For many years, the only way of working with membrane proteins was to optimise the environment to suit the protein, through the use of different detergents, solubilising additives, and other adaptations. However, with innovative protein engineering methodologies, the membrane proteins themselves are now being adapted to suit the environment. This mini-review looks at the types of adaptations which are applied to membrane proteins from a variety of different fields, including water solubilising fusion tags, thermostabilising mutation screening, scaffold proteins, stabilising protein chimeras, and isolating water-soluble domains.

Taxonomic Landscape of the Dark Proteomes: Whole-Proteome Scale Interplay Between Structural Darkness, Intrinsic Disorder, and Crystallization Propensity

Growth rate of the protein sequence universe dramatically exceeds the speed of expansion for the protein structure universe, generating an immense dark proteome that includes proteins with unknown structure. A whole-proteome scale analysis of 5.4 million proteins from 987 proteomes in the three domains of life and viruses to systematically dissect an interplay between structural coverage, degree of putative intrinsic disorder, and predicted propensity for structure determination is performed. It has been found that Archaean and Bacterial proteomes have relatively high structural coverage and low amounts of disorder, whereas Eukaryotic and Viral proteomes are characterized by a broad spread of structural coverage and higher disorder levels. The analysis reveals that dark proteomes (i.e., proteomes containing high fractions of proteins with unknown structure) have significantly elevated amounts of intrinsic disorder and are predicted to be difficult to solve structurally. Although the majority of dark proteomes are of viral origin, many dark viral proteomes have at least modest crystallization propensity and only a handful of them are enriched in the intrinsic disorder. The disorder, structural coverage, and propensity are mapped for structural determination onto a novel proteome-level sequence similarity network to analyze the interplay of these characteristics in the taxonomic landscape.

Role of simple descriptors and applicability domain in predicting change in protein thermostability

The melting temperature (Tm) of a protein is the temperature at which half of the protein population is in a folded state. Therefore, Tm is a measure of the thermostability of a protein. Increasing the Tm of a protein is a critical goal in biotechnology and biomedicine. However, predicting the change in melting temperature (dTm) due to mutations at a single residue is difficult because it depends on an intricate balance of forces. Existing methods for predicting dTm have had similar levels of success using generally complex models. We find that training a machine learning model with a simple set of easy to calculate physicochemical descriptors describing the local environment of the mutation performed as well as more complicated machine learning models and is 2-6 orders of magnitude faster. Importantly, unlike in most previous publications, we perform a blind prospective test on our simple model by designing 96 variants of a protein not in the training set. Results from retrospective and prospective predictions reveal the limited applicability domain of each model. This study highlights the current deficiencies in the available dTm dataset and is a call to the community to systematically design a larger and more diverse experimental dataset of mutants to prospectively predict dTm with greater certainty.

Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets

To investigate large library docking's ability to find molecules with joint activity against on-targets and selectivity versus antitargets, the dopamine D2 and serotonin 5-HT2A receptors were targeted, seeking selectivity against the histamine H1 receptor. In a second campaign, κ-opioid receptor ligands were sought with selectivity versus the μ-opioid receptor. While hit rates ranged from 40% to 63% against the on-targets, they were just as good against the antitargets, even though the molecules were selected for their putative lack of binding to the off-targets. Affinities, too, were often as good or better for the off-targets. Even though it was occasionally possible to find selective molecules, such as a mid-nanomolar D2/5-HT2A ligand with 21-fold selectivity versus the H1 receptor, this was the exception. Whereas false-negatives are tolerable in docking screens against on-targets, they are intolerable against antitargets; addressing this problem may demand new strategies in the field.

Tuning Vector Stability and Integration Frequency Elevates Functional GPCR Production and Homogeneity in Saccharomyces cerevisiae

Membrane proteins play a valuable role in biotechnology, yet the difficulty of producing high yields of functional membrane protein limits their use in synthetic biology. The practical application of G protein-coupled receptors in whole cell biosensors, for example, is restricted to those that are functionally produced at the cell surface in the chosen host, limiting the range of detectable molecules. Here, we present a facile approach to significantly improve the yield and homogeneity of functional membrane proteins in Saccharomyces cerevisiae by altering only the choice of expression vector. Expression of a model GPCR, the human adenosine A₂a receptor, from commonly used centromeric and episomal vectors leads to low yields and cellular heterogeneity due to plasmid loss in 20-90% of the cell population. In contrast, homogeneous production of GPCR is attained using a multisite integrating vector or a novel, modified high copy vector that does not require genomic integration or addition of any selection agents. Finally, we introduce a FACS-based screen, which enables rapid isolation of cells with 4- to 15-fold increases in gene dosage and up to a 9-fold increase in functional protein yield without loss of homogeneity compared to a strain isolated through conventional, low-throughput methods. These results can be extended to improve the cellular homogeneity and yield of other membrane proteins, expanding the repertoire of useful receptors for synthetic biology applications.

Functional and structural characterization of osteocytic MLO-Y4 cell proteins encoded by genes differentially expressed in response to mechanical signals in vitro

The anabolic response of bone to mechanical load is partially the result of osteocyte response to fluid flow-induced shear stress. Understanding signaling pathways activated in osteocytes exposed to fluid flow could identify novel signaling pathways involved in the response of bone to mechanical load. Bioinformatics allows for a unique perspective and provides key first steps in understanding these signaling pathways. We examined proteins encoded by genes differentially expressed in response to fluid flow in murine osteocytic MLO-Y4 cells. We considered structural and functional characteristics including putative intrinsic disorder, evolutionary conservation, interconnectedness in protein-protein interaction networks, and cellular localization. Our analysis suggests that proteins encoded by fluid flow activated genes have lower than expected conservation, are depleted in intrinsic disorder, maintain typical levels of connectivity for the murine proteome, and are found in the cytoplasm and extracellular space. Pathway analyses reveal that these proteins are associated with cellular response to stress, chemokine and cytokine activity, enzyme binding, and osteoclast differentiation. The lower than expected disorder of proteins encoded by flow activated genes suggests they are relatively specialized.

Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode

An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A_2A and A₁ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A_2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A_2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A_2A/A₃, ARs, or A₃AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A_2A/A₃, A₁/A₃, A₁/A_2A/A₃, and A₃). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A₃ receptor.

Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design

The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the "golden age for GPCR structural biology." Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand- and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Molecular Evidence of Adenosine Deaminase Linking Adenosine A2A Receptor and CD26 Proteins

Adenosine is an endogenous purine nucleoside that acts in all living systems as a homeostatic network regulator through many pathways, which are adenosine receptor (AR)-dependent and -independent. From a metabolic point of view, adenosine deaminase (ADA) is an essential protein in the regulation of the total intracellular and extracellular adenosine in a tissue. In addition to its cytosolic localization, ADA is also expressed as an ecto-enzyme on the surface of different cells. Dipeptidyl peptidase IV (CD26) and some ARs act as binding proteins for extracellular ADA in humans. Since CD26 and ARs interact with ADA at opposite sites, we have investigated if ADA can function as a cell-to-cell communication molecule by bridging the anchoring molecules CD26 and A2AR present on the surfaces of the interacting cells. By combining site-directed mutagenesis of ADA amino acids involved in binding to A2AR and a modification of the bioluminescence resonance energy transfer (BRET) technique that allows detection of interactions between two proteins expressed in different cell populations with low steric hindrance (NanoBRET), we show direct evidence of the specific formation of trimeric complexes CD26-ADA-A2AR involving two cells. By dynamic mass redistribution assays and ligand binding experiments, we also demonstrate that A2AR-NanoLuc fusion proteins are functional. The existence of this ternary complex is in good agreement with the hypothesis that ADA could bridge T-cells (expressing CD26) and dendritic cells (expressing A2AR). This is a new metabolic function for ecto-ADA that, being a single chain protein, it has been considered as an example of moonlighting protein, because it performs more than one functional role (as a catalyst, a costimulator, an allosteric modulator and a cell-to-cell connector) without partitioning these functions in different subunits.

Computational Methods Used in Hit-to-Lead and Lead Optimization Stages of Structure-Based Drug Discovery

GPCR modeling approaches are widely used in the hit-to-lead (H2L) and lead optimization (LO) stages of drug discovery. The aims of these modeling approaches are to predict the 3D structures of the receptor-ligand complexes, to explore the key interactions between the receptor and the ligand and to utilize these insights in the design of new molecules with improved binding, selectivity or other pharmacological properties. In this book chapter, we present a brief survey of key computational approaches integrated with hierarchical GPCR modeling protocol (HGMP) used in hit-to-lead (H2L) and in lead optimization (LO) stages of structure-based drug discovery (SBDD). We outline the differences in modeling strategies used in H2L and LO of SBDD and illustrate how these tools have been applied in three drug discovery projects.

Survey of Predictors of Propensity for Protein Production and Crystallization with Application to Predict Resolution of Crystal Structures

Selection of proper targets for the X-ray crystallography will benefit biological research community immensely. Several computational models were proposed to predict propensity of successful protein production and diffraction quality crystallization from protein sequences. We reviewed a comprehensive collection of 22 such predictors that were developed in the last decade. We found that almost all of these models are easily accessible as webservers and/or standalone software and we demonstrated that some of them are widely used by the research community. We empirically evaluated and compared the predictive performance of seven representative methods. The analysis suggests that these methods produce quite accurate propensities for the diffraction-quality crystallization. We also summarized results of the first study of the relation between these predictive propensities and the resolution of the crystallizable proteins. We found that the propensities predicted by several methods are significantly higher for proteins that have high resolution structures compared to those with the low resolution structures. Moreover, we tested a new meta-predictor, MetaXXC, which averages the propensities generated by the three most accurate predictors of the diffraction-quality crystallization. MetaXXC generates putative values of resolution that have modest levels of correlation with the experimental resolutions and it offers the lowest mean absolute error when compared to the seven considered methods. We conclude that protein sequences can be used to fairly accurately predict whether their corresponding protein structures can be solved using X-ray crystallography. Moreover, we also ascertain that sequences can be used to reasonably well predict the resolution of the resulting protein crystals.

Exploring GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity. It is essential for an efficient structure-based drug design (SBDD) process. FMO enables ab initio approaches to be applied to systems that conventional quantum-mechanical (QM) methods would find challenging. The key advantage of the Fragment Molecular Orbital Method (FMO) is that it can reveal atomistic details about the individual contributions and chemical nature of each residue and water molecule toward ligand binding which would otherwise be difficult to detect without using QM methods. In this chapter, we demonstrate the typical use of FMO to analyze 19 crystal structures of β1 and β2 adrenergic receptors with their corresponding agonists and antagonists.

Human Adenosine A2A Receptor: Molecular Mechanism of Ligand Binding and Activation

Adenosine receptors (ARs) comprise the P1 class of purinergic receptors and belong to the largest family of integral membrane proteins in the human genome, the G protein-coupled receptors (GPCRs). ARs are classified into four subtypes, A₁, A_2A, A_2B, and A₃, which are all activated by extracellular adenosine, and play central roles in a broad range of physiological processes, including sleep regulation, angiogenesis and modulation of the immune system. ARs are potential therapeutic targets in a variety of pathophysiological conditions, including sleep disorders, cancer, and dementia, which has made them important targets for structural biology. Over a decade of research and innovation has culminated with the publication of more than 30 crystal structures of the human adenosine A_2A receptor (A_2AR), making it one of the best structurally characterized GPCRs at the atomic level. In this review we analyze the structural data reported for A_2AR that described for the first time the binding of mode of antagonists, including newly developed drug candidates, synthetic and endogenous agonists, sodium ions and an engineered G protein. These structures have revealed the key conformational changes induced upon agonist and G protein binding that are central to signal transduction by A_2AR, and have highlighted both similarities and differences in the activation mechanism of this receptor compared to other class A GPCRs. Finally, comparison of A_2AR with the recently solved structures of A₁R has provided the first structural insight into the molecular determinants of ligand binding specificity in different AR subtypes.

Trends in GPCR drug discovery: new agents, targets and indications

G protein-coupled receptors (GPCRs) are the most intensively studied drug targets, mostly due to their substantial involvement in human pathophysiology and their pharmacological tractability. Here, we report an up-to-date analysis of all GPCR drugs and agents in clinical trials, which reveals current trends across molecule types, drug targets and therapeutic indications, including showing that 475 drugs (~34% of all drugs approved by the US Food and Drug Administration (FDA)) act at 108 unique GPCRs. Approximately 321 agents are currently in clinical trials, of which ~20% target 66 potentially novel GPCR targets without an approved drug, and the number of biological drugs, allosteric modulators and biased agonists has increased. The major disease indications for GPCR modulators show a shift towards diabetes, obesity and Alzheimer disease, although several central nervous system disorders are also highly represented. The 224 (56%) non-olfactory GPCRs that have not yet been explored in clinical trials have broad untapped therapeutic potential, particularly in genetic and immune system disorders. Finally, we provide an interactive online resource to analyse and infer trends in GPCR drug discovery.

Structure-Based Design and Discovery of New M2 Receptor Agonists

Muscarinic receptor agonists are characterized by apparently strict restraints on their tertiary or quaternary amine and their distance to an ester or related center. On the basis of the active state crystal structure of the muscarinic M₂ receptor in complex with iperoxo, we explored potential agonists that lacked the highly conserved functionalities of previously known ligands. Using structure-guided pharmacophore design followed by docking, we found two agonists (compounds 3 and 17), out of 19 docked and synthesized compounds, that fit the receptor well and were predicted to form a hydrogen-bond conserved among known agonists. Structural optimization led to compound 28, which was 4-fold more potent than its parent 3. Fortified by the discovery of this new scaffold, we sought a broader range of chemotypes by docking 2.2 million fragments, which revealed another three micromolar agonists unrelated either to 28 or known muscarinics. Even pockets as tightly defined and as deeply studied as that of the muscarinic reveal opportunities for the structure-based design and the discovery of new chemotypes.

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.

Understanding the common themes and diverse roles of the second extracellular loop (ECL2) of the GPCR super-family

The extracellular loops (ECLs) of G protein-coupled receptors (GPCRs) can bind directly to docked orthosteric or allosteric ligands, they can contain transient contact points for ligand entry into the transmembrane (TM) bundle and they can regulate the activation of the receptor signalling pathways. Of the three ECLs, ECL2 is the largest and most structurally diverse reflecting its functional importance. This has been shown through biochemical techniques and has been supported by the many subsequent crystal structures of GPCRs bound to both agonists and antagonists. ECL2 shares common structural features between (and sometimes across) receptor sub-families and can facilitate ligand entry to the TM core or act directly as a surface of the ligand-binding pocket. Structural similarities seem to underpin common binding mechanisms; however, where these exist, variations in primary sequence ensure ligand-binding specificity. This review will compare current understanding of the structural themes and main functional roles of ECL2 in ligand binding, activation and regulation of the major families of GPCRs.

Rapid and accurate assessment of GPCR-ligand interactions Using the fragment molecular orbital-based density-functional tight-binding method

The reliable and precise evaluation of receptor–ligand interactions and pair‐interaction energy is an essential element of rational drug design. While quantum mechanical (QM) methods have been a promising means by which to achieve this, traditional QM is not applicable for large biological systems due to its high computational cost. Here, the fragment molecular orbital (FMO) method has been used to accelerate QM calculations, and by combining FMO with the density‐functional tight‐binding (DFTB) method we are able to decrease computational cost 1000 times, achieving results in seconds, instead of hours. We have applied FMO‐DFTB to three different GPCR–ligand systems. Our results correlate well with site directed mutagenesis data and findings presented in the published literature, demonstrating that FMO‐DFTB is a rapid and accurate means of GPCR–ligand interactions. © 2017 Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Applying Structure-Based Drug Design Approaches to Allosteric Modulators of GPCRs

Structural insights have been revealed from X-ray co-complexes of a range of G protein-coupled receptors (GPCRs) and their allosteric ligands. The understanding of how small molecules can modulate the function of this important class of receptors by binding to a diverse range of pockets on and inside the proteins has had a profound impact on the structure-based drug design (SBDD) of new classes of therapeutic agents. The types of allosteric pockets and the mode of modulation as well as the advantages and disadvantages of targeting allosteric pockets (as opposed to the natural orthosteric site) are considered in the context of these new structural findings.

Photoaffinity Labeling of the Human A2A Adenosine Receptor and Cross-link Position Analysis by Mass Spectrometry

Photoaffinity labeling (PAL) is widely used for the identification of ligand-binding proteins and elucidation of ligand-binding sites. PAL has also been employed for the characterization of G protein-coupled receptors (GPCRs); however, a limited number of reports has successfully identified their cross-linked amino acids. This report is the first of its kind to determine the cross-link position of the human A_2A adenosine receptor by PAL with the novel diazirine-based photoaffinity probe 9.

Supporting the Identification of Novel Fragment-Based Positive Allosteric Modulators Using a Supervised Molecular Dynamics Approach: A Retrospective Analysis Considering the Human A2A Adenosine Receptor as a Key Example

Structure-driven fragment-based (SDFB) approaches have provided efficient methods for the identification of novel drug candidates. This strategy has been largely applied in discovering several pharmacological ligand classes, including enzyme inhibitors, receptor antagonists and, more recently, also allosteric (positive and negative) modulators. Recently, Siegal and collaborators reported an interesting study, performed on a detergent-solubilized StaR adenosine A_2A receptor, describing the existence of both fragment-like negative allosteric modulators (NAMs), and fragment-like positive allosteric modulators (PAMs). From this retrospective study, our results suggest that Supervised Molecular Dynamics (SuMD) simulations can support, on a reasonable time scale, the identification of fragment-like PAMs following their receptor recognition pathways and characterizing the possible allosteric binding sites.

Using the fragment molecular orbital method to investigate agonist-orexin-2 receptor interactions

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in 'real-time' and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allowsab initioapproaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure-function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available.