Emergence of novel enzyme quasi-species depends on the substrate matrix

Current research on enzyme evolution has shown that many enzymes are promiscuous and have activities with alternative substrates. Mutagenesis tends to relax substrate selectivity, and evolving enzymes can be regarded (summed over evolutionary time) as clusters of enzyme variants, or “quasi-species,” tested against a “substrate matrix” defined by all chemical substances to which the evolvants are exposed.In this investigation, the importance of the substrate matrix for identification of evolvable clusters of enzymes was evaluated by random sampling of variants from a library of glutathione transferase (GST) mutants. The variant GSTs were created by DNA shuffling of homologous Alpha class sequences. The substrate matrix was an array of alternative substrates used under defined experimental conditions. The measured enzyme activities produced a rectangular matrix, in which the rows can be projected as enzyme vectors in substrate-activity space and, reciprocally, the columns can be projected as alternative substrate vectors in enzyme-activity space. Multivariate analysis of the catalytic activities demonstrated that the enzyme vectors formed two primary clusters or functional “molecular quasi-species.” These quasi-species serve as the raw material from which more specialized enzymes eventually could evolve. The substrate vectors similarly formed two major groups. Identification of separate quasi-species of GSTs in a mutant library was critically dependent on the nature of the substrate matrix. When substrates from just one of the two groups were used, only one cluster of enzymes could be recognized. On the other hand, expansion of the substrate matrix to include additional substrates showed the presence of a third quasi-species among the GST variants already analyzed. Thus, the portrayal of the functional quasi-species is intimately linked to the effective substrate matrix. In natural evolution, the substrates actually encountered therefore play a pivotal role in determining whether latent catalytic abilities become manifest in novel enzymes.     

Sequential binding of FurA from Anabaena sp. PCC 7120 to iron boxes: Exploring regulation at the nanoscale

Fur (ferric uptake regulator) proteins are involved in the control of a variety of processes in most prokaryotes. Although it is assumed that this regulator binds its DNA targets as a dimer, the way in which this interaction occurs remains unknown. We have focused on FurA from the cyanobacterium Anabaena sp. PCC 7120. To assess the molecular mechanism by which FurA specifically binds to “iron boxes” in P_furA, we examined the topology arrangement of FurA–DNA complexes by atomic force microscopy. Interestingly, FurA–P_furA complexes exhibit several populations, in which one is the predominant and depends clearly on the regulator/promoter ratio on the environment. Those results together with EMSA and other techniques suggest that FurA binds P_furA using a sequential mechanism: (i) a monomer specifically binds to an “iron box” and bends P_furA; (ii) two situations may occur, that a second FurA monomer covers the free “iron box" or that joins to the previously used forming a dimer which would maintain the DNA kinked; (iii) trimerization in which the DNA is unbent; and (iv) finally undergoes a tetramerization; the next coming molecules cover the DNA strands unspecifically. In summary, the bending appears when an “iron box” is bound to one or two molecules and decreases when both “iron boxes” are covered. These results suggest that DNA bending contributes at the first steps of FurA repression promoting the recruitment of new molecules resulting in a fine regulation in the Fur-dependent cluster associated genes.     

Exploring the Linkage Dependence of Polyubiquitin Conformations Using Molecular Modeling

Posttranslational modification of proteins by covalent attachment of a small protein ubiquitin (Ub) or a polymeric chain of Ub molecules (called polyubiquitin) is involved in controlling a vast variety of processes in eukaryotic cells. The question of how different polyubiquitin signals are recognized is central to understanding the specificity of various types of polyubiquitination. In polyubiquitin, monomers are linked to each other via an isopeptide bond between the C-terminal glycine of one Ub and a lysine of the other. The functional outcome of polyubiquitination depends on the particular lysine involved in chain formation and appears to rely on linkage-dependent conformation of polyubiquitin. Thus, K48-linked chains, a universal signal for proteasomal degradation, under physiological conditions adopt a closed conformation where functionally important residues L8, I44, and V70 are sequestered at the interface between two adjacent Ub monomers. By contrast, K63-linked chains, which act as a nonproteolytic regulatory signal, adopt an extended conformation that lacks hydrophobic interubiquitin contact. Little is known about the functional roles of the so-called “noncanonical” chains (linked via K6, K11, K27, K29, or K33, or linked head-to-tail), and no structural information on these chains is available, except for information on the crystal structure of the head-to-tail-linked diubiquitin (Ub₂). In this study, we use molecular modeling to examine whether any of the noncanonical chains can adopt a closed conformation similar to that in K48-linked polyubiquitin. Our results show that the eight possible Ub₂ chains can be divided into two groups: chains linked via K6, K11, K27, or K48 are predicted to form a closed conformation, whereas chains linked via K29, K33, or K63, or linked head-to-tail are unable to form such a contact due to steric occlusion. These predictions are validated by the known structures of K48-, K63-, and head-to-tail-linked chains. Our study also predicts structural models for Ub₂ chains linked via K6, K11, or K27. The implications of these findings for linkage-selective recognition of noncanonical polyubiquitin signals by various receptors are discussed.     

