Modelling the Efficiency of Codon–tRNA Interactions Based on Codon Usage Bias

The tRNA adaptation index (tAI) is a widely used measure of the efficiency by which a coding sequence is recognized by the intra-cellular tRNA pool. This index includes among others weights that represent wobble interactions between codons and tRNA molecules. Currently, these weights are based only on the gene expression in Saccharomyces cerevisiae. However, the efficiencies of the different codon–tRNA interactions are expected to vary among different organisms. In this study, we suggest a new approach for adjusting the tAI weights to any target model organism without the need for gene expression measurements. Our method is based on optimizing the correlation between the tAI and a measure of codon usage bias. Here, we show that in non-fungal the new tAI weights predict protein abundance significantly better than the traditional tAI weights. The unique tRNA–codon adaptation weights computed for 100 different organisms exhibit a significant correlation with evolutionary distance. The reported results demonstrate the usefulness of the new measure in future genomic studies.     

Comparative analysis of bacterial genomes: identification of divergent regions in mycobacterial strains using an anchor-based approach

Comparative genomic approaches are useful in identifying molecular differences between organisms. Currently available methods fail to identify small changes in genomes, such as expansion of short repetitive motifs and to analyse divergent sequences. In this report, we describe an anchor-based whole genome comparison (ABWGC) method. ABWGC is based on random sampling of anchor sequences from one genome, followed by analysis of sampled and homologous regions from the target genome. The method was applied to compare two strains of Mycobacterium tuberculosis CDC1551 and H37Rv. ABWGC was able to identify a total of 104 indels including 20 expansion of short repetitive sequences and five recombination events. It included 18 new unidentified genomic differences. ABWGC also identified 188 SNPs including eight new ones. The method was also used to compare M. tuberculosis H37Rv and M. avium genomes. ABWGC was able to correctly pick 1002 additional indels (size>100nt) between the two organisms in contrast to MUMmer, a popular tool for comparative genomics. ABWGC was able to identify correctly repeat expansion and indels in a set of simulated sequences. The study also revealed important role of small repeat expansion in the evolution of M. tuberculosis strains.     

Proposing a Radial Basis Function and CSDM Indices to Predict the Traumatic Brain Injury Risk

BackgroundMany head injury indices and finite element (FE) head models have been proposed to predict traumatic brain injury (TBI). Although FE head models are suitable methods with high accuracy, they are computationally intensive. Head motion-based brain injury criteria are usually fast tools with lower accuracy. So, the objective of this study is to propose new criteria along with an artificial neural network model to predict TBI risks, which can be fast and accurate.MethodsFor this purpose, 250 FE head simulations have been carried out at 5 magnitudes and 50 rotational impact directions using the SIMon model. The effects of directions and magnitudes of rotational impacts were assessed for cumulative strain damage measure (CSDM) values. Next, statistical analysis and neural network were applied to predict CSDM values.ResultsThe results of the present research showed that the direction of rotation in the sagittal and frontal planes had a considerable effect on the CSDM values. Furthermore, new brain injury indices and a radial basis function neural network have been proposed to predict CSDM values which having high correlation coefficients with SIMon responses.ConclusionsThe results of this research demonstrated that rotational impact directions should be used to develop new head injury criteria being able to predict CSDM values. However, findings of present research proved that head motion-based brain injury criteria and RBF network can be used to predict FE head model responses with high speed and accuracy.     

Genome-wide identification and characterization of replication origins by deep sequencing

Background

DNA replication initiates at distinct origins in eukaryotic genomes, but the genomic features that define these sites are not well understood.

Results

We have taken a combined experimental and bioinformatic approach to identify and characterize origins of replication in three distantly related fission yeasts: Schizosaccharomyces pombe, Schizosaccharomyces octosporus and Schizosaccharomyces japonicus. Using single-molecule deep sequencing to construct amplification-free high-resolution replication profiles, we located origins and identified sequence motifs that predict origin function. We then mapped nucleosome occupancy by deep sequencing of mononucleosomal DNA from the corresponding species, finding that origins tend to occupy nucleosome-depleted regions.

Conclusions

