Structure of the TRAIL-DR5 complex reveals mechanisms conferring specificity in apoptotic initiation.

TRAIL, an apoptosis inducing ligand, has at least four cell surface receptors including the death receptor DR5. Here we report the crystal structure at 2.2 A resolution of a complex between TRAIL and the extracellular region of DR5. TRAIL forms a central homotrimer around which three DR5 molecules bind. Radical differences in the surface charge of the ligand, together with variation in the alignment of the two receptor domains confer specificity between members of these ligand and receptor families. The existence of a switch mechanism allowing variation in receptor domain alignment may mean that it is possible to engineer receptors with multiple specificities by exploiting contact positions unique to individual receptor-ligand pairs.     

Multiple octamer binding sites in the promoter region of the bovine alpha s2-casein gene.

Using a set of overlapping oligonucleotides from the promoter region of the bovine alpha s2-casein gene we have identified two nuclear factors which probably are involved in expression of this gene and the related calcium sensitive alpha s1- and beta-casein genes. One of these factors which was present in extracts of all tissues that have been tested including Hela cells turned out to be the octamer binding protein OCT-1. Oct-1 binds with different affinity to 4 sites at positions centred around -480, -260, -210 and -50. The strongest of these 4 binding sites, the one around position -50, is highly conserved in all calcium sensitive caseins of mouse, rat, rabbit and cattle. The other nuclear factor (MGF, mammary gland factor) which is specifically expressed in the mammary gland, binds to a site around position -90. This binding site is also highly conserved in all calcium sensitive caseins of mouse, rat, rabbit and cattle.     

A quantitative analysis of interfacial amino acid conservation in protein-protein hetero complexes

A long-standing question in molecular biology is whether interfaces of protein-protein complexes are more conserved than the rest of the protein surfaces. Although it has been reported that conservation can be used as an indicator for predicting interaction sites on proteins, there are recent reports stating that the interface regions are only slightly more conserved than the rest of the protein surfaces, with conservation signals not being statistically significant enough for predicting protein-protein binding sites. In order to properly address these controversial reports we have studied a set of 28 well resolved hetero complex structures of proteins that consists of transient and non-transient complexes. The surface positions were classified into four conservation classes and the conservation index of the surface positions was quantitatively analyzed. The results indicate that the surface density of highly conserved positions is significantly higher in the protein-protein interface regions compared with the other regions of the protein surface. However, the average conservation index of the patches in the interface region is not significantly higher compared with other surface regions of the protein structures. This finding demonstrates that the number of conserved residue positions is a more appropriate indicator for predicting protein-protein binding sites than the average conservation index in the interacting region. We have further validated our findings on a set of 59 benchmark complex structures. Furthermore, an analysis of 19 complexes of antigen-antibody interactions shows that there is no conservation of amino acid positions in the interacting regions of these complexes, as expected, with the variable region of the immunoglobulins interacting mostly with the antigens. Interestingly, antigen interacting regions also have a higher number of non-conserved residue positions in the interacting region than the rest of the protein surface.     

Prediction of the odorant binding site of olfactory receptor proteins by human-mouse comparisons

Olfactory receptors (ORs) are a large family of proteins involved in the recognition and discrimination of numerous odorants. These receptors belong to the G-protein coupled receptor (GPCR) hyperfamily, for which little structural data are available. In this study we predict the binding site residues of OR proteins by analyzing a set of 1441 OR protein sequences from mouse and human. The central insight utilized is that functional contact residues would be conserved among pairs of orthologous receptors, but considerably less conserved among paralogous pairs. Using judiciously selected subsets of 218 ortholog pairs and 518 paralog pairs, we have identified 22 sequence positions that are both highly conserved among the putative orthologs and variable among paralogs. These residues are disposed on transmembrane helices 2 to 7, and on the second extracellular loop of the receptor. Strikingly, although the prediction makes no assumption about the location of the binding site, these amino acid positions are clustered around a pocket in a structural homology model of ORs, mostly facing the inner lumen. We propose that the identified positions constitute the odorant binding site. This conclusion is supported by the observation that all but one of the predicted binding site residues correspond to ligand-contact positions in other rhodopsin-like GPCRs.     

Centripetal modules and ancient introns

