The LassoProt server, http://lassoprot.cent.uw.edu.pl/, enables analysis of biopolymers with entangled configurations called lassos. The server offers various ways of visualizing lasso configurations, as well as their time trajectories, with all the results and plots downloadable. Broad spectrum of applications makes LassoProt a useful tool for biologists, biophysicists, chemists, polymer physicists and mathematicians. The server and our methods have been validated on the whole PDB, and the results constitute the database of proteins with complex lassos, supported with basic biological data. This database can serve as a source of information about protein geometry and entanglement-function correlations, as a reference set in protein modeling, and for many other purposes.

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.AvailabilityGPCRM server and database: http://gpcrm.biomodellab.eu

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.

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.

The CB1 cannabinoid receptor (CB1R) contains one of the longest N termini among class A G protein-coupled receptors. Mutagenesis studies suggest that the allosteric binding site of cannabidiol (CBD) involves residues from the N terminal domain. In order to study the allosteric binding of CBD to CB1R we modeled the whole N-terminus of this receptor using the replica exchange molecular dynamics with solute tempering (REST2) approach. Then, the obtained structures of CB1R with the N terminus were used for ligand docking. A natural cannabinoid receptor agonist, Δ9-THC, was docked to the orthosteric site and a negative allosteric modulator, CBD, to the allosteric site positioned between extracellular ends of helices TM1 and TM2. The molecular dynamics simulations were then performed for CB1R with ligands: (i) CBD together with THC, and (ii) THC-only. Analyses of the differences in the residue-residue interaction patterns between those two cases allowed us to elucidate the allosteric network responsible for the modulation of the CB1R by CBD. In addition, we identified the changes in the orthosteric binding mode of Δ9-THC, as well as the changes in its binding energy, caused by the CBD allosteric binding. We have also found that the presence of a complete N-terminal domain is essential for a stable binding of CBD in the allosteric site of CB1R as well as for the allosteric-orthosteric coupling mechanism.

GPCRsignal (https://gpcrsignal.biomodellab.eu/) is a webserver devoted to signaling complexes of G-protein–coupled receptors (GPCRs). The recent improvement in cryo-electron microscopy resulted in the determination of a large number of high-resolution structures of GPCRs bound to their effector proteins: G proteins or arrestins. Analyzing the interfaces between receptor and an effector protein is of high importance since a selection of proper G protein or specific conformation of arrestin leads to changes of signaling that can significantly affect action of drugs. GPCRsignal provides a possibility of running molecular dynamics simulations of all currently available GPCR-effector protein complexes for curated structures: wild-type, with crystal/cryo-EM mutations, or with mutations introduced by the user. The simulations are performed in an implicit water-membrane environment, so they are rather fast. User can run several simulations to obtain statistically valid results. The simulations can be analyzed separately using dynamic FlarePlots for particular types of interactions. One can also compare groups of simulations in Interaction frequency analysis as HeatMaps and also in interaction frequency difference analysis as sticks, linking the interacting residues, of different color and size proportional to differences in contact frequencies.

No abstract

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.