Measurement of protein surface shape by solid angles

A method for measuring the gross shape of local regions of protein surface is presented and applied to an examination of three proteins. Any point on the protein surface can be assigned a number that measures the degree of convexity or concavity of the surface in the vicinity of the point. This number is computed by centring a sphere at that point and measuring how much of the sphere lies inside the protein. The sphere radius is a parameter chosen according to the scale of the features that are being analysed. The amount of sphere intersecting the protein is interpreted as a solid angle, denoted omega. Three proteins are analysed by this method: lysozyme, Superoxide dismutase and chymotrypsin. The resulting omega values are used to colour code the protein surfaces displayed on a colour raster graphics terminal. The method can be seen to reliably identify protrusions and depressions. The difficulty of developing a generally useful method for measuring protein surface shape is discussed. Possible applications of the solid-angle method include the analysis of shape complementarity at protein interfaces, the development of computer algorithms for predicting complexes between proteins, or between proteins and ligands, the identification of homologous epitopes on different proteins that might be immunologically crossreactive, and the determination of correlations between surface geometry and chemical properties.     

Protein docking using surface matching and supervised machine learning

Computational prediction of protein complex structures through docking offers a means to gain a mechanistic understanding of protein interactions that mediate biological processes. This is particularly important as the number of experimentally determined structures of isolated proteins exceeds the number of structures of complexes. A comprehensive docking procedure is described in which efficient sampling of conformations is achieved by matching surface normal vectors, fast filtering for shape complementarity, clustering by RMSD, and scoring the docked conformations using a supervised machine learning approach. Contacting residue pair frequencies, residue propensities, evolutionary conservation, and shape complementarity score for each docking conformation are used as input data to a Random Forest classifier. The performance of the Random Forest approach for selecting correctly docked conformations was assessed by cross-validation using a nonredundant benchmark set of X-ray structures for 93 heterodimer and 733 homodimer complexes. The single highest rank docking solution was the correct (near-native) structure for slightly more than one third of the complexes. Furthermore, the fraction of highly ranked correct structures was significantly higher than the overall fraction of correct structures, for almost all complexes. A detailed analysis of the difficult to predict complexes revealed that the majority of the homodimer cases were explained by incorrect oligomeric state annotation. Evolutionary conservation and shape complementarity score as well as both underrepresented and overrepresented residue types and residue pairs were found to make the largest contributions to the overall prediction accuracy. Finally, the method was also applied to docking unbound subunit structures from a previously published benchmark set.     

Use of pair potentials across protein interfaces in screening predicted docked complexes.

Empirical residue-residue pair potentials are used to screen possible complexes for protein-protein dockings. A correct docking is defined as a complex with not more than 2.5 A root-mean-square distance from the known experimental structure. The complexes were generated by "ftdock" (Gabb et al. J Mol Biol 1997;272:106-120) that ranks using shape complementarity. The complexes studied were 5 enzyme-inhibitors and 2 antibody-antigens, starting from the unbound crystallographic coordinates, with a further 2 antibody-antigens where the antibody was from the bound crystallographic complex. The pair potential functions tested were derived both from observed intramolecular pairings in a database of nonhomologous protein domains, and from observed intermolecular pairings across the interfaces in sets of nonhomologous heterodimers and homodimers. Out of various alternate strategies, we found the optimal method used a mole-fraction calculated random model from the intramolecular pairings. For all the systems, a correct docking was placed within the top 12% of the pair potential score ranked complexes. A combined strategy was developed that incorporated "multidock," a side-chain refinement algorithm (Jackson et al. J Mol Biol 1998;276:265-285). This placed a correct docking within the top 5 complexes for enzyme-inhibitor systems, and within the top 40 complexes for antibody-antigen systems.     

Protein-protein recognition analyzed by docking simulation.