Genome wide sequence analysis grants unbiased definition of species boundaries in “Candidatus Phytoplasma”

The phytoplasmas are currently named using the Candidatus category, as the inability to grow them in vitro prevented (i) the performance of tests, such as DNA-DNA hybridization, that are regarded as necessary to establish species boundaries, and (ii) the deposition of type strains in culture collections. The recent accession to complete or nearly complete genome sequence information disclosed the opportunity to apply to the uncultivable phytoplasmas the same taxonomic approaches used for other bacteria. In this work, the genomes of 14 strains, belonging to the 16SrI, 16SrIII, 16SrV and 16SrX groups, including the species “Ca. P. asteris”, “Ca. P. mali”, “Ca. P. pyri”, “Ca. P. pruni”, and “Ca. P. australiense” were analyzed along with Acholeplasma laidlawi, to determine their taxonomic relatedness. Average nucleotide index (ANIm), tetranucleotide signature frequency correlation index (Tetra), and multilocus sequence analysis of 107 shared genes using both phylogenetic inference of concatenated (DNA and amino acid) sequences and consensus networks, were carried out. The results were in large agreement with the previously established 16S rDNA based classification schemes. Moreover, the taxonomic relationships within the 16SrI, 16SrIII and 16SrX groups, that represent clusters of strains whose relatedness could not be determined by 16SrDNA analysis, could be comparatively evaluated with non-subjective criteria. “Ca. P. mali” and “Ca. P. pyri” were found to meet the genome characteristics for the retention into two different, yet strictly related species; representatives of subgroups 16SrI-A and 16SrI-B were also found to meet the standards used in other bacteria to distinguish separate species; the genomes of the strains belonging to 16SrIII were found more closely related, suggesting that their subdivision into Candidatus species should be approached with caution.     

New Aspects of the Spontaneous Polymerization of Actin in the Presence of Salts

The mechanism of salt-induced actin polymerization involves the energetically unfavorable nucleation step, followed by filament elongation by the addition of monomers. The use of a bifunctional cross-linker, N,N′-(1,4-phenylene)dimaleimide, revealed rapid formation of the so-called lower dimers (LD) in which actin monomers are arranged in an antiparallel fashion. The filament elongation phase is characterized by a gradual LD decay and an increase in the yield of “upper dimers” (UD) characteristic of F-actin. Here we have used 90° light scattering, electron microscopy, and N,N′-(1,4-phenylene)dimaleimide cross-linking to reinvestigate relationships between changes in filament morphology, LD decay, and increase in the yield of UD during filament growth in a wide range of conditions influencing the rate of the nucleation reaction. The results show irregularity and instability of filaments at early stages of polymerization under all conditions used, and suggest that an earlier documented coassembling of LD with monomeric actin contributes to the initial disordering of the filaments rather than to the nucleation of polymerization. The effects of the type of G-actin-bound divalent cation (Ca²⁺/Mg²⁺), nucleotide (ATP/ADP), and polymerizing salt on the relation between changes in filament morphology and progress in G-actin-to-F-actin transformation show that ligand-dependent alterations in G-actin conformation determine not only the nucleation rate but also the kinetics of ordering of the filament structure in the elongation phase. The time courses of changes in the yield of UD suggest that filament maturation involves cooperative propagation of “proper” interprotomer contacts. Acceleration of this process by the initially bound MgATP supports the view that the filament-destabilizing conformational changes triggered by ATP hydrolysis and P_i liberation during polymerization are constrained by the intermolecular contacts established between MgATP monomers prior to ATP hydrolysis. An important role of contacts involving the DNase-I-binding loop and the C-terminus of actin is proposed.     

