OptGraft: A computational procedure for transferring a binding site onto an existing protein scaffold

One of the many challenging tasks of protein design is the introduction of a completely new function into an existing protein scaffold. In this study, we introduce a new computational procedure OptGraft for placing a novel binding pocket onto a protein structure so as its geometry is minimally perturbed. This is accomplished by introducing a two‐level procedure where we first identify where are the most appropriate locations to graft the new binding pocket into the protein fold by minimizing the departure from a set of geometric restraints using mixed‐integer linear optimization. On identifying the suitable locations that can accommodate the new binding pocket, CHARMM energy calculations are employed to identify what mutations in the neighboring residues, if any, are needed to ensure that the minimum energy conformation of the binding pocket conserves the desired geometry. This computational framework is benchmarked against the results available in the literature for engineering a copper binding site into thioredoxin protein. Subsequently, OptGraft is used to guide the transfer of a calcium‐binding pocket from thermitase protein (PDB: 1thm) into the first domain of CD2 protein (PDB:1hng). Experimental characterization of three de novo redesigned proteins with grafted calcium‐binding centers demonstrated that they all exhibit high affinities for terbium (K_d ～ 22, 38, and 55 μM) and can selectively bind calcium over magnesium.     

Caenorhabditis elegans NONO‐1: Insights into DBHS protein structure,architecture, and function

Members of the Drosophila behavior/human splicing (DBHS) protein family have been characterized in the vertebrates Homo sapiens and Mus musculus, and the invertebrates Drosophila melanogaster and Chironomus tentans. Collectively, both vertebrate and invertebrate DBHS proteins function throughout gene regulation, largely but not always, within the nucleus. In this study, we report a structural and bioinformatic analysis of the DBHS protein family to guide future studies into DBHS protein function. To explore the structural plasticity of the family, we describe the 2.4 Å crystal structure of Caenorhabditis elegans non‐POU domain‐containing octamer‐binding protein 1 (NONO‐1). The structure is dimeric, with a domain arrangement consistent with mammalian DBHS proteins. Comparison with the DBHS structures available from H. sapiens reveals that there is inherent domain flexibility within the homologous DBHS region. Mapping amino acid similarity within the family to the NONO‐1 dimer highlights the dimer interface, coiled‐coil oligomerization motif, and putative RNA binding surfaces. Surprisingly, the interior surface of RNA recognition motif 2 (RRM2) that faces a large internal void is highly variable, but the external β2–β3 loops of RRM2 show remarkable preservation. Overall, the DBHS region is under strong purifying selection, whereas the sequences N‐ and C‐terminal to the DBHS region are less constrained. The findings described in this study provide a molecular basis for further investigation into the mechanistic function of the DBHS protein family in biology.     

Ligand association rates to the inner-variable-domain of a dual-variable-domain immunoglobulin are significantly impacted by linker design

The DVD-Ig^TM protein is a dual-specific immunoglobulin. Each of the two arms of the molecule contains two variable domains, an inner variable domain and an outer variable domain linked in tandem, each with binding specificity for different targets or epitopes. One area of on-going research involves determining how the proximity of the outer variable domain affects the binding of ligands to the inner variable domain. To explore this area, we prepared a series of DVD-Ig proteins with binding specificities toward TNFα and an alternate therapeutic target. Kinetic measurements of TNFα binding to this series of DVD-Ig proteins were used to probe the effects of variable domain position and linker design on ligand on- and off-rates. We found that affinities for TNFα are generally lower when binding to the inner domain than to the outer domain and that this loss of affinity is primarily due to reduced association rate. This effect could be mitigated, to some degree, by linker design. We show several linker sequences that mitigate inner domain affinity losses in this series of DVD-Ig proteins. Moreover, we show that single chain proteolytic cleavage between the inner and outer domains, or complete outer domain removal, can largely restore inner domain TNFα affinity to that approaching the reference antibody. Taken together, these results suggest that a loss of affinity for inner variable domains in this set of DVD-Ig proteins may be largely driven by simple steric hindrance effects and can be reduced by careful linker design.     