The sequences that specify origins are evolutionarily plastic, with low complexity nucleosome-excluding sequences functioning in S. pombe and S. octosporus, and binding sites for trans-acting nucleosome-excluding proteins functioning in S. japonicus. Furthermore, chromosome-scale variation in replication timing is conserved independently of origin location and via a mechanism distinct from known heterochromatic effects on origin function. These results are consistent with a model in which origins are simply the nucleosome-depleted regions of the genome with the highest affinity for the origin recognition complex. This approach provides a general strategy for understanding the mechanisms that define DNA replication origins in eukaryotes.     

Four chromosome replication origins in the archaeon Pyrobaculum calidifontis

Replication origins were mapped in hyperthermophilic crenarchaea, using high‐throughput sequencing‐based marker frequency analysis. We confirm previous origin mapping in Sulfolobus acidocaldarius, and demonstrate that the single chromosome of Pyrobaculum calidifontis contains four replication origins, the highest number detected in a prokaryotic organism. The relative positions of the origins in both organisms coincided with regions enriched in highly conserved (core) archaeal genes. We show that core gene distribution provides a useful tool for origin identification in archaea, and predict multiple replication origins in a range of species. One of the P. calidifontis origins was mapped in detail, and electrophoretic mobility shift assays demonstrated binding of the Cdc6/Orc1 replication initiator protein to a repeated sequence element, denoted Orb‐1, within the origin. The high‐throughput sequencing approach also allowed for an annotation update of both genomes, resulting in the restoration of open reading frames encoding proteins involved in, e.g., sugar, nitrate and energy metabolism, as well as in glycosylation and DNA repair.     

Quantitative analysis of the mechanism of DNA binding by Bacillus DnaA protein

DnaA protein has the sole responsibility of initiating a new round of DNA replication in prokaryotic organisms. It recognizes the origin of DNA replication, and initiates chromosomal DNA replication in the bacterial genome. In Gram-negative Escherichia coli, a large number of DnaA molecules bind to specific DNA sequences (known as DnaA boxes) in the origin of DNA replication, oriC, leading to the activation of the origin. We have cloned, expressed, and purified full-length DnaA protein in large quantity from Gram-positive pathogen Bacillus anthracis (DnaA_BA). DnaA_BA was a highly soluble monomeric protein making it amenable to quantitative analysis of its origin recognition mechanisms. DnaA_BA bound DnaA boxes with widely divergent affinities in sequence and ATP-dependent manner. In the presence of ATP, the K_D ranged from 3.8 × 10⁻⁸ M for a specific DnaA box sequence to 4.1 × 10⁻⁷ M for a non-specific DNA sequence and decreased significantly in the presence of ADP. Thermodynamic analyses of temperature and salt dependence of DNA binding indicated that hydrophobic (entropic) and ionic bonds contributed to the DnaA_BA·DNA complex formation. DnaA_BA had a DNA-dependent ATPase activity. DNA sequences acted as positive effectors and modulated the rate (V_max) of ATP hydrolysis without any significant change in ATP binding affinity.     

Excess SeqA leads to replication arrest and a cell division defect in Vibrio cholerae

Although most bacteria contain a single circular chromosome, some have complex genomes, and all Vibrio species studied so far contain both a large and a small chromosome. In recent years, the divided genome of Vibrio cholerae has proven to be an interesting model system with both parallels to and novel features compared with the genome of Escherichia coli. While factors influencing the replication and segregation of both chromosomes have begun to be elucidated, much remains to be learned about the maintenance of this genome and of complex bacterial genomes generally. An important aspect of replicating any genome is the correct timing of initiation, without which organisms risk aneuploidy. During DNA replication in E. coli, newly replicated origins cannot immediately reinitiate because they undergo sequestration by the SeqA protein, which binds hemimethylated origin DNA. This DNA is already methylated by Dam on the template strand and later becomes fully methylated; aberrant amounts of Dam or the deletion of seqA leads to asynchronous replication. In our study, hemimethylated DNA was detected at both origins of V. cholerae, suggesting that these origins are also subject to sequestration. The overproduction of SeqA led to a loss of viability, the condensation of DNA, and a filamentous morphology. Cells with abnormal DNA content arose in the population, and replication was inhibited as determined by a reduced ratio of origin to terminus DNA in SeqA-overexpressing cells. Thus, excessive SeqA negatively affects replication in V. cholerae and prevents correct progression to downstream cell cycle events such as segregation and cell division.     

