Abnormal open states patterns in the ATXN2 DNA sequence depends on the CAG repeats length.

Hereditary ataxias are one of the «anticipation diseases¼ types. Spinocerebral ataxia type 2 occurs when the number of CAG repeats in the coding region of the ATXN2 gene exceeds 34 or more. In healthy people, the CAG repeat region in the ATXN2 gene usually consists of 22-23 CAG trinucleotides. Mutations that increase the length of CAG repeats can cause severe neurodegenerative and neuromuscular disorders known as trinucleotide repeat expansion diseases. The mechanisms causing such diseases are associated with non-canonical configurations that can be formed in the CAG repeat region during replication, transcription or repair. This makes it relevant to study the zones of open states that arise in the region of CAG repeats under torque. The purpose of this work is to study, using mathematical modeling, zones of open states in the region of CAG repeats of the ATXN2 gene, caused by torque. It has been established that the torque effect on the 1st exon of the ATXN2 gene, in addition to the formation of open states in the promoter region, can lead to the formation of additional various sizes open states zones in the CAG repeats region. Moreover, the frequency of additional large zones genesis increases with increasing number of CAG repeats. The inverse of this frequency correlates with the dependence of the disease onset average age on the CAG repeats length. The obtained results will allow us to get closer to understanding the genetic mechanisms that cause trinucleotide repeat diseases.

Mathematical models describing oxygen binding by hemoglobin.

Despite the fact that the investigation of the structural and functional properties of hemoglobin dates back more than 150 years, the topic has not lost its relevance today. The most important component of these studies is the development of mathematical models that formalize and generalize the mechanisms determining the cooperative binding of ligands based on data on the structural and functional state of the protein. In this work, we review the mathematical relationships describing oxygen binding by hemoglobin, ranging from the classical Hüfner, Hill, and Adair equations to the Szabo-Karplus and tertiary two-state mathematical models based on the Monod-Wyman-Changeux and Koshland-Némethy-Filmer concepts. The generality of the considered equations as mathematical functions, bearing in their basis a power dependence, is demonstrated. The problems and possible solutions related to approximation of experimental data by the oxygenation equations with correlated fitting parameters are noted. Attention is paid to empirical equations, extended versions of the Hill equation, where the coefficient of cooperation is modulated by Gauss and Lorentz distributions as functions of partial oxygen pressure.

The Degradation and Aging of Biological Systems as a Process of Information Loss and Entropy Increase.

The problem of the degradation and aging of bioorganisms is herein considered from the viewpoint of statistical physics. Two typical timescales in biological systems-the time of metabolic processes and the time of the life cycle-are used. A kinetic equation describing the small timescales of the systems' characteristic processes in is proposed. Maintaining a biosystem in a time-stable state requires a constant inflow of negative entropy (negentropy). Ratios are proposed to evaluate the aging and degradation of systems in terms of entropy. As an example, the aging of the epithelium is studied. The connection of our approach to the information theory of aging is discussed, as well as theoretical constructions related to the concept of cooperon and its changing with time.

Theoretical Explanation for the Rarity of Antibody-Dependent Enhancement of Infection (ADE) in COVID-19.

Global vaccination against the SARS-CoV-2 virus has proved to be highly effective. However, the possibility of antibody-dependent enhancement of infection (ADE) upon vaccination remains underinvestigated. Here, we aimed to theoretically determine conditions for the occurrence of ADE in COVID-19. We developed a series of mathematical models of antibody response: model Ab-a model of antibody formation; model Cv-a model of infection spread in the body; and a complete model, which combines the two others. The models describe experimental data on SARS-CoV and SARS-CoV-2 infections in humans and cell cultures, including viral load dynamics, seroconversion times and antibody concentration kinetics. The modelling revealed that a significant proportion of macrophages can become infected only if they bind antibodies with high probability. Thus, a high probability of macrophage infection and a sufficient amount of pre-existing antibodies are necessary for the development of ADE in SARS-CoV-2 infection. However, from the point of view of the dynamics of pneumocyte infection, the two cases where the body has a high concentration of preexisting antibodies and a high probability of macrophage infection and where there is a low concentration of antibodies in the body and no macrophage infection are indistinguishable. This conclusion could explain the lack of confirmed ADE cases for COVID-19.

A New Model of Hemoglobin Oxygenation.