Antibody-lysozyme and protease-inhibitor complexes are reconstituted by docking lysozyme as a rigid body onto the combining site of the antibodies and the inhibitors onto the active site of the proteases. Simplified protein models with one sphere per residue are subjected to simulated annealing using a crude energy function where the attractive component is proportional to the interface area. The procedure finds clusters of orientations in which a steric fit between the two protein components is achieved over a large contact surface. With five out of six complexes, the native structure of the complexes determined by X-ray crystallography is among those retained. Docked complexes are then subjected to conformational energy refinement with full atomic detail. With Fab HyHEL 5 and lysozyme, a native-like complex has the lowest refined energy. It can also be retrieved when starting with the X-ray structure of free lysozyme. However, some non-native complexes cannot be rejected: they form large interfaces, have a large number of H-bonds, and few unpaired polar groups. While these are necessary features of protein-protein recognition, they are not sufficient in determining specificity.     

A simple shape characteristic of protein-protein recognition

MOTIVATION: Observation of co-crystallized protein-protein complexes and low-resolution protein-protein docking studies suggest the existence of a binding-related anisotropic shape characteristic of protein-protein complexes. RESULTS: Our study systematically assessed the global shape of proteins in a non-redundant database of co-crystallized protein-protein complexes by measuring the distance of the surface residues to the protein's center of mass. The results show that on average the binding site residues are closer to the center of mass than the non-binding surface residues. Thus, the study directly detects an important and simple binding-related characteristic of protein shapes. The results provide an insight into one of the fundamental properties of protein structure and association.     

Reaction mechanism of glyoxalase I explored by an X-ray crystallographic analysis of the human enzyme in complex with a transition state analogue.

The structures of human glyoxalase I in complexes with S-(N-hydroxy-N-p-iodophenylcarbamoyl)glutathione (HIPC-GSH) and S-p-nitrobenzyloxycarbonylglutathione (NBC-GSH) have been determined at 2.0 and 1.72 A resolution, respectively. HIPC-GSH is a transition state analogue mimicking the enediolate intermediate that forms along the reaction pathway of glyoxalase I. In the structure, the hydroxycarbamoyl function is directly coordinated to the active site zinc ion. In contrast, the equivalent group in the NBC-GSH complex is approximately 6 A from the metal in a conformation that may resemble the product complex with S-D-lactoylglutathione. In this complex, two water molecules occupy the liganding positions at the zinc ion occupied by the hydroxycarbamoyl function in the enediolate analogue complex. Coordination of the transition state analogue to the metal enables a loop to close down over the active site, relative to its position in the product-like structure, allowing the glycine residue of the glutathione moiety to hydrogen bond with the protein. The structure of the complex with the enediolate analogue supports an "inner sphere mechanism" in which the GSH-methylglyoxal thiohemiacetal substrate is converted to product via a cis-enediolate intermediate. The zinc ion is envisioned to play an electrophilic role in catalysis by directly coordinating this intermediate. In addition, the carboxyl of Glu 172 is proposed to be displaced from the inner coordination sphere of the metal ion during substrate binding, thus allowing this group to facilitate proton transfer between the adjacent carbon atoms of the substrate. This proposal is supported by the observation that in the complex with the enediolate analogue the carboxyl group of Glu 172 is 3.3 A from the metal and is in an ideal position for reprotonation of the transition state intermediate. In contrast, Glu 172 is directly coordinated to the zinc ion in the complexes with S-benzylglutathione and with NBC-GSH.     

Plasmonic Spectral Engineering via Interferometric Illumination of Colloid Sphere Monolayers

