Measures of residue density in protein structures

A hierarchy of residue density assessments and packing properties in protein structures are contrasted, including a regular density, a variety of charge densities, a hydrophobic density, a polar density, and an aromatic density. These densities are investigated by alternative distance measures and also at the interface of multiunit structures. Amino acids are divided into nine structural categories according to three secondary structure states and three solvent accessibility levels. To take account of amino acid abundance differences across protein structures, we normalize the observed density by the expected density defining a density index. Solvent accessibility levels exert the predominant influence in determinations of the regular residue density. Explicitly, the regular density values vary approximately linearly with respect to solvent accessibility levels, the linearity parameters depending on the amino acid. The charge index reveals pronounced inequalities between lysine and arginine in their interactions with acidic residues. The aromatic density calculations in all structural categories parallel the regular density calculations, indicating that the aromatic residues are distributed as a random sample of all residues. Moreover, aromatic residues are found to be over-represented in the neighborhood of all amino acids. This result might be attributed to nucleation sites and protein stability being substantially associated with aromatic residues.

Atom density in protein structures

The residue environment in protein structures is studied with respect to the density of carbon (C), oxygen (O), and nitrogen (N) atoms within a certain distance (say 5 Å) of each residue. Two types of environments are evaluated: one based on side-chain atom contacts (abbreviated S-S) and the other based on all atom (side-chain + backbone) contacts (abbreviated A-A). Different atom counts are observed about nine-residue structural categories defined by three solvent accessibility levels and three secondary structure states. Among the structural categories, the S-S atom count ratios generally vary more than the A-A atom count ratios because of the fact that the backbone (O) and (N) atoms contribute equal counts. Secondary structure affects the (C) density for the A-A contacts whereas secondary structure has little influence on the (C) density for the S-S contacts. For S-S contacts, a greater density of (O) over (N) atom neighbors stands out in the environment of most amino acid types. By contrast, for A-A contacts, independent of the solvent accessibility levels, the ratio (O)/(N) is ≈1 in helical states, consistent with the geometry of α-helical residues whose side-chains tilt oppositely to the amino to carboxy α-helical axis. The highest ratio of neighbor (O)/(N) is achieved under solvent exposed conditions. This (O) vs. (N) prevalence is advantageous at the protein surface that generally exhibits an acid excess that helps to enhance protein solubility in the cell and to avoid nonspecific interactions with phosphate groups of DNA, RNA, and other plasma constituents.

The extended environment of mononuclear metal centers in protein structures

The objectives of this and the following paper are to identify commonalities and disparities of the extended environment of mononuclear metal sites centering on Cu, Fe, Mn, and Zn. The extended environment of a metal site within a protein embodies at least three layers: the metal core, the ligand group, and the second shell, which is defined here to consist of all residues distant less than 3.5 Å from some ligand of the metal core. The ligands and second-shell residues can be characterized in terms of polarity, hydrophobicity, secondary structures, solvent accessibility, hydrogen-bonding interactions, and membership in statistically significant residue clusters of different kinds. Findings include the following: (i) Both histidine ligands of type I copper ions exclusively attach the Nδ1 nitrogen of the histidine imidazole ring to the metal, whereas histidine ligands for all mononuclear iron ions and nearly all type II copper ions are ligated via the Nɛ2 nitrogen. By contrast, multinuclear copper centers are coordinated predominantly by histidine Nɛ2, whereas diiron histidine contacts are predominantly Nδ1. Explanations in terms of steric differences between Nδ1 and Nɛ2 are considered. (ii) Except for blue copper (type I), the second-shell composition favors polar residues. (iii) For blue copper, the second shell generally contains multiple methionine residues, which are elements of a statistically significant histidine–cysteine–methionine cluster. Almost half of the second shell of blue copper consists of solvent-accessible residues, putatively facilitating electron transfer. (iv) Mononuclear copper atoms are never found with acidic carboxylate ligands, whereas single Mn2+ ion ligands are predominantly acidic and the second shell tends to be mostly buried. (v) The extended environment of mononuclear Fe sites often is associated with histidine–tyrosine or histidine–acidic clusters.

Classification of mononuclear zinc metal sites in protein structures

