BACKGROUND AND OBJECTIVE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the causative agent of COVID-19 has significantly challenged the public health landscape in late 2019. After almost 3 years of the first ever SARS-CoV-2 case, the World Health Organization (WHO) declared the end of this global health emergency in May 2023. Although, despite the subsequent drop of COVID-19 cases, the SARS-CoV-2 infection still exhibited multiple waves of infection, primarily attributed to the appearance of new variants. Five of these variants have been classified as Variants of Concern (VOC): Alpha, Beta, Gamma, Delta, and the most recent, Omicron. Therefore, the development of methods for the timely and accurate detection of viral variants remains fundamental, ensuring an ongoing and effective response to the disease. This study aims to evaluate the feasibility of the application of an in-house approach in genomic surveillance for the detection of SARS-CoV-2 variants using in silico designed primers.

METHODS

The primers used for the study were particularly designed based on conserved regions of certain genes in the virus, targeting distinct mutations found in known variants of SARS-CoV-2. Viral RNA extracts from nasopharyngeal samples (n=14) were subjected to quantitative and qualitative tests (Nanodrop and AGE). Selected samples were then analyzed by RT-PCR and amplicons were submitted for sequencing. Sequence alignment analysis was carried out to identify the prevailing COVID-19 variant present in the sample population.

RESULTS

The study findings demonstrated that the in-house method was able to successfully amplify conserved sequences (spike, envelope, membrane, ORF1ab) and enabled identification of the circulating SARS-CoV-2 variant among the samples. Majority of the samples were identified as Omicron variant. Three out of four designed primers effectively bound into the conserved sequence of target genes present in the sample, revealing the specific SARSCoV-2 variant. The detected mutations characterized for Omicron found in the identified lineages included K417N, S477N, and P681H which were also identified as mutations of interest. Furthermore, identification of the B.1.448 lineage which was not classified in any known variant also provided the potential of the developed in-house method in detecting unknown variants of COVID-19.

CONCLUSION

Among the five VOCs, Omicron is the most prevalent and dominant variant. The in-house direct PCR product sequencing surveillance (DPPSS) method provided an alternative platform for SAR-CoV-2 variant analysis which is accessible and affordable than the conventional diagnostic surveillance methods and the whole genome sequencing. Further evaluation and improvements on the oligonucleotide primers may offer significant contribution to the development of a specific and direct PCRbased detection of new emerging COVID-19 variants.

Sars-cov-2 ; Polymerase Chain Reaction ; Dna Primers ; Oligonucleotide Primers ; Computer Simulation ; Conserved Sequence ; Coronavirus ; Covid-19 ; Disease ; Emergencies ; Evaluation Studies As Topic ; Genes ; Genome ; Global Health ; Health ; Identification (psychology) ; Infection ; Infections ; Membranes ; Methods ; Mutation ; Oligonucleotides ; Organizations ; Population ; Public Health ; Rna ; Rna, Viral ; Sars Virus ; Sequence Alignment ; Severe Acute Respiratory Syndrome ; Syndrome ; Viruses ; Whole Genome Sequencing ; World Health Organization

Computational modeling is an invaluable tool for mechanism analysis, directed engineering, and rational design of biological parts, metabolic networks, and even cellular systems. It can provide new technological solutions to address biological challenges at different levels and has become a central focus of research in biomanufacturing. In the computational modeling of proteins, which are the key parts in biological systems, the traditional physics-based methods (computer software and mathematical model) have been widely used to study the physical and chemical processes in the functioning of proteins, and have thus been recognized as a powerful tool for understanding complex biological systems and guiding experimental designs. As the scale of computational modeling continues to expand, traditional modeling techniques face difficulties in balancing computational accuracy and speed. In recent years, the explosive growth of biological data has made it possible to construct high-performance artificial intelligence (AI) models, which brings new opportunities to the computational modeling of proteins, and the AI-enhanced physics-based computational modeling technologies have emerged. This combined strategy not only incorporates the chemical knowledge and established physical principles but also is powerful in data processing and pattern recognition, which greatly improves the computational efficiency and prediction accuracy, as well as possesses stronger interpretation ability, transferability, and robustness. The AI-enhanced physics-based computational modeling technologies have already shown great potential and value in biocatalysis, paving a new way for the future development of biomanufacturing.
Artificial Intelligence ; Proteins/chemistry* ; Computer Simulation ; Software ; Computational Biology/methods*

Transcription factor (TF)-based biosensors have been widely applied in metabolic engineering, synthetic biology, metabolites monitoring, etc. These biosensors are praised for the high orthogonality, modularity, and operability. However, most natural TFs with weak responses and low specificity still demand optimization for desired performance in applications. Herein, we comprehensively summarize the recent advances in the engineering and optimization of TF-based biosensors with the assistance of computational simulation and artificial intelligence. This review includes the regulatory protein engineering aided by protein structure prediction and ligand binding simulation and the regulatory protein responses predicted by a mathematical model obtained from machine learning of mutagenesis data. In comparison with conventional tools, computational simulation and artificial intelligence enable more accurate and rapid design and construction of biosensors. Thus, these technologies will greatly promote the development of novel biosensors for applications.
Biosensing Techniques/methods* ; Transcription Factors/metabolism* ; Artificial Intelligence ; Protein Engineering/methods* ; Computer Simulation ; Synthetic Biology ; Machine Learning

