Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Comparative evolutionary rates of introns and exons in murine rodents

Analysis of DNA sequences of 132 introns and 140 exons from 42 pairs of orthologous genes of mouse and rat was used to compare patterns of evolutionary change between introns and exons. The mean of the absolute difference in length (measured in base pairs) between the two species was nearly five times as high in the case of introns as in the case of exons. The average rate of nucleotide substitution in introns was very similar to the rate of synonymous substitution in exons, and both were about three times the rate of substitution at nonsynonymous sites in exons. G+C content of introns and exons of the same gene were correlated; but mean G+C content at the third positions of exons was significantly higher than that of introns or positions 1–2 of exons from the same gene. G+C content was conserved over evolutionary time, as indicated by strong correlations between mouse and rat; but the change in G+C content was greatest at position 3 of exons, intermediate in introns, and lowest at positions 1–2 in introns. Received: 23 December 1996 / Accepted: 1 April 1997     

The Human Obesity Gene Map: The 1996 Update

An update of the human obesity gene map up to October 1996 is presented. Evidence from Mendelian disorders exhibiting obesity as a clinical feature, single-gene mutation rodent models, quantitative trait loci uncovered in crossbreeding experiments with mouse, rat, and pig models, association and case-control studies with candidate genes, and linkage studies with genes and other markers is reviewed. All chromosomal locations of the animal loci are converted into human genome locations based on syntenic relationships between the genomes. A complete listing of all these loci reveals that only 4 of the 24 human chromosomes are not yet represented, i.e., 9, 18, 21, and Y. Several chromosome arms are characterized by the presence of several putative loci. The following arms include at least three such loci: 1p, 1q, 3p, 4q, 6p, 7q, 8p, 8q, 11p, 11q, 15q, 20q, and Xq. Studies with negative association and linkage results are also reviewed.     

Global stability of an SIR epidemic model with time delays

An SIR disease transmission model is formulated under the assumption that the force of infection at the present time depends on the number of infectives at the past. It is shown that a disease free equilibrium point is globally stable if no endemic equilibrium point exists. Further the endemic point (if it exists) is globally stable with respect to the whole state space except the neighborhood of the disease free state.Research partly supported by the Ministry of Education, Science and Culture, Japan, Grant 05640256     

Optimal local propensities for model proteins

Lattice models of proteins were used to examine the role of local propensities in stabilizing the native state of a protein, using techniques drawn from spin-glass theory to characterize the free-energy landscapes. In the strong evolutionary limit, optimal conditions for folding are achieved when the contributions from local interactions to the stability of the native state is small. Further increasing the local interactions rapidly decreases the foldability. © 1995 Wiley-Liss, Inc.     

Brief overview of control of genetic expression by antisense oligonucleotides and in vivo applications

Over the past several years, the use of synthetic oligonucleotides and functional analogs thereof as a possibly general means of controlling genetic expression has received widespread attention. Following a brief overview of some of the basic principles and strategies for this approach, attention is focused here on summarizing some recent reports of in vitro and, in particular, in vivo investigations in various animal models using phosphorothioate analogs of 2′-deoxyoligo-nucleotides. In view of these findings, which include studies related to neurobiology, this field should find significant utility in applications of the antisense method for controlling genetic expression.     

Structural analysis of zinc substitutions in the active site of thermolysin.

Native thermolysin binds a single catalytically essential zinc ion that is tetrahedrally coordinated by three protein ligands and a water molecule. During catalysis the zinc ligation is thought to change from fourfold to fivefold. Substitution of the active-site zinc with Cd2+, Mn2+, Fe2+, and Co2+ alters the catalytic activity (Holmquist B, Vallee BL, 1974, J Biol Chem 249:4601-4607). Excess zinc inhibits the enzyme. To investigate the structural basis of these changes in activity, we have determined the structures of a series of metal-substituted thermolysins at 1.7-1.9 A resolution. The structure of the Co(2+)-substituted enzyme is shown to be very similar to that of wild type except that two solvent molecules are liganded to the metal at positions that are thought to be occupied by the two oxygens of the hydrated scissile peptide in the transition state. Thus, the enhanced activity toward some substrates of the cobalt-relative to the zinc-substituted enzyme may be due to enhanced stabilization of the transition state. The ability of Zn2+ and Co2+ to accept tetrahedral coordination in the Michaelis complex, as well as fivefold coordination in the transition state, may also contribute to their effectiveness in catalysis. The Cd(2+)- and Mn(2+)-substituted thermolysins display conformational changes that disrupt the active site to varying degrees and could explain the associated reduction of activity. The conformational changes involve not only the essential catalytic residue, Glu 143, but also concerted side-chain rotations in the adjacent residues Met 120 and Leu 144. Some of these side-chain movements are similar to adjustments that have been observed previously in association with the "hinge-bending" motion that is presumed to occur during catalysis by the zinc endoproteases. In the presence of excess zinc, a second zinc ion is observed to bind at His 231 within 3.2 A of the zinc bound to native thermolysin, explaining the inhibitory effect.     

