Study of absorption and emission spectra of substituted terthiophene compounds by methods of theoretical chemistry.

Terthiophene derivatives attract interest due to their prospective applications in optoelectronic or sensor devices. Due to their nontoxicity they can be considered as suitable candidates in biological applications. Supramolecular organization of the matter is one of the most interesting topics in contemporary materials science. Amphiphilic chromophores based on substituted terthiophenes are capable of self-assembly into supramolecular architectures. In this work, we aim at simulation of the spectral properties of terthiophene with oligo(ethylene oxide) substituents by the methods of quantum chemistry and molecular dynamics (MD). The potential energy surface (PES) of this molecule was determined by the methods of density functional theory (DFT) for the ground state and time-dependent density-functional theory (TD-DFT) for the excited state. MD simulations in water than revealed the most frequented molecular conformations in both these states. Absorption and fluorescence spectra were determined for all these conformations, including the surrounding water molecules, using TD-DFT and averaged over the conformation space to obtain the final absorption and fluorescence spectrum. The calculated spectra were compared with their experimental counterparts and the differences were discussed in context of the supramolecular structure revealed by confocal microscopy. In spite of its simplicity, this approach provides a satisfactory approximation of absorption and fluorescent spectra of these molecules obtained by computational methods.

Hyaluronan oligosaccharides form double-helical duplexes in water:1,4-dioxane mixed solvent.

Hyaluronic acid (HA) is a hydrophilic natural polysaccharide consisting of alternating monosaccharide units of glucuronic acid and N-acetyl glucosamine. In aqueous solutions the electrostatic repulsion of the carboxylate groups hampers the formation of supermolecular structures that can be partially stabilized by the addition of salt. Increased permittivity of the mixed water:organic solvents causes better compensation of the negative charge of HA chains by dissolved cations which changes their interactions with other molecules. In this study we simulate interactions of two HA chains in water:1,4-dioxane and water:tert-butanol mixed solvents with varying NaCl concentrations using molecular dynamics (MD). Anti-parallel double-helix-like duplexes are formed in NaCl-containing water:1,4-dioxane mixture and remain stable even when NaCl is removed. Parallel duplexes separate after a short time. In water:tert-butanol analogous duplexes are unstable. Stability of HA duplexes is thus determined by the solvent composition and the ability of its components to separate in the solvation shell of HA molecules, as well as by the mutual orientation of the chains.

Structure and dynamics of the hyaluronan oligosaccharides and their solvation shell in water: organic mixed solvents.

Hyaluronan (HA) is a natural polysaccharide occurring ubiquitously in the connective tissues of vertebrates widely used in the cosmetic and pharmaceutic industries. In numerous applications HA oligosaccharides are being chemically modified using reactions incompatible with aqueous solutions, often carried out in water:organic mixed solvents. We carry out molecular-dynamics (MD) simulations of HA oligosaccharides in water:1,4-dioxane and water:tert-butanol mixtures of different compositions. HA molecule causes a separation of the solvent components in its surroundings, especially in tert-butanol containing solutions, constituting thus a solvation shell enriched by water. Furthermore, interactions with ions are stronger than in pure water and depend on the solvent composition. Consequently, the dynamics of the HA chain varies with the solvent composition and causes observable conformational changes of the HA oligosaccharide. Composition of mixed solvents thus enables us to modify the interaction of HA with other molecules as well as its reactivity.

Salt-dependent intermolecular interactions of hyaluronan molecules mediate the formation of temporary duplex structures.

Hyaluronic acid (HA) is a natural polysaccharide present in the connective tissues of vertebrates, often used in the cosmetics and pharmaceutical industries. HA is a strongly hydrophilic macromolecule forming highly swollen random coils in aqueous solutions. Although some authors reported the secondary and tertiary structures of HA chain, others brought convincing evidence contradicting this hypothesis. This study aims at investigation of the stability and dynamics of the temporary duplex HA structures at different NaCl concentrations by molecular-dynamics (MD) simulations. The tendency to duplex formation grows with NaCl concentration reaching its maximum at 0.6 M. This profile is a result of two counteracting NaCl-concentration dependent phenomena, the growing electrostatic-repulsion screening on one side and the disturbance of hydrogen-bonds formation on the other side. Although the weak intermolecular attraction cannot lead to long-lived secondary and tertiary structures, it may influence the properties of large HA macromolecules and concentrated HA solutions.

The rate and evenness of the substitutions on hyaluronan grafted by dodecanoic acid influenced by the mixed-solvent composition.

