pKa Calculations of GPCRs: Understanding Protonation States in Receptor Activation.

The increase in the available G protein-coupled receptor (GPCR) structures has been pivotal in helping to understand their activation process. However, the role of protonation-conformation coupling in GPCR activation still needs to be clarified. We studied the protonation behavior of the highly conserved Asp2.50 residue in five different class A GPCRs (active and inactive conformations) using a linear response approximation (LRA) pKa calculation protocol. We observed consistent differences (1.3 pK units) for the macroscopic pKa values between the inactive and active states of the A2AR and B2AR receptors, indicating the protonation of Asp2.50 during GPCR activation. This process seems to be specific and not conserved, as no differences were observed in the pKa values of the remaining receptors (CB1R, NT1R, and GHSR).

Current genetic engineering strategies for the production of antihypertensive ACEI peptides.

Hypertension is a major and highly prevalent risk factor for various diseases. Among the most frequently prescribed antihypertensive first-line drugs are synthetic angiotensin I-converting enzyme inhibitors (ACEI). However, since  their use in hypertension therapy has been linked to various side effects, interest in the application of food-derived ACEI peptides (ACEIp) as antihypertensive agents is rapidly growing. Although promising, the industrial production of ACEIp through conventional methods such as chemical synthesis or enzymatic hydrolysis of food proteins has been proven troublesome. We here provide an overview of current antihypertensive therapeutics, focusing on ACEI, and illustrate how biotechnology and bioengineering can overcome the limitations of ACEIp large-scale production. Latest advances in ACEIp research and current genetic engineering-based strategies for heterologous production of ACEIp (and precursors) are also presented. Cloning approaches include tandem repeats of single ACEIp, ACEIp fusion to proteins/polypeptides, joining multivariate ACEIp into bioactive polypeptides, and producing ACEIp-containing modified plant storage proteins. Although bacteria have been privileged ACEIp heterologous hosts, particularly when testing for new genetic engineering strategies, plants and microalgae-based platforms are now emerging. Besides being generally safer, cost-effective and scalable, these "pharming" platforms can perform therelevant posttranslational modifications and produce (and eventually deliver) biologically active protein/peptide-based antihypertensive medicines.

Understanding the Binding Specificity of G-Protein Coupled Receptors toward G-Proteins and Arrestins: Application to the Dopamine Receptor Family.

G-Protein coupled receptors (GPCRs) are involved in a myriad of pathways key for human physiology through the formation of complexes with intracellular partners such as G-proteins and arrestins (Arrs). However, the structural and dynamical determinants of these complexes are still largely unknown. Herein, we developed a computational big-data pipeline that enables the structural characterization of GPCR complexes with no available structure. This pipeline was used to study a well-known group of catecholamine receptors, the human dopamine receptor (DXR) family and its complexes, producing novel insights into the physiological properties of these important drug targets. A detailed description of the protein interfaces of all members of the DXR family (D1R, D2R, D3R, D4R, and D5R) and the corresponding protein interfaces of their binding partners (Arrs: Arr2 and Arr3; G-proteins: Gi1, Gi2, Gi3, Go, Gob, Gq, Gslo, Gssh, Gt2, and Gz) was generated. To produce reliable structures of the DXR family in complex with either G-proteins or Arrs, we performed homology modeling using as templates the structures of the ß2-adrenergic receptor (ß2AR) bound to Gs, the rhodopsin bound to Gi, and the recently acquired neurotensin receptor-1 (NTSR1) and muscarinic 2 receptor (M2R) bound to arrestin (Arr). Among others, the work demonstrated that the three partner groups, Arrs and Gs- and Gi-proteins, are all structurally and dynamically distinct. Additionally, it was revealed the involvement of different structural motifs in G-protein selective coupling between D1- and D2-like receptors. Having constructed and analyzed 50 models involving DXR, this work represents an unprecedented large-scale analysis of GPCR-intracellular partner interface determinants. All data is available at www.moreiralab.com/resources/dxr.

A Complete Assessment of Dopamine Receptor- Ligand Interactions through Computational Methods.

