Capri Research Articles

PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations.

Predicting protein-peptide interactions is crucial for understanding peptide binding processes and designing peptide drugs. However, traditional computational modeling approaches face challenges in accurately predicting peptide-protein binding structures due to the slow dynamics and high flexibility of the peptides. Here, we introduce a new workflow termed "PepBinding" for predicting peptide binding structures, which combines peptide docking, all-atom enhanced sampling simulations using the Peptide Gaussian accelerated Molecular Dynamics (Pep-GaMD) method, and structural clustering. PepBinding has been demonstrated on seven distinct model peptides. In peptide docking using HPEPDOCK, the peptide backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures ranged from 3.8 to 16.0 Å, corresponding to the medium to inaccurate quality models according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. The Pep-GaMD simulations performed for only 200 ns significantly improved the docking models, resulting in five medium and two acceptable quality models. Therefore, PepBinding is an efficient workflow for predicting peptide binding structures and is publicly available at https://github.com/MiaoLab20/PepBinding.

A significance score for protein-protein interaction models through random docking.

Comparing accuracies of structural protein-protein interaction (PPI) models for different complexes on an absolute scale is a challenge, requiring normalization of scores across structures of different sizes and shapes. To help address this challenge, we have developed a statistical significance metric for docking models, called random-docking (RD) p-value. This score evaluates a PPI model based on how likely a random docking process is to produce a model of better or equal accuracy. The binding partners are randomly docked against each other a large number of times, and the probability of sampling a model of equal or greater accuracy from this reference distribution is the RD p-value. Using a subset of top predicted models from CAPRI (Critical Assessment of PRediction of Interactions) rounds over 2017-2020, we find that the ease of achieving a given root mean squared deviation or DOCKQ score varies considerably by target; achieving the same relative metric can be thousands of times easier for one complex compared to another. In contrast, RD p-values inherently normalize scores for models of different complexes, making them globally comparable. Furthermore, one can calculate RD p-values after generating a reference distribution that accounts for prior information about the interface geometry, such as residues involved in binding, by giving the random-docking process access the same information. Thus, one can decouple improvements in prediction accuracy that arise solely from basic modeling constraints from those due to the rest of the method. We provide efficient code for computing RD p-values at https://github.com/Grigoryanlab/RDP.

Computational Methods to Target Protein-Protein Interactions.

The chapter emphasizes the importance of understanding protein-protein interactions in cellular mechanisms and highlights the role of computational modeling in predicting these interactions. It discusses sequence-based approaches such as evolutionary trace (ET), correlated mutation analysis (CMA), and subtractive correlated mutation (SCM) for identifying crucial amino acid residues, considering interface conservation or evolutionary changes. The chapter also explores methods like differential ET, hidden-site class model, and spatial cluster detection (SCD) for interface specificity and spatial clustering. Furthermore, it examines approaches combining structural and sequential methodologies and evaluates modeled predictions through initiatives like critical assessment of prediction of interactions (CAPRI). Additionally, the chapter provides an overview of various software programs used for molecular docking, detailing their search, sampling, refinement and scoring stages, along with innovative techniques and tools like normal mode analysis (NMA) and adaptive Poisson-Boltzmann solver (APBS) for electrostatic calculations. These computational and experimental approaches are crucial for unraveling protein-protein interactions and aid in developing potential therapeutics for various diseases.

Could local foods achieve recommended calcium intakes for nutritionally vulnerable populations in Uganda, Guatemala, and Bangladesh?

Globally, dietary intake of calcium is often insufficient, and it is unclear if adequacy could be achieved by promoting calcium-rich local foods. This study used linear programming and household consumption data from Uganda, Bangladesh, and Guatemala to assess whether local foods could meet calcium population reference intakes (CaPRIs). The most promising food-based approaches to promote dietary calcium adequacy were identified for 12- to 23-month-old breastfed children, 4- to 6-year-old children, 10- to 14-year-old girls, and nonpregnant and nonbreastfeeding (NPNB) women of reproductive age living in two regions of each country. Calcium-optimized diets achieved 75-253% of the Ca PRI, depending on the population, and were <100% for 4- to 6-year-olds in one region of each country and 10- to 14-year-old girls in Sylhet, Bangladesh. The best food sources of calcium were green leafy vegetables and milk, across geographic locations, and species of small fish, nixtamalized (lime-treated) maize products, sesame seeds, and bean varieties, where consumed. Food-based recommendations (FBRs) achieving the minimum calcium threshold were identified for 12- to 23-month-olds and NPNB women across geographic locations, and for 4- to 6-year-olds and 10-to 14-year-old girls in Uganda. However, for 4- to 6-year-olds and 10- to 14-year-old girls in Bangladesh and Guatemala, calcium-adequate FBRs could not be identified, indicating a need for alternative calcium sources or increased access to and consumption of local calcium-rich foods.