A novel method is presented for complex structure fabrication, which is capable of breaking the hexagonal symmetry of the conventional colloid sphere lithography via the interferometric illumination of colloid sphere monolayers (IICSM). It is demonstrated that the perfect lateral synchronization of a linear intensity modulation originating from two-beam interference with respect to a hexagonal colloid sphere monolayer makes it possible to tune four complex structure parameters independently. Based on comparative study of hexagonal and rectangular hole doublet-arrays, which can be generated by linearly polarized light via homogeneous illumination and via IICSM, it is shown that the novel IICSM method enables plasmonic spectral engineering with higher degrees of freedom. The unique spectral properties of the complex patterns attainable via IICSM are more precisely tunable by properly selected azimuthal orientation during illumination and by the surrounding medium. It is shown that coupling phenomena between propagating and localized plasmonic modes on appropriately designed complex structures result in unique charge and near-field distribution accompanied by narrow Fano lines. Optimal configurations of complex plasmonic structures consisting of a rectangular array of hole doublets with different geometrical size parameters are presented, which ensure enhanced sensitivity in bio-detection.     

Multiplexed complexome profiling using tandem mass tags

Complexome profiling is a rapidly spreading, powerful technique to gain insight into the nature of protein complexes. It identifies and quantifies protein complexes separated into multiple fractions of increasing molecular mass using mass spectrometry-based, label-free bottom-up proteomics. Complexome profiling enables a sophisticated and thorough characterization of the composition, molecular mass, assembly, and interactions of protein complexes. However, in practice, its application is limited by the large number of samples it generates and the related time of mass spectrometry analyses. Here, we report an improved process workflow that implements tandem mass tags for multiplexing complexome profiling. This workflow substantially reduces the number of samples and measuring time without compromising protein identification or quantification reliability. In profiles from mitochondrial fractions of cells recovering from chloramphenicol treatment, tandem mass tags-multiplexed complexome profiling exhibited migration patterns of mature ATP synthase (complex V) and assembly intermediates that were consistent in composition and abundance with profiles obtained by the label-free approach. Reporter ion quantifications of proteins and complexes unaffected by the chloramphenicol treatment presented less variation in comparison to the label-free method. Incorporation of tandem mass tags enabled an efficient and robust complexome profiling analysis and may foster broader application for protein complex profiling in biomedical research and diagnostics.     

Cbl: many adaptations to regulate protein tyrosine kinases

Responses to extracellular stimuli are often transduced from cell-surface receptors to protein tyrosine kinases which, when activated, initiate the formation of protein complexes that transmit signals throughout the cell. A prominent component of these complexes is the product of the proto-oncogene c-Cbl, which specifically targets activated protein tyrosine kinases and regulates their signalling. How, then, does this multidomain protein shape the responses generated by these signalling complexes?     

Identifying foldable regions in protein sequence from the hydrophobic signal

Structural genomics initiatives aim to elucidate representative 3D structures for the majority of protein families over the next decade, but many obstacles must be overcome. The correct design of constructs is extremely important since many proteins will be too large or contain unstructured regions and will not be amenable to crystallization. It is therefore essential to identify regions in protein sequences that are likely to be suitable for structural study. Scooby-Domain is a fast and simple method to identify globular domains in protein sequences. Domains are compact units of protein structure and their correct delineation will aid structural elucidation through a divide-and-conquer approach. Scooby-Domain predictions are based on the observed lengths and hydrophobicities of domains from proteins with known tertiary structure. The prediction method employs an A*-search to identify sequence regions that form a globular structure and those that are unstructured. On a test set of 173 proteins with consensus CATH and SCOP domain definitions, Scooby-Domain has a sensitivity of 50% and an accuracy of 29%, which is better than current state-of-the-art methods. The method does not rely on homology searches and, therefore, can identify previously unknown domains.     

SIMS: computation of a smooth invariant molecular surface.

