Trendspotting in the Protein Data Bank

Network medicine approaches for identification of novel prognostic systems biomarkers and drug candidates for papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is one of the most common endocrine carcinomas worldwide and the aetiology of this cancer is still not well understood. Therefore, it remains important to understand the disease mechanism and find prognostic biomarkers and/or drug candidates for PTC. Compared with approaches based on single-gene assessment, network medicine analysis offers great promise to address this need. Accordingly, in the present study, we performed differential co-expressed network analysis using five transcriptome datasets in patients with PTC and healthy controls. Following meta-analysis of the transcriptome datasets, we uncovered common differentially expressed genes (DEGs) for PTC and, using these genes as proxies, found a highly clustered differentially expressed co-expressed module: a 'PTC-module'. Using independent data, we demonstrated the high prognostic capacity of the PTC-module and designated this module as a prognostic systems biomarker. In addition, using the nodes of the PTC-module, we performed drug repurposing and text mining analyzes to identify novel drug candidates for the disease. We performed molecular docking simulations, and identified: 4-demethoxydaunorubicin hydrochloride, AS605240, BRD-A60245366, ER 27319 maleate, sinensetin, and TWS119 as novel drug candidates whose efficacy was also confirmed by in silico analyzes. Consequently, we have highlighted here the need for differential co-expression analysis to gain a systems-level understanding of a complex disease, and we provide candidate prognostic systems biomarker and novel drugs for PTC.

Illuminating the "Twilight Zone": Advances in Difficult Protein Modeling

Homology modeling was long considered a method of choice in tertiary protein structure prediction. However, it used to provide models of acceptable quality only when templates with appreciable sequence identity with a target could be found. The threshold value was long assumed to be around 20-30%. Below this level, obtained sequence identity was getting dangerously close to values that can be obtained by chance, after aligning any random, unrelated sequences. In these cases, other approaches, including ab initio folding simulations or fragment assembly, were usually employed. The most recent editions of the CASP and CAMEO community-wide modeling methods assessment have brought some surprising outcomes, proving that much more clues can be inferred from protein sequence analyses than previously thought. In this chapter, we focus on recent advances in the field of difficult protein modeling, pushing the threshold deep into the "twilight zone", with particular attention devoted to improvements in applications of machine learning and model evaluation.

Crystallography-based exploration of non-covalent interactions for the design and synthesis of coumarin for stronger protein binding

Protein molecules are a good target for the inhibition or promotion of biological processes. Different methods like QSAR and molecular docking have been developed to accurately design small binder molecules for target proteins. An alternative model has been developed wherein a statistical method is used to find the propensity of different non-covalent interactions between small molecules and amino acid residues of the protein. The results give hints as to the choice of substituents required at the SM to strongly bind to a protein. In this case, 75 different types of proteins bound with coumarin derivatives have been investigated and the non-covalent interactions observed between the basic coumarin moiety and amino acids have been analyzed. Density functional theory (DFT) calculations were used to identify the electronic features of coumarin to understand the feasibility of the observed non-covalent interactions and to find appropriate groups that can modulate these interactions. The binding affinity towards a protein (β-lactoglobulin (BLG)) and the stability of the protein complex have been investigated through docking and molecular dynamics of 100 ns, respectively. The modeled compounds were synthesized and investigated with regards to their interactions with the model carrier protein. The thermodynamics of the interactions were also investigated and the binding is governed by the Le Chatelier principle.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Protein oligomer structure prediction using GALAXY in CASP14

Proteins perform their functions by interacting with other biomolecules. For these interactions, proteins often form homo- or hetero-oligomers as well. Thus, oligomer protein structures provide important clues regarding the biological roles of proteins. To this end, computational prediction of oligomer structures may be a useful tool in the absence of experimentally resolved structures. Here, we describe our server and human-expert methods used to predict oligomer structures in the CASP14 experiment. Examples are provided for cases in which manual domain-splitting led to improved oligomeric domain structures by ab initio docking, automated oligomer structure refinement led to improved subunit orientation and terminal structure, and manual oligomer modeling utilizing literature information generated a reasonable oligomer model. We also discussed the results of post-prediction docking calculations with AlphaFold2 monomers as input in comparison to our blind prediction results. Overall, ab initio docking of AlphaFold2 models did not lead to better oligomer structure prediction, which may be attributed to the interfacial structural difference between the AlphaFold2 monomer structures and the crystal oligomer structures. This result poses a next-stage challenge in oligomer structure prediction after the success of AlphaFold2. For successful protein assembly structure prediction, a different approach that exploits further evolutionary information on the interface and/or flexible docking taking the interfacial conformational flexibilities of subunit structures into account is needed.

A crystallography-based investigation of weak interactions for drug design against COVID-19

Interactions between proteins and small molecules play important roles in the inhibition of protein function. However, a lack of proper knowledge about non-covalent interactions can act as a barrier towards gaining a complete understanding of the factors that control these associations. To find effective molecules for COVID-19 inhibition, we have quantitatively investigated 143 X-ray crystal structures of the SARS-CoV-2 M^pro protein of coronavirus with covalently or non-covalently bound small molecules (SMs). Our present study is able to explain ordinary and perceptive aspects relating to protein inhibition. The active site of the protein consists of 21 amino acid residues, but only nine are actively involved in the ligand binding process. The H41, M49, and C145 residues have highest priority with respect to interactions with small molecules through hydrogen bond, CH-π, and van der Waals interactions. At the active site, this ranking of amino acids is clear, based on different spatial orientations of ligands, and consistent with the electronic properties. SMs with aromatic moieties that bind to the active site of the protein play a distinct role in the determination of the following order of interaction frequency with the amino acids: CH-π > H-bonding > polar interactions. This present study revealed that the G143 and C145 residues play crucial roles in the recognition of the carbonyl functionality of SMs through hydrogen bonding. With this knowledge in mind, an effective inhibitor small-molecule for SARS-CoV-2 M^pro was designed: docking studies showed that the designed molecule has strong binding affinity towards the protein. The non-covalent interactions in the protein-ligand complex are in good agreement with the results obtained from X-ray crystallography. Moreover, the present study focused on weak forces and their influence on protein inhibition, henceforth shedding much light on the essential requirements for moieties that should be present in a good inhibitor and their orientations at the ligand binding site.

