AlphaFold2 versus experimental structures: evaluation on G protein-coupled receptors

Physical stability of oil-in-water multi-layered coenzyme Q10 nano-emulsions

As a lipophilic antioxidant, coenzyme Q10 (CoQ10) has limited application owing to its low water solubility and instability. In the present study, potato protein (PP) and soybean soluble polysaccharide (SSPS) were used as carriers to prepare a multilayer SSPS-PP-CoQ10 nano-emulsion using the reversed-phase emulsification method; further, the water solubility, stability, and formation mechanism of the nano-emulsion were analyzed. The results showed that the particle size of SSPS-PP-CoQ10 nano-emulsions was 253-422 nm with good polydispersity. The encapsulation efficiency (EE) could reach up to 88.87 %. When the concentration of SSPS was 0.1 wt%, the decrease in interfacial tension and increase in viscoelasticity indicated that nano-emulsion improved CoQ10 physical stability. SSPS incorporation altered the microscopic environment of the hydrophobic residues, rendering them more hydrophilic and enhancing their water solubility. According to molecular docking results, hydrogen bonds promote binding among SSPS, PP, and CoQ10, and increase emulsion stability.

Leveraging Artificial Intelligence in GPCR Activation Studies: Computational Prediction Methods as Key Drivers of Knowledge

G protein-coupled receptors (GPCRs) are key molecules involved in cellular signaling and are attractive targets for pharmacological intervention. This chapter is designed to explore the range of algorithms used to predict GPCRs' activation states, while also examining the pharmaceutical implications of these predictions. Our primary objective is to show how artificial intelligence (AI) is key in GPCR research to reveal the intricate dynamics of activation and inactivation processes, shedding light on the complex regulatory mechanisms of this vital protein family. We describe several computational strategies that leverage diverse structural data from the Protein Data Bank, molecular dynamic simulations, or ligand-based methods to predict the activation states of GPCRs. We demonstrate how the integration of AI into GPCR research not only enhances our understanding of their dynamic properties but also presents immense potential for driving pharmaceutical research and development, offering promising new avenues in the search for newer, better therapeutic agents.

AlphaFold3 versus experimental structures: assessment of the accuracy in ligand-bound G protein-coupled receptors

G protein-coupled receptors (GPCRs) are critical drug targets involved in numerous physiological processes, yet many of their structures remain unresolved due to inherent flexibility and diverse ligand interactions. This study systematically evaluates the accuracy of AlphaFold3-predicted GPCR structures compared to experimentally determined structures, with a primary focus on ligand-bound states. Our analysis reveals that while AlphaFold3 shows improved performance over AlphaFold2 in predicting overall GPCR backbone architecture, significant discrepancies persist in ligand-binding poses, particularly for ions, peptides, and proteins. Despite advancements, these limitations constrain the utility of AlphaFold3 models in functional studies and structure-based drug design, where high-resolution details of ligand interactions are crucial. We assess the accuracy of predicted structures across various ligand types, quantifying deviations in binding pocket geometries and ligand orientations. Our findings highlight specific challenges in the computational prediction of ligand-bound GPCR structures, emphasizing areas where further refinement is needed. This study provides valuable insights for researchers using AlphaFold3 in GPCR studies, underscores the ongoing necessity for experimental structure determination, and offers direction for improving protein-ligand interaction predictions in future computational models.

Exploring key features of selectivity in somatostatin receptors through molecular dynamics simulations

Somatostatin receptors (SSTRs) are widely distributed throughout the human body and play crucial roles in various physiological processes. They are recognized as key targets for both radiotherapy and radiodiagnosis due to their overexpression in several cancer types. However, the discovery and design of selective drugs for each of the five isoforms have been significantly hindered by the lack of complete structural information. In this study, we conducted a systematic computational analysis of all five SSTRs in complex with the endogenous ligand somatostatin to elucidate their structural and dynamic features. We thoroughly characterized each isoform using available experimental structures for SSTR2 and SSTR4, as well as AlphaFold2 models for SSTR1, SSTR3, and SSTR5. By performing multi-copy μs-long molecular dynamics simulations, we examined the differences and similarities in dynamical behavior and somatostatin binding among all SSTRs. Our analysis focused on understanding the opening and closing movements of the extracellular loop 2, which are crucial for ligand binding and recognition. Interestingly, we observed a unique conformation of somatostatin within the binding pocket of SSTR5 in which the loop can partially close, as compared to the other isoforms. Fingerprint analyses provided distinct interaction patterns of somatostatin with all receptors, thus enabling precise guidelines for the discovery and development of more selective somatostatin-based pharmaceuticals tailored for precision medicine therapies.

Template-based modeling of insect odorant receptors outperforms AlphaFold3 for ligand binding predictions

Insects rely on odorant receptors (ORs) to detect and respond to volatile environmental cues, so the ORs are attracting increasing interest as potential targets for pest control. However, experimental analysis of their structures and functions faces significant challenges. Computational methods such as template-based modeling (TBM) and AlphaFold3 (AF3) could facilitate the structural characterisation of ORs. This study first showed that both models accurately predicted the structural fold of MhOR5, a jumping bristletail OR with known experimental 3D structures, although accuracy was higher in the extracellular region of the protein and binding mode of their cognate ligands with TBM. The two approaches were then compared for their ability to predict the empirical binding evidence available for OR-odorant complexes in two economically important fruit fly species, Bactrocera dorsalis and B. minax. Post-simulation analyses including binding affinities, complex and ligand stability and receptor-ligand interactions (RLIs) revealed that TBM performed better than AF3 in discriminating between binder and non-binder complexes. TBM's superior performance is attributed to hydrophobicity-based helix-wise multiple sequence alignment (MSA) between available insect OR templates and the ORs for which the binding data were generated. This MSA identified conserved residues and motifs which could be used as anchor points for refining the alignments.

Substrate Specificity of ABCB Transporters Predicted by Docking Simulations Can Be Confirmed by Experimental Tests