Exhaustive enumeration of protein conformations using experimental restraints.

We present an efficient new algorithm that enumerates all possible conformations of a protein that satisfy a given set of distance restraints. Rapid growth of all possible self-avoiding conformations on the diamond lattice provides construction of alpha-carbon representations of a protein fold. We investigated the dependence of the number of conformations on pairwise distance restraints for the proteins crambin, pancreatic trypsin inhibitor, and ubiquitin. Knowledge of between one and two contacts per monomer is shown to be sufficient to restrict the number of candidate structures to approximately 1,000 conformations. Pairwise RMS deviations of atomic position comparisons between pairs of these 1,000 structures revealed that these conformations can be grouped into about 25 families of structures. These results suggest a new approach to assessing alternative protein folds given a very limited number of distance restraints. Such restraints are available from several experimental techniques such as NMR, NOESY, energy transfer fluorescence spectroscopy, and crosslinking experiments. This work focuses on exhaustive enumeration of protein structures with emphasis on the possible use of NOESY-determined distance restraints.     

Establishing isostructural metal substitution in metalloproteins using 1H NMR, circular dichroism, and Fourier transform infrared spectroscopy.

Far-UV CD, 1H-NMR, and Fourier transform infrared (FTIR) spectroscopy are three of the most commonly used methods for the determination of protein secondary structure composition. These methods are compared and evaluated as a means of establishing isostructural metal substitution in metalloproteins, using the crystallographically defined rubredoxin from Desulfovibrio gigas and its well-characterized cadmium derivative as a model system. It is concluded that analysis of the FTIR spectrum of the protein amide I resonance represents the most facile and generally applicable method of determining whether the overall structure of a metalloprotein has been altered upon metal reconstitution. This technique requires relatively little biological material (ca. 300 micrograms total protein) and, unlike either CD or 1H-NMR spectroscopy, is unaffected by the presence of different metal ions, thus allowing the direct comparison of FTIR spectra before and after metal substitution.     

Interactions between two sheets of a bilayer membrane and its internal lateral pressure

We have modelled a phospholipid bilayer as two monolayer sheets which interact with each other by a coupling which depends upon the states of the lipid hydrocarbon chains in each sheet. We make use of a model (Georgallas and Pink 1982a) and its parameters, already used to study monolayer phase changes at the LC-LE transition, in order to study the lipid main transition. Although the monolayer coexistence curve can be calculated exactly, we have made use of high-temperature series expansions to calculate the critical point of the bilayer. We also present the results of computer simulations on triangular lattices for the pressure-area isotherms. We find: (i) the interaction between the sheets of a DPPC bilayer is about 1.5–2% of the maximum interaction within the plane of each sheet; (ii) the internal lateral pressure of a DPPC bilayer is about 30.5 dyne/cm; (iii) the bilayer transition enthalpy depends sensitively upon the coupling between the sheets. Should this coupling vary from sample to sample (due, possibly, to its preparation) then very different values of transition enthalpy may be measured. (iv) We present a rough rule-of-thumb for estimating the internal lateral pressure of a bilayer from a knowledge of the corresponding monolayer pressure-area isotherms.Abbreviations LC-LE liquid condensed — liquid expanded - DPPC dipalmitoylphosphatidylcholine - Q transition enthalpy Work supported in part by the Natural Sciences and Engineering Research Council of Canada