Packaging of a unit-length viral genome: the role of nucleotides and the gpD decoration protein in stable nucleocapsid assembly in bacteriophage lambda

The developmental pathways for a variety of eukaryotic and prokaryotic double-stranded DNA viruses include packaging of viral DNA into a preformed procapsid structure, catalyzed by terminase enzymes and fueled by ATP hydrolysis. In most instances, a capsid expansion process accompanies DNA packaging, which significantly increases the volume of the capsid to accommodate the full-length viral genome. “Decoration” proteins add to the surface of the expanded capsid lattice, and the terminase motors tightly package DNA, generating up to ∼ 20 atm of internal capsid pressure. Herein we describe biochemical studies on genome packaging using bacteriophage λ as a model system. Kinetic analysis suggests that the packaging motor possesses at least four ATPase catalytic sites that act cooperatively to effect DNA translocation, and that the motor is highly processive. While not required for DNA translocation into the capsid, the phage λ capsid decoration protein gpD is essential for the packaging of the penultimate 8-10 kb (15-20%) of the viral genome; virtually no DNA is packaged in the absence of gpD when large DNA substrates are used, most likely due to a loss of capsid structural integrity. Finally, we show that ATP hydrolysis is required to retain the genome in a packaged state subsequent to condensation within the capsid. Presumably, the packaging motor continues to “idle” at the genome end and to maintain a positive pressure towards the packaged state. Surprisingly, ADP, guanosine triphosphate, and the nonhydrolyzable ATP analog 5'-adenylyl-beta,gamma-imidodiphosphate (AMP-PNP) similarly stabilize the packaged viral genome despite the fact that they fail to support genome packaging. In contrast, the poorly hydrolyzed ATP analog ATP-γS only partially stabilizes the nucleocapsid, and a DNA is released in “quantized” steps. We interpret the ensemble of data to indicate that (i) the viral procapsid possesses a degree of plasticity that is required to accommodate the packaging of large DNA substrates; (ii) the gpD decoration protein is required to stabilize the fully expanded capsid; and (iii) nucleotides regulate high-affinity DNA binding interactions that are required to maintain DNA in the packaged state.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

DGCA: A comprehensive R package for Differential Gene Correlation Analysis

Background

Dissecting the regulatory relationships between genes is a critical step towards building accurate predictive models of biological systems. A powerful approach towards this end is to systematically study the differences in correlation between gene pairs in more than one distinct condition.

Results

In this study we develop an R package, DGCA (for Differential Gene Correlation Analysis), which offers a suite of tools for computing and analyzing differential correlations between gene pairs across multiple conditions. To minimize parametric assumptions, DGCA computes empirical p-values via permutation testing. To understand differential correlations at a systems level, DGCA performs higher-order analyses such as measuring the average difference in correlation and multiscale clustering analysis of differential correlation networks. Through a simulation study, we show that the straightforward z-score based method that DGCA employs significantly outperforms the existing alternative methods for calculating differential correlation. Application of DGCA to the TCGA RNA-seq data in breast cancer not only identifies key changes in the regulatory relationships between TP53 and PTEN and their target genes in the presence of inactivating mutations, but also reveals an immune-related differential correlation module that is specific to triple negative breast cancer (TNBC).

Conclusions

DGCA is an R package for systematically assessing the difference in gene-gene regulatory relationships under different conditions. This user-friendly, effective, and comprehensive software tool will greatly facilitate the application of differential correlation analysis in many biological studies and thus will help identification of novel signaling pathways, biomarkers, and targets in complex biological systems and diseases.

     

Codon sextets with leading role of serine create “ideal” symmetry classification scheme of the genetic code