ATP-binding cassette (ABC) transporters, particularly those of subfamily B, are involved in cell detoxification, multidrug resistance, drug treatment pharmacodynamics, and also ecological adaptation. In this regard, ABCB transporters may play a decisive role in the co-evolution between plants and herbivores. Cardenolides, toxic steroid glycosides, are secondary plant metabolites that defend plants against herbivores by targeting their sodium-potassium ATPase. Despite their toxicity, several herbivorous insects such as the large milkweed bug (Oncopeltus fasciatus) have evolved adaptations to tolerate cardenolides and sequester them for their own defense. We investigate the role of two ABCB transporters of O. fasciatus for the paracellular transport of cardenolides by docking simulations and ATPase assays. Cardenolide binding of OfABCB1 and OfABCB2 is predicted by docking simulations and calculated binding energies are compared with substrate specificities determined in ATPase assays. Both tested ABCB transporters showed activity upon exposure to cardenolides and Km values that agreed well with the predictions of our docking simulations. We conclude that docking simulations can help identify transporter binding regions and predict substrate specificity, as well as provide deeper insights into the structural basis of ABC transporter function.

GPCR-BSD: a database of binding sites of human G-protein coupled receptors under diverse states

G-protein coupled receptors (GPCRs), the largest family of membrane proteins in human body, involve a great variety of biological processes and thus have become highly valuable drug targets. By binding with ligands (e.g., drugs), GPCRs switch between active and inactive conformational states, thereby performing functions such as signal transmission. The changes in binding pockets under different states are important for a better understanding of drug-target interactions. Therefore it is critical, as well as a practical need, to obtain binding sites in human GPCR structures. We report a database (called GPCR-BSD) that collects 127,990 predicted binding sites of 803 GPCRs under active and inactive states (thus 1,606 structures in total). The binding sites were identified from the predicted GPCR structures by executing three geometric-based pocket prediction methods, fpocket, CavityPlus and GHECOM. The server provides query, visualization, and comparison of the predicted binding sites for both GPCR predicted and experimentally determined structures recorded in PDB. We evaluated the identified pockets of 132 experimentally determined human GPCR structures in terms of pocket residue coverage, pocket center distance and redocking accuracy. The evaluation showed that fpocket and CavityPlus methods performed better and successfully predicted orthosteric binding sites in over 60% of the 132 experimentally determined structures. The GPCR Binding Site database is freely accessible at https://gpcrbs.bigdata.jcmsc.cn . This study not only provides a systematic evaluation of the commonly-used fpocket and CavityPlus methods for the first time but also meets the need for binding site information in GPCR studies.

A deep learning framework combining molecular image and protein structural representations identifies candidate drugs for pain

Artificial intelligence (AI) and deep learning technologies hold promise for identifying effective drugs for human diseases, including pain. Here, we present an interpretable deep-learning-based ligand image- and receptor's three-dimensional (3D)-structure-aware framework to predict compound-protein interactions (LISA-CPI). LISA-CPI integrates an unsupervised deep-learning-based molecular image representation (ImageMol) of ligands and an advanced AlphaFold2-based algorithm (Evoformer). We demonstrated that LISA-CPI achieved ∼20% improvement in the average mean absolute error (MAE) compared to state-of-the-art models on experimental CPIs connecting 104,969 ligands and 33 G-protein-coupled receptors (GPCRs). Using LISA-CPI, we prioritized potential repurposable drugs (e.g., methylergometrine) and identified candidate gut-microbiota-derived metabolites (e.g., citicoline) for potential treatment of pain via specifically targeting human GPCRs. In summary, we presented that the integration of molecular image and protein 3D structural representations using a deep learning framework offers a powerful computational drug discovery tool for treating pain and other complex diseases if broadly applied.

AI Prediction of Structural Stability of Nanoproteins Based on Structures and Residue Properties by Mean Pooled Dual Graph Convolutional Network

The structural stability of proteins is an important topic in various fields such as biotechnology, pharmaceuticals, and enzymology. Specifically, understanding the structural stability of protein is crucial for protein design. Artificial design, while pursuing high thermodynamic stability and rigidity of proteins, inevitably sacrifices biological functions closely related to protein flexibility. The thermodynamic stability of proteins is not always optimal when they are highest to perfectly perform their biological functions. Extensive theoretical and experimental screening is often required to obtain stable protein structures. Thus, it becomes critically important to develop a stability prediction model based on the balance between protein stability and bioactivity. To design protein drugs with better functionality in a broader structural space, a novel protein structural stability predictor called PSSP has been developed in this study. PSSP is a mean pooled dual graph convolutional network (GCN) model based on sequence characteristics and secondary structure, distance matrix, graph, and residue properties of a nanoprotein to provide rapid prediction and judgment. This model exhibits excellent robustness in predicting the structural stability of nanoproteins. Comparing with previous artificial intelligence algorithms, the results indicate this model can provide a rapid and accurate assessment of the structural stability of artificially designed proteins, which shows the great promises for promoting the robust development of protein design.

Utilization of an optimized AlphaFold protein model for structure-based design of a selective HDAC11 inhibitor with anti-neuroblastoma activity

AlphaFold is an artificial intelligence approach for predicting the three-dimensional (3D) structures of proteins with atomic accuracy. One challenge that limits the use of AlphaFold models for drug discovery is the correct prediction of folding in the absence of ligands and cofactors, which compromises their direct use. We have previously described the optimization and use of the histone deacetylase 11 (HDAC11) AlphaFold model for the docking of selective inhibitors such as FT895 and SIS17. Based on the predicted binding mode of FT895 in the optimized HDAC11 AlphaFold model, a new scaffold for HDAC11 inhibitors was designed, and the resulting compounds were tested in vitro against various HDAC isoforms. Compound 5a proved to be the most active compound with an IC50 of 365 nM and was able to selectively inhibit HDAC11. Furthermore, docking of 5a showed a binding mode comparable to FT895 but could not adopt any reasonable poses in other HDAC isoforms. We further supported the docking results with molecular dynamics simulations that confirmed the predicted binding mode. 5a also showed promising activity with an EC50 of 3.6 µM on neuroblastoma cells.

