Citius, Altius, Fortius.

2020 is a leap year. That means that we have one day extra and, if the Olympic games had survived the corona crisis, we would all be watching television and ask the eternal question whether Olympic records will for ever be broken and broken again, or that there are limits to human biology1 . In this article we ask the same question, but rather than discussing aspects of Citius, Altius, and Fortius of athletes we will discuss them for macromolecules. It is remarkable how many parallels can be found between Olympic records in these two seemingly different worlds. People involved in structure validation and re-refinement try to make us believe that most aspects of macromolecular structures can be caught by a number that has some constant value with little variation around it. We will show here that the PDB2 databank proves this idea to be wrong. In the protein structure world, it holds for many that "participating is more important than winning", but some, fortunately, still go for the record books.

Bio-knowledge-based filters improve residue-residue contact prediction accuracy.

Motivation: Residue-residue contact prediction through direct coupling analysis has reached impressive accuracy, but yet higher accuracy will be needed to allow for routine modelling of protein structures. One way to improve the prediction accuracy is to filter predicted contacts using knowledge about the particular protein of interest or knowledge about protein structures in general. Results: We focus on the latter and discuss a set of filters that can be used to remove false positive contact predictions. Each filter depends on one or a few cut-off parameters for which the filter performance was investigated. Combining all filters while using default parameters resulted for a test set of 851 protein domains in the removal of 29% of the predictions of which 92% were indeed false positives. Availability and implementation: All data and scripts are available at http://comprec-lin.iiar.pwr.edu.pl/FPfilter/. Supplementary information: Supplementary data are available at Bioinformatics online.

Forecasting residue-residue contact prediction accuracy.

MOTIVATION: Apart from meta-predictors, most of today's methods for residue-residue contact prediction are based entirely on Direct Coupling Analysis (DCA) of correlated mutations in multiple sequence alignments (MSAs). These methods are on average ∼40% correct for the 100 strongest predicted contacts in each protein. The end-user who works on a single protein of interest will not know if predictions are either much more or much less correct than 40%, which is especially a problem if contacts are predicted to steer experimental research on that protein. RESULTS: We designed a regression model that forecasts the accuracy of residue-residue contact prediction for individual proteins with an average error of 7 percentage points. Contacts were predicted with two DCA methods (gplmDCA and PSICOV). The models were built on parameters that describe the MSA, the predicted secondary structure, the predicted solvent accessibility and the contact prediction scores for the target protein. Results show that our models can be also applied to the meta-methods, which was tested on RaptorX. AVAILABILITY AND IMPLEMENTATION: All data and scripts are available from http://comprec-lin.iiar.pwr.edu.pl/dcaQ/. CONTACT: malgorzata.kotulska@pwr.edu.pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Correlated mutations select misfolded from properly folded proteins.

MOTIVATION: The recently developed direct coupling analysis (DCA) method has greatly improved the accuracy with which residue-residue contacts can be predicted from multiple sequence alignments. Contact prediction accuracy, though, is still often not sufficient for complete ab initio protein structure prediction. DCA can, however, support protein structure studies in several ways. RESULTS: We show that DCA can select the better structure from among properly folded and misfolded variants. This idea was tested by comparing obsolete PDB files with their more correctly folded successors and by the comparison of structures with deliberately misfolded decoy models from the Decoys 'R' Us database. The DCA method systematically predicts more contacts for properly folded structures than for misfolded ones. The method works much better for X-ray structures than for NMR structures. AVAILABILITY AND IMPLEMENTATION: All data are available from http://comprec-lin.iiar.pwr.edu.pl/dcaVSmisfolds/ and http://swift.cmbi.ru.nl/dcaVSmisfolds/ . CONTACT: malgorzata.kotulska@pwr.edu.pl . SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel IRF6 Mutations Detected in Orofacial Cleft Patients by Targeted Massively Parallel Sequencing.

Common variants in interferon regulatory factor 6 ( IRF6) have been associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) as well as with tooth agenesis (TA). These variants contribute a small risk towards the 2 congenital conditions and explain only a small percentage of heritability. On the other hand, many IRF6 mutations are known to be a monogenic cause of disease for syndromic orofacial clefting (OFC). We hypothesize that IRF6 mutations in some rare instances could also cause nonsyndromic OFC. To find novel rare variants in IRF6 responsible for nonsyndromic OFC and TA, we performed targeted multiplex sequencing using molecular inversion probes (MIPs) in 1,072 OFC patients, 67 TA patients, and 706 controls. We identified 3 potentially pathogenic de novo mutations in OFC patients. In addition, 3 rare missense variants were identified, for which pathogenicity could not unequivocally be shown, as all variants were either inherited from an unaffected parent or the parental DNA was not available. Retrospective investigation of the patients with these variants revealed the presence of lip pits in one of the patients with a de novo mutation suggesting a Van der Woude syndrome (VWS) phenotype, whereas, in other patients, no lip pits were identified.