Beyond history and "on a roll": The list of the most well-studied human protein structures and overall trends in the protein data bank

Of the roughly 20,000 canonical human protein sequences, as of January 20, 2021, 7,077 proteins have had their full or partial, medium- to high-resolution structures determined by x-ray crystallography or other methods. Which of these proteins dominate the protein data bank (the PDB) and why? In this paper, we list the 273 top human protein structures based on the number of their PDB entries. This set of proteins accounts for more than 40% of all available human PDB entries and represent past trends as well as current status for protein structural biology. We briefly discuss the relationship which some of the prominent protein structures have with protein research as a whole and mention their relevance to human diseases. The top-10 soluble and membrane proteins are all well-known (most of their first structures being deposited more than 30 years ago). Overall, there is no dramatic change in recent trends in the PDB. Remarkably, the number of structure depositions has grown nearly exponentially over the last 10 or more years (with a doubling time of 7 years for proteins, obtained from any organism). Growth in human protein structures is slightly faster (at 5.9 years). The information in this paper may be informative to senior scientists but also inspire researchers who are new to protein science, providing the year 2021 snap-shot for the state of protein structural biology.

Structure prediction of biological assemblies using GALAXY in CAPRI rounds 38-45

We participated in CARPI rounds 38-45 both as a server predictor and a human predictor. These CAPRI rounds provided excellent opportunities for testing prediction methods for three classes of protein interactions, that is, protein-protein, protein-peptide, and protein-oligosaccharide interactions. Both template-based methods (GalaxyTBM for monomer protein, GalaxyHomomer for homo-oligomer protein, GalaxyPepDock for protein-peptide complex) and ab initio docking methods (GalaxyTongDock and GalaxyPPDock for protein oligomer, GalaxyPepDock-ab-initio for protein-peptide complex, GalaxyDock2 and Galaxy7TM for protein-oligosaccharide complex) have been tested. Template-based methods depend heavily on the availability of proper templates and template-target similarity, and template-target difference is responsible for inaccuracy of template-based models. Inaccurate template-based models could be improved by our structure refinement and loop modeling methods based on physics-based energy optimization (GalaxyRefineComplex and GalaxyLoop) for several CAPRI targets. Current ab initio docking methods require accurate protein structures as input. Small conformational changes from input structure could be accounted for by our docking methods, producing one of the best models for several CAPRI targets. However, predicting large conformational changes involving protein backbone is still challenging, and full exploration of physics-based methods for such problems is still to come.

Blind prediction of homo- and hetero-protein complexes: The CASP13-CAPRI experiment

We present the results for CAPRI Round 46, the third joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of 20 targets including 14 homo-oligomers and 6 heterocomplexes. Eight of the homo-oligomer targets and one heterodimer comprised proteins that could be readily modeled using templates from the Protein Data Bank, often available for the full assembly. The remaining 11 targets comprised 5 homodimers, 3 heterodimers, and two higher-order assemblies. These were more difficult to model, as their prediction mainly involved "ab-initio" docking of subunit models derived from distantly related templates. A total of ~30 CAPRI groups, including 9 automatic servers, submitted on average ~2000 models per target. About 17 groups participated in the CAPRI scoring rounds, offered for most targets, submitting ~170 models per target. The prediction performance, measured by the fraction of models of acceptable quality or higher submitted across all predictors groups, was very good to excellent for the nine easy targets. Poorer performance was achieved by predictors for the 11 difficult targets, with medium and high quality models submitted for only 3 of these targets. A similar performance "gap" was displayed by scorer groups, highlighting yet again the unmet challenge of modeling the conformational changes of the protein components that occur upon binding or that must be accounted for in template-based modeling. Our analysis also indicates that residues in binding interfaces were less well predicted in this set of targets than in previous Rounds, providing useful insights for directions of future improvements.

Automatic recognition of ligands in electron density by machine learning

Motivation

The correct identification of ligands in crystal structures of protein complexes is the cornerstone of structure-guided drug design. However, cognitive bias can sometimes mislead investigators into modeling fictitious compounds without solid support from the electron density maps. Ligand identification can be aided by automatic methods, but existing approaches are based on time-consuming iterative fitting.

Results

Here we report a new machine learning algorithm called CheckMyBlob that identifies ligands from experimental electron density maps. In benchmark tests on portfolios of up to 219 931 ligand binding sites containing the 200 most popular ligands found in the Protein Data Bank, CheckMyBlob markedly outperforms the existing automatic methods for ligand identification, in some cases doubling the recognition rates, while requiring significantly less time. Our work shows that machine learning can improve the automation of structure modeling and significantly accelerate the drug screening process of macromolecule-ligand complexes.

Availability and implementation

Code and data are available on GitHub at https://github.com/dabrze/CheckMyBlob.

Supplementary information

Supplementary data are available at Bioinformatics online.

A comparative anatomy of protein crystals: lessons from the automatic processing of 56 000 samples

The fully automatic processing of crystals of macromolecules has presented a unique opportunity to gather information on the samples that is not usually recorded. This has proved invaluable in improving sample-location, characterization and data-collection algorithms. After operating for four years, MASSIF-1 has now processed over 56 000 samples, gathering information at each stage, from the volume of the crystal to the unit-cell dimensions, the space group, the quality of the data collected and the reasoning behind the decisions made in data collection. This provides an unprecedented opportunity to analyse these data together, providing a detailed landscape of macromolecular crystals, intimate details of their contents and, importantly, how the two are related. The data show that mosaic spread is unrelated to the size or shape of crystals and demonstrate experimentally that diffraction intensities scale in proportion to crystal volume and molecular weight. It is also shown that crystal volume scales inversely with molecular weight. The results set the scene for the development of X-ray crystallography in a changing environment for structural biology.

