The Development of the SARS-CoV-2 Epidemic in Different Regions of Siberia in the 2020-2022 Period.

The comparison of the development of the SARS-CoV-2 epidemic in several neighboring regions can help researchers to assess the risks and develop more effective strategies and approaches in the field of preventive medicine. We analyzed the infection and mortality statistics for the 2020-2022 period in ten individual regions of the Siberian Federal District of Russia. We also sequenced complete genomes, which allowed us to analyze the genetic diversity of SARS-CoV-2 circulated in each of the ten regions and to build a phylogenetic dendrogram for the virus variants. The ParSeq v.1.0 software was developed to automate and speed up the processing and analysis of viral genomes. At the beginning of the pandemic, in the first two waves, the B.1.1 variant (20B) dominated in all regions of the Siberian Federal District. The third and fourth waves were caused by the Delta variant. Mortality during this period was at a maximum; the incidence was quite high, but the number of deposited genomes with GISAID during this period was extremely low. The maximum incidence was at the beginning of 2022, which corresponds to the arrival of the Omicron variant in the region. The BA.5.2 variant became the dominant one. In addition, by using NextClade, we identified three recombinants in the most densely populated areas.

From decision to action: Detailed modelling of frog tadpoles reveals neuronal mechanisms of decision-making and reproduces unpredictable swimming movements in response to sensory signals.

How does the brain process sensory stimuli, and decide whether to initiate locomotor behaviour? To investigate this question we develop two whole body computer models of a tadpole. The "Central Nervous System" (CNS) model uses evidence from whole-cell recording to define 2300 neurons in 12 classes to study how sensory signals from the skin initiate and stop swimming. In response to skin stimulation, it generates realistic sensory pathway spiking and shows how hindbrain sensory memory populations on each side can compete to initiate reticulospinal neuron firing and start swimming. The 3-D "Virtual Tadpole" (VT) biomechanical model with realistic muscle innervation, body flexion, body-water interaction, and movement is then used to evaluate if motor nerve outputs from the CNS model can produce swimming-like movements in a volume of "water". We find that the whole tadpole VT model generates reliable and realistic swimming. Combining these two models opens new perspectives for experiments.

Three-dimensional simulation of the Caenorhabditis elegans body and muscle cells in liquid and gel environments for behavioural analysis.

To better understand how a nervous system controls the movements of an organism, we have created a three-dimensional computational biomechanical model of the Caenorhabditis elegans body based on real anatomical structure. The body model is created with a particle system-based simulation engine known as Sibernetic, which implements the smoothed particle-hydrodynamics algorithm. The model includes an elastic body-wall cuticle subject to hydrostatic pressure. This cuticle is then driven by body-wall muscle cells that contract and relax, whose positions and shape are mapped from C. elegans anatomy, and determined from light microscopy and electron micrograph data. We show that by using different muscle activation patterns, this model is capable of producing C. elegans-like behaviours, including crawling and swimming locomotion in environments with different viscosities, while fitting multiple additional known biomechanical properties of the animal. This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

OpenWorm: overview and recent advances in integrative biological simulation of Caenorhabditis elegans.

The adoption of powerful software tools and computational methods from the software industry by the scientific research community has resulted in a renewed interest in integrative, large-scale biological simulations. These typically involve the development of computational platforms to combine diverse, process-specific models into a coherent whole. The OpenWorm Foundation is an independent research organization working towards an integrative simulation of the nematode Caenorhabditis elegans, with the aim of providing a powerful new tool to understand how the organism's behaviour arises from its fundamental biology. In this perspective, we give an overview of the history and philosophy of OpenWorm, descriptions of the constituent sub-projects and corresponding open-science management practices, and discuss current achievements of the project and future directions.This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

Hydrodynamic description of protein folding: the decrease of the probability fluxes as an indicator of transition states in two-state folders.

Using hydrodynamic description of protein folding, the process of the first-passage folding of ubiquitin has been studied. Since a large number of folding trajectories were required to obtain converged folding flows, a coarse-grained representation of the protein in the form of a C-bead Go-model was employed, and discrete molecular dynamics was used to perform simulations. It has been found that the free energy surface has a maximum width in the transition state region, so that the densities of folding flows (probability fluxes) decrease to minimum when the system passes through the transition state. There are indications that the increasing number of different protein conformations in the transition state region compared with those in the neighboring regions of semi-compact and native-like states is responsible for the present phenomena. It has also been shown that if the free energy is projected onto a single reaction coordinate, the low populations of the transition states can be compensated by the increasing number of states, which can lead to a considerable decrease or even disappearance of the free energy barrier in the transition state.

OpenWorm: an open-science approach to modeling Caenorhabditis elegans.

OpenWorm is an international collaboration with the aim of understanding how the behavior of Caenorhabditis elegans (C. elegans) emerges from its underlying physiological processes. The project has developed a modular simulation engine to create computational models of the worm. The modularity of the engine makes it possible to easily modify the model, incorporate new experimental data and test hypotheses. The modeling framework incorporates both biophysical neuronal simulations and a novel fluid-dynamics-based soft-tissue simulation for physical environment-body interactions. The project's open-science approach is aimed at overcoming the difficulties of integrative modeling within a traditional academic environment. In this article the rationale is presented for creating the OpenWorm collaboration, the tools and resources developed thus far are outlined and the unique challenges associated with the project are discussed.

Insights from the HuR-interacting transcriptome: ncRNAs, ubiquitin pathways, and patterns of secondary structure dependent RNA interactions.