SIMS, a new method of calculating a smooth invariant molecular dot surface, is presented. The SIMS method generates the smooth molecular surface by rolling two probe spheres. A solvent probe sphere is rolled over the molecule and produces a Richards-Connolly molecular surface (MS), which envelops the solvent-excluded volume of the molecule. In deep crevices, Connolly's method of calculating the MS has two deficiencies. First, it produces self-intersecting parts of the molecular surface, which must be removed to obtain the correct MS. Second, the correct MS is not smooth, i.e., the direction of the normal vector of the MS is not continuous, and some points of the MS are singular. We present an exact method for removing self-intersecting parts and smoothing the singular regions of the MS. The singular MS is smoothed by rolling a smoothing probe sphere over the inward side of the singular MS. The MS in the vicinity of singularities is replaced with the reentrant surface of the smoothing probe sphere. The smoothing method does not disturb the topology of a singular MS, and the smooth MS is a better approximation of the dielectric border between high dielectric solvent and the low dielectric molecular interior. The SIMS method generates a smooth molecular dot surface, which has a quasi-uniform dot distribution in two orthogonal directions on the molecular surface, which is invariant with molecular rotation and stable under changes in the molecular conformation, and which can be used in a variety of implicit methods of modeling solvent effects. The SIMS program is faster than the Connolly MS program, and in a matter of seconds generates a smooth dot MS of a 200-residue protein. The program is available from the authors on request (see http:@femto.med.unc.edu/SIMS).     

Electron spin echo modulation and nuclear relaxation studies of staphylococcal nuclease and its metal-coordinating mutants

Electron spin echo envelope modulation (ESEEM) spectroscopy has been applied to the determination of the number of water molecules coordinated to the metal in the binary complex of staphylococcal nuclease with Mn2+, to the ternary enzyme-Mn2+-3',5'-pdTp complex, and to ternary complexes of a number of mutant enzymes in which metal-binding ligands have been individually altered. Quantitation of coordinated water is based on ESEEM spectral comparisons of Mn2+-EDTA and Mn2+-DTPA, which differ by a single inner sphere water, and with Mn2+-(H2O)6. It was found that Mn2+ in the ternary complex of the wild-type enzyme has a single additional coordinated water, as compared to Mn2+ in the binary complex, confirming earlier findings based on T1 measurements of bound water [Serpersu, E. H., Shortle, D. L., & Mildvan, A. S. (1987) Biochemistry 26, 1289-1300]. Ternary complexes of the mutant proteins D40E, D40G, and D21Y have the same number of water ligands as the ternary complex of the wild-type enzyme, while the D21E mutant has one less water ligand. In order to maintain octahedral coordination geometry, these findings require two additional ligands to Mn2+ from the protein in the binary complex of the wild-type enzyme, probably Glu 43 and Asp 19, and one additional ligand from the protein in the ternary D40G and D21E complexes. Other ESEEM studies of ternary Mn2+ complexes of wild-type, D21E, and D21Y mutants indicate the coordination by Mn2+ of a phosphate of 3',5'-pdTp, as demonstrated by a 31P contact interaction of 3.9 +/- 0.3 MHz. Magnetic interaction of Mn2+ with 31P could not be demonstrated with the D40G and D40E mutants, suggesting that metal-phosphate distances are greater in these mutants than in the wild-type protein. In a parallel NMR study of the paramagnetic effects of enzyme-bound Co2+ (which occupies the Mn2+ site on the enzyme) on the T1 of 31P from enzyme-bound 3',5'-pdTp and 5'-TMP, it was found that metal to 5'-phosphate distances are 0.9-1.6 A shorter in ternary complexes of the wild-type enzyme and of the D21E mutant than in ternary complexes of the D40G mutant. In all cases, the 5'-phosphate of pdTp is in the inner coordination sphere of Co2+ and the 3'-phosphate is predominantly in the second coordination sphere.     

Residue conservation information for generating near-native structures in protein-protein docking