Evaluation of template-based modeling in CASP13

Performance in the template‐based modeling (TBM) category of CASP13 is assessed here, using a variety of metrics. Performance of the predictor groups that participated is ranked using the primary ranking score that was developed by the assessors for CASP12. This reveals that the best results are obtained by groups that include contact predictions or inter‐residue distance predictions derived from deep multiple sequence alignments. In cases where there is a good homolog in the wwPDB (TBM‐easy category), the best results are obtained by modifying a template. However, for cases with poorer homologs (TBM‐hard), very good results can be obtained without using an explicit template, by deep learning algorithms trained on the wwPDB. Alternative metrics are introduced, to allow testing of aspects of structural models that are not addressed by traditional CASP metrics. These include comparisons to the main‐chain and side‐chain torsion angles of the target, and the utility of models for solving crystal structures by the molecular replacement method. The alternative metrics are poorly correlated with the traditional metrics, and it is proposed that modeling has reached a sufficient level of maturity that the best models should be expected to satisfy this wider range of criteria.

Structural and dynamical description of the enzymatic reaction of a phosphohexomutase

Enzymes are known to adopt various conformations at different points along their catalytic cycles. Here, we present a comprehensive analysis of 15 isomorphous, high resolution crystal structures of the enzyme phosphoglucomutase from the bacterium Xanthomonas citri. The protein was captured in distinct states critical to function, including enzyme-substrate, enzyme-product, and enzyme-intermediate complexes. Key residues in ligand recognition and regions undergoing conformational change are identified and correlated with the various steps of the catalytic reaction. In addition, we use principal component analysis to examine various subsets of these structures with two goals: (1) identifying sites of conformational heterogeneity through a comparison of room temperature and cryogenic structures of the apo-enzyme and (2) a priori clustering of the enzyme-ligand complexes into functionally related groups, showing sensitivity of this method to structural features difficult to detect by traditional methods. This study captures, in a single system, the structural basis of diverse substrate recognition, the subtle impact of covalent modification, and the role of ligand-induced conformational change in this representative enzyme of the α-D-phosphohexomutase superfamily.

Improving inverse docking target identification with Z-score selection

The utilization of inverse docking methods for target identification has been driven by an increasing demand for efficient tools for detecting potential drug side-effects. Despite impressive achievements in the field of inverse docking, identifying true positives from a pool of potential targets still remains challenging. Notably, most of the developed techniques have low accuracies, limit the pool of possible targets that can be investigated or are not easy to use for non-experts due to a lack of available scripts or webserver. Guided by our finding that the absolute docking score was a poor indication of a ligand's protein target, we developed a novel "combined Z-score" method that used a weighted fraction of ligand and receptor-based Z-scores to identify the most likely binding target of a ligand. With our combined Z-score method, an additional 14%, 3.6%, and 6.3% of all ligand-protein pairs of the Astex, DUD, and DUD-E databases, respectively, were correctly predicted compared to a docking score-based selection. The combined Z-score had the highest area under the curve in a ROC curve analysis of all three datasets and the enrichment factor for the top 1% predictions using the combined Z-score analysis was the highest for the Astex and DUD-E datasets. Additionally, we developed a user-friendly python script (compatible with both Python2 and Python3) that enables users to employ the combined Z-score analysis for target identification using a user-defined list of ligands and targets. We are providing this python script and a user tutorial as part of the supplemental information.

Taxonomic Landscape of the Dark Proteomes: Whole-Proteome Scale Interplay Between Structural Darkness, Intrinsic Disorder, and Crystallization Propensity

Growth rate of the protein sequence universe dramatically exceeds the speed of expansion for the protein structure universe, generating an immense dark proteome that includes proteins with unknown structure. A whole-proteome scale analysis of 5.4 million proteins from 987 proteomes in the three domains of life and viruses to systematically dissect an interplay between structural coverage, degree of putative intrinsic disorder, and predicted propensity for structure determination is performed. It has been found that Archaean and Bacterial proteomes have relatively high structural coverage and low amounts of disorder, whereas Eukaryotic and Viral proteomes are characterized by a broad spread of structural coverage and higher disorder levels. The analysis reveals that dark proteomes (i.e., proteomes containing high fractions of proteins with unknown structure) have significantly elevated amounts of intrinsic disorder and are predicted to be difficult to solve structurally. Although the majority of dark proteomes are of viral origin, many dark viral proteomes have at least modest crystallization propensity and only a handful of them are enriched in the intrinsic disorder. The disorder, structural coverage, and propensity are mapped for structural determination onto a novel proteome-level sequence similarity network to analyze the interplay of these characteristics in the taxonomic landscape.

A pipeline for computational design of novel RNA-like topologies

Designing novel RNA topologies is a challenge, with important therapeutic and industrial applications. We describe a computational pipeline for design of novel RNA topologies based on our coarse-grained RNA-As-Graphs (RAG) framework. RAG represents RNA structures as tree graphs and describes RNA secondary (2D) structure topologies (currently up to 13 vertices, ≈260 nucleotides). We have previously identified novel graph topologies that are RNA-like among these. Here we describe a systematic design pipeline and illustrate design for six broad design problems using recently developed tools for graph-partitioning and fragment assembly (F-RAG). Following partitioning of the target graph, corresponding atomic fragments from our RAG-3D database are combined using F-RAG, and the candidate atomic models are scored using a knowledge-based potential developed for 3D structure prediction. The sequences of the top scoring models are screened further using available tools for 2D structure prediction. The results indicate that our modular approach based on RNA-like topologies rather than specific 2D structures allows for greater flexibility in the design process, and generates a large number of candidate sequences quickly. Experimental structure probing using SHAPE-MaP for two sequences agree with our predictions and suggest that our combined tools yield excellent candidates for further sequence and experimental screening.

