BiGGER: A new (soft) docking algorithm for predicting protein interactions

Approaches to Backbone Flexibility in Protein-Protein Docking

Proteins are the fundamental organic macromolecules in living systems that play a key role in a variety of biological functions including immunological detection, intracellular trafficking, and signal transduction. The docking of proteins has greatly advanced during recent decades and has become a crucial complement to experimental methods. Protein-protein docking is a helpful method for simulating protein complexes whose structures have not yet been solved experimentally. This chapter focuses on major search tactics along with various docking programs used in protein-protein docking algorithms, which include: direct search, exhaustive global search, local shape feature matching, randomized search, and broad category of post-docking approaches. As backbone flexibility predictions and interactions in high-resolution protein-protein docking remain important issues in the overall optimization context, we have put forward several methods and solutions used to handle backbone flexibility. In addition, various docking methods that are utilized for flexible backbone docking, including ATTRACT, FlexDock, FLIPDock, HADDOCK, RosettaDock, FiberDock, etc., along with their scoring functions, algorithms, advantages, and limitations are discussed. Moreover, what progress in search technology is expected, including not only the creation of new search algorithms but also the enhancement of existing ones, has been debated. As conformational flexibility is one of the most crucial factors affecting docking success, more work should be put into evaluating the conformational flexibility upon binding for a particular case in addition to developing new algorithms to replace the rigid body docking and scoring approach.

Fine tuning rigid body docking results using the Dreiding force field: A computational study of 36 known nanobody-protein complexes

This paper aims to understand the binding strategies of a nanobody-protein pair by studying known complexes. Rigid body protein-ligand docking programs produce several complexes, called decoys, which are good candidates with high scores of shape complementarity, electrostatic interactions, desolvation, buried surface area, and Lennard-Jones potentials. However, the decoy that corresponds to the native structure is not known. We studied 36 nanobody-protein complexes from the single domain antibody database, sd-Ab DB, http://www.sdab-db.ca/. For each structure, a large number of decoys are generated using the Fast Fourier Transform algorithm of the software ZDOCK. The decoys were ranked according to their target protein-nanobody interaction energies, calculated by using the Dreiding Force Field, with rank 1 having the lowest interaction energy. Out of 36 protein data bank (PDB) structures, 25 true structures were predicted as rank 1. Eleven of the remaining structures required Ångstrom size rigid body translations of the nanobody relative to the protein to match the given PDB structure. After the translation, the Dreiding interaction (DI) energies of all complexes decreased and became rank 1. In one case, rigid body rotations as well as translations of the nanobody were required for matching the crystal structure. We used a Monte Carlo algorithm that randomly translates and rotates the nanobody of a decoy and calculates the DI energy. Results show that rigid body translations and the DI energy are sufficient for determining the correct binding location and pose of ZDOCK created decoys. A survey of the sd-Ab DB showed that each nanobody makes at least one salt bridge with its partner protein, indicating that salt bridge formation is an essential strategy in nanobody-protein recognition. Based on the analysis of the 36 crystal structures and evidence from existing literature, we propose a set of principles that could be used in the design of nanobodies.

The History of Desulfovibrio gigas Aldehyde Oxidoreductase-A Personal View

A story going back almost 40 years is presented in this manuscript. This is a different and more challenging way of reporting my research and I hope it will be useful to and target a wide-ranging audience. When preparing the manuscript and collecting references on the subject of this paper-aldehyde oxidoreductase from Desulfovibrio gigas-I felt like I was travelling back in time (and space), bringing together the people that have contributed most to this area of research. I sincerely hope that I can give my collaborators the credit they deserve. This study is not presented as a chronologic narrative but as a grouping of topics, the development of which occurred over many years.

Protein-protein interfaces in molecular glue-induced ternary complexes: classification, characterization, and prediction

Molecular glues are a class of small molecules that stabilize the interactions between proteins. Naturally occurring molecular glues are present in many areas of biology where they serve as central regulators of signaling pathways. Importantly, several clinical compounds act as molecular glue degraders that stabilize interactions between E3 ubiquitin ligases and target proteins, leading to their degradation. Molecular glues hold promise as a new generation of therapeutic agents, including those molecular glue degraders that can redirect the protein degradation machinery in a precise way. However, rational discovery of molecular glues is difficult in part due to the lack of understanding of the protein-protein interactions they stabilize. In this review, we summarize the structures of known molecular glue-induced ternary complexes and the interface properties. Detailed analysis shows different mechanisms of ternary structure formation. Additionally, we also review computational approaches for predicting protein-protein interfaces and highlight the promises and challenges. This information will ultimately help inform future approaches for rational molecular glue discovery.