As green, sustainable, and environmentally friendly material processing processes using biological cells or enzymes to achieve substance conversion, bioprocesses play an increasingly important role in biomanufacturing. It is difficult to optimize bioprocesses because of the complex relationship at multiple levels and multiple scales. The knowledge of mesoscale behaviors is the key to understanding the dynamics of bioprocesses and to sort out the complex relationships of parameter variations in the spatial-temporal domain. Mesoscale numerical simulation paves a way for understanding these phenomena, and the integration of artificial intelligence (AI) and mesoscale simulation offers new vitality into the optimization of bioprocesses. This article reviews the progress in mesoscale simulation and AI optimization of bioprocesses and discusses the possible development directions, aiming to promote the development of this field.
Artificial Intelligence ; Biotechnology/trends* ; Computer Simulation

Transcranial temporal interference stimulation (tTIS) is a novel non-invasive transcranial electrical stimulation technique that achieves deep brain stimulation through multiple electrodes applying electric fields of different frequencies. Current studies on the mechanism of tTIS effects are primarily based on rodents, but experimental outcomes are often significantly influenced by electrode configurations. To enhance the performance of tTIS within the limited cranial space of rodents, we proposed various electrode configurations for tTIS and conducted finite element simulations using a realistic mouse model. Results demonstrated that ventral-dorsal, four-channel bipolar, and two-channel configurations performed best in terms of focality, diffusion of activated brain regions, and scalp impact, respectively. Compared to traditional transcranial direct current stimulation (tDCS), these configurations improved by 94.83%, 50.59%, and 3 514.58% in the respective evaluation metrics. This study provides a reference for selecting electrode configurations in future tTIS research on rodents.
Animals ; Transcranial Direct Current Stimulation/instrumentation* ; Electrodes ; Mice ; Computer Simulation ; Finite Element Analysis ; Brain/physiology*

Clinical trials have demonstrated that kilohertz-frequency transcutaneous spinal cord stimulation (TSCS) can be used to facilitate the recovery of sensory-motor function for patients with spinal cord injury, whereas the neural mechanism of TSCS is still undetermined so that the choice of stimulation parameters is largely dependent on the clinical experience. In this paper, a finite element model of transcutaneous spinal cord stimulation was used to calculate the electric field distribution of human spinal cord segments T ₁₂ to L ₂, whereas the activation thresholds of spinal fibers were determined by using a double-cable neuron model. Then the variation of activation thresholds was obtained by varying the carrier waveform, the interphase delay, the modulating frequency, and the modulating pulse width. Compared with the sinusoidal carrier, the usage of square carrier could significantly reduce the activation threshold of dorsal root (DR) fibers. Moreover, the variation of activation thresholds was no more than 1 V due to the varied modulating frequency and decreases with the increased modulating pulse width. For a square carrier at 10 kHz modulated by rectangular pulse with the frequency of 50 Hz and the pulse width of 1 ms, the lowest activation thresholds of DR fibers and dorsal column fibers were 27.6 V and 55.8 V, respectively. An interphase delay of 5 μs was able to reduce the activation thresholds of the DR fibers to 20.1 V. The simulation results can lay a theoretical foundation on the selection of TSCS parameters in clinical trials.
Humans ; Spinal Cord Stimulation/methods* ; Nerve Fibers/physiology* ; Finite Element Analysis ; Spinal Cord/physiology* ; Computer Simulation ; Spinal Cord Injuries/physiopathology* ; Lumbosacral Region ; Lumbar Vertebrae ; Transcutaneous Electric Nerve Stimulation/methods* ; Models, Neurological

Alveolar bone reconstruction simulation is an effective means for quantifying orthodontics, but currently, it is not possible to directly obtain human alveolar bone material models for simulation. This study introduces a prediction method for the equivalent shear modulus of three-dimensional random porous materials, integrating the first-order Ogden hyperelastic model to construct a computed tomography (CT) based porous hyperelastic Ogden model (CT-PHO) for human alveolar bone. Model parameters are derived by combining results from micro-CT, nanoindentation experiments, and uniaxial compression tests. Compared to previous predictive models, the CT-PHO model shows a lower root mean square error (RMSE) under all bone density conditions. Simulation results using the CT-PHO model parameters in uniaxial compression experiments demonstrate more accurate prediction of the mechanical behavior of alveolar bone under compression. Further prediction and validation with different individual human alveolar bone samples yield accurate results, confirming the generality of the CT-PHO model. The study suggests that the CT-PHO model proposed in this paper can estimate the material properties of human alveolar bone and may eventually be used for bone reconstruction simulations to guide clinical treatment.
Humans ; Tomography, X-Ray Computed/methods* ; Porosity ; Alveolar Process/physiology* ; Bone Density ; Computer Simulation ; Elasticity ; X-Ray Microtomography ; Stress, Mechanical ; Finite Element Analysis ; Models, Biological