GPCRdb: the G protein-coupled receptor database - an introduction.

GPCRs make up the largest family of human membrane proteins and of drug targets. Recent advances in GPCR pharmacology and crystallography have shed new light on signal transduction, allosteric modulation and biased signalling, translating into new mechanisms and principles for drug design. The GPCR database, GPCRdb, has served the community for over 20 years and has recently been extended to include a more multidisciplinary audience. This review is intended to introduce new users to the services in GPCRdb, which meets three overall purposes: firstly, to provide reference data in an integrated, annotated and structured fashion, with a focus on sequences, structures, single-point mutations and ligand interactions. Secondly, to equip the community with a suite of web tools for swift analysis of structures, sequence similarities, receptor relationships, and ligand target profiles. Thirdly, to facilitate dissemination through interactive diagrams of, for example, receptor residue topologies, phylogenetic relationships and crystal structure statistics. Herein, these services are described for the first time; visitors and guides are provided with good practices for their utilization. Finally, we describe complementary databases cross-referenced by GPCRdb and web servers with corresponding functionality.

The role of the second and third extracellular loops of the adenosine A1 receptor in activation and allosteric modulation.

The adenosine A1 receptor is a member of the large membrane protein family that signals through G proteins, the G protein-coupled receptors (GPCRs). GPCRs consist of seven transmembrane domains connected by three intracellular and three extracellular loops. Their N-terminus is extracellular, the C-terminal tail is in the cytoplasm. The transmembrane domains in receptor subfamilies that bind the same endogenous ligand, such as dopamine or adenosine, tend to be highly similar. In contrast, the loop regions can vary greatly, both in sequence and in length, and the role these loops have in the activation mechanism of the receptors remains unclear. Here, we investigated the activating role of the second and third extracellular loop of the human adenosine A1 receptor. By means of an (Ala)3 mutagenic scan in which consecutive sets of three amino acids were mutated into alanine residues in EL2 and a classical alanine scan in EL3, we revealed a strong regulatory role for the second extracellular loop (EL2) of the human adenosine A1 receptor. Besides many residues in the second and the third extracellular loops important for adenosine A1 receptor activation, we also identified two residues in EL2, a tryptophan and a glutamate, that affect the influence of the allosteric modulator PD81,723. These results, combined with a comparison of the different receptor loop regions, provide insight in the activation mechanism of this typical class A GPCR and further emphasize the unique pharmacological profile the loops can provide to individual receptors, even within subfamilies of GPCRs.

AsteriX: a Web server to automatically extract ligand coordinates from figures in PDF articles.

Coordinates describing the chemical structures of small molecules that are potential ligands for pharmaceutical targets are used at many stages of the drug design process. The coordinates of the vast majority of ligands can be obtained from either publicly accessible or commercial databases. However, interesting ligands sometimes are only available from the scientific literature, in which case their coordinates need to be reconstructed manually--a process that consists of a series of time-consuming steps. We present a Web server that helps reconstruct the three-dimensional (3D) coordinates of ligands for which a two-dimensional (2D) picture is available in a PDF file. The software, called AsteriX, analyses every picture contained in the PDF file and attempts to determine automatically whether or not it contains ligands. Areas in pictures that may contain molecular structures are processed to extract connectivity and atom type information that allow coordinates to be subsequently reconstructed. The AsteriX Web server was tested on a series of articles containing a large diversity in graphical representations. In total, 88% of 3249 ligand structures present in the test set were identified as chemical diagrams. Of these, about half were interpreted correctly as 3D structures, and a further one-third required only minor manual corrections. It is principally impossible to always correctly reconstruct 3D coordinates from pictures because there are many different protocols for drawing a 2D image of a ligand, but more importantly a wide variety of semantic annotations are possible. The AsteriX Web server therefore includes facilities that allow the users to augment partial or partially correct 3D reconstructions. All 3D reconstructions are submitted, checked, and corrected by the users domain at the server and are freely available for everybody. The coordinates of the reconstructed ligands are made available in a series of formats commonly used in drug design research. The AsteriX Web server is freely available at http://swift.cmbi.ru.nl/bitmapb/.

Drug design for ever, from hype to hope.

In its first 25 years JCAMD has been disseminating a large number of techniques aimed at finding better medicines faster. These include genetic algorithms, COMFA, QSAR, structure based techniques, homology modelling, high throughput screening, combichem, and dozens more that were a hype in their time and that now are just a useful addition to the drug-designers toolbox. Despite massive efforts throughout academic and industrial drug design research departments, the number of FDA-approved new molecular entities per year stagnates, and the pharmaceutical industry is reorganising accordingly. The recent spate of industrial consolidations and the concomitant move towards outsourcing of research activities requires better integration of all activities along the chain from bench to bedside. The next 25 years will undoubtedly show a series of translational science activities that are aimed at a better communication between all parties involved, from quantum chemistry to bedside and from academia to industry. This will above all include understanding the underlying biological problem and optimal use of all available data.