Multivariate tools to investigate the spatial contaminant distribution in a highly anthropized area (Gulf of Naples, Italy)

The Gulf of Naples located in a high anthropized coastal area is subjected to an infrastructural intervention for the installation of a submarine power pipeline. In order to evaluate the distribution of contaminants in the seafloor sediments, a preliminary study has been conducted in the area using multivariate techniques. The statistic approach was performed to gain insights on the occurrence of organic and inorganic contaminants within the area, aiming to identify the relevant hot spots. Three geographical sub-areas influenced by different contaminant association were recognized: Torre Annunziata (TA), Capri (CA), and middle offshore (MO). TA and CA resulted marked by a severe contamination pattern due to anthropogenic pressures. In addition, the influence of the depositional basin in governing the contamination trend has been pointed out. The supervised technique PLS_DA resulted to be a powerful tool in addressing the complexity of the huge dataset acquired during the marine survey, highlighting the main trends in the variability of quality indicators, orienting thus the deeper investigations during follow-up monitoring activities.

An editorial perspective on Barbieri and Capri article: Is vaccination good value for money? A review of cost-utility analyses of vaccination strategies in eight European countries

An editorial perspective on Barbieri and Capri article: Is vaccination good value for money? A review of cost-utility analyses of vaccination strategies in eight European countries

Challenge Accepted – Paratope States in Solution Improve Structure Prediction and Docking

The rise of antibodies as biotherapeutic proteins has attracted the interest in characterizing and understanding the antibody binding interface. Structure-based antibody design and accurate predictions of antibody-antigen interactions remain major challenges in computational biology. By using molecular dynamics simulations in combination with enhanced sampling techniques, we show that a single static X-ray structure is not sufficient to identify determinants of antibody recognition and to characterize functionally relevant conformational rearrangements. Thus, we suggest that antibodies should be described as ensembles of paratope states in solution, which are defined by a combination of correlated loop movements, which favor distinct interface orientations and elbow angles. In this study, we focus our attention on characterizing paratope states in solution of antibodies which undergo substantial conformational changes upon antigen binding. Because of their conformational changes, they have been classified as “difficult” cases for antibody-antigen docking in the Critical Assessment of PRediction of Interactions (CAPRI) blind docking challenge. Here, we present transition kinetics and thermodynamics of these conformational rearrangements of the paratope and show that kinetically relevant states can be used to improve antibody-antigen docking. We use the Rosetta SnugDock docking algorithm to show that by even using the unbound antibody X-ray structure as starting structures for molecular dynamics simulations, a binding competent conformation substantially different to the unbound antibody X-ray structure, can be retained. We also observe that the kinetically dominant antibody conformations are chosen by the bound antigen conformation with the highest probability. Thus, we show that kinetically relevant antibody conformations obtained by a combination of molecular dynamics simulations and Markov-state models can improve antibody-antigen docking and structure prediction.

Template-based modeling and ab-initio docking using CoDock in CAPRI.

Integration of template-based modeling, global sampling and precise scoring is crucial for the development of molecular docking programs with improved accuracy. We combined template-based modeling and ab-initio docking protocol as hybrid docking strategy called CoDock for the docking and scoring experiments of the seventh CAPRI edition. For CAPRI rounds 38-45, we obtained acceptable or better models in the top 10 submissions for eight out of the 16 evaluated targets as predictors, nine out of the 16 targets as scorers. Especially, we submitted acceptable models for all of the evaluated protein-oligosaccharide targets. For the CASP13-CAPRI experiment (round 46), we obtained acceptable or better models in the top 5 submissions for 10 out of the 20 evaluated targets as predictors, 11 out of the 20 targets as scorers. The failed cases for our group were mainly the difficult targets and the protein-peptide systems in CAPRI and CASP13-CAPRI experiments. In summary, this CAPRI edition showed that our hybrid docking strategy can be efficiently adapted to the increasing variety of challenges in the field of molecular interactions.

