Dynamical structure of transfer RNA studied by normal mode analysis

The internal motion of yeast phenylalanine transfer RNA is studied by normal mode analysis in extended dihedral angle space, in which the flexibility of five-membered ribose rings is treated faithfully by introducing a variable for its pseudo-rotational motion. Analysis of global molecular motions reveals that the molecule is very soft. We show that this softness comes not from the property of the “material” comprising the molecule but from its slender shape. Analysis of thermal distance fluctuations reveals that this molecule can be regarded as consisting dynamically of three blocks. Thermal fluctuations of the mainchain dihedral angles show rigidity of the anticodon region. They also show flexibility of regions around nonstacking bases. Base-stacking interactions cause suppression of the correlated functions of mainchain dihedral angles beyond a ribose ring. We analyze the thermal fluctuation of parameters describing the positions of two adjacent bases. Fluctuations of relative translational parameters in the anticodon and acceptor stem regions are found to be larger than those in other stem regions. The relative translational motions cause the two stem regions to undergo global twisting and bending motions. We show that the role of pseudo-rotational motion of sugars is important in regions around bases which are involved in nonregular interactions. Received: 29 October 1998 / Revised version: 8 February 1999 / Accepted: 11 February 1999     

Domain structure of neurofilament subunits as revealed by monoclonal antibodies

Monoclonal antibodies have been prepared against purified neurofilament (NF) subunits (NF68, NF150, and NF200). From 25 fusions, several hundred strongly positive antibodies have been obtained. Among them are antibodies against the specific subunits as well as antibodies recognizing common antigenic determinants. These have all been characterized according to the following properties: ELISA (enzyme-linked immunosorbant assay) testing against each subunit, immunoblots against enriched neurofilament preparation, immunoblots of cyanogen bromide or chymotrypsin-treated neurofilaments, immunofluorescence with PC12 cells, and immunohistochemistry of cerebellum. Whereas the antibodies against the NF68 and NF150 appear to react with single cyanogen bromide fragments, the antibodies against the NF200 react with multiple cyanogen bromide fragments. These data are consistent with the hypothesis that the NF200 is partially composed of several repeated structural determinants. Furthermore, all of the antibodies that react with the NF200 recognize the solubilized "sidearm" domain from limited chymotryptic digestions. The locations of the common and variable domains of the three subunits are discussed in light of these results.     

The relationship between periodic dinucleotides and the nucleosomal DNA deformation revealed by normal mode analysis

Nucleosomes, which contain DNA and proteins, are the basic unit of eukaryotic chromatins. Polymers such as DNA and proteins are dynamic, and their conformational changes can lead to functional changes. Periodic dinucleotide patterns exist in nucleosomal DNA chains and play an important role in the nucleosome structure. In this paper, we use normal mode analysis to detect significant structural deformations of nucleosomal DNA and investigate the relationship between periodic dinucleotides and DNA motions. We have found that periodic dinucleotides are usually located at the peaks or valleys of DNA and protein motions, revealing that they dominate the nucleosome dynamics. Also, a specific dinucleotide pattern CA/TG appears most frequently.     

Dynamic structure of subtilisin-eglin c complex studied by normal mode analysis

Normal mode analysis of subtilisin-eglin c complex was performed to investigate the dynamics at the interface between the enzyme and the inhibitor. The internal motions of the complex calculated from the normal modes were divided into three parts: the internal motions changing the shape of each molecule, the external rigid-body motions changing their mutual dispositions, and the coupling between the internal and external motions. From the results of the analysis, the following characteristic features were found in the dynamics at the interface regions: 1) negative correlation between the internal and external motions within each molecule, and 2) positive correlation between the external motions of the two molecules. The former decreases the apparent amplitudes of motions at the interface. The latter minimizes the interference between individual motions of the two molecules. These dynamic characteristics allow the enzyme and the inhibitor to move as freely as possible. This finding suggests that the experimental evidence of the large entropy gain on binding should be attributed not only to strong hydrophobic interactions, but also to the dynamic structure of the complex, which is found to minimize an unavoidable loss of the conformational entropy on binding. Proteins 32:324–333, 1998. © 1998 Wiley-Liss, Inc.     