Ligand association rates to the inner-variable-domain of a dual-variable-domain immunoglobulin are significantly impacted by linker design

The DVD-Ig™ protein is a dual-specific immunoglobulin. Each of the two arms of the molecule contains two variable domains, an inner variable domain and an outer variable domain linked in tandem, each with binding specificity for different targets or epitopes. One area of on-going research involves determining how the proximity of the outer variable domain affects the binding of ligands to the inner variable domain. To explore this area, we prepared a series of DVD-Ig proteins with binding specificities toward TNFα and an alternate therapeutic target. Kinetic measurements of TNFα binding to this series of DVD-Ig proteins were used to probe the effects of variable domain position and linker design on ligand on- and off-rates. We found that affinities for TNFα are generally lower when binding to the inner domain than to the outer domain and that this loss of affinity is primarily due to reduced association rate. This effect could be mitigated, to some degree, by linker design. We show several linker sequences that mitigate inner domain affinity losses in this series of DVD-Ig proteins. Moreover, we show that single chain proteolytic cleavage between the inner and outer domains, or complete outer domain removal, can largely restore inner domain TNFα affinity to that approaching the reference antibody. Taken together, these results suggest that a loss of affinity for inner variable domains in this set of DVD-Ig proteins may be largely driven by simple steric hindrance effects and can be reduced by careful linker design.Key words: dual variable domain immunoglobulin, DVD-Ig, immunotherapy, variable domain, antibody engineering, dual-specific, linker     

Score_set: A CAPRI benchmark for scoring protein complexes

Critical Assessment of PRedicted Interactions (CAPRI) has proven to be a catalyst for the development of docking algorithms. An essential step in docking is the scoring of predicted binding modes in order to identify stable complexes. In 2005, CAPRI introduced the scoring experiment, where upon completion of a prediction round, a larger set of models predicted by different groups and comprising both correct and incorrect binding modes, is made available to all participants for testing new scoring functions independently from docking calculations. Here we present an expanded benchmark data set for testing scoring functions, which comprises the consolidated ensemble of predicted complexes made available in the CAPRI scoring experiment since its inception. This consolidated scoring benchmark contains predicted complexes for 15 published CAPRI targets. These targets were subjected to 23 CAPRI assessments, due to existence of multiple binding modes for some targets. The benchmark contains more than 19,000 protein complexes. About 10% of the complexes represent docking predictions of acceptable quality or better, the remainder represent incorrect solutions (decoys). The benchmark set contains models predicted by 47 different predictor groups including web servers, which use different docking and scoring procedures, and is arguably as diverse as one may expect, representing the state of the art in protein docking. The data set is publicly available at the following URL: http://cb.iri.univ‐lille1.fr/Users/lensink/Score_set . Proteins 2014; 82:3163–3169. © 2014 Wiley Periodicals, Inc.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Mutational analysis of the KIX domain of CBP reveals residues critical for SREBP binding

Structure-based mutagenesis was used to probe the binding surface for the activation domain of sterol-responsive element binding protein (SREBP) in the KIX domain of CREB binding protein. A set of conserved residues scattering in the alpha2 helix and the extended C-terminal region of alpha 3 helix in the KIX domain including two arginines previously characterized as a hot spot for cofactor-mediated methylation was shown to be crucial for SREBP-KIX interaction, and was not essential for phosphorylated KID recognition. Therefore, our results suggest the existence of a SREBP binding site formed by positively charged residues in the C-terminal part of the extended alpha 3 helix of the KIX domain distinct from the previously identified phosphorylated KID binding site.     

Domain assignment for protein structures using a consensus approach: characterization and analysis.