We have created an algorithm which instantiates the centripetal definition of modules, compact regions of protein structure, as introduced by Go and Nosaka (M. Go and M. Nosaka, 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Bio. 52, 915-924). That definition seeks the minima of a function that sums the squares of C-alpha carbon distances over a window around each amino acid residue in a three-dimensional protein structure and identifies such minima with module boundaries. We analyze a set of 44 ancient conserved proteins, with known three-dimensional structures, which have intronless homologues in bacteria and intron-containing homologues in the eukaryotes, with a corresponding set of 988 intron positions. We show that the phase zero intron positions are significantly correlated with the module boundaries (p = 0.0002), while the intron positions that lie within codons, in phase one and phase two, are not correlated with these 'centripetal' module boundaries. Furthermore, we analyze the phylogenetic distribution of intron positions and identify a subset of putatively 'ancient' intron positions: phase zero positions in one phylogenetic kingdom which have an associated intron either in an identical position or within three codons in another phylogenetic kingdom (a notion of intron sliding). This subset of 120 'ancient' introns lies closer to the module boundaries than does the full set of phase zero introns with high significance, a p-value of 0.008. We conclude that the behavior of this set of introns supports the prediction of a mixed theory: that some introns are very old and were used for exon shuffling in the progenote, while many introns have been lost and added since.     

Branching in Polysphondylium whorls: Two-dimensional patterning in a three-dimensional system

During morphogenesis in the slime mold Polysphondylium pallidum, spherical masses of cells called whorls pinch off from the slug at regular intervals. Soon afterward, branches form at equidistant positions around the whorl equator. We have quantified the relationship between the number of cells in a whorl, and the number of branches that form. We find that the number of branches produced is proportional to the surface area of the whorl, suggesting that the patterning process is confined to the whorl surface. This observation is consistent with theoretical arguments that mechanisms for pattern formation would likely operate in one or two dimensions, not three.     

Identification of biochemically neutral positions in liver pyruvate kinase

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions—high sequence entropy, high surface exposure, and alanine scanning—lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a “neutrality” score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.     

Conformational Flexibility of Ubiquitin-Modified and SUMO-Modified PCNA Shown by Full-Ensemble Hybrid Methods

Ubiquitin-modified proliferating cell nuclear antigen (PCNA) and small ubiquitin-like modifier (SUMO)-modified PCNA regulate DNA damage tolerance pathways. X-ray crystal structures of these proteins suggested that they do not have much conformational flexibility because the modifiers have preferred binding sites on the surface of PCNA. By contrast, small-angle X-ray scattering analyses of these proteins suggested that they have different degrees of conformational flexibility, with SUMO-modified PCNA being more flexible. These conclusions were based on minimal-ensemble hybrid approaches, which produce unrealistic models by representing flexible proteins with only a few static structures. To overcome the limitations of minimal-ensemble hybrid approaches and to determine the degree of conformational flexibility of ubiquitin-modified PCNA and SUMO-modified PCNA, we utilized a novel full-ensemble hybrid approach. We carried out molecular simulations and small-angle X-ray scattering analyses of both proteins and obtained outstanding agreement between the full ensembles generated by the simulations and the experimental data. We found that both proteins have a high degree of conformational flexibility. The modifiers occupy many positions around the back and side of the PCNA ring. Moreover, we found no preferred ubiquitin-binding or SUMO-binding sites on PCNA. This conformational flexibility likely facilitates the recognition of downstream effector proteins and the formation of PCNA tool belts.     

The relative intranuclear positions of Barr bodies in XXX non-transformed human fibroblasts

The extent to which chromosomes exist in an ordered three-dimensional arrangement within the interphase nucleus remains unknown. As a means of tackling this issue we have developed simple statistical methods to test whether the positions of any intranuclear markers are mutually correlated. Applying these methods to non-transformed XXX human fibroblasts as a test case, we have examined the relative nuclear locations of the two Barr bodies in 2C cells of this type. Our results show that while individual Barr bodies have a highly non-random intranuclear distribution, segregating preferentially at the nuclear periphery in a plane parallel to the plane of cell growth (i.e., the substrate surface) and bisecting the nucleus, no evidence was found of any significant correlation between the positions of the two Barr bodies within a given nucleus.     

Modelling of solvent positions around polar groups in proteins

Previous analysis of the distribution of experimental solvent molecule positions around amino acid side chains showed that distinct clustering occurred close to polar or charged atoms in proteins. We have used those data to predict likely solvent positions around proteins not used in our initial analysis. We envisage that this algorithm, AQUARIUS, will be useful for finding solvent positions in electron density maps generated by protein crystallography and as useful starting positions for solvent molecules in computer simulation studies of macromolecules.     

Mechanism of product chain length determination for heptaprenyl diphosphate synthase from Bacillus stearothermophilus.