Reliability of AlphaFold2 Models in Virtual Drug Screening: A Focus on Selected Class A GPCRs

Protein three-dimensional (3D) structure prediction is one of the most challenging issues in the field of computational biochemistry, which has overwhelmed scientists for almost half a century. A significant breakthrough in structural biology has been established by developing the artificial intelligence (AI) system AlphaFold2 (AF2). The AF2 system provides a state-of-the-art prediction of protein structures from nearly all known protein sequences with high accuracy. This study examined the reliability of AF2 models compared to the experimental structures in drug discovery, focusing on one of the most common protein drug-targeted classes known as G protein-coupled receptors (GPCRs) class A. A total of 32 representative protein targets were selected, including experimental structures of X-ray crystallographic and Cryo-EM structures and their corresponding AF2 models. The quality of AF2 models was assessed using different structure validation tools, including the pLDDT score, RMSD value, MolProbity score, percentage of Ramachandran favored, QMEAN Z-score, and QMEANDisCo Global. The molecular docking was performed using the Genetic Optimization for Ligand Docking (GOLD) software. The AF2 models' reliability in virtual drug screening was determined by their ability to predict the ligand binding poses closest to the native binding pose by assessing the Root Mean Square Deviation (RMSD) metric and docking scoring function. The quality of the docking and scoring function was evaluated using the enrichment factor (EF). Furthermore, the capability of using AF2 models in molecular docking to identify hits with key protein-ligand interactions was analyzed. The posing power results showed that the AF2 models successfully predicted ligand binding poses (RMSD < 2 Å). However, they exhibited lower screening power, with average EF values of 2.24, 2.42, and 1.82 for X-ray, Cryo-EM, and AF2 structures, respectively. Moreover, our study revealed that molecular docking using AF2 models can identify competitive inhibitors. In conclusion, this study found that AF2 models provided docking results comparable to experimental structures, particularly for certain GPCR targets, and could potentially significantly impact drug discovery.

Exploring the constitutive activation mechanism of the class A orphan GPR20

GPR20, an orphan G protein-coupled receptor (GPCR), shows significant expression in intestinal tissue and represents a potential therapeutic target to treat gastrointestinal stromal tumors. GPR20 performs high constitutive activity when coupling with Gi. Despite the pharmacological importance of GPCR constitutive activation, determining the mechanism has long remained unclear. In this study, we explored the constitutive activation mechanism of GPR20 through large-scale unbiased molecular dynamics simulations. Our results unveil the allosteric nature of constitutively activated GPCR signal transduction involving extracellular and intracellular domains. Moreover, the constitutively active state of the GPR20 requires both the N-terminal cap and Gi protein. The N-terminal cap of GPR20 functions like an agonist and mediates long-range activated conformational shift. Together with the previous study, this study enhances our knowledge of the self-activation mechanism of the orphan receptor, facilitates the drug discovery efforts that target GPR20.

Interpreting forces as deep learning gradients improves quality of predicted protein structures

Protein structure predictions from deep learning models like AlphaFold2, despite their remarkable accuracy, are likely insufficient for direct use in downstream tasks like molecular docking. The functionality of such models could be improved with a combination of increased accuracy and physical intuition. We propose a new method to train deep learning protein structure prediction models using molecular dynamics force fields to work toward these goals. Our custom PyTorch loss function, OpenMM-Loss, represents the potential energy of a predicted structure. OpenMM-Loss can be applied to any all-atom representation of a protein structure capable of mapping into our software package, SidechainNet. We demonstrate our method's efficacy by finetuning OpenFold. We show that subsequently predicted protein structures, both before and after a relaxation procedure, exhibit comparable accuracy while displaying lower potential energy and improved structural quality as assessed by MolProbity metrics.

The path to the G protein-coupled receptor structural landscape: Major milestones and future directions

G protein-coupled receptors (GPCRs) play a crucial role in cell function by transducing signals from the extracellular environment to the inside of the cell. They mediate the effects of various stimuli, including hormones, neurotransmitters, ions, photons, food tastants and odorants, and are renowned drug targets. Advancements in structural biology techniques, including X-ray crystallography and cryo-electron microscopy (cryo-EM), have driven the elucidation of an increasing number of GPCR structures. These structures reveal novel features that shed light on receptor activation, dimerization and oligomerization, dichotomy between orthosteric and allosteric modulation, and the intricate interactions underlying signal transduction, providing insights into diverse ligand-binding modes and signalling pathways. However, a substantial portion of the GPCR repertoire and their activation states remain structurally unexplored. Future efforts should prioritize capturing the full structural diversity of GPCRs across multiple dimensions. To do so, the integration of structural biology with biophysical and computational techniques will be essential. We describe in this review the progress of nuclear magnetic resonance (NMR) to examine GPCR plasticity and conformational dynamics, of atomic force microscopy (AFM) to explore the spatial-temporal dynamics and kinetic aspects of GPCRs, and the recent breakthroughs in artificial intelligence for protein structure prediction to characterize the structures of the entire GPCRome. In summary, the journey through GPCR structural biology provided in this review illustrates how far we have come in decoding these essential proteins architecture and function. Looking ahead, integrating cutting-edge biophysics and computational tools offers a path to navigating the GPCR structural landscape, ultimately advancing GPCR-based applications.

AlphaFold accelerated discovery of psychotropic agonists targeting the trace amine-associated receptor 1