MOTIVATION: Protein-protein docking algorithms typically generate large numbers of possible complex structures with only a few of them resembling the native structure. Recently (Duan et al., Protein Sci, 14:316-218, 2005), it was observed that the surface density of conserved residue positions is high at the interface regions of interacting protein surfaces, except for antibody-antigen complexes, where a lesser number of conserved positions than average is observed at the interface regions. Using this observation, we identified putative interacting regions on the surface of interacting partners and significantly improved docking results by assigning top ranks to near-native complex structures. In this paper, we combine the residue conservation information with a widely used shape complementarity algorithm to generate candidate complex structures with a higher percentage of near-native structures (hits). What is new in this work is that the conservation information is used early in the generation stage and not only in the ranking stage of the docking algorithm. This results in a significantly larger number of generated hits and an improved predictive ability in identifying the native structure of protein-protein complexes. RESULTS: We report on results from 48 well-characterized protein complexes, which have enough residue conservation information from the same 59 benchmark complexes used in our previous work. We compute conservation indices of residue positions on the surfaces of interacting proteins using available homologous sequences from UNIPROT and calculate the solvent accessible surface area. We combine this information with shape-complementarity scores to generate candidate protein-protein complex structures. When compared with pure shape-complementarity algorithms, performed by FTDock, our method results in significantly more hits, with the improvement being over 100% in many instances. We demonstrate that residue conservation information is useful not only in refinement and scoring of docking solutions, but also helpful in enrichment of near-native-structures during the generation of candidate geometries of complex structures.     

Coordination scheme and stereochemical configuration of manganese(II) adenosine 5'-diphosphate at the active site of 3-phosphoglycerate kinase

The structure of the MnIIADP complex at the active site of 3-phosphoglycerate kinase from yeast has been investigated by electron paramagnetic resonance (EPR) spectroscopy. Inhomogeneous broadening in the EPR signals for Mn(II) resulting from unresolved superhyperfine coupling to 17O regiospecifically incorporated into ADP shows that Mn(II) is coordinated to the alpha- and beta-phosphate groups of ADP at the active site of the enzyme. The EPR pattern for the enzyme-MnIIADP complex is characteristic of a predominantly axially symmetric zero-field splitting tensor. The symmetry and magnitude of the zero-field splitting interaction suggest that there is an additional negatively charged oxygen ligand in the coordination sphere of Mn(II). EPR measurements for solutions of the enzyme-MnIIADP complex in 17O-enriched water indicate that there are also two or three water molecules in the coordination sphere of the metal ion. EPR data for complexes with the two epimers of [alpha-17O]ADP have been used to determine the stereochemical configuration of the MnIIADP complex at the active site. EPR spectra for Mn(II) in the enzymic complex with (Rp)-[alpha-17O]ADP show an inhomogeneous broadening due to superhyperfine coupling with 17O whereas spectra for (Sp)-[alpha-17O]ADP complexes are indistinguishable from those for matched samples with unlabeled ADP. These results show that 3-phosphoglycerate kinase selectivity binds the alpha configuration of the alpha, beta chelate of MnIIADP. Addition of 3-phosphoglycerate to form the dead-end complex (enzyme-MnIIADP-3-phosphoglycerate) does not alter the EPR spectrum, but addition of vanadate to this complex causes marked changes in the spectral parameters.(ABSTRACT TRUNCATED AT 250 WORDS)     

Kinetics and mechanisms of the oxidation of myoglobin by Fe(III) and Cu(II) complexes

Two distinct mechanisms by which sperm whale myoglobin reduces, respectively, complexes of Fe(III) and Cu(II) and, in turn, is oxidized to metmyoglobin have been characterized. For both mechanisms, deoxymyoglobin is the active reductant. An outer sphere electron transfer, probably at the edge of the heme, is involved for Fe(III)NTA (NTA is nitrilotriacetic acid). This pathway does not involve ionic binding of the Fe(III) complex to the protein. The most reactive species of Fe(III)NTA is uncharged. No inhibition is observed with Ni(II) or Zn(II). An outer sphere site specific electron transfer is operative for reduction of Cu(II) complexes. The site has been characterized using NMR spectroscopy and involves one or more histidines. There is an initial binding of the Cu(II) chelate. The ternary complex of chelator-Cu(II)-deoxymyoglobin is a mandatory intermediate. Ni(II) and Zn(II) compete with Cu(II) for the binding site. A scheme for the participation of either or both of these mechanisms in reduction reactions of heme proteins is proposed. Both the overall redox potential, delta E0, and the stability constant for the ternary complex, K, govern the pathway and the reaction rate.     

