A maximum likelihood algorithm for genome mapping of cytogenetic loci from meiotic configuration data.

Frequencies of meiotic configurations in cytogenetic stocks are dependent on chiasma frequencies in segments defined by centromeres, breakpoints, and telomeres. The expectation maximization algorithm is proposed as a general method to perform maximum likelihood estimations of the chiasma frequencies in the intervals between such locations. The estimates can be translated via mapping functions into genetic maps of cytogenetic landmarks. One set of observational data was analyzed to exemplify application of these methods, results of which were largely concordant with other comparable data. The method was also tested by Monte Carlo simulation of frequencies of meiotic configurations from a monotelodisomic translocation heterozygote, assuming six different sample sizes. The estimate averages were always close to the values given initially to the parameters. The maximum likelihood estimation procedures can be extended readily to other kinds of cytogenetic stocks and allow the pooling of diverse cytogenetic data to collectively estimate lengths of segments, arms, and chromosomes.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

A peptide that inhibits hydroxyapatite growth is in an extended conformation on the crystal surface

Proteins play an important role in the biological mechanisms controlling hard tissue development, but the details of molecular recognition at inorganic crystal interfaces remain poorly characterized. We have applied a recently developed homonuclear dipolar recoupling solid-state NMR technique, dipolar recoupling with a windowless sequence (DRAWS), to directly probe the conformation of an acidic peptide adsorbed to hydroxyapatite (HAP) crystals. The phosphorylated hexapeptide, DpSpSEEK (N6, where pS denotes phosphorylated serine), was derived from the N terminus of the salivary protein statherin. Constant-composition kinetic characterization demonstrated that, like the native statherin, this peptide inhibits the growth of HAP seed crystals when preadsorbed to the crystal surface. The DRAWS technique was used to measure the internuclear distance between two 13C labels at the carbonyl positions of the adjacent phosphoserine residues. Dipolar dephasing measured at short mixing times yielded a mean separation distance of 3.2 ± 0.1 Å. Data obtained by using longer mixing times suggest a broad distribution of conformations about this average distance. Using a more complex model with discrete α-helical and extended conformations did not yield a better fit to the data and was not consistent with chemical shift analysis. These results suggest that the peptide is predominantly in an extended conformation rather than an α-helical state on the HAP surface. Solid-state NMR approaches can thus be used to determine directly the conformation of biologically relevant peptides on HAP surfaces. A better understanding of peptide and protein conformation on biomineral surfaces may provide design principles useful for the modification of orthopedic and dental implants with coatings and biological growth factors that are designed to enhance biocompatibility with surrounding tissue.

Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

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.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

RANK is the intrinsic hematopoietic cell surface receptor that controls osteoclastogenesis and regulation of bone mass and calcium metabolism

We have generated RANK (receptor activator of NF-κB) nullizygous mice to determine the molecular genetic interactions between osteoprotegerin, osteoprotegerin ligand, and RANK during bone resorption and remodeling processes. RANK−/− mice lack osteoclasts and have a profound defect in bone resorption and remodeling and in the development of the cartilaginous growth plates of endochondral bone. The osteopetrosis observed in these mice can be reversed by transplantation of bone marrow from rag1−/− (recombinase activating gene 1) mice, indicating that RANK−/− mice have an intrinsic defect in osteoclast function. Calciotropic hormones and proresorptive cytokines that are known to induce bone resorption in mice and human were administered to RANK−/− mice without inducing hypercalcemia, although tumor necrosis factor α treatment leads to the rare appearance of osteoclast-like cells near the site of injection. Osteoclastogenesis can be initiated in RANK−/− mice by transfer of the RANK cDNA back into hematopoietic precursors, suggesting a means to critically evaluate RANK structural features required for bone resorption. Together these data indicate that RANK is the intrinsic cell surface determinant that mediates osteoprotegerin ligand effects on bone resorption and remodeling as well as the physiological and pathological effects of calciotropic hormones and proresorptive cytokines.

Using rare mutations to estimate population divergence times: A maximum likelihood approach

In this paper we propose a method to estimate by maximum likelihood the divergence time between two populations, specifically designed for the analysis of nonrecurrent rare mutations. Given the rapidly growing amount of data, rare disease mutations affecting humans seem the most suitable candidates for this method. The estimator RD, and its conditional version RDc, were derived, assuming that the population dynamics of rare alleles can be described by using a birth–death process approximation and that each mutation arose before the split of a common ancestral population into the two diverging populations. The RD estimator seems more suitable for large sample sizes and few alleles, whose age can be approximated, whereas the RDc estimator appears preferable when this is not the case. When applied to three cystic fibrosis mutations, the estimator RD could not exclude a very recent time of divergence among three Mediterranean populations. On the other hand, the divergence time between these populations and the Danish population was estimated to be, on the average, 4,500 or 15,000 years, assuming or not a selective advantage for cystic fibrosis carriers, respectively. Confidence intervals are large, however, and can probably be reduced only by analyzing more alleles or loci.

Factors affecting likelihood of applicants being offered a place in medical schools in the United Kingdom in 1996 and 1997: retrospective study

Objective To assess the relation between a range of measures and the likelihood of applicants to medical schools in the United Kingdom being offered a place overall and at each medical school, with particular emphasis on ethnic minority applicants.

Galactose-extended glycans of antibodies produced by transgenic plants