A two-domain structure for the two subunits of NAD(P)H:quinone acceptor oxidoreductase.

NAD(P)H:quinone acceptor oxidoreductase (EC 1.6.99.2) (DT-diaphorase) is a FAD-containing reductase that catalyzes a unique 2-electron reduction of quinones. It consists of 2 identical subunits. In this study, it was found that the carboxyl-terminal portion of the 2 subunits can be cleaved by various proteases, whereas the amino-terminal portion cannot. It was also found that proteolytic digestion of the enzyme can be blocked by the prosthetic group FAD, substrates NAD(P)H and menadione, and inhibitors dicoumarol and phenindione. Interestingly, chrysin and Cibacron blue, 2 additional inhibitors, cannot protect the enzyme from proteolytic digestion. The results obtained from this study indicate that the subunit of the quinone reductase has a 2-domain structure, i.e., an amino-terminal compact domain and a carboxyl-terminal flexible domain. A structural model of the quinone reductase is generated based on results obtained from amino-terminal and carboxyl-terminal protein sequence analyses and electrospray mass spectral analyses of hydrolytic products of the enzyme generated by trypsin, chymotrypsin, and Staphylococcus aureus protease. Furthermore, based on the data, it is suggested that the binding of substrates involves an interaction between 2 structural domains.     

A novel two-domain architecture within the amino acid kinase enzyme family revealed by the crystal structure of Escherichia coli glutamate 5-kinase

Glutamate 5-kinase (G5K) makes the highly unstable product glutamyl 5-phosphate (G5P) in the initial, controlling step of proline/ornithine synthesis, being feedback-inhibited by proline or ornithine, and causing, when defective, clinical hyperammonaemia. We determined two crystal structures of G5K from Escherichia coli, at 2.9 A and 2.5 A resolution, complexed with glutamate and sulphate, or with G5P, sulphate and the proline analogue 5-oxoproline. E. coli G5K presents a novel tetrameric (dimer of dimers) architecture. Each subunit contains a 257 residue AAK domain, typical of acylphosphate-forming enzymes, with characteristic alpha(3)beta(8)alpha(4) sandwich topology. This domain is responsible for catalysis and proline inhibition, and has a crater on the beta sheet C-edge that hosts the active centre and bound 5-oxoproline. Each subunit contains a 93 residue C-terminal PUA domain, typical of RNA-modifying enzymes, which presents the characteristic beta(5)beta(4) sandwich fold and three alpha helices. The AAK and PUA domains of one subunit associate non-canonically in the dimer with the same domains of the other subunit, leaving a negatively charged hole between them that hosts two Mg ions in one crystal, in line with the G5K requirement for free Mg. The tetramer, formed by two dimers interacting exclusively through their AAK domains, is flat and elongated, and has in each face, pericentrically, two exposed active centres in alternate subunits. This would permit the close apposition of two active centres of bacterial glutamate-5-phosphate reductase (the next enzyme in the proline/ornithine-synthesising route), supporting the postulated channelling of G5P. The structures clarify substrate binding and catalysis, justify the high glutamate specificity, explain the effects of known point mutations, and support the binding of proline near glutamate. Proline binding may trigger the movement of a loop that encircles glutamate, and which participates in a hydrogen bond network connecting active centres, which is possibly involved in the cooperativity for glutamate.     

Understanding ribosome assembly: the structure of in vivo assembled immature 30S subunits revealed by cryo-electron microscopy