In this work, low molecular weight (17 kDa) hyaluronan was modified by dodecanoyl substituents. The activation of dodecanoic acid was mediated by benzoyl chloride towards the preparation of a mixed anhydride, which reacts in a second step with HA in water mixed with an organic solvent. The effect of the cosolvent was studied and showed an even distribution of substituents and higher reaction rate in water: 1,4-dioxane compared to water:tert-butanol where substituents occupy adjacent positions. The chemical characterization of the prepared derivatives was elucidated by NMR, FTIR spectroscopy, thermal analyses, and gas chromatography, while the distribution of substituents was evaluated by enzymatic degradation. Molecular-dynamics simulations reveal opposite solvent separations around HA and dodecanoyl chains, that is stronger in water:tert-butanol solution. The resulting incompatibility of solvation-shells of the two entities repels the reaction intermediates from the HA chain and drives them towards the already bound substituents, explaining the observed differences in the distribution evenness. Thus, the influence of the solvent on the reaction selectivity is observed by shielding reactive sites around HA. Therefore, a control of the distribution of the substituents was obtained by defining the concentration of HA and used cosolvent.

Effect of solvent and ions on the structure and dynamics of a hyaluronan molecule.

Hyaluronic acid (hyaluronan, HA) is a negatively charged polysaccharide forming highly swollen random coils in aqueous solutions. Their size decreases along with growing salt concentration, but the mechanism of this phenomenon remains unclear. We carry out molecular-dynamics simulations of a 48-monosaccharide HA oligomer in varying salt concentration and temperature. They identify the interaction points of Na+ ions with the HA chain and reveal their influence on the HA solvation-shell structure. The salt-dependent variation of the molecular size does not consist in the distribution of the dihedral angles of the glycosidic connections but is driven by the random flips of individual dihedral angles, which cause the formation of temporary hairpin-like structures effectively shortening the chain. They are induced by the frequency of cation-chain interactions that grow with the salt concentration, but is reduced by the simultaneous decrease of ions' activities. This leads to an anomalous random-coil shrinkage at 0.6 M salt concentration.

Hyaluronan random coils in electrolyte solutions-a molecular dynamics study.

A computational method of modeling random coils of hyaluronan was developed based on the molecular-dynamics simulations. An oligosaccharide of 48 monosaccharide units was equilibrated within a 70-100ns simulation and randomly chosen pieces of this molecule from different simulation frames were combined to constitute a long polysaccharide chain, both for hyaluronan and its non-ionic analog containing glucose instead of glucuronic acid. The dihedral angles of the glycoside connections of the pieces obeyed the statistics deduced from the simulation. The simulations were performed at various concentrations of NaCl and MgCl2. The calculated radii of gyration show a striking agreement with experimental data from the literature and indicate a key importance of the polymer-ion interactions for the random-coil conformation, but a low influence of the excluded volume of the chain and the carboxylate-groups repulsion. The method has thus the potential to become a versatile tool of modeling macromolecules of various semirigid polymers.

Equilibria of oligomeric proteins under high pressure - A theoretical description.

High pressure methods have become a useful tool for studying protein structure and stability. Using them, various physico-chemical processes including protein unfolding, aggregation, oligomer dissociation or enzyme-activity decrease were studied on many different proteins. Oligomeric protein dissociation is a process that can perfectly utilize the potential of high-pressure techniques, as the high pressure shifts the equilibria to higher concentrations making them better observable by spectroscopic methods. This can be especially useful when the oligomeric form is highly stable at atmospheric pressure. These applications may be, however, hindered by less intensive experimental response as well as interference of the oligomerization equilibria with unfolding or aggregation of the subunits, but also by more complex theoretical description. In this study we develop mathematical models describing different kinds of oligomerization equilibria, both closed (equilibrium of monomer and the highest possible oligomer without any intermediates) and consecutive. Closed homooligomer equilibria are discussed for any oligomerization degree, while the more complex heterooligomer equilibria and the consecutive equilibria in both homo- and heterooligomers are taken into account only for dimers and trimers. In all the cases, fractions of all the relevant forms are evaluated as functions of pressure and concentration. Significant points (inflection points and extremes) of the resulting transition curves, that can be determined experimentally, are evaluated as functions of pressure and/or concentration. These functions can be further used in order to evaluate the thermodynamic parameters of the system, i.e. atmospheric-pressure equilibrium constants and volume changes of the individual steps of the oligomer-dissociation processes.

Pressure induced structural changes and dimer destabilization of HIV-1 protease studied by molecular dynamics simulations.

High-pressure methods have become attractive tools for investigation of the structural stability of proteins. Besides protein unfolding, dimerization can be studied this way, too. HIV-1 protease is a convenient target of experimental and theoretical high-pressure studies. In this study molecular-dynamics simulations are used to predict the response of HIV-1 protease to the pressure of 0.1 to 600 MPa. The protease conformation of both the monomer and the dimer is highly rigid changing insignificantly with growing pressure. Hydrophobicity of the protease decreases with increasing pressure. Water density inside the active-site cavity grows from 87% to 100% of the bulk water density within the pressure range. The dimer-dissociation volume change is negative for most of the pressure ranges with the minimum of -105 ml mol(-1), except for a short interval of positive values at low pressures. The dimer is thus slightly stabilized up to 160 MPa, but strongly destabilized by higher pressures.