Flexible structural comparison allowing hinge-bending, swiveling motions

We present an efficient method for flexible comparison of protein structures, allowing swiveling motions. In all currently available methodologies developed and applied to the comparisons of protein structures, the molecules are considered to be rigid objects. The method described here extends and generalizes current approaches to searches for structural similarity between molecules by viewing proteins as objects consisting of rigid parts connected by rotary joints. During the matching, the rigid subparts are allowed to be rotated with respect to each other around swiveling points in one of the molecules. This technique straightforwardly detects structural motifs having hinge(s) between their domains. Whereas other existing methods detect hinge-bent motifs by initially finding the matching rigid parts and subsequently merging these together, our method automatically detects recurring substructures, allowing full 3 dimensional rotations about their swiveling points. Yet the method is extremely fast, avoiding the time-consuming full conformational space search. Comparison of two protein structures, without a predefinition of the motif, takes only seconds to one minute on a workstation per hinge. Hence, the molecule can be scanned for many potential hinge sites, allowing practically all C(alpha) atoms to be tried as swiveling points. This algorithm provides a highly efficient, fully automated tool. Its complexity is only O(n2), where n is the number of C(alpha) atoms in the compared molecules. As in our previous methodologies, the matching is independent of the order of the amino acids in the polypeptide chain. Here we illustrate the performance of this highly powerful tool on a large number of proteins exhibiting hinge-bending domain movements. Despite the motions, known hinge-bent domains/motifs which have been assembled and classified, are correctly identified. Additional matches are detected as well. This approach has been motivated by a technique for model based recognition of articulated objects originating in computer vision and robotics.     

Stability domains, substrate-induced conformational changes, and hinge-bending motions in a psychrophilic phosphoglycerate kinase. A microcalorimetric study

The cold-active phosphoglycerate kinase from the Antarctic bacterium Pseudomonas sp. TACII18 exhibits two distinct stability domains in the free, open conformation. It is shown that these stability domains do not match the structural N- and C-domains as the heat-stable domain corresponds to about 80 residues of the C-domain, including the nucleotide binding site, whereas the remaining of the protein contributes to the main heat-labile domain. This was demonstrated by spectroscopic and microcalorimetric analyses of the native enzyme, of its mutants, and of the isolated recombinant structural domains. It is proposed that the heat-stable domain provides a compact structure improving the binding affinity of the nucleotide, therefore increasing the catalytic efficiency at low temperatures. Upon substrate binding, the enzyme adopts a uniformly more stable closed conformation. Substrate-induced stability changes suggest that the free energy of ligand binding is converted into an increased conformational stability used to drive the hinge-bending motions and domain closure.     

Effects of ligand binding and oxidation on hinge-bending motions in S-adenosyl-L-homocysteine hydrolase

Domain motions of S-adenosyl-l-homocysteine (AdoHcy) hydrolase have been detected by time-resolved fluorescence anisotropy measurements. Time constants for reorientational motions in the native enzyme were compared with those for enzymes where key residues were altered by site-directed mutation. Mutations M351P, H353A, and P354A were selected in a hinge region for motion between the open and closed forms of the enzyme, as identified in a previous normal-mode study [Wang et al. (2005) Domain motions and the open-to-closed conformational transition of an enzyme: A normal-mode analysis of S-adenosyl-l-homocysteine hydrolase, Biochemistry 44, 7228-7239]. In wild-type, substrate-free AdoHcy hydrolase (NAD(+) cofactor in each subunit), reorientational motions were detected on time scales of 10-20 and 80-90 ns. The faster motion is attributed to the domain motion, and the slower motion is attributed to the tumbling of the enzyme. The domain motion was also detected for the enzyme complexes E(NADH/3'-keto-adenosine) and E(NAD(+)/3'-deoxyadenosine) but was absent for the complex E(NADH/3'-keto-neplanocin A). The results indicate that AdoHcy hydrolase exists in equilibrium of open and closed structures, with the equilibrium shifted toward the more mobile open form for the substrate-free enzyme, E(NAD(+)), and for intermediates formed early in the catalytic cycle after substrate binding or formed late prior to product release, E(NAD(+)/ligand). However, the strong inhibitor neplanocin A upon binding undergoes oxidation, forming the complex E(NADH/3'-keto-neplanocin). For this complex, which is analogous to the enzyme complex with the central catalytic intermediate, the equilibrium was shifted toward the more rigid closed form. A similar pattern was observed for M351P and P354A mutants. In contrast, the domain motion could not be detected, either in the absence or presence of ligands or with the cofactor in either the oxidized or reduced state, for the H353A protein, suggesting that this mutation changes the hinge-bending dynamics of the enzyme.     