Antibody engineering reveals the important role of J segments in the production efficiency of llama single-domain antibodies in Saccharomyces cerevisiae.

Variable domains of llama heavy-chain antibodies (VHH) are becoming a potent tool for a wide range of biotechnological and medical applications. Because of structural features typical of their single-domain nature, they are relatively easy to produce in lower eukaryotes, but it is not uncommon that some molecules have poor secretion efficiency. We therefore set out to study the production yield of VHH. We computationally identified five key residues that are crucial for folding and secretion, and we validated their importance with systematic site-directed mutations. The observation that all key residues were localised in the V segment, in proximity of the J segment of VHH, led us to study the importance of J segment in secretion efficiency. Intriguingly, we found that the use of specific J segments in VHH could strongly influence the production yield. Sequence analysis and expression experiments strongly suggested that interactions with chaperones, especially with the J segment, are a crucial aspect of the production yield of VHH.

The role of calcium ions in the stability and instability of a thermolysin-like protease.

Thermolysin and other secreted broad-specificity proteases, such as subtilisin or alpha-lytic protease, are produced as pre-pro-proteins that stay at least partially unfolded while in the cytosol. After secretion, the pro-proteases fold to their active conformations in a process that includes the autolytic removal of the pro-peptide. We review the life cycle of the thermolysin-like protease from Bacillus stearothermophilus in light of the calcium dependent stability and instability of the N-terminal domain. The protease binds calcium ions in the regions that are involved in the autolytic maturation process. It is generally assumed that the calcium ions contribute to the extreme stability of the protease, but experimental evidence for TLP-ste indicates that at least one of the calcium ions plays a regulatory role. We hypothesize that this calcium ion plays an important role as a switch that modulates the protease between stable and unstable states as appropriate to the biological need.

Systematic generation of in vivo G protein-coupled receptor mutants in the rat.

G-protein-coupled receptors (GPCRs) constitute a large family of cell surface receptors that are involved in a wide range of physiological and pathological processes, and are targets for many therapeutic interventions. However, genetic models in the rat, one of the most widely used model organisms in physiological and pharmacological research, are largely lacking. Here, we applied N-ethyl-N-nitrosourea (ENU)-driven target-selected mutagenesis to generate an in vivo GPCR mutant collection in the rat. A pre-selected panel of 250 human GPCR homologs was screened for mutations in 813 rats, resulting in the identification of 131 non-synonymous mutations. From these, seven novel potential rat gene knockouts were established as well as 45 lines carrying missense mutations in various genes associated with or involved in human diseases. We provide extensive in silico modeling results of the missense mutations and show experimental data, suggesting loss-of-function phenotypes for several models, including Mc4r and Lpar1. Taken together, the approach used resulted not only in a set of novel gene knockouts, but also in allelic series of more subtle amino acid variants, similar as commonly observed in human disease. The mutants presented here may greatly benefit studies to understand specific GPCR function and support the development of novel therapeutic strategies.

WIWS: a protein structure bioinformatics Web service collection.

The WHAT IF molecular-modelling and drug design program is widely distributed in the world of protein structure bioinformatics. Although originally designed as an interactive application, its highly modular design and inbuilt control language have recently enabled its deployment as a collection of programmatically accessible web services. We report here a collection of WHAT IF-based protein structure bioinformatics web services: these relate to structure quality, the use of symmetry in crystal structures, structure correction and optimization, adding hydrogens and optimizing hydrogen bonds and a series of geometric calculations. The freely accessible web services are based on the industry standard WS-I profile and the EMBRACE technical guidelines, and are available via both REST and SOAP paradigms. The web services run on a dedicated computational cluster; their function and availability is monitored daily.

An active registry for bioinformatics web services.

SUMMARY: The EMBRACE Registry is a web portal that collects and monitors web services according to test scripts provided by the their administrators. Users are able to search for, rank and annotate services, enabling them to select the most appropriate working service for inclusion in their bioinformatics analysis tasks. AVAILABILITY AND IMPLEMENTATION: Web site implemented with PHP, Python, MySQL and Apache, with all major browsers supported. (www.embraceregistry.net).

Modeling GPCRs.

Many GPCR models have been built over the years for many different purposes, of which drug-design undoubtedly has been the most frequent one. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn things for the models we build today. We conclude that the GPCR modeling field is riddled with "common knowledge". Several characteristics of the bovine rhodopsin structure came as a big surprise, and had obviously not been predicted, which led to large errors in the models. Some of these surprises, however, could have been predicted if the modelers had more rigidly stuck to the rule that holds for all models, namely that a model should explain all experimental facts, and not just those facts that agree with the modeler's preconceptions.