The standard classification scheme of the genetic code is organized for alphabetic ordering of nucleotides. Here we introduce the new, “ideal” classification scheme in compact form, for the first time generated by codon sextets encoding Ser, Arg and Leu amino acids. The new scheme creates the known purine/pyrimidine, codon–anticodon, and amino/keto type symmetries and a novel A + U rich/C + G rich symmetry. This scheme is built from “leading” and “nonleading” groups of 32 codons each. In the ensuing 4 × 16 scheme, based on trinucleotide quadruplets, Ser has a central role as initial generator. Six codons encoding Ser and six encoding Arg extend continuously along a linear array in the “leading” group, and together with four of six Leu codons uniquely define construction of the “leading” group. The remaining two Leu codons enable construction of the “nonleading” group. The “ideal” genetic code suggests the evolution of genetic code with serine as an initiator.     

Ampelozyziphus amazonicus Ducke (Rhamnaceae), a medicinal plant used to prevent malaria in the Amazon Region, hampers the development of Plasmodium berghei sporozoites

Most medicinal plants used against malaria in endemic areas aim to treat the acute symptoms of the disease such as high temperature fevers with periodicity and chills. In some endemic areas of the Brazilian Amazon region one medicinal plant seems to be an exception: Ampelozyziphus amazonicus, locally named “Indian beer” or “Saracura-mira”, used to prevent the disease when taken daily as a cold suspension of powdered dried roots. In previous work we found no activity of the plant extracts against malaria blood parasites in experimentally infected animals (mice and chickens) or in cultures of Plasmodium falciparum. However, in infections induced by sporozoites, chickens treated with plant extracts were partially protected against Plasmodium gallinaceum and showed reduced numbers of exoerythrocytic forms in the brain. We now present stronger evidence that the ethanolic extract of “Indian beer” roots hampers in vitro and in vivo development of Plasmodium berghei sporozoites, a rodent malaria parasite. Some mice treated with high doses of the plant extract did not become infected after sporozoite inoculation, whereas others had a delayed prepatent period and lower parasitemia. Our data validates the use of “Indian beer” as a remedy for malaria prophylaxis in the Amazon, where the plant exists and the disease represents an important problem which is difficult to control. Studies aiming to identify the active compounds responsible for the herein described causal prophylactic activity are needed and may lead to a new antimalarial prophylactic.     

The Capsid Proteins of a Large, Icosahedral dsDNA Virus

Chilo iridescent virus (CIV) is a large (∼ 1850 Å diameter) insect virus with an icosahedral, T = 147 capsid, a double-stranded DNA (dsDNA) genome, and an internal lipid membrane. The structure of CIV was determined to 13 Å resolution by means of cryoelectron microscopy (cryoEM) and three-dimensional image reconstruction. A homology model of P50, the CIV major capsid protein (MCP), was built based on its amino acid sequence and the structure of the homologous Paramecium bursaria chlorella virus 1 Vp54 MCP. This model was fitted into the cryoEM density for each of the 25 trimeric CIV capsomers per icosahedral asymmetric unit. A difference map, in which the fitted CIV MCP capsomers were subtracted from the CIV cryoEM reconstruction, showed that there are at least three different types of minor capsid proteins associated with the capsomers outside the lipid membrane. “Finger” proteins are situated at many, but not all, of the spaces between three adjacent capsomers within each trisymmetron, and “zip” proteins are situated between sets of three adjacent capsomers at the boundary between neighboring trisymmetrons and pentasymmetrons. Based on the results of segmentation and density correlations, there are at least eight finger proteins and three dimeric and two monomeric zip proteins in one asymmetric unit of the CIV capsid. These minor proteins appear to stabilize the virus by acting as intercapsomer cross-links. One transmembrane “anchor” protein per icosahedral asymmetric unit, which extends from beneath one of the capsomers in the pentasymmetron to the internal leaflet of the lipid membrane, may provide additional stabilization for the capsid. These results are consistent with the observations for other large, icosahedral dsDNA viruses that also utilize minor capsid proteins for stabilization and for determining their assembly.     

The Crystal Structure of the Hinge Domain of the Escherichia coli Structural Maintenance of Chromosomes Protein MukB

