Efficient silencing of endogenous microRNAs using artificial microRNAs in Arabidopsis thaliana

We report here that the expression of endogenous microRNAs (miRNAs) can be efficiently silenced in Arabidopsis thaliana (Arabidopsis) using artificial miRNA (amiRNA) technology. We demonstrate that an amiRNA designed to target a mature miRNA directs silencing against all miRNA family members, whereas an amiRNA designed to target the stem-loop region of a miRNA precursor transcript directs silencing against only the individual family member targeted. Furthermore, our results indicate that amiRNAs targeting both the mature miRNA and stem-loop sequence direct RNA silencing through cleavage of the miRNA precursor transcript, which presumably occurs in the nucleus of a plant cell during the initial stages of miRNA biogenesis. This suggests that small RNA (sRNA)-guided RNA cleavage in plants occurs not only in the cytoplasm, but also in the nucleus. Many plant miRNA gene families have been identified via sequencing and bioinformatic analysis, but, to date, only a small tranche of these have been functionally characterized due to a lack of effective forward or reverse genetic tools. Our findings therefore provide a new and powerful reverse-genetic tool for the analysis of miRNA function in plants. © The Author 2010. Published by the Molecular Plant Shanghai Editorial Office in association with Oxford University Press on behalf of CSPP and IPPE, SIBS, CAS.

Alteration of the C-terminal ligand specificity of the Erbin PDZ domain by allosteric mutational effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.

Arsenic trioxide promotes mitochondrial DNA mutation and cell apoptosis in primary APL cells and NB4 cell line

This study aimed to investigate the effects of arsenic trioxide (As2O3) on the mitochondrial DNA (mtDNA) of acute promyelocytic leukemia (APL) cells. The NB4 cell line was treated with 2.0 μmol/L As2O3in vitro, and the primary APL cells were treated with 2.0 μmol/L As2O3in vitro and 0.16 mg kg-1 d-1 As2O3in vivo. The mitochondrial DNA of all the cells above was amplified by PCR, directly sequenced and analyzed by Sequence Navigatore and Factura software. The apoptosis rates were assayed by flow cytometry. Mitochondrial DNA mutation in the D-loop region was found in NB4 and APL cells before As2O3 use, but the mutation spots were remarkably increased after As2O3 treatment, which was positively correlated to the rates of cellular apoptosis, the correlation coefficient: rNB4-As2O3=0.973818, and rAPL-As2O3=0.934703. The mutation types include transition, transversion, codon insertion or deletion, and the mutation spots in all samples were not constant and regular. It is revealed that As2O3 aggravates mtDNA mutation in the D-loop region of acute promyelocytic leukemia cells both in vitro and in vivo. Mitochondrial DNA might be one of the targets of As2O3 in APL treatment.

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Computer modeling studies on the binding of 2',5'-linked dinucleoside phosphates to ribonuclease T1-influence of subsite interactions on the substrate specificity

The modes of binding of Gp(2',5')A, Gp(2',5')C, Gp(2',5')G and Gp(2',5')U to RNase T1 have been determined by computer modelling studies. All these dinucleoside phosphates assume extended conformations in the active site leading to better interactions with the enzyme. The 5'-terminal guanine of all these ligands is placed in the primary base binding site of the enzyme in an orientation similar to that of 2'-GMP in the RNase T1-2'-GMP complex. The 2'-terminal purines are placed close to the hydrophobic pocket formed by the residues Gly71, Ser72, Pro73 and Gly74 which occur in a loop region. However, the orientation of the 2'-terminal pyrimidines is different from that of 2'-terminal purines. This perhaps explains the higher binding affinity of the 2',5'-linked guanine dinucleoside phosphates with 2'-terminal purines than those with 2'-terminal pyrimidines. A comparison of the binding of the guanine dinucleoside phosphates with 2',5'- and 3',5'-linkages suggests significant differences in the ribose pucker and hydrogen bonding interactions between the catalytic residues and the bound nucleoside phosphate implying that 2',5'-linked dinucleoside phosphates may not be the ideal ligands to probe the role of the catalytic amino acid residues. A change in the amino acid sequence in the surface loop region formed by the residues Gly71 to Gly74 drastically affects the conformation of the base binding subsite, and this may account for the inactivity of the enzyme with altered sequence i.e., with Pro, Gly and Ser at positions 71 to 73 respectively. These results thus suggest that in addition to recognition and catalytic sites, interactions at the loop regions which constitute the subsite for base binding are also crucial in determining the substrate specificity.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Intraloop and Interloop Interactions in the Folded-G Quartet Structure of Oxytricha Telomeric Sesquence

The antiparallel intramolecular G quartet structure for the 3.5 copy Oxytricha telomeric sequence d(G(4)T(4))(3)G4 has been established using a combination of spectroscopic and chemical probing methods. In the presence of Naf ions, this sequence exhibits a circular dichroism spectrum with a positive band at 295 nm and a negative band around 265 nm, characteristic of an antiparallel G quartet structure. Further, we show that d(G(4)T(4))(3)G(4) adopts an antiparallel intramolecular G quartet structure even in K+ unlike d(G(4)T(4)G(4)). KMnO4 probing experiments indicated the existence of intra and interloop interactions in the Na+ induced structure. We have found that K+ not only increases the thermal stability of,G quartet structure but also binds to the loop region and disrupts stacking and interloop interactions. Biological consequences of such cation-dependent conformational micro-heterogeneity in the loop region of G quartet structures is also discussed.

