@article{WohlgemuthMitric2016,
  author    = {Wohlgemuth, Matthias and Mitric, Roland},
  title     = {Photochemical Chiral Symmetry Breaking in Alanine},
  series = {Journal of Physical Chemistry A},
  volume    = {45},
  journal   = {Journal of Physical Chemistry A},
  number    = {120},
  doi       = {10.1021/acs.jpca.6b07611},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-158557},
  pages     = {8976-8982},
  year      = {2016},
  abstract  = {We introduce a general theoretical approach for the simulation of photochemical dynamics under the influence of circularly polarized light to explore the possibility of generating enantiomeric enrichment through polarized-light-selective photochemistry. The method is applied to the simulation of the photolysis of alanine, a prototype chiral amino acid. We show that a systematic enantiomeric enrichment can be obtained depending on the helicity of the circularly polarized light that induces the excited-state photochemistry of alanine. By analyzing the patterns of the photoinduced fragmentation of alanine we find an inducible enantiomeric enrichment up to 1.7\%, which is also in good correspondence to the experimental findings. Our method is generally applicable to complex systems and might serve to systematically explore the photochemical origin of homochirality.},
  language  = {en}
}
@article{ShityakovSalvadorPastorinetal.2015,
  author    = {Shityakov, Sergey and Salvador, Ellaine and Pastorin, Giorgia and F{\"o}rster, Carola},
  title     = {Blood-brain barrier transport studies, aggregation, and molecular dynamics simulation of multiwalled carbon nanotube functionalized with fluorescein isothiocyanate},
  series = {International Journal of Nanomedicine},
  volume    = {10},
  journal   = {International Journal of Nanomedicine},
  doi       = {10.2147/IJN.S68429},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-149233},
  pages     = {1703-1713},
  year      = {2015},
  abstract  = {In this study, the ability of a multiwalled carbon nanotube functionalized with fluorescein isothiocyanate (MWCNT-FITC) was assessed as a prospective central nervous system-targeting drug delivery system to permeate the blood-brain barrier. The results indicated that the MWCNT-FITC conjugate is able to penetrate microvascular cerebral endothelial monolayers; its concentrations in the Transwell® system were fully equilibrated after 48 hours. Cell viability test, together with phase-contrast and fluorescence microscopies, did not detect any signs of MWCNT-FITC toxicity on the cerebral endothelial cells. These microscopic techniques also revealed presumably the intracellular localization of fluorescent MWCNT-FITCs apart from their massive nonfluorescent accumulation on the cellular surface due to nanotube lipophilic properties. In addition, the 1,000 ps molecular dynamics simulation in vacuo discovered the phenomenon of carbon nanotube aggregation driven by van der Waals forces via MWCN-TFITC rapid dissociation as an intermediate phase.},
  language  = {en}
}
@article{SarukhanyanShityakovDandekar2020,
  author    = {Sarukhanyan, Edita and Shityakov, Sergey and Dandekar, Thomas},
  title     = {Rational drug design of Axl tyrosine kinase type I inhibitors as promising candidates against cancer},
  series = {Frontiers in Chemistry},
  volume    = {7},
  journal   = {Frontiers in Chemistry},
  number    = {920},
  issn      = {2296-2646},
  doi       = {10.3389/fchem.2019.00920},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-199505},
  year      = {2020},
  abstract  = {The high level of Axl tyrosine kinase expression in various cancer cell lines makes it an attractive target for the development of anti-cancer drugs. In this study, we carried out several sets of in silico screening for the ATP-competitive Axl kinase inhibitors based on different molecular docking protocols. The best drug-like candidates were identified, after parental structure modifications, by their highest affinity to the target protein. We found that our newly designed compound R5, a derivative of the R428 patented analog, is the most promising inhibitor of the Axl kinase according to the three molecular docking algorithms applied in the study. The molecular docking results are in agreement with the molecular dynamics simulations using the MM-PBSA/GBSA implicit solvation models, which confirm the high affinity of R5 toward the protein receptor. Additionally, the selectivity test against other kinases also reveals a high affinity of R5 toward ABL1 and Tyro3 kinases, emphasizing its promising potential for the treatment of malignant tumors.},
  language  = {en}
}
@article{SarukhanyanShityakovDandekar2018,
  author    = {Sarukhanyan, Edita and Shityakov, Sergey and Dandekar, Thomas},
  title     = {In silico designed Axl receptor blocking drug candidates against Zika virus infection},
  series = {ACS Omega},
  volume    = {3},
  journal   = {ACS Omega},
  number    = {5},
  doi       = {10.1021/acsomega.8b00223},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-176739},
  pages     = {5281-5290},
  year      = {2018},
  abstract  = {After a large outbreak in Brazil, novel drugs against Zika virus became extremely necessary. Evaluation of virus-based pharmacological strategies concerning essential host factors brought us to the idea that targeting the Axl receptor by blocking its dimerization function could be critical for virus entry. Starting from experimentally validated compounds, such as RU-301, RU-302, warfarin, and R428, we identified a novel compound 2′ (R428 derivative) to be the most potent for this task amongst a number of alternative compounds and leads. The improved affinity of compound 2′ was confirmed by molecular docking as well as molecular dynamics simulation techniques using implicit solvation models. The current study summarizes a new possibility for inhibition of the Axl function as a potential target for future antiviral therapies.},
  language  = {en}
}
@article{QuastGescheidtSpichty2020,
  author    = {Quast, Helmut and Gescheidt, Georg and Spichty, Martin},
  title     = {Topological dynamics of a radical ion pair: Experimental and computational assessment at the relevant nanosecond timescale},
  series = {Chemistry},
  volume    = {2},
  journal   = {Chemistry},
  number    = {2},
  issn      = {2624-8549},
  doi       = {10.3390/chemistry2020014},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-285195},
  pages     = {219 -- 230},
  year      = {2020},
  abstract  = {Chemical processes mostly happen in fluid environments where reaction partners encounter via diffusion. The bimolecular encounters take place at a nanosecond time scale. The chemical environment (e.g., solvent molecules, (counter)ions) has a decisive influence on the reactivity as it determines the contact time between two molecules and affects the energetics. For understanding reactivity at an atomic level and at the appropriate dynamic time scale, it is crucial to combine matching experimental and theoretical data. Here, we have utilized all-atom molecular-dynamics simulations for accessing the key time scale (nanoseconds) using a QM/MM-Hamiltonian. Ion pairs consisting of a radical ion and its counterion are ideal systems to assess the theoretical predictions because they reflect dynamics at an appropriate time scale when studied by temperature-dependent EPR spectroscopy. We have investigated a diketone radical anion with its tetra-ethylammonium counterion. We have established a funnel-like transition path connecting two (equivalent) complexation sites. The agreement between the molecular-dynamics simulation and the experimental data presents a new paradigm for ion-ion interactions. This study exemplarily demonstrates the impact of the molecular environment on the topological states of reaction intermediates and how these states can be consistently elucidated through the combination of theory and experiment. We anticipate that our findings will contribute to the prediction of bimolecular transformations in the condensed phase with relevance to chemical synthesis, polymers, and biological activity.},
  language  = {en}
}
@phdthesis{Merget2015,
  author    = {Merget, Benjamin},
  title     = {Computational methods for assessing drug-target residence times in bacterial enoyl-ACP reductases and predicting small-molecule permeability for the \(Mycobacterium\) \(tuberculosis\) cell wall},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-127386},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2015},
  abstract  = {\textbf{Molecular Determinants of Drug-Target Residence Times of Bacterial Enoyl-ACP Reductases.} Whereas optimization processes of early drug discovery campaigns are often affinity-driven, the drug-target residence time \$t_R\$ should also be considered due to an often strong correlation with \textit{in vivo} efficacy of compounds. However, rational optimization of \$t_R\$ is not straightforward and generally hampered by the lack of structural information about the transition states of ligand association and dissociation. The enoyl-ACP reductase FabI of the fatty acid synthesis (FAS) type II is an important drug-target in antibiotic research. InhA is the FabI enzyme of \textit{Mycobacterium tuberculosis}, which is known to be inhibited by various compound classes. Slow-onset inhibition of InhA is assumed to be associated with the ordering of the most flexible protein region, the substrate binding loop (SBL). Diphenylethers are one class of InhA inhibitors that can promote such SBL ordering, resulting in long drug-target residence times. Although these inhibitors are energetically and kinetically well characterized, it is still unclear how the structural features of a ligand affect \$t_R\$. Using classical molecular dynamics (MD) simulations, recurring conformational families of InhA protein-ligand complexes were detected and structural determinants of drug-target residence time of diphenyl\-ethers with different kinetic profiles were described. This information was used to deduce guidelines for efficacy improvement of InhA inhibitors, including 5'-substitution on the diphenylether B-ring. The validity of this suggestion was then analyzed by means of MD simulations. Moreover, Steered MD (SMD) simulations were employed to analyze ligand dissociation of diphenylethers from the FabI enzyme of \textit{Staphylococcus aureus}. This approach resulted in a very accurate and quantitative linear regression model of the experimental \$ln(t_R)\$ of these inhibitors as a function of the calculated maximum free energy change of induced ligand extraction. This model can be used to predict the residence times of new potential inhibitors from crystal structures or valid docking poses. Since correct structural characterization of the intermediate enzyme-inhibitor state (EI) and the final state (EI*) of two-step slow-onset inhibition is crucial for rational residence time optimization, the current view of the EI and EI* states of InhA was revisited by means of crystal structure analysis, MD and SMD simulations. Overall, the analyses affirmed that the EI* state is a conformation resembling the 2X23 crystal structure (with slow-onset inhibitor \textbf{PT70}), whereas a twist of residues Ile202 and Val203 with a further opened helix \$\alpha 6\$ corresponds to the EI state. Furthermore, MD simulations emphasized the influence of close contacts to symmetry mates in the SBL region on SBL stability, underlined by the observation that an MD simulation of \textbf{PT155} chain A with chain B' of a symmetry mate in close proximity of the SBL region showed significantly more stable loops, than a simulation of the tetrameric assembly. Closing Part I, SMD simulations were employed which allow the delimitation of slow-onset InhA inhibitors from rapid reversible ligands. \textbf{Prediction of \textit{Mycobacterium tuberculosis} Cell Wall Permeability.} The cell wall of \textit{M. tuberculosis} hampers antimycobacterial drug design due to its unique composition, providing intrinsic antibiotic resistance against lipophilic and hydrophilic compounds. To assess the druggability space of this pathogen, a large-scale data mining endeavor was conducted, based on multivariate statistical analysis of differences in the physico-chemical composition of a normally distributed drug-like chemical space and a database of antimycobacterial--and thus very likely permeable--compounds. The approach resulted in the logistic regression model MycPermCheck, which is able to predict the permeability probability of small organic molecules based on their physico-chemical properties. Evaluation of MycPermCheck suggests a high predictive power. The model was implemented as a freely accessible online service and as a local stand-alone command-line version. Methodologies and findings from both parts of this thesis were combined to conduct a virtual screening for antimycobacterial substances. MycPermCheck was employed to screen the chemical permeability space of \textit{M. tuberculosis} from the entire ZINC12 drug-like database. After subsequent filtering steps regarding ADMET properties, InhA was chosen as an exemplary target. Docking to InhA led to a principal hit compound, which was further optimized. The quality of the interaction of selected derivatives with InhA was subsequently evaluated using MD and SMD simulations in terms of protein and ligand stability, as well as maximum free energy change of induced ligand egress. The results of the presented computational experiments suggest that compounds with an indole-3-acethydrazide scaffold might constitute a novel class of InhA inhibitors, worthwhile of further investigation.},
  subject      = {Computational chemistry},
  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}
}
@article{CapraBusnelliPerennaetal.2013,
  author    = {Capra, Val{\´e}rie and Busnelli, Marta and Perenna, Alessandro and Ambrosio, Manuela and Accomazzo, Maria Rosa and Gal{\´e}s, Celine and Chini, Bice and Rovati, G. Enrico},
  title     = {Full and Partial Agonists of Thromboxane Prostanoid Receptor Unveil Fine Tuning of Receptor Superactive Conformation and G Protein Activation},
  series = {PLoS ONE},
  volume    = {8},
  journal   = {PLoS ONE},
  number    = {3},
  doi       = {10.1371/journal.pone.0060475},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-131013},
  pages     = {e60475},
  year      = {2013},
  abstract  = {The intrahelical salt bridge between \(E/D^{3.49}\) and \(R^{3.50}\) within the E/DRY motif on helix 3 (H3) and the interhelical hydrogen bonding between the E/DRY and residues on H6 are thought to be critical in stabilizing the class A G protein-coupled receptors in their inactive state. Removal of these interactions is expected to generate constitutively active receptors. This study examines how neutralization of \(E^{3.49/6.30}\) in the thromboxane prostanoid (TP) receptor alters ligand binding, basal, and agonist-induced activity and investigates the molecular mechanisms of G protein activation. We demonstrate here that a panel of full and partial agonists showed an increase in affinity and potency for E129V and E240V mutants. Yet, even augmenting the sensitivity to detect constitutive activity (CA) with overexpression of the receptor or the G protein revealed resistance to an increase in basal activity, while retaining fully the ability to cause agonist-induced signaling. However, direct G protein activation measured through bioluminescence resonance energy transfer (BRET) indicates that these mutants more efficiently communicate and/or activate their cognate G proteins. These results suggest the existence of additional constrains governing the shift of TP receptor to its active state, together with an increase propensity of these mutants to agonist-induced signaling, corroborating their definition as superactive mutants. The particular nature of the TP receptor as somehow "resistant" to CA should be examined in the context of its pathophysiological role in the cardiovascular system. Evolutionary forces may have favored regulation mechanisms leading to low basal activity and selected against more highly active phenotypes.},
  language  = {en}
}