Artificial intelligence is revolutionizing protein structure prediction, providing unprecedented opportunities for drug design. To assess the potential impact on ligand discovery, we compared virtual screens using protein structures generated by the AlphaFold machine learning method and traditional homology modeling. More than 16 million compounds were docked to models of the trace amine-associated receptor 1 (TAAR1), a G protein-coupled receptor of unknown structure and target for treating neuropsychiatric disorders. Sets of 30 and 32 highly ranked compounds from the AlphaFold and homology model screens, respectively, were experimentally evaluated. Of these, 25 were TAAR1 agonists with potencies ranging from 12 to 0.03 μM. The AlphaFold screen yielded a more than twofold higher hit rate (60%) than the homology model and discovered the most potent agonists. A TAAR1 agonist with a promising selectivity profile and drug-like properties showed physiological and antipsychotic-like effects in wild-type but not in TAAR1 knockout mice. These results demonstrate that AlphaFold structures can accelerate drug discovery.

AlphaFold2 in biomedical research: facilitating the development of diagnostic strategies for disease

Proteins, as the primary executors of physiological activity, serve as a key factor in disease diagnosis and treatment. Research into their structures, functions, and interactions is essential to better understand disease mechanisms and potential therapies. DeepMind's AlphaFold2, a deep-learning protein structure prediction model, has proven to be remarkably accurate, and it is widely employed in various aspects of diagnostic research, such as the study of disease biomarkers, microorganism pathogenicity, antigen-antibody structures, and missense mutations. Thus, AlphaFold2 serves as an exceptional tool to bridge fundamental protein research with breakthroughs in disease diagnosis, developments in diagnostic strategies, and the design of novel therapeutic approaches and enhancements in precision medicine. This review outlines the architecture, highlights, and limitations of AlphaFold2, placing particular emphasis on its applications within diagnostic research grounded in disciplines such as immunology, biochemistry, molecular biology, and microbiology.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

OPUS-Rota5: A highly accurate protein side-chain modeling method with 3D-Unet and RotaFormer

Accurate protein side-chain modeling is crucial for protein folding and design. This is particularly true for molecular docking as ligands primarily interact with side chains. In this study, we introduce a two-stage side-chain modeling approach called OPUS-Rota5. It leverages a modified 3D-Unet to capture the local environmental features, including ligand information of each residue, and then employs the RotaFormer module to aggregate various types of features. Evaluation on three test sets, including recently released targets from CAMEO and CASP15, shows that OPUS-Rota5 significantly outperforms some other leading side-chain modeling methods. We also employ OPUS-Rota5 to refine the side chains of 25 G protein-coupled receptor targets predicted by AlphaFold2 and achieve a significantly improved success rate in a subsequent "back" docking of their natural ligands. Therefore, OPUS-Rota5 is a useful and effective tool for molecular docking, particularly for targets with relatively accurate predicted backbones but not side chains such as high-homology targets.

The application of artificial intelligence to accelerate G protein-coupled receptor drug discovery

The application of artificial intelligence (AI) approaches to drug discovery for G protein-coupled receptors (GPCRs) is a rapidly expanding area. Artificial intelligence can be used at multiple stages during the drug discovery process, from aiding our understanding of the fundamental actions of GPCRs to the discovery of new ligand-GPCR interactions or the prediction of clinical responses. Here, we provide an overview of the concepts behind artificial intelligence, including the subfields of machine learning and deep learning. We summarise the published applications of artificial intelligence to different stages of the GPCR drug discovery process. Finally, we reflect on the benefits and limitations of artificial intelligence and share our vision for the exciting potential for further development of applications to aid GPCR drug discovery. In addition to making the drug discovery process "faster, smarter and cheaper," we anticipate that the application of artificial intelligence will create exciting new opportunities for GPCR drug discovery. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

AlphaFold2 structures guide prospective ligand discovery

AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.

Highly accurate carbohydrate-binding site prediction with DeepGlycanSite

As the most abundant organic substances in nature, carbohydrates are essential for life. Understanding how carbohydrates regulate proteins in the physiological and pathological processes presents opportunities to address crucial biological problems and develop new therapeutics. However, the diversity and complexity of carbohydrates pose a challenge in experimentally identifying the sites where carbohydrates bind to and act on proteins. Here, we introduce a deep learning model, DeepGlycanSite, capable of accurately predicting carbohydrate-binding sites on a given protein structure. Incorporating geometric and evolutionary features of proteins into a deep equivariant graph neural network with the transformer architecture, DeepGlycanSite remarkably outperforms previous state-of-the-art methods and effectively predicts binding sites for diverse carbohydrates. Integrating with a mutagenesis study, DeepGlycanSite reveals the guanosine-5'-diphosphate-sugar-recognition site of an important G-protein coupled receptor. These findings demonstrate DeepGlycanSite is invaluable for carbohydrate-binding site prediction and could provide insights into molecular mechanisms underlying carbohydrate-regulation of therapeutically important proteins.

Overview of AlphaFold2 and breakthroughs in overcoming its limitations

Predicting three-dimensional (3D) protein structures has been challenging for decades. The emergence of AlphaFold2 (AF2), a deep learning-based machine learning method developed by DeepMind, became a game changer in the protein folding community. AF2 can predict a protein's three-dimensional structure with high confidence based on its amino acid sequence. Accurate prediction of protein structures can dramatically accelerate our understanding of biological mechanisms and provide a solid foundation for reliable drug design. Although AF2 breaks through the barriers in predicting protein structures, many rooms remain to be further studied. This review provides a brief historical overview of the development of protein structure prediction, covering template-based, template-free, and machine learning-based methods. In addition to reviewing the potential benefits (Pros) and considerations (Cons) of using AF2, this review summarizes the diverse applications, including protein structure predictions, dynamic changes, point mutation, integration of language model and experimental data, protein complex, and protein-peptide interaction. It underscores recent advancements in efficiency, reliability, and broad application of AF2. This comprehensive review offers valuable insights into the applications of AF2 and AF2-inspired AI methods in structural biology and its potential for clinically significant drug target discovery.

Revolutionizing GPCR-ligand predictions: DeepGPCR with experimental validation for high-precision drug discovery

