Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

Cloning and characterization of cDNAs for matrix metalloproteinases of regenerating newt limbs.

Matrix metalloproteinases (MMPs) of regenerating urodele limbs have been suggested to play crucial roles in the process of the dedifferentiation of cells in the damaged tissues and the ensuing blastema formation because the activation of MMPs is an early and conspicuous event occurring in the amputated limb. MMP cDNAs were cloned as products of the reverse transcription-PCR from cDNA libraries of newt limbs, and their structures were characterized. Three cDNAs encoding newt MMPs (2D-1, 2D-19, and 2D-24) have been cloned from second day postamputation regenerating limbs, and a cDNA (EB-1) was cloned from early bud-stage regenerating limbs. These cDNAs included the full-length coding regions. The deduced amino acid sequences of 2D-1, 2D-19, 2D-24, and EB-1 had a homology with mammalian MMP9, MMP3/10, MMP3/10, and MMP13, respectively. The basic motif of these newt MMP genes was similar to mammalian counterparts and contained regions encoding a putative signal sequence, a propeptide, an active site with three zinc-binding histidine residues, a calcium-binding domain, a hemopexin region, and three key cysteine residues. However, some unique molecular evolutionary features were also found in the newt MMPs. cDNAs of 2D-19 and 2D-24 contained a specific insertion and deletion, respectively. The insertion of 2D-19 is threonine-rich, similar to the threonine cluster found in the collagenase-like sea urchin hatching enzyme. Northern blot analysis showed that the expression levels of the newt MMPs were dramatically increased after amputation, suggesting that they play an important role(s) in tissue remodeling of the regenerating limb.

Recruitment and replacement of hippocampal neurons in young and adult chickadees: an addition to the theory of hippocampal learning.

We used [3H]thymidine to document the birth of neurons and their recruitment into the hippocampal complex (HC) of juvenile (4.5 months old) and adult blackcapped chickadees (Parus atricapillus) living in their natural surroundings. Birds received a single dose of [3H]thymidine in August and were recaptured and killed 6 weeks later, in early October. All brains were stained with Cresyl violet, a Nissl stain. The boundaries of the HC were defined by reference to the ventricular wall, the brain surface, or differences in neuronal packing density. The HC of juveniles was as large as or larger than that of adults and packing density of HC neurons was 31% higher in juveniles than in adults. Almost all of the 3H-labeled HC neurons were found in a 350-m-wide layer of tissue adjacent to the lateral ventricle. Within this layer the fraction of 3H-labeled neurons was 50% higher in juveniles than in adults. We conclude that the HC of juvenile chickadees recruits more neurons and has more neurons than that of adults. We speculate that juveniles encounter greater environmental novelty than adults and that the greater number of HC neurons found in juveniles allows them to learn more than adults. At a more general level, we suggest that (i) long-term learning alters HC neurons irreversibly; (ii) sustained hippocampal learning requires the periodic replacement of HC neurons; (iii) memories coded by hippocampal neurons are transferred elsewhere before the neurons are replaced.

The tobacco wounding-activated mitogen-activated protein kinase is encoded by SIPK

It has been demonstrated that both salicylic acid and fungal elicitors activate a 48-kDa mitogen-activated protein kinase termed salicylic acid-induced protein kinase (SIPK) in tobacco suspension cells. Here, we show that infiltration of these agents into tobacco leaves also activates SIPK. Of particular interest, infiltration of water alone activated a kinase of the same size, possibly because of wounding and/or osmotic stresses. The kinetics of kinase activation, however, differ for these different treatments. Various mechanical stresses, including cutting and wounding by abrasion, also activated a 48-kDa kinase. By using an immune-complex kinase assay with antibodies specific for SIPK or wounding-induced protein kinase, we demonstrate that this wounding-activated 48-kDa kinase is SIPK, rather than wounding-induced protein kinase, as reported [Seo, S., Okamoto, M., Seto, H., Ishizuka, K., Sano, H. & Ohashi, Y. (1995) Science 270, 1988–1992]. Activation of SIPK after wounding was associated with tyrosine phosphorylation but not with increases in SIPK mRNA or protein levels. Thus, the same mitogen-activated protein kinase, SIPK, appears to facilitate signaling for two distinct pathways that lead to disease resistance responses and wounding responses.

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.