Background: Selectively targeting dopamine receptors (DRs) has been a persistent challenge in the last years for the development of new treatments to combat the large variety of diseases involving these receptors. Although, several drugs have been successfully brought to market, the subtype-specific binding mode on a molecular basis has not been fully elucidated. Methods: Homology modeling and molecular dynamics were applied to construct robust conformational models of all dopamine receptor subtypes (D1-like and D2-like). Fifteen structurally diverse ligands were docked. Contacts at the binding pocket were fully described in order to reveal new structural findings responsible for selective binding to DR subtypes. Results: Residues of the aromatic microdomain were shown to be responsible for the majority of ligand interactions established to all DRs. Hydrophobic contacts involved a huge network of conserved and non-conserved residues between three transmembrane domains (TMs), TM2-TM3-TM7. Hydrogen bonds were mostly mediated by the serine microdomain. TM1 and TM2 residues were main contributors for the coupling of large ligands. Some amino acid groups form electrostatic interactions of particular importance for D1R-like selective ligands binding. Conclusions: This in silico approach was successful in showing known receptor-ligand interactions as well as in determining unique combinations of interactions, which will support mutagenesis studies to improve the design of subtype-specific ligands.

MUG: A mutation overview of GPCR subfamily A17 receptors.

G protein-coupled receptors (GPCRs) mediate several signaling pathways through a general mechanism that involves their activation, upholding a chain of events that lead to the release of molecules responsible for cytoplasmic action and further regulation. These physiological functions can be severely altered by mutations in GPCR genes. GPCRs subfamily A17 (dopamine, serotonin, adrenergic and trace amine receptors) are directly related with neurodegenerative diseases, and as such it is crucial to explore known mutations on these systems and their impact in structure and function. A comprehensive and detailed computational framework - MUG (Mutations Understanding GPCRs) - was constructed, illustrating key reported mutations and their effect on receptors of the subfamily A17 of GPCRs. We explored the type of mutations occurring overall and in the different families of subfamily A17, as well their localization within the receptor and potential effects on receptor functionality. The mutated residues were further analyzed considering their pathogenicity. The results reveal a high diversity of mutations in the GPCR subfamily A17 structures, drawing attention to the considerable number of mutations in conserved residues and domains. Mutated residues were typically hydrophobic residues enriched at the ligand binding pocket and known activating microdomains, which may lead to disruption of receptor function. MUG as an interactive web application is available for the management and visualization of this dataset. We expect that this interactive database helps the exploration of GPCR mutations, their influence, and their familywise and receptor-specific effects, constituting the first step in elucidating their structures and molecules at the atomic level.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements.

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the ß-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or ß-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

Decoding Partner Specificity of Opioid Receptor Family.

This paper describes an exciting big data analysis compiled in a freely available database, which can be applied to characterize the coupling of different G-Protein coupled receptors (GPCRs) families with their intracellular partners. Opioid receptor (OR) family was used as case study in order to gain further insights into the physiological properties of these important drug targets, known to be associated with the opioid crisis, a huge socio-economic issue directly related to drug abuse. An extensive characterization of all members of the ORs family (µ (MOR), δ (DOR), κ (KOR), nociceptin (NOP)) and their corresponding binding partners (ARRs: Arr2, Arr3; G-protein: Gi1, Gi2, Gi3, Go, Gob, Gz, Gq, G11, G14, G15, G12, Gssh, Gslo) was carried out. A multi-step approach including models' construction (multiple sequence alignment, homology modeling), complex assembling (protein complex refinement with HADDOCK and complex equilibration), and protein-protein interface characterization (including both structural and dynamics analysis) were performed. Our database can be easily applied to several GPCR sub-families, to determine the key structural and dynamical determinants involved in GPCR coupling selectivity.

In Silico End-to-End Protein-Ligand Interaction Characterization Pipeline: The Case of SARS-CoV-2.

SARS-CoV-2 triggered a worldwide pandemic disease, COVID-19, for which an effective treatment has not yet been settled. Among the most promising targets to fight this disease is SARS-CoV-2 main protease (Mpro), which has been extensively studied in the last few months. There is an urgency for developing effective computational protocols that can help us tackle these key viral proteins. Hence, we have put together a robust and thorough pipeline of in silico protein-ligand characterization methods to address one of the biggest biological problems currently plaguing our world. These methodologies were used to characterize the interaction of SARS-CoV-2 Mpro with an α-ketoamide inhibitor and include details on how to upload, visualize, and manage the three-dimensional structure of the complex and acquire high-quality figures for scientific publications using PyMOL (Protocol 1); perform homology modeling with MODELLER (Protocol 2); perform protein-ligand docking calculations using HADDOCK (Protocol 3); run a virtual screening protocol of a small compound database of SARS-CoV-2 candidate inhibitors with AutoDock 4 and AutoDock Vina (Protocol 4); and, finally, sample the conformational space at the atomic level between SARS-CoV-2 Mpro and the α-ketoamide inhibitor with Molecular Dynamics simulations using GROMACS (Protocol 5). Guidelines for careful data analysis and interpretation are also provided for each Protocol.

Prediction and targeting of GPCR oligomer interfaces.

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.