@misc{Helfmann2018, type = {Master Thesis}, author = {Helfmann, Luzie}, title = {Stochastic Modeling of Interacting Agent Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-71887}, pages = {68}, year = {2018}, language = {en} } @phdthesis{Reuter2018, author = {Reuter, Bernhard}, title = {Generalisierte Markov-Modellierung von Nichtgleichgewichtssystemen - Simulation und Modellierung der Amyloid-beta(1-40)-Konformationsdynamik unter Mikrowelleneinfluss}, year = {2018}, language = {de} } @phdthesis{Quer2018, author = {Quer, Jannes}, title = {Importance Sampling for metastable dynamical systems in molecular dynamics}, year = {2018}, language = {en} } @phdthesis{Klimm2018, author = {Klimm, Martina}, title = {New Strategies in Conformation Dynamics}, year = {2018}, language = {en} } @phdthesis{Villatoro2018, author = {Villatoro, Jos{\´e}}, title = {A combined approach for the analysis of biomolecules using IR-MALDI ion mobility spectrometry and molecular dynamics simulations of peptide ions in the gas phase}, year = {2018}, language = {en} } @article{ReidelbachWeberImhof2018, author = {Reidelbach, Marco and Weber, Marcus and Imhof, Petra}, title = {Prediction of perturbed proton transfer networks}, volume = {13}, journal = {PLoS ONE}, number = {12}, doi = {https://doi.org/10.1371/journal.pone.0207718}, pages = {e0207718 -- e0207718}, year = {2018}, abstract = {The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80\%, while important network properties are maintained in most predictions.}, language = {en} } @article{QuerDonatiKelleretal.2018, author = {Quer, Jannes and Donati, Luca and Keller, Bettina and Weber, Marcus}, title = {An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates}, volume = {40}, journal = {SIAM Journal on Scientific Computing}, number = {2}, doi = {10.1137/17m1124772}, pages = {A653 -- A670}, year = {2018}, abstract = {In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities.}, language = {en} } @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} } @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} } @article{WitzigBeckenbachEifleretal.2018, author = {Witzig, Jakob and Beckenbach, Isabel and Eifler, Leon and Fackeldey, Konstantin and Gleixner, Ambros and Grever, Andreas and Weber, Marcus}, title = {Mixed-Integer Programming for Cycle Detection in Non-reversible Markov Processes}, volume = {16}, journal = {Multiscale Modeling and Simulation}, number = {1}, issn = {1438-0064}, doi = {10.1137/16M1091162}, pages = {248 -- 265}, year = {2018}, abstract = {In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.}, language = {en} }