Due to the specificity of their structure, protein systems are adapted to carry various ligands. The structure of many proteins potentially allows for two types of immobilization of a therapeutic agent, either on the outer surface of the protein or within the protein structure. The existence of two active sites in BSA’s structure, the so-called Sudlow I and II, was confirmed. The conducted research involved determining the effectiveness of BSA as a potential carrier of 5-fluorouracil (5FU). 5-fluorouracil is a broad-spectrum anticancer drug targeting solid tumors. The research was carried out to estimate the physicochemical properties of the system using complementary measurement techniques. The optimization of the complex formation conditions made it possible to obtain significant correlations between the form of the drug and the effective localization of the active substance in the structure of the protein molecule. The presence of two amino groups in the 5FU structure contributes to the deprotonation of the molecule at high pH values (pH > 8) and the transition to the anionic form (AN1 and AN3). To investigate the binding affinity of the tautomeric form with BSA, UV-vis absorption, fluorescence quenching, zeta potential, QCM-D, and CD spectroscopic studies were performed. The experimental research was supported by molecular dynamics (MD) simulations and molecular docking. The simulations confirm the potential location of 5FU tautomers inside the BSA structure and on its surface.

The unique structure of G4.0 PAMAM dendrimers allows a drug to be enclosed in internal spaces or immobilized on the surface. In the conducted research, the conditions for the formation of the active G4.0 PAMAM complex with doxorubicin hydrochloride (DOX) were optimized. The physicochemical properties of the system were monitored using dynamic light scattering (DLS), circular dichroism (CD), and fluorescence spectroscopy. The Quartz Crystal Microbalance with Dissipation Monitoring (QCM-D) method was chosen to determine the preferential conditions for the complex formation. The highest binding efficiency of the drug to the cationic dendrimer was observed under basic conditions when the DOX molecule was deprotonated. The decrease in the zeta potential of the complex confirms that DOX immobilizes through electrostatic interaction with the carrier’s surface amine groups. The binding constants were determined from the fluorescence quenching of the DOX molecule in the presence of G4.0 PAMAM. The two-fold way of binding doxorubicin in the structure of dendrimers was visible in the Isothermal calorimetry (ITC) isotherm. Fluorescence spectra and release curves identified the reversible binding of DOX to the nanocarrier. Among the selected cancer cells, the most promising anticancer activity of the G4.0-DOX complex was observed in A375 malignant melanoma cells. Moreover, the preferred intracellular location of the complexes concerning the free drug was found, which is essential from a therapeutic point of view.

Understanding the mechanism of adsorption of Rhodamine 6G (R6G) to various crystal structures of silica nanoparticles (SNPs) is important to elucidate the impact of dye size when measuring the size of the dye-SNP complex via the time-resolved fluorescence anisotropy method. In this work, molecular dynamics (MD) simulations were used to get an insight into the R6G adsorption process, which cannot be observed using experimental methods. It was found that at low pH, α-Cristobalite structured SNPs have a strong affinity to R6G; however, at high pH, more surface silanol groups undergo ionization when compared with α-Quartz, preventing the adsorption. Therefore, α-Quartz structured SNPs are more suitable for R6G adsorption at high pH than the α-Cristobalite ones. Furthermore, it was found that stable adsorption can occur only when the R6G xanthene core is oriented flat with respect to the SNP surface, indicating that the dye size does not contribute significantly to the measured size of the dye-SNP complex. The requirement of correct dipole moment orientation indicates that only one R6G molecule can adsorb on any sized SNP, and the R6G layer formation on SNP is not possible. Moreover, the dimerization process of R6G and its competition with the adsorption has been explored. It has been shown that the highest stable R6G aggregate is a dimer, and in this form, R6G does not adsorb to SNPs. Finally, using steered molecular dynamics (SMD) with constant-velocity pulling, the binding energies of R6G dimers and R6G complexes with both α-Quartz and α-Cristobalite SNPs of 40 Å diameter were estimated. These confirm that R6G adsorption is most stable on 40 Å α-Quartz at pH 7, although dimerization is equally possible.

π–π stacking are omnipresent interactions, crucial in many areas of chemistry, and often studied using quantum chemical methods. Here, we report a simple and computationally efficient method of estimating the binding energies of stacked polycyclic aromatic hydrocarbons based on steered molecular dynamics. This method leverages the force field parameters for accurate calculation. The presented results show good agreement with those obtained through DFT at the ωB97X-D3/cc-pVQZ level of theory. It is demonstrated that this force field-driven SMD method can be applied to other aromatic molecules, allowing insight into the complexity of the stacking interactions and, more importantly, reporting π–π stacking energy values with reasonable precision.

With the ever-growing need to study systems of increased size and complexity, modern density functional theory (DFT) methods often encounter problems arising from the growing computational demands. In this work, we have presented a comprehensive DFT validation of the steered molecular dynamics (SMD) approach in estimating the binding energies of aromatic dimers. By performing DFT calculations on optimized and unoptimized anthracene and rhodamine 6G (R6G) dimers using functionals with progressively enhanced exchange-correlation energy description and comparing the obtained results with SMD-predicted values, it was found that SMD predictions are in good agreement with the results obtained from hybrid DFT calculations. The average binding energies for optimized anthracene dimers were found to be 6.46 kcal/mol using DFT at ωB97X-D4/def2-QZVPP and 7.64 ± 1.61 kcal/mol as predicted by the SMD. For the R6G H-type dimer, the binding energies were 17.48 and 19.02 ± 2.22 kcal/mol, respectively. The study also revealed that due to the lack of explicit terms accounting for electron-electron interactions in MD force fields, the proposed method tends to overbind dimers. It is anticipated that the presented method can be applied to more complex dimers, potentially accelerating the calculations of binding energies. Moreover, this study further validates the accuracy of the CHARMM36 FF.

We investigate the nanometrology of sub-nanometre particle sizes in industrially manufactured sodium silicate liquors at high pH using time-resolved fluorescence anisotropy. Rather than the previous approach of using a single dye label, we investigate and quantify the advantages and limitations of multiplexing two fluorescent dye labels. Rotational times of the non-binding rhodamine B and adsorbing rhodamine 6G dyes are used to independently determine the medium microviscosity and the silicate particle radius, respectively. The anisotropy measurements were performed on the range of samples prepared by diluting the stock solution of silicate to concentrations ranging between 0.2 M and 2 M of NaOH and on the stock solution at different temperatures. Additionally, it was shown that the particle size can also be measured using a single excitation wavelength when both dyes are present in the sample. The recovered average particle size has an upper limit of 7.0 ± 1.2 Å. The obtained results were further verified using small-angle X-ray scattering, with the recovered particle size equal to 6.50 ± 0.08 Å. To disclose the impact of the dye label on the measured complex size, we further investigated the adsorption state of rhodamine 6G on silica nanoparticles using molecular dynamics simulations, which showed that the size contribution is strongly impacted by the size of the nanoparticle of interest. In the case of the higher radius of curvature (less curved) of larger particles, the size contribution of the dye label is below 10%, while in the case of smaller and more curved particles, the contribution increases significantly, which also suggests that the particles of interest might not be perfectly spherical.