Motivation and Aim: Intracellular accumulation of organic solutes is one of the strategies used by living organisms to manage the variety of extreme conditions that they encounter in nature. This is also the case for luminous bacteria that inhabit the ocean and can adapt to the variety of temperatures, salinity, concentration of nutrients and toxins. The cells accumulate small molecules to counter the osmotic stress. These low molecular weight protective compounds called osmolytes include polyols, sugars, amino acid derivatives etc. [1]. They are assumed to increase protein stability and prevent protein aggregation. Nevertheless, small organic compounds can cause structural and dynamic changes of proteins and so influence its function, i.e. catalytic activity. For Photobacterium leiognathi luciferase it was revealed that integral intensity of the reaction changes in the presence of different cosolvents as follows: sucrose > sorbitol > glucose > glycerol > ethylene glycol that correlates with the relation of hydrodynamic radii of the cosolvent molecules: 5.2 > 3.9 > 3.6 > 3.1 > 2.6 (Å) [2]. Moreover, sucrose increases catalytic constant of excited state intermediate formation up to 5-fold unlike other cosolvents. Atomic-level molecular dynamics (MD) is useful tool to examine effects of osmolytes on P. leiognathi structure. Thus, the aim of this study was to reveal the structural and dynamic effects caused by osmolytes on the bacterial luciferase using molecular modelling approach. Methods and Algorithms: Three-dimension structure of luciferase from P. leiognathi was obtained by homology modelling using Swiss-model server. Crystal structure of luciferase from Vibrio harveyi (PDB ID: 3FGC) was used as a template. MD simulation of structure of P. leiognathi luciferase in water and in water-cosolvent mixtures (40 w/w %) was performed with GROMACS 2020.4 software. Sucrose, sorbitol, glucose, glycerol and ethylene glycol molecules were used as cosolvents. MD of each modeled system was calculated for 100 ns with three independent runs. Results: To estimate the penetration of cosolvents into active site gorge of P. leiognathi luciferase a radial distribution function (RDF) was calculated between centre of mass of αGLU43 and cosolvent molecules. αGLU43 was chosen as reference point because it locates in the deep of luciferase active site. RDF demonstrates that ethylene glycol and glycerol molecules penetrate deeper the active site of luciferase than other cosolvents that may explain the decrease of integral intensity in the presence of these two osmolytes. Minimum distance of RDF for ethylene glycol is 1,5 Å and for glycerol -1,6 Å, which is close to the value of water (1,45 Å). Sorbitol molecules can be found at 1,7 Å, but this case is extremely rare and the RDF remains close to zero up to the distance of 8 Å relative to the reference point. For glucose and sucrose the RDF function starts to rise from the distance of 4,3 and 5,5 Å respectively, and reaches significant values at >10 Å, i.e. closer to the entrance of the ac...