Acetylation increases access of remodelling complexes to their nucleosome targets to enhance initiation of V(D)J recombination

Targeted chromatin remodelling is essential for many nuclear processes, including the regulation of V(D)J recombination. ATP-dependent nucleosome remodelling complexes are important players in this process whose activity must be tightly regulated. We show here that histone acetylation regulates nucleosome remodelling complex activity to boost RAG cutting during the initiation of V(D)J recombination. RAG cutting requires nucleosome mobilization from recombination signal sequences. Histone acetylation does not stimulate nucleosome mobilization per se by CHRAC, ACF or their catalytic subunit, ISWI. Instead, we find the more open structure of acetylated chromatin regulates the ability of nucleosome remodelling complexes to access their nucleosome templates. We also find that bromodomain/acetylated histone tail interactions can contribute to this targeting at limited concentrations of remodelling complex. We therefore propose that the changes in higher order chromatin structure associated with histone acetylation contribute to the correct targeting of nucleosome remodelling complexes and this is a novel way in which histone acetylation can modulate remodelling complex activity.     

Modeling of protein binary complexes using structural mass spectrometry data

In this article, we describe a general approach to modeling the structure of binary protein complexes using structural mass spectrometry data combined with molecular docking. In the first step, hydroxyl radical mediated oxidative protein footprinting is used to identify residues that experience conformational reorganization due to binding or participate in the binding interface. In the second step, a three-dimensional atomic structure of the complex is derived by computational modeling. Homology modeling approaches are used to define the structures of the individual proteins if footprinting detects significant conformational reorganization as a function of complex formation. A three-dimensional model of the complex is constructed from these binary partners using the ClusPro program, which is composed of docking, energy filtering, and clustering steps. Footprinting data are used to incorporate constraints-positive and/or negative-in the docking step and are also used to decide the type of energy filter-electrostatics or desolvation-in the successive energy-filtering step. By using this approach, we examine the structure of a number of binary complexes of monomeric actin and compare the results to crystallographic data. Based on docking alone, a number of competing models with widely varying structures are observed, one of which is likely to agree with crystallographic data. When the docking steps are guided by footprinting data, accurate models emerge as top scoring. We demonstrate this method with the actin/gelsolin segment-1 complex. We also provide a structural model for the actin/cofilin complex using this approach which does not have a crystal or NMR structure.     

Structures and mechanisms of formation of liprotides

Many proteins form complexes called liprotides with oleic acid and other cis-fatty acids under conditions where the protein is partially unfolded. The complexes vary in structure depending on the ratio of protein and lipid, but the most common structural organization is the core-shell structure, in which a layer of dynamic, partially unfolded and extended proteins surrounds a micelle-like fatty acid core. This structure, first reported for α-lactalbumin together with OA, resembles complexes formed between proteins and anionic surfactants like SDS. Liprotides first rose to fame through their anti-carcinogenic properties which still remains promising for topical applications though not yet implemented in the clinic. In addition, liprotides show potential in drug delivery thanks to the ability of the micelle core to solubilize and stabilize hydrophobic compounds, though applications are challenged by their sensitivity to acidic pH and dynamic exchange of lipids which makes them easy prey for serum “hoovers” such as albumin. However, liprotides are also of fundamental interest as a generic “protein complex structure”, demonstrating the many and varied structural consequences of protein-lipid interactions. Here we provide an overview of the different types of liprotide complexes, ranging from quasi-native complexes via core-shell structures to multi-layer structures, and discuss the many conditions under which they form. Given the many variable types of complexes that can form, rigorous biophysical analysis (stoichiometry, shape and structure of the complexes) remains crucial for a complete understanding of the mechanisms of action of this fascinating group of protein-lipid complexes both in vitro and in vivo.