Search

Homeostasis of the host immune system is regulated by white blood cells with a variety of cell surface receptors for cytokines. Chemotactic cytokines (chemokines) activate their receptors to evoke the chemotaxis of immune cells in homeostatic migrations or inflammatory conditions towards inflamed tissue or pathogens. Dysregulation of the immune system leading to disorders such as allergies, autoimmune diseases, or cancer requires efficient, fast-acting drugs to minimize the long-term effects of chronic inflammation. Here, we performed structure-based virtual screening (SBVS) assisted by the Keras/TensorFlow neural network (NN) to find novel compound scaffolds acting on three chemokine receptors: CCR2, CCR3, and one CXC receptor, CXCR3. Keras/TensorFlow NN was used here not as a typically used binary classifier but as an efficient multi-class classifier that can discard not only inactive compounds but also low- or medium-activity compounds. Several compounds proposed by SBVS and NN were tested in 100 ns all-atom molecular dynamics simulations to confirm their binding affinity. To improve the basic binding affinity of the compounds, new chemical modifications were proposed. The modified compounds were compared with known antagonists of these three chemokine receptors. Known CXCR3 compounds were among the top predicted compounds; thus, the benefits of using Keras/TensorFlow in drug discovery have been shown in addition to structure-based approaches. Furthermore, we showed that Keras/TensorFlow NN can accurately predict the receptor subtype selectivity of compounds, for which SBVS often fails. We cross-tested chemokine receptor datasets retrieved from ChEMBL and curated datasets for cannabinoid receptors. The NN model trained on the cannabinoid receptor datasets retrieved from ChEMBL was the most accurate in the receptor subtype selectivity prediction. Among NN models trained on the chemokine receptor datasets, the CXCR3 model showed the highest accuracy in differentiating the receptor subtype for a given compound dataset.

Chemokines modulate the immune response by regulating the migration of immune cells. They are also known to participate in such processes as cell-cell adhesion, allograft rejection, and angiogenesis. Chemokines interact with two different subfamilies of G protein-coupled receptors: conventional chemokine receptors and atypical chemokine receptors. Here, we focused on the former one which has been linked to many inflammatory diseases, including: multiple sclerosis, asthma, nephritis, and rheumatoid arthritis. Available crystal and cryo-EM structures and homology models of six chemokine receptors (CCR1 to CCR6) were described and tested in terms of their usefulness in structure-based drug design. As a result of structure-based virtual screening for CCR2 and CCR3, several new active compounds were proposed. Known inhibitors of CCR1 to CCR6, acquired from ChEMBL, were used as training sets for two machine learning algorithms in ligand-based drug design. Performance of LightGBM was compared with a sequential Keras/TensorFlow model of neural network for these diverse datasets. A combination of structure-based virtual screening with machine learning allowed to propose several active ligands for CCR2 and CCR3 with two distinct compounds predicted as CCR3 actives by all three tested methods: Glide, Keras/TensorFlow NN, and LightGBM. In addition, the performance of these three methods in the prediction of the CCR2/CCR3 receptor subtype selectivity was assessed.

Bacteriophages effectively counteract diverse bacterial infections, and their ability to treat most types of cancer has been explored using phage engineering or phage-virus hybrid platforms. In the present study, it was demonstrated that the bacteriophage MS2 can affect the expression of genes associated with the proliferation and survival of LNCaP prostate epithelial cells. LNCaP cells were exposed to bacteriophage MS2 at a concentration of 1×10 7 plaque forming units/ml for 24-48 h. After exposure, various cellular parameters, including cell viability, morphology, and changes in gene expression, were examined. MS2 affected cell viability adversely, reducing viability by 25% in the first 4 h of treatment; however, cell viability recovered within 24-48 h. Similarly, the AKT , androgen receptor, integrin α5, integrin β1, MAPK1, MAPK3, STAT3 , and peroxisome proliferator-activated receptor-γ coactivator 1α genes, which are involved in various normal cellular processes and tumor progression, were significantly upregulated, whereas the expression levels of HSP90, ITGB5, ITGB3, HSP27, ITGAV , and PI3K genes were unchanged. Therefore, based on viability and gene expression changes, bacteriophage MS2 severely impaired LNCaP cells by reducing anchorage-dependent survival and androgen signaling. A caveolin-mediated endocytosis mechanism for MS2-mediated signaling in prostate cancer cells was proposed based on reports involving bacteriophages T4, M13, and MS2, and their interactions with LNCaP and PC3 cell lines., Competing Interests: The authors declare that they have no competing interests., (Copyright: © Sanmukh et al.)