A Guide to In Silico Drug Design

The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.

Classification and prediction of protein-protein interaction interface using machine learning algorithm

Structural insight of the protein-protein interaction (PPI) interface can provide knowledge about the kinetics, thermodynamics and molecular functions of the complex while elucidating its role in diseases and further enabling it as a potential therapeutic target. However, owing to experimental lag in solving protein-protein complex structures, three-dimensional (3D) knowledge of the PPI interfaces can be gained via computational approaches like molecular docking and post-docking analyses. Despite development of numerous docking tools and techniques, success in identification of native like interfaces based on docking score functions is limited. Hence, we employed an in-depth investigation of the structural features of the interface that might successfully delineate native complexes from non-native ones. We identify interface properties, which show statistically significant difference between native and non-native interfaces belonging to homo and hetero, protein-protein complexes. Utilizing these properties, a support vector machine (SVM) based classification scheme has been implemented to differentiate native and non-native like complexes generated using docking decoys. Benchmarking and comparative analyses suggest very good performance of our SVM classifiers. Further, protein interactions, which are proven via experimental findings but not resolved structurally, were subjected to this approach where 3D-models of the complexes were generated and most likely interfaces were predicted. A web server called Protein Complex Prediction by Interface Properties (PCPIP) is developed to predict whether interface of a given protein-protein dimer complex resembles known protein interfaces. The server is freely available at http://www.hpppi.iicb.res.in/pcpip/ .

Recent Advances in Molecular Docking for the Research and Discovery of Potential Marine Drugs

Marine drugs have long been used and exhibit unique advantages in clinical practices. Among the marine drugs that have been approved by the Food and Drug Administration (FDA), the protein–ligand interactions, such as cytarabine–DNA polymerase, vidarabine–adenylyl cyclase, and eribulin–tubulin complexes, are the important mechanisms of action for their efficacy. However, the complex and multi-targeted components in marine medicinal resources, their bio-active chemical basis, and mechanisms of action have posed huge challenges in the discovery and development of marine drugs so far, which need to be systematically investigated in-depth. Molecular docking could effectively predict the binding mode and binding energy of the protein–ligand complexes and has become a major method of computer-aided drug design (CADD), hence this powerful tool has been widely used in many aspects of the research on marine drugs. This review introduces the basic principles and software of the molecular docking and further summarizes the applications of this method in marine drug discovery and design, including the early virtual screening in the drug discovery stage, drug target discovery, potential mechanisms of action, and the prediction of drug metabolism. In addition, this review would also discuss and prospect the problems of molecular docking, in order to provide more theoretical basis for clinical practices and new marine drug research and development.

Prediction of Protein-Protein Interactions with Local Weight-Sharing Mechanism in Deep Learning

Protein-protein interactions (PPIs) are important for almost all cellular processes, including metabolic cycles, DNA transcription and replication, and signaling cascades. The experimental methods for identifying PPIs are always time-consuming and expensive. Therefore, it is important to develop computational approaches for predicting PPIs. In this paper, an improved model is proposed to use a machine learning method in the study of protein-protein interactions. With the consideration of the factors affecting the prediction of the PPIs, a method of feature extraction and fusion is proposed to improve the variety of the features to be considered in the prediction. Besides, with the consideration of the effect affected by the different input order of the two proteins, we propose a "Y-type" Bi-RNN model and train the network by using a method which both needs backward and forward training. In order to insure the training time caused on the extra training either a backward one or a forward one, this paper proposes a weight-sharing policy to minimize the parameters in the training. The experimental results show that the proposed method can achieve an accuracy of 99.57%, recall of 99.36%, sensitivity of 99.76%, precision of 99.74%, MCC of 99.14%, and AUC of 99.56% under the benchmark dataset.

Molecular characterization, computational analysis and expression profiling of Dmrt1 gene in Indian major carp, Labeo rohita (Hamilton 1822)

Sexual dimorphism of fish morphology, physiology and behavior is diverse and complex in nature. Doublesex and mab-3 related transcription factor (Dmrt) is a large protein family whose function is sexual development and differentiation in vertebrates. Here, we report a full-length cDNA sequence of Labeo rohita (rohu) Dmrt1 of 907 bp length having 798 bp of open reading frame encoding 265 amino acids. The molecular weight of rohu DMRT1 protein was found to be 28.74 KDa and isoelectric point was 7.53. DMRT1 protein contains 23 positively and 24 negatively charged amino acids with a GRAVY score of -0.618. A characteristic DM domain was found in DMRT1 protein, which is a novel DNA-binding domain. Phylogenetic analysis showed maximum similarity with Cyprinus carpio when compared with DMRT1 of other vertebrates. Molecular docking study identified active sites to be targeted for drug designing. Rohu DMRT1 was observed to interact with other proteins such as FOXL2, CYP19a1a, AMH and SOX9a. Differential expression study revealed higher expression in testis tissue implying its role in male sex differentiation and testicular development. The information generated in the present work could facilitate further research to resolve the issues related to gonadal maturation and reproduction of commercially important aquaculture species.