Our study of the extended metal environment, particularly of the second shell, focuses in this paper on zinc sites. Key findings include: (i) The second shell of mononuclear zinc centers is generally more polar than hydrophobic and prominently features charged residues engaged in an abundance of hydrogen bonding with histidine ligands. Histidine–acidic or histidine–tyrosine clusters commonly overlap the environment of zinc ions. (ii) Histidine tautomeric metal bonding patterns in ligating zinc ions are mixed. For example, carboxypeptidase A, thermolysin, and sonic hedgehog possess the same ligand group (two histidines, one unibidentate acidic ligand, and a bound water), but their histidine tautomeric geometries markedly differ such that the carboxypeptidase A makes only Nδ1 contacts, thermolysin makes only Nɛ2 contacts, and sonic hedgehog uses one of each. Thus the presence of a similar ligand cohort does not necessarily imply the same topology or function at the active site. (iii) Two close histidine ligands HXmH, m ≤ 5, rarely both coordinate a single metal ion in the Nδ1 tautomeric conformation, presumably to avoid steric conflicts. Mononuclear zinc sites can be classified into six types depending on the ligand composition and geometry. Implications of the results are discussed in terms of divergent and convergent evolution.

Structural code for DNA recognition revealed in crystal structures of papillomavirus E2-DNA targets

Transcriptional regulation in papillomaviruses depends on sequence-specific binding of the regulatory protein E2 to several sites in the viral genome. Crystal structures of bovine papillomavirus E2 DNA targets reveal a conformational variant of B-DNA characterized by a roll-induced writhe and helical repeat of 10.5 bp per turn. A comparison between the free and the protein-bound DNA demonstrates that the intrinsic structure of the DNA regions contacted directly by the protein and the deformability of the DNA region that is not contacted by the protein are critical for sequence-specific protein/DNA recognition and hence for gene-regulatory signals in the viral system. We show that the selection of dinucleotide or longer segments with appropriate conformational characteristics, when positioned at correct intervals along the DNA helix, can constitute a structural code for DNA recognition by regulatory proteins. This structural code facilitates the formation of a complementary protein–DNA interface that can be further specified by hydrogen bonds and nonpolar interactions between the protein amino acids and the DNA bases.

From local to global in quasiconformal structures.

We exhibit a large class of metric spaces whose infinitesimal quasiconformal structure is strong enough to capture the global quasiconformal structure. A sufficient condition for this to happen is described in terms of a Poincaré-type inequality.

Trapping of branched DNA in microfabricated structures.

We have observed electrostatic trapping of tribranched DNA molecules undergoing electrophoresis in a microfabricated pseudo-two-dimensional array of posts. Trapping occurs in a unique transport regimen in which the electrophoretic mobility is extremely sensitive to polymer topology. The arrest of branched polymers is explained by considering their center-of-mass motion; in certain conformations, owing to the constraints imposed by the obstacles a molecule cannot advance without the center of mass first moving a short distance backwards. The depth of the resulting local potential well can be much greater than the thermal energy so that escape of an immobilized molecule can be extremely slow. We summarize the expected behavior of the mobility as a function of field strength and topology and point out that the microfabricated arrays are highly suitable for detecting an extremely small number of branched molecules in a very large population of linear molecules.

How are close residues of protein structures distributed in primary sequence?

Structurally neighboring residues are categorized according to their separation in the primary sequence as proximal (1-4 positions apart) and otherwise distal, which in turn is divided into near (5-20 positions), far (21-50 positions), very far ( > 50 positions), and interchain (from different chains of the same structure). These categories describe the linear distance histogram (LDH) for three-dimensional neighboring residue types. Among the main results are the following: (i) nearest-neighbor hydrophobic residues tend to be increasingly distally separated in the linear sequence, thus most often connecting distinct secondary structure units. (ii) The LDHs of oppositely charged nearest-neighbors emphasize proximal positions with a subsidiary maximum for very far positions. (iii) Cysteine-cysteine structural interactions rarely involve proximal positions. (iv) The greatest numbers of interchain specific nearest-neighbors in protein structures are composed of oppositely charged residues. (v) The largest fraction of side-chain neighboring residues from beta-strands involves near positions, emphasizing associations between consecutive strands. (vi) Exposed residue pairs are predominantly located in proximal linear positions, while buried residue pairs principally correspond to far or very far distal positions. The results are principally invariant to protein sizes, amino acid usages, linear distance normalizations, and over- and underrepresentations among nearest-neighbor types. Interpretations and hypotheses concerning the LDHs, particularly those of hydrophobic and charged pairings, are discussed with respect to protein stability and functionality. The pronounced occurrence of oppositely charged interchain contacts is consistent with many observations on protein complexes where multichain stabilization is facilitated by electrostatic interactions.