A consensus approach for the assignment of structural domains in proteins is presented. The approach combines a number of previously published algorithms, and takes advantage of the elevated accuracy obtained when assignments from the individual algorithms are in agreement. The consensus approach is tested on a data set of 55 protein chains, for which domain assignments from four automated methods were known, and for which crystallographers assignments had been reported in the literature. Accuracy was found to increase in this test from 72% using individual algorithms to 100% when all four methods were in agreement. However a consensus prediction using all four methods was only possible for 52% of the dataset. The consensus approach [using three publicly available domain assignment algorithms (PUU, DETECTIVE, DOMAK)] was then used to make domain assignments for a data set of 787 protein chains from the Protein Data Bank. Analysis of the assignments showed 55.7% of assignments could be made automatically, and of these, 13.5% were multi-domain proteins. Of the remaining 44.3% that could not be assigned by the consensus procedure 90.4% had their domain boundaries assigned correctly by at least one of the algorithms. Once identified, these domains were analyzed for trends in their size and secondary structure class. In addition, the discontinuity of each domain along the protein chain was considered.     

A new, structurally nonredundant, diverse data set of protein-protein interfaces and its implications

Here, we present a diverse, structurally nonredundant data set of two-chain protein-protein interfaces derived from the PDB. Using a sequence order-independent structural comparison algorithm and hierarchical clustering, 3799 interface clusters are obtained. These yield 103 clusters with at least five nonhomologous members. We divide the clusters into three types. In Type I clusters, the global structures of the chains from which the interfaces are derived are also similar. This cluster type is expected because, in general, related proteins associate in similar ways. In Type II, the interfaces are similar; however, remarkably, the overall structures and functions of the chains are different. The functional spectrum is broad, from enzymes/inhibitors to immunoglobulins and toxins. The fact that structurally different monomers associate in similar ways, suggests "good" binding architectures. This observation extends a paradigm in protein science: It has been well known that proteins with similar structures may have different functions. Here, we show that it extends to interfaces. In Type III clusters, only one side of the interface is similar across the cluster. This structurally nonredundant data set provides rich data for studies of protein-protein interactions and recognition, cellular networks and drug design. In particular, it may be useful in addressing the difficult question of what are the favorable ways for proteins to interact. (The data set is available at http://protein3d.ncifcrf.gov/~keskino/ and http://home.ku.edu.tr/~okeskin/INTERFACE/INTERFACES.html.)     

Recursive domains in proteins.

The domain is a fundamental unit of protein structure. Numerous studies have analyzed folding patterns in protein domains of known structure to gain insight into the underlying protein folding process. Are such patterns a haphazard assortment or are they similar to sentences in a language, which can be generated by an underlying grammar? Specifically, can a small number of intuitively sensible rules generate a large class of folds, including feasible new folds? In this paper, we explore the extent to which four simple rules can generate the known all‐β folds, using tools from graph theory. As a control, an exhaustive set of β‐sandwiches was tested and found to be largely incompatible with such a grammar. The existence of a protein grammar has potential implications for both the mechanism of folding and the evolution of domains.     

A de novo designed protein protein interface

As an approach to both explore the physical/chemical parameters that drive molecular self-assembly and to generate novel protein oligomers, we have developed a procedure to generate protein dimers from monomeric proteins using computational protein docking and amino acid sequence design. A fast Fourier transform-based docking algorithm was used to generate a model for a dimeric version of the 56-amino-acid beta1 domain of streptococcal protein G. Computational amino acid sequence design of 24 residues at the dimer interface resulted in a heterodimer comprised of 12-fold and eightfold variants of the wild-type protein. The designed proteins were expressed, purified, and characterized using analytical ultracentrifugation and heteronuclear NMR techniques. Although the measured dissociation constant was modest ( approximately 300 microM), 2D-[(1)H,(15)N]-HSQC NMR spectra of one of the designed proteins in the absence and presence of its binding partner showed clear evidence of specific dimer formation.     

Highly covarying residues have a functional role in antibody constant domains

The ability to generate and design antibodies recognizing specific targets has revolutionized the pharmaceutical industry and medical imaging. Engineering antibody therapeutics in some cases requires modifying their constant domains to enable new and altered interactions. Engineering novel specificities into antibody constant domains has proved challenging due to the complexity of inter‐domain interactions. Covarying networks of residues that tend to cluster on the protein surface and near binding sites have been identified in some proteins. However, the underlying role these networks play in the protein resulting in their conservation remains unclear in most cases. Resolving their role is crucial, because residues in these networks are not viable design targets if their role is to maintain the fold of the protein. Conversely, these networks of residues are ideal candidates for manipulating specificity if they are primarily involved in binding, such as the myriad interdomain interactions maintained within antibodies. Here, we identify networks of evolutionarily‐related residues in C‐class antibody domains by evaluating covariation, a measure of propensity with which residue pairs vary dependently during evolution. We computationally test whether mutation of residues in these networks affects stability of the folded antibody domain, determining their viability as design candidates. We find that members of covarying networks cluster at domain‐domain interfaces, and that mutations to these residues are diverse and frequent during evolution, precluding their importance to domain stability. These results indicate that networks of covarying residues exist in antibody domains for functional reasons unrelated to thermodynamic stability, making them ideal targets for antibody design. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Solution structure of villin 14T, a domain conserved among actin-severing proteins.

The solution structure of the N-terminal domain of the actin-severing protein villin has been determined by multidimensional heteronuclear resonance spectroscopy. Villin is a member of a family of actin-severing proteins that regulate the organization of actin in the eukaryotic cytoskeleton. Members of this family are built from 3 or 6 homologous repeats of a structural domain of approximately 130 amino acids that is unrelated to any previously known structure. The N-terminal domain of villin (14T) contains a central beta-sheet with 4 antiparallel strands and a fifth parallel strand at one edge. This sheet is sandwiched between 2 helices on one side and a 2-stranded parallel beta-sheet with another helix on the other side. The strongly conserved sequence characteristic of the protein family corresponds to internal hydrophobic residues. Calcium titration experiments suggest that there are 2 binding sites for Ca2+, a stronger site near the N-terminal end of the longest helix, with a Kd of 1.8 +/- 0.4 mM, and a weaker site near the C-terminal end of the same helix, with a Kd of 11 +/- 2 mM. Mutational and biochemical studies of this domain in several members of the family suggest that the actin monomer binding site is near the parallel strand at the edge of the central beta-sheet.     

SLIDE, the Protein Interacting Domain of Imitation Switch Remodelers, Binds DDT-Domain Proteins of Different Subfamilies in Chromatin Remodeling Complexes

The Imitation Switch (ISWI) type adenosine triphosphate (ATP)-dependent chromatin remodeling factors are conserved proteins in eukaryotes, and some of them are known to form stable remodeling complexes...     

Prediction of protein conformational freedom from distance constraints

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints. Proteins 29:240–251, 1997. © 1997 Wiley-Liss, Inc.     

Survey of the geometric association of domain-domain interfaces

Considering the limited success of the most sophisticated docking methods available and the amount of computation required for systematic docking, cataloging all the known interfaces may be an alternative basis for the prediction of protein tertiary and quaternary structures. We classify domain interfaces according to the geometry of domain-domain association. By applying a simple and efficient method called "interface tag clustering," more than 4,000 distinct types of domain interfaces are collected from Protein Quaternary Structure Server and Protein Data Bank. Given a pair of interacting domains, we define "face" as the set of interacting residues in each single domain and the pair of interacting faces as an "interface." We investigate how the geometry of interfaces relates to a network of interacting protein families, such as how many different binding orientations are possible between two families or whether a family uses distinct surfaces or the same surface when the family has diverse interaction partners from various families. We show there are, on average, 1.2-1.9 different types of interfaces between interacting domains and a significant number of family pairs associate in multiple orientations. In general, a family tends to use distinct faces for each partner when the family has diverse interaction partners. Each face is highly specific to its interaction partner and the binding orientation. The relative positions of interface residues are generally well conserved within the same type of interface even between remote homologs. The classification result is available at http://www.biotec.tu-dresden.de/~wkim/supplement.     

Protein kinase CK2 interacts with a multi-protein binding domain of p53

p53 is one of the most powerful negative regulators of growth. To manage this in an efficient way it has to interact with a set of different cellular proteins. Most contacts with the cellular environment occur in the N- or the C-terminal domain of the protein. Since we previously found that p53 binds to the regulatory -subunit of CK2 we now analyzed N- and C-terminal domains of p53 separately for the binding of protein kinase CK2, an enzyme which seems to have a certain importance for proliferation processes. With different overlay assays we could map the binding domain of protein kinase CK2 to a sequence between amino acids 325-344, a region which coincides with the interaction domain of some other p53 binding proteins. We also found that the regulatory -subunit of protein kinase CK2 binds independent of the catalytic -subunit to this C-terminal domain of p53.     

Differential activities of cellular and viral macro domain proteins in binding of ADP-ribose metabolites

Macro domain is a highly conserved protein domain found in both eukaryotes and prokaryotes. Macro domains are also encoded by a set of positive-strand RNA viruses that replicate in the cytoplasm of animal cells, including coronaviruses and alphaviruses. The functions of the macro domain are poorly understood, but it has been suggested to be an ADP-ribose-binding module. We have here characterized three novel human macro domain proteins that were found to reside either in the cytoplasm and nucleus [macro domain protein 2 (MDO2) and ganglioside-induced differentiation-associated protein 2] or in mitochondria [macro domain protein 1 (MDO1)], and compared them with viral macro domains from Semliki Forest virus, hepatitis E virus, and severe acute respiratory syndrome coronavirus, and with a yeast macro protein, Poa1p. MDO2 specifically bound monomeric ADP-ribose with a high affinity (K_d = 0.15 μM), but did not bind poly(ADP-ribose) efficiently. MDO2 also hydrolyzed ADP-ribose-1″ phosphate, resembling Poa1p in all these properties. Ganglioside-induced differentiation-associated protein 2 did not show affinity for ADP-ribose or its derivatives, but instead bound poly(A). MDO1 was generally active in these reactions, including poly(A) binding. Individual point mutations in MDO1 abolished monomeric ADP-ribose binding, but not poly(ADP-ribose) binding; in poly(ADP-ribose) binding assays, the monomer did not compete against polymer binding. The viral macro proteins bound poly(ADP-ribose) and poly(A), but had a low affinity for monomeric ADP-ribose. Thus, the viral proteins do not closely resemble any of the human proteins in their biochemical functions. The differential activity profiles of the human proteins implicate them in different cellular pathways, some of which may involve RNA rather than ADP-ribose derivatives.     

Consensus design of a NOD receptor leucine rich repeat domain with binding affinity for a muramyl dipeptide,a bacterial cell wall fragment

Repeat proteins have recently emerged as especially well‐suited alternative binding scaffolds due to their modular architecture and biophysical properties. Here we present the design of a scaffold based on the consensus sequence of the leucine rich repeat (LRR) domain of the NOD family of cytoplasmic innate immune system receptors. Consensus sequence design has emerged as a protein design tool to create de novo proteins that capture sequence‐structure relationships and interactions present in nature. The multiple sequence alignment of 311 individual LRRs, which are the putative ligand‐recognition domain in NOD proteins, resulted in a consensus sequence protein containing two internal and N‐ and C‐capping repeats named CLRR2. CLRR2 protein is a stable, monomeric, and cysteine free scaffold that without any affinity maturation displays micromolar binding to muramyl dipeptide, a bacterial cell wall fragment. To our knowledge, this is the first report of direct interaction of a NOD LRR with a physiologically relevant ligand.