MRS: a fast and compact retrieval system for biological data.

The biological data explosion of the 'omics' era requires fast access to many data types in rapidly growing data banks. The MRS server allows for very rapid queries in a large number of flat-file data banks, such as EMBL, UniProt, OMIM, dbEST, PDB, KEGG, etc. This server combines a fast and reliable backend with a very user-friendly implementation of all the commonly used information retrieval facilities. The MRS server is freely accessible at http://mrs.cmbi.ru.nl/. Moreover, the MRS software is freely available at http://mrs.cmbi.ru.nl/download/ for those interested in making their own data banks available via a web-based server.

Heavier-than-air flying machines are impossible.

Many G protein-coupled receptor (GPCR) models have been built over the years. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn more about the models we build today. We conclude that the GPCR modelling field is riddled with 'common knowledge' similar to Lord Kelvin's remark in 1895 that "heavier-than-air flying machines are impossible", and we summarize what we think are the (im)possibilities of modelling GPCRs using the coordinates of bovine rhodopsin as a template. Associated WWW pages: www.gpcr.org/articles/2003_mod

Correlated mutation analyses on very large sequence families.

The 'omics era' (the era of genomics, proteomics, and so forth) is marked by a flood of data that need to be interpreted to become useful information. Thanks to genome sequencing projects, large numbers of sequence families with more than a thousand members each are now available. Novel analytical techniques are needed to deal with this avalanche of sequence data. Sequence entropy is a measure of the information present in an alignment, whereas sequence variability represents the mutational flexibility at a particular position. Entropy versus variability plots can reveal the roles of groups of residues in the overall function of a protein. Such roles can be as part of the main active site, part of a modulator binding site, or transduction of a signal between those sites. Residues that are involved in a common function tend to stay conserved as a group, but when they mutate, they tend to mutate together. Correlated mutation analysis can detect groups of residue positions that show this behaviour. The combination of entropy, variability and correlation is a powerful tool to convert sequence data into useful information. This analysis can, for example, detect the key residues involved in cooperativity in globins, the switch regions in ras-like proteins and the calcium binding and signalling residues in serine proteases. We have extrapolated from these three classes of structurally and functionally well-described proteins to G-protein-coupled receptors (GPCRs). We can detect the residues in the main functional site in GPCRs that are responsible for G-protein coupling, the residues in the endogenous agonist binding site, and the residues in between that transduce the signal to and fro between these sites. The results are discussed in the light of a simple two-step evolutionary model for the development of functional proteins.

Completion and refinement of 3-D homology models with restricted molecular dynamics: application to targets 47, 58, and 111 in the CASP modeling competition and posterior analysis.

A method is presented to refine models built by homology by the use of restricted molecular dynamics (MD) techniques. The basic idea behind this method is the use of structure validation software to determine for each residue the likelihood that it is modeled correctly. This information is used to determine constraints and restraints in an MD simulation including explicit solvent molecules, which is used for model refinement. The procedure is based on the idea that residues that the validation software identifies as correctly positioned should be strongly constrained or restrained in the MD simulations, whereas residues that are likely to be positioned wrongly should move freely. Two different protocols are compared: one (applied to CASP3 target T58) using full structural constraints with separate optimization of each short fragment and the other (applied to T47) allowing some freedom using harmonic restraining potentials, with automatic optimization of the whole molecule. Structures along the MD trajectory that scored best in structural checks were selected for the construction of models that appeared to be successful in the CASP3 competition. Model refinement with MD in general leads to a model that is less like the experimental structure (Levitt et al. Nature Struct Biol 1999;6:108-111). Actually, refined T47 was slightly improved compared to the starting model; changes in model T58 led not to further enhancement. After the X-ray structure of the modeled proteins became known, the procedure was evaluated for two targets (T47 and the CASP4 target T111) by comparing a long simulation in water with the experimental target structures. It was found that structural improvements could be obtained on a nanosecond time scale by allowing appropriate freedom in the simulation. Structural checks applied to fast fluctuations do not appear to be informative for the correctness of the structure. However, both a simple hydrogen bond count and a simple compactness measure, if averaged over times of typically 300 ps, correlate well with structural correctness and we suggest that criteria based on these properties may be used in computational folding strategies.

Four-helix bundle topology re-engineered: monomeric Rop protein variants with different loop arrangements.

We converted the small homodimeric four-helix bundle repressor of primer protein (Rop) into a monomeric four-helix bundle by introduction of connecting loops. Both left- and right-handed four-helix bundles were produced. The left-handed bundles were more stable and were used to introduce biologically interesting peptides in one of the loops.