Crystal structure of Taq DNA polymerase in complex with an inhibitory Fab: The Fab is directed against an intermediate in the helix-coil dynamics of the enzyme

We report the crystal structure of Thermus aquaticus DNA polymerase I in complex with an inhibitory Fab, TP7, directed against the native enzyme. Some of the residues present in a helical conformation in the native enzyme have adopted a γ turn conformation in the complex. Taken together, structural information that describes alteration of helical structure and solution studies that demonstrate the ability of TP7 to inhibit 100% of the polymerase activity of the enzyme suggest that the change in conformation is probably caused by trapping of an intermediate in the helix-coil dynamics of this helix by the Fab. Antibodies directed against modified helices in proteins have long been anticipated. The present structure provides direct crystallographic evidence. The Fab binds within the DNA binding cleft of the polymerase domain, interacting with several residues that are used by the enzyme in binding the primer:template complex. This result unequivocally corroborates inferences drawn from binding experiments and modeling calculations that the inhibitory activity of this Fab is directly attributable to its interference with DNA binding by the polymerase domain of the enzyme. The combination of interactions made by the Fab residues in both the polymerase and the vestigial editing nuclease domain of the enzyme reveal the structural basis of its preference for binding to DNA polymerases of the Thermus species. The orientation of the structure-specific nuclease domain with respect to the polymerase domain is significantly different from that seen in other structures of this polymerase. This reorientation does not appear to be antibody-induced and implies remarkably high relative mobility between these two domains.

A possible role for kinase-associated protein phosphatase in the Arabidopsis CLAVATA1 signaling pathway

Continuous growth and development in plants are accomplished by meristems, groups of undifferentiated cells that persist as stem cells and initiate organs. While the structures of the apical and floral meristems in dicotyledonous plants have been well described, little is known about the underlying molecular mechanisms controlling cell proliferation and differentiation in these structures. We have shown previously that the CLAVATA1 (CLV1) gene in Arabidopsis encodes a receptor kinase-like protein that controls the size of the apical and floral meristems. Here, we show that KAPP, a gene encoding a kinase-associated protein phosphatase, is expressed in apical and young floral meristems, along with CLV1. Overexpression of KAPP mimics the clv1 mutant phenotype. Furthermore, CLV1 has kinase activity: it phosphorylates both itself and KAPP. Finally, KAPP binds and dephosphorylates CLV1. We present a model where KAPP functions as a negative regulator of the CLAVATA1 signal transduction pathway.

Gestational exposure to ethanol suppresses msx2 expression in developing mouse embryos

Ethanol acts as a teratogen in developing fetuses causing abnormalities of the brain, heart, craniofacial bones, and limb skeletal elements. To assess whether some teratogenic actions of ethanol might occur via dysregulation of msx2 expression, we examined msx2 expression in developing mouse embryos exposed to ethanol on embryonic day (E) 8 of gestation and subjected to whole mount in situ hybridization on E11–11.5 using a riboprobe for mouse msx2. Control mice exhibited expression of msx2 in developing brain, the developing limb buds and apical ectodermal ridge, the lateral and nasal processes, olfactory pit, palatal shelf of the maxilla, the eye, the lens of the eye, otic vesicle, prevertebral bodies (notochord), and endocardial cushion. Embryos exposed to ethanol in utero were significantly smaller than their normal counterparts and did not exhibit expression of msx2 in any structures. Similarly, msx2 expression, as determined by reverse transcription–PCR and Northern blot hybridization, was reduced ≈40–50% in fetal mouse calvarial osteoblastic cells exposed to 1% ethanol for 48 hr while alkaline phosphatase was increased by 2-fold and bone morphogenetic protein showed essentially no change. Transcriptional activity of the msx2 promoter was specifically suppressed by alcohol in MC3T3-E1 osteoblasts. Taken together, these data demonstrate that fetal alcohol exposure decreases msx2 expression, a known regulator of osteoblast and myoblast differentiation, and suggest that one of the “putative” mechanisms for fetal alcohol syndrome is the inhibition of msx2 expression during key developmental periods leading to developmental retardation, altered craniofacial morphogenesis, and cardiac defects.