Plant-specific N-glycosylation can represent an important limitation for the use of recombinant glycoproteins of mammalian origin produced by transgenic plants. Comparison of plant and mammalian N-glycan biosynthesis indicates that β1,4-galactosyltransferase is the most important enzyme that is missing for conversion of typical plant N-glycans into mammalian-like N-glycans. Here, the stable expression of human β1,4-galactosyltransferase in tobacco plants is described. Proteins isolated from transgenic tobacco plants expressing the mammalian enzyme bear N-glycans, of which about 15% exhibit terminal β1,4-galactose residues in addition to the specific plant N-glycan epitopes. The results indicate that the human enzyme is fully functional and localizes correctly in the Golgi apparatus. Despite the fact that through the modified glycosylation machinery numerous proteins have acquired unusual N-glycans with terminal β1,4-galactose residues, no obvious changes in the physiology of the transgenic plants are observed, and the feature is inheritable. The crossing of a tobacco plant expressing human β1,4-galactosyltransferase with a plant expressing the heavy and light chains of a mouse antibody results in the expression of a plantibody that exhibits partially galactosylated N-glycans (30%), which is approximately as abundant as when the same antibody is produced by hybridoma cells. These results are a major step in the in planta engineering of the N-glycosylation of recombinant antibodies.

Crystal structure of the ectodomain of Methuselah, a Drosophila G protein-coupled receptor associated with extended lifespan

The Drosophila mutant methuselah (mth) was identified from a screen for single gene mutations that extended average lifespan. Mth mutants have a 35% increase in average lifespan and increased resistance to several forms of stress, including heat, starvation, and oxidative damage. The protein affected by this mutation is related to G protein-coupled receptors of the secretin receptor family. Mth, like secretin receptor family members, has a large N-terminal ectodomain, which may constitute the ligand binding site. Here we report the 2.3-Å resolution crystal structure of the Mth extracellular region, revealing a folding topology in which three primarily β-structure-containing domains meet to form a shallow interdomain groove containing a solvent-exposed tryptophan that may represent a ligand binding site. The Mth structure is analyzed in relation to predicted Mth homologs and potential ligand binding features.

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

Induction of a Carbon-Starvation-Related Proteolysis in Whole Maize Plants Submitted to Light/Dark Cycles and to Extended Darkness1

Three-week-old maize (Zea mays L.) plants were submitted to light/dark cycles and to prolonged darkness to investigate the occurrence of sugar-limitation effects in different parts of the whole plant. Soluble sugars fluctuated with light/dark cycles and dropped sharply during extended darkness. Significant decreases in protein level were observed after prolonged darkness in mature roots, root tips, and young leaves. Glutamine and asparagine (Asn) changed in opposite ways, with Asn increasing in the dark. After prolonged darkness the increase in Asn accounted for most of the nitrogen released by protein breakdown. Using polyclonal antibodies against a vacuolar root protease previously described (F. James, R. Brouquisse, C. Suire, A. Pradet, P. Raymond [1996] Biochem J 320: 283–292) or the 20S proteasome, we showed that the increase in proteolytic activities was related to an enrichment of roots in the vacuolar protease, with no change in the amount of 20S proteasome in either roots or leaves. Our results show that no significant net proteolysis is induced in any part of the plant during normal light/dark cycles, although changes in metabolism and growth appear soon after the beginning of the dark period, and starvation-related proteolysis probably appears in prolonged darkness earlier in sink than in mature tissues.

Cell Adhesion Molecule L1 in Folded (Horseshoe) and Extended Conformations

We have investigated the structure of the cell adhesion molecule L1 by electron microscopy. We were particularly interested in the conformation of the four N-terminal immunoglobulin domains, because x-ray diffraction showed that these domains are bent into a horseshoe shape in the related molecules hemolin and axonin-1. Surprisingly, rotary-shadowed specimens showed the molecules to be elongated, with no indication of the horseshoe shape. However, sedimentation data suggested that these domains of L1 were folded into a compact shape in solution; therefore, this prompted us to look at the molecules by an alternative technique, negative stain. The negative stain images showed a compact shape consistent with the expected horseshoe conformation. We speculate that in rotary shadowing the contact with the mica caused a distortion of the protein, weakening the bonds forming the horseshoe and permitting the molecule to extend. We have thus confirmed that the L1 molecule is primarily in the horseshoe conformation in solution, and we have visualized for the first time its opening into an extended conformation. Our study resolves conflicting interpretations from previous electron microscopy studies of L1.

Cellular aging, destabilization, and cancer.

Three major characteristics of aging in animals are a slowdown of cell proliferation, an increase in residual bodies associated with age pigments, and a marked increase in the likelihood of neoplastic transformation. The 28 L subline of the NIH 3T3 line of mouse embryo fibroblasts exhibits all these characteristics when held at confluence for extended periods. The impairment of proliferation is the first behavioral characteristic detected in low density subcultures from the confluent cultures, and it persists through many cell generations of exponential multiplication. There is an equal degree of growth impairment among replicate cultures (lineages) recovered after each of 2 successive rounds of confluence, although heterogeneity appears after the third round. The growth impairment pervades the entire cell population of each lineage. The degree and duration of impairment increase with repeated rounds of confluence. A marked increase of residual bodies characteristic of age pigments occurs in the cytoplasm of all the cells kept under prolonged confluence. Neoplastic transformation first appears as foci of multilayered cells on a monolayered background of nontransformed cells. The transformed cells arise at different times in the lineages and originate from a very small fraction of the population. The transformed cells selectively overgrow the entire population in successive rounds of confluence leading to an increase in saturation density of each lineage at different times. Under cloning conditions, isolated colonies of transformed cells develop more slowly than colonies of nontransformed cells but eventually reach a higher population density. The regularity of persistent growth impairment among the lineages and the appearance of large numbers of residual bodies in all the cells of each population are more characteristic of an epigenetic process than of specific local mutations. although random chromosomal lesions cannot be ruled out. By contrast, the low frequency and stochastic character of neoplastic transformation are consistent with a conventional genetic origin. The advent in long-term confluent NIH 3T3 cultures of three cardinal characteristics of cellular aging in vivo recommends it as a model for aging cells.