Pushing the accuracy limit of shape complementarity for protein-protein docking

Background

Protein-protein docking is a valuable computational approach for investigating protein-protein interactions. Shape complementarity is the most basic component of a scoring function and plays an important role in protein-protein docking. Despite significant progresses, shape representation remains an open question in the development of protein-protein docking algorithms, especially for grid-based docking approaches.

Results

We have proposed a new pairwise shape-based scoring function (LSC) for protein-protein docking which adopts an exponential form to take into account long-range interactions between protein atoms. The LSC scoring function was incorporated into our FFT-based docking program and evaluated for both bound and unbound docking on the protein docking benchmark 4.0. It was shown that our LSC achieved a significantly better performance than four other similar docking methods, ZDOCK 2.1, MolFit/G, GRAMM, and FTDock/G, in both success rate and number of hits. When considering the top 10 predictions, LSC obtained a success rate of 51.71% and 6.82% for bound and unbound docking, respectively, compared to 42.61% and 4.55% for the second-best program ZDOCK 2.1. LSC also yielded an average of 8.38 and 3.94 hits per complex in the top 1000 predictions for bound and unbound docking, respectively, followed by 6.38 and 2.96 hits for the second-best ZDOCK 2.1.

Conclusions

The present LSC method will not only provide an initial-stage docking approach for post-docking processes but also have a general implementation for accurate representation of other energy terms on grids in protein-protein docking. The software has been implemented in our HDOCK web server at http://hdock.phys.hust.edu.cn/.

RapidRMSD: rapid determination of RMSDs corresponding to motions of flexible molecules

Motivation

The root mean square deviation (RMSD) is one of the most used similarity criteria in structural biology and bioinformatics. Standard computation of the RMSD has a linear complexity with respect to the number of atoms in a molecule, making RMSD calculations time-consuming for the large-scale modeling applications, such as assessment of molecular docking predictions or clustering of spatially proximate molecular conformations. Previously, we introduced the RigidRMSD algorithm to compute the RMSD corresponding to the rigid-body motion of a molecule. In this study, we go beyond the limits of the rigid-body approximation by taking into account conformational flexibility of the molecule. We model the flexibility with a reduced set of collective motions computed with e.g. normal modes or principal component analysis.

Results

The initialization of our algorithm is linear in the number of atoms and all the subsequent evaluations of RMSD values between flexible molecular conformations depend only on the number of collective motions that are selected to model the flexibility. Therefore, our algorithm is much faster compared to the standard RMSD computation for large-scale modeling applications. We demonstrate the efficiency of our method on several clustering examples, including clustering of flexible docking results and molecular dynamics (MD) trajectories. We also demonstrate how to use the presented formalism to generate pseudo-random constant-RMSD structural molecular ensembles and how to use these in cross-docking.

Availability and implementation

We provide the algorithm written in C++ as the open-source RapidRMSD library governed by the BSD-compatible license, which is available at http://team.inria.fr/nano-d/software/RapidRMSD/. The constant-RMSD structural ensemble application and clustering of MD trajectories is available at http://team.inria.fr/nano-d/software/nolb-normal-modes/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Interaction between Neisseria gonorrhoeae bacterial peroxidase and its electron donor, the lipid-modified azurin

The Neisseria gonorrhoeae bacterial cytochrome c peroxidase plays a key role in detoxifying the cells from H₂ O₂ by reducing it to water using the lipid-modified azurin, LAz, a small type 1 copper protein, as electron donor. Here, the interaction between these two proteins was characterized by steady-state kinetics, two-dimensional NMR and molecular docking simulations. LAz is an efficient electron donor capable of activating this enzyme. This electron transfer complex is weak with a hydrophobic character, with LAz binding close to the electron transferring heme of the enzyme. The high catalytic rate (39 ± 0.03 s^-1 ) is explained by the LAz pre-orientation, due to a positive dipole moment, and by the fast-dynamic ensemble of orientations, suggested by the small chemical shifts.

The NMR2 Method to Determine Rapidly the Structure of the Binding Pocket of a Protein–Ligand Complex with High Accuracy