Accurately Predicting Disordered Regions of Proteins Using Rosetta ResidueDisorder Application

Although many proteins necessitate well-folded structures to properly instigate their biological functions, a large fraction of functioning proteins contain regions-known as intrinsically disordered protein regions-where stable structures are not likely to form. Notable functional roles of intrinsically disordered proteins are in transcriptional regulation, translation, and cellular signal transduction. Moreover, intrinsically disordered protein regions are highly abundant in many proteins associated with various human diseases, therefore these segments have become attractive drug targets for potential therapeutics. Over the past decades, numerous computational methods have been developed to accurately predict disordered regions of proteins. Here we introduce a user-friendly and reliable approach for the prediction of disordered protein regions using the structure prediction software Rosetta. Using 245 proteins from a benchmark data set (16 DisProt database proteins) and a test data set (229 proteins with NMR data), we use Rosetta to predict the global protein structures and then show that there is a statistically significant difference between Rosetta scores in disordered and ordered regions, with scores being less favorable in disordered regions. Furthermore, the difference in scores between ordered and disordered protein regions is sufficient to accurately identify disordered protein regions. As a result, our Rosetta ResidueDisorder method (benchmark data set prediction accuracy of 71.77% and independent test data set prediction accuracy of 65.37%) outperformed other established disorder prediction tools and did not exhibit a biased prediction toward either ordered or disordered regions. To facilitate usage, a Rosetta application has been developed for the Rosetta ResidueDisorder method.

Structural database resources for biological macromolecules

This Briefing reviews the widely used, currently active, up-to-date databases derived from the worldwide Protein Data Bank (PDB) to facilitate browsing, finding and exploring its entries. These databases contain visualization and analysis tools tailored to specific kinds of molecules and interactions, often including also complex metrics precomputed by experts or external programs, and connections to sequence and functional annotation databases. Importantly, updates of most of these databases involves steps of curation and error checks based on specific expertise about the subject molecules or interactions, and removal of sequence redundancy, both leading to better data sets for mining studies compared with the full list of raw PDB entries. The article presents the databases in groups such as those aimed to facilitate browsing through PDB entries, their molecules and their general information, those built to link protein structure with sequence and dynamics, those specific for transmembrane proteins, nucleic acids, interactions of biomacromolecules with each other and with small molecules or metal ions, and those concerning specific structural features or specific protein families. A few webservers directly connected to active databases, and a few databases that have been discontinued but would be important to have back, are also briefly commented on. Along the Briefing, sample cases where these databases have been used to aid structural studies or advance our knowledge about biological macromolecules are referenced. A few specific examples are also given where using these databases is easier and more informative than using raw PDB data.

The PPI3D web server for searching, analyzing and modeling protein-protein interactions in the context of 3D structures

Summary

The PPI3D web server is focused on searching and analyzing the structural data on protein-protein interactions. Reducing the data redundancy by clustering and analyzing the properties of interaction interfaces using Voronoi tessellation makes this software a highly effective tool for addressing different questions related to protein interactions.

Availability and Implementation

The server is freely accessible at http://bioinformatics.lt/software/ppi3d/ .

Contact

ceslovas.venclovas@bti.vu.lt.

Supplementary information

Supplementary data are available at Bioinformatics online.

The challenge of modeling protein assemblies: the CASP12-CAPRI experiment

We present the quality assessment of 5613 models submitted by predictor groups from both CAPRI and CASP for the total of 15 most tractable targets from the second joint CASP-CAPRI protein assembly prediction experiment. These targets comprised 12 homo-oligomers and 3 hetero-complexes. The bulk of the analysis focuses on 10 targets (of CAPRI Round 37), which included all 3 hetero-complexes, and whose protein chains or the full assembly could be readily modeled from structural templates in the PDB. On average, 28 CAPRI groups and 10 CASP groups (including automatic servers), submitted models for each of these 10 targets. Additionally, about 16 groups participated in the CAPRI scoring experiments. A range of acceptable to high quality models were obtained for 6 of the 10 Round 37 targets, for which templates were available for the full assembly. Poorer results were achieved for the remaining targets due to the lower quality of the templates available for the full complex or the individual protein chains, highlighting the unmet challenge of modeling the structural adjustments of the protein components that occur upon binding or which must be accounted for in template-based modeling. On the other hand, our analysis indicated that residues in binding interfaces were correctly predicted in a sizable fraction of otherwise poorly modeled assemblies and this with higher accuracy than published methods that do not use information on the binding partner. Lastly, the strengths and weaknesses of the assessment methods are evaluated and improvements suggested.

Iterative Molecular Dynamics-Rosetta Membrane Protein Structure Refinement Guided by Cryo-EM Densities

Knowing atomistic details of proteins is essential not only for the understanding of protein function but also for the development of drugs. Experimental methods such as X-ray crystallography, NMR, and cryo-electron microscopy (cryo-EM) are the preferred forms of protein structure determination and have achieved great success over the most recent decades. Computational methods may be an alternative when experimental techniques fail. However, computational methods are severely limited when it comes to predicting larger macromolecule structures with little sequence similarity to known structures. The incorporation of experimental restraints in computational methods is becoming increasingly important to more reliably predict protein structure. One such experimental input used in structure prediction and refinement is cryo-EM densities. Recent advances in cryo-EM have arguably revolutionized the field of structural biology. Our previously developed cryo-EM-guided Rosetta-MD protocol has shown great promise in the refinement of soluble protein structures. In this study, we extended cryo-EM density-guided iterative Rosetta-MD to membrane proteins. We also improved the methodology in general by picking models based on a combination of their score and fit-to-density during the Rosetta model selection. By doing so, we have been able to pick models superior to those with the previous selection based on Rosetta score only and we have been able to further improve our previously refined models of soluble proteins. The method was tested with five membrane spanning protein structures. By applying density-guided Rosetta-MD iteratively we were able to refine the predicted structures of these membrane proteins to atomic resolutions. We also showed that the resolution of the density maps determines the improvement and quality of the refined models. By incorporating high-resolution density maps (∼4 Å), we were able to more significantly improve the quality of the models than when medium-resolution maps (6.9 Å) were used. Beginning from an average starting structure root mean square deviation (RMSD) to native of 4.66 Å, our protocol was able to refine the structures to bring the average refined structure RMSD to 1.66 Å when 4 Å density maps were used. The protocol also successfully refined the HIV-1 CTD guided by an experimental 5 Å density map.