Four decades after early in vitro assembly studies demonstrated that ribosome assembly is a controlled process, our understanding of ribosome assembly is still incomplete. Just as structure determination has been so important to understanding ribosome function, so too will it be critical to sorting out the assembly process. Here, we used a viable deletion in the yjeQ gene, a recognized ribosome assembly factor, to isolate and structurally characterize immature 30S subunits assembled in vivo. These small ribosome subunits contained unprocessed 17S rRNA and lacked some late ribosomal proteins. Cryo-electron microscopy reconstructions revealed that the presence of precursor sequences in the rRNA induces a severe distortion in the 3' minor domain of the subunit involved in the decoding of mRNA and interaction with the large ribosome subunit. These findings suggest that rRNA processing events induce key local conformational changes directing the structure toward the mature assembly. We concluded that rRNA processing, folding, and the entry of tertiary r-proteins are interdependent events in the late stages of 30S subunit assembly. In addition, we demonstrate how studies of emerging assembly factors in ribosome biogenesis can help to elucidate the path of subunit assembly in vivo.     

Community structure in a methanotroph biofilter as revealed by phospholipid fatty acid analysis

The microbial community structure of two biofilters used for the oxidation of methane and organic trace gases generated in landfills was analysed by phospholipid fatty acid composition. Community structure varied with biofilter depth, reflecting varying conditions of substrate supply as well as of organic carbon content, nutrient status and osmotic stress determined by the different materials used for the individual biofilter layers. Both biofilters were dominated by type II methanotrophs. In the biofilter charged with landfill gas containing significant amounts of trace organics, fatty acid 18:1omega7c constituted 87% of the methanotrophic PLFA, while the recognised signature fatty acids 16:1omega8 and 18:1omega8, which were well represented in the other biofilter, were entirely absent. This indicates the development of a highly specific methanotrophic population, presumably as a result of the adaption to continuous organic trace gas exposure.     

A novel mode of Gleevec binding is revealed by the structure of spleen tyrosine kinase

Spleen tyrosine kinase (Syk) is a non-receptor tyrosine kinase required for signaling from immunoreceptors in various hematopoietic cells. Phosphorylation of two tyrosine residues in the activation loop of the Syk kinase catalytic domain is necessary for signaling, a phenomenon typical of tyrosine kinase family members. Syk in vitro enzyme activity, however, does not depend on phosphorylation (activation loop tyrosine --> phenylalanine mutants retain catalytic activity). We have determined the x-ray structure of the unphosphorylated form of the kinase catalytic domain of Syk. The enzyme adopts a conformation of the activation loop typically seen only in activated, phosphorylated tyrosine kinases, explaining why Syk does not require phosphorylation for activation. We also demonstrate that Gleevec (STI-571, Imatinib) inhibits the isolated kinase domains of both unphosphorylated Syk and phosphorylated Abl with comparable potency. Gleevec binds Syk in a novel, compact cis-conformation that differs dramatically from the binding mode observed with unphosphorylated Abl, the more Gleevec-sensitive form of Abl. This finding suggests the existence of two distinct Gleevec binding modes: an extended, trans-conformation characteristic of tight binding to the inactive conformation of a protein kinase and a second compact, cis-conformation characteristic of weaker binding to the active conformation. Finally, the Syk-bound cis-conformation of Gleevec bears a striking resemblance to the rigid structure of the nonspecific, natural product kinase inhibitor staurosporine.     

Population structure in African Drosophila melanogaster revealed by microsatellite analysis

Tropical sub-Saharan regions are considered to be the geographical origin of Drosophila melanogaster. Starting from there, the species colonized the rest of the world after the last glaciation about 10 000 years ago. Consistent with this demographic scenario, African populations have been shown to harbour higher levels of microsatellite and sequence variation than cosmopolitan populations. Nevertheless, limited information is available on the genetic structure of African populations. We used X chromosomal microsatellite variation to study the population structure of D. melanogaster populations using 13 sampling sites in North, West and East Africa. These populations were compared to six European and one North American population. Significant population structure was found among African D. melanogaster populations. Using a Bayesian method for inferring population structure we detected two distinct groups of populations among African D. melanogaster. Interestingly, the comparison to cosmopolitan D. melanogaster populations indicated that one of the divergent African groups is closely related to cosmopolitan flies. Low, but significant levels of differentiation were observed for sub-Saharan D. melanogaster populations from West and East Africa.     