Structural characterization of complexes is crucial for a better understanding of biological processes and structure-based drug design. However, many protein–ligand structures are not solvable by X-ray crystallography, for example those with low affinity binders or dynamic binding sites. Such complexes are usually targeted by solution-state NMR spectroscopy. Unfortunately, structure calculation by NMR is very time consuming since all atoms in the complex need to be assigned to their respective chemical shifts. To circumvent this problem, we recently developed the Nuclear Magnetic Resonance Molecular Replacement (NMR2) method. NMR2 very quickly provides the complex structure of a binding pocket as measured by solution-state NMR. NMR2 circumvents the assignment of the protein by using previously determined structures and therefore speeds up the whole process from a couple of months to a couple of days. Here, we recall the main aspects of the method, show how to apply it, discuss its advantages over other methods and outline its limitations and future directions.

Topography of human cytochrome b5/cytochrome b5 reductase interacting domain and redox alterations upon complex formation

Cytochrome b₅ is the main electron acceptor of cytochrome b₅ reductase. The interacting domain between both human proteins has been unidentified up to date and very little is known about its redox properties modulation upon complex formation. In this article, we characterized the protein/protein interacting interface by solution NMR and molecular docking. In addition, upon complex formation, we measured an increase of cytochrome b₅ reductase flavin autofluorescence that was dependent upon the presence of cytochrome b_5. Data analysis of these results allowed us to calculate a dissociation constant value between proteins of 0.5±0.1μM and a 1:1 stoichiometry for the complex formation. In addition, a 30mV negative shift of cytochrome b₅ reductase redox potential in presence of cytochrome b₅ was also measured. These experiments suggest that the FAD group of cytochrome b₅ reductase increase its solvent exposition upon complex formation promoting an efficient electron transfer between the proteins.

Histone chaperone activity of Arabidopsis thaliana NRP1 is blocked by cytochrome c

Higher-order plants and mammals use similar mechanisms to repair and tolerate oxidative DNA damage. Most studies on the DNA repair process have focused on yeast and mammals, in which histone chaperone-mediated nucleosome disassembly/reassembly is essential for DNA to be accessible to repair machinery. However, little is known about the specific role and modulation of histone chaperones in the context of DNA damage in plants. Here, the histone chaperone NRP1, which is closely related to human SET/TAF-Iβ, was found to exhibit nucleosome assembly activity in vitro and to accumulate in the chromatin of Arabidopsis thaliana after DNA breaks. In addition, this work establishes that NRP1 binds to cytochrome c, thereby preventing the former from binding to histones. Since NRP1 interacts with cytochrome c at its earmuff domain, that is, its histone-binding domain, cytochrome c thus competes with core histones and hampers the activity of NRP1 as a histone chaperone. Altogether, the results obtained indicate that the underlying molecular mechanisms in nucleosome disassembly/reassembly are highly conserved throughout evolution, as inferred from the similar inhibition of plant NRP1 and human SET/TAF-Iβ by cytochrome c during DNA damage response.

Molecular Simulations of Disulfide-Rich Venom Peptides with Ion Channels and Membranes

Disulfide-rich peptides isolated from the venom of arthropods and marine animals are a rich source of potent and selective modulators of ion channels. This makes these peptides valuable lead molecules for the development of new drugs to treat neurological disorders. Consequently, much effort goes into understanding their mechanism of action. This paper presents an overview of how molecular simulations have been used to study the interactions of disulfide-rich venom peptides with ion channels and membranes. The review is focused on the use of docking, molecular dynamics simulations, and free energy calculations to (i) predict the structure of peptide-channel complexes; (ii) calculate binding free energies including the effect of peptide modifications; and (iii) study the membrane-binding properties of disulfide-rich venom peptides. The review concludes with a summary and outlook.

Software for molecular docking: a review

Molecular docking methodology explores the behavior of small molecules in the binding site of a target protein. As more protein structures are determined experimentally using X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy, molecular docking is increasingly used as a tool in drug discovery. Docking against homology-modeled targets also becomes possible for proteins whose structures are not known. With the docking strategies, the druggability of the compounds and their specificity against a particular target can be calculated for further lead optimization processes. Molecular docking programs perform a search algorithm in which the conformation of the ligand is evaluated recursively until the convergence to the minimum energy is reached. Finally, an affinity scoring function, ΔG [U total in kcal/mol], is employed to rank the candidate poses as the sum of the electrostatic and van der Waals energies. The driving forces for these specific interactions in biological systems aim toward complementarities between the shape and electrostatics of the binding site surfaces and the ligand or substrate.

Specificity and stability of transient protein-protein interactions