“Capri Sun v. Riha”

“Capri Sun v. Riha”

Template-Based Modeling of Protein Complexes Using the PPI3D Web Server.

There is a large gap between the numbers of known protein-protein interactions and the corresponding experimentally solved structures of protein complexes. Fortunately, this gap can be in part bridged by computational structure modeling methods. Currently, template-based modeling is the most accurate means to predict both individual protein structures and protein complexes. One of the major issues in template-based modeling is to identify homologous structures that could be utilized as templates. To simplify this task, we have developed the PPI3D web server. The server is not only able to search for homologous protein complexes, but also provides means to analyze identified interactions and to model protein complexes. In recent CASP and CAPRI experiments, PPI3D proved to be a useful tool for homology modeling of multimeric proteins. In this chapter, we provide a brief description of the PPI3D web server capabilities and how to use the server for modeling of protein complexes.

Novel sampling strategies and a coarse-grained score function for docking homomers, flexible heteromers, and oligosaccharides using Rosetta in CAPRI rounds 37-45.

Critical Assessment of PRediction of Interactions (CAPRI) rounds 37 through 45 introduced larger complexes, new macromolecules, and multistage assemblies. For these rounds, we used and expanded docking methods in Rosetta to model 23 target complexes. We successfully predicted 14 target complexes and recognized and refined near-native models generated by other groups for two further targets. Notably, for targets T110 and T136, we achieved the closest prediction of any CAPRI participant. We created several innovative approaches during these rounds. Since round 39 (target 122), we have used the new RosettaDock 4.0, which has a revamped coarse-grained energy function and the ability to perform conformer selection during docking with hundreds of pregenerated protein backbones. Ten of the complexes had some degree of symmetry in their interactions, so we tested Rosetta SymDock, realized its shortcomings, and developed the next-generation symmetric docking protocol, SymDock2, which includes docking of multiple backbones and induced-fit refinement. Since the last CAPRI assessment, we also developed methods for modeling and designing carbohydrates in Rosetta, and we used them to successfully model oligosaccharide-protein complexes in round 41. Although the results were broadly encouraging, they also highlighted the pressing need to invest in (a) flexible docking algorithms with the ability to model loop and linker motions and in (b) new sampling and scoring methods for oligosaccharide-protein interactions.

Docking proteins and peptides under evolutionary constraints in Critical Assessment of PRediction of Interactions rounds 38 to 45.

Computational structural prediction of macromolecular interactions is a fundamental tool toward the global understanding of cellular processes. The Critical Assessment of PRediction of Interactions (CAPRI) community-wide experiment provides excellent opportunities for blind testing computational docking methods and includes original targets, thus widening the range of docking applications. Our participation in CAPRI rounds 38 to 45 enabled us to expand the way we include evolutionary information in structural predictions beyond our standard free docking InterEvDock pipeline. InterEvDock integrates a coarse-grained potential that accounts for interface coevolution based on joint multiple sequence alignments of two protein partners (co-alignments). However, even though such co-alignments could be built for none of the CAPRI targets in rounds 38 to 45, including host-pathogen and protein-oligosaccharide complexes and a redesigned interface, we identified multiple strategies that can be used to incorporate evolutionary constraints, which helped us to identify the most likely macromolecular binding modes. These strategies include template-based modeling where only local adjustments should be applied when query-template sequence identity is above 30% and larger perturbations are needed below this threshold; covariation-based structure prediction for individual protein partners; and the identification of evolutionarily conserved and structurally recurrent anchoring interface motifs. Overall, we submitted correct predictions among the top 5 models for 12 out of 19 interface challenges, including four High- and five Medium-quality predictions. Our top 20 models included correct predictions for three out of the five targets we missed in the top 5, including two targets for which misleading biological data led us to downgrade correct free docking models.

Performance and enhancement of the LZerD protein assembly pipeline in CAPRI 38-46.