Augmenting the anisotropic network model with torsional potentials improves PATH performance, enabling detailed comparison with experimental rate data

PATH algorithms for identifying conformational transition states provide computational parameters-time to the transition state, conformational free energy differences, and transition state activation energies-for comparison to experimental data and can be carried out sufficiently rapidly to use in the "high throughput" mode. These advantages are especially useful for interpreting results from combinatorial mutagenesis experiments. This report updates the previously published algorithm with enhancements that improve correlations between PATH convergence parameters derived from virtual variant structures generated by RosettaBackrub and previously published kinetic data for a complete, four-way combinatorial mutagenesis of a conformational switch in Tryptophanyl-tRNA synthetase.

A close-up view of codon selection in eukaryotic initiation

When given an option to choose among a set of alternatives and only one selection is right, one might stop and reflect over which one is best. However, the ribosome has no time to stop and make such reflections, proteins need to be produced and very fast. Eukaryotic translation initiation is an example of such a conundrum. Here, scanning for the correct codon match must be fast, efficient and accurate. We highlight our recent computational findings, which show how the initiation machinery manages to recognize one specific codon among many possible challengers, by fine-tuning the energetic landscape of base-pairing with the aid of the initiation factors eIF1 and eIF1A. Using a recent 3-dimensional structure of the eukaryotic initiation complex we have performed simulations of codon recognition in atomic detail. These calculations provide an in-depth energetic and structural view of how discrimination against near-cognate codons is achieved by the initiation complex.

Modeling protein-protein and protein-peptide complexes: CAPRI 6th edition

We present the sixth report evaluating the performance of methods for predicting the atomic resolution structures of protein complexes offered as targets to the community-wide initiative on the Critical Assessment of Predicted Interactions (CAPRI). The evaluation is based on a total of 20,670 predicted models for 8 protein-peptide complexes, a novel category of targets in CAPRI, and 12 protein-protein targets in CAPRI prediction Rounds held during the years 2013-2016. For two of the protein-protein targets, the focus was on the prediction of side-chain conformation and positions of interfacial water molecules. Seven of the protein-protein targets were particularly challenging owing to their multicomponent nature, to conformational changes at the binding site, or to a combination of both. Encouragingly, the very large multiprotein complex with the nucleosome was correctly predicted, and correct models were submitted for the protein-peptide targets, but not for some of the challenging protein-protein targets. Models of acceptable quality or better were obtained for 14 of the 20 targets, including medium quality models for 13 targets and high quality models for 8 targets, indicating tangible progress of present-day computational methods in modeling protein complexes with increased accuracy. Our evaluation suggests that the progress stems from better integration of different modeling tools with docking procedures, as well as the use of more sophisticated evolutionary information to score models. Nonetheless, adequate modeling of conformational flexibility in interacting proteins remains an important area with a crucial need for improvement. Proteins 2017; 85:359-377. © 2016 Wiley Periodicals, Inc.

Synchrotron-based macromolecular crystallography module for an undergraduate biochemistry laboratory course

This paper describes the introduction of synchrotron-based macromolecular crystallography (MX) into an undergraduate laboratory class. An introductory 2 week experimental module on MX, consisting of four laboratory sessions and two classroom lectures, was incorporated into a senior-level biochemistry class focused on a survey of biochemical techniques, including the experimental characterization of proteins. Students purified recombinant protein samples, set up crystallization plates and flash-cooled crystals for shipping to a synchrotron. Students then collected X-ray diffraction data sets from their crystals via the remote interface of the Molecular Biology Consortium beamline (4.2.2) at the Advanced Light Source in Berkeley, CA, USA. Processed diffraction data sets were transferred back to the laboratory and used in conjunction with partial protein models provided to the students for refinement and model building. The laboratory component was supplemented by up to 2 h of lectures by faculty with expertise in MX. This module can be easily adapted for implementation into other similar undergraduate classes, assuming the availability of local crystallographic expertise and access to remote data collection at a synchrotron source.

The Intersection of Structural and Chemical Biology - An Essential Synergy

The continual improvement in our ability to generate high resolution structural models of biological molecules has stimulated and supported innovative chemical biology projects that target increasingly challenging ligand interaction sites. In this review we outline some of the recent developments in chemical biology and rational ligand design and show selected examples that illustrate the synergy between these research areas.

The RCSB protein data bank: integrative view of protein, gene and 3D structural information

The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB, http://rcsb.org), the US data center for the global PDB archive, makes PDB data freely available to all users, from structural biologists to computational biologists and beyond. New tools and resources have been added to the RCSB PDB web portal in support of a ‘Structural View of Biology.’ Recent developments have improved the User experience, including the high-speed NGL Viewer that provides 3D molecular visualization in any web browser, improved support for data file download and enhanced organization of website pages for query, reporting and individual structure exploration. Structure validation information is now visible for all archival entries. PDB data have been integrated with external biological resources, including chromosomal position within the human genome; protein modifications; and metabolic pathways. PDB-101 educational materials have been reorganized into a searchable website and expanded to include new features such as the Geis Digital Archive.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

Biological membranes