Remarkable features that are achieved in a protein-protein complex to precise levels are stability and specificity. Deviation from the normal levels of specificity and stability, which is often caused by mutations, could result in disease conditions. Chemical nature, 3-D arrangement and dynamics of interface residues code for both specificity and stability. This article reviews roles of interfacial residues in transient protein-protein complexes. It is proposed that aside from hotspot residues conferring stability to the complex, a small set of 'rigid' residues at the interface that maintain conformation between complexed and uncomplexed forms, play a major role in conferring specificity. Exceptionally, 'super hotspot' residues, which confer both stability and specificity, are attractive sites for interaction with small molecule inhibitors.

Molecular Docking Challenges and Limitations

Today, the development of new drugs is a challenging task of science. Researchers already applied molecular docking in the drug design field to simulate ligand- receptor interactions. Docking is a term used for computational schemes that attempt to find the “best” matching between two molecules in a complex formed from constituent molecules. It has a wide range of uses and applications in drug discovery. However, some defects still exist; the accuracy and speed of docking calculation is a challenge to explore and these methods can be enhanced as a solution to docking problem. The molecular docking problem can be defined as follows: Given the atomic coordinates of two molecules, predict their “correct” bound association. The chapter discusses common challenges critical aspects of docking method such as ligand- and receptor- conformation, flexibility and cavity detection, etc. It emphasis to the challenges and inadequacies with the theories behind as well as the examples.

Protein-Protein Docking

Protein-protein docking algorithms are powerful computational tools, capable of analyzing the protein-protein interactions at the atomic-level. In this chapter, we will review the theoretical concepts behind different protein-protein docking algorithms, highlighting their strengths as well as their limitations and pointing to important case studies for each method. The methods we intend to cover in this chapter include various search strategies and scoring techniques. This includes exhaustive global search, fast Fourier transform search, spherical Fourier transform-based search, direct search in Cartesian space, local shape feature matching, geometric hashing, genetic algorithm, randomized search, and Monte Carlo search. We will also discuss the different ways that have been used to incorporate protein flexibility within the docking procedure and some other future directions in this field, suggesting possible ways to improve the different methods.

Charge-Disproportionation Symmetry Breaking Creates a Heterodimeric Myoglobin Complex with Enhanced Affinity and Rapid Intracomplex Electron Transfer

We report rapid photoinitiated intracomplex electron transfer (ET) within a "charge-disproportionated" myoglobin (Mb) dimer with greatly enhanced affinity. Two mutually supportive Brownian Dynamics (BD) interface redesign strategies, one a new "heme-filtering" approach, were employed to "break the symmetry" of a Mb homodimer by pairing Mb constructs with complementary highly positive and highly negative net surface charges, introduced through D/E → K and K → E mutations, respectively. BD simulations using a previously developed positive mutant, Mb(+6) = Mb(D44K/D60K/E85K), led to construction of the complementary negative mutant Mb(-6) = Mb(K45E, K63E, K95E). Simulations predict the pair will form a well-defined complex comprising a tight ensemble of conformations with nearly parallel hemes, at a metal-metal distance ∼18-19 Å. Upon expression and X-ray characterization of the partners, BD predictions were verified through ET photocycle measurements enabled by Zn-deuteroporphyrin substitution, forming the [ZnMb(-6), Fe(3+)Mb(+6)] complex. Triplet ET quenching shows charge disproportionation increases the binding constant by no less than ∼5 orders of magnitude relative to wild-type Mb values. All progress curves for charge separation (CS) and charge recombination (CR) are reproduced by a generalized kinetic model for the interprotein ET photocycle. The intracomplex ET rate constants for both CS and CR are increased by over 5 orders of magnitude, and their viscosity independence is indicative of true interprotein ET, rather than dynamic gating as seen in previous studies. The complex displays an unprecedented timecourse for CR of the CS intermediate I. After a laser flash, I forms through photoinduced CS, accumulates to a maximum concentration, then dies away through CR. However, before completely disappearing, I reappears without another flash and reaches a second maximum before disappearing completely.

Predicting Protein-Protein Interactions Using BiGGER: Case Studies

The importance of understanding interactomes makes preeminent the study of protein interactions and protein complexes. Traditionally, protein interactions have been elucidated by experimental methods or, with lower impact, by simulation with protein docking algorithms. This article describes features and applications of the BiGGER docking algorithm, which stands at the interface of these two approaches. BiGGER is a user-friendly docking algorithm that was specifically designed to incorporate experimental data at different stages of the simulation, to either guide the search for correct structures or help evaluate the results, in order to combine the reliability of hard data with the convenience of simulations. Herein, the applications of BiGGER are described by illustrative applications divided in three Case Studies: (Case Study A) in which no specific contact data is available; (Case Study B) when different experimental data (e.g., site-directed mutagenesis, properties of the complex, NMR chemical shift perturbation mapping, electron tunneling) on one of the partners is available; and (Case Study C) when experimental data are available for both interacting surfaces, which are used during the search and/or evaluation stage of the docking. This algorithm has been extensively used, evidencing its usefulness in a wide range of different biological research fields.