Host response to infection involves the activation of the complement system leading to the production of anaphylatoxins C3a and C5a. Complement factor C5a exerts its effect through the activation of C5aR1, chemotactic receptor 1, and triggers the G protein-coupled signaling cascade. Orthosteric and allosteric antagonists of C5aR1 are a novel strategy for anti-inflammatory therapies. Here, we discuss recent crystal structures of inactive C5aR1 in terms of an inverted orientation of helix H8, unobserved in other GPCR structures. An analysis of mutual interactions of subunits in the C5aR1-G protein complex has provided new insights into the activation mechanism of this distinct receptor. By comparing two C5aR receptors C5aR1 and C5aR2 we explained differences between their signaling pathways on the molecular level. By means of molecular dynamics we explained why C5aR2 cannot transduce signal through the G protein pathway but instead recruits beta-arrestin. A comparison of microsecond MD trajectories started from active and inactive C5aR1 receptor conformations has provided insights into details of local and global changes in the transmembrane domain induced by interactions with the Gα subunit and explained the impact of inverted H8 on the C5aR1 activation., Competing Interests: The authors have declared that no competing interests exist.

Homeostasis of the human immune system is regulated by many cellular components, including two neuropeptides, VIP and PACAP, primary stimuli for three class B G protein-coupled receptors, VPAC1, VPAC2, and PAC1. Vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating polypeptide (PACAP) regulate intestinal motility and secretion and influence the functioning of the endocrine and immune systems. Inhibition of VIP and PACAP receptors is an emerging concept for new pharmacotherapies for chronic inflammation and cancer, while activation of their receptors provides neuroprotection. A small number of known active compounds for these receptors still impose limitations on their use in therapeutics. Recent cryo-EM structures of VPAC1 and PAC1 receptors in their agonist-bound active state have provided insights regarding their mechanism of activation. Here, we describe major molecular switches of VPAC1, VPAC2, and PAC1 that may act as triggers for receptor activation and compare them with similar non-covalent interactions changing upon activation that were observed for other GPCRs. Interhelical interactions in VIP and PACAP receptors that are important for agonist binding and/or activation provide a molecular basis for the design of novel selective drugs demonstrating anti-inflammatory, anti-cancer, and neuroprotective effects. The impact of genetic variants of VIP, PACAP, and their receptors on signalling mediated by endogenous agonists is also described. This sequence diversity resulting from gene splicing has a significant impact on agonist selectivity and potency as well as on the signalling properties of VIP and PACAP receptors.

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating polypeptide (PACAP) are two neuropeptides that contribute to the regulation of intestinal motility and secretion, exocrine and endocrine secretions, and homeostasis of the immune system. Their biological effects are mediated by three receptors named VPAC1, VPAC2 and PAC1 that belong to class B GPCRs. VIP and PACAP receptors have been identified as potential therapeutic targets for the treatment of chronic inflammation, neurodegenerative diseases and cancer. However, pharmacological use of endogenous ligands for these receptors is limited by their lack of specificity (PACAP binds with high affinity to VPAC1, VPAC2 and PAC1 receptors while VIP recognizes both VPAC1 and VPAC2 receptors), their poor oral bioavailability (VIP and PACAP are 27- to 38-amino acid peptides) and their short half-life. Therefore, the development of non-peptidic small molecules or specific stabilized peptidic ligands is of high interest. Structural similarities between VIP and PACAP receptors are major causes of difficulties in the design of efficient and selective compounds that could be used as therapeutics. In this study we performed structure-based virtual screening against the subset of the ZINC15 drug library. This drug repositioning screen provided new applications for a known drug: ticagrelor, a P2Y12 purinergic receptor antagonist. Ticagrelor inhibits both VPAC1 and VPAC2 receptors which was confirmed in VIP-binding and calcium mobilization assays. A following analysis of detailed ticagrelor binding modes to all three VIP and PACAP receptors with molecular dynamics revealed its allosteric mechanism of action. Using a validated homology model of inactive VPAC1 and a recently released cryo-EM structure of active VPAC1 we described how ticagrelor could block conformational changes in the region of 'tyrosine toggle switch' required for the receptor activation. We also discuss possible modifications of ticagrelor comparing other P2Y12 antagonist - cangrelor, closely related to ticagrelor but not active for VPAC1/VPAC2. This comparison with inactive cangrelor could lead to further improvement of the ticagrelor activity and selectivity for VIP and PACAP receptor sub-types., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Langer and Latek.)

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q 2 > 0.63 and Q 2 > 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