Biological membranes allow life as we know it to exist. They form cells and enable separation between the inside and outside of an organism, controlling by means of their selective permeability which substances enter and leave. By allowing gradients of ions to be created across them, membranes also enable living organisms to generate energy. In addition, they control the flow of messages between cells by sending, receiving and processing information in the form of chemical and electrical signals. This essay summarizes the structure and function of membranes and the proteins within them, and describes their role in trafficking and transport, and their involvement in health and disease. Techniques for studying membranes are also discussed.

Modulating Mobility: a Paradigm for Protein Engineering?

Proteins are highly mobile structures. In addition to gross conformational changes occurring on, for example, ligand binding, they are also subject to constant thermal motion. The mobility of a protein varies through its structure and can be modulated by ligand binding and other events. It is becoming increasingly clear that this mobility plays an important role in key functions of proteins including catalysis, allostery, cooperativity, and regulation. Thus, in addition to an optimum structure, proteins most likely also require an optimal dynamic state. Alteration of this dynamic state through protein engineering will affect protein function. A dramatic example of this is seen in some inherited metabolic diseases where alternation of residues distant from the active site affects the mobility of the protein and impairs function. We postulate that using molecular dynamics simulations, experimental data or a combination of the two, it should be possible to engineer the mobility of active sites. This may be useful in, for example, increasing the promiscuity of enzymes. Thus, a paradigm for protein engineering is suggested in which the mobility of the active site is rationally modified. This might be combined with more "traditional" approaches such as altering functional groups in the active site.

Application of advanced X-ray methods in life sciences

BACKGROUND

Synchrotron radiation (SR) sources provide diverse X-ray methods for the investigation of structure-function relationships in biological macromolecules.

SCOPE OF REVIEW

Recent developments in SR sources and in the X-ray tools they offer for life sciences are reviewed. Specifically, advances in macromolecular crystallography, small angle X-ray solution scattering, X-ray absorption and fluorescence spectroscopy, and imaging are discussed with examples.

MAJOR CONCLUSIONS

SR sources offer a range of X-ray techniques that can be used in a complementary fashion in studies of biological systems at a wide range of resolutions from atomic to cellular scale. Emerging applications of X-ray techniques include the characterization of disordered proteins, noncrystalline and nonequilibrium systems, elemental imaging of tissues, cells and organs, and detection of time-resolved changes in molecular structures.

GENERAL SIGNIFICANCE

X-ray techniques are in the center of hybrid approaches that are used to gain insight into complex problems relating to biomolecular mechanisms, disease and possible therapeutic solutions. This article is part of a Special Issue entitled "Science for Life". Guest Editors: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Fitmunk: improving protein structures by accurate, automatic modeling of side-chain conformations