G-protein coupled receptors (GPCRs), crucial in various diseases, are targeted of over 40% of approved drugs. However, the reliable acquisition of experimental GPCRs structures is hindered by their lipid-embedded conformations. Traditional protein-ligand interaction models falter in GPCR-drug interactions, caused by limited and low-quality structures. Generalized models, trained on soluble protein-ligand pairs, are also inadequate. To address these issues, we developed two models, DeepGPCR_BC for binary classification and DeepGPCR_RG for affinity prediction. These models use non-structural GPCR-ligand interaction data, leveraging graph convolutional networks and mol2vec techniques to represent binding pockets and ligands as graphs. This approach significantly speeds up predictions while preserving critical physical-chemical and spatial information. In independent tests, DeepGPCR_BC surpassed Autodock Vina and Schrödinger Dock with an area under the curve of 0.72, accuracy of 0.68 and true positive rate of 0.73, whereas DeepGPCR_RG demonstrated a Pearson correlation of 0.39 and root mean squared error of 1.34. We applied these models to screen drug candidates for GPR35 (Q9HC97), yielding promising results with three (F545-1970, K297-0698, S948-0241) out of eight candidates. Furthermore, we also successfully obtained six active inhibitors for GLP-1R. Our GPCR-specific models pave the way for efficient and accurate large-scale virtual screening, potentially revolutionizing drug discovery in the GPCR field.

A Novel Peptidomimetic Insecticide: Dippu-AstR-Based Rational Design and Biological Activity of Allatostatin Analogs

Insect neuropeptides play an essential role in regulating growth, development, reproduction, nerve conduction, metabolism, and behavior in insects; therefore, G protein-coupled receptors of neuropeptides are considered important targets for designing green insecticides. Cockroach-type allatostatins (ASTs) (FGLamides allatostatins) are important insect neuropeptides in Diploptera punctata that inhibit juvenile hormone (JH) synthesis in the corpora allata and affect growth, development, and reproduction of insects. Therefore, the pursuit of novel insecticides targeting the allatostatin receptor (AstR) holds significant importance. Previously, we identified an AST analogue, H17, as a promising candidate for pest control. Herein, we first modeled the 3D structure of AstR in D. punctata (Dippu-AstR) and predicted the binding mode of H17 with Dippu-AstR to study the critical interactions and residues favorable to its bioactivity. Based on this binding mode, we designed and synthesized a series of H17 derivatives and assessed their insecticidal activity against D. punctata. Among them, compound Q6 showed higher insecticidal activity than H17 against D. punctata by inhibiting JH biosynthesis, indicating that Q6 is a potential candidate for a novel insect growth regulator (IGR)-based insecticide. Moreover, Q6 exhibited insecticidal activity against Plutella xylostella, indicating that these AST analogs may have a wider insecticidal spectrum. The underlying mechanisms and molecular conformations mediating the interactions of Q6 with Dippu-AstR were explored to understand its effects on the bioactivity. The present work clarifies how a target-based strategy facilitates the discovery of new peptide mimics with better bioactivity, enabling improved IGR-based insecticide potency in sustainable agriculture.

Accelerating Inference in Molecular Diffusion Models with Latent Representations of Protein Structure

Diffusion generative models have emerged as a powerful framework for addressing problems in structural biology and structure-based drug design. These models operate directly on 3D molecular structures. Due to the unfavorable scaling of graph neural networks (GNNs) with graph size as well as the relatively slow inference speeds inherent to diffusion models, many existing molecular diffusion models rely on coarse-grained representations of protein structure to make training and inference feasible. However, such coarse-grained representations discard essential information for modeling molecular interactions and impair the quality of generated structures. In this work, we present a novel GNN-based architecture for learning latent representations of molecular structure. When trained end-to-end with a diffusion model for de novo ligand design, our model achieves comparable performance to one with an all-atom protein representation while exhibiting a 3-fold reduction in inference time.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Identifying Residues for Substrate Recognition in Human GPAT4 by Molecular Dynamics Simulations

Glycerol-3-phosphate acyltransferase (GPAT) catalyzes the first step in triacylglycerol synthesis. Understanding its substrate recognition mechanism may help to design drugs to regulate the production of glycerol lipids in cells. In this work, we investigate how the native substrate, glycerol-3-phosphate (G3P), and palmitoyl-coenzyme A (CoA) bind to the human GPAT isoform GPAT4 via molecular dynamics simulations (MD). As no experimentally resolved GPAT4 structure is available, the AlphaFold model is employed to construct the GPAT4-substrate complex model. Using another isoform, GPAT1, we demonstrate that once the ligand binding is properly addressed, the AlphaFold complex model can deliver similar results to the experimentally resolved structure in MD simulations. Following the validated protocol of complex construction, we perform MD simulations using the GPAT4-substrate complex. Our simulations reveal that R427 is an important residue in recognizing G3P via a stable salt bridge, but its motion can bring the ligand to different binding hotspots on GPAT4. Such high flexibility can be attributed to the flexible region that exists only on GPAT4 and not on GPAT1. Our study reveals the substrate recognition mechanism of GPAT4 and hence paves the way towards designing GPAT4 inhibitors.

Molecular dynamics of the human RhD and RhAG blood group proteins

Introduction: Blood group antigens of the RH system (formerly known as "Rhesus") play an important role in transfusion medicine because of the severe haemolytic consequences of antibodies to these antigens. No crystal structure is available for RhD proteins with its partner RhAG, and the precise stoichiometry of the trimer complex remains unknown. Methods: To analyse their structural properties, the trimers formed by RhD and/or RhAG subunits were generated by protein modelling and molecular dynamics simulations were performed. Results: No major differences in structural behaviour were found between trimers of different compositions. The conformation of the subunits is relatively constant during molecular dynamics simulations, except for three large disordered loops. Discussion: This work makes it possible to propose a reasonable stoichiometry and demonstrates the potential of studying the structural behaviour of these proteins to investigate the hundreds of genetic variants relevant to transfusion medicine.

AlphaFold2 structures template ligand discovery