The Kinetics of Mannose 6-Phosphate Receptor Trafficking in the Endocytic Pathway in HEp-2 Cells: The Receptor Enters and Rapidly Leaves Multivesicular Endosomes without Accumulating in a Prelysosomal Compartment

We have previously shown that in HEp-2 cells, multivesicular bodies (MVBs) processing internalized epidermal growth factor–epidermal growth factor receptor complexes mature and fuse directly with lysosomes in which the complexes are degraded. The MVBs do not fuse with a prelysosomal compartment enriched in mannose 6-phosphate receptor (M6PR) as has been described in other cell types. Here we show that the cation-independent M6PR does not become enriched in the endocytic pathway en route to the lysosome, but if a pulse of M6PR or an M6PR ligand, cathepsin D, is followed, a significant fraction of these proteins are routed from the trans-Golgi to MVBs. Accumulation of M6PR does not occur because when the ligand dissociates, the receptor rapidly leaves the MVB. At steady state, most M6PR are distributed within the trans-Golgi and trans-Golgi network and in vacuolar structures distributed in the peripheral cytoplasm. We suggest that these M6PR-rich vacuoles are on the return route from MVBs to the trans-Golgi network and that a separate stable M6PR-rich compartment equivalent to the late endosome/prelysosome stage does not exist on the endosome–lysosome pathway in these cells.

Cytoplasmically sequestered wild-type p53 protein in neuroblastoma is relocated to the nucleus by a C-terminal peptide

Cytoplasmic sequestration of wild-type p53 protein occurs in a subset of primary human tumors including breast cancer, colon cancer, and neuroblastoma (NB). The sequestered p53 localizes to punctate cytoplasmic structures that represent large protein aggregates. One functional consequence of this blocked nuclear access is impairment of the p53-mediated G1 checkpoint after DNA damage. Here we show that cytoplasmic p53 from NB cells is incompetent for specific DNA binding, probably due to its sequestration. Importantly, the C-terminal domain of sequestered p53 is masked, as indicated by the failure of a C-terminally directed antibody to detect p53 in these structures. To determine (i) which domain of p53 is involved in the aggregation and (ii) whether this phenotype is potentially reversible, we generated stable NB sublines that coexpress the soluble C-terminal mouse p53 peptide DD1 (amino acids 302–390). A dramatic phenotypic reversion occurred in five of five lines. The presence of DD1 blocked the sequestration of wild-type p53 and relocated it to the nucleus, where it accumulated. The nuclear translocation is due to shuttling of wild-type p53 by heteroligomerization to DD1, as shown by coimmunoprecipitation. As expected, the nuclear heterocomplexes were functionally inactive, since DD1 is a dominant negative inhibitor of wild-type p53. In summary, we show that nuclear access of p53 can be restored in NB cells.

Statistical methods for characterizing similarities and differences between semantic structures

This paper describes a variety of statistical methods for obtaining precise quantitative estimates of the similarities and differences in the structures of semantic domains in different languages. The methods include comparing mean correlations within and between groups, principal components analysis of interspeaker correlations, and analysis of variance of speaker by question data. Methods for graphical displays of the results are also presented. The methods give convergent results that are mutually supportive and equivalent under suitable interpretation. The methods are illustrated on the semantic domain of emotion terms in a comparison of the semantic structures of native English and native Japanese speaking subjects. We suggest that, in comparative studies concerning the extent to which semantic structures are universally shared or culture-specific, both similarities and differences should be measured and compared rather than placing total emphasis on one or the other polar position.