Structural basis for the functional and inhibitory mechanisms of beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The beta-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of ``U''. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms. (C) 2011 Elsevier Inc. All rights reserved.

X-ray structure of a galactose-specific lectin from Spatholobous parviflorous

A galactose-specific seed lectin from Spatholobous parviflorus (SPL) has been purified, crystallized and its X-ray structure solved. It is the first lectin purified and crystallized from the genus Spatholobus (family: Fabaceae). The crystals belong to the space group P1, with a = 60.792 angstrom, b = 60.998 angstrom, c = 78.179 angstrom, alpha = 78.68 degrees, beta = 88.62 degrees, gamma = 104.32 degrees. The data were collected at 2.04 angstrom resolution under cryocondition, on a MAR image-plate detector system, mounted on a rotating anode X-ray generator. The coordinates of Dolichos biflorus lectin (1lu1) were successfully used for the structure solution by molecular replacement method. The primary structure of the SPL was not known earlier and it was unambiguously visible in the electron density. S. parviflorus lectin is a hetero-dimeric-tetramer with two alpha and two beta chains of 251 and 239 residues respectively. SPL has two metal ions, Ca(2+) and Mn(2+), bound to a loop region of each chain. The SPL monomers are in jelly roll form. (C) 2011 Elsevier B.V. All rights reserved.

Autoinhibitory mechanism and activity-related structural changes in a mycobacterial adenylyl cyclase

An adenylyl cyclase from Mycobacterium avium, Mal 120, is a functional orthologue of a pseudogene Rv1120c from Mycobacterium tuberculosis. We report the crystal structure of Mal 120 in a monomeric form and its truncated construct as a dimer. Mal 120 exists as a monomer in solution and crystallized as a monomer in the absence of substrate or inhibitor. An additional alpha-helix present at the N-terminus of the monomeric structure blocks the active site by interacting with the substrate binding residues and occupying the dimer interface region. However, the enzyme has been found to be active in solution, indicating the movement of the helix away from the interface to facilitate the formation of active dimers in conditions favourable for catalysis. Thus, the N-terminal helix of Ma1120 keeps the enzyme in an autoinhibited state when it is not active. Deletion of this helix enabled us to crystallize the molecule as an active homodimer in the presence of a P-site inhibitor 2',5'-dideoxy-3'-ATP, or pyrophosphate along with metal ions. The substrate specifying lysine residue plays a dual role of interacting with the substrate and stabilizing the dimer. The dimerization loop region harbouring the second substrate specifying residue, an aspartate, shows significant differences in conformation and position between the monomeric and dimeric structures. Thus, this study has not only revealed that significant structural transitions are required for the interconversion of the inactive and the active forms of the enzyme, but also provided precise nature of these transitions. (C) 2015 Elsevier Inc. All rights reserved.

Genetic diversity and population history of the red panda (Ailurus fulgens) as inferred from mitochondrial DNA sequence variations

The red panda (Ailurus fulgens) is one of the flagship species in worldwide conservation and is of special interest in evolutionary studies due to its taxonomic uniqueness. We sequenced a 236-bp fragment of the mitochondrial D-loop region in a sample of 5

Sequence variation and genetic diversity in the giant panda

About 336-444 bp mitochondrial D-loop region and tRNA gene were sequenced for 40 individuals of the giant panda which were collected from Mabian, Meigu, Yuexi, Baoxing, Pingwu, Qingchuan, Nanping and Baishuijiang, respectively. 9 haplotypes were found in 21 founders. The results showed that the giant panda has low genetic variations, and that there is no notable genetic isolation among geographical populations. The ancestor of the living giant panda population perhaps appeared in the late Pleistocene, and unfortunately, might have suffered bottle-neck attacks. Afterwards, its genetic diversity seemed to recover to same extent.

The origin and genetic diversity of Chinese native chicken breeds

The first 539 bases of mitochondrial DNA D-loop region of six Chinese native chicken breeds (Gallus gallus domesticus) were sequenced and compared to those of the red junglefowl (Gallus gallus), the gray junglefowl (Gallus sonneratii), the green junglefow

Sequence variations of the S7 ribosomal protein gene in primitive cyprinid fishes: Implication on phylogenetic analysis

Cyprinidae is the largest fish family in the world and contains about 210 genera and 2010 species. Appropriate DNA markers must be selected for the phylogenetic analyses of Cyprinidae. In present study, the 1st intron of the S7 ribosomal protein (r-protein) gene is first used to examine the relationships among cyprinid fishes. The length of the 1st intron obtained by PCR amplification ranges from 655 to 859 by in the 16 cyprinid species investigated, and is 602 by in Myxocyprinus asiaticus. Out of the alignment of 925 nucleotide sites obtained, the parsimony informative sites are 499 and occupy 54% of the total sites. The results indicate that the 1st intron sequences of the S7 r-protein gene in cyprinids are rich in informative sites and vary remarkably in sequence divergence from 2.3% between close species to 66.6% between distant species. The bootstrap values of the interior nodes in the NJ (neighbor-joining) and MP (most-parsimony) trees based on the present S7 r-protein gene data are higher than those based on cytochrome b and the d-loop region respectively. Therefore, the 1st intron sequences of the S7 r-protein gene in cyprinids are sensitive enough for phylogenetic analyses, and the 1st intron is an appropriate genetic marker for the phylogenetic reconstruction of the taxa in different cyprinid subfamilies. However, attempts to discuss whether the present S7 r-protein gene data can be applied to the phylogeny of the taxa at the level of the family or the higher categories in Cypriniformes need further studies.