@article{Weber2018, author = {Weber, Marcus}, title = {Transformationsprodukte im Kl{\"a}rwerk: Mathematische Ans{\"a}tze der Bewertung}, journal = {KA Korrespondenz Abwasser, Abfall}, year = {2018}, language = {de} } @misc{Breitbach2018, type = {Master Thesis}, author = {Breitbach, Paula}, title = {Molecular Binding Kinetics of CYP P450 by Using the Infinitesimal Generator Approach}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-69076}, year = {2018}, language = {en} } @article{ReuterWeberFackeldeyetal.2018, author = {Reuter, Bernhard and Weber, Marcus and Fackeldey, Konstantin and R{\"o}blitz, Susanna and Garcia, Martin E.}, title = {Generalized Markov State Modeling Method for Nonequilibrium Biomolecular Dynamics: Exemplified on Amyloid β Conformational Dynamics Driven by an Oscillating Electric Field}, volume = {14}, journal = {Journal of Chemical Theory and Computation}, number = {7}, doi = {10.1021/acs.jctc.8b00079}, pages = {3579 -- 3594}, year = {2018}, abstract = {Markov state models (MSMs) have received an unabated increase in popularity in recent years, as they are very well suited for the identification and analysis of metastable states and related kinetics. However, the state-of-the-art Markov state modeling methods and tools enforce the fulfillment of a detailed balance condition, restricting their applicability to equilibrium MSMs. To date, they are unsuitable to deal with general dominant data structures including cyclic processes, which are essentially associated with nonequilibrium systems. To overcome this limitation, we developed a generalization of the common robust Perron Cluster Cluster Analysis (PCCA+) method, termed generalized PCCA (G-PCCA). This method handles equilibrium and nonequilibrium simulation data, utilizing Schur vectors instead of eigenvectors. G-PCCA is not limited to the detection of metastable states but enables the identification of dominant structures in a general sense, unraveling cyclic processes. This is exemplified by application of G-PCCA on nonequilibrium molecular dynamics data of the Amyloid β (1-40) peptide, periodically driven by an oscillating electric field.}, language = {en} } @article{WagnerZapataWanetal.2018, author = {Wagner, Sabine and Zapata, Carlos and Wan, Wei and Gawlitza, Kornelia and Weber, Marcus and Rurack, Knut}, title = {Role of Counterions in Molecularly Imprinted Polymers for Anionic Species}, volume = {34}, journal = {Langmuir}, number = {23}, doi = {10.1021/acs.langmuir.8b00500}, pages = {6963 -- 6975}, year = {2018}, abstract = {Small-molecule oxoanions are often imprinted noncovalently as carboxylates into molecularly imprinted polymers (MIPs), requiring the use of an organic counterion. Popular species are either pentamethylpiperidine (PMP) as a protonatable cation or tetraalkylammonium (TXA) ions as permanent cations. The present work explores the influence of the TXA as a function of their alkyl chain length, from methyl to octyl, using UV/vis absorption, fluorescence titrations, and HPLC as well as MD simulations. Protected phenylalanines (Z-L/D-Phe) served as templates/analytes. While the influence of the counterion on the complex stability constants and anion-induced spectral changes shows a monotonous trend with increasing alkyl chain length at the prepolymerization stage, the cross-imprinting/rebinding studies showed a unique pattern that suggested the presence of adaptive cavities in the MIP matrix, related to the concept of induced fit of enzyme-substrate interaction. Larger cavities formed in the presence of larger counterions can take up pairs of Z-x-Phe and smaller TXA, eventually escaping spectroscopic detection. Correlation of the experimental data with the MD simulations revealed that counterion mobility, the relative distances between the three partners, and the hydrogen bond lifetimes are more decisive for the response features observed than actual distances between interacting atoms in a complex or the orientation of binding moieties. TBA has been found to yield the highest imprinting factor, also showing a unique dual behavior regarding the interaction with template and fluorescent monomer. Finally, interesting differences between both enantiomers have been observed in both theory and experiment, suggesting true control of enantioselectivity. The contribution concludes with suggestions for translating the findings into actual MIP development.}, language = {en} } @article{SpahnDelVecchioRodriguezGaztelumendietal.2018, author = {Spahn, Viola and Del Vecchio, Giovanna and Rodriguez-Gaztelumendi, Antonio and Temp, Julia and Labuz, Dominika and Kloner, Michael and Reidelbach, Marco and Machelska, Halina and Weber, Marcus and Stein, Christoph}, title = {Opioid receptor signaling, analgesic and side effects induced by a computationally designed pH-dependent agonist}, volume = {8}, journal = {Scientific Reports}, publisher = {Springer Nature}, pages = {8965}, year = {2018}, abstract = {Novel pain killers without adverse effects are urgently needed.}, language = {en} } @article{ErlekamIgdeRoeblitzetal.2019, author = {Erlekam, Franziska and Igde, Sinaida and R{\"o}blitz, Susanna and Hartmann, Laura and Weber, Marcus}, title = {Modeling of Multivalent Ligand-Receptor Binding Measured by kinITC}, volume = {7}, journal = {Computation}, number = {3}, doi = {10.3390/computation7030046}, pages = {46}, year = {2019}, abstract = {In addition to the conventional Isothermal Titration Calorimetry (ITC), kinetic ITC (kinITC) not only gains thermodynamic information, but also kinetic data from a biochemical binding process. Moreover, kinITC gives insights into reactions consisting of two separate kinetic steps, such as protein folding or sequential binding processes. The ITC method alone cannot deliver kinetic parameters, especially not for multivalent bindings. This paper describes how to solve the problem using kinITC and an invariant subspace projection. The algorithm is tested for multivalent systems with different valencies.}, language = {en} } @article{Weber2019, author = {Weber, Marcus}, title = {Transformationsprodukte im Kl{\"a}rwerk: Mathematische Ans{\"a}tze der Bewertung}, volume = {7}, journal = {KA Korrespondenz Abwasser, Abfall}, pages = {551 -- 557}, year = {2019}, language = {de} } @article{VenkatareddyWilkeErnstetal.2019, author = {Venkatareddy, Narendra Lagumaddepalli and Wilke, Patrick and Ernst, Natalia and Horch, Justus and Weber, Marcus and Dallmann, Andre and B{\"o}rner, Hans G.}, title = {Mussel-glue inspired adhesives: A study on the relevance of L-Dopa and the function of the sequence at nanomaterial-peptide interfaces}, volume = {6}, journal = {Advanced Materials Interfaces}, number = {13}, doi = {10.1002/admi.201900501}, pages = {1900501}, year = {2019}, abstract = {Mussel glue-proteins undergo structural transitions at material interfaces to optimize adhesive surface contacts. Those intriguing structure responses are mimicked by a mussel-glue mimetic peptide (HSY*SGWSPY*RSG (Y* = l-Dopa)) that was previously selected by phage-display to adhere to Al2O3 after enzymatic activation. Molecular level insights into the full-length adhesion domain at Al2O3 surfaces are provided by a divergent-convergent analysis, combining nuclear Overhauser enhancement based 2D NOESY and saturation transfer difference NMR analysis of submotifs along with molecular dynamics simulations of the full-length peptide. The peptide is divided into two submotifs, each containing one Dopa "anchor" (Motif-1 and 2). The analysis proves Motif-1 to constitute a dynamic Al2O3 binder and adopting an "M"-structure with multiple surface contacts. Motif-2 binds stronger by two surface contacts, forming a compact "C"-structure. Taking these datasets as constraints enables to predict the structure and propose a binding process model of the full-length peptide adhering to Al2O3.}, language = {en} } @article{DelVecchioLabuzTempetal.2019, author = {Del Vecchio, Giovanna and Labuz, Dominika and Temp, Julia and Seitz, Viola and Kloner, Michael and Negrete, Roger and Rodriguez-Gaztelumendi, Antonio and Weber, Marcus and Machelska, Halina and Stein, Christoph}, title = {pKa of opioid ligands as a discriminating factor for side effects}, volume = {9}, journal = {Scientific Reports}, doi = {10.1038/s41598-019-55886-1}, pages = {19344}, year = {2019}, abstract = {The non-selective activation of central and peripheral opioid receptors is a major shortcoming of currently available opioids. Targeting peripheral opioid receptors is a promising strategy to preclude side effects. Recently, we showed that fentanyl-derived μ-opioid receptor (MOR) agonists with reduced acid dissociation constants (pKa) due to introducing single fluorine atoms produced injury-restricted antinociception in rat models of inflammatory, postoperative and neuropathic pain. Here, we report that a new double-fluorinated compound (FF6) and fentanyl show similar pKa, MOR affinity and [35S]-GTPγS binding at low and physiological pH values. In vivo, FF6 produced antinociception in injured and non-injured tissue, and induced sedation and constipation. The comparison of several fentanyl derivatives revealed a correlation between pKa values and pH-dependent MOR activation, antinociception and side effects. An opioid ligand's pKa value may be used as discriminating factor to design safer analgesics.}, language = {en} } @misc{Bender2019, type = {Master Thesis}, author = {Bender, Wilhelm}, title = {On a Novel Approach for Global Optimization of Non-Convex Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-77074}, year = {2019}, abstract = {In this work new approaches for global optimization are developed. We move from theoretical connections in complexity theory to very practical results in protein-ligand docking. In regard to the efficiency of the suggested optimizer, continuity methods are motivated, that yield an improved optimization method.}, language = {en} } @book{Reuter2020, author = {Reuter, Bernhard}, title = {Generalisierte Markov-Modellierung}, publisher = {Springer Spektrum, Wiesbaden}, isbn = {978-3-658-29711-4}, doi = {10.1007/978-3-658-29712-1}, year = {2020}, abstract = {Markov State Models (MSM) sind der Goldstandard zur Modellierung biomolekularer Dynamik, da sie die Identifizierung und Analyse metastabiler Zust{\"a}nde erm{\"o}glichen. Die robuste Perron-Cluster-Cluster-Analyse (PCCA+) ist ein verbreiteter Spectral-Clustering-Algorithmus, der f{\"u}r das Clustering hochdimensionaler MSM verwendet wird. Da die PCCA+ auf reversible Prozesse beschr{\"a}nkt ist, wird sie zur Generalisierten PCCA+ (G-PCCA) verallgemeinert, die geeignet ist, nichtreversible Prozesse aufzukl{\"a}ren. Bernhard Reuter untersucht hier mittels G-PCCA die nichtthermischen Auswirkungen von Mikrowellen auf die Proteindynamik. Dazu f{\"u}hrt er molekulardynamische Nichtgleichgewichtssimulationen des Amyloid-β-(1-40)-Peptids durch und modelliert diese.}, language = {de} } @misc{Wulkow2017, type = {Master Thesis}, author = {Wulkow, Niklas}, title = {Modelling the Spread of Innovations by a Markov Process in a Bayesian Framework}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66869}, year = {2017}, language = {en} } @phdthesis{Durmaz2016, author = {Durmaz, Vedat}, title = {Atomistic Binding Free Energy Estimations for Biological Host-Guest Systems}, publisher = {FU Dissertationen Online}, pages = {216}, year = {2016}, abstract = {Accurate quantifications of protein-ligand binding affinities by means of in silico methods increasingly gain importance in various scientific branches including toxicology and pharmacology. In silico techniques not only are generally less demanding than laboratory experiments regarding time as well as cost, in particular, if binding assays or synthesis protocols need to be developed in advance. At times, they also provide the only access to risk assessments on novel chemical compounds arising from biotic or abiotic degradation of anthropogenic substances. However, despite the continuous technological and algorithmic progress over the past decades, binding free energy estimations through molecular dynamics simulations still pose an enormous computational challenge owed to the mathematical complexity of solvated macromolecular systems often consisting of hundreds of thousands of atoms. The goals of this thesis can roughly be divided into two categories dealing with different aspects of host-guest binding quantification. On the one side algorithmic strategies for a comprehensive exploration and decomposition of conformational space in conjunction with an automated selection of representative molecular geometries and binding poses have been elaborated providing initial structures for free energy calculations. In light of the dreaded trapping problem typically associated with molecular dynamics simulations, the focus was laid on a particularly systematic generation of representatives covering a broad range of physically accessible molecular conformations and interaction modes. On the other side and ensuing from these input geometries, binding affinity models based on the linear interaction energy (LIE) method have been developed for a couple of (bio)molecular systems. The applications included a successful prediction of the liquid-chromatographic elution order as well as retention times of highly similar hexabromocyclododecane (HBCD) stereoisomers, a novel empirical LIE-QSAR hybrid binding affinity model related to the human estrogen receptor α (ERα), and, finally, the (eco)toxicological prioritization of transformation products originating from the antibiotic sulfamethoxazole with respect to their binding affinities to the bacterial enzyme dihydropteroate synthase. Altogether, a fully automated approach to binding mode and affinity estimation has been presented that is content with an arbitrary geometry of a small molecule under observation and a spatial vector specifying the binding site of a potential target molecule. According to our studies, it is superior to conventional docking and thermodynamic average methods and primarily suggesting binding free energy calculation on the basis of several heavily distinct complex geometries. Both chromatographic retention times of HBCD and binding affinities to ERα yielded squared coefficients of correlation with experimental results significantly higher than 0.8. Approximately 85 \% (100 \%) of predicted receptor-ligand binding modes deviated less than 1.53 {\AA} (2.05 {\AA}) from available crystallographic structures.}, language = {en} } @misc{WeberDurmazSabrietal.2017, author = {Weber, Marcus and Durmaz, Vedat and Sabri, Peggy and Reidelbach, Marco}, title = {Supplementary simulation data for Science Manuscript ai8636}, doi = {10.12752/5.MWB.1.0}, year = {2017}, abstract = {The simulation data has been produced by Vedat Durmaz, Peggy Sabri and Marco Reidelbach inside the "Computational Molecular Design" Group headed by Marcus Weber at Zuse-Institut Berlin, Takustr. 7, D-14195 Berlin, Germany. The file contains classical simulation data for different fentanyl derivates in the MOR binding pocket at different pHs. It also includes instruction files for quantum-chemical pKa-value estimations and a description of how we derived the pKa-values from the Gaussian09 log-files.}, language = {en} } @misc{Weber2018, author = {Weber, Marcus}, title = {Supplementary: Implications of PCCA+ in Molecular Simulation}, year = {2018}, abstract = {Matlab-software and data sets to recapitulate the presented results in M. Weber: Implications of PCCA+ in Molecular Simulation. Computation, 6(1):20, 2018.}, language = {en} } @article{Weber2018, author = {Weber, Marcus}, title = {Implications of PCCA+ in Molecular Simulation}, volume = {6}, journal = {Computation}, number = {1}, doi = {10.3390/computation6010020}, pages = {20}, year = {2018}, abstract = {Upon ligand binding or during chemical reactions the state of a molecular system changes in time. Usually we consider a finite set of (macro-) states of the system (e.g., 'bound' vs. 'unbound'), although the process itself takes place in a continuous space. In this context, the formula chi=XA connects the micro-dynamics of the molecular system to its macro-dynamics. Chi can be understood as a clustering of micro-states of a molecular system into a few macro-states. X is a basis of an invariant subspace of a transfer operator describing the micro-dynamics of the system. The formula claims that there is an unknown linear relation A between these two objects. With the aid of this formula we can understand rebinding effects, the electron flux in pericyclic reactions, and systematic changes of binding rates in kinetic ITC experiments. We can also analyze sequential spectroscopy experiments and rare event systems more easily. This article provides an explanation of the formula and an overview of some of its consequences.}, language = {en} } @article{SchradeTroegerEldashanetal.2018, author = {Schrade, Katharina and Tr{\"o}ger, Jessica and Eldashan, Adeep and Z{\"u}hlke, Kerstin and Abdul Azees, Kamal R. and Elkins, Jonathan M. and Neuenschwander, Martin and Oder, Andreas and Elkewedi, Mohamed and Jaksch, Sarah and Andrae, Karsten and Li, Jinliang and Fernandes, Jaoa and M{\"u}ller, Paul Markus and Grunwald, Stephan and Marino, Stephen F. and Vukicevic, Tanja and Eichhorst, Jenny and Wiesner, Burkhard and Weber, Marcus and Kapiloff, Michael and Rocks, Oliver and Daumke, Oliver and Wieland, Thomas and Knapp, Stefan and von Kries, Jens Peter and Klussmann, Enno}, title = {An AKAP-Lbc-RhoA interaction inhibitor promotes the translocation of aquaporin-2 to the plasma membrane of renal collecting duct principal cells}, volume = {13}, journal = {PLOS ONE}, number = {1}, doi = {10.1371/journal.pone.0191423}, pages = {e0191423 -- e0191423}, year = {2018}, abstract = {Stimulation of renal collecting duct principal cells with antidiuretic hormone (arginine-vasopressin, AVP) results in inhibition of the small GTPase RhoA and the enrichment of the water channel aquaporin-2 (AQP2) in the plasma membrane. The membrane insertion facilitates water reabsorption from primary urine and fine-tuning of body water homeostasis. Rho guanine nucleotide exchange factors (GEFs) interact with RhoA, catalyze the exchange of GDP for GTP and thereby activate the GTPase. However, GEFs involved in the control of AQP2 in renal principal cells are unknown. The A-kinase anchoring protein, AKAP-Lbc, possesses GEF activity, specifically activates RhoA, and is expressed in primary renal inner medullary collecting duct principal (IMCD) cells. Through screening of 18,431 small molecules and synthesis of a focused library around one of the hits, we identified an inhibitor of the interaction of AKAP-Lbc and RhoA. This molecule, Scaff10-8, bound to RhoA, inhibited the AKAP-Lbc-mediated RhoA activation but did not interfere with RhoA activation through other GEFs or activities of other members of the Rho family of small GTPases, Rac1 and Cdc42. Scaff10-8 promoted the redistribution of AQP2 from intracellular vesicles to the periphery of IMCD cells. Thus, our data demonstrate an involvement of AKAP-Lbc-mediated RhoA activation in the control of AQP2 trafficking.}, language = {en} } @misc{HelfmannDjurdjevacConradDjurdjevacetal.2019, author = {Helfmann, Luzie and Djurdjevac Conrad, Natasa and Djurdjevac, Ana and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {From interacting agents to density-based modeling with stochastic PDEs}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-73456}, year = {2019}, abstract = {Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models.}, language = {en} } @article{HelfmannDjurdjevacConradDjurdjevacetal.2021, author = {Helfmann, Luzie and Djurdjevac Conrad, Natasa and Djurdjevac, Ana and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {From interacting agents to density-based modeling with stochastic PDEs}, volume = {16}, journal = {Communications in Applied Mathematics and Computational Science}, number = {1}, arxiv = {http://arxiv.org/abs/1905.13525}, doi = {10.2140/camcos.2021.16.1}, pages = {1 -- 32}, year = {2021}, abstract = {Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models.}, language = {en} } @article{ChewleEmmerlingWeber2020, author = {Chewle, Surahit and Emmerling, Franziska and Weber, Marcus}, title = {Effect of choice of solvent on crystallization pathway of Paracetamol: An experimental and theoretical case study}, volume = {10}, journal = {Crystals}, number = {12}, doi = {10.3390/cryst10121107}, pages = {1107}, year = {2020}, abstract = {The choice of solvents influences crystalline solid formed during the crystallization of active pharmaceutical ingredients (API). The underlying effects are not always well understood because of the complexity of the systems. Theoretical models are often insufficient to describe this phenomenon. In this study, the crystallization behavior of the model drug paracetamol in different solvents was studied based on experimental and molecular dynamics data. The crystallization process was followed in situ using time-resolved Raman spectroscopy. Molecular dynamics with simulated annealing algorithm was used for an atomistic understanding of the underlying processes. The experimental and theoretical data indicate that paracetamol molecules adopt a particular geometry in a given solvent predefining the crystallization of certain polymorphs.}, language = {en} } @article{AbendrothBujotzekShanetal.2011, author = {Abendroth, Frank and Bujotzek, Alexander and Shan, Min and Haag, Rainer and Weber, Marcus and Seitz, Oliver}, title = {DNA-controlled bivalent presentation of ligands for the estrogen receptor}, journal = {Angew. Chem. Int. Ed.}, year = {2011}, language = {en} } @article{AndraeDinhReinhardetal.2006, author = {Andrae, Karsten and Dinh, P. and Reinhard, P.-G. and Suraud, E.}, title = {Pump and probe analysis of metal cluster dynamics}, volume = {35}, journal = {Computational Materials Science}, number = {3}, doi = {DOI: 10.1016/j.commatsci.2004.09.062}, pages = {169 -- 173}, year = {2006}, language = {en} } @article{AndraeReinhardSuraud2004, author = {Andrae, Karsten and Reinhard, P.-G. and Suraud, E.}, title = {Crossed Beam Pump and Probe Dynamics in Metal Clusters}, volume = {92}, journal = {Phys. Rev. Lett.}, number = {17}, publisher = {American Physical Society}, doi = {10.1103/PhysRevLett.92.173402}, pages = {173402}, year = {2004}, language = {en} } @article{AndraeReinhardSuraud2007, author = {Andrae, Karsten and Reinhard, P.-G. and Suraud, E.}, title = {Systematics of spin-polarized small Na clusters}, journal = {Annals of the European Academy of Science}, arxiv = {http://arxiv.org/abs/arXiv:physics/0701028}, year = {2007}, language = {en} } @article{BujotzekShanHaagetal.2011, author = {Bujotzek, Alexander and Shan, Min and Haag, Rainer and Weber, Marcus}, title = {Towards a rational spacer design for bivalent inhibition of estrogen receptor}, volume = {25(3)}, journal = {J. Comput.-Aided Mol. Des.}, pages = {253 -- 262}, year = {2011}, language = {en} } @incollection{DeuflhardWeber2005, author = {Deuflhard, Peter and Weber, Marcus}, title = {Robust Perron Cluster Analysis in Conformation Dynamics}, volume = {398}, booktitle = {Lin. Alg. Appl. - Special Issue on Matrices and Mathematical Biology}, publisher = {Elsevier Journals}, address = {Germany}, pages = {161 -- 184}, year = {2005}, language = {en} } @article{KoschekDurmazKrylovaetal.2015, author = {Koschek, and Durmaz, Vedat and Krylova, and Wieczorek, and Gupta, Pooja and Richter, and Bujotzek, Alexander and Fischer, and Haag, Rainer and Freund, and Weber, Marcus and Rademann,}, title = {Peptide polymer ligands for a tandem WW-domain, a soft multivalent protein-protein interaction: lessons on the thermodynamic fitness of flexible ligands}, volume = {11}, journal = {Beilstein J. Org. Chem.}, pages = {837 -- 847}, year = {2015}, language = {en} } @article{DurmazWeberMeyeretal.2015, author = {Durmaz, Vedat and Weber, Marcus and Meyer, and M{\"u}ckter,}, title = {Computergest{\"u}tzte Simulationen zur Absch{\"a}tzung gesundheitlicher Risiken durch anthropogene Spurenstoffe der Wassermatrix}, volume = {3/15}, journal = {KA Korrespondenz Abwasser, Abfall}, pages = {264 -- 267}, year = {2015}, language = {de} } @misc{WeberFackeldey2015, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {G-PCCA: Spectral Clustering for Non-reversible Markov Chains}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55505}, year = {2015}, abstract = {Spectral clustering methods are based on solving eigenvalue problems for the identification of clusters, e.g., the identification of metastable subsets of a Markov chain. Usually, real-valued eigenvectors are mandatory for this type of algorithms. The Perron Cluster Analysis (PCCA+) is a well-known spectral clustering method of Markov chains. It is applicable for reversible Markov chains, because reversibility implies a real-valued spectrum. We extend this spectral clustering method also to non-reversible Markov chains and give some illustrative examples. The main idea is to replace the eigenvalue problem by a real-valued Schur decomposition. By this extension, non-reversible Markov chains can be analyzed. Furthermore, the chains need not have a positive stationary distribution. And additionally to metastabilities, dominant cycles and sinks can be identified, too.}, language = {en} } @article{AbendrothSollederWelkeretal.2015, author = {Abendroth, Frank and Solleder, Marthe and Welker, Pia and Licha, Kai and Weber, Marcus and Seitz, Oliver and Mangoldt, Dorothea}, title = {High affinity flourescence labelled ligands for the estrogen receptor}, volume = {2015}, journal = {Eur. J. Org. Chem.}, number = {10}, pages = {2157 -- 2166}, year = {2015}, language = {en} } @article{WeberZoschkeSedighietal.2014, author = {Weber, Marcus and Zoschke, Christian and Sedighi, Amir and Fleige, Emanuel and Haag, Rainer and Sch{\"a}fer-Korting, Monika}, title = {Free Energy Simulations of Drug loading for Core-Multishell Nanotransporters}, volume = {5}, journal = {J Nanomed Nanotechnol}, number = {5}, doi = {10.4172/2157-7439.1000234}, pages = {234}, year = {2014}, language = {en} } @misc{SteinWeberZoellneretal.2013, author = {Stein, Christoph and Weber, Marcus and Z{\"o}llner, Christian and Scharkoi, Olga}, title = {Fentanyl derivatives as pH-dependent opioid receptor agonists}, journal = {European Patent Application, Bulletin 2013/08}, year = {2013}, language = {en} } @misc{SteinWeberScharkoietal.2013, author = {Stein, Christoph and Weber, Marcus and Scharkoi, Olga and Deuflhard, Peter}, title = {Method and system for identifying compounds that bind and preferably activate a target opioid receptor in a pH-dependent manner}, journal = {European Patent Application, Bulletin 2013/28}, year = {2013}, language = {en} } @article{ZhangWangHartmannetal.2014, author = {Zhang, Wei and Wang, Han and Hartmann, Carsten and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Applications of the cross-entropy method to importance sampling and optimal control of diffusions}, volume = {36}, journal = {Siam Journal on Scientific Computing}, number = {6}, doi = {10.1137/14096493X}, pages = {A2654 -- A2672}, year = {2014}, language = {en} } @article{AndraeMerkelDurmazetal.2014, author = {Andrae, Karsten and Merkel, Stefan and Durmaz, Vedat and Fackeldey, Konstantin and K{\"o}ppen, Robert and Weber, Marcus and Koch, Matthias}, title = {Investigation of the Ergopeptide Epimerization Process}, volume = {2}, journal = {Computation}, number = {3}, doi = {10.3390/computation2030102}, pages = {102 -- 111}, year = {2014}, abstract = {Ergopeptides, like ergocornine and a-ergocryptine, exist in an S- and in an R-configuration. Kinetic experiments imply that certain configurations are preferred depending on the solvent. The experimental methods are explained in this article. Furthermore, computational methods are used to understand this configurational preference. Standard quantum chemical methods can predict the favored configurations by using minimum energy calculations on the potential energy landscape. However, the explicit role of the solvent is not revealed by this type of methods. In order to better understand its influence, classical mechanical molecular simulations are applied. It appears from our research that "folding" the ergopeptide molecules into an intermediate state (between the S- and the R-configuration) is mechanically hindered for the preferred configurations.}, language = {en} } @misc{DeuflhardWeber2014, author = {Deuflhard, Peter and Weber, Marcus}, title = {Mathematics without pain}, volume = {1}, journal = {MATHEON-Mathematics for Key Technologies}, editor = {Deuflhard, Peter and Gr{\"o}tschel, Martin and H{\"o}mberg, Dietmar and Horst, Ulrich and Kramer, J{\"u}rg and Mehrmann, Volker and Polthier, Konrad and Schmidt, Frank and Skutella, Martin and Sprekels, J{\"u}rgen}, publisher = {European Mathematical Society}, doi = {10.4171/137}, pages = {26 -- 28}, year = {2014}, language = {en} } @article{SchuetteNielsenWeber2015, author = {Sch{\"u}tte, Christof and Nielsen, Adam and Weber, Marcus}, title = {Markov State Models and Molecular Alchemy}, volume = {113}, journal = {Molecular Physics}, number = {1}, doi = {10.1080/00268976.2014.944597}, pages = {69 -- 78}, year = {2015}, abstract = {In recent years Markov State Models (MSMs) have attracted a consid- erable amount of attention with regard to modelling conformation changes and associated function of biomolecular systems. They have been used successfully, e.g., for peptides including time-resolved spectroscopic experiments, protein function and protein folding , DNA and RNA, and ligand-receptor interaction in drug design and more complicated multivalent scenarios. In this article a novel reweighting scheme is introduced that allows to construct an MSM for certain molecular system out of an MSM for a similar system. This permits studying how molecular properties on long timescales differ between similar molecular systems without performing full molecular dynamics simulations for each system under con- sideration. The performance of the reweighting scheme is illustrated for simple test cases including one where the main wells of the respective energy landscapes are located differently and an alchemical transformation of butane to pentane where the dimension of the state space is changed.}, language = {en} } @article{WeberFackeldey2014, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the Minimal Rebinding Effect Included in a Given Kinetics}, volume = {12}, journal = {Multiscale Model. Simul.}, number = {1}, doi = {10.1137/13091124X}, pages = {318 -- 334}, year = {2014}, abstract = {The rebinding effect is a phenomenon which occurs when observing a ligand-receptor binding process. On the macro scale this process comprises the Markov property. This Makovian view is spoiled when switching to the atomistic scale of a binding process. We therefore suggest a model which accurately describes the rebinding effect on the atomistic scale by allowing ''intermediate'' bound states. This allows us to define an indicator for the magnitude of rebinding and to formulate an optimization problem. The results form our examples show good agreement with data form laboratory.}, language = {en} } @masterthesis{Solleder2014, type = {Bachelor Thesis}, author = {Solleder, Marthe}, title = {Bewertung von Transformationsprodukten im Wasserkreislauf durch Computersimulationen}, year = {2014}, language = {de} } @masterthesis{Schmiedel2013, type = {Bachelor Thesis}, author = {Schmiedel, Jessica}, title = {Singul{\"a}rwertzerlegung mit Zufallsalgorithmen}, year = {2013}, language = {de} } @masterthesis{Toudic2013, type = {Bachelor Thesis}, author = {Toudic, Remi}, title = {Confjump: A method for biomolecular sampling}, year = {2013}, language = {de} } @misc{Pommer2013, type = {Master Thesis}, author = {Pommer, Daniel}, title = {Konditionsverbesserung von Markov-Modellen}, year = {2013}, language = {de} } @misc{Riesland2012, type = {Master Thesis}, author = {Riesland, Mark Daniel}, title = {Anwendbarkeit von Thermostaten in der molek{\"u}ldynamischen Simulation}, year = {2012}, language = {de} } @masterthesis{Jeschke2012, type = {Bachelor Thesis}, author = {Jeschke, Bj{\"o}rn-Marcel}, title = {Spezielle Eigenschaften von Simulated Annealing}, year = {2012}, language = {de} } @misc{Josten2012, type = {Master Thesis}, author = {Josten, Jan-Martin}, title = {{\"O}strogenit{\"a}ts-Approximation mittels physikochemischem Oberfl{\"a}chenvergleich}, year = {2012}, language = {de} } @misc{Lie2012, type = {Master Thesis}, author = {Lie, Han Cheng}, title = {A Meshless Method for Conformational Analysis}, year = {2012}, language = {en} } @misc{Haack2009, type = {Master Thesis}, author = {Haack, Fiete}, title = {Representative Spectral Clustering for Large Data Sets applied to Gene Expression Data}, year = {2009}, language = {en} } @misc{NielsenWeber2014, author = {Nielsen, Adam and Weber, Marcus}, title = {Computing the nearest reversible Markov chain}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53292}, year = {2014}, abstract = {Reversible Markov chains are the basis of many applications. However, computing transition probabilities by a finite sampling of a Markov chain can lead to truncation errors. Even if the original Markov chain is reversible, the approximated Markov chain might be non-reversible and will lose important properties, like the real valued spectrum. In this paper, we show how to find the closest reversible Markov chain to a given transition matrix. It turns out that this matrix can be computed by solving a convex minimization problem.}, language = {en} } @article{NielsenWeber2015, author = {Nielsen, Adam and Weber, Marcus}, title = {Computing the nearest reversible Markov chain}, volume = {22}, journal = {Numerical Linear Algebra with Applications}, number = {3}, doi = {10.1002/nla.1967}, pages = {483 -- 499}, year = {2015}, abstract = {Reversible Markov chains are the basis of many applications. However, computing transition probabilities by a finite sampling of a Markov chain can lead to truncation errors. Even if the original Markov chain is reversible, the approximated Markov chain might be non-reversible and will lose important properties, like the real valued spectrum. In this paper, we show how to find the closest reversible Markov chain to a given transition matrix. It turns out that this matrix can be computed by solving a convex minimization problem.}, language = {en} } @article{WeberFackeldeySchuette2017, author = {Weber, Marcus and Fackeldey, Konstantin and Sch{\"u}tte, Christof}, title = {Set-Free Markov State Model Building}, volume = {146}, journal = {Journal of Chemical Physics}, number = {12}, doi = {10.1063/1.4978501}, year = {2017}, language = {en} } @article{QuerLie2017, author = {Quer, Jannes and Lie, Han Cheng}, title = {Some connections between importance sampling and enhanced sampling methods in molecular dynamics}, journal = {Journal of Chemical Physics}, year = {2017}, abstract = {Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Sch{\"u}tte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling.}, language = {en} }