Search

Research has been conducted on the Anabaena apoflavodoxin. The folding reaction of this apoflavodoxin has been investigated via the site-directed mutagenesis and stopped-flow kinetics and the results indicate that a transient monomeric intermediate accumulates in the cause of the reaction.

Examines the folding reaction of apoflavodoxin by stopped-flow kinetics and site-directed mutagenesis. Accumulation of transient intermediates during folding reaction; Detection of intermediates at low urea concentration; Role of the intermediates in the folding reaction.

RNA-binding proteins are involved in many important regulatory processes in cells and their study is essential for a complete understanding of living organisms. They show a large variability from both structural and functional points of view. However, several recent studies performed on protein-RNA crystal structures have revealed interesting common properties. RNA-binding sites usually constitute patches of positively charged or polar residues that make most of the specific and non-specific contacts with RNA. Negatively charged or aliphatic residues are less frequent at protein-RNA interfaces, although they can also be found either forming aliphatic and positive-negative pairs in protein RNA-binding sites or contacting RNA through their main chains. Aromatic residues found within these interfaces are usually involved in specific base recognition at RNA single-strand regions. This specific recognition, in combination with structural complementarity, represents the key source for specificity in protein-RNA association. From all this knowledge, a variety of computational methods for prediction of RNA-binding sites have been developed based either on protein sequence or on protein structure. Some reported methods are really successful in the identification of RNA-binding proteins or the prediction of RNA-binding sites. Given the growing interest in the field, all these studies and prediction methods will undoubtedly contribute to the identification and comprehension of protein-RNA interactions.

Abstract Background Protein-protein interactions are involved in most cellular processes, and their detailed physico-chemical and structural characterization is needed in order to understand their function at the molecular level. In-silico docking tools can complement experimental techniques, providing three-dimensional structural models of such interactions at atomic resolution. In several recent studies, protein structures have been modeled as networks (or graphs), where the nodes represent residues and the connecting edges their interactions. From such networks, it is possible to calculate different topology-based values for each of the nodes, and to identify protein regions with high centrality scores, which are known to positively correlate with key functional residues, hot spots, and protein-protein interfaces. Results Here we show that this correlation can be efficiently used for the scoring of rigid-body docking poses. When integrated into the pyDock energy-based docking method, the new combined scoring function significantly improved the results of the individual components as shown on a standard docking benchmark. This improvement was particularly remarkable for specific protein complexes, depending on the shape, size, type, or flexibility of the proteins involved. Conclusions The network-based representation of protein structures can be used to identify protein-protein binding regions and to efficiently score docking poses, complementing energy-based approaches.

Abstract Background Protein-protein interactions are fundamental for the majority of cellular processes and their study is of enormous biotechnological and therapeutic interest. In recent years, a variety of computational approaches to the protein-protein docking problem have been reported, with encouraging results. Most of the currently available protein-protein docking algorithms are composed of two clearly defined parts: the sampling of the rotational and translational space of the interacting molecules, and the scoring and clustering of the resulting orientations. Although this kind of strategy has shown some of the most successful results in the CAPRI blind test http://www.ebi.ac.uk/msd-srv/capri, more efforts need to be applied. Thus, the sampling protocol should generate a pool of conformations that include a sufficient number of near-native ones, while the scoring function should discriminate between near-native and non-near-native proposed conformations. On the other hand, protocols to efficiently include full flexibility on the protein structures are increasingly needed. Results In these work we present new computational tools for protein-protein docking. We describe here the RotBUS (Rotation-Based Uniform Sampling) method to generate uniformly distributed sets of rigid-body docking poses, with a new fast calculation of the optimal contacting distance between molecules. We have tested the method on a standard benchmark of unbound structures and we can find near-native solutions in 100% of the cases. After applying a new fast filtering scheme based on residue-based desolvation, in combination with FTDock plus pyDock scoring, near-native solutions are found with rank ≤ 50 in 39% of the cases. Knowledge-based experimental restraints can be easily included to reduce computational times during sampling and improve success rates, and the method can be extended in the future to include flexibility of the side-chains. Conclusions This new sampling algorithm has the advantage of its high speed achieved by fast computing of the intermolecular distance based on a coarse representation of the interacting surfaces. In addition, a fast desolvation scoring permits the screening of millions of conformations at low computational cost, without compromising accuracy. The protocol presented here can be used as a framework to include restraints, flexibility and ensemble docking approaches.

