@phdthesis{Sauer2021,
  author    = {Sauer, Susanne},
  title     = {Implementation and Application of QM/MM Hybrid Methods},
  doi       = {10.25972/OPUS-24321},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-243213},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2021},
  abstract  = {Within this work, an additive and a subtractive QM/MM interface were implemented into CAST. The interactions between QM and MM system are described via electrostatic embedding. Link atoms are used to saturate dangling bonds originating from the separation of QM and MM system. Available energy evaluation methods to be combined include force fields (OPLSAA and AMBER), semi-empirical programs (Mopac and DFTB+), and quantum-chemical methods (from Gaussian, Orca, and Psi4). Both the additive and the subtractive interface can deal with periodic boundary conditions. The subtractive scheme was extended to enable QM/QM, three-layer, and multi-center calculations. Another feature only available within the subtractive interface is the microiteration procedure for local optimizations. The novel QM/MM methods were applied to the investigation of the reaction path for the complex formation between rhodesain and K11777. Benchmark calculations show a very good agreement with results from Gaussian-ONIOM. When comparing the relative energies obtained with different options to a computation where the whole system was treated with the "QM method" DFTB3, the electrostatic embedding scheme with option "delM3" gives the best results. "delM3" means that atoms with up to three bonds distance to the QM region are ignored when creating the external charges. This is done in order to avoid a double counting of Coulomb interactions between QM and MM system. The embedding scheme for the inner system in a three-layer calculation, however, does not have a significant influence on the energies. The same is true for the choice of the coupling scheme: Whether the additive or the subtractive QM/MM interface is applied does not alter the results significantly. The choice of the QM region, though, proved to be an important factor. As can be seen from the comparison of two QM systems of different size, bigger is not always better here. Instead, one has to make sure not to separate important (polar) interactions by the QM/MM border. After this benchmark study with singlepoint calculations, the various possibilities of CAST were used to approximate the solution of a remaining problem: The predicted reaction energy for the formation of the rhodesain-K11777 complex differs significantly depending on the starting point of the reaction path. The reason for this is assumed to be an inadequate adjustment of the environment during the scans, which leads to a better stabilization of the starting structure in comparison to the final structure. The first approach to improve this adjustment was performing the relaxed scan with a bigger QM region instead of the minimal QM system used before. While the paths starting from the covalent complex do not change significantly, those starting from the non-covalent complex become more exothermic, leading to a higher similarity of the two paths. Nevertheless, the difference of the reaction energy is still around 15 kcal/mol, which is far from a perfect agreement. For this reason, Umbrella Samplings were run. Here, the adjustment of the environment is not done by local optimizations like in the scans, but by MD simulations. This has the advantage that the system can cross barriers and reach different local minima. The relative free energies obtained by Umbrella Samplings with suitable QM regions are nearly identical, independently of the starting point of the calculation. Thus, \(\Delta A\) evaluated by these computations can be assumed to reproduce the real energy change best. An MD simulation that was started from the transition state in order to mimic a "real-time" reaction indicates a very fast adjustment of the environment during the formation of the complex. This confirms that Umbrella Sampling is probably better suitable to describe the reaction path than a scan, where the environment can never move strong enough to leave the current local minimum.},
  subject      = {Quantenmechanik},
  language  = {en}
}
@phdthesis{Le2020,
  author    = {Le, Thien Anh},
  title     = {Theoretical investigations of proton transfer and interactions or reactions of covalent and non-covalent inhibitors in different proteins},
  doi       = {10.25972/OPUS-17051},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-170511},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2020},
  abstract  = {Nowadays, computational-aided investigations become an essential part in the chemical, biochemical or pharmaceutical research. With increasing computing power, the calculation of larger biological systems becomes feasible. In this work molecular mechanical (MM) and quantum mechanical approaches (QM) and the combination of both (QM/MM) have been applied to study several questions which arose from different working groups. Thus, this work comprises eight different subjects which deals with chemical reactions or proton transfer in enzymes, conformational changes of ligands or proteins and verification of experimental data. This work firstly deals with reaction mechanisms of aromatic inhibitors of cysteine proteases which can be found in many organisms. These enzymes are responsible for various cancer or diseases as for example Human African Trypanosomiasis (HAT) or the Chagas disease. Aromatic SNAr-type electrophiles might offer a new possibility to covalently modify these proteases. Quantum mechanical calculations have been performed to gain insights into the energetics and possible mechanisms. The next chapter also deals with Trypanosomiasis but the focus was set on a different enzyme. The particularity of Trypanosomiasis is the thiol metabolism which can also be modified by covalent inhibitors. In this context, the wild type and point mutations of the enzyme tryparedoxin have been investigated via molecular dynamic (MD) simulations to examine the influence of specific amino acids in regard to the inhibitor. Experimental data showed that a dimerization of the enzyme occurs if the inhibitor is present. Simulations revealed that the stability of the dimer decreases in absence of the inhibitor and thus confirms these experiments. Further investigations concerning cysteine proteases such as cruzain and rhodesain have been conducted with respect to experimental kinetic data of covalent vinylsulfone inhibitors. Several approaches such as QM or QM/MM calculations and docking, MD or MMPBSA/MMGBSA simulations have been applied to reproduce these data. The utilization of force field approaches resulted in a qualitatively accurate prediction. The kinase AKT is involved in a range of diseases and plays an important role in the formation of cancer. Novel covalent-allosteric inhibitors have been developed and crystallized in complex with AKT. It was shown that depending on the inhibitor a different cysteine residue is modified. To investigate these differences in covalent modification computational simulations have been applied. Enoyl-(acyl carrier) (ENR) proteins are essential in the last step of the fatty acid biosynthesis II (FAS) and represent a good target for inhibition. The diphenylether inhibitor SKTS1 which was originally designed to target the ENR's of Staphylococcus aureus was also crystallized in InhA, the ENR of Mycobacterium tuberculosis (TB). Crystal structures indicate a change of the inhibitor's tautomeric form. This subject was investigated via MD simulations. Results of these simulations confirmed the tautomerization of the inhibitor. This work also deals with the development of a covalent inhibitor originating from a non-covalent ligand. The target FadA5 is an essential enzyme for the degradation of steroids in TB and is responsible for chronic tuberculosis. This enzyme was crystallized in complex with a non-covalent ligand which served as starting point for this study. Computations on QM or QM/MM level and docking and MD simulations have been applied to evaluate potential candidates. The next chapter focuses on the modification of the product spectrum of Bacillus megaterium levansucrase, a polymerase which catalyzes the biosynthesis of fructans. The covalent modification of the wild type or mutants of the enzyme lead to an accumulation of oligosaccharides but also to polymers with higher polymerization degree. To understand these changes in product spectra MD simulations have been performed. Finally, the proton transfer in catalytic cysteine histidine dyads was investigated. The focus was set on the influence of the relaxation of the protein environment to the reaction. Calculations of the enzymes FadA5 and rhodesain revealed that the preferred protonation state of the dyade depends on the protein environment and has an impact on the reaction barrier. Furthermore, the adaptation of the environment to a fixed protonation state was analyzed via MD simulations.},
  subject      = {Computational chemistry},
  language  = {en}
}
@phdthesis{Bellinger2016,
  author    = {Bellinger, Daniel},
  title     = {Implementation of new reaction pathway determining methods and study of solvent effects on the excited state nature of perylene based dyes},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-144435},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2016},
  abstract  = {Two thematic complexes were addressed within this work. One part is related to improvements and new implementations into the CAST program package. Thereby the main focus laid on the delivery of a tool which can be used to characterize complex reactions and their mechanisms. But also within the new force ﬁeld (FF) method (SAPT-FF) within the CAST program, several improvements were made. The second topic is related to the description of dye molecules and their spectral properties. The main focus within these studies was set on the inﬂuence of the environment on these properties. In the ﬁrst topic improvements of the local acting NEB (nudged elastic band) methods were included and the number of available methods was extended. The initial pathway generation was improved by implementing the IDPP (image dependent pair potential) method and a new method was implemented for describing temperature dependent pathways. Additionally, improvements have been made to the optimization routines (global NEB). As a second part the Pathopt (PO) method was considerably improved. In the beginning of the work the original PO idea was used. In this approach one starts with a global optimization on one n-1 dimensional hyperplane which divides the reaction into two sub-areas for obtaining guesses of TSs (transition states). These found TS guesses were used to optimize to the "true" TS. Starting from the optimized ones a relaxation to the next connected minima is done. This idea has been automatically implemented and extended to several number of hyperplanes. In this manner a group of pathsegments is obtained which needs to be connected, but within this work it was realized that such a procedure might be not very eﬃcient. Therefore, a new strategy was implemented which is founded on the same constrained global optimization scheme (MCM) for which the user deﬁnes the number of hyperplanes generated. The number of such generated hyperplanes should be large enough 134 to describe the space between the concerning reactants in a suﬃcient way. The found minima are directly used to built up the reaction pathway. For this purpose a RMSD (root mean square deviation) criterion is used to walk along ways of minimal change from one to another hyperplane. To prove the implementations various test calculations were carried out and extensions included to prove the capabilities of the new strategy. Related to these tests a new strategy for applying the move steps in MCM (Monte Carlo with minimization) was realized which is also related to the question of the coordinates representation. We were able to show that the hopping steps in MCM can be improved by applying Cartesian steps in combination of random dihedral moves with respect to the constraint. In this way it was possible to show that a large variety of systems can be treated. An additional chapter shows the improvements of the SAPT-FF implementation and related test cases. It was possible to treat benzene dimer and cluster systems of diﬀerent sizes consistently also in accordance with high level ab initio based approaches. Furthermore, we showed that the SAPT-FF with the right parameters outperforms the standard AMOEBA implementation which is the basis of the SAPT-FF implementation. In the last three chapters deal with the description of perlyene-based dyes. In the ﬁrst smaller chapter ground state chemistry description of macro cycles of PBI (perylene bisimide) derivatives were investigated. Therefore, AFM (atomic force microscopy) based pictures were explained within our study. The methods to explain aggregation behavior in dependency of the ring size were MD simulations and conﬁguration studies. The last two chapters deal with opto-electronic or photo-physical properties of PBI and PTCDA (perylene-3,4,9,10-tetracarboxylic dianhydride). In detail, we investigated the role of the environment and the aggregate or crystal surrounding by applying diﬀerent models. In that way implicit and explicit solvation models, the size of aggregates and vibration motions were used. In the case of PBI the recent work is found on preliminary studies related to my bachelor thesis and extends it. It was shown that the direct inﬂuence of a polarizable surrounding, as well as explicit inclusion of solvent molecules on the overall description of the excitations and nature of the excited states is weaker as one might expect. However the inclusion of intra-molecular degrees of freedom showed a stronger inﬂuence on the state characteristics and can induce a change of the order of states within the dimer picture. For the PTCDA molecule the main focus was set on the description of the absorption spectrum of crystalline thin ﬁlms. Related to this older works exist which already gave a description and assignment of the absorption band, but are based on diﬀerent approaches compared to the one used in this work. We used the supermolecule ansatz, whereas the environment and diﬀerent aggregate sizes were investigated. Within the dimer based approach we were able to show that using continuum solvation (IEFPCM/COSMO) based description for the environment the relative order of states remains unchanged. Similar to the PBI calculations the inﬂuence of the vibrational motions /distortions is larger. The simulation of the crystal environment by using QM/MM (quantum mechanics/molecular mechanics) approaches delivered that an asymmetric charge distribution might induce a localization of the excitation and a stronger mixing of states. For obtaining further insights we go beyond the dimer picture and aggregates of diﬀerent sizes were used, whereas the simulations up to the octadecamer mono- and even dual-layer stack were carried out. Within these calculations it was shown that the H-coupling is dominating over a weaker J-coupling between diﬀerent stacks. Additionally the calculations based on DFT (density functional theory) and semi-empirics showed that the lowest state in terms of energy are mostly of Frenkel type, whereas the higher lying states are CT ones which mix with embedded Frenkel type states. The ﬁrst band of the absorption spectrum was explained by inclusion of vibrational motions within the stacks which induce an intensity gain of the ﬁrst excited state. This intensity was not explainable by using the undistorted stacks. Also relaxations at the crystal surface might play a role, but are experimentally not explainable.},
  subject      = {Globale Optimierung},
  language  = {en}
}
@phdthesis{Lee2013,
  author    = {Lee, Wook},
  title     = {Computational study on the catalytic mechanism of mtKasA},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-83989},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2013},
  abstract  = {Das Enzym KasA spielt eine entscheidende Rolle in der Biosynthese von Mykols{\"a}uren, den Bausteinen der Zellw{\"a}nde von Mycobacteriumtuberculosis. Dessen essentielle Notwendigkeit zeigt sich bei Abwesenheit von KasA in einer Zelllyse (Aufl{\"o}sung von Zellen) bei Mycobacteriumtuberculosis. Durch seine Bedeutung f{\"u}r Mycobacteriumtuberculosis, dem Erreger von Tuberkulose und damit der zweith{\"a}ufigsten Todesursache durch Infektionskrankheiten, stellt KasA ein vielversprechendes Ziel f{\"u}r die Entwicklung neuer Medikamente gegen Tuberkulose dar. Durch das Auftreten von extensiv resistenten St{\"a}mmen welche die meisten bekannten Antibiotika zur Bek{\"a}mpfung von Tuberkulose inaktivieren wird es dringend notwendig neue Medikamente gegen Tuberkulose zu entwickeln. In Kapitel 3.1 wird der Protonierungszustand der katalytischen Reste im Ruhezustand untersucht. F{\"u}r diese Untersuchungen wurden Free Energy Perturbation (FEP) Rechnungen und MD Simulationen verwendet. Die Ergebnisse zeigten, dass der zwitterionische Zustand am wahrscheinlichsten ist. Um diese Aussage mit weiteren handfesten Daten zu untermauern wurden Potential(hyper)fl{\"a}chen (PES) f{\"u}r den Protonentransfer zwischen neutralen und zwitterionischen Zustand mit Hilfe von QM/MM Methoden berechnet. Durch die starke Abh{\"a}ngigkeit der QM/MM Optimierung von der Ausgangsstruktur war es nicht m{\"o}glich konsistente Ergebnisse f{\"u}r diese Berechnungen zu bekommen. Um dieses Problem zu umgehen wurde ein auf QM/MM basierendes Umbrella Sampling mit Semiempirischen Methoden (RM1) durchgef{\"u}hrt. Die sich daraus ergebende PMF Fl{\"a}che zeigt das der zwitterionische Zustand stabiler ist als der neutrale Zustand. In Kapitel 3.2 wurde der Protonierungszustand der entsprechenden Reste im Acyl-Enzym Zustand untersucht. Im Unterschied zu anderen katalytischen Resten ist der Protonierungszustand von His311 ist nicht eindeutig im Acyl-Enzym Zustand und es ergeben sich aus den verschiedenen Protonierungszust{\"a}nden verschiedene Decarboxylierungsmechanismen. Um den wahrscheinlichsten Protonierungszustand bez{\"u}glich der freien Energie zu bestimmen wurden FEP Rechnungen durchgef{\"u}hrt. Die Ergebnisse zeigen, dass der pKa Wert an Nδ betr{\"a}chtlich durch die Enzymumgebung verringert wird, w{\"a}hrend dies f{\"u}r Nε nicht der Fall ist. Zus{\"a}tzlich dazu wurden die PMF Profile f{\"u}r den Protonentransfer zwischen Lys340 und Glu354 mit der QM/MM basierten Umbrella Sampling Methode berechnet. Die Ergebnisse zeigen, dass das Lys340/Glu354 Paar eher neutral als ionisch ist, wenn His311 an Nε protoniert ist. Ein relativ hoher ionischer Charakter des Lys340/Glu354 Paares, wenn His311 doppelt protoniert ist, gibt einen wertvollen Einblick in die Rolle welche das Lys340/Glu354 Paar beim verschieben des Protonierungszustandes von Nδ zu Nε im His311 nach dem Acyltransferschritt spielt. Die Ergebnisse zeigen, dass His311 neutral und an Nε protoniert ist. Ebenso ist das Lys340/Glu354 Paar neutral im Acyl-Enzym Zustand. Diese berechneten Ergebnisse f{\"u}hren zu dem Schluss, dass die Decarboxylierung durch ein Oxyanion Loch erleichtert wird welches aus zwei katalytischen Histidin Resten besteht. In Kapitel 3.3 wurde der Protonierungszustand der katalytischen Reste im Ruhezustand erneut untersucht da eine aktuelle Benchmarkstudie zeigte, dass die verwendete Semiempirische Methode (RM1) in Kapitel 3.1 dazu tendiert die Stabilisation des zwitterionischen Zustandes zu {\"u}bersch{\"a}tzen. Auch wurde in Kapitel 3.1 das Lys340/Glu354 Paar als rein ionisch angesehen, w{\"a}hrend sich in Kapitel 3.2 herausstellte, dass es sich um eine Mischung aus neutralen und ionischen Charakter handelt. Die neuen Untersuchungen beinhalten eine gr{\"o}ßere QM Region inklusive des Lys340/Glu354 Paares. Der daf{\"u}r verwendete BLYP/6-31G** Ansatz ist ausreichend akkurat f{\"u}r die aktuelle Fragestellung, was durch Vergleichsrechnungen bewiesen wurde. Die neuen Ergebnisse der QM/MM MD und FEP Rechnungen deuten an, dass die katalytischen Reste im Ruhezustand h{\"o}chst wahrscheinlich neutral vorliegen. Dies wiederum f{\"u}hrt zu der Frage wie KasA aktiviert werden kann um die katalytische Reaktion zu initiieren. Auf der Basis der Ergebnisse der MD Simulationen und FEP Rechnungen f{\"u}r den His311Ala Mutanten in Kapitel 3.1 stellten wir die Hypothese auf, dass die offene Konformation von Phe404 die Aktivierung der katalytischen Reste durch die (Aus)bildung einer starken Wasserstoffbindung einleitet. Die QM/MM MD Simulation best{\"a}tigt dass diese Aktivierung der katalytischen Reste durch die offene Konformation des Phe404 bewerkstelligt werden kann. Das entsprechende auf Kraftfeld basierende PMF Profil zeigt auch, dass dieser Konformationswechsel energetisch realisierbar ist. Die Verteilung der hydrophilen und hydrophoben Reste in der Malonyl Bindungstasche in Verbindung mit unseren berechneten Ergebnissen geben einen Einblick in den detaillierten},
  subject      = {Tuberkelbakterium},
  language  = {en}
}