AlphaFold2 (AF2) and RosettaFold have greatly expanded the number of structures available for structure-based ligand discovery, even though retrospective studies have cast doubt on their direct usefulness for that goal. Here, we tested unrefined AF2 models prospectively, comparing experimental hit-rates and affinities from large library docking against AF2 models vs the same screens targeting experimental structures of the same receptors. In retrospective docking screens against the σ2 and the 5-HT2A receptors, the AF2 structures struggled to recapitulate ligands that we had previously found docking against the receptors' experimental structures, consistent with published results. Prospective large library docking against the AF2 models, however, yielded similar hit rates for both receptors versus docking against experimentally-derived structures; hundreds of molecules were prioritized and tested against each model and each structure of each receptor. The success of the AF2 models was achieved despite differences in orthosteric pocket residue conformations for both targets versus the experimental structures. Intriguingly, against the 5-HT2A receptor the most potent, subtype-selective agonists were discovered via docking against the AF2 model, not the experimental structure. To understand this from a molecular perspective, a cryoEM structure was determined for one of the more potent and selective ligands to emerge from docking against the AF2 model of the 5-HT2A receptor. Our findings suggest that AF2 models may sample conformations that are relevant for ligand discovery, much extending the domain of applicability of structure-based ligand discovery.

Drugging the entire human proteome: Are we there yet?

Each of the ∼20,000 proteins in the human proteome is a potential target for compounds that bind to it and modify its function. The 3D structures of most of these proteins are now available. Here, we discuss the prospects for using these structures to perform proteome-wide virtual HTS (VHTS). We compare physics-based (docking) and AI VHTS approaches, some of which are now being applied with large databases of compounds to thousands of targets. Although preliminary proteome-wide screens are now within our grasp, further methodological developments are expected to improve the accuracy of the results.

Exploring Novel Antidepressants Targeting G Protein-Coupled Receptors and Key Membrane Receptors Based on Molecular Structures

Major Depressive Disorder (MDD) is a complex mental disorder that involves alterations in signal transmission across multiple scales and structural abnormalities. The development of effective antidepressants (ADs) has been hindered by the dominance of monoamine hypothesis, resulting in slow progress. Traditional ADs have undesirable traits like delayed onset of action, limited efficacy, and severe side effects. Recently, two categories of fast-acting antidepressant compounds have surfaced, dissociative anesthetics S-ketamine and its metabolites, as well as psychedelics such as lysergic acid diethylamide (LSD). This has led to structural research and drug development of the receptors that they target. This review provides breakthroughs and achievements in the structure of depression-related receptors and novel ADs based on these. Cryo-electron microscopy (cryo-EM) has enabled researchers to identify the structures of membrane receptors, including the N-methyl-D-aspartate receptor (NMDAR) and the 5-hydroxytryptamine 2A (5-HT2A) receptor. These high-resolution structures can be used for the development of novel ADs using virtual drug screening (VDS). Moreover, the unique antidepressant effects of 5-HT1A receptors in various brain regions, and the pivotal roles of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) and tyrosine kinase receptor 2 (TrkB) in regulating synaptic plasticity, emphasize their potential as therapeutic targets. Using structural information, a series of highly selective ADs were designed based on the different role of receptors in MDD. These molecules have the favorable characteristics of rapid onset and low adverse drug reactions. This review offers researchers guidance and a methodological framework for the structure-based design of ADs.

Advances in structure-based drug design: The potential for precision therapeutics in psychiatric disorders

Over the years, the field of GPCR drug design has undergone a remarkable evolution, fueled by advancements in science and technology. This evolution has given rise to a diverse range of ideas and approaches in structure-based drug design, bolstering the versatility and strength of the GPCR drug design toolbox. This review encapsulates the iterative development process, navigating challenges and opportunities in structure-based drug design within GPCRs. With a focused emphasis on its impact on psychiatric disorders, the review accentuates recent advancements and delves into the potentials unlocked by emerging technologies. The review explores the intricate interplay between scientific progress and iterative refinement, offering profound insights into the potential pathways that lie ahead for GPCR drug design.

Computational drug development for membrane protein targets

The application of computational biology in drug development for membrane protein targets has experienced a boost from recent developments in deep learning-driven structure prediction, increased speed and resolution of structure elucidation, machine learning structure-based design and the evaluation of big data. Recent protein structure predictions based on machine learning tools have delivered surprisingly reliable results for water-soluble and membrane proteins but have limitations for development of drugs that target membrane proteins. Structural transitions of membrane proteins have a central role during transmembrane signaling and are often influenced by therapeutic compounds. Resolving the structural and functional basis of dynamic transmembrane signaling networks, especially within the native membrane or cellular environment, remains a central challenge for drug development. Tackling this challenge will require an interplay between experimental and computational tools, such as super-resolution optical microscopy for quantification of the molecular interactions of cellular signaling networks and their modulation by potential drugs, cryo-electron microscopy for determination of the structural transitions of proteins in native cell membranes and entire cells, and computational tools for data analysis and prediction of the structure and function of cellular signaling networks, as well as generation of promising drug candidates.

Comparative Structure-Based Virtual Screening Utilizing Optimized AlphaFold Model Identifies Selective HDAC11 Inhibitor

HDAC11 is a class IV histone deacylase with no crystal structure reported so far. The catalytic domain of HDAC11 shares low sequence identity with other HDAC isoforms, which makes conventional homology modeling less reliable. AlphaFold is a machine learning approach that can predict the 3D structure of proteins with high accuracy even in absence of similar structures. However, the fact that AlphaFold models are predicted in the absence of small molecules and ions/cofactors complicates their utilization for drug design. Previously, we optimized an HDAC11 AlphaFold model by adding the catalytic zinc ion and minimization in the presence of reported HDAC11 inhibitors. In the current study, we implement a comparative structure-based virtual screening approach utilizing the previously optimized HDAC11 AlphaFold model to identify novel and selective HDAC11 inhibitors. The stepwise virtual screening approach was successful in identifying a hit that was subsequently tested using an in vitro enzymatic assay. The hit compound showed an IC50 value of 3.5 µM for HDAC11 and could selectively inhibit HDAC11 over other HDAC subtypes at 10 µM concentration. In addition, we carried out molecular dynamics simulations to further confirm the binding hypothesis obtained by the docking study. These results reinforce the previously presented AlphaFold optimization approach and confirm the applicability of AlphaFold models in the search for novel inhibitors for drug discovery.

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24β1, -A24β2, -A34α2, -A34β1, and -A34β2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded β-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34β1/β2 had an additional two-stranded β-sheet. We hypothesize that this second β-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