The study of hemoglobin oxygenation, starting from the classical works of Hill, has laid the foundation for molecular biophysics. The cooperative nature of oxygen binding to hemoglobin has been variously described in different models. In the Adair model, which better fits the experimental data, the constants of oxygen binding at various stages differ. However, the physical meaning of the parameters in this model remains unclear. In this work, we applied Hill's approach, extending its interpretation; we obtained a good agreement between the theory and the experiment. The equation in which the Hill coefficient is modulated by the Lorentz distribution for oxygen partial pressure approximates the experimental data better than not only the classical Hill equation, but also the Adair equation.

The Use of the Statistical Entropy in Some New Approaches for the Description of Biosystems.

Some problems of describing biological systems with the use of entropy as a measure of the complexity of these systems are considered. Entropy is studied both for the organism as a whole and for its parts down to the molecular level. Correlation of actions of various parts of the whole organism, intercellular interactions and control, as well as cooperativity on the microlevel lead to a more complex structure and lower statistical entropy. For a multicellular organism, entropy is much lower than entropy for the same mass of a colony of unicellular organisms. Cooperativity always reduces the entropy of the system; a simple example of ligand binding to a macromolecule carrying two reaction centers shows how entropy is consistent with the ambiguity of the result in the Bernoulli test scheme. Particular attention is paid to the qualitative and quantitative relationship between the entropy of the system and the cooperativity of ligand binding to macromolecules. A kinetic model of metabolism. corresponding to Schrödinger's concept of the maintenance biosystems by "negentropy feeding", is proposed. This model allows calculating the nonequilibrium local entropy and comparing it with the local equilibrium entropy inherent in non-living matter.

The Role of Acidosis in the Pathogenesis of Severe Forms of COVID-19.

COVID-19 has specific characteristics that distinguish this disease from many other infections. We suggest that the pathogenesis of severe forms of COVID-19 can be associated with acidosis. This review article discusses several mechanisms potentially linking the damaging effects of COVID-19 with acidosis and shows the existence of a vicious cycle between the development of hypoxia and acidosis in COVID-19 patients. At the early stages of the disease, inflammation, difficulty in gas exchange in the lungs and thrombosis collectively contribute to the onset of acidosis. In accordance with the Verigo-Bohr effect, a decrease in blood pH leads to a decrease in oxygen saturation, which contributes to the exacerbation of acidosis and results in a deterioration of the patient's condition. A decrease in pH can also cause conformational changes in the S-protein of the virus and thus lead to a decrease in the affinity and avidity of protective antibodies. Hypoxia and acidosis lead to dysregulation of the immune system and multidirectional pro- and anti-inflammatory reactions, resulting in the development of a "cytokine storm". In this review, we highlight the potential importance of supporting normal blood pH as an approach to COVID-19 therapy.

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA.

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC+ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.

A Method for Identification of the Methylation Level of CpG Islands From NGS Data.

In the course of sample preparation for Next Generation Sequencing (NGS), DNA is fragmented by various methods. Fragmentation shows a persistent bias with regard to the cleavage rates of various dinucleotides. With the exception of CpG dinucleotides the previously described biases were consistent with results of the DNA cleavage in solution. Here we computed cleavage rates of all dinucleotides including the methylated CpG and unmethylated CpG dinucleotides using data of the Whole Genome Sequencing datasets of the 1000 Genomes project. We found that the cleavage rate of CpG is significantly higher for the methylated CpG dinucleotides. Using this information, we developed a classifier for distinguishing cancer and healthy tissues based on their CpG islands statuses of the fragmentation. A simple Support Vector Machine classifier based on this algorithm shows an accuracy of 84%. The proposed method allows the detection of epigenetic markers purely based on mechanochemical DNA fragmentation, which can be detected by a simple analysis of the NGS sequencing data.

Structural features of DNA that determine RNA polymerase II core promoter.