The selective targeting of the cannabinoid receptor 2 (CB2) is crucial for the development of peripheral system-acting cannabinoid analgesics. This work aimed at computer-assisted identification of prospective CB2-selective compounds among the constituents of Cannabis Sativa . The molecular structures and corresponding binding affinities to CB1 and CB2 receptors were collected from ChEMBL. The molecular structures of Cannabis Sativa constituents were collected from a phytochemical database. The collected records were curated and applied for the development of quantitative structure-activity relationship (QSAR) models with a machine learning approach. The validated models predicted the affinities of Cannabis Sativa constituents. Four structures of CB2 were acquired from the Protein Data Bank (PDB) and the discriminatory ability of CB2-selective ligands and two sets of decoys were tested. We succeeded in developing the QSAR model by achieving Q 2 5-CV > 0.62. The QSAR models helped to identify three prospective CB2-selective molecules that are dissimilar to already tested compounds. In a complementary structure-based virtual screening study that used available PDB structures of CB2, the agonist-bound, Cryogenic Electron Microscopy structure of CB2 showed the best statistical performance in discriminating between CB2-active and non-active ligands. The same structure also performed best in discriminating between CB2-selective ligands from non-selective ligands.

Vasoactive intestinal peptide receptor 1 (VPAC1) is a member of a secretin-like subfamily of G protein-coupled receptors. Its endogenous neuropeptide (VIP), secreted by neurons and immune cells, modulates various physiological functions such as exocrine and endocrine secretions, immune response, smooth muscles relaxation, vasodilation, and fetal development. As a drug target, VPAC1 has been selected for therapy of inflammatory diseases but drug discovery is still hampered by lack of its crystal structure. In this study we presented the homology model of this receptor constructed with the well-known web service GPCRM. The VPAC1 model is composed of extracellular and transmembrane domains that form a complex with an endogenous hormone VIP. Using the homology model of VPAC1 the mechanism of action of potential drug candidates for VPAC1 was described. Only two series of small-molecule antagonists of confirmed biological activity for VPAC1 have been described thus far. Molecular docking and a series of molecular dynamics simulations were performed to elucidate their binding to VPAC1 and resulting antagonist effect. The presented work provides the basis for the possible binding mode of VPAC1 antagonists and determinants of their molecular recognition in the context of other class B GPCRs. Until the crystal structure of VPAC1 will be released, the presented homology model of VPAC1 can serve as a scaffold for drug discovery studies and is available from the author upon request.

