Evolution of the 28S ribosomal RNA gene in anurans: regions of variability and their phylogenetic implications

Fifteen restriction sites were mapped to the 28S ribosomal RNA gene of individuals representing 54 species of frogs, two species of salamanders, a caecilian, and a lungfish. Eight of these sites were present in all species examined, and two were found in all but one species. Alignment of these conserved restriction sites revealed, among anuran 28S rRNA genes, five regions of major length variation that correspond to four of 12 previously identified divergent domains of this gene. One of the divergent domains (DD8) consists of two regions of length variation separated by a short segment that is conserved at least throughout tetrapods. Most of the insertions, deletions, and restriction-site variations identified in the 28S gene will require sequence-level analysis for a detailed reconstruction of their history. However, an insertion in DD9 that is coextensive with frogs in the suborder Neobatrachia, a BstEII site that is limited to representatives of two leptodactylid subfamilies, and a deletion in DD10 that is found only in three ranoid genera are probably synapomorphies.      

Characterization of cross-bridge elasticity and kinetics of cross-bridge cycling during force development in single smooth muscle cells

Force development in smooth muscle, as in skeletal muscle, is believed to reflect recruitment of force-generating myosin cross-bridges. However, little is known about the events underlying cross-bridge recruitment as the muscle cell approaches peak isometric force and then enters a period of tension maintenance. In the present studies on single smooth muscle cells isolated from the toad (Bufo marinus) stomach muscularis, active muscle stiffness, calculated from the force response to small sinusoidal length changes (0.5% cell length, 250 Hz), was utilized to estimate the relative number of attached cross-bridges. By comparing stiffness during initial force development to stiffness during force redevelopment immediately after a quick release imposed at peak force, we propose that the instantaneous active stiffness of the cell reflects both a linearly elastic cross-bridge element having 1.5 times the compliance of the cross-bridge in frog skeletal muscle and a series elastic component having an exponential length-force relationship. At the onset of force development, the ratio of stiffness to force was 2.5 times greater than at peak isometric force. These data suggest that, upon activation, cross-bridges attach in at least two states (i.e., low-force-producing and high-force-producing) and redistribute to a steady state distribution at peak isometric force. The possibility that the cross-bridge cycling rate was modulated with time was also investigated by analyzing the time course of tension recovery to small, rapid step length changes (0.5% cell length in 2.5 ms) imposed during initial force development, at peak force, and after 15 s of tension maintenance. The rate of tension recovery slowed continuously throughout force development following activation and slowed further as force was maintained. Our results suggest that the kinetics of force production in smooth muscle may involve a redistribution of cross-bridge populations between two attached states and that the average cycling rate of these cross-bridges becomes slower with time during contraction.     

Turn of promotor DNA by cAMP receptor protein characterized by bead model simulation of rotational diffusion

The rotation diffusion of DNA double helices and their complexes with the cAMP receptor protein (CRP) has been simulated by bead models, in order to derive information on their structure in solution by comparison with results obtained from dichroism decay measurements. Straight DNA double helices are simulated by linear, rigid strings of overlapping beads. The radius of the beads and the length of the string are increased simultaneously by the same increments from initial outer dimensions derived from crystallographic data to final values, which are fitted to experimental rotation time constants observed for short DNA fragments (less than 100 bp). The final values reflect the solvated structure with the same 'solvation layer' added in all three dimensions. The protein is simulated by overlapping beads, which are assembled to a structure very similar to that found by x-ray crystallography. Complexes of the protein with DNA are formed with the centres of palindromic DNA sites at the centre of the two helix-turn-helix-motifs of the protein with some overlap of the two components. Simulation of the experimental data obtained for CRP complexes with specific DNA in the presence of cAMP requires strong bending of the double helices. According to our simulation the DNA is almost completely wrapped around the protein both in the complexes with a 62 bp fragment containing the standard CRP site and with a 80 bp fragment containing the second binding site of the lac operon. Simulations of the data obtained for a 203 bp fragment with both binding sites suggest that the two bound CRP proteins are in contact with each other and that the DNA is wrapped around the two protein dimers. A stereochemical model is suggested with a tetrahedral arrangement of the four protein subunits, which provides the advantage that two binding sites of the protein formed by two subunits each are located favorable for tight contacts to two binding sites on bent DNA, provided that the DNA sites are separated by an integer number of helix turns. In summary, the simulations demonstrate strong bending, which can be reflected by an arc radius in the range around 50 A. According to these data the overall bending angle of our longest DNA fragment is approximately 180 degrees, and thus the protruding ends are sufficiently close to each other such that RNA polymerase, for example, could contact both helical segments.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Arabidopsis TCH4, regulated by hormones and the environment, encodes a xyloglucan endotransglycosylase.

Adaptation of plants to environmental conditions requires that sensing of external stimuli be linked to mechanisms of morphogenesis. The Arabidopsis TCH (for touch) genes are rapidly upregulated in expression in response to environmental stimuli, but a connection between this molecular response and developmental alterations has not been established. We identified TCH4 as a xyloglucan endotransglycosylase by sequence similarity and enzyme activity. Xyloglucan endotransglycosylases most likely modify cell walls, a fundamental determinant of plant form. We determined that TCH4 expression is regulated by auxin and brassinosteroids, by environmental stimuli, and during development, by a 1-kb region. Expression was restricted to expanding tissues and organs that undergo cell wall modification. Regulation of genes encoding cell wall-modifying enzymes, such as TCH4, may underlie plant morphogenetic responses to the environment.     