Genetic structure of Mediterranean populations revealed by Y-chromosome haplotype analysis

The allelic variability at six Y-chromosome-specific polymorphisms (YAP, DYS19, DYS389-I, DYS390, DYS391, and DYS392) was used to generate male-specific haplotypes in 333 males representing 12 population samples from the region around the Mediterranean sea. Extreme interindividual variation was observed, as more than 160 distinct Y-chromosome variants could be defined as six-locus haplotypes. Concomitant with this high variability, low levels of population genetic structure were observed. In particular, a "core" of populations directly facing the north and the east of the Mediterranean basin, from the Middle East to the Italian Peninsula, was found to be genetically undifferentiated. This observation, supported by a reanalysis of Y-specific binary polymorphisms in the same populations, suggests that at least part of the male-specific gene pools of these populations has either a very recent common origin (that could be related with the Neolithic demic diffusion hypothesis), and/or that gene flow has played a significant role in shaping the patterns of genetic variability in this region. In agreement with both hypotheses, we found that the spatial distribution of DYS392 alleles revealed a marked differentiation between the East and the West of the Mediterranean area. Through the analysis of microsatellite variation, the time to the most recent common ancestor (TMRCA) of the YAP(+) sublineage 4 has been estimated. The estimations, based on two different data sets, turn out to be quite recent (7,000-11,000 YBP), suggesting that this lineage may have been first introduced into Southern Europe through Neolithic migrations from the Middle East.     

A novel mode of DNA recognition by a beta-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA.

The 3D solution structure of the GCC-box binding domain of a protein from Arabidopsis thaliana in complex with its target DNA fragment has been determined by heteronuclear multidimensional NMR in combination with simulated annealing and restrained molecular dynamic calculation. The domain consists of a three-stranded anti-parallel beta-sheet and an alpha-helix packed approximately parallel to the beta-sheet. Arginine and tryptophan residues in the beta-sheet are identified to contact eight of the nine consecutive base pairs in the major groove, and at the same time bind to the sugar phosphate backbones. The target DNA bends slightly at the central CG step, thereby allowing the DNA to follow the curvature of the beta-sheet.     

The roles of the extrinsic subunits in Photosystem II as revealed by EPR

The effects of selective removal of extrinsic proteins on donor side electron transport in oxygen-evolving PS II particles were examined by monitoring the decay time of the EPR signal from the oxidized secondary donor, Z⁺, and the amplitude of the multiline manganese EPR signal. Removal of the 16 and 24 kDa proteins by washing with 1 M NaCl inhibits oxygen evolution, but rapid electron transfer to Z⁺ still occurs as evidenced by the near absence of Signal II_f. The absence of a multiline EPR signal shows that NaCl washing induces a modification of the oxygen-evolving complex which prevents the formation of the S₂ state. This modification is different from the one induced by chloride depletion of PS II particles, since in these a large multiline EPR signal is found. After removal of the 33 kDa protein with 1 M MgCl₂, Signal II_f is generated after a light flash. Readdition of the 33 kDa component to the depleted membranes accelerates the reduction of Z⁺. Added calcium ions show a similar effect. These findings suggest that partial advancement through the oxygen-evolving cycle can occur in the absence of the 16 and 24 kDa proteins. The 33 kDa protein, on the other hand, may be necessary for such reactions to take place.     

Coarse-grained normal mode analysis in structural biology

The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.     

Compartmentalized structure of the plasma membrane for receptor movements as revealed by a nanometer-level motion analysis