MukB, a divergent structural maintenance of chromosomes (SMC) protein, is important for chromosomal segregation and condensation in γ-proteobacteria. MukB and canonical SMC proteins share a characteristic five-domain structure. Globular N- and C-terminal domains interact to form an ATP-binding cassette-like ATPase or “head” domain, which is connected to a smaller dimerization or “hinge” domain by a long, antiparallel coiled coil. In addition to mediating dimerization, this hinge region has been implicated in both conformational flexibility and dynamic protein-DNA interactions. We report here the first crystallographic model of the MukB hinge domain. This model also contains approximately 20% of the coiled-coil domain, including an unusual coiled-coil deviation. These results will facilitate studies to clarify the roles of both the hinge and the coiled-coil domains in MukB function.     

Evolutionary genomics of Culex pipiens: global and local adaptations associated with climate,life-history traits and anthropogenic factors

We present the first genome-wide study of recent evolution in Culex pipiens species complex focusing on the genomic extent, functional targets and likely causes of global and local adaptations. We resequenced pooled samples of six populations of C. pipiens and two populations of the outgroup Culex torrentium. We used principal component analysis to systematically study differential natural selection across populations and developed a phylogenetic scanning method to analyse admixture without haplotype data. We found evidence for the prominent role of geographical distribution in shaping population structure and specifying patterns of genomic selection. Multiple adaptive events, involving genes implicated with autogeny, diapause and insecticide resistance were limited to specific populations. We estimate that about 5–20% of the genes (including several histone genes) and almost half of the annotated pathways were undergoing selective sweeps in each population. The high occurrence of sweeps in non-genic regions and in chromatin remodelling genes indicated the adaptive importance of gene expression changes. We hypothesize that global adaptive processes in the C. pipiens complex are potentially associated with South to North range expansion, requiring adjustments in chromatin conformation. Strong local signature of adaptation and emergence of hybrid bridge vectors necessitate genomic assessment of populations before specifying control agents.     

Cytochrome c1 exhibits two binding sites for cytochrome c in plants

In plants, channeling of cytochrome c molecules between complexes III and IV has been purported to shuttle electrons within the supercomplexes instead of carrying electrons by random diffusion across the intermembrane bulk phase. However, the mode plant cytochrome c behaves inside a supercomplex such as the respirasome, formed by complexes I, III and IV, remains obscure from a structural point of view. Here, we report ab-initio Brownian dynamics calculations and nuclear magnetic resonance-driven docking computations showing two binding sites for plant cytochrome c at the head soluble domain of plant cytochrome c₁, namely a non-productive (or distal) site with a long heme-to-heme distance and a functional (or proximal) site with the two heme groups close enough as to allow electron transfer. As inferred from isothermal titration calorimetry experiments, the two binding sites exhibit different equilibrium dissociation constants, for both reduced and oxidized species, that are all within the micromolar range, thus revealing the transient nature of such a respiratory complex. Although the docking of cytochrome c at the distal site occurs at the interface between cytochrome c₁ and the Rieske subunit, it is fully compatible with the complex III structure. In our model, the extra distal site in complex III could indeed facilitate the functional cytochrome c channeling towards complex IV by building a “floating boat bridge” of cytochrome c molecules (between complexes III and IV) in plant respirasome.     

A novel fluorescent timer based on bicistronic expression strategy in Caenorhabditis elegans

Fluorescent timers are useful tools for studying the spatial and temporal cellular or molecular events. Based on the trans-splicing mechanism in Caenorhabditis elegans, we constructed a “fluorescent timer” through bicistronic expression of two fluorescent proteins with different maturation times. When used in vivo, this “timer” changes its color over time and therefore can be used to monitor the activity of the targeted promoters in C. elegans. Using this “timer”, we have successfully traced the time-dependent activity of myo-3 promoter which drives expression in body wall muscle and vulval muscle. We found that the myo-3 promoter started to be active about 7 h after egg-laying and sustained its activity in the following hatching process. We have also determined the myo-3 promoter activity during larval development by this “timer”. We anticipate that more new “fluorescent timers” with variable time-resolution could be designed by bicistronic expression of different fluorescent protein pairs.     

A systematic comparison of quantitative high-resolution DNA methylation analysis and methylation-specific PCR

