This paper numerically investigates a 3D mesh independency study of a straight blade horizontal axis tidal turbine modelled using Computational Fluid Dynamics (CFD). The solution was produced by employing two turbulence models, the standard k-ε model and Shear Stress Transport (SST) in ANSYS CFX. Three parameters were investigated: mesh resolution, turbulence model, and power coefficient in the initial CFD, analysis. It was found that the mesh resolution and the turbulence model affect the power coefficient results. The power coefficients obtained from the standard k-ε model are 15% to 20% lower than the accuracy of the SST model. It can also be demonstrated that the torque coefficient increases with the increasing Tip Speed Ratio (TSR), but drops drastically after TSR = 5 and k-ε model failing to capture the non-linearity in the torque coefficient with the increasing TSR.

Starting from December 2019, coronavirus disease 2019 (COVID-19) has emerged as a once-in-a-century pandemic with deadly consequences, which urgently calls for new treatments, cures and supporting apparatuses. Remdesivir was reported by World Health Organization (WHO) as the most promising drug currently available for the treatment of COVID-19. Here, we use molecular dynamics simulations and free energy perturbation methods to study the inhibition mechanism of remdesivir to its target SARS-CoV-2 virus RNA-dependent RNA polymerase (RdRp). In the absence of a crystal structure of the SARS-CoV-2 RdRp, we first construct the homology model of this polymerase based on a previously available structure of SARS-CoV NSP12 RdRp (with a sequence identify of 95.8%). We then build the putative binding mode by aligning the remdesivir + RdRp complex to the ATP bound poliovirus RdRp. The putative binding structure is further optimized with molecular dynamics simulations and demonstrated to be stable, indicating a reasonable binding mode for remdesivir. The relative binding free energy of remdesivir is calculated to be -8.28 ± 0.65 kcal/mol, much stronger than the natural substrate ATP (-4.14 ± 0.89 kcal/mol) which is needed for the polymerization. The ~800-fold improvement in the K_d from remdesivir over ATP indicates an effective replacement of APT in blocking of the RdRp binding pocket. Key residues D618, S549 and R555 are found to be the contributors to the binding affinity of remdesivir. These findings demonstrate that remdesivir can potentially act as a SARS-CoV-2 RNA-chain terminator, effectively stopping its RNA reproduction, with key residues also identified for future lead optimization and/or drug resistance studies.

In this paper, we introduce a simple yet powerful and working version of the molecular dynamics code using the Python 3.9 language. The code contents are published in the link given in the appendix 1. The structure and components of the program is given in detail using flowcharts and code snippets. The program consists of major features like velocity verlet integrator, thermostats, COM removal, input and output modules, virial, pressure, and other thermodynamic quantities estimation etc. The author believes that this program will be helpful to graduate students who perform research in molecular dynamics simulations who intend to write their own code instead of the sophisticated open source packages.

Measurements of optical properties have been used for decades to study particle distributions in the ocean. They have been found useful to constrain suspended mass concentration as well as particle-related properties such as size, composition, packing (particle porosity or density) and settling velocity. Optical properties, however, provide measurements that are biased, as certain particles (based on size, composition, shape or packing) contribute to a specific property more than others. Here we study this issue both theoretically as well as by contrasting different optical properties collected simultaneously in a bottom boundary layer, to highlight the utility of such measurements as well as the biases we are likely to encounter using different optical properties to study suspended particles. In particular, we investigate the possibility to infer settling velocity from vertical profiles of optical measurements, finding that the effects of aggregation dynamics can seldom be ignored.

In this paper, the current technologies of the regenerative shock absorber systems have been categorized and evaluated. Three drive modes of the regenerative shock absorber systems, namely the direct drive mode, the indirect drive mode and hybrid drive mode are reviewed for their readiness to be implemented. The damping performances of the three different modes are listed and compared. Electrical circuit and control algorithms have also been evaluated to maximize the power output and to deliver the premium ride comfort and handling performance. Different types of parameterized road excitations have been applied to vehicle suspension systems to investigate the performance of the regenerative shock absorbers including that of the nonlinear regenerative shock absorber. The research gaps for comparison of the different drive modes and the nonlinearity analysis of the regenerative shock absorbers are identified and, the corresponding research questions have been proposed for future work.

Biomolecular structure drives function, and computational capabilities have progressed such that the prediction and computational design of biomolecular structures is increasingly feasible. Because computational biophysics attracts students from many different backgrounds and with different levels of resources, teaching the subject can be challenging. One strategy to teach diverse learners is with interactive multimedia material that promotes self-paced, active learning. We have created a hands-on education strategy with a set of fifteen modules that teach topics in biomolecular structure and design, from fundamentals of conformational sampling and energy evaluation to applications like protein docking, antibody design, and RNA structure prediction. Our modules are based on PyRosetta, a Python library that encapsulates all computational modules and methods in the Rosetta software package. The workshop-style modules are implemented as Jupyter Notebooks that can be executed in the Google Colaboratory, allowing learners access with just a web browser. The digital format of Jupyter Notebooks allows us to embed images, molecular visualization movies, and interactive coding exercises. This multimodal approach may better reach students from different disciplines and experience levels as well as attract more researchers from smaller labs and cognate backgrounds to leverage PyRosetta in their science and engineering research. All materials are freely available at https://github.com/RosettaCommons/PyRosetta.notebooks.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is pandemic. Prevention and control strategies require an improved understanding of SARS-CoV-2 dynamics. We did a rapid review of the literature on SARS-CoV-2 viral dynamics with a focus on infective dose. We sought comparisons of SARS-CoV-2 with other respiratory viruses including SARS-CoV-1 and MERS-CoV. We examined laboratory animal, and human studies. The literature on infective dose, transmission, and routes of exposure was limited specially in humans, and varying endpoints were used for measurement of infection. We propose the minimum infective dose of COVID-19 in humans, is higher than 100 particles, possibly slightly lower than the 700 particles estimated for H1N1 influenza. Despite variability in animal studies, there was some evidence that increased dose at exposure correlated with higher viral load clinically, and severer symptoms. Higher viral load measures did not reflect COVID-19 severity. Aerosol transmission seemed to raise the risk of more severe respiratory complications in animals. An accurate quantitative estimate of the infective dose of SARS-CoV-2 in humans is not currently feasible and needs further research. Further work is also required on the relationship between routes of transmission, infective dose, co-infection, and outcomes.

Alzheimer's disease (AD) is a progressive neurodegenerative brain disorder. One of the important therapeutic approaches of AD is the inhibition of β‐site APP cleaving enzyme‐1 (BACE1). This enzyme plays a central role in the synthesis of the pathogenic β-amyloid peptides (Aβ) in Alzheimer's disease. A group of potent BACE1 inhibitors with known x-ray structures (PDB ID 5i3X, 5i3Y, 5iE1, 5i3V, 5i3W, 4LC7, 3TPP) were studied by molecular dynamics simulation and binding energy calculation employing MM_GB(PB)SA. The calculated binding energies gave Kd values 0.139 µM, 1.39 nM, 4.39 mM, 24.3 nM, 1.39 mM, 29.13 mM and 193.07 nM, respectively. These inhibitors showed potent inhibitory activities in enzymatic and cell assays. The Kd values were compared with experimental values, the structures were discussed in view of the energy contributions to binding. Drug likeness of these inhibitors is also discussed. Accommodation of ligands in the catalytic site of BACE1 is discussed depending on the type of fragment involved in each structure. Molecular dynamics (MD) simulations and energy studies were used to explore the recognition of the selected BACE1 inhibitors by Asp 32, Asp228 and the hydrophobic flap. The results show that selective BACE1 inhibition may be due to the formation of strong electrostatic interactions with Asp32 and Asp228 and a large number of hydrogen bonds, π-π and Van der Waals interactions with the amino acid residues located inside the catalytic cavity. Interactions with the ligands show a similar binding mode with BACE1. These results help to rationalize the design of selective BACE1 inhibitors.

The global pandemic caused by the new SARS-COV-2 coronavirus makes it necessary to search for drugs for its control. Within of this research it has been known that the ivermectin drug, a FDA-approved drugs which is formulated as an 80:20 mixture of ivermectin B1a and B1b and used commonly for parasitic infections, has an inhibitory effect on viruses, includes SARS-COV-2 at in vitro level. In the particular case of SARS-COV-2 its mechanism of action remains elusive and controversial. Interestingly, the energy of interaction of ivermectin with any of the proteins the SARS-CoV-2 and the possible structural alterations at the protein level that this drug can cause have not been reported. In this sense, we carried out a bioinformatics study with docking strategies and molecular dynamics to predict the binding and disturbance induced by ivermectin in proteins associated with SARS-CoV-2. We use DockThor and Molegro docking scores. The DockThor server and myPresto software were used to build complexes and dynamics studies, respectively. The results obtained suggested that ivermectin is capable of docking with the 3CL protease and the HR2 domain, and may promote structural changes in these proteins by inducing unfolding/folding. Specifically, ivermectin brings protease to a significantly more deployed conformational state and the HR2 domain to a more compact state compared to the native state. Finally, it is shown that B1a and B1b macrocyclic lactones have a behavior different from to each target protein. These results suggest a possible inhibitory effect against SARS-CoV-2 due to a synergistic role of this drug to spontaneously bind to two important proteins involve in the proliferation of this virus. However, more studies are required on this possible mechanism of action.

A computational simulation of a PT6A gas turbine engine operating at off-design conditions is described in the present article. The model consists of a 0-dimensional thermo-fluidic description of the engine by applying the mass, linear momentum, angular momentum, and energy balances in each engine's component. The transient behavior of the engine is simulated for different blends of the original JET-A1 fuel with bio-diesel. Simulated thermodynamic variables of the air at each engine's component, as well as some performance parameters correspond to the reported experimental measurements. The numerical results also demonstrate the ability of the computational simulation to predict acceptable fuel blends, such that the efficiency of the engine is maximized and its structural integrity is maintained.

The SARS-CoV-2 was confirmed to cause the regional outbreak of coronavirus disease 2019 (COVID-19) in Wuhan, China. The 3C-like protease (3CL^pro), an essential enzyme for viral replication, is a valid target to compacts SARS-CoV and MERS-CoV. In this research, an integrated library consisting of 1000 compounds from Asinex Focused Covalent (AFCL) library and 16 FDA-approved protease inhibitors were screened against SARS-CoV-2 3CL^pro. Top compounds with significant docking scores and making stable interactions with catalytic dyad residues were obtained. The screening results in identification of compound 621 from AFCL library as well as Paritaprevir and Simeprevir from FDA-approved protease inhibitors as potential inhibitors of SARS-CoV-2 3CL^pro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CL^pro were characterized using 50 nanoseconds (ns) molecular dynamic (MD) simulation approach. The identified compounds are potential inhibitors worthy of further development as SARS-CoV-2 3CLpro inhibitors/drugs. Importantly, the identified FDA-approved therapeutics could be ready for clinical trials to treat infected patients and help to curb the COVID-19.

Copper demand for a country's copper industry has a greater pull effect. China's copper consumption in 2015 has accounted for 50% of the world. The scientific forecast of China's copper demands trend is also an important basis for analyzing its future environmental impact. This paper assumes that China's economy will be developing high, medium and low scenarios, and forecasts economic and social indicators such as total GDP, population and per capita GDP in China from 2016 to 2030. Then, predicted the demand of copper resources in China from 2016 to 2030 by the combination of system dynamics model, vector autoregressive moving average model and inverted U-type empirical model. The results show that: (1) in 2020, 2025 and 2030, China's refined copper demand will be 13 Mt, 15 Mt and 15.5 Mt. (2) China's copper demand growth slowed down significantly from 2016-2030. (3) 2025-2030, China's copper resource demand is stable, into the platform of demand growth, the highest peak value in 2027 will be 15.5 Mt. (4) 2030 years later, China's copper resource demand will enter a slow decline.