The HuR protein regulates the expression of thousands of cellular transcripts by modulating mRNA splicing, trafficking, translation, and stability. Although it serves as a model of RNA-protein interactions, many features of HuR's interactions with RNAs remain unknown. In this report, we deployed the cryogenic RNA immunoprecipitation technique to analyze HuR-interacting RNAs with the Affymetrix all-exon microarray platform. We revealed several thousand novel HuR-interacting RNAs, including hundreds of non-coding RNAs such as natural antisense transcripts from stress responsive loci. To gain insight into the mechanisms of specificity and sensitivity of HuR's interaction with its target RNAs, we searched HuR-interacting RNAs for composite patterns of primary sequence and secondary structure. We provide evidence that secondary structures of 66-75 nucleotides enhance HuR's recognition of its specific RNA targets composed of short primary sequence patterns. We validated thousands of these RNAs by analysis of overlap with recently published findings, including HuR's interaction with RNAs in the pathways of RNA splicing and stability. Finally, we observed a striking enrichment for members of ubiquitin ligase pathways among the HuR-interacting mRNAs, suggesting a new role for HuR in the regulation of protein degradation to mirror its known function in protein translation.

Towards a virtual C. elegans: a framework for simulation and visualization of the neuromuscular system in a 3D physical environment.

The nematode C. elegans is the only animal with a known neuronal wiring diagram, or "connectome". During the last three decades, extensive studies of the C. elegans have provided wide-ranging data about it, but few systematic ways of integrating these data into a dynamic model have been put forward. Here we present a detailed demonstration of a virtual C. elegans aimed at integrating these data in the form of a 3D dynamic model operating in a simulated physical environment. Our current demonstration includes a realistic flexible worm body model, muscular system and a partially implemented ventral neural cord. Our virtual C. elegans demonstrates successful forward and backward locomotion when sending sinusoidal patterns of neuronal activity to groups of motor neurons. To account for the relatively slow propagation velocity and the attenuation of neuronal signals, we introduced "pseudo neurons" into our model to simulate simplified neuronal dynamics. The pseudo neurons also provide a good way of visualizing the nervous system's structure and activity dynamics.

Hydrodynamic description of protein folding.

A hydrodynamic description of protein folding is proposed and illustrated with a lattice protein model, which has a free energy surface (FES) typical of proteins with two-state folding kinetics. The flows from the unfolded to the native state are concentrated in a limited region of the FES. The rest is occupied by a flow "vortex", which does not lead to the native state. In contrast with intermediates that are associated with local minima, the vortex is not visible on the FES. The hydrodynamic interpretation thus provides new insights into the mechanism of protein folding and can be a useful complement to standard analyses.

AUG_hairpin: prediction of a downstream secondary structure influencing the recognition of a translation start site.

BACKGROUND: The translation start site plays an important role in the control of translation efficiency of eukaryotic mRNAs. The recognition of the start AUG codon by eukaryotic ribosomes is considered to depend on its nucleotide context. However, the fraction of eukaryotic mRNAs with the start codon in a suboptimal context is relatively large. It may be expected that mRNA should possess some features providing efficient translation, including the proper recognition of a translation start site. It has been experimentally shown that a downstream hairpin located in certain positions with respect to start codon can compensate in part for the suboptimal AUG context and also increases translation from non-AUG initiation codons. Prediction of such a compensatory hairpin may be useful in the evaluation of eukaryotic mRNA translation properties. RESULTS: We evaluated interdependency between the start codon context and mRNA secondary structure at the CDS beginning: it was found that a suboptimal start codon context significantly correlated with higher base pairing probabilities at positions 13 - 17 of CDS of human and mouse mRNAs. It is likely that the downstream hairpins are used to enhance translation of some mammalian mRNAs in vivo. Thus, we have developed a tool, AUG_hairpin, to predict local stem-loop structures located within the defined region at the beginning of mRNA coding part. The implemented algorithm is based on the available published experimental data on the CDS-located stem-loop structures influencing the recognition of upstream start codons. CONCLUSION: An occurrence of a potential secondary structure downstream of start AUG codon in a suboptimal context (or downstream of a potential non-AUG start codon) may provide researchers with a testable assumption on the presence of additional regulatory signal influencing mRNA translation initiation rate and the start codon choice. AUG_hairpin, which has a convenient Web-interface with adjustable parameters, will make such an evaluation easy and efficient.

A lattice protein with an amyloidogenic latent state: stability and folding kinetics.

We have designed a model lattice protein that has two stable folded states, the lower free energy native state and a latent state of somewhat higher energy. The two states have a sizable part of their structures in common (two "alpha-helices") and differ in the content of "alpha-helices" and "beta-strands" in the rest of their structures; i.e. for the native state, this part is alpha-helical, and for the latent state it is composed of beta-strands. Thus, the lattice protein free energy surface mimics that of amyloidogenic proteins that form well organized fibrils under appropriate conditions. A Go-like potential was used and the folding process was simulated with a Monte Carlo method. To gain insight into the equilibrium free energy surface and the folding kinetics, we have combined standard approaches (reduced free energy surfaces, contact maps, time-dependent populations of the characteristic states, and folding time distributions) with a new approach. The latter is based on a principal coordinate analysis of the entire set of contacts, which makes possible the introduction of unbiased reaction coordinates and the construction of a kinetic network for the folding process. The system is found to have four characteristic basins, namely a semicompact globule, an on-pathway intermediate (the bifurcation basin), and the native and latent states. The bifurcation basin is shallow and consists of the structure common to the native and latent states, with the rest disorganized. On the basis of the simulation results, a simple kinetic model describing the transitions between the characteristic states was developed, and the rate constants for the essential transitions were estimated. During the folding process the system dwells in the bifurcation basin for a relatively short time before it proceeds to the native or latent state. We suggest that such a bifurcation may occur generally for proteins in which native and latent states have a sizable part of their structures in common. Moreover, there is the possibility of introducing changes in the system (e.g., mutations), which guide the system toward the native or misfolded state.