Creating a Novel Origin of Replication through Modulating DNA-Protein Interfaces

Background

While the molecular mechanisms of DNA-protein specificity at the origin of replication have been determined in many model organisms, these interactions remain unknown in the majority of higher eukaryotes and numerous vertebrate viruses. Similar to many viral origins of replication, adeno-associated virus (AAV) utilizes a cis-acting origin of replication and a virus specific Replication protein (Rep) to faithfully carry out self-priming replication. The mechanisms of AAV DNA replication are generally well understood. However, the molecular basis of specificity between the Rep protein and the viral origin of replication between different AAV serotypes remains uncharacterized.

Methodology/Principal Findings

By generating a panel of chimeric and mutant origins between two AAV serotypes, we have mapped two independent DNA-Protein interfaces involved in replicative specificity. In vivo replication assays and structural modeling demonstrated that three residues in the AAV2 Rep active site are necessary to cleave its cognate origin. An analogous origin (AAV5) possesses a unique interaction between an extended Rep binding element and a 49 aa region of Rep containing two DNA binding interfaces.

Conclusions/Significance

The elucidation of these structure-function relationships at the AAV origin led to the creation of a unique recombinant origin and compatible Rep protein with properties independent of either parent serotype. This novel origin may impact the safety and efficacy of AAV as a gene delivery tool. This work may also explain the unique ability of certain AAV serotypes to achieve site-directed integration into the human chromosome. Finally, this result impacts the study of conserved DNA viruses which employ rolling circle mechanisms of replication.     

Functional centromeres determine the activation time of pericentric origins of DNA replication in Saccharomyces cerevisiae

The centromeric regions of all Saccharomyces cerevisiae chromosomes are found in early replicating domains, a property conserved among centromeres in fungi and some higher eukaryotes. Surprisingly, little is known about the biological significance or the mechanism of early centromere replication; however, the extensive conservation suggests that it is important for chromosome maintenance. Do centromeres ensure their early replication by promoting early activation of nearby origins, or have they migrated over evolutionary time to reside in early replicating regions? In Candida albicans, a neocentromere contains an early firing origin, supporting the first hypothesis but not addressing whether the new origin is intrinsically early firing or whether the centromere influences replication time. Because the activation time of individual origins is not an intrinsic property of S. cerevisiae origins, but is influenced by surrounding sequences, we sought to test the hypothesis that centromeres influence replication time by moving a centromere to a late replication domain. We used a modified Meselson-Stahl density transfer assay to measure the kinetics of replication for regions of chromosome XIV in which either the functional centromere or a point-mutated version had been moved near origins that reside in a late replication region. We show that a functional centromere acts in cis over a distance as great as 19 kb to advance the initiation time of origins. Our results constitute a direct link between establishment of the kinetochore and the replication initiation machinery, and suggest that the proposed higher-order structure of the pericentric chromatin influences replication initiation.     

Universal Temporal Profile of Replication Origin Activation in Eukaryotes

Although replication proteins are conserved among eukaryotes, the sequence requirements for replication initiation differ between species. In all species, however, replication origins fire asynchronously throughout S phase. The temporal program of origin firing is reproducible in cell populations but largely probabilistic at the single-cell level. The mechanisms and the significance of this program are unclear. Replication timing has been correlated with gene activity in metazoans but not in yeast. One potential role for a temporal regulation of origin firing is to minimize fluctuations in replication end time and avoid persistence of unreplicated DNA in mitosis. Here, we have extracted the population-averaged temporal profiles of replication initiation rates for S. cerevisiae, S. pombe, D. melanogaster, X. laevis and H. sapiens from genome-wide replication timing and DNA combing data. All the profiles have a strikingly similar shape, increasing during the first half of S phase then decreasing before its end. A previously proposed minimal model of stochastic initiation modulated by accumulation of a recyclable, limiting replication-fork factor and fork-promoted initiation of new origins, quantitatively described the observed profiles without requiring new implementations.The selective pressure for timely completion of genome replication and optimal usage of replication proteins that must be imported into the cell nucleus can explain the generic shape of the profiles. We have identified a universal behavior of eukaryotic replication initiation that transcends the mechanisms of origin specification. The population-averaged efficiency of replication origin usage changes during S phase in a strikingly similar manner in a highly diverse set of eukaryotes. The quantitative model previously proposed for origin activation in X. laevis can be generalized to explain this evolutionary conservation.     