This article attempts to formulate a mathematical model for a potential explanation regarding the unavoidable impact of a rigid body's peculiar shape on the seamless flow over it. The solid body completely immersed in a Newtonian fluid and respectively has a relative open circuit flow on it will typically experience various observable phenomena like flow separation, flow transition, down-wash, stalling at the higher angle of attack, stalling velocity and how cambered airfoil can typically generate lift at a zero incidence angle. This article respectively represents an understanding of the laminar flow over a rigid body's external surface with due respect to its distinctive shape and size. This working paper formulates a more realistic and simplified mathematical model for open circuit laminar flow over a body, based on the historical data of aerodynamics and theoretical mechanics. This is intended to properly estimate forces on the continuous surface of the body in a laminar flow, to properly explain, understand and predict mentioned phenomena. Most of all the mechanism of streamline formation and its deformation with due regards to flow, shape and size of the body in an open-circuit laminar are formulated mathematically to enhance better design theory which can reduce experimentation while designing a streamlined body.

In this paper, the boiling flow inside a nanochannel with 700000 argon particle has been simulated by molecular dynamics (MD) simulation. This approach has been employed to analysis the superheated flow and its heat transfer pattern as well. For all simulations an external thrust force varying from 1 PN to 12 PN is exerted on inlet nanoparticles along the channel to have the forced annular boiling flow. Computations reveal that saturation condition and superheat degree have significant impacts on the liquid-vapor interface. Furthermore, because of the major influence of surface tension throughout a nanochannel, the x-velocity of liquid film and vapor core has not considerable fluctuations and stay smooth. All provided results show the behaviors completely similar to the available outcomes in the literature.

Caveolin-1 (Cav-1) is 22 kDa caveolae protein, acts as a scaffold within caveolar membranes. It interacts with alpha subunits of G-protein and thereby regulates their activity. Earlier studies reported elevated or up-regulated levels of caveolin-1 in the serum of prostate cancer patients. Secreted Cav-1 promotes angiogenesis, cell proliferation and anti-apoptotic activities in prostate cancer patients. Cav-1 upregulation is mainly related to prostate cancer metastasis. Keeping above facts in view, the present study was designed to explore Cav-1 as a target for prostate cancer therapy using computational approach. Molecular docking, structural base molecular modelling and molecular dynamics simulations were performed to investigate Cav-1 inhibitors. A predictive model was generated and validated to establish a stable structure. ZINC database of biogenic compounds was used for induced fit docking (IFD) and high throughput virtual screening. The H-bond interactions of the compounds with active site residues of Cav-1 were estimated by IFD and 100 ns long molecular dynamic simulations. The reported compounds showed significant binding and thus can be considered as potent therapeutic inhibitors of Cav-1. This study provides a valuable insight into biochemical interactions of Cav-1 for therapeutic applications and warrants for experimental validation of the predicted ‘active(s)’.

