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.

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.

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.

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.

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.

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.

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).

The site of attachment in human rhinovirus 14 for antiviral agents that inhibit uncoating.

WIN 51711 and WIN 52084 are structurally related, antiviral compounds that inhibit the replication of rhino (common cold) viruses and related picornaviruses. They prevent the pH-mediated uncoating of the viral RNA. The compounds consist of a 3-methylisoxazole group that inserts itself into the hydrophobic interior of the VP1 beta-barrel, a connecting seven-membered aliphatic chain, and a 4-oxazolinylphenoxy group (OP) that covers the entrance to an ion channel in the floor of the "canyon." Viral disassembly may be inhibited by preventing the collapse of the VP1 hydrophobic pocket or by blocking the flow of ions into the virus interior.

The atomic structure of Mengo virus at 3.0 A resolution.

The structure of Mengo virus, a representative member of the cardio picornaviruses, is substantially different from the structures of rhino- and polioviruses. The structure of Mengo virus was solved with the use of human rhinovirus 14 as an 8 A resolution structural approximation. Phase information was then extended to 3 A resolution by use of the icosahedral symmetry. This procedure gives promise that many other virus structures also can be determined without the use of the isomorphous replacement technique. Although the organization of the major capsid proteins VP1, VP2, and VP3 of Mengo virus is essentially the same as in rhino- and polioviruses, large insertions and deletions, mostly in VP1, radically alter the surface features. In particular, the putative receptor binding "canyon" of human rhinovirus 14 becomes a deep "pit" in Mengo virus because of polypeptide insertions in VP1 that fill part of the canyon. The minor capsid peptide, VP4, is completely internal in Mengo virus, but its association with the other capsid proteins is substantially different from that in rhino- or poliovirus. However, its carboxyl terminus is located at a position similar to that in human rhinovirus 14 and poliovirus, suggesting the same autocatalytic cleavage of VP0 to VP4 and VP2 takes place during assembly in all these picornaviruses.

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.

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.

Collecting and harvesting biological data: the GPCRDB and NucleaRDB information systems.

The amount of genomic and proteomic data that is entered each day into databases and the experimental literature is outstripping the ability of experimental scientists to keep pace. While generic databases derived from automated curation efforts are useful, most biological scientists tend to focus on a class or family of molecules and their biological impact. Consequently, there is a need for molecular class-specific or other specialized databases. Such databases collect and organize data around a single topic or class of molecules. If curated well, such systems are extremely useful as they allow experimental scientists to obtain a large portion of the available data most relevant to their needs from a single source. We are involved in the development of two such databases with substantial pharmacological relevance. These are the GPCRDB and NucleaRDB information systems, which collect and disseminate data related to G protein-coupled receptors and intra-nuclear hormone receptors, respectively. The GPCRDB was a pilot project aimed at building a generic molecular class-specific database capable of dealing with highly heterogeneous data. A first version of the GPCRDB project has been completed and it is routinely used by thousands of scientists. The NucleaRDB was started recently as an application of the concept for the generalization of this technology. The GPCRDB is available via the WWW at http://www.gpcr.org/7tm/ and the NucleaRDB at http://www.receptors.org/NR/.

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.

Role of the N-terminal part of the coat protein in the assembly of cowpea chlorotic mottle virus. A 500 MHz proton nuclear magnetic resonance study and structural calculations.

The interaction of the oligonucleotides (Ap)8A and (A-T)5 with empty capsids of the coat protein of cowpea chlorotic mottle virus (CCMV) has been studied with 500 MHz 1H nuclear magnetic resonance. It is found that these oligonucleotides specifically bind to the arginine and lysine residues of the N-terminal arm of the protein. Upon this binding, immobilization of part of the N-terminal arm occurs. In addition, secondary structure predictions and energy calculations have been performed on the N-terminal arm. These calculations were carried out as a function of the charges on the arginine and lysine side-chains. For free coat protein, where the arginine and lysine side-chains are charged, the arm is found in a random-coil conformation. In the neutralized state, as for the coat protein in the virus, the arm adopts an alpha-helical conformation. The results support a previously published model for the assembly of CCMV, in which a random-coil to alpha-helix conformational transition, induced by neutralizing the arginine and lysine side-chains, plays an essential role.

Rubredoxin reductase of Pseudomonas oleovorans. Structural relationship to other flavoprotein oxidoreductases based on one NAD and two FAD fingerprints.

The oxidation of alkanes to alkanols by Pseudomonas oleovorans involves a three-component enzyme system: alkane hydroxylase, rubredoxin and rubredoxin reductase. Alkane hydroxylase and rubredoxin are encoded by the alkBFGHJKL operon, while previous studies indicated that rubredoxin reductase is most likely encoded on the second alk cluster: the alkST operon. In this study we show that alkT encodes the 41 x 10(3) Mr rubredoxin reductase, on the basis of a comparison of the expected amino acid composition of AlkT and the previously established amino acid composition of the purified rubredoxin reductase. The alkT sequence revealed significant similarities between AlkT and several NAD(P)H and FAD-containing reductases and dehydrogenases. All of these enzymes contain two ADP binding sites, which can be recognized by a common beta alpha beta-fold or fingerprint, derived from known structures of cofactor binding enzymes. By means of this amino acid fingerprint we were able to determine that one ADP binding site in rubredoxin reductase (AlkT) is located at the N terminus and is involved in FAD binding, while the second site is located in the middle of the sequence and is involved in the binding of NAD or NADP. In addition, we derived from the sequences of FAD binding reductases a second amino acid fingerprint for FAD binding, and we used this fingerprint to identify a third amino acid sequence in AlkT near the carboxy terminus for binding of the flavin moiety of FAD. On the basis of the known architecture and relative spatial orientations of the NAD and FAD binding sites in related dehydrogenases, a model for part of the tertiary structure of AlkT was developed.