A model of self-thinning through local competition.

Explanations of self-thinning in plant populations have focused on plant shape and packing. A dynamic model based on the structure of local interactions successfully reproduces the pattern and can be approximated to identify key parameters and relationships. The approach generates testable new explanations for differences between species and populations, unifies self-thinning with other patterns in plant population dynamics, and indicates why organisms other than plants can follow the law.

The effects of sequence context on DNA curvature.

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.

Structure and function in rhodopsin: Packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain are coupled*

A previous study of the retinitis pigmentosa mutation L125R and two designed mutations at this site, L125A and L125F, showed that these mutations cause partial or total misfolding of the opsins expressed in COS cells from the corresponding mutant opsin genes. We now report on expression and characterization of the opsins from the following retinitis pigmentosa mutants in the transmembrane domain of rhodopsin that correspond to six of the seven helices: G51A and G51V (helix A), G89D (helix B), A164V (helix D), H211P (helix E), P267L and P267R (helix F), and T297R (helix G). All the mutations caused partial misfolding of the opsins as observed by the UV/visible absorption characteristics and by separation of the expressed opsins into fractions that bound 11-cis-retinal to form the corresponding mutant rhodopsins and those that did not bind 11-cis-retinal. Further, all the mutant rhodopsins prepared from the above mutants, except for G51A, showed strikingly abnormal bleaching behavior with abnormal metarhodopsin II photointermediates. The results show that retinitis pigmentosa mutations in every one of the transmembrane helices can cause misfolding of the opsin. Therefore, on the basis of these and previous results, we conclude that defects in the packing of the transmembrane helices resulting from these mutations are relayed to the intradiscal domain, where they cause misfolding of the opsin by inducing the formation of a disulfide bond other than the native Cys-110—Cys-187 disulfide bond. Thus, there is coupling between packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain.

A general method for determining helix packing in membrane proteins in situ: Helices I and II are close to helix VII in the lactose permease of Escherichia coli

It was previously shown that coexpression of the lactose permease of Escherichia coli in two contiguous fragments leads to functional complementation. We demonstrate here that site-directed thiol crosslinking of coexpressed permease fragments can be used to determine helix proximity in situ without the necessity of purifying the permease. After coexpression of the six N-terminal (N6) and six C-terminal (C6) transmembrane helices, each with a single Cys residue, crosslinking was carried out in native membranes and assessed by the mobility of anti-C-terminal-reactive polypeptides on immunoblots. A Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 28 or 29 (helix I), but not with a Cys residue at position 27, which is on the opposite face of helix I, thereby indicating that the face of helix I containing Pro-28 and Phe-29 is close to helix VII. Similarly, a Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 52 or 53 (helix II), but not with a Cys residue at position 54. Furthermore, low-efficiency crosslinking is observed between a Cys residue at position 52 or 53 and a Cys residue at position 361 (helix XI). The results indicate that helix VII lies in close proximity to both helices I and II and that helix II is also close to helix XI. The method should be applicable to a number of different polytopic membrane proteins.

Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

A cooperative model for receptor recognition and cell adhesion: evidence from the molecular packing in the 1.6-A crystal structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The crystal structure of the pheromone Er-1 from the unicellular eukaryotic organism Euplotes raikovi was determined at 1.6 A resolution and refined to a crystallographic R factor of 19.9%. In the tightly packed crystal, two extensive intermolecular helix-helix interactions arrange the Er-1 molecules into layers. Since the putative receptor of the pheromone is a membrane-bound protein, whose extracellular C-terminal domain is identical in amino acid sequence to the soluble pheromone, the interactions found in the crystal may mimic the pheromone-receptor interactions as they occur on a cell surface. Based on this, we propose a model for the interaction between soluble pheromone molecules and their receptors. In this model, strong pheromone-receptor binding emerges as a consequence of the cooperative utilization of several weak interactions. The model offers an explanation for the results of binding studies and may also explain the adhesion between cells that occurs during mating.