Key aspects of modern GPCR drug discovery

G-protein-coupled receptors (GPCRs) are the largest and most versatile cell surface receptor family with a broad repertoire of ligands and functions. We've learned an enormous amount about discovering drugs of this receptor class since the first GPCR was cloned and expressed in 1986, such that it's now well-recognized that GPCRs are the most successful target class for approved drugs. Here we take the reader through a GPCR drug discovery journey from target to the clinic, highlighting the key learnings, best practices, challenges, trends and insights on discovering drugs that ultimately modulate GPCR function therapeutically in patients. The future of GPCR drug discovery is inspiring, with more desirable drug mechanisms and new technologies enabling the delivery of better and more successful drugs.

How accurately can one predict drug binding modes using AlphaFold models?

Computational prediction of protein structure has been pursued intensely for decades, motivated largely by the goal of using structural models for drug discovery. Recently developed machine-learning methods such as AlphaFold 2 (AF2) have dramatically improved protein structure prediction, with reported accuracy approaching that of experimentally determined structures. To what extent do these advances translate to an ability to predict more accurately how drugs and drug candidates bind to their target proteins? Here, we carefully examine the utility of AF2 protein structure models for predicting binding poses of drug-like molecules at the largest class of drug targets, the G-protein-coupled receptors. We find that AF2 models capture binding pocket structures much more accurately than traditional homology models, with errors nearly as small as differences between structures of the same protein determined experimentally with different ligands bound. Strikingly, however, the accuracy of ligand-binding poses predicted by computational docking to AF2 models is not significantly higher than when docking to traditional homology models and is much lower than when docking to structures determined experimentally without these ligands bound. These results have important implications for all those who might use predicted protein structures for drug discovery.

Artificial Intelligence-based database for prediction of protein structure and their alterations in ocular diseases

The aim of the study is to establish an online database for predicting protein structures altered in ocular diseases by Alphafold2 and RoseTTAFold algorithms. Totally, 726 genes of multiple ocular diseases were collected for protein structure prediction. Both Alphafold2 and RoseTTAFold algorithms were built locally using the open-source codebases. A dataset with 48 protein structures from Protein Data Bank (PDB) was adopted for algorithm set-up validation. A website was built to match ocular genes with the corresponding predicted tertiary protein structures for each amino acid sequence. The predicted local distance difference test-Cα (pLDDT) and template modeling (TM) scores of the validation protein structure and the selected ocular genes were evaluated. Molecular dynamics and molecular docking simulations were performed to demonstrate the applications of the predicted structures. For the validation dataset, 70.8% of the predicted protein structures showed pLDDT greater than 90. Compared to the PDB structures, 100% of the AlphaFold2-predicted structures and 97.9% of the RoseTTAFold-predicted structure showed TM score greater than 0.5. Totally, 1329 amino acid sequences of 430 ocular disease-related genes have been predicted, of which 75.9% showed pLDDT greater than 70 for the wildtype sequences and 76.1% for the variant sequences. Small molecule docking and molecular dynamics simulations revealed that the predicted protein structures with higher confidence scores showed similar molecular characteristics with the structures from PDB. We have developed an ocular protein structure database (EyeProdb) for ocular disease, which is released for the public and will facilitate the biological investigations and structure-based drug development for ocular diseases. Database URL:  http://eyeprodb.jsiec.org.

An overview of recent advances and challenges in predicting compound-protein interaction (CPI)

Compound-protein interactions (CPIs) are critical in drug discovery for identifying therapeutic targets, drug side effects, and repurposing existing drugs. Machine learning (ML) algorithms have emerged as powerful tools for CPI prediction, offering notable advantages in cost-effectiveness and efficiency. This review provides an overview of recent advances in both structure-based and non-structure-based CPI prediction ML models, highlighting their performance and achievements. It also offers insights into CPI prediction-related datasets and evaluation benchmarks. Lastly, the article presents a comprehensive assessment of the current landscape of CPI prediction, elucidating the challenges faced and outlining emerging trends to advance the field.

AlphaFold2-RAVE: From Sequence to Boltzmann Ranking

While AlphaFold2 is rapidly being adopted as a new standard in protein structure predictions, it is limited to single structures. This can be insufficient for the inherently dynamic world of biomolecules. In this Letter, we propose AlphaFold2-RAVE, an efficient protocol for obtaining Boltzmann-ranked ensembles from sequence. The method uses structural outputs from AlphaFold2 as initializations for artificial intelligence-augmented molecular dynamics. We release the method as an open-source code and demonstrate results on different proteins.

Promiscuity in drug discovery on the verge of the structural revolution: recent advances and future chances

INTRODUCTION

Promiscuity denotes the ability of ligands and targets to specifically interact with multiple binding partners. Despite negative aspects like side effects, promiscuity is receiving increasing attention in drug discovery as it can enhance drug efficacy and provides a molecular basis for drug repositioning. The three-dimensional structure of ligand-target complexes delivers exclusive insights into the molecular mechanisms of promiscuity and structure-based methods enable the identification of promiscuous interactions. With the recent breakthrough in protein structure prediction, novel possibilities open up to reveal unknown connections in ligand-target interaction networks.

AREAS COVERED

This review highlights the significance of structure in the identification and characterization of promiscuity and evaluates the potential of protein structure prediction to advance our knowledge of drug-target interaction networks. It discusses the definition and relevance of promiscuity in drug discovery and explores different approaches to detecting promiscuous ligands and targets.

EXPERT OPINION