Abstract Background Modelling proteins with multiple domains is one of the central challenges in Structural Biology. Although homology modelling has successfully been applied for prediction of protein structures, very often domain-domain interactions cannot be inferred from the structures of homologues and their prediction requires ab initio methods. Here we present a new structural prediction approach for modelling two-domain proteins based on rigid-body domain-domain docking. Results Here we focus on interacting domain pairs that are part of the same peptide chain and thus have an inter-domain peptide region (so called linker). We have developed a method called pyDockTET (tethered-docking), which uses rigid-body docking to generate domain-domain poses that are further scored by binding energy and a pseudo-energy term based on restraints derived from linker end-to-end distances. The method has been benchmarked on a set of 77 non-redundant pairs of domains with available X-ray structure. We have evaluated the docking method ZDOCK, which is able to generate acceptable domain-domain orientations in 51 out of the 77 cases. Among them, our method pyDockTET finds the correct assembly within the top 10 solutions in over 60% of the cases. As a further test, on a subset of 20 pairs where domains were built by homology modelling, ZDOCK generates acceptable orientations in 13 out of the 20 cases, among which the correct assembly is ranked lower than 10 in around 70% of the cases by our pyDockTET method. Conclusion Our results show that rigid-body docking approach plus energy scoring and linker-based restraints are useful for modelling domain-domain interactions. These positive results will encourage development of new methods for structural prediction of macromolecules with multiple (more than two) domains.

Abstract Background The study of protein-protein interactions is becoming increasingly important for biotechnological and therapeutic reasons. We can define two major areas therein: the structural prediction of protein-protein binding mode, and the identification of the relevant residues for the interaction (so called 'hot-spots'). These hot-spot residues have high interest since they are considered one of the possible ways of disrupting a protein-protein interaction. Unfortunately, large-scale experimental measurement of residue contribution to the binding energy, based on alanine-scanning experiments, is costly and thus data is fairly limited. Recent computational approaches for hot-spot prediction have been reported, but they usually require the structure of the complex. Results We have applied here normalized interface propensity (NIP) values derived from rigid-body docking with electrostatics and desolvation scoring for the prediction of interaction hot-spots. This parameter identifies hot-spot residues on interacting proteins with predictive rates that are comparable to other existing methods (up to 80% positive predictive value), and the advantage of not requiring any prior structural knowledge of the complex. Conclusion The NIP values derived from rigid-body docking can reliably identify a number of hot-spot residues whose contribution to the interaction arises from electrostatics and desolvation effects. Our method can propose residues to guide experiments in complexes of biological or therapeutic interest, even in cases with no available 3D structure of the complex.