Tesla blood pumps demonstrate a reduced propensity for hemolysis and thrombosis compared with vane blood pumps. Considering the restricted driving force within the secondary flow channel of vane blood pumps, along with the low hydraulic efficiency of conventional Tesla blood pumps and their internal flow characteristics that significantly contribute to hemolysis and thrombosis, this study introduces a set of vanes atop the rotor of the Tesla blood pump. This forms a dual-fluid domain rotor, and an axial dual-outlet volute shell structure is adopted to realize the separation of the fluid domains. Through numerical simulations of the new structure, a comparative analysis was conducted in this study on the internal flow characteristics of double-outlet and single-outlet volute shells, and symmetric and asymmetric cross-sections of the same rotor. The results indicate that the flow field distribution is more uniform under the double-outlet volute shell structure, and overall energy dissipation is decreased. After implementing the double-outlet design, in the asymmetric cross-section, compared with the symmetric cross-section, the fluid velocity gradient and turbulent kinetic energy at the tongue of the septum are reduced, and the fluid velocity gradient at the convergence of the diffuser tube outlets are also decreased. The maximum scalar stress is lower, and the decline in head and efficiency is mitigated. Moreover, compared with the single-outlet volute shell, the hemolysis index in the asymmetric cross-section is reduced. In summary, this paper proposes a novel dual-outlet centrifugal disk blood pumps, which can provide a reference for the structural design and performance optimization of magnetically levitated centrifugal blood pumps.
Heart-Assist Devices ; Humans ; Equipment Design ; Hemolysis ; Computer Simulation

Glaucoma is the leading cause of irreversible blindness worldwide, with its primary risk factor arising from elevated intraocular pressure (IOP) due to an imbalance between aqueous humor production and outflow. This study aims to establish quantitative correlations among IOP, iris mechanical properties, channel microstructures, and aqueous humor dynamics through three-dimensional modeling and finite element analysis, overcoming the limitations of conventional experimental techniques in studying aqueous flow within the trabecular meshwork (TM) outflow pathway. A three-dimensional fluid-structure interaction (FSI) model incorporating the layered TM structure, Schlemm's canal (SC), iris, and other anterior segment tissues was developed based on human ocular anatomy. FSI simulations were performed to quantify the effects of IOP variations and iris Young's modulus on tissue morphology and aqueous humor dynamics parameters. The computational results demonstrated that axial iris deformation showed significant correlations with IOP and iris Young's modulus. Although elevated IOP exhibited minimal effects on hydrodynamic parameters in the anterior and posterior chambers, it markedly suppressed aqueous flow velocity in the TM region. Additionally, wall shear stress in SC and collector channels displayed high sensitivity to IOP variations. These findings reveal that the tissue mechanics-FSI mechanism modulates outflow resistance by regulating aqueous humor dynamics, offering valuable references for developing clinical therapies targeting IOP reduction in glaucoma management.
Humans ; Trabecular Meshwork/anatomy & histology* ; Finite Element Analysis ; Aqueous Humor/metabolism* ; Intraocular Pressure/physiology* ; Glaucoma/physiopathology* ; Iris/anatomy & histology* ; Computer Simulation ; Models, Biological

In order to solve the problems of difficult test, high cost and long cycle in the development of large-scale airborne negative pressure isolation system, the simulation analysis of negative pressure response characteristics is carried out around various aviation conditions such as aircraft ascending, leveling and descending, especially rapid decompression, based on the computational fluid dynamics (CFD) method. The results showed that the isolation cabin could achieve -50 Pa pressure difference environment and form a certain pressure gradient. The exhaust air volume reached the maximum value in the early stage of the aircraft's ascent, and gradually decreased with the increase of altitude until it was level flying. In the process of aircraft descent, the exhaust fan could theoretically maintain a pressure difference far below -50 Pa without working; Under the special condition of rapid pressure loss, it was difficult to deal with the rapid change of low pressure only by the exhaust fan, so it was necessary to design safety valve and other anti-leakage measures in the isolation cabin structure. Therefore, the initial stage of aircraft ascent is the key stage for the adjustment and control of the negative pressure isolation system. By controlling the exhaust air volume and adjusting parameters, it can adapt to the change of low pressure under normal flight conditions, form a relatively stable negative pressure environment, and meet the needs of biological control, isolation and transport.
Aircraft ; Computer Simulation ; Aviation/instrumentation* ; Humans ; Hydrodynamics ; Air Pressure ; Equipment Design ; Pressure