BACKGROUND: The general structure and action of all eukaryotic and archaeal RNA polymerases machinery have an astonishing similarity despite the diversity of core promoter sequences in different species. The goal of our work is to find common characteristics of DNA region that define it as a promoter for the RNA polymerase II (Pol II). RESULTS: The profiles of a large number of physical and structural characteristics, averaged over representative sets of the Pol II minimal core promoters of the evolutionary divergent species from animals, plants and unicellular fungi were analysed. In addition to the characteristics defined at the base-pair steps, we, for the first time, use profiles of the ultrasonic cleavage and DNase I cleavage indexes, informative for internal properties of each complementary strand. CONCLUSIONS: DNA of the core promoters of metazoans and Schizosaccharomyces pombe has similar structural organization. Its mechanical and 3D structural characteristics have singular properties at the positions of TATA-box. The minor groove is broadened and conformational motion is decreased in that region. Special characteristics of conformational behavior are revealed in metazoans at the region, which connects the end of TATA-box and the transcription start site (TSS). The intensities of conformational motions in the complementary strands are periodically changed in opposite phases. They are noticeable, best of all, in mammals. Such conformational features are lacking in the core promoters of S. pombe. The profiles of Saccharomyces cerevisiae core promoters significantly differ: their singular region is shifted down thus pointing to the uniqueness of their structural organization. Obtained results may be useful in genetic engineering for artificial modulation of the promoter strength.

Optimization of signal-to-noise ratio for efficient microarray probe design.

MOTIVATION: Target-specific hybridization depends on oligo-probe characteristics that improve hybridization specificity and minimize genome-wide cross-hybridization. Interplay between specific hybridization and genome-wide cross-hybridization has been insufficiently studied, despite its crucial role in efficient probe design and in data analysis. RESULTS: In this study, we defined hybridization specificity as a ratio between oligo target-specific hybridization and oligo genome-wide cross-hybridization. A microarray database, derived from the Genomic Comparison Hybridization (GCH) experiment and performed using the Affymetrix platform, contains two different types of probes. The first type of oligo-probes does not have a specific target on the genome and their hybridization signals are derived from genome-wide cross-hybridization alone. The second type includes oligonucleotides that have a specific target on the genomic DNA and their signals are derived from specific and cross-hybridization components combined together in a total signal. A comparative analysis of hybridization specificity of oligo-probes, as well as their nucleotide sequences and thermodynamic features was performed on the database. The comparison has revealed that hybridization specificity was negatively affected by low stability of the fully-paired oligo-target duplex, stable probe self-folding, G-rich content, including GGG motifs, low sequence complexity and nucleotide composition symmetry. CONCLUSION: Filtering out the probes with defined 'negative' characteristics significantly increases specific hybridization and dramatically decreasing genome-wide cross-hybridization. Selected oligo-probes have two times higher hybridization specificity on average, compared to the probes that were filtered from the analysis by applying suggested cutoff thresholds to the described parameters. A new approach for efficient oligo-probe design is described in our study. CONTACT: shabalin@ncbi.nlm.nih.gov or olga.matveeva@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Non-random DNA fragmentation in next-generation sequencing.

Next Generation Sequencing (NGS) technology is based on cutting DNA into small fragments, and their massive parallel sequencing. The multiple overlapping segments termed "reads" are assembled into a contiguous sequence. To reduce sequencing errors, every genome region should be sequenced several dozen times. This sequencing approach is based on the assumption that genomic DNA breaks are random and sequence-independent. However, previously we showed that for the sonicated restriction DNA fragments the rates of double-stranded breaks depend on the nucleotide sequence. In this work we analyzed genomic reads from NGS data and discovered that fragmentation methods based on the action of the hydrodynamic forces on DNA, produce similar bias. Consideration of this non-random DNA fragmentation may allow one to unravel what factors and to what extent influence the non-uniform coverage of various genomic regions.

Sequence-specific ultrasonic cleavage of DNA.

We investigated the phenomenon of ultrasonic cleavage of DNA by analyzing a large set of cleavage patterns of DNA restriction fragments using polyacrylamide gel electrophoresis. The cleavage intensity of individual phosphodiester bonds was found to depend on the nucleotide sequence and the position of the bond with respect to the ends of the fragment. The relative intensities of cleavage of the central phosphodiester bond in 16 dinucleotides and 256 tetranucleotides were determined by multivariate statistical analysis. We observed a remarkable enhancement of the mean values of the relative intensities of cleavage (cleavage rates) in phosphodiester bonds following deoxycytidine, which diminished in the row of dinucleotides: d(CpG) > d(CpA) > d(CpT) >> d(CpC). The cleavage rates for all pairs of complementary dinucleotides were significantly different from each other. The effect of flanking nucleotides in tetranucleotides on cleavage rates of all 16 types of central dinucleotides was also statistically significant. The sequence-dependent ultrasonic cleavage rates of dinucleotides are consistent with reported data on the intensity of the conformational motion of their 5'-deoxyribose. As a measure of local conformational dynamics, cleavage rates may be useful for characterizing functional regions of the genome.