A member of the medium-chain prenyl diphosphate synthases, Bacillus stearothermophilus heptaprenyl diphosphate synthase, catalyzes the consecutive condensation of isopentenyl diphosphate with allylic diphosphate to produce (all-E)-C35 prenyl diphosphate as the ultimate product. We previously showed that the product specificity of short-chain prenyl diphosphate synthases is regulated by the structure around the first aspartate-rich motif (FARM). The FARM is also conserved in a subunit of heptaprenyl diphosphate synthase, component II', which suggests that the structure around the FARM of component II' regulates the elongation. To determine whether component II' regulates the product chain length by a mode similar to that of the short-chain prenyl diphosphate synthases, we replaced a bulky amino acid at the eighth position before the FARM of component II', isoleucine 76, by glycine and analyzed the product specificity. The mutated enzyme, I76G, can catalyze condensations of isopentenyl diphosphate beyond the native chain length of C35. Moreover, two mutated enzymes of A79Y and S80F, which have a single replacement to the aromatic residue at the fourth or the fifth position before the FARM, mainly yielded a C20 product. These results strongly suggest that a common mechanism controls the product chain length of both short-chain and medium-chain prenyl diphosphate synthases and that, in wild-type heptaprenyl diphosphate synthase, the prenyl chain can grow on the surface of the small residues at positions 79 and 80, and the elongation is precisely blocked at the length of C35 by isoleucine 76.     

Industrial Ecology as a Cultural Phenomenon: On Objectivity as a Normative Position

Taking a recent column by Braden Allenby in this journal as a starting point, we argue the need for researchers in the field of industrial ecology to reflect upon its normative aspects. We argue that the field will advance through an explicit discourse on such issues as epistemological positions, the inherent normative nature of using metaphors, and the way in which the field of study relates to the field of practice.
In "Culture and Industrial Ecology," Allenby raises, and begins to address, some critical questions that define the field of industrial ecology. We suggest that his arguments dismiss too lightly issues that are central to this developing field. The purpose of this reply is to open up for further discussion what industrial ecology is, and what it is not.
We begin by summarizing Allenby's line of reasoning. We explore his arguments around the notion of objectivity and the questions he raises about the relationship among objective science, normative positions, opinions, and fads. While Allenby makes a strong plea for industrial ecology to maintain objectivity and avoid normative positions, we contend that any positions on what should or should not be included in a field of study carry normative intent. Although science might seek to be objective in its method, scientists are rarely free of such normative positions. These positions need to be explicit and open to debate. Also, Allenby's position implies that industrial ecology as a field of study should be separated from industrial ecology as a field of practice. We argue that the interrelation of these fields provides important advantages.
Our main concern is that Allenby's position forestalls discussion on important issues at the heart of the field of study. We examine some key issues around the concept of objectivity and propose an agenda for research on the cultural and ideological aspects of industrial ecology.     

Kinetic stabilization of Bacillus licheniformis alpha-amylase through introduction of hydrophobic residues at the surface

It is generally assumed that in proteins hydrophobic residues are not favorable at solvent-exposed sites, and that amino acid substitutions on the surface have little effect on protein thermostability. Contrary to these assumptions, we have identified hyperthermostable variants of Bacillus licheniformis alpha-amylase (BLA) that result from the incorporation of hydrophobic residues at the surface. Under highly destabilizing conditions, a variant combining five stabilizing mutations unfolds 32 times more slowly and at a temperature 13 degrees C higher than the wild-type. Crystal structure analysis at 1.7 A resolution suggests that stabilization is achieved through (a) extension of the concept of increased hydrophobic packing, usually applied to cavities, to surface indentations, (b) introduction of favorable aromatic-aromatic interactions on the surface, (c) specific stabilization of intrinsic metal binding sites, and (d) stabilization of a beta-sheet by introducing a residue with high beta-sheet forming propensity. All mutated residues are involved in forming complex, cooperative interaction networks that extend from the interior of the protein to its surface and which may therefore constitute "weak points" where BLA unfolding is initiated. This might explain the unexpectedly large effect induced by some of the substitutions on the kinetic stability of BLA. Our study shows that substantial protein stabilization can be achieved by stabilizing surface positions that participate in underlying cooperatively formed substructures. At such positions, even the apparently thermodynamically unfavorable introduction of hydrophobic residues should be explored.     

Evolutionary trace analysis of scorpion toxins specific for K-channels