Acetylcholinesterase: diffusional encounter rate constants for dumbbell models of ligand.

For some enzymes, virtually every substrate molecule that encounters the entrance to the active site proceeds to reaction, at low substrate concentrations. Such diffusion-limited enzymes display high apparent bimolecular rate constants ((kcat/KM)), which depend strongly upon solvent viscosity. Some experimental studies provide evidence that acetylcholinesterase falls into this category. Interestingly, the asymmetric charge distribution of acetylcholinesterase, apparent from the crystallographic structure, suggests that its electrostatic field accelerates the encounter of its cationic substrate, acetylcholine, with the entrance to the active site. Here we report simulations of the diffusion of substrate in the electrostatic field of acetylcholinesterase. We find that the field indeed guides the substrate to the mouth of the active site. The computed encounter rate constants depend upon the particular relative geometries of substrate and enzyme that are considered to represent successful encounters. With loose reaction criteria, the computed rates exceed those measured experimentally, but the rate constants vary appropriately with ionic strength. Although more restrictive reaction criteria lower the computed rates, they also lead to unrealistic variation of the rate constants with ionic strength. That these simulations do not agree well with experiment suggests that the simple diffusion model is incomplete. Structural fluctuations in the enzyme or events after the encounter may well contribute to rate limitation.     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Evidence for "twisted plane" undulations in golden hamster sperm tails

Motile spermatozoa from the golden hamster have been arrested by rapid freezing and then fixed with glutaraldehyde at low temperature after substitution with ethylene glycol. As far as can be judged, the flagellar waveforms thus stabilized are similar to those seen in living sperm; in contrast, fixation in glutaraldehyde, without prior freezing, induces agonal changes in flagellar conformation. The characteristics waveform after freeze substitution contains three bends. Approx. half of these flagella are entirely planar. The rest are three dimensional, with the third bend displaced in a regular way from the plane containing the proximal two bends. From the geometry of these flagella, it is concluded that the plane of action of a given bending cycle undergoes a clockwise twist (from a forward viewpoint) as the cycle is succeeded by new bending cycles. This "twisted plane" undulation is quite different from helical movement. The twisting seems to occur abruptly, between cycles, as if each bending cycle has a preferred plane of action. The mechanism underlying the twisting is uncertain. However, on the basis of the angular displacements between the preferred planes, and the findings from electron microscopy, the following idea is presented as a working hypothesis: that, if the most proximal plane of bending is topographically determined by peripheral doublet 1, then successive distal planes of action are influenced predominantly by doublets 2, 3, etc., in clockwise sequence. The merits and weaknesses of this hypothesis are discussed.     

Expression of HvHMA2 in tobacco modifies Zn–Fe–Cd homeostasis

HvHMA2 is a plasma membrane P_1B-ATPase from barley that functions in Zn/Cd root-to-shoot transport. To assess the usefulness of HvHMA2 for modifying the metal content in aerial plant parts, it was expressed in tobacco under the CaMV35S promoter. Transformation with HvHMA2 did not produce one unique pattern of Zn and Cd accumulation; instead it depended on external metal supply. Thus Zn and Cd root-to-shoot translocation was facilitated, but not at all applied Zn/Cd concentrations. Metal uptake was restricted in HvHMA2-transformed plants and the level in the shoot was not enhanced. It was shown that HvHMA2 localizes to the plasma membrane of tobacco cells, and overloads the apoplast with Zn, which could explain the overall decrease in metal uptake observed. Despite the lower levels in the shoot, HvHMA2 transformants showed increased Zn sensitivity. Moreover, introduction of HvHMA2 into tobacco interfered with Fe metabolism and Fe accumulation was modified in HvHMA2-transformants in a Zn- and Cd-concentration dependent manner. The results indicate that ectopic expression of the export protein HvHMA2 in tobacco interferes with tobacco metal Zn–Cd–Fe cross-homeostasis, inducing internal mechanisms regulating metal uptake and tolerance.     

Impact of human activities on the reproduction of Hooded Vultures Necrosyrtes monachus in Burkina Faso

During the last decades, the critically endangered Hooded Vulture Necrosyrtes monachus has strongly declined across its African range. Although direct persecution has been suggested as a major cause of this decline, little is known about the impact of humans on reproductive output in West Africa. We studied the impact of human activities on the reproductive output of Hooded Vultures in the Garango area of Burkina Faso. Twenty and 56 nesting attempts were monitored, respectively, during the breeding season in 2013/14 and 2014/15, to determine reproductive success and identify causes of nest failure. Annual breeding success varied between 0.68 and 0.71 chicks fledged per breeding pair per year and productivity was assessed at 0.57 chicks fledged per territorial pair in 2014/15. The main threats imposed by humans were poaching of eggs, chicks and collection of nest materials, leading to 20% (13 out of 64 breeding attempts) of nest failures over the two years. An additional important reason for nest failure was the pruning and (partial) cutting of nest trees. Despite this high level of human interference, we found that Hooded Vulture nest success increased with proximity to human settlements, probably because breeding vultures benefit from protection by people against persecution and disturbance.