Protein-protein interactions are at the heart of all cellular processes, with the ribosome emerging as a platform, orchestrating the nascent-chain interplay dynamics. Here, to study the characteristics governing co-translational protein folding and complex assembly, we combine selective ribosome profiling, imaging, and N-terminomics with all-atoms molecular dynamics. Focusing on conserved N-terminal acetyltransferases (NATs), we uncover diverging co-translational assembly pathways, where highly homologous subunits serve opposite functions. We find that only a few residues serve as "hotspots," initiating co-translational assembly interactions upon exposure at the ribosome exit tunnel. These hotspots are characterized by high binding energy, anchoring the entire interface assembly. Alpha-helices harboring hotspots are highly thermolabile, folding and unfolding during simulations, depending on their partner subunit to avoid misfolding. In vivo hotspot mutations disrupted co-translational complexation, leading to aggregation. Accordingly, conservation analysis reveals that missense NATs variants, causing neurodevelopmental and neurodegenerative diseases, disrupt putative hotspot clusters. Expanding our study to include phosphofructokinase, anthranilate synthase, and nucleoporin subcomplex, we employ AlphaFold-Multimer to model the complexes' complete structures. Computing MD-derived interface energy profiles, we find similar trends. Here, we propose a model based on the distribution of interface energy as a strong predictor of co-translational assembly., (© 2024. The Author(s).)

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem., (© 2023 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Reliably scoring and ranking candidate models of protein complexes and assigning their oligomeric state from the structure of the crystal lattice represent outstanding challenges. A community-wide effort was launched to tackle these challenges. The latest resources on protein complexes and interfaces were exploited to derive a benchmark dataset consisting of 1677 homodimer protein crystal structures, including a balanced mix of physiological and non-physiological complexes. The non-physiological complexes in the benchmark were selected to bury a similar or larger interface area than their physiological counterparts, making it more difficult for scoring functions to differentiate between them. Next, 252 functions for scoring protein-protein interfaces previously developed by 13 groups were collected and evaluated for their ability to discriminate between physiological and non-physiological complexes. A simple consensus score generated using the best performing score of each of the 13 groups, and a cross-validated Random Forest (RF) classifier were created. Both approaches showed excellent performance, with an area under the Receiver Operating Characteristic (ROC) curve of 0.93 and 0.94, respectively, outperforming individual scores developed by different groups. Additionally, AlphaFold2 engines recalled the physiological dimers with significantly higher accuracy than the non-physiological set, lending support to the reliability of our benchmark dataset annotations. Optimizing the combined power of interface scoring functions and evaluating it on challenging benchmark datasets appears to be a promising strategy., (© 2023 Wiley-VCH GmbH.)

Recycling of membrane proteins enables the reuse of receptors, ion channels and transporters. A key component of the recycling machinery is the endosomal sorting complex for promoting exit 1 (ESCPE-1), which rescues transmembrane proteins from the endolysosomal pathway for transport to the trans-Golgi network and the plasma membrane. This rescue entails the formation of recycling tubules through ESCPE-1 recruitment, cargo capture, coat assembly and membrane sculpting by mechanisms that remain largely unknown. Herein, we show that ESCPE-1 has a single-layer coat organization and suggest how synergistic interactions between ESCPE-1 protomers, phosphoinositides and cargo molecules result in a global arrangement of amphipathic helices to drive tubule formation. Our results thus define a key process of tubule-based endosomal sorting., (© 2023. The Author(s).)

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., (© 2021 Wiley Periodicals LLC.)

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., (© 2019 Wiley Periodicals, Inc.)

MAPK/ERK kinase (MEK) 1/2 are central signaling proteins that serve as specificity determinants of the MAPK/ERK cascade. More than twenty activating mutations have been reported for MEK1/2, and many of them are known to cause diseases such as cancers, arteriovenous malformation and RASopathies. Changes in their intrinsic activity do not seem to correlate with the severity of the diseases. Here we studied four MEK1/2 mutations using biochemical and molecular dynamic methods. Although the studied mutants elevated the activating phosphorylation of MEK they had no effect on the stimulated ERK1/2 phosphorylation. Studying the regulatory mechanism that may explain this lack of effect, we found that one type of mutation affects MEK stability and two types of mutations demonstrate a reduced sensitivity to PP2A. Together, our results indicate that some MEK mutations exert their function not only by their elevated intrinsic activity, but also by modulation of regulatory elements such as protein stability or dephosphorylation.

Metallocarboxypeptidase D (CPD) is a membrane-bound component of the trans-Golgi network that cycles to the cell surface through exocytic and endocytic pathways. Unlike other members of the metallocarboxypeptidase family, CPD is a multicatalytic enzyme with three carboxypeptidase-like domains, although only the first two domains are predicted to be enzymatically active. To investigate the enzymatic properties of each domain in human CPD, a critical active site Glu in domain I and/or II was mutated to Gln and the protein expressed, purified, and assayed with a wide variety of peptide substrates. CPD with all three domains intact displays >50% activity from pH 5.0 to 7.5 with a maximum at pH 6.5, as does CPD with mutation of domain I. In contrast, the domain II mutant displayed >50% activity from pH 6.5-7.5. CPD with mutations in both domains I and II was completely inactive towards all substrates and at all pH values. A quantitative peptidomics approach was used to compare the activities of CPD domains I and II towards a large number of peptides. CPD cleaved C-terminal Lys or Arg from a subset of the peptides. Most of the identified substrates of domain I contained C-terminal Arg, whereas comparable numbers of Lys- and Arg-containing peptides were substrates of domain II. We also report that some peptides with C-terminal basic residues were not cleaved by either domain I or II, showing the importance of the P1 position for CPD activity. Finally, the preference of domain I for C-terminal Arg was validated through molecular docking experiments. Together with the differences in pH optima, the different substrate specificities of CPD domains I and II allow the enzyme to perform distinct functions in the various locations within the cell.

We present an updated and integrated version of our widely used protein-protein docking and binding affinity benchmarks. The benchmarks consist of non-redundant, high-quality structures of protein-protein complexes along with the unbound structures of their components. Fifty-five new complexes were added to the docking benchmark, 35 of which have experimentally measured binding affinities. These updated docking and affinity benchmarks now contain 230 and 179 entries, respectively. In particular, the number of antibody-antigen complexes has increased significantly, by 67% and 74% in the docking and affinity benchmarks, respectively. We tested previously developed docking and affinity prediction algorithms on the new cases. Considering only the top 10 docking predictions per benchmark case, a prediction accuracy of 38% is achieved on all 55 cases and up to 50% for the 32 rigid-body cases only. Predicted affinity scores are found to correlate with experimental binding energies up to r=0.52 overall and r=0.72 for the rigid complexes., (Copyright © 2015. Published by Elsevier Ltd.)

Protein S functions as a cofactor for tissue factor pathway inhibitor (TFPI) and activated protein C (APC). The sex hormone binding globulin (SHBG)-like region of protein S, consisting of two laminin G-like domains (LG1 and LG2), contains the binding site for C4b-binding protein (C4BP) and TFPI. Furthermore, the LG-domains are essential for the TFPI-cofactor function and for expression of full APC-cofactor function. The aim of the current study was to localise functionally important interaction sites in the protein S LG-domains using amino acid substitutions. Four protein S variants were created in which clusters of surface-exposed amino acid residues within the LG-domains were substituted. All variants bound normally to C4BP and were fully functional as cofactors for APC in plasma and in pure component assays. Two variants, SHBG2 (E612A, I614A, F265A, V393A, H453A), involving residues from both LG-domains, and SHBG3 (K317A, I330A, V336A, D365A) where residues in LG1 were substituted, showed 50-60 % reduction in enhancement of TFPI in FXa inhibition assays. For SHBG3 the decreased TFPI cofactor function was confirmed in plasma based thrombin generation assays. Both SHBG variants bound to TFPI with decreased affinity in surface plasmon resonance experiments. The TFPI Kunitz 3 domain is known to contain the interaction site for protein S. Using in silico analysis and protein docking exercises, preliminary models of the protein S SHBG/TFPI Kunitz domain 3 complex were created. Based on a combination of experimental and in silico data we propose a binding site for TFPI on protein S, involving both LG-domains.

EMMPRIN/CD147 is mainly known for its protease inducing function but a role in promoting tumor angiogenesis has also been demonstrated. This study provides evidence that EMMPRIN is a new coreceptor for the VEGFR-2 tyrosine kinase receptor in both endothelial and tumor cells, as it directly interacts with it and regulates its activation by its VEGF ligand, signalling and functional consequences both in vitro and in vivo. Computational docking analyses and mutagenesis studies identified a molecular binding site in the extracellular domain of EMMPRIN located close to the cell membrane and containing the amino acids 195/199. EMMPRIN is overexpressed in cancer and hence is able to further potentiate VEGFR-2 activation, suggesting that a combinatory therapy of an antiangiogenic drug together with an inhibitor of EMMPRIN/VEGFR-2 interaction may have a greater impact on inhibiting angiogenesis and malignancy.

The tumor suppressor p53 regulates the expression of genes involved in cell cycle progression, senescence and apoptosis. Here, we investigated the effect of single point mutations in the oligomerization domain (OD) on tetramerization, transcription, ubiquitylation and stability of p53. As predicted by docking and molecular dynamics simulations, p53 OD mutants show functional defects on transcription, Mdm2-dependent ubiquitylation and 26S proteasome-mediated degradation. However, mutants unable to form tetramers are well degraded by the 20S proteasome. Unexpectedly, despite the lower structural stability compared to WT p53, p53 OD mutants form heterotetramers with WT p53 when expressed transiently or stably in cells wild type or null for p53. In consequence, p53 OD mutants interfere with the capacity of WT p53 tetramers to be properly ubiquitylated and result in changes of p53-dependent protein expression patterns, including the pro-apoptotic proteins Bax and PUMA under basal and adriamycin-induced conditions. Importantly, the patient derived p53 OD mutant L330R (OD1) showed the more severe changes in p53-dependent gene expression. Thus, in addition to the well-known effects on p53 stability, ubiquitylation defects promote changes in p53-dependent gene expression with implications on some of its functions., (Copyright © 2014 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

In addition to protein-protein docking, this CAPRI edition included new challenges, like protein-water and protein-sugar interactions, or the prediction of binding affinities and ΔΔG changes upon mutation. Regarding the standard protein-protein docking cases, our approach, mostly based on the pyDock scheme, submitted correct models as predictors and as scorers for 67% and 57% of the evaluated targets, respectively. In this edition, available information on known interface residues hardly made any difference for our predictions. In one of the targets, the inclusion of available experimental small-angle X-ray scattering (SAXS) data using our pyDockSAXS approach slightly improved the predictions. In addition to the standard protein-protein docking assessment, new challenges were proposed. One of the new problems was predicting the position of the interface water molecules, for which we submitted models with 20% and 43% of the water-mediated native contacts predicted as predictors and scorers, respectively. Another new problem was the prediction of protein-carbohydrate binding, where our submitted model was very close to being acceptable. A set of targets were related to the prediction of binding affinities, in which our pyDock scheme was able to discriminate between natural and designed complexes with area under the curve = 83%. It was also proposed to estimate the effect of point mutations on binding affinity. Our approach, based on machine learning methods, showed high rates of correctly classified mutations for all cases. The overall results were highly rewarding, and show that the field is ready to move forward and face new interesting challenges in interactomics., (Copyright © 2013 Wiley Periodicals, Inc.)

Community-wide blind prediction experiments such as CAPRI and CASP provide an objective measure of the current state of predictive methodology. Here we describe a community-wide assessment of methods to predict the effects of mutations on protein-protein interactions. Twenty-two groups predicted the effects of comprehensive saturation mutagenesis for two designed influenza hemagglutinin binders and the results were compared with experimental yeast display enrichment data obtained using deep sequencing. The most successful methods explicitly considered the effects of mutation on monomer stability in addition to binding affinity, carried out explicit side-chain sampling and backbone relaxation, evaluated packing, electrostatic, and solvation effects, and correctly identified around a third of the beneficial mutations. Much room for improvement remains for even the best techniques, and large-scale fitness landscapes should continue to provide an excellent test bed for continued evaluation of both existing and new prediction methodologies., (© 2013 Wiley Periodicals, Inc.)

Understanding and predicting the energetics of protein-protein interactions is fundamental to the structural modeling of protein complexes. Binding free energy can be approximated as a sum of pairwise atomic or residue contact energies, which are commonly inferred from contact frequencies observed in experimental protein structures. However, such statistically inferred potentials require certain assumptions and approximation. Here, we explore the possibility of deriving atomic and residue contact potentials directly from experimental binding free energy changes following mutation and present a number of such potentials. The first set of potentials is obtained by unweighted least-squares fitting and bootsrap aggregating. The second set is calculated using a weighting scheme optimized against absolute binding affinity data, so as to account for the over-representation of certain complexes, residues, and families of interactions. The congruence of the potentials with known physical chemistry is investigated. The potentials are further validated by ranking and clustering protein-protein docking poses.

Protein-protein interactions are central to almost all biological functions, and the atomic details of such interactions can yield insights into the mechanisms that underlie these functions. We present a web server that wraps and extends the SwarmDock flexible protein-protein docking algorithm. After uploading PDB files of the binding partners, the server generates low energy conformations and returns a ranked list of clustered docking poses and their corresponding structures. The user can perform full global docking, or focus on particular residues that are implicated in binding. The server is validated in the CAPRI blind docking experiment, against the most current docking benchmark, and against the ClusPro docking server, the highest performing server currently available.

Conformational fluctuations in proteins play key roles in their functions and interactions. In this work, validated conformational ensembles for ubiquitin have been used in docking trials. The ensembles were used in a systematic predictive study of known ubiquitin complexes by applying a cross-docking strategy against the bound structure of each partner. The global docking predictions obtained with the complete ubiquitin ensembles were significantly better than those obtained with the crystallographic structure of free ubiquitin. Importantly, in all cases we identified an individual ensemble member that performed equally well, or even better, than the bound structure of ubiquitin. These results unequivocally demonstrate that, for proteins that recognize binding partners by conformational selection, the availability of conformational ensembles can greatly improve the performance of automatic docking predictions. Our results highlight the need for docking methodologies to capitalize on validated ensemble representations of biomacromolecules.

Protein-protein interactions are responsible for the transfer of information inside the cell and represent one of the most interesting research fields in structural biology. Unfortunately, after decades of intense research, experimental approaches still have difficulties in providing 3D structures for the hundreds of thousands of interactions formed between the different proteins in a living organism. The use of theoretical approaches like docking aims to complement experimental efforts to represent the structure of the protein interactome. However, we cannot ignore that current methods have limitations due to problems of sampling of the protein-protein conformational space and the lack of accuracy of available force fields. Cases that are especially difficult for prediction are those in which complex formation implies a non-negligible change in the conformation of the interacting proteins, i.e., those cases where protein flexibility plays a key role in protein-protein docking. In this work, we present a new approach to treat flexibility in docking by global structural relaxation based on ultrafast discrete molecular dynamics. On a standard benchmark of protein complexes, the method provides a general improvement over the results obtained by rigid docking. The method is especially efficient in cases with large conformational changes upon binding, in which structure relaxation with discrete molecular dynamics leads to a predictive success rate double that obtained with state-of-the-art rigid-body docking.

Protein-protein interactions are fundamental for the majority of biological processes, so their structural, functional, and energetic characterization is of enormous biotechnological and therapeutic interest. In recent years, a variety of computational docking approaches to the structural prediction of protein-protein complexes have been reported, with encouraging results. However, a major bottleneck is found in cases with conformational movements upon binding, for which docking algorithms have to be extended beyond the rigid-body framework by introducing flexibility. Given the high computational cost of flexible docking, coarse-grained models offer an efficient alternative to full-atom descriptions. This work describes pyDockCG, a new coarse-grained potential for protein-protein docking scoring and refinement, based on the known UNRES model for polypeptide chains. The main novelty is the inclusion of two new terms accounting for the Coulomb electrostatics and the solvation energy. The latter has been devised by adapting the EEF1 model to the coarse-grained approach, with optimal parameters for protein-protein docking. The coarse-grained potential yielded highly similar values to the full-atom scoring function pyDock when applied to the rigid body docking sets, but at much lower computational cost. This efficiency makes it suitable for the treatment of flexibility during docking., (© 2011 American Chemical Society)

A detailed and complete structural knowledge of the interactome is one of the grand challenges in Biology, and a variety of computational docking approaches have been developed to complement experimental efforts and help in the characterization of protein-protein interactions. Among the different docking scoring methods, those based on physicochemical considerations can give the maximum accuracy at the atomic level, but they are usually computationally demanding and necessarily noisy when implemented in rigid-body approaches. Coarser-grained knowledge-based potentials are less sensitive to details of atomic arrangements, thus providing an efficient alternative for scoring of rigid-body docking poses. In this study, we have extracted new statistical potentials from intermolecular pairs of exposed residues in known complex structures, which were then used to score protein-protein docking poses. The new method, called SIPPER (scoring by intermolecular pairwise propensities of exposed residues), combines the value of residue desolvation based on solvent-exposed area with the propensity-based contribution of intermolecular residue pairs. This new scoring function found a near-native orientation within the top 10 predictions in nearly one-third of the cases of a standard docking benchmark and proved to be also useful as a filtering step, drastically reducing the number of docking candidates needed by energy-based methods like pyDock.

We describe here our results in the last CAPRI edition. We have participated in all targets, both as predictors and as scorers, using our pyDock docking methodology. The new challenges (homology-based modeling of the interacting subunits, domain-domain assembling, and protein-RNA interactions) have pushed our computer tools to the limits and have encouraged us to devise new docking approaches. Overall, the results have been quite successful, in line with previous editions, especially considering the high difficulty of some of the targets. Our docking approaches succeeded in five targets as predictors or as scorers (T29, T34, T35, T41, and T42). Moreover, with the inclusion of available information on the residues expected to be involved in the interaction, our protocol would have also succeeded in two additional cases (T32 and T40). In the remaining targets (except T37), results were equally poor for most of the groups. We submitted the best model (in ligand RMSD) among scorers for the unbound-bound target T29, the second best model among scorers for the protein-RNA target T34, and the only correct model among predictors for the domain assembly target T35. In summary, our excellent results for the new proposed challenges in this CAPRI edition showed the limitations and applicability of our approaches and encouraged us to continue developing methodologies for automated biomolecular docking., (© 2010 Wiley-Liss, Inc.)

Motivation: Prediction of protein-protein complexes from the coordinates of their unbound components usually starts by generating many potential predictions from a rigid-body 6D search followed by a second stage that aims to refine such predictions. Here, we present and evaluate a new method to effectively address the complexity and sampling requirements of the initial exhaustive search. In this approach we combine the projection of the interaction terms into 3D grid-based potentials with the efficiency of spherical harmonics approximations to accelerate the search. The binding energy upon complex formation is approximated as a correlation function composed of van der Waals, electrostatics and desolvation potential terms. The interaction-energy minima are identified by a novel, fast and exhaustive rotational docking search combined with a simple translational scanning. Results obtained on standard protein-protein benchmarks demonstrate its general applicability and robustness. The accuracy is comparable to that of existing state-of-the-art initial exhaustive rigid-body docking tools, but achieving superior efficiency. Moreover, a parallel version of the method performs the docking search in just a few minutes, opening new application opportunities in the current 'omics' world., Availability: http://sbg.cib.csic.es/Software/FRODOCK/

Plant immune responses often depend on leucine-rich repeat receptors that recognize microbe-associated molecular patterns or pathogen-specific virulence proteins, either directly or indirectly. When the recognition is direct, a molecular arms race takes place where plant receptors continually and rapidly evolve in response to virulence factor evolution. A useful model system to study ligand-receptor coevolution dynamics at the protein level is represented by the interaction between pathogen-derived polygalacturonases (PGs) and plant polygalacturonase-inhibiting proteins (PGIPs). We have applied codon substitution models to PGIP sequences of different eudicotyledonous families to identify putative positively selected sites and then compared these sites with the propensity of protein surface residues to interact with protein partners, based on desolvation energy calculations. The 2 approaches remarkably correlated in pinpointing several residues in the concave face of the leucine-rich repeat domain. These residues were mutated into alanine and their effect on the recognition of several PGs was tested, leading to the identification of unique hotspots for the PGIP-PG interaction. The combined approach used in this work can be of general utility in cases where structural information about a pattern-recognition receptor or resistance-gene product is available.

Desolvation property is used here to predict protein-protein binding sites exploiting the fact that lower-valued 'optimal docking area' ODA (Fernandez-Recio et al., 2005) points form cluster at the interface. The proposed method involves two steps; clustering the ODA points and representing ODA points by average ODA values. On 51 nonredundant proteins, results show the success rate improved considerably. Considering only significant ODA, the previous ODA method has obtained a success rate of 65% with overall success rate of 39%. The proposed method improved the overall success rate to 61%. Further, comparable results were found for X-ray and NMR structures.

The interactions between redox proteins are transient in nature. Therefore, very few crystal structures are available for the complexes formed between these proteins. Computational docking simulations thus provide a useful alternative method for studying the interactions between electron transfer proteins. In this paper, we have studied the interactions between the aa(3)-type cytochrome c oxidase of the cyanobacterium Phormidium laminosum and its redox partners plastocyanin and cytochrome c(6) using a combination of comparative modelling techniques and docking simulations. Rigid-body docking orientations were scored with a combined energy function that accounts for electrostatics and desolvation. These simulations have identified two plausible docking sites, one of which appears to be unique to the binding of plastocyanin to the oxidase. This unique binding site may be due to the presence of a long loop region in the subunit II of cyanobacterial oxidases. Control simulations were performed with the ba(3)-type cytochrome c oxidase and its redox partner cytochrome c(552) from Thermus thermophilus. The docking between cytochrome c oxidase and its redox partners plastocyanin and cytochrome c(6) is dominated by hydrophobic residues, a feature already observed from kinetic and structural studies in other complexes of P. laminosum (e.g. plastocyanin or cytochrome c(6) with cytochrome f and photosystem I).

Ferredoxin (Fd) interacts with ferredoxin-NADP(+) reductase (FNR) to transfer two electrons to the latter, one by one, which will finally be used to reduce NADP(+) to NADPH. The formation of a transient complex between Fd and FNR is required for the electron transfer (ET), and extensive mutational and crystallographic studies have been reported to characterize such protein-protein interaction. However, some aspects of the association mechanism still remain unclear. Moreover, in spite of their structural differences, flavodoxin (Fld) can replace Fd in its function and interact with FNR to transfer electrons with only slightly lower efficiency. Although crystallographic structures for the FNR:Fd association have been reported, experimental structural data for the FNR:Fld interaction are highly elusive. We have modeled here the interactions between FNR and both of its protein partners, Fd and Fld, using surface energy analysis, computational rigid-body docking simulations, and interface side-chain refinement. The results, consistent with previous experimental data, suggest the existence of alternative binding modes in these ET proteins.