Long-term prediction on Baltic fish stocks based on periodicity of solar activity

The methods presently in wide use for the assessment of marine ecosystems and fish stocks cannot provide the high-quality long-term prognoses urgently needed for improved management of marine ecosystems and their fishery resources. A novel method of forecasting the long-term qualitative composition of ecosystem and fish stocks in the Baltic Sea based on the periodicity and succession of the processes: extraterrestrial factors????climate changes????regime forming factors (chiefly salinity and temperature)????ecosystems????fish stocks was used in 1989 to predict the approximate year class abundances of cod (Gadus morhua callarias), sprat (Sprattus sprattus balticus), sea herring and gulf herring (Clupea harengus membras) in the Baltic Sea for the years 1989?C2008. A comparison with actual development up to 2008 show that this prognosis was confirmed to a considerable extent. The method based on periodical fluctuations of freshwater input and the thermal regime is described. It predicts a new regime shift in the 2020s to a higher salinity and more acceptable conditions for the organisms of marine origin.     

Deep learning allows genome-scale prediction of Michaelis constants from structural features

The Michaelis constant K_M describes the affinity of an enzyme for a specific substrate and is a central parameter in studies of enzyme kinetics and cellular physiology. As measurements of K_M are often difficult and time-consuming, experimental estimates exist for only a minority of enzyme–substrate combinations even in model organisms. Here, we build and train an organism-independent model that successfully predicts K_M values for natural enzyme–substrate combinations using machine and deep learning methods. Predictions are based on a task-specific molecular fingerprint of the substrate, generated using a graph neural network, and on a deep numerical representation of the enzyme’s amino acid sequence. We provide genome-scale K_M predictions for 47 model organisms, which can be used to approximately relate metabolite concentrations to cellular physiology and to aid in the parameterization of kinetic models of cellular metabolism.

To understand the action of an enzyme, we need to know its affinity for its substrates, quantified by Michaelis constants, but these are difficult to measure experimentally. This study shows that a deep learning model that can predict them from structural features of the enzyme and substrate, providing KM predictions for all enzymes across 47 model organisms.     

Functional Characterization of the Origin of Replication of pRN1 from Sulfolobus islandicus REN1H1

Plasmid pRN1 from Sulfolobus islandicus REN1H1 is believed to replicate by a rolling circle mechanism but its origin and mechanism of replication are not well understood. We sought to create minimal expression vectors based on pRN1 that would be useful for heterologous gene expression in S. acidocaldarius, and in the process improve our understanding of the mechanism of replication. We constructed and transformed shuttle vectors that harbored different contiguous stretches of DNA from pRN1 into S. acidocaldarius E4-39, a uracil auxotroph. A 232-bp region 3’ of orf904 was found to be critical for pRN1 replication and is therefore proposed to be the putative origin of replication. This 232-bp region contains a 100-bp stem-loop structure believed to be the double-strand origin of replication. The loop of the 100-bp structure contains a GTG tri-nucleotide motif, a feature that was previously reported to be important for the primase activity of Orf904. This putative origin and the associated orf56 and orf904 were identified as the minimal replicon of pRN1 because transformants of plasmids lacking any of these three features were not recovered. Plasmids lacking orf904 and orf56 but harboring the putative origin were transformable when orf904 and orf56 were provided in-trans; a 75-bp region 5’ of the orf904 start codon was found to be essential for this complementation. Detailed knowledge of the pRN1 origin of replication will broaden the application of the plasmid as a genetic tool for Sulfolobus species.     

Bacteriophage SPP1 DNA replication strategies promote viral and disable host replication in vitro

Complex viruses that encode their own initiation proteins and subvert the host’s elongation apparatus have provided valuable insights into DNA replication. Using purified bacteriophage SPP1 and Bacillus subtilis proteins, we have reconstituted a rolling circle replication system that recapitulates genetically defined protein requirements. Eleven proteins are required: phage-encoded helicase (G40P), helicase loader (G39P), origin binding protein (G38P) and G36P single-stranded DNA-binding protein (SSB); and host-encoded PolC and DnaE polymerases, processivity factor (β₂), clamp loader (τ-δ-δ′) and primase (DnaG). This study revealed a new role for the SPP1 origin binding protein. In the presence of SSB, it is required for initiation on replication forks that lack origin sequences, mimicking the activity of the PriA replication restart protein in bacteria. The SPP1 replisome is supported by both host and viral SSBs, but phage SSB is unable to support B. subtilis replication, likely owing to its inability to stimulate the PolC holoenzyme in the B. subtilis context. Moreover, phage SSB inhibits host replication, defining a new mechanism by which bacterial replication could be regulated by a viral factor.     