Scorpion alpha-K(+) channel toxins are a large family of polypeptides with a similar structure but diverse pharmacological activities. Despite many structural and functional data available at present, little progress has been made in understanding the toxin's molecular basis responsible for the functional diversification. In this paper, we report the first complete cDNA sequences of toxins belonging to subfamily 6 and identify five new members, called alpha-KTx 6.6-6.10. By analyzing the rates of mutations that occurred in the corresponding cDNAs, we suggest that accelerated evolution in toxin-coding regions may be associated with the functional diversification of this subfamily. To pinpoint sites probably involved in the functional diversity of alpha-KTx family, we analyzed this family of sequences using the evolutionary trace method. This analysis highlighted one channel-binding surface common for all the members. This surface is composed of one conserved lysine residue at position 29 assisted by other residues at positions 10, 26, 27, 32, 34, and 36. Of them, the positions 29, 32, and 34 have been reported to be the most major determinants of channel specificity. Interestingly, another contrary surface was also observed at a higher evolutionary time cut-off value, which may be involved in the binding of ERG (ether-a-go-go-related gene) channel-specific toxins. The good match between the trace residues and the functional epitopes of the toxins suggested that the evolutionary trace results reported here can be applied to predict channel-binding sites of the toxins. Because, the side-chain variation in the trace positions is strongly linked with the functional alteration and channel-binding surface transfer of alpha-KTx family, we conclude that our findings should also be important for the rational design of new toxins targeting a given potassium channel with high selectivity.     

Approximation and visualization of large-scale motion of protein surfaces

We present a method for the approximation and real-time visualization of large-scale motion of protein surfaces. A molecular surface is represented by an expansion of spherical harmonic functions, and the motion of protein atoms around their equilibrium positions is computed by normal mode analysis. The motion of the surface is approximated by projecting the normal mode vectors of the solvent-accessible atoms to the spherical harmonic representation of the molecular surface. These surface motion vectors are represented by a separated spherical harmonic expansion. Representing the surface geometry and the surface motion vectors by spherical harmonic expansions allows variable-resolution analysis and real-time display of the large-scale surface motion. This technique has been applied to interactive visualization, interactive surface manipulation, and animation.     

Side-chain contributions to membrane protein structure and stability

The molecular forces that stabilize membrane protein structure are poorly understood. To investigate these forces we introduced alanine substitutions at 24 positions in the B helix of bacteriorhodopsin and examined their effects on structure and stability. Although most of the results can be rationalized in terms of the folded structure, there are a number of surprises. (1) We find a remarkably high frequency of stabilizing mutations (17%), indicating that membrane proteins are not highly optimized for stability. (2) Helix B is kinked, with the kink centered around Pro50. The P50A mutation has no effect on stability, however, and a crystal structure reveals that the helix remains bent, indicating that tertiary contacts dominate in the distortion of this helix. (3) We find that the protein is stabilized by about 1kcal/mol for every 38A(2) of surface area buried, which is quite similar to soluble proteins in spite of their dramatically different environments. (4) We find little energetic difference, on average, in the burial of apolar surface or polar surface area, implying that van der Waals packing is the dominant force that drives membrane protein folding.     

Measuring Actin Flow in 3D Cell Protrusions

Actin dynamics is important in determining cell shape, tension, and migration. Methods such as fluorescent speckle microscopy and spatial temporal image correlation spectroscopy have been used to capture high-resolution actin turnover dynamics within cells in two dimensions. However, these methods are not directly applicable in 3D due to lower resolution and poor contrast. Here, we propose to capture actin flow in 3D with high spatial-temporal resolution by combining nanoscale precise imaging by rapid beam oscillation and fluctuation spectroscopy techniques. To measure the actin flow along cell protrusions in cell expressing actin-eGFP cultured in a type I collagen matrix, the laser was orbited around the protrusion and its trajectory was modulated in a clover-shaped pattern perpendicularly to the protrusion. Orbits were also alternated at two positions closely spaced along the protrusion axis. The pair cross-correlation function was applied to the fluorescence fluctuation from these two positions to capture the flow of actin. Measurements done on nonmoving cellular protrusion tips showed no pair-correlation at two orbital positions indicating a lack of flow of F-actin bundles. However, in some protrusions, the pair-correlation approach revealed directional flow of F-actin bundles near the protrusion surface with flow rates in the range of ∼1 μm/min, comparable to results in two dimensions using fluorescent speckle microscopy. Furthermore, we found that the actin flow rate is related to the distance to the protrusion tip. We also observed collagen deformation by concomitantly detecting collagen fibers with reflectance detection during these actin motions. The implementation of the nanoscale precise imaging by rapid beam oscillation method with a cloverleaf-shaped trajectory in conjunction with the pair cross-correlation function method provides a quantitative way of capturing dynamic flows and organization of proteins during cell migration in 3D in conditions of poor contrast.