Optimization of duplex stability and terminal asymmetry for shRNA design.

Prediction of efficient oligonucleotides for RNA interference presents a serious challenge, especially for the development of genome-wide RNAi libraries which encounter difficulties and limitations due to ambiguities in the results and the requirement for significant computational resources. Here we present a fast and practical algorithm for shRNA design based on the thermodynamic parameters. In order to identify shRNA and siRNA features universally associated with high silencing efficiency, we analyzed structure-activity relationships in thousands of individual RNAi experiments from publicly available databases (ftp://ftp.ncbi.nlm.nih.gov/pub/shabalin/siRNA/si_shRNA_selector/). Using this statistical analysis, we found free energy ranges for the terminal duplex asymmetry and for fully paired duplex stability, such that shRNAs or siRNAs falling in both ranges have a high probability of being efficient. When combined, these two parameters yield a approximately 72% success rate on shRNAs from the siRecords database, with the target RNA levels reduced to below 20% of the control. Two other parameters correlate well with silencing efficiency: the stability of target RNA and the antisense strand secondary structure. Both parameters also correlate with the short RNA duplex stability; as a consequence, adding these parameters to our prediction scheme did not substantially improve classification accuracy. To test the validity of our predictions, we designed 83 shRNAs with optimal terminal asymmetry, and experimentally verified that small shifts in duplex stability strongly affected silencing efficiency. We showed that shRNAs with short fully paired stems could be successfully selected by optimizing only two parameters: terminal duplex asymmetry and duplex stability of the hypothetical cleavage product, which also relates to the specificity of mRNA target recognition. Our approach performs at the level of the best currently utilized algorithms that take into account prediction of the secondary structure of the target and antisense RNAs, but at significantly lower computational costs. Based on this study, we created the si-shRNA Selector program that predicts both highly efficient shRNAs and functional siRNAs (ftp://ftp.ncbi.nlm.nih.gov/pub/shabalin/siRNA/si_shRNA_selector/).

Distribution of unselectively bound ligands along DNA.

Unselective and reversible adsorption of ligands on DNA for a model of binding proposed by Zasedatelev, Gursky, and Volkenshtein is considered. In this model, the interaction between neighboring ligands located at the distance of i binding centers is characterized by the statistical weight ai. Each ligand covers L binding centers. For this model, expressions for binding averages are represented in a new simple form. This representation is convenient for the calculation of the fraction of inter-ligand distances of i binding centers fd(i) and the fraction of binding centers included in the distances of i binding centers fbc(i) for various types of interaction between bound ligands. It is shown that, for non-cooperative binding, contact cooperativity and long-range cooperativity, the fraction of the zero inter-ligand distance fd(0) is maximal at any relative concentration of bound ligands (r). Calculations demonstrate that, at low r, fd(0) approximately r.ao, and fd(i) approximately r at 11/r-L, then fd(i) rapidly decreases with i at any r for all types of inter-ligand interaction. At high ligand concentration (r is close to rmax=L(-1)), fd(0) is close to unity and fd(i) rapidly decreases with i for any type of inter-ligand interaction. For strong contact cooperativity, fd(0) is close to unity in a much lager r interval ((0.5-1).rmax), and fd(1) approximately ao(-1) at r approximately 0.5.rmax. In the case of long-range interaction between bound ligands, the dependence fd(i) is more complex and has a maximum at i approximately (1/r-L)1/2 for anti-cooperative binding. fbc(i) is maximal at i approximately 1/r-L for all types of binding except the contact cooperativity. A strong asymmetry in the influence of contact cooperativity and anticooperativity on the ligand distribution along DNA is demonstrated.

Nanoconstructions based on double-stranded nucleic acids.

We describe the formation and properties of nanoconstruction that consists of the double-stranded DNA molecules located at distance of 35-50 A in the spatial structure of particles of their cholesteric liquid-crystalline dispersions and cross-linked by artificial nanobridges. The resulting nanostructures possess the peculiar spatial and optical properties.