Progress in understanding DNA replication control

Completion of genome duplication during the S-phase of the cell cycle is crucial for the maintenance of genomic integrity. In eukaryotes, chromosomal DNA replication is accomplished by the activity of multiple origins of DNA replication scattered across the genome. Origin specification, selection and activity as well as the availability of replication factors and the regulation of DNA replication licensing, have unique and common features among eukaryotes. Although the initial studies on the semiconservative nature of chromosome duplication were carried out in the mid 1950s in Vicia faba, since then plant DNA replication studies have been scarce. However, they have received an unprecedented drive in the last decade after the completion of sequencing the Arabidopsis thaliana genome, and more recently of other plant genomes. In particular, the past year has witnessed major advances with the use of genomic approaches to study chromosomal replication timing, DNA replication origins and licensing control mechanisms. In this minireview article we discuss these recent discoveries in plants in the context of what is known at the genomic level in other eukaryotes. These studies constitute the basis for addressing in the future key questions about replication origin specification and function that will be of relevance not only for plants but also for the rest of multicellular organisms.     

Deoxyribonucleic acid replication time in Mycobacterium tuberculosis H37 Rv

The DNA increment method, designed for measuring the increment in the amount of DNA after inhibition of initiation of fresh rounds of replication initiation was employed to measure the rate of deoxyribonucleic acid (DNA) chain growth in Mycobacterium tuberculosis H37Rv growing in Youman and Karlson's medium at 37°C with a generation time of 24 h and also in relatively fast growing species like Mycobacterium smegmatis and Escherichia coli. From the results obtained, the time required for a DNA replication fork to traverse the chromosome from origin to terminus (C period) was calculated. The chain elongation rates of DNA of the three organisms was determined from the C period and the known genome sizes assuming that all these genomes have a single replication origin and bidirectional replication fork. The rate for M. tuberculosis was 3,200 nucleotides per min about 11 times slower than that of M. smegmatis and about 13–18 times slower than that of E. coli.Abbreviations DNA deoxyribonucleic acid - td delay in initiation - OD optical density - CAM chloramphenicol - RIF rifampicin     

Assembly of Helicobacter pylori Initiation Complex Is Determined by Sequence-Specific and Topology-Sensitive DnaA–oriC Interactions

In bacteria, chromosome replication is initiated by binding of the DnaA initiator protein to DnaA boxes located in the origin of chromosomal replication (oriC). This leads to DNA helix opening within the DNA-unwinding element. Helicobacter pylori oriC, the first bipartite origin identified in Gram-negative bacteria, contains two subregions, oriC1 and oriC2, flanking the dnaA gene. The DNA-unwinding element region is localized in the oriC2 subregion downstream of dnaA. Surprisingly, oriC2–DnaA interactions were shown to depend on DNA topology, which is unusual in bacteria but is similar to initiator–origin interactions observed in higher organisms. In this work, we identified three DnaA boxes in the oriC2 subregion, two of which were bound only as supercoiled DNA. We found that all three DnaA boxes play important roles in orisome assembly and subsequent DNA unwinding, but different functions can be assigned to individual boxes. This suggests that the H. pylori oriC may be functionally divided, similar to what was described recently for Escherichia coli oriC. On the basis of these results, we propose a model of initiation complex formation in H. pylori.     

In silico prediction of human pathogenicity in the γ-proteobacteria

Background

Although the majority of bacteria are innocuous or even beneficial for their host, others are highly infectious pathogens that can cause widespread and deadly diseases. When investigating the relationships between bacteria and other living organisms, it is therefore essential to be able to separate pathogenic organisms from non-pathogenic ones. Using traditional experimental methods for this purpose can be very costly and time-consuming, and also uncertain since animal models are not always good predictors for pathogenicity in humans. Bioinformatics-based methods are therefore strongly needed to mine the fast growing number of genome sequences and assess in a rapid and reliable way the pathogenicity of novel bacteria.