We report the performance of the protein docking prediction pipeline of our group and the results for Critical Assessment of Prediction of Interactions (CAPRI) rounds 38-46. The pipeline integrates programs developed in our group as well as other existing scoring functions. The core of the pipeline is the LZerD protein-protein docking algorithm. If templates of the target complex are not found in PDB, the first step of our docking prediction pipeline is to run LZerD for a query protein pair. Meanwhile, in the case of human group prediction, we survey the literature to find information that can guide the modeling, such as protein-protein interface information. In addition to any literature information and binding residue prediction, generated docking decoys were selected by a rank aggregation of statistical scoring functions. The top 10 decoys were relaxed by a short molecular dynamics simulation before submission to remove atom clashes and improve side-chain conformations. In these CAPRI rounds, our group, particularly the LZerD server, showed robust performance. On the other hand, there are failed cases where some other groups were successful. To understand weaknesses of our pipeline, we analyzed sources of errors for failed targets. Since we noted that structure refinement is a step that needs improvement, we newly performed a comparative study of several refinement approaches. Finally, we show several examples that illustrate successful and unsuccessful cases by our group.

Template‐based modeling of diverse protein interactions in CAPRI rounds 38‐45

Structures of proteins complexed with other proteins, peptides, or ligands are essential for investigation of molecular mechanisms. However, the experimental structures of protein complexes of interest are often not available. Therefore, computational methods are widely used to predict these structures, and, of those methods, template-based modeling is the most successful. In the rounds 38-45 of the Critical Assessment of PRediction of Interactions (CAPRI), we applied template-based modeling for 9 of 11 protein-protein and protein-peptide interaction targets, resulting in medium and high-quality models for six targets. For the protein-oligosaccharide docking targets, we used constraints derived from template structures, and generated models of at least acceptable quality for most of the targets. Apparently, high flexibility of oligosaccharide molecules was the main cause preventing us from obtaining models of higher quality. We also participated in the CAPRI scoring challenge, the goal of which was to identify the highest quality models from a large pool of decoys. In this experiment, we tested VoroMQA, a scoring method based on interatomic contact areas. The results showed VoroMQA to be quite effective in scoring strongly binding and obligatory protein complexes, but less successful in the case of transient interactions. We extensively used manual intervention in both CAPRI modeling and scoring experiments. This oftentimes allowed us to select the correct templates from available alternatives and to limit the search space during the model scoring.

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations.

Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3–4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6–2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.

VoroMQA web server for assessing three-dimensional structures of proteins and protein complexes.

The VoroMQA (Voronoi tessellation-based Model Quality Assessment) web server is dedicated to the estimation of protein structure quality, a common step in selecting realistic and most accurate computational models and in validating experimental structures. As an input, the VoroMQA web server accepts one or more protein structures in PDB format. Input structures may be either monomeric proteins or multimeric protein complexes. For every input structure, the server provides both global and local (per-residue) scores. Visualization of the local scores along the protein chain is enhanced by providing secondary structure assignment and information on solvent accessibility. A unique feature of the VoroMQA server is the ability to directly assess protein-protein interaction interfaces. If this type of assessment is requested, the web server provides interface quality scores, interface energy estimates, and local scores for residues involved in inter-chain interfaces. VoroMQA, the underlying method of the web server, was extensively tested in recent community-wide CASP and CAPRI experiments. During these experiments VoroMQA showed outstanding performance both in model selection and in estimation of accuracy of local structural regions. The VoroMQA web server is available at http://bioinformatics.ibt.lt/wtsam/voromqa.

Modeling of protein complexes in CAPRI Round 37 using template-based approach combined with model selection.

We participated in Round 37 of the Critical Assessment of PRediction of Interactions (CAPRI), held jointly with the 12th edition of the Critical Assessment of protein Structure Prediction (CASP12), having two major objectives. First, we intended to test the utility of our PPI3D web server in finding and selecting templates for comparative modeling of structures of protein complexes. Our second aim was to evaluate the ability of our model accuracy estimation method VoroMQA to score and rank structural models for protein-protein interactions. Using sequence search in PPI3D and HHpred servers we identified multimeric templates for 7 of 11 CAPRI targets, and models of at least acceptable quality were constructed for 6 of them. The clustering and visual analysis features implemented in the PPI3D software were instrumental in detecting alternative protein-protein interaction interfaces among the identified templates. When a single binding mode was observed for homologous proteins, the structural modeling of the protein complex was fairly straightforward, whereas choosing the correct interaction template from several alternatives turned out to be a difficult task requiring manual intervention. The combination of full structure and interaction interface VoroMQA scores effectively ranked structural models of protein complexes and selected models of better quality from the CAPRI Scoring sets. The overall results show possible uses of PPI3D and VoroMQA in structural modeling of protein-protein interactions and suggest ways for further improvements of both methods.