Movements of transferrin and alpha 2-macroglobulin receptor molecules in the plasma membrane of cultured normal rat kidney (NRK) fibroblastic cells were investigated by video-enhanced contrast optical microscopy with 1.8 nm spatial precision and 33 ms temporal resolution by labeling the receptors with the ligand-coated nanometer-sized colloidal gold particles. For both receptor species, most of the movement trajectories are of the confined diffusion type, within domains of approximately 0.25 microns2 (500-700 nm in diagonal length). Movement within the domains is random with a diffusion coefficient approximately 10(-9) cm2/s, which is consistent with that expected for free Brownian diffusion of proteins in the plasma membrane. The receptor molecules move from one domain to one of the adjacent domains at an average frequency of 0.034 s-1 (the residence time within a domain approximately 29 s), indicating that the plasma membrane is compartmentalized for diffusion of membrane receptors and that long- range diffusion is the result of successive intercompartmental jumps. The macroscopic diffusion coefficients for these two receptor molecules calculated on the basis of the compartment size and the intercompartmental jump rate are approximately 2.4 x 10(-11) cm2/s, which is consistent with those determined by averaging the long-term movements of many particles. Partial destruction of the cytoskeleton decreased the confined diffusion mode, increased the simple diffusion mode, and induced the directed diffusion (transport) mode. These results suggest that the boundaries between compartments are made of dynamically fluctuating membrane skeletons (membrane-skeleton fence model).     

Reduced genetic structure of north Ethiopian cattle revealed by Y-chromosome analysis

Ethiopia is considered to be a putative migratory corridor for both Near-East Bos taurine and Arabian and Indian B. indicus cattle into East Africa. African pastoralism, which is associated with adaptation to specific habitats and farming systems, has contributed to the composite constitution of Ethiopian cattle. We analyse, for the first time, five Y-chromosome microsatellite markers from seven north Ethiopian cattle populations, using a European Holstein-Friesian population as a reference, to assess the paternal gene pool and to explore the mechanisms behind the genetic structure. Our results reveal that the indicine alleles predominate in the present populations, with only one animal in the Arado carrying the taurine alleles. The north Ethiopian cattle populations with one exception (Abergelle) are characterized by a general low Y-chromosome haplotype diversity, as well as by a reduced interpopulation variance (Phi(ST)=4.0%), which can be a result of strong male-mediated selective sweeps. Population structure revealed by multidimensional-scaling analysis differentiates two populations (Arado and Abergelle) from the rest. Analysis of molecular variance does not lend support to the traditional classification for the populations, which is mainly based on physical characteristics. A network analysis indicates two closely related founding haplotypes accounting for a large proportion (50.0% in Abergelle and 85.0-94.7% in others) of north Ethiopian cattle Y-chromosomes. Our findings point to a common, but limited, paternal origin of the north Ethiopian cattle populations and strong male-mediated gene flow among them. The findings also provide insight into the historical immigration of cattle into East Africa.     

Relatedness structure in Rhododendron metternichii var. hondoense revealed by microsatellite analysis

The relatedness structure of Rhododendron metternichii Sieb. et Zucc. var. hondoense Nakai was analysed in a 150 x 70-m quadrat in Hiroshima Prefecture, western Japan. The population of R. metternichii occurred as three subpopulations at the study site (A1-A3). Pairwise relatedness based on microsatellite genotypes at eight loci and Mantel tests revealed a hierarchical structure of relatedness within and among subpopulations: (i) relatedness between individuals within 10 m of one another was significantly positive; (ii) relatedness between individuals in the same subpopulation was significantly positive, but negative between individuals in distant subpopulations; and (iii) relatedness was not significantly different from zero among neighbouring subpopulations. In detail, however, relatedness within each subpopulation was significantly positive in subpopulation A1, relatively weak but significantly positive in subpopulation A2, and not significantly different from zero in subpopulation A3. Relatedness within each subpopulation was inversely related to correlations between interindividual distance and relatedness. The aggregation of related individuals at short distances from one another may lead to decreasing relatedness within subpopulations as a whole. Moreover, negative correlations between interindividual distance and relatedness corresponded to high flowering densities at less than 10-m distance, implying that high flowering densities reduce pollinator foraging distance and lead to stronger genetic structure within subpopulations. Small individuals (< 2.0 m in height) showed stronger genetic structure compared with that of large individuals (> or = 2.0 m in height). The different relatedness structure within and among subpopulations may be caused by various degrees of gene flow affected by distribution patterns of individuals and population density.