Methodology/Principal Findings

We describe a new in silico method for the prediction of bacterial pathogenicity, based on the identification in microbial genomes of features that appear to correlate with virulence. The method does not rely on identifying genes known to be involved in pathogenicity (for instance virulence factors), but rather it inherently builds families of proteins that, irrespective of their function, are consistently present in only one of the two kinds of organisms, pathogens or non-pathogens. Whether a new bacterium carries proteins contained in these families determines its prediction as pathogenic or non-pathogenic. The application of the method on a set of known genomes correctly classified the virulence potential of 86% of the organisms tested. An additional validation on an independent test-set assigned correctly 22 out of 24 bacteria.

Conclusions

The proposed approach was demonstrated to go beyond the species bias imposed by evolutionary relatedness, and performs better than predictors based solely on taxonomy or sequence similarity. A set of protein families that differentiate pathogenic and non-pathogenic strains were identified, including families of yet uncharacterized proteins that are suggested to be involved in bacterial pathogenicity.     

Singlet Oxygen-Mediated Oxidation during UVA Radiation Alters the Dynamic of Genomic DNA Replication

UVA radiation (320–400 nm) is a major environmental agent that can exert its deleterious action on living organisms through absorption of the UVA photons by endogenous or exogenous photosensitizers. This leads to the production of reactive oxygen species (ROS), such as singlet oxygen (¹O₂) and hydrogen peroxide (H₂O₂), which in turn can modify reversibly or irreversibly biomolecules, such as lipids, proteins and nucleic acids. We have previously reported that UVA-induced ROS strongly inhibit DNA replication in a dose-dependent manner, but independently of the cell cycle checkpoints activation. Here, we report that the production of ¹O₂ by UVA radiation leads to a transient inhibition of replication fork velocity, a transient decrease in the dNTP pool, a quickly reversible GSH-dependent oxidation of the RRM1 subunit of ribonucleotide reductase and sustained inhibition of origin firing. The time of recovery post irradiation for each of these events can last from few minutes (reduction of oxidized RRM1) to several hours (replication fork velocity and origin firing). The quenching of ¹O₂ by sodium azide prevents the delay of DNA replication, the decrease in the dNTP pool and the oxidation of RRM1, while inhibition of Chk1 does not prevent the inhibition of origin firing. Although the molecular mechanism remains elusive, our data demonstrate that the dynamic of replication is altered by UVA photosensitization of vitamins via the production of singlet oxygen.     

Characterization of a Drosophila Ortholog of the Cdc7 Kinase: A ROLE FOR Cdc7 IN ENDOREPLICATION INDEPENDENT OF CHIFFON*

Cdc7 is a serine-threonine kinase that phosphorylates components of the pre-replication complex during DNA replication initiation. Cdc7 is highly conserved, and Cdc7 orthologs have been characterized in organisms ranging from yeast to humans. Cdc7 is activated specifically during late G₁/S phase by binding to its regulatory subunit, Dbf4. Drosophila melanogaster contains a Dbf4 ortholog, Chiffon, which is essential for chorion amplification in Drosophila egg chambers. However, no Drosophila ortholog of Cdc7 has yet been characterized. Here, we report the functional and biochemical characterization of a Drosophila ortholog of Cdc7. Co-expression of Drosophila Cdc7 and Chiffon is able to complement a growth defect in yeast containing a temperature-sensitive Cdc7 mutant. Cdc7 and Chiffon physically interact and can be co-purified from insect cells. Cdc7 phosphorylates the known Cdc7 substrates Mcm2 and histone H3 in vitro, and Cdc7 kinase activity is stimulated by Chiffon and inhibited by the Cdc7-specific inhibitor XL413. Drosophila egg chamber follicle cells deficient for Cdc7 have a defect in two types of DNA replication, endoreplication and chorion gene amplification. However, follicle cells deficient for Chiffon have a defect in chorion gene amplification but still undergo endocycling. Our results show that Cdc7 interacts with Chiffon to form a functional Dbf4-dependent kinase complex and that Cdc7 is necessary for DNA replication in Drosophila egg chamber follicle cells. Additionally, we show that Chiffon is a member of an expanding subset of DNA replication initiation factors that are not strictly required for endoreplication in Drosophila.