Improvements in crystallographic hardware and software have allowed automated structure-solution pipelines to approach a near-`one-click' experience for the initial determination of macromolecular structures. However, in many cases the resulting initial model requires a laborious, iterative process of refinement and validation. A new method has been developed for the automatic modeling of side-chain conformations that takes advantage of rotamer-prediction methods in a crystallographic context. The algorithm, which is based on deterministic dead-end elimination (DEE) theory, uses new dense conformer libraries and a hybrid energy function derived from experimental data and prior information about rotamer frequencies to find the optimal conformation of each side chain. In contrast to existing methods, which incorporate the electron-density term into protein-modeling frameworks, the proposed algorithm is designed to take advantage of the highly discriminatory nature of electron-density maps. This method has been implemented in the program Fitmunk, which uses extensive conformational sampling. This improves the accuracy of the modeling and makes it a versatile tool for crystallographic model building, refinement and validation. Fitmunk was extensively tested on over 115 new structures, as well as a subset of 1100 structures from the PDB. It is demonstrated that the ability of Fitmunk to model more than 95% of side chains accurately is beneficial for improving the quality of crystallographic protein models, especially at medium and low resolutions. Fitmunk can be used for model validation of existing structures and as a tool to assess whether side chains are modeled optimally or could be better fitted into electron density. Fitmunk is available as a web service at http://kniahini.med.virginia.edu/fitmunk/server/ or at http://fitmunk.bitbucket.org/.

Carbohydrate-Aromatic Interactions in Proteins

Protein-carbohydrate interactions play pivotal roles in health and disease. However, defining and manipulating these interactions has been hindered by an incomplete understanding of the underlying fundamental forces. To elucidate common and discriminating features in carbohydrate recognition, we have analyzed quantitatively X-ray crystal structures of proteins with noncovalently bound carbohydrates. Within the carbohydrate-binding pockets, aliphatic hydrophobic residues are disfavored, whereas aromatic side chains are enriched. The greatest preference is for tryptophan with an increased prevalence of 9-fold. Variations in the spatial orientation of amino acids around different monosaccharides indicate specific carbohydrate C-H bonds interact preferentially with aromatic residues. These preferences are consistent with the electronic properties of both the carbohydrate C-H bonds and the aromatic residues. Those carbohydrates that present patches of electropositive saccharide C-H bonds engage more often in CH-π interactions involving electron-rich aromatic partners. These electronic effects are also manifested when carbohydrate-aromatic interactions are monitored in solution: NMR analysis indicates that indole favorably binds to electron-poor C-H bonds of model carbohydrates, and a clear linear free energy relationships with substituted indoles supports the importance of complementary electronic effects in driving protein-carbohydrate interactions. Together, our data indicate that electrostatic and electronic complementarity between carbohydrates and aromatic residues play key roles in driving protein-carbohydrate complexation. Moreover, these weak noncovalent interactions influence which saccharide residues bind to proteins, and how they are positioned within carbohydrate-binding sites.

Fully automatic characterization and data collection from crystals of biological macromolecules

Considerable effort is dedicated to evaluating macromolecular crystals at synchrotron sources, even for well established and robust systems. Much of this work is repetitive, and the time spent could be better invested in the interpretation of the results. In order to decrease the need for manual intervention in the most repetitive steps of structural biology projects, initial screening and data collection, a fully automatic system has been developed to mount, locate, centre to the optimal diffraction volume, characterize and, if possible, collect data from multiple cryocooled crystals. Using the capabilities of pixel-array detectors, the system is as fast as a human operator, taking an average of 6 min per sample depending on the sample size and the level of characterization required. Using a fast X-ray-based routine, samples are located and centred systematically at the position of highest diffraction signal and important parameters for sample characterization, such as flux, beam size and crystal volume, are automatically taken into account, ensuring the calculation of optimal data-collection strategies. The system is now in operation at the new ESRF beamline MASSIF-1 and has been used by both industrial and academic users for many different sample types, including crystals of less than 20 µm in the smallest dimension. To date, over 8000 samples have been evaluated on MASSIF-1 without any human intervention.

Structural Perspectives on the Evolutionary Expansion of Unique Protein-Protein Binding Sites

Structures of protein complexes provide atomistic insights into protein interactions. Human proteins represent a quarter of all structures in the Protein Data Bank; however, available protein complexes cover less than 10% of the human proteome. Although it is theoretically possible to infer interactions in human proteins based on structures of homologous protein complexes, it is still unclear to what extent protein interactions and binding sites are conserved, and whether protein complexes from remotely related species can be used to infer interactions and binding sites. We considered biological units of protein complexes and clustered protein-protein binding sites into similarity groups based on their structure and sequence, which allowed us to identify unique binding sites. We showed that the growth rate of the number of unique binding sites in the Protein Data Bank was much slower than the growth rate of the number of structural complexes. Next, we investigated the evolutionary roots of unique binding sites and identified the major phyletic branches with the largest expansion in the number of novel binding sites. We found that many binding sites could be traced to the universal common ancestor of all cellular organisms, whereas relatively few binding sites emerged at the major evolutionary branching points. We analyzed the physicochemical properties of unique binding sites and found that the most ancient sites were the largest in size, involved many salt bridges, and were the most compact and least planar. In contrast, binding sites that appeared more recently in the evolution of eukaryotes were characterized by a larger fraction of polar and aromatic residues, and were less compact and more planar, possibly due to their more transient nature and roles in signaling processes.

X-ray crystallography over the past decade for novel drug discovery - where are we heading next?

INTRODUCTION

Macromolecular X-ray crystallography has been the primary methodology for determining the three-dimensional structures of proteins, nucleic acids and viruses. Structural information has paved the way for structure-guided drug discovery and laid the foundations for structural bioinformatics. However, X-ray crystallography still has a few fundamental limitations, some of which may be overcome and complemented using emerging methods and technologies in other areas of structural biology.

AREAS COVERED

This review describes how structural knowledge gained from X-ray crystallography has been used to advance other biophysical methods for structure determination (and vice versa). This article also covers current practices for integrating data generated by other biochemical and biophysical methods with those obtained from X-ray crystallography. Finally, the authors articulate their vision about how a combination of structural and biochemical/biophysical methods may improve our understanding of biological processes and interactions.

EXPERT OPINION

X-ray crystallography has been, and will continue to serve as, the central source of experimental structural biology data used in the discovery of new drugs. However, other structural biology techniques are useful not only to overcome the major limitation of X-ray crystallography, but also to provide complementary structural data that is useful in drug discovery. The use of recent advancements in biochemical, spectroscopy and bioinformatics methods may revolutionize drug discovery, albeit only when these data are combined and analyzed with effective data management systems. Accurate and complete data management is crucial for developing experimental procedures that are robust and reproducible.

Ligand-target prediction by structural network biology using nAnnoLyze

Target identification is essential for drug design, drug-drug interaction prediction, dosage adjustment and side effect anticipation. Specifically, the knowledge of structural details is essential for understanding the mode of action of a compound on a target protein. Here, we present nAnnoLyze, a method for target identification that relies on the hypothesis that structurally similar binding sites bind similar ligands. nAnnoLyze integrates structural information into a bipartite network of interactions and similarities to predict structurally detailed compound-protein interactions at proteome scale. The method was benchmarked on a dataset of 6,282 pairs of known interacting ligand-target pairs reaching a 0.96 of area under the Receiver Operating Characteristic curve (AUC) when using the drug names as an input feature for the classifier, and a 0.70 of AUC for “anonymous” compounds or compounds not present in the training set. nAnnoLyze resulted in higher accuracies than its predecessor, AnnoLyze. We applied the method to predict interactions for all the compounds in the DrugBank database with each human protein structure and provide examples of target identification for known drugs against human diseases. The accuracy and applicability of our method to any compound indicate that a comparative docking approach such as nAnnoLyze enables large-scale annotation and analysis of compound–protein interactions and thus may benefit drug development.

Description of the “mode-of-action” of a small chemical compound against a protein target is essential for the drug discovery process. Such description relies on three main steps: i) the identification of the target protein within the thousands of proteins in an organism, ii) the localization of the binding interaction site in the identified target protein, and iii) the molecular characterization of the compound’s binding mode in the binding site of the target protein. Here, we introduce a new computational method, called nAnnoLyze, which uses graph theory principles to relate compounds and target proteins based on comparative principles. nAnnoLyze aims at addressing two of the three previous steps, that is, target identification and binding site localization. Our results suggest that the nAnnoLyze accuracy and proteome-wide applicability enables the large-scale annotation and analysis of compound–protein interaction and thus may benefit drug development.

Quantitative and Systems-Based Approaches for Deciphering Bacterial Membrane Interactome and Gene Function

High-throughput genomic and proteomic methods provide a concise description of the molecular constituents of a cell, whereas systems biology strives to understand the way these components function as a whole. Recent developments, such as genome editing technologies and protein epitope-tagging coupled with high-sensitivity mass-spectrometry, allow systemic studies to be performed at an unprecedented scale. Available methods can be successfully applied to various goals, both expanding fundamental knowledge and solving applied problems. In this review, we discuss the present state and future of bacterial cell envelope interactomics, with a specific focus on host-pathogen interactions and drug target discovery. Both experimental and computational methods will be outlined together with examples of their practical implementation.

Plant-PrAS: a database of physicochemical and structural properties and novel functional regions in plant proteomes

Arabidopsis thaliana is an important model species for studies of plant gene functions. Research on Arabidopsis has resulted in the generation of high-quality genome sequences, annotations and related post-genomic studies. The amount of annotation, such as gene-coding regions and structures, is steadily growing in the field of plant research. In contrast to the genomics resource of animals and microorganisms, there are still some difficulties with characterization of some gene functions in plant genomics studies. The acquisition of information on protein structure can help elucidate the corresponding gene function because proteins encoded in the genome possess highly specific structures and functions. In this study, we calculated multiple physicochemical and secondary structural parameters of protein sequences, including length, hydrophobicity, the amount of secondary structure, the number of intrinsically disordered regions (IDRs) and the predicted presence of transmembrane helices and signal peptides, using a total of 208,333 protein sequences from the genomes of six representative plant species, Arabidopsis thaliana, Glycine max (soybean), Populus trichocarpa (poplar), Oryza sativa (rice), Physcomitrella patens (moss) and Cyanidioschyzon merolae (alga). Using the PASS tool and the Rosetta Stone method, we annotated the presence of novel functional regions in 1,732 protein sequences that included unannotated sequences from the Arabidopsis and rice proteomes. These results were organized into the Plant Protein Annotation Suite database (Plant-PrAS), which can be freely accessed online at http://plant-pras.riken.jp/.

Single nucleotide variations: biological impact and theoretical interpretation

Genome-wide association studies (GWAS) and whole-exome sequencing (WES) generate massive amounts of genomic variant information, and a major challenge is to identify which variations drive disease or contribute to phenotypic traits. Because the majority of known disease-causing mutations are exonic non-synonymous single nucleotide variations (nsSNVs), most studies focus on whether these nsSNVs affect protein function. Computational studies show that the impact of nsSNVs on protein function reflects sequence homology and structural information and predict the impact through statistical methods, machine learning techniques, or models of protein evolution. Here, we review impact prediction methods and discuss their underlying principles, their advantages and limitations, and how they compare to and complement one another. Finally, we present current applications and future directions for these methods in biological research and medical genetics.

Structural models for the design of novel antiviral agents against Greek Goat Encephalitis

The Greek Goat Encephalitis virus (GGE) belongs to the Flaviviridae family of the genus Flavivirus. The GGE virus constitutes an important pathogen of livestock that infects the goat’s central nervous system. The viral enzymes of GGE, helicase and RNA-dependent RNA polymerase (RdRP), are ideal targets for inhibitor design, since those enzymes are crucial for the virus’ survival, proliferation and transmission. In an effort to understand the molecular structure underlying the functions of those viral enzymes, the three dimensional structures of GGE NS3 helicase and NS5 RdRP have been modelled. The models were constructed in silico using conventional homology modelling techniques and the known 3D crystal structures of solved proteins from closely related species as templates. The established structural models of the GGE NS3 helicase and NS5 RdRP have been evaluated for their viability using a repertoire of in silico tools. The goal of this study is to present the 3D conformations of the GGE viral enzymes as reliable structural models that could provide the platform for the design of novel anti-GGE agents.

New structural insights into the apelin receptor: identification of key residues for apelin binding

Apelin is the endogenous ligand of the orphan 7-transmembrane domain GPCR APJ, now named the apelin receptor (ApelinR). Apelin plays a prominent role in body fluid and cardiovascular homeostasis. To better understand the structural organization of the ApelinR, we built 3 homology 3-dimensional (3D) models of the human ApelinR using the validated cholecystokinin receptor-1 3D model or the X-ray structures of the β2-adrenergic and CXCR4 receptors as templates. Docking of the pyroglutamyl form of apelin 13 (pE13F) into these models revealed the conservation at the bottom of the binding site of a hydrophobic cavity in which the C-terminal Phe of pE13F was embedded. In contrast, at the top of the binding site, depending on the model, different interactions were visualized between acidic residues of the ApelinR and the basic residues of pE13F. Using site-directed mutagenesis, we showed that Asp 92, Glu 172, and Asp 282 of rat ApelinR are key residues in apelin binding by interacting with Lys 8, Arg 2, and Arg 4 of pE13F, respectively. These residues are only seen in the CXCR4-based ApelinR 3D model, further validating this model. These findings bring new insights into the structural organization of the ApelinR and the mode of apelin binding.

Redundancy-weighting for better inference of protein structural features

MOTIVATION

Structural knowledge, extracted from the Protein Data Bank (PDB), underlies numerous potential functions and prediction methods. The PDB, however, is highly biased: many proteins have more than one entry, while entire protein families are represented by a single structure, or even not at all. The standard solution to this problem is to limit the studies to non-redundant subsets of the PDB. While alleviating biases, this solution hides the many-to-many relations between sequences and structures. That is, non-redundant datasets conceal the diversity of sequences that share the same fold and the existence of multiple conformations for the same protein. A particularly disturbing aspect of non-redundant subsets is that they hardly benefit from the rapid pace of protein structure determination, as most newly solved structures fall within existing families.

RESULTS

In this study we explore the concept of redundancy-weighted datasets, originally suggested by Miyazawa and Jernigan. Redundancy-weighted datasets include all available structures and associate them (or features thereof) with weights that are inversely proportional to the number of their homologs. Here, we provide the first systematic comparison of redundancy-weighted datasets with non-redundant ones. We test three weighting schemes and show that the distributions of structural features that they produce are smoother (having higher entropy) compared with the distributions inferred from non-redundant datasets. We further show that these smoothed distributions are both more robust and more correct than their non-redundant counterparts. We suggest that the better distributions, inferred using redundancy-weighting, may improve the accuracy of knowledge-based potentials and increase the power of protein structure prediction methods. Consequently, they may enhance model-driven molecular biology.