Conformational changes in the chaperonin GroEL: new insights into the allosteric mechanism.

Conformational changes are known to play a crucial role in the function of the bacterial GroE chaperonin system. Here, results are presented from an essential dynamics analysis of known experimental structures and from computer simulations of GroEL using the CONCOORD method. The results indicate a possible direct form of inter-ring communication associated with internal fluctuations in the nucleotide-binding domains upon nucleotide and GroES binding that are involved in the allosteric mechanism of GroEL. At the level of conformational transitions in entire GroEL rings, nucleotide-induced structural changes were found to be distinct and in principle uncoupled from changes occurring upon GroES binding. However, a coupling is found between nucleotide-induced conformational changes and GroES-mediated transitions, but only in simulations of GroEL double rings, and not in simulations of single rings. This provides another explanation for the fact that GroEL functions a double ring system.

Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex.

Triosephosphate isomerase (TIM) is a dimeric glycolytic enzyme. TIM from Trypanosoma brucei brucei has been crystallized at pH 7.0 in 2.4 M-ammonium sulphate. The well-diffracting crystals have one dimer per asymmetric unit. The structure has been refined at 1.83 A resolution with an R-factor of 18.3% for all data between 6 A and 1.83 A (37,568 reflections). The model consists of 3778 protein atoms and 297 solvent atoms. Subunit 1 is involved in considerably more crystal contacts than subunit 2. Correlated with these differences in crystal packing is the observation that only in the active site of subunit 2 is a sulphate ion bound. Furthermore, significant differences with respect to structure and flexibility are observed in three loops near the active site. In particular, there is a 7 A positional difference of the tip of the flexible loop (loop 6) when comparing subunit 1 and subunit 2. Also, the neighbouring loops (loop 5 and loop 7) have significantly different conformations and flexibility. In subunit 1, loop 6 is in an "open" conformation, in subunit 2, loop 6 is in an "almost closed" conformation. Only in the presence of a phosphate-containing ligand, such as glycerol-3-phosphate, does loop 6 take up the "closed" conformation. Loop 6 and loop 7 (and also to some extent loop 5) are rather flexible in the almost closed conformation, but well defined in the open and closed conformations. The closing of loop 6 (167 to 180), as observed in the almost closed conformation, slightly changes the main-chain conformation of the catalytic glutamate, Glu167, leading to a change of the chi 1 angle of this residue from approximately -60 degrees to approximately 60 degrees and the weakening of the hydrogen bonds between its polar side-chain atoms and Ser96. In the closed conformation, in the presence of glycerol-3-phosphate, the main-chain atoms of Glu167 remain in the same position as in the almost closed conformation, but the side-chain has rotated around the CA-CB bond changing chi 1 from approximately 60 degrees to approximately -60 degrees. In this new position the hydrogen bonding to Ser96 is completely lost and also a water-mediated salt bridge between OE2(Glu167) and NE(Arg99) is lost. Comparison of the two independently refined subunits, showed that the root-mean-square deviation for all 249 CA atoms is 0.9 A; for the CA atoms of the beta-strands this is only 0.2 A. The average B-factor for all subunit 1 and subunit 2 atoms is 20 A2 and 25 A2, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)

Crystal structure of the p-hydroxybenzoate hydroxylase-substrate complex refined at 1.9 A resolution. Analysis of the enzyme-substrate and enzyme-product complexes.

Using synchrotron radiation, the X-ray diffraction intensities of crystals of p-hydroxy-benzoate hydroxylase, complexed with the substrate p-hydroxybenzoate, were measured to a resolution of 1.9 A. Restrained least-squares refinement alternated with rebuilding in electron density maps yielded an atom model of the enzyme-substrate complex with a crystallographic R-factor of 15.6% for 31,148 reflections between 6.0 and 1.9 A. A total of 330 solvent molecules was located. In the final model, only three residues have deviating phi-psi angle combinations. One of them, the active site residue Arg44, has a well-defined electron density and may be strained to adopt this conformation for efficient catalysis. The mode of binding of FAD is distinctly different for the different components of the coenzyme. The adenine ring is engaged in three water-mediated hydrogen bonds with the protein, while making only one direct hydrogen bond with the enzyme. The pyrophosphate moiety makes five water-mediated versus three direct hydrogen bonds. The ribityl and ribose moieties make only direct hydrogen bonds, in all cases, except one, with side-chain atoms. The isoalloxazine ring also makes only direct hydrogen bonds, but virtually only with main-chain atoms. The conformation of FAD in p-hydroxybenzoate hydroxylase is strikingly similar to that in glutathione reductase, while the riboflavin-binding parts of these two enzymes have no structural similarity at all. The refined 1.9 A structure of the p-hydroxybenzoate hydroxylase-substrate complex was the basis of further refinement of the 2.3 A structure of the enzyme-product complex. The result was a final R-factor of 16.7% for 14,339 reflections between 6.0 and 2.3 A and an improved geometry. Comparison between the complexes indicated only small differences in the active site region, where the product molecule is rotated by 14 degrees compared with the substrate in the enzyme-substrate complex. During the refinements of the enzyme-substrate and enzyme-product complexes, the flavin ring was allowed to bend or twist by imposing planarity restraints on the benzene and pyrimidine ring, but not on the flavin ring as a whole. The observed angle between the benzene ring and the pyrimidine ring was 10 degrees for the enzyme-substrate complex and 19 degrees for the enzyme-product complex. Because of the high temperature factors of the flavin ring in the enzyme-product complex, the latter value should be treated with caution. Six out of eight peptide residues near the flavin ring are oriented with their nitrogen atom pointing towards the ring.(ABSTRACT TRUNCATED AT 400 WORDS)