Molecular Docking Challenges and Limitations

Today, the development of new drugs is a challenging task of science. Researchers already applied molecular docking in the drug design field to simulate ligand- receptor interactions. Docking is a term used for computational schemes that attempt to find the “best” matching between two molecules in a complex formed from constituent molecules. It has a wide range of uses and applications in drug discovery. However, some defects still exist; the accuracy and speed of docking calculation is a challenge to explore and these methods can be enhanced as a solution to docking problem. The molecular docking problem can be defined as follows: Given the atomic coordinates of two molecules, predict their “correct” bound association. The chapter discusses common challenges critical aspects of docking method such as ligand- and receptor- conformation, flexibility and cavity detection, etc. It emphasis to the challenges and inadequacies with the theories behind as well as the examples.

Surfing the Protein-Protein Interaction Surface Using Docking Methods: Application to the Design of PPI Inhibitors

Blocking protein-protein interactions (PPI) using small molecules or peptides modulates biochemical pathways and has therapeutic significance. PPI inhibition for designing drug-like molecules is a new area that has been explored extensively during the last decade. Considering the number of available PPI inhibitor databases and the limited number of 3D structures available for proteins, docking and scoring methods play a major role in designing PPI inhibitors as well as stabilizers. Docking methods are used in the design of PPI inhibitors at several stages of finding a lead compound, including modeling the protein complex, screening for hot spots on the protein-protein interaction interface and screening small molecules or peptides that bind to the PPI interface. There are three major challenges to the use of docking on the relatively flat surfaces of PPI. In this review we will provide some examples of the use of docking in PPI inhibitor design as well as its limitations. The combination of experimental and docking methods with improved scoring function has thus far resulted in few success stories of PPI inhibitors for therapeutic purposes. Docking algorithms used for PPI are in the early stages, however, and as more data are available docking will become a highly promising area in the design of PPI inhibitors or stabilizers.

Protein docking with predicted constraints

This paper presents a constraint-based method for improving protein docking results. Efficient constraint propagation cuts over 95% of the search time for finding the configurations with the largest contact surface, provided a contact is specified between two amino acid residues. This makes it possible to scan a large number of potentially correct constraints, lowering the requirements for useful contact predictions. While other approaches are very dependent on accurate contact predictions, ours requires only that at least one correct contact be retained in a set of, for example, one hundred constraints to test. It is this feature that makes it feasible to use readily available sequence data to predict specific potential contacts. Although such prediction is too inaccurate for most purposes, we demonstrate with a Naïve Bayes Classifier that it is accurate enough to more than double the average number of acceptable models retained during the crucial filtering stage of protein docking when combined with our constrained docking algorithm. All software developed in this work is freely available as part of the Open Chemera Library.

Structural and functional analysis of novel human cytochrome C targets in apoptosis

Since the first description of apoptosis four decades ago, great efforts have been made to elucidate, both in vivo and in vitro, the molecular mechanisms involved in its regulation. Although the role of cytochrome c during apoptosis is well established, relatively little is known about its participation in signaling pathways in vivo due to its essential role during respiration. To obtain a better understanding of the role of cytochrome c in the onset of apoptosis, we used a proteomic approach based on affinity chromatography with cytochrome c as bait in this study. In this approach, novel cytochrome c interaction partners were identified whose in vivo interaction and cellular localization were facilitated through bimolecular fluorescence complementation. Modeling of the complex interface between cytochrome c and its counterparts indicated the involvement of the surface surrounding the heme crevice of cytochrome c, in agreement with the vast majority of known redox adducts of cytochrome c. However, in contrast to the high turnover rate of the mitochondrial cytochrome c redox adducts, those occurring under apoptosis led to the formation of stable nucleo-cytoplasmic ensembles, as inferred mainly from surface plasmon resonance and nuclear magnetic resonance measurements, which permitted us to corroborate the formation of such complexes in vitro. The results obtained suggest that human cytochrome c interacts with pro-survival, anti-apoptotic proteins following its release into the cytoplasm. Thus, cytochrome c may interfere with cell survival pathways and unlock apoptosis in order to prevent the spatial and temporal coexistence of antagonist signals.

A knowledge-based scoring function for protein-RNA interactions derived from a statistical mechanics-based iterative method