Examination of structural data is essential for understanding and quantifying promiscuity. The recent advancements in structure prediction have resulted in an abundance of targets that are well-suited for structure-based methods like docking. In silico approaches may eventually completely transform our understanding of drug-target networks by complementing the millions of predicted protein structures with billions of predicted drug-target interactions.

Computational approaches streamlining drug discovery

Computer-aided drug discovery has been around for decades, although the past few years have seen a tectonic shift towards embracing computational technologies in both academia and pharma. This shift is largely defined by the flood of data on ligand properties and binding to therapeutic targets and their 3D structures, abundant computing capacities and the advent of on-demand virtual libraries of drug-like small molecules in their billions. Taking full advantage of these resources requires fast computational methods for effective ligand screening. This includes structure-based virtual screening of gigascale chemical spaces, further facilitated by fast iterative screening approaches. Highly synergistic are developments in deep learning predictions of ligand properties and target activities in lieu of receptor structure. Here we review recent advances in ligand discovery technologies, their potential for reshaping the whole process of drug discovery and development, as well as the challenges they encounter. We also discuss how the rapid identification of highly diverse, potent, target-selective and drug-like ligands to protein targets can democratize the drug discovery process, presenting new opportunities for the cost-effective development of safer and more effective small-molecule treatments.

AlphaFold2 and its applications in the fields of biology and medicine

AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.

Developing novel antifungals: lessons from G protein-coupled receptors

Up to 1.5 million people die yearly from fungal disease, but the repertoire of antifungal drug classes is minimal and the incidence of drug resistance is rising rapidly. This dilemma was recently declared by the World Health Organization as a global health emergency, but the discovery of new antifungal drug classes remains excruciatingly slow. This process could be accelerated by focusing on novel targets, such as G protein-coupled receptor (GPCR)-like proteins, that have a high likelihood of being druggable and have well-defined biology and roles in disease. We discuss recent successes in understanding the biology of virulence and in structure determination of yeast GPCRs, and highlight new approaches that might pay significant dividends in the urgent search for novel antifungal drugs.

Exploring AlphaFold2's Performance on Predicting Amino Acid Side-Chain Conformations and Its Utility in Crystal Structure Determination of B318L Protein

Recent technological breakthroughs in machine-learning-based AlphaFold2 (AF2) are pushing the prediction accuracy of protein structures to an unprecedented level that is on par with experimental structural quality. Despite its outstanding structural modeling capability, further experimental validations and performance assessments of AF2 predictions are still required, thus necessitating the development of integrative structural biology in synergy with both computational and experimental methods. Focusing on the B318L protein that plays an essential role in the African swine fever virus (ASFV) for viral replication, we experimentally demonstrate the high quality of the AF2 predicted model and its practical utility in crystal structural determination. Structural alignment implies that the AF2 model shares nearly the same atomic arrangement as the B318L crystal structure except for some flexible and disordered regions. More importantly, side-chain-based analysis at the individual residue level reveals that AF2's performance is likely dependent on the specific amino acid type and that hydrophobic residues tend to be more accurately predicted by AF2 than hydrophilic residues. Quantitative per-residue RMSD comparisons and further molecular replacement trials suggest that AF2 has a large potential to outperform other computational modeling methods in terms of structural determination. Additionally, it is numerically confirmed that the AF2 model is accurate enough so that it may well potentially withstand experimental data quality to a large extent for structural determination. Finally, an overall structural analysis and molecular docking simulation of the B318L protein are performed. Taken together, our study not only provides new insights into AF2's performance in predicting side-chain conformations but also sheds light upon the significance of AF2 in promoting crystal structural determination, especially when the experimental data quality of the protein crystal is poor.

Modeling the Orthosteric Binding Site of the G Protein-Coupled Odorant Receptor OR5K1

With approximately 400 encoding genes in humans, odorant receptors (ORs) are the largest subfamily of class A G protein-coupled receptors (GPCRs). Despite its high relevance and representation, the odorant-GPCRome is structurally poorly characterized: no experimental structures are available, and the low sequence identity of ORs to experimentally solved GPCRs is a significant challenge for their modeling. Moreover, the receptive range of most ORs is unknown. The odorant receptor OR5K1 was recently and comprehensively characterized in terms of cognate agonists. Here, we report two additional agonists and functional data of the most potent compound on two mutants, L104^3.32 and L255^6.51. Experimental data was used to guide the investigation of the binding modes of OR5K1 ligands into the orthosteric binding site using structural information from AI-driven modeling, as recently released in the AlphaFold Protein Structure Database, and from homology modeling. Induced-fit docking simulations were used to sample the binding site conformational space for ensemble docking. Mutagenesis data guided side chain residue sampling and model selection. We obtained models that could better rationalize the different activity of active (agonist) versus inactive molecules with respect to starting models and also capture differences in activity related to minor structural differences. Therefore, we provide a model refinement protocol that can be applied to model the orthosteric binding site of ORs as well as that of GPCRs with low sequence identity to available templates.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Evaluating GPCR modeling and docking strategies in the era of deep learning-based protein structure prediction

While deep learning (DL) has brought a revolution in the protein structure prediction field, still an important question remains how the revolution can be transferred to advances in structure-based drug discovery. Because the lessons from the recent GPCR dock challenge were inconclusive primarily due to the size of the dataset, in this work we further elaborated on 70 diverse GPCR complexes bound to either small molecules or peptides to investigate the best-practice modeling and docking strategies for GPCR drug discovery. From our quantitative analysis, it is shown that substantial improvements in docking and virtual screening have been possible by the advance in DL-based protein structure predictions with respect to the expected results from the combination of best pre-DL tools. The success rate of docking on DL-based model structures approaches that of cross-docking on experimental structures, showing over 30% improvement from the best pre-DL protocols. This amount of performance could be achieved only when two modeling points were considered properly: 1) correct functional-state modeling of receptors and 2) receptor-flexible docking. Best-practice modeling strategies and the model confidence estimation metric suggested in this work may serve as a guideline for future computer-aided GPCR drug discovery scenarios.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.