Interdomain interactions in hinge-bending transitions.

BACKGROUND: The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions. RESULTS: We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface. CONCLUSIONS: Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.     

Engineering photocycle dynamics. Crystal structures and kinetics of three photoactive yellow protein hinge-bending mutants.

Crystallographic and spectroscopic analyses of three hinge-bending mutants of the photoactive yellow protein are described. Previous studies have identified Gly(47) and Gly(51) as possible hinge points in the structure of the protein, allowing backbone segments around the chromophore to undergo large concerted motions. We have designed, crystallized, and solved the structures of three mutants: G47S, G51S, and G47S/G51S. The protein dynamics of these mutants are significantly affected. Transitions in the photocycle, measured with laser induced transient absorption spectroscopy, show rates up to 6-fold different from the wild type protein and show an additive effect in the double mutant. Compared with the native structure, no significant conformational differences were observed in the structures of the mutant proteins. We conclude that the structural and dynamic integrity of the region around these mutations is of crucial importance to the photocycle and suggest that the hinge-bending properties of Gly(51) may also play a role in PAS domain proteins where it is one of the few conserved residues.     

Simulating protein motions with rigidity analysis.

Protein motions, ranging from molecular flexibility to large-scale conformational change, play an essential role in many biochemical processes. Despite the explosion in our knowledge of structural and functional data, our understanding of protein movement is still very limited. In previous work, we developed and validated a motion planning based method for mapping protein folding pathways from unstructured conformations to the native state. In this paper, we propose a novel method based on rigidity theory to sample conformation space more effectively, and we describe extensions of our framework to automate the process and to map transitions between specified conformations. Our results show that these additions both improve the accuracy of our maps and enable us to study a broader range of motions for larger proteins. For example, we show that rigidity-based sampling results in maps that capture subtle folding differences between protein G and its mutants, NuG1 and NuG2, and we illustrate how our technique can be used to study large-scale conformational changes in calmodulin, a 148 residue signaling protein known to undergo conformational changes when binding to Ca(2+). Finally, we announce our web-based protein folding server which includes a publicly available archive of protein motions: (http://parasol.tamu.edu/foldingserver/).     

Low-frequency motions in protein molecules. Beta-sheet and beta-barrel.

Low-frequency internal motions in protein molecules play a key role in biological functions. Based on previous work with alpha-helical structure, the quasi-continuum model is extended to the beta-structure, another fundamental element in protein molecules. In terms of the equations derived here, one can easily calculate the low-frequency wave number of a beta-sheet in an accordionlike motion, and the low-frequency wave number of a beta-barrel in a breathing motion. The calculated results for immunoglobulin G and concanavalin A agree well with the observations. These findings further verify that the observed low-frequency motion (or the so-called dominant low-frequency mode) in a protein molecule is essentially governed by the collective fluctuations of its weak bonds, especially hydrogen bonds, and the internal displacement of the massive atoms therein, as described by the quasi-continuum model.     

Motif-based searching in TOPS protein topology databases.

MOTIVATION: TOPS cartoons are a schematic ion of protein three-dimensional structures in two dimensions, and are used for understanding and manual comparison of protein folds. Recently, an algorithm that produces the cartoons automatically from protein structures has been devised and cartoons have been generated to represent all the structures in the structural databank. There is now a need to be able to define target topological patterns and to search the database for matching domains. RESULTS: We have devised a formal language for describing TOPS diagrams and patterns, and have designed an efficient algorithm to match a pattern to a set of diagrams. A pattern-matching system has been implemented, and tested on a database derived from all the current entries in the Protein Data Bank (15,000 domains). Users can search on patterns selected from a library of motifs or, alternatively, they can define their own search patterns. AVAILABILITY: The system is accessible over the Web at http://tops.ebi.ac.uk/tops     

Membrane localization and topology of a viral assembly protein.

The gene I protein (pI) of the filamentous bacteriophage f1 is required for the assembly of this virus. Antibodies specific to either the amino or carboxyl terminus of this protein were used to determine the location and topology of the gene I protein in f1-infected bacteria. pI is anchored in the inner membrane of Escherichia coli cells via a 20-amino-acid hydrophobic stretch, with its carboxyl-terminal 75 residues located in the periplasm and its amino-terminal 253 amino acids residing in the cytoplasm. By using the carboxyl-terminal pI antibody, a smaller protein, pI*, is also detected in f1-infected cells at a ratio of one to two molecules per molecule of pI. Analysis of proteins produced from a gene I amber mutant plasmid or bacteriophage suggests that pI* is most likely the result of an in-frame internal translational initiation event at methionine 241 of the 348-amino-acid pI. pI* is shown to be an integral inner membrane protein inserted in the same orientation as pI. The relation of the cellular locations of pI and pI* to some of the proposed functions of pI is discussed.     

Building blocks for plant gene assembly

The MultiSite Gateway cloning system, based on site-specific recombination, enables the assembly of multiple DNA fragments in predefined order, orientation, and frame register. To streamline the construction of recombinant genes for functional analysis in plants, we have built a collection of 36 reference Gateway entry clones carrying promoters, terminators, and reporter genes, as well as elements of the LhG4/LhGR two-component system. This collection obeys simple engineering rules. The genetic elements (parts) are designed in a standard format. They are interchangeable, fully documented, and can be combined at will according to the desired output. We also took advantage of the MultiSite Gateway recombination sites to create vectors in which two or three genes can be cloned simultaneously in separate expression cassettes. To illustrate the flexibility of these core resources for the construction of a wide variety of plant transformation vectors, we generated various transgenes encoding fluorescent proteins and tested their activity in plant cells. The structure and sequence of all described plasmids are accessible online at http://www.psb.ugent.be/gateway/. All accessions can be requested via the same Web site.     

The cellular prion protein: biochemistry,topology, and physiologic functions

Studies on the transmission from man to animals of Creutzfeld-Jacob disease (CJD) led Prusiner to identify a proteinaceous infectious particle lacking nucleic acid, which was called prion. The identification of the infectious prion (PrPsc) then led to the discovery of the normal cellular counterpart (PrPc). One of the still enigmatic aspects regarding prion diseases is actually how, where, and when the transformation PrPc/PrPsc is occurring, this being due to the result of a large extent to the fact that so far most studies have been dedicated to the formation and transmission of PrPsc, whereas the understanding of physiologic roles of PrPc are in their infancy. In this review, we hope to identify the most reliable hypotheses for future experiments on PrPc. This is relevant not only for the understanding of PrPc functions but also to unravel the enigmatic nature of PrPc/PrPsc conversion.     

Internal protein motions, concentrated glycerol, and hydrogen exchange studied in myoglobin

Experiments were carried out to measure the effect of concentrations of glycerol on H-exchange (HX) rates by using myoglobin as a test protein. Concentrated glycerol has only a small slowing effect on the HX kinetics of freely exposed amides, studied in a small molecule model (acetamide). Larger effects occur in structured proteins. The effect of solvent glycerol on different parts of the HX curve of myoglobin was studied by use of a selective "kinetic labeling" approach. Concentrated glycerol exerts an apparently reverse effect on protein H exchange; the faster exchanging "surface" protons are least affected, while the slower and slower amide NH is further slowed by larger and larger factors. These results seem inconsistent with solvent penetration models which generally visualize slower and slower protons as being placed, and undergoing exchange, farther and farther from the solvent-protein interface. On the other hand, the results are as expected for the local unfolding model for protein H exchange since concentrated glycerol is known to stabilize proteins against unfolding. In the local unfolding model, slower exchanging protons are released by way of higher energy and therefore generally larger, unfolding reactions. Larger unfoldings must be more inhibited by the glycerol effect.     

Structure of a hinge-bending bacteriophage T4 lysozyme mutant, Ile3-->Pro.

The mutant T4 phage lysozyme in which isoleucine 3 is replaced by proline (I3P) crystallizes in an orthorhombic form with two independent molecules in the asymmetric unit. Relative to wild-type lysozyme, which crystallizes in a trigonal form, the two I3P molecules undergo large hinge-bending displacements with the alignments of the amino-terminal and carboxy-terminal domains changed by 28.9 degrees and 32.9 degrees, respectively. The introduction of the mutation, together with the hinge-bending displacement, is associated with repacking of the side-chains of Phe4, Phe67 and Phe104. These aromatic residues are clustered close to the site of the mutation and are at the junction between the amino and carboxyl-terminal domains. As a result of this structural rearrangement the side-chain of Phe4 moves from a relatively solvent-exposed conformation to one that is largely buried. Mutant I3P also crystallizes in the same trigonal form as wild-type and, in this case, the observed structural changes are restricted to the immediate vicinity of the replacement. The main change is a shift of 0.3 to 0.5 A in the backbone of residues 1 to 5. The ability to crystallize I3P under similar conditions but in substantially different conformations suggests that the molecule undergoes large-scale hinge-bending displacements in solution. It is also likely that these conformational excursions are associated with repacking at the junction of the N-terminal and C-terminal domains. On the other hand, the analysis is complicated by possible effects of crystal packing. The different I3P crystal structures show substantial differences in the binding of solvent, both at the site of the Ile3-->Pro replacement and at other internal sites.     

GuHC1 induced unfolding-folding transition of a hinge-bending protein: horse muscle phosphoglycerate kinase

The antineoplastic compound N2-methyl-9-hydroxyellipticinium (9-OH-NME) is able to bind to different biological molecules after an oxidative activation by horseradish peroxidase and hydrogen peroxide. In this study, the efficient covalent binding in vitro of 9-OH-NME onto RNA and poly A is described. The phenomenon is analyzed by different HPLC methods and the yield of binding is determined using [3H]9-OH-NME. For an initial ratio drug per nucleotide of 0.07, the rb obtained (ratio of drug bound per nucleotide) of 0.026 for RNA and 0.044 for poly A, which represent respectively a yield of 40% and 60% for the drug fixation onto these macromolecules. These facts demonstrate the high electrophilicity of para-quinone-imine derivatives in ellipticinium series.     

Membrane topology of the Escherichia coli ExbD protein.

The ExbD protein is involved in the energy-coupled transport of ferric siderophores, vitamin B12, and B-group colicins across the outer membrane of Escherichia coli. In order to study ExbD membrane topology, ExbD-beta-lactamase fusion proteins were constructed. Cells expressing beta-lactamase fusions to residues 53, 57, 70, 76, 78, 80, 92, 121, and 134 of ExbD displayed high levels of ampicillin resistance, whereas fusions to residues 9 and 19 conferred no ampicillin resistance. It is concluded that the only hydrophobic segment of ExbD, encompassing residues 23 to 43, forms a transmembrane domain and that residues 1 to 22 are located in the cytoplasm and residues 44 to 141 are located in the periplasm.     

Control of protein topology at the endoplasmic reticulum.

I have described recent work that supports several conclusions that might not have been previously expected: first, that stop transfer, like the initiation of translocation, is receptor-mediated; second, that at least some of the topology-determining events at the ER membrane can be regulated (an example is provided where regulation may occur developmentally [PrP] and a possible example where receptor interactions for stop transfer seem to have been dissociated from those of integration in the membrane, in the course of evolution [apo B]); third, that these variations on the universal mechanism of eukaryotic secretory and transmembrane protein biogenesis can occur either through the variations in sequences presented to the common machinery of translocation or through variations in the machinery with which these sequences interact. Thus, on the one hand, at least some of these variations are directed by signal and stop transfer sequence subtypes and, on the other hand, in at least one case, a special cytoplasmic factor distinct from the core machinery for chain translocation also seems to be involved (RRL cytosolic factor effect on PrP topology) in the special handling of the STE stop transfer sequence subtype. In another case, the conserved universal machinery is engaged by a protein (apo B) to carry out an unusual, if not unique, mechanism presumably related to the lipid carrying role of this soluble secretory protein. Whether stop transfer sequence subtypes are involved here remains to be demonstrated, but it is a tempting hypothesis. Taken together, these findings suggest that the ER is more than a barrier to be overcome in protein export. In some cases, it plays a regulatory role in gene expression (e.g., alternate fates of PrP), and in other cases, it plays a role as a specialized assembly line for biogenesis of proteins with unusual properties. It seems likely that many other examples of proteins using these two mechanisms will be found, as well as entirely different variations on the mechanisms of protein biogenesis. A common conceptual theme is likely to be that they are all directed by discrete sequences within the particular newly synthesized proteins engaging both/either the common and/or distinctive component of the cellular machinery for protein biogenesis.     

Rapid membrane protein topology prediction

State-of-the-art methods for topology of α-helical membrane proteins are based on the use of time-consuming multiple sequence alignments obtained from PSI-BLAST or other sources. Here, we examine if it is possible to use the consensus of topology prediction methods that are based on single sequences to obtain a similar accuracy as the more accurate multiple sequence-based methods. Here, we show that TOPCONS-single performs better than any of the other topology prediction methods tested here, but ~6% worse than the best method that is utilizing multiple sequence alignments. AVAILABILITY AND IMPLEMENTATION: TOPCONS-single is available as a web server from http://single.topcons.net/ and is also included for local installation from the web site. In addition, consensus-based topology predictions for the entire international protein index (IPI) is available from the web server and will be updated at regular intervals.