The application of commercial herbicide restricts weed growth and significantly improves control over crop vitality and yield. Despite their utility in the agriculture sector, herbicides have the potential to contaminate local water sources. To minimize environmental impacts, the development of efficient separation processes to clean-up contaminated water bodies is necessary. However, complex speciation and conformational flexibility in the condensed phase poses a significant challenge. In this work, we investigate structure and speciation of three common organic herbicides (glyphosate, atrazine, and metolachlor) in aqueous solution. We employ the PBE-D3 density functional to perform ab initio molecular dynamics (MD) simulations in the canonical and isothermal-isobaric ensembles. We analyze MD trajectories to understand hydrogen bonding dynamics and lifetime as well as diffusional and vibrational characteristics. To enhance configurational sampling, we conduct metadynamics simulations to obtain the free energies of dissociation and intramolecular proton transfer of glyphosate.

Title: Raman optical activity and conformational flexibility of peptides in solution Author: Jana Hrudíková Department: Institute of Physics of Charles University Supervisor: Doc. RNDr. Vladimír Baumruk, DrSc. Supervisor's e-mail address: baumruk@karlov.mff.cuni.cz Abstract: Molecular flexibility can significantly modify Raman and ROA spectral intensities, band positions and the ROA signs. Taking into account dynamic aspects of behavior of studied molecules in solution via conformational averaging therefore seems to be crucial for spectral interpretation. The first of studied models, histidine, plays an important role in metallo-enzymatic reactions and peptide folding, due to its imidazole ring. ROA spectra of His at different pH, His complexed with Cu2+ and dipeptides His- Gly and Gly-His were recorded on the spectrometer built at the Institute of Physics of the Charles University as a first step of the subsequent study. The second studied system, a cyclic hexapeptide c-(Phe-D-Pro-Gly-Arg-Gly-Asp), serves as a convenient model for β- hairpin and anti-parallel β-sheet. It was previously studied by means of VCD and IR. From molecular dynamics simulations 10 peptide geometries were selected for spectral modeling. The Raman and ROA spectra were calculated ab initio. For a model fragment Phe-D-Pro, which...

Métodos computacionais são usados para gerar estruturas de complexo proteína-ligante e estimar suas afinidades. Esse trabalho investigou como as diferentes representações da flexibilidade proteica afetam as poses obtidas por ancoragem molecular e as afinidades atribuídas a essas poses. Os mutantes L99A e L99A/M102Q da lisozima T4 foram escolhidos como sistemas modelo. Um descritor para predição de afinidades baseado na aproximação de energia de interação linear (LIE) foi parametrizado especificamente para ligantes da lisozima e foi usado para estimar as afinidades. A proteína foi representada como um grupo de estruturas cristalográficas ou de estruturas de trajetória de dinâmica molecular. O campo de força OPLS-AA para modelar a proteína e os ligantes e a aproximação de Born generalizada para modelar o solvente foram empregados. O descritor de afinidades parametrizado resultou em desvios médios entre afinidades experimentais e calculadas de 1,8 kcal/mol para um conjunto de testes. O descritor teve desempenho satisfatório na separação entre poses cristalográficas e poses falso-positivo e na identificação de poses falso-positivo. Experimentos de agrupamento de complexos realizados com o objetivo de reduzir o custo computacional para estimar afinidades apresentaram resultados insatisfatórios. As melhores aproximações da teoria do ligante implícito propostas aqui para estimar afinidades consideram conjuntos de estruturas de receptor com o mesmo peso. Configurações de ligante também apresentam o mesmo peso ou são dominadas por uma única configuração. A representação da flexibilidade requer um tratamento estatístico adequado para estimativa de afinidades. Aqui, a associação entre LIE e a teoria do ligante implícito mostrou-se frutífera. / Computational methods are used to generate protein-ligand complex structures and estimate their binding affinities. This work investigated how different representations of protein flexibility affect poses obtained by molecular docking and the affinities attributed to these poses. T4 lysozyme mutants L99A and L99A/M102Q were chosen as model systems. A descriptor for prediction of affinities based on linear interaction energy (LIE) approximation was parametrized specifically to lysozyme ligands and was used to estimate affinities. The protein was represented as a group of crystal structures or as structures from a molecular dynamics trajectory. OPLS-AA force field was used to model protein and ligands and the Generalized Born approximation was used to model solvent. The parametrized affinity descriptor resulted in average deviations between experimental and calculated affinities of 1.8 kcal/mol for a test set. Descriptor performance was satisfactory in the separation between crystal poses and false-positive ones and in the identification of false-positive poses. Clustering of complexes was tried out to reduce computational cost to estimate affinities, but results were poor. The best approximations to the implicit ligand theory proposed here in order to estimate affinities consider groups of receptor structures with the same weight. Ligand configurations also have the same weight or are dominated by only one configuration. The representation of protein flexibility requires an adequate statistical treatment when used to estimate affinities. Here, the linking between LIE and the implicit ligand theory proved itself useful.