Assessment of DNA methylation has become a critical factor for the identification, development and application of methylation based biomarkers. Here we describe a systematic comparison of a quantitative high-resolution mass spectrometry-based approach (MassARRAY), pyrosequencing and the broadly used methylation-specific PCR (MSP) technique analyzing clinically relevant epigenetically silenced genes in acute myeloid leukemia (AML). By MassARRAY and pyrosequencing, we identified significant DNA methylation differences at the ID4 gene promoter and in the 5′ region of members of the SFRP gene family in 62 AML patients compared with healthy controls. We found a good correlation between data obtained by MassARRAY and pyrosequencing (correlation coefficient R² = 0.88). MSP-based assessment of the identical samples showed less pronounced differences between AML patients and controls. By direct comparison of MSP-derived and MassARRAY-based methylation data as well as pyrosequencing, we could determine overestimation of DNA methylation data by MSP. We found sequence-context dependent highly variable cut-off values of quantitative DNA methylation values serving as discriminator for the two MSP methylation categories. Moreover, good agreements between quantitative methods and MSP could not be achieved for all investigated loci. Significant correlation of the quantitative assessment but not of MSP-derived methylation data with clinically important characteristics in our patient cohort demonstrated clinical relevance of quantitative DNA methylation assessment. Taken together, while MSP is still the most commonly applied technique for DNA methylation assessment, our data highlight advantages of quantitative approaches for precise characterization and reliable biomarker use of aberrant DNA methylation in primary patient samples, particularly.     

The role of epibotic bacteria from the surface of the soft coral Dendronephthya sp. in the inhibition of larval settlement

It has been suggested that bacteria associated with soft-bodied organisms are suggested to produce bioactive compounds against the attachment of invertebrate larvae and bacteria onto the surface of these organisms. Our recent study has demonstrated that epibiotic bacteria from the surface of the soft coral Dendronephthya sp. (Coelenterata: Octocoralia, Alcyonacea) inhibit the growth of bacteria commonly found in marine natural biofilms. In the present study, the effect of 11 epibiotic bacteria isolated from the surface of Dendronephthya sp. on larval settlement of the tubeworms Hydroides elegans was examined using laboratory bioassay. Among 11 bacterial isolates, 2 strains (18%) inhibited the larval settlement of H. elegans (Haswell), 4 strains (36%) were “inductive” to larvae and the remaining 5 strains (46%) were “non-inductive”. There was no correlation between the antifouling activities of bacterial isolates and their phylogenetic origin, i.e. closely related bacterial strains showed different effects on larval settlement of H. elegans. When all “inductive”, “non-inductive” and “inhibitive” bacterial isolates were mixed in a 1:1:1 ratio, the effect of the resultant multispecies film on larval settlement became “inhibitive”. Waterborne compounds of Vibrio sp. and an unidentified α-Proteobacterium, which suppressed the settlement of H. elegans and Bugula neritina (L.) larvae, were further investigated using size fractionation and bioassay-guided enzymatic analysis. It was found that antilarval settlement compounds from these bacteria were heat-stable polysaccharides with a molecular weight >100 kDa. The results indicate that the bacteria associated with the soft coral Dendronephthya sp. may contribute to the antifouling mechanisms of the soft-bodied organisms by producing compounds that are against bacterial growth and settlement of macrofoulers on the surface of their host.     

Phase diagram of stigmasterol-dipalmitoylphosphatidylcholine mixtures dispersed in excess water

As a simple model of rafts in plant cells, the effect of stigmasterol, one of the predominant sterols in plant plasma membranes, on the phase behavior of dipalmitoylphosphatidylcholine (DPPC) multilayers has been studied by X-ray diffraction (XRD), differential scanning calorimetry (DSC), and freeze-fracture electron microscopy (FFEM) techniques. A partial phase diagram of the binary system has been constructed. Particularly, the stigmasterol concentrations of the “left endpoint” and “right endpoint” of the three-phase line have been determined using the newly developed linear and nonlinear fitting method. They are 6.2 and 23.7 mol%, respectively. Furthermore, the resemblance and difference of phase diagrams of DPPC/stigmasterol, DPPC/cholesterol, and DPPC/ergosterol have been compared and the efficiency of these sterols in promoting the formation of the liquid-ordered domains (rafts) have also been discussed.