Improved performance in CAPRI round 37 using LZerD docking and template-based modeling with combined scoring functions.

We report our group's performance for protein-protein complex structure prediction and scoring in Round 37 of the Critical Assessment of PRediction of Interactions (CAPRI), an objective assessment of protein-protein complex modeling. We demonstrated noticeable improvement in both prediction and scoring compared to previous rounds of CAPRI, with our human predictor group near the top of the rankings and our server scorer group at the top. This is the first time in CAPRI that a server has been the top scorer group. To predict protein-protein complex structures, we used both multi-chain template-based modeling (TBM) and our protein-protein docking program, LZerD. LZerD represents protein surfaces using 3D Zernike descriptors (3DZD), which are based on a mathematical series expansion of a 3D function. Because 3DZD are a soft representation of the protein surface, LZerD is tolerant to small conformational changes, making it well suited to docking unbound and TBM structures. The key to our improved performance in CAPRI Round 37 was to combine multi-chain TBM and docking. As opposed to our previous strategy of performing docking for all target complexes, we used TBM when multi-chain templates were available and docking otherwise. We also describe the combination of multiple scoring functions used by our server scorer group, which achieved the top rank for the scorer phase.

Lessons from (co-)evolution in the docking of proteins and peptides for CAPRI Rounds 28-35.

Computational protein-protein docking is of great importance for understanding protein interactions at the structural level. Critical assessment of prediction of interactions (CAPRI) experiments provide the protein docking community with a unique opportunity to blindly test methods based on real-life cases and help accelerate methodology development. For CAPRI Rounds 28-35, we used an automatic docking pipeline integrating the coarse-grained co-evolution-based potential InterEvScore. This score was developed to exploit the information contained in the multiple sequence alignments of binding partners and selectively recognize co-evolved interfaces. Together with Zdock/Frodock for rigid-body docking, SOAP-PP for atomic potential and Rosetta applications for structural refinement, this pipeline reached high performance on a majority of targets. For protein-peptide docking and interfacial water position predictions, we also explored different means of taking evolutionary information into account. Overall, our group ranked 1st by correctly predicting 10 targets, composed of 1 High, 7 Medium and 2 Acceptable predictions. Excellent and Outstanding levels of accuracy were reached for each of the two water prediction targets, respectively. Altogether, in 15 out of 18 targets in total, evolutionary information, either through co-evolution or conservation analyses, could provide key constraints to guide modeling towards the most likely assemblies. These results open promising perspectives regarding the way evolutionary information can be valuable to improve docking prediction accuracy. Proteins 2017; 85:378-390. © 2016 Wiley Periodicals, Inc.

Human and server docking prediction for CAPRI round 30-35 using LZerD with combined scoring functions.

We report the performance of protein-protein docking predictions by our group for recent rounds of the Critical Assessment of Prediction of Interactions (CAPRI), a community-wide assessment of state-of-the-art docking methods. Our prediction procedure uses a protein-protein docking program named LZerD developed in our group. LZerD represents a protein surface with 3D Zernike descriptors (3DZD), which are based on a mathematical series expansion of a 3D function. The appropriate soft representation of protein surface with 3DZD makes the method more tolerant to conformational change of proteins upon docking, which adds an advantage for unbound docking. Docking was guided by interface residue prediction performed with BindML and cons-PPISP as well as literature information when available. The generated docking models were ranked by a combination of scoring functions, including PRESCO, which evaluates the native-likeness of residues' spatial environments in structure models. First, we discuss the overall performance of our group in the CAPRI prediction rounds and investigate the reasons for unsuccessful cases. Then, we examine the performance of several knowledge-based scoring functions and their combinations for ranking docking models. It was found that the quality of a pool of docking models generated by LZerD, that is whether or not the pool includes near-native models, can be predicted by the correlation of multiple scores. Although the current analysis used docking models generated by LZerD, findings on scoring functions are expected to be universally applicable to other docking methods. Proteins 2017; 85:513-527. © 2016 Wiley Periodicals, Inc.