Protein-RNA interactions play important roles in many biological processes. Given the high cost and technique difficulties in experimental methods, computationally predicting the binding complexes from individual protein and RNA structures is pressingly needed, in which a reliable scoring function is one of the critical components. Here, we have developed a knowledge-based scoring function, referred to as ITScore-PR, for protein-RNA binding mode prediction by using a statistical mechanics-based iterative method. The pairwise distance-dependent atomic interaction potentials of ITScore-PR were derived from experimentally determined protein–RNA complex structures. For validation, we have compared ITScore-PR with 10 other scoring methods on four diverse test sets. For bound docking, ITScore-PR achieved a success rate of up to 86% if the top prediction was considered and up to 94% if the top 10 predictions were considered, respectively. For truly unbound docking, the respective success rates of ITScore-PR were up to 24 and 46%. ITScore-PR can be used stand-alone or easily implemented in other docking programs for protein–RNA recognition.

Intuitive, but not simple: including explicit water molecules in protein-protein docking simulations improves model quality

Characterizing the nature of interaction between proteins that have not been experimentally cocrystallized requires a computational docking approach that can successfully predict the spatial conformation adopted in the complex. In this work, the Hydropathic INTeractions (HINT) force field model was used for scoring docked models in a data set of 30 high-resolution crystallographically characterized "dry" protein-protein complexes and was shown to reliably identify native-like models. However, most current protein-protein docking algorithms fail to explicitly account for water molecules involved in bridging interactions that mediate and stabilize the association of the protein partners, so we used HINT to illuminate the physical and chemical properties of bridging waters and account for their energetic stabilizing contributions. The HINT water Relevance metric identified the "truly" bridging waters at the 30 protein-protein interfaces and we utilized them in "solvated" docking by manually inserting them into the input files for the rigid body ZDOCK program. By accounting for these interfacial waters, a statistically significant improvement of ∼24% in the average hit-count within the top-10 predictions the protein-protein dataset was seen, compared to standard "dry" docking. The results also show scoring improvement, with medium and high accuracy models ranking much better than incorrect ones. These improvements can be attributed to the physical presence of water molecules that alter surface properties and better represent native shape and hydropathic complementarity between interacting partners, with concomitantly more accurate native-like structure predictions.

How to use not-always-reliable binding site information in protein-protein docking prediction

In many protein-protein docking algorithms, binding site information is used to help predicting the protein complex structures. Using correct and accurate binding site information can increase protein-protein docking success rate significantly. On the other hand, using wrong binding sites information should lead to a failed prediction, or, at least decrease the success rate. Recently, various successful theoretical methods have been proposed to predict the binding sites of proteins. However, the predicted binding site information is not always reliable, sometimes wrong binding site information could be given. Hence there is a high risk to use the predicted binding site information in current docking algorithms. In this paper, a softly restricting method (SRM) is developed to solve this problem. By utilizing predicted binding site information in a proper way, the SRM algorithm is sensitive to the correct binding site information but insensitive to wrong information, which decreases the risk of using predicted binding site information. This SRM is tested on benchmark 3.0 using purely predicted binding site information. The result shows that when the predicted information is correct, SRM increases the success rate significantly; however, even if the predicted information is completely wrong, SRM only decreases success rate slightly, which indicates that the SRM is suitable for utilizing predicted binding site information.

The scoring of poses in protein-protein docking: current capabilities and future directions

BACKGROUND

Protein-protein docking, which aims to predict the structure of a protein-protein complex from its unbound components, remains an unresolved challenge in structural bioinformatics. An important step is the ranking of docked poses using a scoring function, for which many methods have been developed. There is a need to explore the differences and commonalities of these methods with each other, as well as with functions developed in the fields of molecular dynamics and homology modelling.

RESULTS

We present an evaluation of 115 scoring functions on an unbound docking decoy benchmark covering 118 complexes for which a near-native solution can be found, yielding top 10 success rates of up to 58%. Hierarchical clustering is performed, so as to group together functions which identify near-natives in similar subsets of complexes. Three set theoretic approaches are used to identify pairs of scoring functions capable of correctly scoring different complexes. This shows that functions in different clusters capture different aspects of binding and are likely to work together synergistically.

CONCLUSIONS

All functions designed specifically for docking perform well, indicating that functions are transferable between sampling methods. We also identify promising methods from the field of homology modelling. Further, differential success rates by docking difficulty and solution quality suggest a need for flexibility-dependent scoring. Investigating pairs of scoring functions, the set theoretic measures identify known scoring strategies as well as a number of novel approaches, indicating promising augmentations of traditional scoring methods. Such augmentation and parameter combination strategies are discussed in the context of the learning-to-rank paradigm.