The prolonged use of many currently available drugs results in the severe side effect of the disruption of glucose metabolism leading to type 2 diabetes mellitus (T2DM. Gut hormone receptors including glucagon receptor (GCGR) and the incretin hormone receptors: glucagon-like peptide 1 receptor (GLP1R) and gastric inhibitory polypeptide receptor (GIPR) are important drug targets for the treatment of T2DM, as they play roles in the regulation of glucose and insulin levels and of food intake. In this study, we hypothesized that we could compensate for the negative influences of specific drugs on glucose metabolism by the positive incretin effect enhanced by the off-target interactions with incretin GPCR receptors. As a test case, we chose to examine beta-blockers because beta-adrenergic receptors and incretin receptors are expressed in a similar location, making off-target interactions possible. The binding affinity of drugs for incretin receptors was approximated by using two docking scoring functions of Autodock VINA (GUT-DOCK) and Glide (Schrodinger) and juxtaposing these values with the medical information on drug-induced T2DM. We observed that beta-blockers with the highest theoretical binding affinities for gut hormone receptors were reported as the least harmful to glucose homeostasis in clinical trials. Notably, a recently discovered beta-blocker compound 15 ([4-((2S)-3-(((S)-3-(3-bromophenyl)-1-(methylamino)-1-oxopropan-2-yl)amino)-2-(2-cyclohexyl-2-phenylacetamido)-3-oxopropyl)benzamide was among the top-scoring drugs, potentially supporting its use in the treatment of hypertension in diabetic patients. Our recently developed web service GUT-DOCK (gut-dock.miningmembrane.com) allows for the execution of similar studies for any drug-like molecule. Specifically, users can compute the binding affinities for various class B GPCRs, gut hormone receptors, VIPR1 and PAC1R., Competing Interests: The authors have declared that no competing interests exist.

A disturbance of glucose homeostasis leading to type 2 diabetes mellitus (T2DM) is one of the severe side effects that may occur during a prolonged use of many drugs currently available on the market. In this manuscript we describe the most common cases of drug-induced T2DM, discuss available pharmacotherapies and propose new ones. Among various pharmacotherapies of T2DM, incretin therapies have recently focused attention due to the newly determined crystal structure of incretin hormone receptor GLP1R. Incretin hormone receptors: GLP1R and GIPR together with the glucagon receptor GCGR regulate food intake and insulin and glucose secretion. Our study showed that incretin hormone receptors, named also gut hormone receptors as they are expressed in the gastrointestinal tract, could potentially act as unintended targets (off-targets) for orally administrated drugs. Such off-target interactions, depending on their effect on the receptor (stimulation or inhibition), could be beneficial, like in the case of incretin mimetics, or unwanted if they cause, e.g., decreased insulin secretion. In this in silico study we examined which well-known pharmaceuticals could potentially interact with gut hormone receptors in the off-target way. We observed that drugs with the strongest binding affinity for gut hormone receptors were also reported in the medical information resources as the least disturbing the glucose homeostasis among all drugs in their class. We suggested that those strongly binding molecules could potentially stimulate GIPR and GLP1R and/or inhibit GCGR which could lead to increased insulin secretion and decreased hepatic glucose production. Such positive effect on the glucose homeostasis could compensate for other, adverse effects of pharmacotherapy which lead to drug-induced T2DM. In addition, we also described several top hits as potential substitutes of peptidic incretin mimetics which were discovered in the drug repositioning screen using gut hormone receptors structures against the ZINC15 compounds subset., Competing Interests: The authors have declared that no competing interests exist.

Due to the involvement of G protein-coupled receptors (GPCRs) in most of the physiological and pathological processes in humans they have been attracting a lot of attention from pharmaceutical industry as well as from scientific community. Therefore, the need for new, high quality structures of GPCRs is enormous. The updated homology modeling service GPCRM (http://gpcrm.biomodellab.eu/) meets those expectations by greatly reducing the execution time of submissions (from days to hours/minutes) with nearly the same average quality of obtained models. Additionally, due to three different scoring functions (Rosetta, Rosetta-MP, BCL::Score) it is possible to select accurate models for the required purposes: the structure of the binding site, the transmembrane domain or the overall shape of the receptor. Currently, no other web service for GPCR modeling provides this possibility. GPCRM is continually upgraded in a semi-automatic way and the number of template structures has increased from 20 in 2013 to over 90 including structures the same receptor with different ligands which can influence the structure not only in the on/off manner. Two types of protein viewers can be used for visual inspection of obtained models. The extended sortable tables with available templates provide links to external databases and display ligand-receptor interactions in visual form.

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

Background: Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family., Results: The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100 ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes., Conclusions: The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated.

The recent GPCR Dock 2013 assessment of serotonin receptor 5-HT1B and 5-HT2B, and smoothened receptor SMO targets, exposed the strengths and weaknesses of the currently used computational approaches. The test cases of 5-HT1B and 5-HT2B demonstrated that both the receptor structure and the ligand binding mode can be predicted with the atomic-detail accuracy, as long as the target-template sequence similarity is relatively high. On the other hand, the observation of a low target-template sequence similarity, e.g., between SMO from the frizzled GPCR family and members of the rhodopsin family, hampers the GPCR structure prediction and ligand docking. Indeed, in GPCR Dock 2013, accurate prediction of the SMO target was still beyond the capabilities of most research groups. Another bottleneck in the current GPCR research, as demonstrated by the 5-HT2B target, is the reliable prediction of global conformational changes induced by activation of GPCRs. In this work, we report details of our protocol used during GPCR Dock 2013. Our structure prediction and ligand docking protocol was especially successful in the case of 5-HT1B and 5-HT2B-ergotamine complexes for which we provide one of the most accurate predictions. In addition to a description of the GPCR Dock 2013 results, we propose a novel hybrid computational methodology to improve GPCR structure and function prediction. This computational methodology employs two separate rankings for filtering GPCR models. The first ranking is ligand-based while the second is based on the scoring scheme of the recently published BCL method. In this work, we prove that the use of knowledge-based potentials implemented in BCL is an efficient way to cope with major bottlenecks in the GPCR structure prediction. Thereby, we also demonstrate that the knowledge-based potentials for membrane proteins were significantly improved, because of the recent surge in available experimental structures.

A key component to success in structure-based drug design is reliable information on protein-ligand interactions. Recent development in NMR techniques has accelerated this process by overcoming some of the limitations of X-ray crystallography and computational protein-ligand docking. In this work we present a new scoring protocol based on NMR-derived interligand INPHARMA NOEs to guide the selection of computationally generated docking modes. We demonstrate the performance in a range of scenarios, encompassing traditionally difficult cases such as docking to homology models and ligand dependent domain rearrangements. Ambiguities associated with sparse experimental information are lifted by searching a consensus solution based on simultaneously fitting multiple ligand pairs. This study provides a previously unexplored integration between molecular modeling and experimental data, in which interligand NOEs represent the key element in the rescoring algorithm. The presented protocol should be widely applicable for protein-ligand docking also in a different context from drug design and highlights the important role of NMR-based approaches to describe intermolecular ligand-receptor interactions.

Unlabelled: G-protein coupled receptors (GPCRs) are targets of nearly one third of the drugs at the current pharmaceutical market. Despite their importance in many cellular processes the crystal structures are available for less than 20 unique GPCRs of the Rhodopsin-like class. Fortunately, even though involved in different signaling cascades, this large group of membrane proteins has preserved a uniform structure comprising seven transmembrane helices that allows quite reliable comparative modeling. Nevertheless, low sequence similarity between the GPCR family members is still a serious obstacle not only in template selection but also in providing theoretical models of acceptable quality. An additional level of difficulty is the prediction of kinks and bulges in transmembrane helices. Usage of multiple templates and generation of alignments based on sequence profiles may increase the rate of success in difficult cases of comparative modeling in which the sequence similarity between GPCRs is exceptionally low. Here, we present GPCRM, a novel method for fast and accurate generation of GPCR models using averaging of multiple template structures and profile-profile comparison. In particular, GPCRM is the first GPCR structure predictor incorporating two distinct loop modeling techniques: Modeller and Rosetta together with the filtering of models based on the Z-coordinate. We tested our approach on all unique GPCR structures determined to date and report its performance in comparison with other computational methods targeting the Rhodopsin-like class. We also provide a database of precomputed GPCR models of the human receptors from that class., Availability: GPCRM SERVER AND DATABASE: http://gpcrm.biomodellab.eu.

Sphingosine 1-phosphate (S1P) is a lysophospholipid mediator which activates G protein-coupled sphingosine 1-phosphate receptors and thus evokes a variety of cell and tissue responses including lymphocyte trafficking, endothelial development, integrity, and maturation. We performed five all-atom 700 ns molecular dynamics simulations of the sphingosine 1-phosphate receptor 1 (S1P₁) based on recently released crystal structure of that receptor with an antagonist. We found that the initial movements of amino acid residues occurred in the area of highly conserved W269⁶·⁴⁸ in TM6 which is close to the ligand binding location. Those residues located in the central part of the receptor and adjacent to kinks of TM helices comprise of a transmission switch. Side chains movements of those residues were coupled to the movements of water molecules inside the receptor which helped in the gradual opening of intracellular part of the receptor. The most stable parts of the protein were helices TM1 and TM2, while the largest movement was observed for TM7, possibly due to the short intracellular part starting with a helix kink at P⁷·⁵⁰, which might be the first helix to move at the intracellular side. We show for the first time the detailed view of the concerted action of the transmission switch and Trp (W⁶·⁴⁸) rotamer toggle switch leading to redirection of water molecules flow in the central part of the receptor. That event is a prerequisite for subsequent changes in intracellular part of the receptor involving water influx and opening of the receptor structure.

Vacuolar ATPases are a potential therapeutic target because of their involvement in a variety of severe diseases such as osteoporosis or cancer. Archazolide A (1) and related analogs have been previously identified as selective inhibitors of V-ATPases with potency down to the subnanomolar range. Herein we report on the determination of the ligand binding mode by a combination of molecular docking, molecular dynamics simulations, and biochemical experiments, resulting in a sound model for the inhibitory mechanism of this class of putative anticancer agents. The binding site of archazolides was confirmed to be located in the equatorial region of the membrane-embedded V(O)-rotor, as recently proposed on the basis of site-directed mutagenesis. Quantification of the bioactivity of a series of archazolide derivatives, together with the docking-derived binding mode of archazolides to the V-ATPase, revealed favorable ligand profiles, which can guide the development of a simplified archazolide analog with potential therapeutic relevance.

The Formyl Peptide Receptor 1 (FPR1) is an important chemotaxis receptor involved in various aspects of host defense and inflammatory processes. We constructed a model of FPR1 using as a novel template the chemokine receptor CXCR4 from the same branch of the phylogenetic tree of G-protein-coupled receptors. The previously employed template of rhodopsin contained a bulge at the extracellular part of TM2 which directly influenced binding of ligands. We also conducted molecular dynamics (MD) simulations of FPR1 in the apo form as well as in a form complexed with the agonist fMLF and the antagonist tBocMLF in the model membrane. During all MD simulation of the fMLF-FPR1 complex a water molecule transiently bridged the hydrogen bond between W254(6.48) and N108(3.35) in the middle of the receptor. We also observed a change in the cytoplasmic part of FPR1 of a rotamer of the Y301(7.53) residue (tyrosine rotamer switch). This effect facilitated movement of more water molecules toward the receptor center. Such rotamer of Y301(7.53) was not observed in any crystal structures of GPCRs which can suggest that this state is temporarily formed to pass the water molecules during the activation process. The presence of a distance between agonist and residues R201(5.38) and R205(5.42) on helix TM5 may suggest that the activation of FPR1 is similar to the activation of β-adrenergic receptors since their agonists are separated from serine residues on helix TM5. The removal of water molecules bridging these interactions in FPR1 can result in shrinking of the binding site during activation similarly to the shrinking observed in β-ARs. The number of GPCR crystal structures with agonists is still scarce so the designing of new ligands with agonistic properties is hampered, therefore homology modeling and docking can provide suitable models. Additionally, the MD simulations can be beneficial to outline the mechanisms of receptor activation and the agonist/antagonist sensing.

The years 2000 and 2007 witnessed milestones in current understanding of G protein-coupled receptor (GPCR) structural biology. In 2000 the first GPCR, bovine rhodopsin, was crystallized and the structure was solved, while in 2007 the structure of β(2)-adrenergic receptor, the first GPCR with diffusible ligands, was determined owing to advances in microcrystallization and an insertion of the fast-folding lysozyme into the receptor. In parallel with those crystallographic studies, the biological and biochemical characterization of GPCRs has advanced considerably because those receptors are molecular targets for many of currently used drugs. Therefore, the mechanisms of activation and signal transduction to the cell interior deduced from known GPCRs structures are of the highest importance for drug discovery. These proteins are the most diversified membrane receptors encoded by hundreds of genes in our genome. They participate in processes responsible for vision, smell, taste and neuronal transmission in response to photons or binding of ions, hormones, peptides, chemokines and other factors. Although the GPCRs share a common seven-transmembrane α-helical bundle structure their binding sites can accommodate thousands of different ligands. The ligands, including agonists, antagonists or inverse agonists change the structure of the receptor. With bound agonists they can form a complex with a suitable G protein, be phosphorylated by kinases or bind arrestin. The discovered signaling cascades invoked by arrestin independently of G proteins makes the GPCR activating scheme more complex such that a ligand acting as an antagonist for G protein signaling can also act as an agonist in arrestin-dependent signaling. Additionally, the existence of multiple ligand-dependent partial activation states as well as dimerization of GPCRs result in a 'microprocessor-like' action of these receptors rather than an 'on-off' switch as was commonly believed only a decade ago.

Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) G protein coupled receptors. Docking of agonists and antagonists to CB(1) and CB(2) cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch and its possible participation in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB(1) and CB(2) receptor models were constructed based on the adenosine A(2A) receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β(2)AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.

Recent development of nuclear magnetic resonance (NMR) techniques provided new types of structural restraints that can be successfully used in fast and low-cost global protein fold determination. Here, we present CABS-NMR, an efficient protein modeling tool, which takes advantage of such structural restraints. The restraints are converted from original NMR data to fit the coarse grained protein representation of the C-Alpha-Beta-Side-group (CABS) algorithm. CABS is a Monte Carlo search algorithm that uses a knowledge-based force field. Its versatile structure enables a variety of protein-modeling protocols, including purely de novo folding, folding guided by restraints derived from template structures or, structure assembly based on experimental data. In particular, CABS-NMR uses the distance and angular restraints set derived from various NMR experiments. This new modeling technique was successfully tested in structure determination of 10 globular proteins of size up to 216 residues, for which sparse NMR data were available. Additional detailed analysis was performed for a S100A1 protein. Namely, we successfully predicted Nuclear Overhauser Effect signals on the basis of low-energy structures obtained from chemical shifts by CABS-NMR. It has been observed that utility of chemical shifts and other types of experimental data (i.e. residual dipolar couplings and methyl-methyl Nuclear Overhauser Effect signals) in the presented modeling pipeline depends mainly on size of a protein and complexity of its topology. In this work, we have provided tools for either post-experiment processing of various kinds of NMR data or fast and low-cost structural analysis in the still challenging field of new fold predictions., (Copyright © 2010 Wiley Periodicals, Inc.)

Background: Several different methods for contact prediction succeeded within the Sixth Critical Assessment of Techniques for Protein Structure Prediction (CASP6). The most relevant were non-local contact predictions for targets from the most difficult categories: fold recognition-analogy and new fold. Such contacts could provide valuable structural information in case a template structure cannot be found in the PDB., Results: We described comprehensive tests of the effectiveness of contact data in various aspects of de novo modeling with CABS, an algorithm which was used successfully in CASP6 by the Kolinski-Bujnicki group. We used the predicted contacts in a simple scoring function for the post-simulation ranking of protein models and as a soft bias in the folding simulations and in the fold-refinement procedure. The latter approach turned out to be the most successful. The CABS force field used in the Replica Exchange Monte Carlo simulations cooperated with the true contacts and discriminated the false ones, which resulted in an improvement of the majority of Kolinski-Bujnicki's protein models. In the modeling we tested different sets of predicted contact data submitted to the CASP6 server. According to our results, the best performing were the contacts with the accuracy balanced with the coverage, obtained either from the best two predictors only or by a consensus from as many predictors as possible., Conclusion: Our tests have shown that theoretically predicted contacts can be very beneficial for protein structure prediction. Depending on the protein modeling method, a contact data set applied should be prepared with differently balanced coverage and accuracy of predicted contacts. Namely, high coverage of contact data is important for the model ranking and high accuracy for the folding simulations.

Routine structure prediction of new folds is still a challenging task for computational biology. The challenge is not only in the proper determination of overall fold but also in building models of acceptable resolution, useful for modeling the drug interactions and protein-protein complexes. In this work we propose and test a comprehensive approach to protein structure modeling supported by sparse, and relatively easy to obtain, experimental data. We focus on chemical shift-based restraints from NMR, although other sparse restraints could be easily included. In particular, we demonstrate that combining the typical NMR software with artificial intelligence-based prediction of secondary structure enhances significantly the accuracy of the restraints for molecular modeling. The computational procedure is based on the reduced representation approach implemented in the CABS modeling software, which proved to be a versatile tool for protein structure prediction during the CASP (CASP stands for critical assessment of techniques for protein structure prediction) experiments (see http://predictioncenter/CASP6/org). The method is successfully tested on a small set of representative globular proteins of different size and topology, including the two CASP6 targets, for which the required NMR data already exist. The method is implemented in a semi-automated pipeline applicable to a large scale structural annotation of genomic data. Here, we limit the computations to relatively small set. This enabled, without a loss of generality, a detailed discussion of various factors determining accuracy of the proposed approach to the protein structure prediction., ((c) 2007 Wiley Periodicals, Inc.)