Métodos computacionais são usados para gerar estruturas de complexo proteína-ligante e estimar suas afinidades. Esse trabalho investigou como as diferentes representações da flexibilidade proteica afetam as poses obtidas por ancoragem molecular e as afinidades atribuídas a essas poses. Os mutantes L99A e L99A/M102Q da lisozima T4 foram escolhidos como sistemas modelo. Um descritor para predição de afinidades baseado na aproximação de energia de interação linear (LIE) foi parametrizado especificamente para ligantes da lisozima e foi usado para estimar as afinidades. A proteína foi representada como um grupo de estruturas cristalográficas ou de estruturas de trajetória de dinâmica molecular. O campo de força OPLS-AA para modelar a proteína e os ligantes e a aproximação de Born generalizada para modelar o solvente foram empregados. O descritor de afinidades parametrizado resultou em desvios médios entre afinidades experimentais e calculadas de 1,8 kcal/mol para um conjunto de testes. O descritor teve desempenho satisfatório na separação entre poses cristalográficas e poses falso-positivo e na identificação de poses falso-positivo. Experimentos de agrupamento de complexos realizados com o objetivo de reduzir o custo computacional para estimar afinidades apresentaram resultados insatisfatórios. As melhores aproximações da teoria do ligante implícito propostas aqui para estimar afinidades consideram conjuntos de estruturas de receptor com o mesmo peso. Configurações de ligante também apresentam o mesmo peso ou são dominadas por uma única configuração. A representação da flexibilidade requer um tratamento estatístico adequado para estimativa de afinidades. Aqui, a associação entre LIE e a teoria do ligante implícito mostrou-se frutífera. / Computational methods are used to generate protein-ligand complex structures and estimate their binding affinities. This work investigated how different representations of protein flexibility affect poses obtained by molecular docking and the affinities attributed to these poses. T4 lysozyme mutants L99A and L99A/M102Q were chosen as model systems. A descriptor for prediction of affinities based on linear interaction energy (LIE) approximation was parametrized specifically to lysozyme ligands and was used to estimate affinities. The protein was represented as a group of crystal structures or as structures from a molecular dynamics trajectory. OPLS-AA force field was used to model protein and ligands and the Generalized Born approximation was used to model solvent. The parametrized affinity descriptor resulted in average deviations between experimental and calculated affinities of 1.8 kcal/mol for a test set. Descriptor performance was satisfactory in the separation between crystal poses and false-positive ones and in the identification of false-positive poses. Clustering of complexes was tried out to reduce computational cost to estimate affinities, but results were poor. The best approximations to the implicit ligand theory proposed here in order to estimate affinities consider groups of receptor structures with the same weight. Ligand configurations also have the same weight or are dominated by only one configuration. The representation of protein flexibility requires an adequate statistical treatment when used to estimate affinities. Here, the linking between LIE and the implicit ligand theory proved itself useful.

Title: Theory and application of optical spectroscopic methods for structural molecular studies Author: RNDr. Jana Hudecová Department / Institute: Institute of Organic Chemistry and Biochemistry Supervisor of the doctoral thesis: Prof. RNDr. Petr Bouř, DSc. Abstract: In the thesis, methods of the chiroptical spectroscopy (Raman optical activity, electronic and vibrational circular dichroism, circularly polarized luminescence) were utilized to obtain information on structure of chiral molecules. In four main projects, we focused on improving accuracy of quantum-chemical computations used for interpretation of experimental spectra by including anharmonic effects, solvent, molecular flexibility and dynamics. In the first project, the normal mode geometry optimization method was investigated and a suitable frequency limit providing realistic vibrational band broadening was found. Then the ability of harmonic and anharmonic computational approaches to describe the C-H stretching vibrations was explored for three terpene molecules and four spectroscopic methods. In the third project, we estimated the role of dispersion forces and different organic solvents for conformer equilibria and dynamics of cyclic dipeptides containing tryptophan. In the last project, circularly polarized luminiscence spectra, which were...