TagDock: an efficient rigid body docking algorithm for oligomeric protein complex model construction and experiment planning

We report here new computational tools and strategies to efficiently generate three-dimensional models for oligomeric biomolecular complexes in cases where there is limited experimental restraint data to guide the docking calculations. Our computational tools are designed to rapidly and exhaustively enumerate all geometrically possible docking poses for an oligomeric complex, rather than generate detailed, atomic-resolution models. Experimental data, such as interatomic distance measurements, are then used to select and refine docking poses that are consistent with the experimental restraints. Our computational toolkit is designed for use with sparse data sets to generate intermediate-resolution docking models, and utilizes distance difference matrix analysis to identify further restraint measurements that will provide maximum additional structural refinement. Thus, these tools can be used to help plan optimal residue positions for probe incorporation in labor-intensive biophysical experiments such as chemical cross-linking, electron paramagnetic resonance, or Förster resonance energy transfer spectroscopy studies. We present benchmark results for docking the collection of all 176 heterodimer protein complexes from the ZDOCK database, as well as a protein homodimer with recently collected experimental distance restraints, to illustrate the toolkit's capabilities and performance, and to demonstrate how distance difference matrix analysis can automatically identify and prioritize additional restraint measurements that allow us to rapidly optimize docking poses.

Using chemical shift perturbation to characterise ligand binding

Chemical shift perturbation (CSP, chemical shift mapping or complexation-induced changes in chemical shift, CIS) follows changes in the chemical shifts of a protein when a ligand is added, and uses these to determine the location of the binding site, the affinity of the ligand, and/or possibly the structure of the complex. A key factor in determining the appearance of spectra during a titration is the exchange rate between free and bound, or more specifically the off-rate koff. When koff is greater than the chemical shift difference between free and bound, which typically equates to an affinity Kd weaker than about 3μM, then exchange is fast on the chemical shift timescale. Under these circumstances, the observed shift is the population-weighted average of free and bound, which allows Kd to be determined from measurement of peak positions, provided the measurements are made appropriately. (1)H shifts are influenced to a large extent by through-space interactions, whereas (13)Cα and (13)Cβ shifts are influenced more by through-bond effects. (15)N and (13)C' shifts are influenced both by through-bond and by through-space (hydrogen bonding) interactions. For determining the location of a bound ligand on the basis of shift change, the most appropriate method is therefore usually to measure (15)N HSQC spectra, calculate the geometrical distance moved by the peak, weighting (15)N shifts by a factor of about 0.14 compared to (1)H shifts, and select those residues for which the weighted shift change is larger than the standard deviation of the shift for all residues. Other methods are discussed, in particular the measurement of (13)CH3 signals. Slow to intermediate exchange rates lead to line broadening, and make Kd values very difficult to obtain. There is no good way to distinguish changes in chemical shift due to direct binding of the ligand from changes in chemical shift due to allosteric change. Ligand binding at multiple sites can often be characterised, by simultaneous fitting of many measured shift changes, or more simply by adding substoichiometric amounts of ligand. The chemical shift changes can be used as restraints for docking ligand onto protein. By use of quantitative calculations of ligand-induced chemical shift changes, it is becoming possible to determine not just the position but also the orientation of ligands.

Antimalarial activity of cupredoxins: the interaction of Plasmodium merozoite surface protein 119 (MSP119) and rusticyanin

The discovery of effective new antimalarial agents is urgently needed. One of the most frequently studied molecules anchored to the parasite surface is the merozoite surface protein-1 (MSP1). At red blood cell invasion MSP1 is proteolytically processed, and the 19-kDa C-terminal fragment (MSP119) remains on the surface and is taken into the red blood cell, where it is transferred to the food vacuole and persists until the end of the intracellular cycle. Because a number of specific antibodies inhibit erythrocyte invasion and parasite growth, MSP119 is therefore a promising target against malaria. Given the structural homology of cupredoxins with the Fab domain of monoclonal antibodies, an approach combining NMR and isothermal titration calorimetry (ITC) measurements with docking calculations based on BiGGER is employed on MSP119-cupredoxin complexes. Among the cupredoxins tested, rusticyanin forms a well defined complex with MSP119 at a site that overlaps with the surface recognized by the inhibitory antibodies. The addition of holo-rusticyanin to infected cells results in parasitemia inhibition, but negligible effects on parasite growth can be observed for apo-rusticyanin and other proteins of the cupredoxin family. These findings point to rusticyanin as an excellent therapeutic tool for malaria treatment and provide valuable information for drug design.