Distant Retrograde Orbits in the Earth-Moon system are gaining in popularity as stable `parking' orbits for various conceptual missions. To investigate the stability of potential Distant Retrograde Orbits, simulations were executed, with propagation running over a thirty-year period. Initial conditions for the vehicle state were limited such that the position and velocity vectors were in the Earth-Moon orbital plane, with the velocity oriented such that it would produce retrograde motion about Moon. The resulting trajectories were investigated for stability against the eccentric relative orbits of Earth and Moon in an environment that also included gravitational perturbations from Sun, Jupiter, and Venus, and the effects of radiation pressure. The results appear to indicate that stability is enhanced for certain resonant states within the Earth-Moon system.

This paper focuses on the problem of finite-time stabilization of homogeneous, non-Lipschitz systems with dilations. A key contribution of this paper is the design of a virtual recursive Hölder, non-Lipschitz state feedback which renders the non-Lipschitz systems in the special case dominated by a lower-triangular nonlinear system, finite-time stable. The proof is based on a recursive design algorithm developed recently, to construct the virtual Hölder continuous, finite-time stabilizer as well as a C¹ positive definite and proper Lyapunov function that guarantees finite-time stability of the non-Lipschitz non-linear systems.

After UV excitation, gas phase thymine returns to ground state in 5 to 7 ps, showing multiple time constants. There is no consensus on the assignment of these processes, with a dispute between models claiming that thymine is trapped either in the first (S1) or in the second (S2) excited states. In the present study, nonadiabatic dynamics simulation of thymine is performed on the basis of ADC(2) surfaces, to understand the role of dynamic electron correlation on the deactivation pathways. The results show that trapping in S2 is strongly reduced in comparison to previous simulations considering only non-dynamic electron correlation on CASSCF surfaces. The reason for the difference is traced back to the energetic cost for formation of a CO p bond in S2.

Assessment of the trends of land cover and vegetation dynamics (VD) using remote sensing (RS) and indicators such as anthropogenic activities and the socio-demographic information is essential in order to make proper planning for sustainable management. This paper attempts to evaluate land cover change (LCC) and VD in Kebbi State, Nigeria using historical Landsat data from 1986-2016 by means of remote sensing. The Driver-Pressure-State-Impact-Response (DPSIR) framework was later employed using both primary and secondary data for a better understanding of the drivers, the state of the environmental condition, the causes as well as the impact of the change. The images were classified into five thematic land cover classes as Dense Vegetation, shrubs/built area, farmland, bare/grassland and water body by means of Maximum likelihood supervised classification technique in accordance with Anderson classification scheme level 1, with acceptable accuracy. Pre-classification and post-classification change detection (CD) methodologies were executed using Normalized difference vegetation index (NDVI) and Image differencing respectively. The study illustrates a steady decline in dense vegetation and shrubs/build areas while farmland and bare/grassland increases, however, water bodies remain unchanged. The DPSIR pin-point that the major drivers of change in the study area have been the pressing need for farming land as the population grows and socioeconomic demands including fuelwood consumption and endemic poverty. Expansion of Farming land, fuelwood consumption and the need for construction materials are identified as the main key elements exerting pressure for the change. The state of the condition indicates a steady decline in dense vegetation and shrubs areas while farmland and bare/grassland are increasing significantly. The impacts include land degradation, the decline in the provision ecosystem goods and services, biodiversity loss through loss of habitats. The study, however, noted that many international and national policies in response to land degradation are channelled toward land restoration and remediating of the environment, through afforestation programs and improving the livelihood of the rural people through providing alternative income sources since they depend heavily on land for sustenance. However, the state governments, communities and individual commonly organized annual tree planting campaign with the main purpose of environmental protection.

Cancers are complex, adaptive ecosystems. They remain the leading cause of disease-related death among children in North America. As we approach computational oncology and Deep Learning Healthcare, our mathematical models of cancer dynamics must be revised. Recent findings support the perspective that cancer-microenvironment interactions consist of turbulent flows. As such, cancer pattern formation, protein-folding and metastatic invasion are discussed herein as processes driven by chemical turbulence within the framework of complex systems theory. Current state-of-the-art quantitative approaches used in reconstructing cancer stem cell networks are reviewed. To conclude, cancer stem cells are presented as strange attractors of the Waddington landscape.

This paper presents an overview of the principal structural and dynamics characteristics of reverse micelles (RMs) in order to highlight their structural flexibility and versatility, along with the possibility to modulate their parameters in a controlled-manner. The multifunctionality in a large range of different scientific fields is exemplified in two distinct directions: a theoretical model for mimicry of biological microenvironment and practical application in the field of nanotechnology and nano-based sensors. RMs represents a convenient experimental approach that limits the drawbacks of the conventionally biological studies in vitro, while the particular structure confers them the status of simplified mimics of cells by reproducing a complex supramolecular organization in an artificial system. The biological relevance of RMs is discussed in some particular cases referring to the confinement and crowding environment, molecular dynamics of water and cell membrane structure. The use of RMs in different range of applications seems to be more promising due to their structural and compositional flexibility, a high efficiency and selectivity being achieved. The advance in nanotechnology is based on developing new methods of nanomaterials synthesis and deposition. This review highlighting the advantages of using RMs in synthesis of nanoparticles with specific properties and in nano (bio)sensors design.

Plant life-long organogenesis involves sequential, time and tissue specific expression of developmental genes. This requires activities of Polycomb Group (PcG) and trithorax Group complexes, respectively responsible for repressive Histone 3 trimethylation at lysine 27 (H3K27me3) and activation-related H3K4me3. However, the genome-wide dynamics in histone modifications that occur during developmental processes have remained elusive. Here, we report the distributions of H3K27me3 and H3K4me3 along with transcriptional changes, in a developmental series including Arabidopsis leaf and three stages of flower development. We found that chromatin mark levels are highly dynamic over the time series on nearly half of all Arabidopsis genes. Moreover, during early flower morphogenesis, changes in H3K4me3 prime over changes in H3K27me3 and quantitatively correlate with transcription changes, while H3K27me3 changes occur after prolonged expression changes. Notably, early activation of PcG target genes is dominated by increases in H3K4me3 while H3K27me3 remains present at the locus. Our results reveal H3K4me3 as greater predictor over H3K27me3 for transcription dynamics, unveil unexpected chromatin mechanisms at gene activation and underline the relevance of tissue-specific temporal epigenomics.

Food systems will need to undergo considerable transformation. To be better prepared for and resilient to uncertainty and disturbances in the future, resource users and managers need to further develop knowledge about the food and farming system, with its dominating feedback structures and complexities, and to test robust and integrated system-based solutions. This paper investigates how participatory system dynamics modeling can be adapted to groups at the community level with low or no formal educational background. The paper also analyzes the refinement of workshop participants’ mental models as a consequence of a participatory system dynamics intervention. For this purpose, we ran two workshops with small-scale farmers in Zambia. Analysis of workshop data and post-workshop interviews shows that participatory system dynamics is well adaptable to support an audience-specific learning-by-doing approach. The use of pictures, objects and water glasses in combination with the basic aspects of causal loop diagramming makes for a well-balanced toolbox. Participants acquire understanding that is also relevant beyond systems thinking in that is offers a range of practical insights such as a critical evaluation of common food security strategies.

A tailored model for the assessment of environmental benefits achievable by “light-weighting” in the automotive field is presented. The model is based on the Fuel Reduction Value (FRV) coefficient, which expresses the Fuel Consumption (FC) saving involved by a 100 kg mass reduction. The work is composed of two main sections: simulation and environmental modelling. Simulation modelling performs an in-depth calculation of weight-induced FC whose outcome is the FRV evaluated for a wide range of Diesel Turbocharged (DT) vehicle case studies. Environmental modelling converts fuel saving to impact reduction basing on the FRVs obtained by simulations. Results show that for the considered case studies, FRV is within the range 0.115–0.143 and 0.142–0.388 L/100 km × 100 kg, respectively, for mass reduction only and powertrain adaptation (secondary effects). The implementation of FRVs within the environmental modelling represents the added value of the research and makes the model a valuable tool for application to real case studies of automotive lightweight LCA.

Waddington’s epigenetic landscape, a famous metaphor in developmental biology, depicts how a stem cell progresses from an undifferentiated phenotype to a differentiated one. The concept of “landscape” in the context of dynamical system theory represents a high-dimensional cell state space, in which each cell phenotype is considered as an “attractor” that is determined by interactions among multiple variables (molecular players), and is buffered against environmental fluctuations. Further, biological noise is thought to play an important role during these cell-fate decisions and in fact controls transitions between different phenotypes. Here, we discuss these phenotypic transitions in cancer from a dynamical systems perspective and invoke the concept of “cancer attractors” – hidden stable states of the underlying regulatory network that are not occupied by normal cells. Using epithelial-to-mesenchymal transition (EMT), cancer stem-like properties, metabolic reprogramming and the emergence of drug/hormone resistance as examples, we illustrate how phenotypic plasticity in cancer cells enables them to acquire hybrid phenotypes (such as hybrid epithelial/mesenchymal and hybrid metabolic phenotypes) that tend to be more aggressive and notoriously resilient to drug/hormone treatment. Furthermore, we highlight multiple factors that may give rise to phenotypic plasticity in cancer cells, such as (a) multi-stability or oscillatory behaviors governed by underlying regulatory networks involved in cell-fate decisions in cancer cells, and (b) network rewiring due to conformational dynamics of intrinsically disordered proteins (IDPs) that are highly enriched in cancer cells. We conclude by discussing why a therapeutic approach that promotes ‘recanalization’, i.e. the exit from “cancer attractors” and re-entry into “normal attractors”, is more likely to succeed rather than a conventional approach that targets individual molecules/pathways.

Terms related to gradients of scalar fields are introduced as scalar products into the formulation of entropy. A Lagrange density is then formulated by adding constraints based on known conservation laws. Applying the Lagrange formalism to the resulting Lagrange density leads to the Poisson equation of gravitation and also includes terms being related to curvature of space. The formalism further leads to terms possibly explaining nonlinear extensions known from modified Newtonian dynamics approaches. The article concludes with a short discussion of the presented methodology and provides an outlook on other phenomena, which might be tackled using this new approach.

We propose a Gauss model (GM), a map from time to the bell-shaped Gauss function to model the deaths per day and country, as a quick and simple model to make predictions on the coronavirus epidemic. Justified by the sigmoidal nature of a pandemic, i.e. initial exponential spread to eventual saturation, we apply the GM to existing data, as of April 2, 2020, from 25 countries during first corona pandemic wave and study the model's predictions. We find that logarithmic daily fatalities caused by Covid-19 are well described by a quadratic function in time. By fitting the data to second order polynomials from a statistical chi2-fit with 95\% confidence, we are able to obtain the characteristic parameters of the GM, i.e. a width, peak height and time of peak, for each country separately, with which we extrapolate to future times to make predictions. We provide evidence that this supposedly oversimplifying model might still have predictive power and use it to forecast the further course of the fatalities caused by Covid-19 per country, including peak number of deaths per day, date of peak, and duration within most deaths occur. While our main goal is to present the general idea of the simple modeling process using GMs, we also describe possible estimates for the number of required respiratory machines and the duration left until the number of infected will be significantly reduced.

Leigh Van Valen was an American evolutionary biologist who made major contributions to evolutionary theory and is particularly remembered by his groundbreaking paper "A New Evolutionary Law" (1973) where he provided evidence from fossil record data that this law maintains that the probability of extinction within any group remains es­sentially constant through time. In order to explain such unexpected result, Van Valen formulated a very influential idea that he dubbed the "Red Queen hypothesis". It states that the constant decay must be a consequence of evolutionary interactions among connected species within ecological networks. In Van Valen's picture, species do not merely evolve: they also coevolve with other species. As a consequence, when thinking in adaptation to an external environment, the other species must be considered as part (may be a major part) of such external world. Van Valen's law provided the first complex systems theory of coevolutionary dynamics and inspired a whole range of theoretical and experimental developments and scholars from very diverse fields, from economics to physics. In that respect, Leigh Van Valen's contribution percolated far beyond its original formulation. Red Queen arms races are nowadays considered a widespread feature of complex adaptive systems.

The computing power of smartphones has not received considerable attention in the mainstream education system. Most of the education-oriented smartphone applications (apps) are limited to general purpose services like Massive Open Online Courses (MOOCs), language learning, and calculators (performing basic mathematical calculations). Greater potential of smartphones lies in educators and researchers developing their customized apps for learners in highly specific domains. In line with this, we present Fluid Dynamics, a highly accurate Android application for measuring flow properties in compressible flows. This app can determine properties across the stationary normal and oblique shock, moving normal shock and across Prandtl $-$ Meyer expansion fan. This app can also measure isentropic flows, Fanno flows, and Rayleigh flows. The functionality of this app is also extended to calculate properties in the atmosphere by assuming the International Standard Atmosphere (ISA) relations and also flows across the Pitot tube. Such measurements are complicated and time-consuming since the relations are implicit and hence require the use of numerical methods, which give rise to repetitive calculations. The app is an efficient semi-implicit solver for gas dynamics formulations and uses underlying numerical methods for the computations in a graphical user interface (GUI), thereby easing and quickening the learning of concerned users. The app is designed for the Android operating system, the most ubiquitous and capable surveillance platform, and its calculations are based on JAVA based code methodology. In order to check its accuracy, the app's results are validated against the existing data given in the literature.

Coronavirus disease 2019 (COVID-19) is an emerging infectious disease that was first reported in Wuhan, China and has subsequently spread worldwide. In the absence of any antiviral or immunomodulatory therapies, the disease is spreading at an alarming rate. 5 to 10% of recovered patients in Wuhan test positive again; this suggest that for controlling COVID-19, vaccines may be better option than drugs. A clinical trial to evaluate an anti-COVID-19 vaccine has started recently. However, its efficacy and potency have to be evaluated and validated. As an alternative, we are presenting a first-of-its-kind, designed multi-peptide subunit based epitope vaccine against COVID-19. The vaccine construct comprise an adjuvant, CTL, HTL, and B-cell epitopes joined by linkers. The vaccine is non-toxic, non-allergenic, thermostable and immunogenic with the capability to elicit a humoral and cell-mediated immune response. The findings are validated with high-end computation-based methods. This unique vaccine is made up of 33 highly antigenic epitopes from three proteins that have a prominent role in host receptor recognition, viral entry, and pathogenicity. We advocate this vaccine must be synthesized and tested urgently as public health priority.

The intracellular calcium dynamics in vascular endothelial cells (VECs) in response to wall shear stress (WSS) and/or adenosine triphosphate (ATP) have been commonly regarded as an important factor in regulating VEC function and behavior including proliferation, migration and apoptosis. However, the effects of time-varying ATP signals have been usually neglected in the past investigations in the field of VEC mechanobiology. In order to investigate the combined effects of WSS and dynamic ATP signals on the intracellular calcium dynamic in VECs, a Y-shaped microfluidic device, which can provide the cultured cells on the bottom of its mixing micro-channel with stimuli of WSS signal alone and different combinations of WSS and ATP signals in one single micro-channel, is proposed. Both numerical simulation and experimental studies verify the feasibility of its application. Cellular experimental results also suggest that a combination of WSS and ATP signals rather than a WSS signal alone might play a more significant role in VEC Ca²⁺ signal transduction induced by blood flow.

The current outbreak of coronavirus disease (COVID-19) caused by SARS-CoV-2 in Wuhan, China has killed more than 2600 people since December 2019. Currently there is no effective treatment for this epidemic. Drug for anti SARS-CoV-2 are urgently needed. In this study we evaluated two compound libraries containing launched drugs and compounds from 300 kinds of Traditional Chinese Medicine in order to find anti SARS-CoV-2 drugs. Docking and then calculating binding free energy were performed as workflow against four key anti-SARS-CoV-2 drug targets, 3CLpro, PLpro and RdRp from SARSCoV-2, and AAK1 from human as well. As a result, drugs launched with potential for antiviral usage were selected in the hope of providing some knowledge for future drug discovery.

The emergence of SARS/MERS drug resistant COVID-19 with high transmission and mortality has recently been declared a deadly pandemic causing economic chaos and significant health problems. Like all coronaviruses, SARS-CoV-2 is a large virus that has many druggable components within its proteome. In this study, we focused on repurposing approved and investigational drugs by identifying potential drugs that are predicted to effectively inhibit critical enzymes within SARS-CoV-2. We shortlisted seven target proteins with enzymatic activities known to be essential at different stages of the virus life cycle. For virtual screening, the energy minimization of a crystal structure or modeled protein was carried out using Protein Preparation Wizard (Schrödinger LLC, 2020-1). Following active site selection based on data mining and COACH predictions, we performed a high-throughput virtual screen of drugs (n=5903) that are already approved by worldwide regulatory bodies including the FDA, using the ZINC database. Screening was performed against viral targets using three sequential docking modes (i.e. HTVS, SP and XP). Our in-silico virtual screening identified ~290 potential drugs based on the criteria of energy, docking parameters, ligand and binding site strain and score. Drugs specific to each target protein were further analyzed for binding free energy perturbation by molecular mechanics (prime MM-GBSA) and pruning the hits to the top 32 candidates. A top lead from each target group was further subjected to molecular dynamics simulation (MDS) using the Desmond module to validate the efficacy of the screening pipeline. All of the simulated hit-target complexes were predicted to strongly interact and with highly stable binding. Thus, we have identified a number of approved and investigational drugs with high likelihood of inhibiting a variety of key SARS-CoV-2 proteins. Follow-up studies will continue to identify inhibitors suitable for combination therapy based on drug-drug synergy to thwart resistance. In addition, the screening hits that we have identified provide excellent probes for understanding the binding properties of the active sites of all seven targets, further enabling us to derive consensus molecules through computer-aided drug design (CADD). While infections are expanding at a rampant pace, it must be recognized that resistance will grow commensurately through either genetic shift and/or genetic drift to all small molecule drugs identified. Vaccines should provide a more permanent solution through prevention, but resistivity is still a possible scenario. Nevertheless, a persistent multi-target drug development program is essential to curb this ongoing pandemic and to keep reemergence in check.

Arylamine N-acetyltransferases (NATs) are cytosolic enzymes, highly polymorphic, present in both eukaryotes and prokaryotes. These enzymes play an important role in the detoxification and activation of xenobiotics as well as in the synthesis of endogenous compounds. Specific NATs have been pointed out in the literature as possible therapeutic targets. In particular, the human NAT1, for the treatment of certain cancers, and the NAT from M. tuberculosis (TBNAT), for the treatment of tuberculosis. This paper describes an in silico approach to prospect and select potentially inhibitors of NAT1 and TBNAT from the Traditional Chinese Medicine (TCM) using free available tools. A library with ligands from TCM was previously screened in order to select only compounds with optimal pharmacological properties. The affinity of the selected ligands with respect to NAT enzymes was then evaluated by virtual screening (VS). Subsequently, the complexes with the best ligands were submitted to molecular dynamics (MD) simulations aiming to obtain better quality information on affinity and selectivity. The results for one specific ligand, ZINC14690579, indicated its potential for affinity and selectivity. ZINC14690579 structure may represent the discovery of a new scaffold for future development of NAT inhibitors.

Osteoarthritis (OA) is a common joint disorder found mostly in elderly people. The role of mechanical behavior in the progression of OA is complex and remains unclear. The stress-relaxation behavior of human articular cartilage in clinically defined osteoarthritic stages may have importance in diagnosis and prognosis of OA. In this study we investigated differences in the biomechanical responses among human cartilage of ICRS grades I, II and III using polymer dynamics theory. We collected 24 explants of human articular cartilage (eight each of ICRS grade I, II and III) and acquired stress-relaxation data applying a continuous load on the articular surface of each cartilage explant for 1180 s. We observed a significant decrease in Young’s modulus, stress-relaxation time, and stretching exponent in advanced stages of OA (ICRS grade III). The stretch exponential model indicated that significant loss in hyaluronic acid polymer might be the reason for the loss of proteoglycan in advanced OA. This work encourages further biomechanical modelling of osteoarthritic cartilage utilizing these data as input parameters to enhance the fidelity of computational models aimed at revealing how mechanical behaviors play a role in pathogenesis of OA.

Sustainable consumption became the community’s attention as the respond of worrying consumption growth direction that tends to excessively exploit nature without considering the continuity of the next generation. Take into account the growing coffee consumption in Indonesia, this article tries to prove whether the connoisseur consumers (CCs) are capable to mediate dynamics capabilities (DCs) of single origin coffee shops (SOCSs) and to encourage the consumption sustainability. In-depth interview to 30 SOCS managers, 60 baristas and 450 consumers, it is found that CCs are capable of supporting the business continuity of SOCSs because they are loyal consumers with the contribution of more than 20% from the total consumers. Correlation testing between connoisseurs’ attributes and the number of visits also shows positive value. It means there is a significant relationship. In addition, related to the sustainability attribute, the existence of CCs has encouraged the practice of consumption sustainability. At the same time, CCs are capable to moderate SOCSs to improve sensing, seizing, and transforming shop management to put up with the competition. Based on the above findings, CCs are moderating SOCSs in improving their DCs to be more potential in improving the sustainable consumption of coffee commodity in the future.

The concept of smart morphing blades, which can control themselves to reduce or eliminate the need for active control systems, is a highly attractive solution in blade technology. In this paper an innovative passive control system based on Shape Memory Alloys (SMAs) is proposed. On the basis of previous thermal and shape characterization of a single morphing blade for a heavy-duty automotive cooling axial fan, this study deals with the numerical analysis of the aerodynamic loads acting on the fan. By coupling CFD and FEM approaches it is possible to analyze the actual blade shape resulting from both the aerodynamic and centrifugal loads. The numerical results indicate that the polymeric blade structure ensures proper resistance and enables shape variation due to the action of the SMA strips.

After a local pollution control plant significantly reduced phosphorus loading into a phytoplankton-dominated tributary of the Potomac River in the early 1980’s, water quality and biological communities were monitored bi-weekly from April-September. After a 10-year time-lag, submerged aquatic vegetation (SAV), once abundant in this freshwater tidal embayment, returned to the area in 1993. After additional reductions in nitrogen load starting in 2000, the system completely switched to a SAV-dominated state in 2005. Fish abundance didn’t change during these distinct phase changes, but the fish community structure did. Increases in SAV provided refuge and additional spawning substrate for species with adhesive eggs such as Banded Killifish (Fundulus diaphanus), which is now the most abundant species in the embayment. Other changes seen were a decrease in the relative contribution of open water dwelling species such as White Perch (Morone americana), and an increase of visual predators such as Largemouth Bass (Micropterus salmoides). The 30-year record of data from this Potomac River tributary has revealed many important long-term trends that validate the effectiveness of initiatives to improve water treatment, and will aid in the continued management of the watershed and point source inputs.

First, this paper broaches the epistemological status of scientific tenets and approaches: phenomenological (descriptive only), well-founded (solid first principles, conducive to deep understanding), provisional (can be falsified if universal and verified if existential), and imaginary (fictitious entities or processes, conducive to empirically unsupported beliefs). The ΛCDM “concordance model” involves such beliefs: the emanation of the universe out of a non-physical stage, cosmic inflation (invented ad hoc), Λ (fictitious energy), and exotic dark matter. Big Bang cosmology further faces conceptual and pragmatic problems in delimiting what expands from what does not. The problems dissolve after untying inertia from space. The cosmology that emerges appears immediately compatible with the considered observations and the ‘perfect cosmological principle’. Waves and field perturbations that propagate at c expand exponentially with distance (a gravitational effect). The cosmic web of galaxies does not. Potential -Φ varies as H/(cz) instead of 1/r. Inertial forces arise from the gravitational action of the rest of the universe. Due to dilatation, they are reduced disproportionately at low accelerations. A cut-off value a₀ = 0.168 cH is deduced. This explains the successful description of galaxy rotation curves by MoND. A fully elaborated physical theory is still pending. Wider implications are briefly discussed.

An important stage in SARS-CoV-2 life cycle is the fusion of spike(S) protein with the ACE2 host-cell receptor. Therefore, to explore conserved features in S protein dynamics and to identify potentially novel regions for drugging, we measured variability derived from 791 viral genomes and studied its properties by MD simulation. The findings indicated that S2 subunit (HR1, CH, and CD domains) showed low variability, low fluctuations in MD, and displayed a trimer cavity. By contrast, the RBD domain, which is typically targeted in drug discovery programmes, exhibits more sequence variability and flexibility. Interpretations from MD suggest that the monomer is in constant motion showing transitions up-to-down state, and the trimer cavity may function as a 'bouncing spring' that may facilitates S protein interactions with ACE2. Feasibility of trimer cavity for potential drug target was examined by SBVS screening. Several hits that have already been validated or suggested to inhibit the SARS-CoV-2 virus in cell systems were identified; in particular, the data suggest an action mechanism for such molecules including Chitosan and macrolide types. These findings identify a novel binding-site formed by the S protein, that might assist in future drug discovery programmes aimed at targeting the CoV family of viruses.

As the interface between human and machine becomes blurred, hydrogel incorporated electronics and devices have emerged to be a new class of flexible/stretchable electronic and ionic devices due to their extraordinary properties, such as soft, mechanically robust and biocompatible. However, heat dissipation in these devices could be a critical issue and remains unexplored. Here, we report the experimental measurements and equilibrium molecular dynamics simulations of thermal conduction in polyacrylamide (PAAm) hydrogels. The thermal conductivity of PAAm hydrogels can be modulated by both the crosslinking density and water content in hydrogels. The crosslinking density dependent thermal conductivity in hydrogels varies from 0.33 to 0.51 Wm^-1K^-1, giving a 54% enhancement. We attribute the crosslinking effect to the competition between the increased conduction pathways and the enhanced phonon scattering effect. Moreover，water content can act as filler in polymers which lead to nearly 40% enhancement in thermal conductivity in PAAm hydrogels with water content vary from 23 to 88 wt%. Furthermore，we find the thermal conductivity of PAAm hydrogel is insensitive to temperature in the range of 25 ^oC – 40 ^oC. Our study offers fundamental understanding of thermal transport in soft materials and provides design guidance for hydrogel-based devices.

The asymmetry in the dynamics of an electro-pneumatic actuating device consisting of an electro-pneumatic transducer and a single-action pneumatic actuator was unexpectedly found experimentally. This asymmetry manifests in response to large step excitations. The dynamic asymmetry effect is understood as a change in the shape of the response of an actuator depending on the direction of the actuators stem movement. The questions appears: How to explain this effect? Does this phenomenon reflect thermodynamic air processes? Is it connected with air-to-mechanical energy conversion? Together, six working hypotheses explaining this effect were formulated. The asymmetry was studied in detail using analytical and simulation modeling, as well as experimental research. In this respect, a nonlinear analytical model was developed, tuned, and later solved using simulations. The simulation model was verified based on the experiment data. In addition, the problem of the efficiency in the energy conversion of a single-action actuator was discussed and, in result, the maximum theoretical energy efficiency was determined. Subsequently, all six working hypotheses were verified. Finally, the hypothesis explaining asymmetry as an effect of the different thermodynamic air processes in both actuator’s stem travel directions was confirmed.

This paper deals with the analytical modeling and control of rectilinear snake robots. During recent times snake robots have created much interest among researchers. The rectilinear pattern gait is one of the four biological snake locomotion modes. Rectilinear snakes have been widely used in rescue operations especially in rough terrains especially in narrow spaces where human intervention is not easy. Computational analysis of rectilinear motion is done using MATLAB.

Full-order state-space models represent the starting point for the development of advanced control methods for wind turbine systems (WTSs). Regarding existing control-oriented WTS models, two research gaps must be noted: (i) There exists no full-order WTS model in form of one overall ordinary differential equation that considers all dynamical effects which significantly influence the electrical power output; (ii) all existing reduced-order WTS models are subject to rather arbitrary simplifications and are not validated against a full-order model. Therefore, in this paper, two full-order nonlinear state-space models (of 11th and 9th-order in the (a, b, c)- and (d, q)-reference frame, resp.) for variable-speed variable-pitch permanent magnet synchronous generator WTSs are derived. The full-order models cover all relevant dynamical effects with significant impact on the system’s power output, including the switching behavior of the power electronic devices. Based on the full-order models, by a step-by-step model reduction procedure, two reduced-order WTS models are deduced: A non-switching (averaging) 7th-order WTS model and a non-switching 3rd-order WTS model. Comparative simulation results reveal that all models capture the dominant system dynamics properly. The full-order models allow for a detailed analysis covering the high frequency oscillations in the instantaneous power output due to the switching in the power converters. The reduced-order models provide a time-averaged instantaneous power output (which still correctly reflects the energy produced by the WTS) and come with a drastically reduced complexity making those models appropriate for large-scale power grid controller design.

In this paper, a two-phase model of air shock wave induced by rock-fall was described. The model was made up of the uniform motion phase (velocity was close to 0 m·s^-1) and the acceleration movement phase. The uniform motion phase was determined by experience, meanwhile the acceleration movement phase was derived by the theoretical analysis. A series of experiments were performed to verify the two-phase model and obtained the law of the uniform motion phase. The acceleration movement phase was taking a larger portion when height of rock-fall was higher with the observations. Experimental results of different falling heights showed good agreements with theoretical analysis values. Computational fluid dynamics (CFD) numerical simulation had been carried out to study the variation velocity with different falling height. As a result of this, the two-phase model could accurately and convenient estimating the velocity of air shock wave induced by rock-fall. The two-phase model could provide a reference and basis for estimating the air shock waves' velocity and designing the protective measures.

In this paper we develop from first principles a unique law pertaining to the flow of fluids through closed conduits. This law, which we call “Quinn’s Law”, may be described as follows: When fluids are forced to flow through closed conduits under the driving force of a pressure gradient, there is a linear relationship between the normalized dimensionless pressure gradient, P_Q, and the normalized dimensionless fluid current, C_Q. The relationship is expressed mathematically as: P_Q = k₁ +k₂C_Q. This linear relationship remains the same whether the conduit is filled with or devoid of solid obstacles. The law differentiates, however, between a packed and an empty conduit by virtue of the tortuosity of the fluid path, which is seamlessly accommodated within the normalization framework of the law itself. When movement of the fluid is very close to being at rest, i.e., very slow, this relationship has the unique minimum constant value of k₁, and as the fluid acceleration increases, it varies with a slope of k₂ as a function of normalized fluid current. Quinn’s Law is validated herein by applying it to the data from published classical studies of measured permeability in both packed and empty conduits, as well as to the data generated by home grown experiments performed in the author’s own laboratory.

In this paper, the nanoscale dissipative mechanisms of a Cu pad in a Ball Grid Array (BGA) packaging structure during isothermal ageing and uniaxial tension were investigated by the molecular dynamics (MD) method and experiments. From the result of the isothermal ageing test, a nonuniform consumption of Cu and large amount of Kirkendall voids were observed at the interface of Cu and Cu₃Sn. To study the effect of pressure and orientation on this phenomenon, MD simulations were conducted on four types of Cu-Cu₃Sn interface structures with different orientations of Cu. By comparing the diffusion coefficients of atoms in those cases, it was found that the tensile stress would inhibit the diffusion of atoms, whereas compressive stress would accelerate it, and this would be more significant under a larger magnitude of stress and temperature. Note that, in the model with the (101) surface Cu at the interface, both Cu and Cu₃Sn have a higher diffusion coefficient compared with the model with (001) surface Cu. Thus, the orientation of Cu will also contribute to the uniform consumption of the pad. Uniaxial tension simulation combined with DXA and CSP analyses on those models also shows the model with (001) surface Cu has a greater mechanical reliability in our simulations and related experiments.

Zinc finger proteins (ZFP) play important roles in cellular processes. The DNA binding region of ZFP consists of 3 zinc finger DNA binding domains connected by amino acid linkers, the sequence TGQKP connects ZF1 and ZF2, and TGEKP connects ZF2 with ZF3. Linkers act to tune the zinc finger protein in the right position to bind its DNA target, the type of amino acid residues and length of linkers reflect on ZF1-ZF2-ZF3 interactions and contribute to the search and recognition process of ZF protein to its DNA target. Linker mutations and the affinity of the resulting mutants to specific and nonspecific DNA targets were studied by MD simulations and MM_GB(PB)SA. The affinity of mutants to DNA varied with type and position of amino acid residue. Mutation of K in TGQKP resulted in loss in affinity due to the loss of positive K interaction with phosphates, mutation of G showed loss in affinity to DNA, WT protein and all linker mutants showed loss in affinity to a nonspecific DNA target, this finding confirms previous reports which interpreted this loss in affinity as due to ZF1 having an anchoring role, and ZF3 playing an explorer role in the binding mechanism. The change in ZFP-DNA affinity with linker mutations is discussed in view of protein structure and role of linker residues in binding.

The active power-assist function greatly expands the potential applications of exoskeleton robots, yet the motion intention estimation (MIE) for active power-assist strategy is quite problematic. Through the analysis of the conduction path and the different stage manifestations of motion intention in human body, we confirmed that the joint torque of human body meets the basic requirements of MIE for the active power-assist that we suggest, namely: (i) it reflects the direction and intensity of the wearer’s efforts; (ii) it precedes the human limb motion; (iii) it generates real-time and continuous output. Thus, an online calculation method of human joint torque was proposed. The sensing system integrated in exoskeleton robots was designed to perceive motion data and foot contact force of a human body. A special inverse dynamics with a parameterized model of the human body was proposed. Contrast experiments were carried out with the motion capture system, which results’ accuracy and similarity were evaluated via the root mean square error and correlation coefficient. The comparative analysis of two synchronous results shows good accuracy of the proposed MIE method, which lays the foundation for the realization of active power-assist.

Identifying or authenticating a computer user are necessary steps to keep systems secure on the network and to prevent fraudulent users from accessing accounts. Keystroke dynamics authentication can be used as an additional authentication method. Keystroke dynamics involves in-depth analysis of how you type on the keyboard, analysis of how long a key is pressed or the time between two consecutive keys. This field has seen a continuous growth in scientific research. In the last five years alone, about 10,000 scientific researches in this field have been published. One of the main problems facing researchers is the small number of public data sets that include how users type on the keyboard. This paper aims to provide researchers with a data set that includes how to type free text on the keyboard by 80 users. The data were collected in a single session via a web platform. The dataset contains 410,633 key-events collected in a total time interval of almost 24 hours. In similar research, most datasets are with texts written by users in English. The language in which the users wrote for this research is Romanian. This paper also provides an extensive analysis of the data set collected and presents relevant information for the analysis of the data set in future research.

The airborne inoculum of Fusarium circinatum, the fungal pathogen causing Pine Pitch Canker (PPC), is one of the main means of spread of the disease in forest stands and forest nurseries. Since this world-wide known pathogen was introduced in Europe, its biology in this newly infected area still remains scarcely known. To shed more light on this topic, we set an experiment on a naturally PPC infected forest of Monterey pine in Galicia (NW Spain) with the following two goals: (i) to describe the seasonal spore dispersal pattern during one year of regular sampling and (ii) to assess the spatial spore dispersal pattern around the infested plot. Portable rotating arm spore traps were used and complemented with meteorological measurements. The abundance of F. circinatum spores in the samples was evaluated by quantitative PCR (qPCR) with hydrolysis probe. The results showed almost permanent occurrence of the air inoculum throughout the whole year, being detected in 27 of the 30 samplings. No clear temporal trends were observed, but higher air inoculum was favoured by previous lower air temperatures and lower leaf wetness. Conversely, neither rainfall nor air humidity seemed to have any significant importance. The spatial spread of the inoculum was noted to be successful up to a distance of 1000 m in the wind direction, even with winds of just 5 m s^-1. Our study shows that rotating arm spore traps combined with qPCR may be an efficient tool for F. circinatum detection.

The novel coronavirus pneumonia epidemic is still raging around the world. Therefore, during the period of a global pandemic, maintaining the normal circulation of medicine and health fields is the key to fight against the epidemic. This paper based on relevant data such as the export volume of cross-border e-commerce in Hangzhou, China from 2015 to 2019, the article uses system dynamics modeling to explore the impact of overseas warehouse services on cross-border e-commerce export transactions. Studies have shown that the introduction of overseas warehouses has significantly improved customer satisfaction and improved the quality of cross-border logistics services, which has a positive impact on cross-border e-commerce export transactions. Through sensitivity analysis, it can be seen that under the same conditions, the efficiency of customs clearance has the greatest impact on overseas warehouse services. Therefore, when allocating various investments in overseas warehouses, you can favor the customs clearance process to enhance the competitiveness of the cross-border e-commerce market and provide a feasible experience for healthcare supplies export from cross-border e-commerce.

In the promptness of the COVID-19 outbreak, it would be very important to observe and estimate the pattern of diseases to reduce the contagious infection. To study this effect, we developed a COVID-19 analytical epidemic framework that combines with isolation and lockdown effect by incorporating five various groups of individuals. Then we analyze the model by evaluating the equilibrium points and analyzing their stability as well as determining the basic reproduction number. The extensive numerical simulations show the dynamics of a different group of the population over time. Thus, our findings based on the sensitivity analysis and the reproduction number highlight the role of outbreak of the virus that can be useful to avoid a massive collapse in Bangladesh and rest of the world. The outcome of this study concludes that outbreak will be in control which ensures the social and economic stability.

An evolutionary game is introduced which considers game-theoretic strategies in the context of non-linear population matrix models. This game considers the states and actions of the organisms of the evolving population, and a notion of dynamic equilibrium between strategies is described. The game’s formalism is expounded and a proof about equilibrium is given; specifically that any stable equilibrium can be described by proportions of pure strategies; particularly when population matrices are not defective.

Purification of biologically-derived therapeutics is a major cost contributor to the production of this rapidly growing class of pharmaceuticals. Monoclonal antibodies comprise a large percentage of these products therefore new antibody purification tools are needed. Small peptides, as opposed to traditional antibody affinity ligands such as Protein A, may have advantages in stability and production costs. Multiple heptapeptides that demonstrate Fc binding behavior that have been identified from a combinatorial peptide library using M13 Phage Display are presented herein. Seven unique peptide sequences of diverse hydrophobicity and charge were identified. All seven peptides showed strong binding to the four major human IgG isotypes, human IgM, as well as binding to canine, rat, and mouse IgG. These seven peptides were also shown to bind human IgG4 from DMEM cell culture media with 5% FCS and 5 g/L ovalbumin present. These peptides may be useful as surface ligands for antibody detection and purification purposes. Molecular docking and classical molecular dynamics (MD) simulations were conducted to elucidate the mechanisms and energetics for the binding of these peptides to the Fc region. The binding site was found to be located between the two glycan chains inside the Fc fragment. Both hydrogen bonding and hydrophobic interactions were found to be crucial for the binding interactions. Excellent agreement for the binding strength was obtained between experimental results and simulations.

Cancers remain the lead cause of disease-related, pediatric death in North America. The emerging field of complex systems has redefined cancer networks as a computational system with intractable algorithmic complexity. Herein, a tumor and its heterogeneous phenotypes are discussed as dynamical systems having multiple, strange attractors. Machine learning, network science and algorithmic information dynamics are discussed as current tools for cancer network reconstruction. Deep Learning architectures and computational fluid models are proposed for better forecasting gene expression patterns in cancer ecosystems. Cancer cell decision-making is investigated within the framework of complex systems and complexity theory.

Since crystals are made of periodic structures in space, predicting their three period vectors starting from any values based on the inside interactions is a basic theoretical physics problem. For the general situation where crystals are under constant external stress, we derived dynamical equations of the period vectors in the framework of Newtonian dynamics, for pair potentials recently (doi:/10.1139/cjp-2014-0518). The derived dynamical equations show that the period vectors are driven by the imbalance between the internal and external stresses. This presents a physical process where when the external stress changes, the crystal structure changes accordingly, since the original internal stress can not balance the external stress. The internal stress has both a full kinetic energy term and a full interaction term. It is extended to many-body interactions in this paper. As a result, all conclusions in the pair-potential case also apply for many-body potentials.

Given the high number of vehicle-crash victims, it has been established as a priority to reduce this figure in the transportation sector. For this reason, many of the recent researches are focused on including control systems in existing vehicles, to improve their stability, comfort and handling. These systems need to know in every moment the behavior of the vehicle (state variables), among others, when the different maneuvers are performed, to actuate by means of the systems in the vehicle (brakes, steering, suspension) and, in this way, to achieve a good behavior. The main problem arises from the lack of ability to directly capture several required dynamic vehicle variables, such as roll angle, from low-cost sensors. Previous studies demonstrate that low-cost sensors can provide data in real-time with the required precision and reliability. Even more, other research works indicate that neural networks are efficient mechanisms to estimate roll angle. Nevertheless, it is necessary to assess that the fusion of data coming from low-cost devices and estimations provided by neural networks can fulfill the reliability and appropriateness requirements for using these technologies to improve overall safety in production vehicles. Because of the increasing of computing power, the reduction of consumption and electric devices size, along with the high variety of communication technologies and networking protocols using Internet have yield to Internet of Things (IoT) development. In order to address this issue, this study has two main goals: 1) Determine the appropriateness and performance of neural networks embedded in low-cost sensors kits to estimate roll angle required to evaluate rollover risk situations. 2) Compare the low-cost control unit devices (Intel Edison and Raspberry Pi 3 Model B), to provide the roll angle estimation with this artificial neural network-based approach. To fulfil these objectives an experimental environment has been set up composed of a van with two set of low-cost kits, one including a Raspberry Pi 3 Model B, low cost Inertial Measurement Unit (BNO055 - 37€) and GPS (Mtk3339 - 53€) and the other having an Intel Edison System on Chip linked to a SparkFun 9 Degrees of Freedom module. This experimental environment will be tested in different maneuvers for comparison purposes. Neural networks embedded in low-cost sensor kits provide roll angle estimations very approximated to real values. Even more, Intel Edison and Raspberry Pi 3 Model B have enough computing capabilities to successfully run roll angle estimation based on neural networks to determine rollover risks situation fulfilling real-time operation restrictions stated for this problem.

Borneol and menthol were two terpenes wildly used as penetrate enhancer in transdermal drug delivery. To explore their penetration enhancement effect towards hydrophilic drug, 5-FU was selected as model drug. A method combined vitro permeation studies and coarse grain molecular dynamics was used to investigate their penetration enhancement effect towards 5-FU. As a result, although both borneol and menthol showed a penetration enhancement effect towards 5-FU, they differed a lot in the penetration enhancement mechanism, which was also thought to account for their different penetration enhancement effect. As for menthol, SC bilayer disrupting effect seemed to be its main mechanism. While for borneol, its mechanism seemed to be more complicated. Except for disrupting the SC bilayer, it could also increase the permeation of 5-FU by enhancing the diffusion rate of 5-FU or inducing the formation of transient pore. All of this enable us a molecular understanding of borneol and menthol’s penetration enhancement effect towards hydrophilic drug, which might provide some guidance in the latter research and application.

The current need to develop sustainable power sources has led to the development of ocean-based conversion systems. Wells turbine is a widely used converter in such systems which suffers from a lack of operational range and power production capacity under operational conditions. The profile named IFS which is concave in the post-mid-chord region, can produce significantly larger lift forces and show better separation behavior than the NACA profiles. In the present study, we tested this profile for the first time in a Wells turbine. The performance of six different blade designs with IFS and NACA profiles were evaluated and compared using a validated computational fluid dynamic model. Although the substitution of the NACA profile with the IFS profile in all cases increased the torque generated, the most efficient power generation and the largest efficient range were achieved in the design with varying thickness from the hub with a 0.15 thickness ratio reaching to the ratio of 0.2 at the tip. The operational span of this design with the IFS profile was 24.1% greater and the maximum torque generation was 71% higher than the case with the NACA profile. Therefore, the use of the IFS profile is suggested for further study and practical trials.

Complex blade responses such as a rotating stall or simultaneous resonances are common in modern engines and their observation can be a challenge even for state-of-the-art tip-timing systems and trained operators. This paper analyses forced vibrations of axial compressor blades, measured during the bench tests of the SO-3 turbojet. In relation to earlier studies conducted in ITWL with a small number of sensors, a multichannel tip-timing system let us observe simultaneous responses or higher-order modes. To find possible symptoms of a failure, blade responses in a healthy and unhealthy engine configuration with an inlet blocker were studied. The used analysis methods covered all-blade spectrum and the circumferential fitting of blade deflections to the harmonic oscillator model. The proposed modal solver can track the vibration frequency and adjust the engine order on the fly. That way, synchronous and asynchronous vibrations are observed and analysed together with an extended variant of least squares. The proposed approach helps to avoid common mistakes and saves a lot of work related to configuring the conventional solver.

Filament wound composite pipes are frequently used in the field were transmission of high pressured chemical fluids, disposal of industrial wastes, oil and natural gas transmission takes place. In oil and gas industry, the pipelines transporting heavy crude oil are subjected to variable pressure waves causing fluctuating stress levels in the pipes. Computational Fluid Dynamics Analysis was performed using Ansys 15.0 Fluent software to study the effects of these pressure waves on some specified joints in the pipes. Depending on the type of heavy crude oil being used, the flow behavior indicated a considerable degree of stress levels in certain connecting joints, causing the joints to become weak over a prolonged period of use. In this research comparison of various pipe joints was done by using different material and the output result of the stress levels of the pipe joints were checked so that the life of the pipe joints can be optimized by the change of material.

The coordination planning between land use and transportation system is an important premise of solving urban transportation problems and realizing land use integration. This study investigates the interactive and feedback relationship between land use and transportation system from the perspective of access management. By integrating the land use and traffic data from Las Vegas Metropolitan area with the system dynamics model, the causal relationship and causal loop diagrams (CLDs) are introduced to analyze the cause-and-effect relationship and quantitative relationship between the factors of the combined system of land use and transportation, and then sub-models partition and system simulation are performed. The systems dynamics model is established by analyzing the relationship between a series of access management techniques, traffic characteristics, and land use features. The results show that system dynamics model can be used as an effective alternative to model the symbiosis relationship of land use and transportation system for urban planning and construction.

Various distance metrics and their induced norms are employed in the quantitative modeling of evolutionary dynamics. Minimization of these distance metrics when applied to evolutionary optimization are hypothesized to result in different outcomes. Here, we apply the different distance metrics to the evolutionary trait dynamics brought about by the interaction between two competing species infected by parasites (exploiters). We present deterministic cases showing the distinctive selection outcomes under the Manhattan, Euclidean and Chebyshev norms. Specifically, we show how they differ in the time of convergence to the desired optima (e.g., no disease), and in the egalitarian sharing of carrying capacity between the competing species. However, when randomness is introduced to the population dynamics of parasites and to the trait dynamics of the competing species, the distinctive characteristics of the outcomes under the three norms become indistinguishable. Our results provide theoretical cases when evolutionary dynamics using different distance metrics exhibit similar outcomes.

Contemporary climate change is exposing plant populations to novel combinations of temperatures, drought stress, [CO₂] and other abiotic and biotic conditions. These changes are rapidly disrupting the evolutionary dynamics of plants. Despite the multifactorial nature of climate change, most studies typically manipulate only one climatic factor. In this opinion piece, we seek to explore how climate change factors interact with each other and with biotic pressures to alter evolutionary processes. We first explore the ramifications of climate change for key life history stages (germination, growth and reproduction). We then examine how mating system variation influences population persistence under rapid environmental change and propose that mixed mating could be advantageous in future climates. Furthermore, we discuss how spatial and temporal mismatches between plants and their mutualists and antagonists could promote or constrain adaptive responses to climate change. For example, plant-virus interactions vary from highly pathogenic to mildly facilitative, and are partly mediated by temperature, moisture availability and [CO₂]. Will host plants exposed to novel, stressful abiotic conditions be more susceptible to viral pathogens? Finally, we propose novel experimental approaches that could illuminate how plants will cope with unprecedented global change, such as resurrection studies combined with experimental evolution, genomics or epigenetics.

In this work, a tool for estimating the contact angle from the molecular dynamics simulations is developed and presented. The tool (Achilles) can detect water droplet on hydrophobic and hydrophilic surfaces. The tool can reconstruct the droplets broken across the periodic boundaries. Further a neighbor density based accurate filter is used to find the droplet liquid vapor interface and a circle is fitted using it after removing the dense layers of water next to solid surface. This fitted circle is solved for contact angle and results are outputted in the form of graphical images and text. The entire content of the internal computations of the tool is broken down into 4 phases and users can monitor the outcomes at every phase through output images. The tool is tested using sample molecular dynamics results of water droplet on hydrophobic and hydrophilic surfaces. We believe this tool can be a good addition to the molecular dynamics simulation community who work on the interfacial physics, droplet evaporation, super hydrophobic surfaces, and wettability etc.

Relationships among the population abundance of four predator groups for rice insect pests, namely: carabid beetles, staphylinid beetles, green mirid bugs, and spiders in three landscape categories were evaluated. Both rice plots and the associated bund margins of these rice plots found among three Bangladesh landscape categories were sampled by sweep net. The results revealed that the abundance significantly varied across landscapes. The rice landscape of one location harbored higher numbers of a specific predator than other location in other regions of Bangladesh. The results also showed a dependency on the width of the rice bund margins of the rice plots, where spiders populations increased with increased bund widths, but the population abundance of these predators did not depend on the diversity of the number of weed species found on the rice bund margins. The relative abundance of predator populations also significantly differed among the three landscapes, with the green mirid bug having the highest number among the four predators. This study indicates that predators of rice insect pests are highly landscape specific. In order to design integrated pest management systems for different Bangladeshi rice production locales, considerations unique to the characteristics of each locale are necessary. Preliminary efforts to apply variography analyses to the RED spectral band of LANDSAT 8 imagery from December 2016 are presented as first step toward learning a suite of methods which describe useful local characteristics affecting rice pest predators.

In this study, we present case studies to explore the dynamics of implementation and maintenance of health interventions. We analyze how specific interventions are built and eroded, how the building and erosion mechanisms are interconnected, and why we can see significantly different erosion rates across otherwise similar organizations. We use multiple comparative obesity prevention case studies to provide empirical information on the mechanisms of interest, and use qualitative systems modeling to integrate our evolving understanding into an internally consistent and transparent theory of the phenomenon. Our preliminary results identify reinforcing feedback mechanisms, including design of organizational processes, motivation of stakeholders, and communication among stakeholders, which influence implementation and maintenance of intervention components. Over time, these feedback mechanisms may drive a wedge between otherwise similar organizations, leading to distinct configurations of implementation and maintenance processes.

Due to poor ventilation conditions in the tunnel, if gas pipeline leaks, the consequence of the accident will be more serious. Therefore, before the emergency repair, gas in the tunnel needs to be excharged so as not to explode during the repair process. Therefore, it is necessary to study the ventilation of gas in the tunnel. Based on the computational fluid dynamics (CFD) theory and taking the Yanyingshan tunnel section of China-Myanmar pipeline as an example, this paper uses Fluent software to establish the leakage model of the gas pipeline and fan model in the tunnel and analyzes the influence of different fan locations and number of fans on gas concentration. It can be concluded that: (1) the use of press-in method makes it more efficient to discharge gas out of the tunnel. (2) In order to make ventilation efficient, the fan should be arranged in a higher position and needs to be at a distance from the top of the tunnel. (3) Parallel use of two fans has better ventilation effect than single fan.

This paper proposes the use of non-uniform extended surfaces installed externally to the tubes of the radiation section of fired heaters, in order to obtain a better heat flux distribution to the coils. To this end, the heat transfer mechanisms present in such equipment were studied through computational fluid dynamics (CFD), using simplified geometries that represent typical sizes of fired heaters. Also, a simplified model for the combustion was considered. Although this model oversimplifies the physics of the problem, it was able to give satisfactory results for the parameters being optimized, considering the main objective of this paper, that is to minimize the non-uniformity of heat flux in the tubes of the radiant section of fired heaters. It was possible to obtain optimized geometric parameters for different types of extended surfaces evaluated, coupling the results of these models with the Particle Swarm optimization method through the use of a response surface technique,. The results indicate a significant improvement in the uniformity of the heat flux distribution to the tubes through the use of the proposed extended surfaces. Thus, this solution reveals to be an interesting alternative to reduce the risks of fluid degradation and coking formation. Future studies must investigate the non-uniformity of the heat flux due to the presence of the flame and consider the interaction between the reactive flow and the participating medium. Nevertheless, this paper presents some results that justify the optimization of such extended surfaces taking into consideration thermal radiation.

Bimetallic nanoparticles are of interest since they lead to many interesting electrical, chemical, catalytic, and optical properties. They are particularly important in the field of catalysis since they show superior catalytic properties than their monometallic counterparts. The structures of bimetallic nanoparticles depend mainly on the synthesis conditions and the miscibility of the two components. In this work, PdPt alloyed-bimetallic nanoparticles (NPs) were synthesized through the polyol method, and characterized using spherical aberration (Cs) corrected scanning transmission electron microscopy (STEM). High-angle annular dark-field (HAADF)-STEM images of bimetallic nanoparticles were obtained. The contrast of images shows that nanoparticles have an alloy structure with an average size of 8.2 nm. Together with the characterization of nanoparticles, a systematic molecular dynamics simulations study, focused on the structural stability and atomic surface segregation trends in 923-atom PdPt alloyed-bimetallic NPs was carried out.

Amorphous solid dispersions are considered a promising formulation strategy for the oral delivery of poorly soluble drugs. The limiting factor for the applicability of this approach is the physical (in)stability of the amorphous phase in solid samples. Minimizing the risk of reduced shelf life for a new drug by establishing a suitable excipient/polymer-type from first principles would be desirable to accelerate formulation development. Here we perform Molecular Dynamics simulations to determine properties of blends of eight different polymer-small molecule drug combinations for which stability data is available from a consistent set of literature data. We calculate thermodynamic factors (mixing energies) as well as mobilities (diffusion rates and roto-vibrational fluctuations). We find that either of the two factors, mobility and energetics, can determine the relative stability of the amorphous form for a given drug. Which factor is rate limiting depends on physico-chemical properties of the drug and the excipients/polymers. The methods outlined here can be readily employed for an in-silico pre-screening of different excipients for a given drug to establish a qualitative ranking of the expected relative stabilities, thereby accelerating and streamlining formulation development.

Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators involve a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function.

Starting with a hypothetical geostrophic zonal current in an unbounded ocean, the investigation points out the response of this simple system to a thermal forcing, applied to the free surface and consistent with the maintenance of the geostrophic balance. The main result is the formation of a meridional component of the current, according to the Sverdrup relation, such that the full velocity vector rotates clockwise for heating and anticlockwise for cooling to adjust eventually in the initial zonal direction for large depths.

Dietary factors are one of the risk factors that can impact the brain chemistry, which leads to mental distress. Based on our data mining approach, we found that mental distress in men is associated with eating unhealthy food. Our aim in this paper is to apply results from our big data analytics approach to inform system dynamics (SD) modeling to investigate the causal relationships between brain structures, nutrients from food and dietary supplements, and mental health. We perform descriptive analysis based on a large data set to estimate the SD modeling parameters. Finally, we calibrate the model towards a time series data collected for individuals on their dietary and distress patterns. The results reveal that bridging these different methodologies leads to further insights from the SD model and decreases the error of calibrated parameter values. Future research is needed to validate our initial results for investigating the relationship between mental distress and dietary intake.

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can provide predictions, reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Molecular dynamics(MD) simulations are carried out to characterize the mechanicalproperties and failure behavior of van der Waals heterostructures composed ofgraphene and hexagonal boron nitride(hBN) single layer sheets. The graphene–hBNand hBN–graphene–hBN heterostructures simulations are carried out under tensileand shear deformation. Accordingly, stress versus strain curves of each systemare plotted and various properties of heterostructures, namely elastic modulusand shear modulus as well as failure stress and failure strain, are evaluatedand compared with one another as well as with the pristine graphene and hBNsheets. Subsequently, the failure mechanism/characterization of each sheet in theheterostructures is described. Alternatively, the elastic and shear modulus of eachheterostructure are computed by means of rule of mixture(ROM) which are in goodagreement with results that are obtained from MD simulations.

Molecular dynamics(MD) simulations are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals the solvent structure and preferential interaction sites (hot spots) on the protein surface. This information is largely transferable across all classes of protein ligands (from water to drugs) and can be used very effectively to understand ligand recognition and improve the predictive capability of well-established methods such as molecular docking. MD simulations for protein and drug or drug-like compounds are now being used but are still computationally expensive and can only be applied to specific cases. On the other hand, MDmix simulations can now be used in SBDD and we will describe the latest developments and implementations. We expect to see an increase in the application of these techniques to a plethora of protein targets to identify new drug candidates with the advent of new tools and faster computers.

A new form of coronavirus called severe acute respiratory disease coronavirus type 2 (SARS-CoV-2) is currently causing a pandemic. A six-month evolutionary history of SARS-CoV-2 is witnessed by characterising the total genome of 821 samples using comparative phylogenomic approaches. Our analyses produced striking inclusive results that may guide scientists/professionals for the past/future of pandemic. Phylogenetic and time estimation analyses suggest the proximate origin of pandemic strain as Guangdong and the origin time as first half of September 2019, not Wuhan and December 2019, respectively. The viral genome experienced a substitution rate similar to other RNA viruses, but it is particularly high in some of the peptides encoding sequences such as leader protein, E gene, orf8, orf10, nsp10, N gene, S gene and M gene and nsp4, while low in nsp11, orf7a, 3C-like proteinase, nsp9, nsp8 and endoRNase. Most strikingly, the divergence rate of amino acid sequences is high proportional to nucleotide divergence. Additionally, specific non-synonymous mutations in nsp3 and nsp6 evolved under positive selection. The exponential growth rate (r), doubling time (Td) and R₀ were estimated to be 47.43 per year, 5.39 days and 2.72, respectively. Comparison of synapomorphies distinguishing the SARS-CoV-2 and the candidate ancestor bat coronavirus indicates that mutation pattern in nsp3 and S gene enabled the new strain to invade human and become a pandemic strain. We arrive at the following main conclusions: (i) six months evolution of viral genome is nearly neutral, (ii) origin of pandemic is not Wuhan and predates formal reports, (iii) although viral population is ongoing an exponential growth, the doubling time is evolving towards shortening, and (iv) divergence rate of total genome is similar to other RNA viruses, but it is prominently high in some genes while low in some others and evolution in these genes should be closely monitored as their protein products intervening to pathogenicity, virulence and immune response.

In this paper an integrated methodology for the coupling between 1D- and 3D-CFD simulation codes is presented, which has been developed to support the design and calibration of new diesel engines. The aim of the proposed methodology is to couple 1D engine models, which may be available in the early-stage engine development phases, with 3D predictive combustion simulations, in order to obtain reliable estimates of engine performance and emissions for newly designed automotive diesel engines. The coupling procedure features simulations performed in 1D-CFD by means of GT-SUITE and in 3D-CFD by means of Converge, executed within a specifically designed calculation methodology. An assessment of the coupling procedure has been performed by comparing its results with experimental data acquired on an automotive Diesel engine, considering different working points including both part load and full load conditions. Different multiple injection schedules have been evaluated for part-load operation, including pre and post injections. The proposed methodology, featuring detailed 3D chemistry modeling, was proven to be capable to properly assess pollutant formation, specifically to estimate NOx concentrations. Soot formation trend was also well-matched for most of the explored working points. The proposed procedure can therefore be considered as a suitable methodology to support the design and calibration of new Diesel engines, thanks to its ability to provide reliable engine performance and emissions estimations from the early-stage of a new engine development.

Aquaporin-3 (AQP3) is one of the aquaglyceroporins, which is expressed in the basolateral layer of the skin membrane. Studies have reported that human skin squamous cell carcinoma overexpresses AQP3 and inhibition of its function may alleviate skin tumorigenesis. In the present study, we have applied a virtual screening method that encompasses filters for physicochemical properties and molecular docking to select potential hit compounds that bind to the Aquaporin-3 protein. Based on molecular docking results, the top 20 hit compounds were analyzed for stability in the binding pocket using unconstrained molecular dynamics simulations and further evaluated for binding free energy. Furthermore, examined the ligand-unbinding pathway of the inhibitor from its bound form to explore possible routes for inhibitor approach to the ligand-binding site. With a good docking score, stability in the binding pocket, and free energy of binding, these hit compounds can be developed as Aquaporin-3 inhibitors in the near future.

Nuclear export of viral mRNAs, is an essential step in the HIV replication cycle. This role is played by a small regulatory protein of HIV-1 called Rev.The N-terminal region of Rev contains an arginine-rich sequence. The arginine-rich motif (ARM) is located between amino acids 38-50 and forms an alpha-helical secondary structure. Expression of the structural proteins of human immunodeficiency virus type 1 requires the direct interaction of multiple copies of the viral Rev protein with its highly structured RNA target sequence, the Rev Response element (RRE). The major viral proteins are not produced if this transport of RNA is stopped. Therefore, knowledge of Rev structure is essential for understanding of its cooperative binding to the RRE, for understanding the mechanism of HIV infection and for the development of antiviral drugs that interfere with Rev’s essential functions and for acknowledgment of good candidate drugs for treatment of AIDS. To understand how REV interact with RRE element of HIV-RNA and its formation of oligomeric complex it is better to characterize the domain wise structure of REV with regard in function of each domain. Due to lack of structural data on Rev no single compound is reported as inhibitor of REV expect antiviral drugs. Identification of a high-affinity RNA-binding site for the human immunodeficiency virus type 1 Rev Protein is much more important. The ARM is a highly specific sequence which allows for the multimerization of Rev and also binding of REV with RNA. Here we are first time exploring the structural characteristics of REV protein both in free form and in complex with RNA at domain function level especially explore the role of ARM motif in REV HIV-1 protein as RNA binding sites by molecular dynamics (MD) simulation and homology modeling studies. Results indicate that the arginine-rich motif (ARM) is crucial in stability of this complex. The residues ARG38, 39, 41, 43, 44, 48, 50, and ASN40 are most interacting with nucleobases of RRE in Crystal structure of Rev and Rev-response-element RNA complex. Our study plays a major role in elaboration of binding of RNA with REV and pave the way for further investigation for therapeutically agent for HIV.

Many species of plants are found in regions to which they are alien. Their global distributions are characterised by a family of exponential functions of the kind that arise in elementary statistical mechanics (an example in ecology is MacArthur's broken stick). We show here that all these functions are quantitatively reproduced by a model containing a single parameter – some global resource partitioned at random on the two axes of species number and site number. A dynamical model generating this equilibrium is a two fold stochastic process and suggests a curious and interesting biological interpretation in terms of niche structures fluctuating with time and productivity; with sites and species highly idiosyncratic. Idiosyncrasy implies that attempts to identify a priori those species likely to become naturalized are unlikely to be successful. Although this paper is primarily concerned with a particular problem in population biology, the two fold stochastic process may be of more general interest.

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. There is evidence that some ligands that bind to the GPCRs activate different downstream signaling pathways. G protein activation or -arrestin biased signaling involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. Molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions between receptors and ligands. Here, we highlight recent simulation studies and methods used to study biased G protein-coupled receptor signaling and their conformational dynamics. We also highlight applications of MD simulations to drug discovery.

Background: Since March 2020, Spain is severely hit by the ongoing pandemic of coronavirus disease 19 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Understanding and disrupting the early transmission dynamics of the infection is crucial for impeding sustained transmission. Methods: We recorded all COVID-19 cases and traced their contacts in an isolated rural community. We also sampled 10 households, 6 public service sites and the wastewater from the village sewage for environmental SARS-CoV-2 RNA. Results: The first village patient diagnosed with COVID-19-compatible symptoms occurred on March 3, 2020, twelve days before lockdown. A peak of 39 cases occurred on March 30. By May 15, the accumulated number of symptomatic cases was 53 (6% of the population), of which only 22 (41%) had been tested and confirmed by RT-PCR as SARS-CoV-2 infected, including 16 hospitalized patients. Contacts (n=144) were six times more likely to develop symptoms. Environmental sampling detected SARS-CoV-2 RNA in two households with known active cases and in two public service sites: the petrol station and the pharmacy. Samples from other sites and the wastewater tested negative. Conclusions: The low proportion of patients tested by RT-PCR calls for urgent changes in disease management. We propose that early testing of all cases and their close contacts would reduce infection spread, reducing the disease burden and fatalities. In a context of restricted testing, environmental RNA surveillance might prove useful for early warning and to identify high-risk settings enabling a targeted resource deployment.

This study developed a system dynamics model to understand the spread of coronavirus disease (COVID-19) in Japan. The model is built on the generic SEIR (Susceptible, Exposed, Infected, and Recovered) framework adopted by Ghaffarzadegan and Rahmandad (2020) to build a system dynamics model for the spread of COVID-19 in Iran. Japan seems to be successful in containing the spread compared with other countries, and its first peak has passed. However, because this also leaves a large population still susceptible to the virus, it could cause a second and potentially higher peak of infection after the state of emergency aimed at reducing contact rate is lifted. As the government has proposed the “New Lifestyle,” it is critical to behave cautiously so as not to be infected. While the model focuses on the SEIR structure, the reflection of other sub-structures such as economy and hospital capacity that have tradeoffs with reducing contact rate should be beneficial. However, before extending the model, it is also critical to conduct estimates using confidence intervals rather than point estimates to better reflect uncertainties.

Three cases of Day-Night Band (DNB) observations of the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) are explored for applications to assess the ocean environment and monitor ocean dynamics. An approach to use the DNB radiance ratio was developed in order to better continuously assess the ocean diurnal and short-term environmental changes with VIIRS DNB observations. In the La Plata River Estuary, the sediment fronts showed 20–25 km diurnal inshore-offshore movements on March 13, 2017. In the waters off the Argentina coast in the South Atlantic, VIIRS DNB measurements provided continuous observations and monitoring of the algae bloom development and migration between 24–26 March 2016. This algae bloom generally kept the same spatial patterns, but moved nearly 20 km eastward in the three-day period. In the Yangtze River Estuary and Hangzhou Bay region in China’s east coast, VIIRS DNB observations also revealed the complicated coastal dynamic changes between 12–14 April 2017. Even though there are still some challenges and limitations for monitoring the ocean environment with VIIRS DNB observations, this study shows that satellite DNB observations can provide additional data sources for ocean observations, especially observations during the nighttime.

This study aims to determine the optimal configuration (position and operation duration) for wall mounted mechanical mixers based on the comparison of three-dimensional computational fluid dynamics (CFD) modelling results and physical data collected from the treatment plant. A three-dimensional model of anoxic zone 1, 2 and 3 of Northern Wastewater Treatment Plant (WWTP) located at Cairns Regional Council, Cairns, Queensland, Australia was developed and validated. The model was used to simulate the flow pattern of the WWTP and the simulation results are in good agreement with the physical data varying between 0% to 15% in key locations. The anoxic zones were subject to velocities less than the desired 0.3 metres per second however results for suspended solids concentration indicate that good mixing is being achieved. Results for suspended solids concentrations suggest that the anoxic zones are towards the upper limits recommended by literature for specific power dissipation. The duration for operation of mechanical mixers was investigated and identified that the duration could be reduced from 900 seconds down to 150 seconds. Alternative mixer positioning was also investigated and identified positioning which would increase the average flow velocity with decreased duration (150 seconds). The study identified that Council may achieve savings of $24,000 per year through optimisation of the mechanical mixers.

We consider a random dynamical system arising in a model of associative memory. This system can be seen as a small (stochastic and deterministic) perturbation of a determinstic system having two weak attractors which are destroyed after the perturbation. We show, with a computer aided proof, that the system has a kind of chaotic itineracy. Typical orbits are globally chaotic, while they spend relatively long time visiting attractor's ruins.

To build a representation of what we see, the human brain recruits regions throughout the visual cortex in cascading sequence. Recently, an approach was proposed to evaluate the dynamics of visual perception in high spatiotemporal resolution at the scale of the whole brain. This method combined functional magnetic resonance imaging (fMRI) data with magnetoencephalography (MEG) data using representational similarity analysis and revealed a hierarchical progression from primary visual cortex through the dorsal and ventral streams. To assess the replicability of this method, here we present results of a visual recognition neuro-imaging fusion experiment, and compare them within and across experimental settings. We evaluated the reliability of this method by assessing the consistency of the results under similar test conditions, showing high agreement within participants. We then generalized these results to a separate group of individuals and visual input by comparing them to the fMRI-MEG fusion data of Cichy et al. (2016), revealing a highly similar temporal progression recruiting both the dorsal and ventral streams. Together these results are a testament to the reproducibility of the fMRI-MEG fusion approach and allows for the interpretation of these spatiotemporal dynamic in a broader context.

The article presents the stages of modeling and simulation of a new design of a wheelchair with the option of moving up and down stairs. These analyzes were aimed at the synthesis of the de-sign parameters and parameters of the sensor and control systems. The simulation results were verified by experimentally testing the prototype.

Among the best-known capillarity phenomena is a capillary rise, the understanding of which is essential in fluidics. Some capillary flows rise monotonically whereas others oscillate, but until now no criteria have been formulated for this scenario. In this paper, the Levine's capillary rise modelling is computed numerically, then the critical radius of the capillary tube is formulated by using the dimensional method and data fitting for identification of exponent index. The phase space diagram of capillary velocity versus height is obtained for the first time and shows that the phase transition from oscillating to monotonic rising happens when the phase trajectory decreases exponentially to somewhere other than the "attractor." Two general Maple codes of the problem are provided as an essential part of this paper.

Histone Lysine Specific Demethylase 1 (LSD1) is overexpressed in many cancers and become a new target for anticancer drugs. In recent years, the small molecule inhibitors with various structures targeting LSD1 have been reported. Here we report the binding interaction modes of a series of thieno[3,2-b]pyrrole-5-carboxamides LSD1 inhibitors using molecular docking, three dimensional quantitative structure-activity relationship (3D-QSAR). Comparative molecular field analysis (CoMFA q²=0.783, r²=0.944, r²_pred=0.851) and Comparative molecular similarity indices analysis (CoMSIA q²=0.728, r²=0.982, r²_pred=0.814) were used to establish 3D-QSAR models, which had good verification and prediction capabilities. Based on the contour maps and the information of molecular docking, 8 novel small molecules were designed in silico, among which compounds D4, D5 and D8 with high predictive activity were subjected to further molecular dynamics simulations (MD), and their possible binding modes were explored. It was found that Asn535 plays a crucial role in stabilizing the inhibitors. Furthermore, the ADME and bioavailability prediction for D4, D5 and D8 were carried out. The results would provide valuable guidance for designing new reversible LSD1 inhibitors in the future.

Flavonoids are widely recognized as natural polydrugs, given their anti-inflammatory, antioxidant, sedative and antineoplastic activities. Recently, different studies have shown that flavonoid have the potential to inhibit BET bromodomains. Previous reports suggest that flavonoids are putative inhibitors of the ZA channel due to their orientation and interactions with P86, V87, L92, L94 and N140. Herein, a comprehensive characterization of the binding mode of the biflavonoid amentoflavone and fisetin is discussed. To this end, both compounds were docked with BRD4 using four docking programs. Results were post-processed with protein-ligand interaction fingerprints. To gain further insights into the binding mode of the two natural products, docking results were further analyzed with molecular dynamics. Results showed that amentoflavone makes numerous contacts in the ZA channel, as previously described for flavonoids and kinase inhibitors. It was also found that amentoflavone can potentially make contacts with non-canonical residues for BET inhibition. Most of these contacts were not observed with fisetin. Based on these results, amentoflavone was tested for BRD4 inhibition, showing activity in the micromolar range. This work may serve as basis for scaffold optimization and further characterization of flavonoids as BET inhibitors.

One of the biggest queries in cognitive sciences is the emergence of consciousness from matter. Modern neurobiological theories of consciousness propose that conscious experience is the result of interactions between large-scale neuronal networks in the brain, traditionally described within the realm of classical physics. Here, we propose a generalized connectionist framework in which the emergence of “conscious networks” is not exclusive of large brain areas, but can be identified in sub-cellular networks exhibiting non-trivial quantum phenomena. The essential feature of such networks is the existence of strong correlations in the system (classical or quantum coherence) and the presence of an optimal point at which the system’s complexity is maximized, expressed either by maximization of the information content in large scale functional networks or by achieving optimal efficiency through the quantum Goldilock effect.

Effect of abrasive shape on the three-body abrasion behaviors of monocrystalline silicon was investigated by Molecular dynamics modelling. The axial ratio of abrasive particle was varied from 1.00, i.e., a complete sphere, to 0.40 to mimic abrasive shape. The particle’s movement turns toward sliding from rolling when the axial ratio becomes less than a critical value 0.46. In the abrasion process, the friction force and normal force showed an approximately sinusoid-like fluctuation for the rolling ellipsoidal particles, while the front cutting of particle caused that friction force increased and became larger than normal force for sliding particles. The phase transformation process was tracked under different particle’ movement patterns. The Si-II and Bct5 phase producing in loading process can partially transform to Si-III/Si-XII phase and backtrack to original crystal silicon under pressure release, which also occurred in the abrasion process. The secondary phase transformation showed difference for particles’ rolling and sliding movements after three-body abrasion. The rolling of particle induced the periodical and inhomogeneous deformation of substrates, while the sliding benefited producing high-quality surface in CMP process. This study aiming to construct more precise model to understand the wear mechanism benefits evaluating the MEMS wear and CMP process of crystal materials.

The Autopoiesis and Cognition Theory (ACT), by Maturana and Varela, based on the notions of Biological Closure and Structural Coupling, is a well-known theory on how to understand biological organization [1, 2, 3]. Although, for example, the Free Energy Principle framework evokes some entailments of autopoiesis in a more formal setting [4, 5]; and ACT has been used in many fields, its impact has been restricted because it lacks quantitative analysis. Here we present a theoretical framework grounded in accepted and well-developed ideas from Mathematics and Physics which advance the understanding of the Principles of Biological Organization under the guidance of Biological Closure and Structural Coupling. The disciplines of Differential Geometry/Topology, Mechanics and Complex Dynamical Systems provide a powerful, elegant, and well-established body of knowledge to support our Biological Organization Principles (BOP) framework. In particular, Stochastic Mechanics and KAM theory (from Kolmogorov, Arnold and Moser theorem) allow us to develop, using the notions of Biological Closure and Structural Coupling, a central core of BOP termed Dynamical Closure Mechanism. Under the proposed framework, a wide variety of bio- logical phenomena can be understood, shedding new light on biological explanations. However, an understanding of biological organization may require the re-evaluation of dogmas on how we think on biology as it seems inescapable that what is needed is an integration of analysis and notions derived from mathematics, physics, and biology to generate a new landscape of ideas.

Microfluidic devices or microdevices refer to systems with a characteristic length in the micrometer range. Systems in this size allow handling small quantities of reagents and samples, with reduced residence time, better control of chemical species concentration, high heat and mass transfers, and high surface/volume ratio. These characteristics led to the application of these microdevices in several areas, such as biological systems, energy, liquid-liquid extraction, food, agricultural sectors, pharmaceuticals, flow chemistry, microreactors, and biodiesel synthesis. Microreactors are devices that have interconnected microchannels, in which small amounts of reagents are manipulated and react for a certain period of time. The traditional characteristics of microreactors are less quantities of reagents and samples, high surface area in relation to volume (10000 m² m^-3), reduction of resistance to heat and mass transfer, reduced reaction times, and narrower residence time distributions. In recent years, several studies have been carried out on biodiesel production in microreactors that explore the influence of operating conditions, mixing and reaction yield, numbering, and especially the microdevices design. Despite all the advantages of microreactors, the literature shows that there are only a few applications on an industrial scale. Two main reasons that hinder the adoption of this technology are the scale-up to a large enough volume to deliver the necessary production capacity and the costs related to industrial manufacturing microreactors. It is often stated that large-scale production of microreactors can be easily achieved by numbering-up. However, researches show that an incredibly high number of microdevices would be needed, which results in a technical unfeasibility and a strong impact on the construction costs of the industrial system. The present review aims to show whether microreactors can replace conventional biodiesel production processes and how this replacement technology could be carried out. The current chapter was divided into the following sections: Introduction, Synthesis and Purification of Biodiesel in Microreactors, Fundamentals of CFD, and Fundamentals of Scale-up. Finally, conclusions and future perspectives are exposed.

Nanoparticles -- regularly patterned or randomly dispersed -- are a key ingredient for emerging technologies in photonics. Of particular interest are scattering and field enhancement effects of metal nanoparticles for energy harvesting and converting systems. An often neglected aspect in the modeling of nanoparticles are light interaction effects beyond classical electrodynamics stemming from electron dynamics in confined and accelerated systems. We give a detailed account on free electron phenomena in metal nanoparticles and discuss analytic expressions, stemming from microscopic and semi-classical theories. These can improve standard computational schemes to produce more reliable results on the optical properties of metal nanoparticles. We combine these solutions into a single framework and study their joint impact on isolated Au, Ag, and Al nanoparticles as well as dimer structures.

Single nucleotide polymorphisms (SNP) are associated with diseases and drug response variabilities in humans. Elucidating the damaging and disease-associated SNPs using wet-laboratory approaches can be challenging and resource-demanding due to the large number of SNPs in the human genome. Due to the growth in the field of computational biology and bioinformatics, algorithms have been developed to help screen and filter out the most deleterious SNPs that are worth considering for wet-laboratory studies. Here we review the existing in-silico based methods used to predict and characterize the effects of SNPs on protein structure and function. This cutting-edge approach will facilitate the search for novel therapeutics, help understand the etiology of diseases and fast-track the personalized medicine agenda.