Exploring theoretically ligands and polymeric linkers to make multivalent ligands

Abstract

The indirect drug design approach outlined in the first section provides an alternative way to model drugs efficiently in case structural knowledge about the target protein is absent. The procedure employs a specialized energy function to be used with molecular dynamics (MD) simulations of ligand pairs binding to same target, allowing the ligands to exchange the information they have about the target. This yields ligand conformers that are most similar to bound ligand conformers. The proposition in the second section is about successful application of an optimized ether force field for polyethylene glycol (PEG), which is used to connect ligands serving as drugs to form dimeric or even multimeric ligands and thus enhancing drug activity by the multivalent binding effect. The proposed ether force field, specifically optimized for implicit solvent (IS) simulation allows to perform molecular dynamics of ligands connected by PEG that are as realistic as explicit solvent (ES) simulation but at much reduced cost. Structure based drug design involves structural knowledge on the target protein and its binding pocket. However, in absence of that and if the structural information of a set of ligands binding to the same target is available, one can still use indirect drug-design approach proposed in the present study. However, this requires a set of ligands with chemically different architecture binding to the same target protein. There exist several indirect drug design methods. The present approach yields coordinates of the ligand that is the bound ligand conformer in binding pocket (BLC), which can then be taken further ahead where different strategies of drug design can be pursued. The assumption that a set of ligands with chemically different architecture binding to the same target protein carry hidden information of their corresponding true BLCs has been exploited in the indirect drug-design approach. Pairwise flexible structure alignment (FSA) using molecular dynamics simulations with specialized energy function employed for this procedure extract the hidden information derived from the molecular similarity of the ligands that allows the ligands to adopt the same space. The conformations of different ligand pairs obtained from FSA once combined yield consensus ligand conformers (CLCs) that should be similar to BLCs. This proposed procedure has been validated on 44 HIV-1 protease (HIV-P) bound ligands of sufficiently diverse chemical composition, the crystal structures of which are available in PDB. The described four different cluster of HIV-P BLCs based on similarity measures are consistent with the H-bond patterns of the ligands bound to HIV-P in the crystal structures exhibiting four different binding modes. MD simulations of biomolecules with realistic representations of cellular environments remain challenging. Explicitly considering water molecules during molecular modeling studies is the most realistic theoretical approach but is computationally expensive. Implicit solvent (IS) approaches, which include the effect of water in a potential of mean force have a clear advantage over explicit solvent (ES) in efficiency. For the attractive features of the IS methodology one has to pay a price in making certain approximations, which may not be sufficient to provide agreement between ES and IS approaches. In our study, an accurate polyether force field for IS MD simulations was developed that matches local and global conformations of 1,2-dimethoxy-ethane (DME) and polyethylene glycol (PEG), respectively. In the IS model, since solvent molecules are missing, the competition of attractive van der Waals (vdW) interactions between solute- solute and solute-solvent atom pairs is absent. As a consequence the remaining attractive vdW interactions between solute atom pairs yield too compact structures. To regain balance of attractive interactions in the IS model, surface energy terms are turned off and attractive vdW interactions are reduced to 90% or as an alternative even slightly negative surface energies with no change in the vdW attractive interactions are employed. However, to obtain quantitatively the same local and global distributions of PEG conformers as in ES, additional force field adjustments involving torsion potentials and 1-4 and 1-5 atom pair Coulomb interactions are also suggested. The proposed ether force field, specifically optimized for IS simulation conditions is equally valid for monomeric and polymeric ethylene glycol.

Der im ersten Absatz beschriebene indirekte Ansatz erlaubt eine effiziente Alternative für das drug design im Falle nicht vorhandener Strukturinformation. Die Methode nutzt eine spezialisierte Energiefunktion für Molekulardynamik (MD) Simulationen von Ligandenpaaren die das selbe Zielprotein binden, die es den Liganden erlaubt ihre Informationen über das Zielprotein auszutauschen. Dies führt zu Ligandenkonformationen die extreme Ähnlichkeit mit denen gebundener Liganden aufweisen. Der Vorschlag im zweiten Absatz beschreibt die erfolgreiche Anwendung optimierter Etherkraftfelder für polyethylenglykol (PEG ). PEG findet im drug design als Linker zwischen Liganden Anwendung, um sie zu Dimeren oder Multimeren zusammenzufassen und dadurch die Wirkung über multivalente Bindungen zu erhöhen. Das vorgestellte Etherkraftfeld, erlaubt MD Simulationen mit implizitem Lösemittel (implicit solvent: IS) für Liganden die über PEG verknüpft sind mit der selben Güte wie mit explizitem Lösemitteln (explicit solvent: ES) aber mit deutlich verringerter benötigter Rechnerleistung. Für strukturbasiertes design von Medikamenten (drugs) benötigt man Strukturinformationen über das Zielprotein und dessen Bindungstasche. Falls diese Informationen nicht verfügbar sind, kann man immer noch eine indirekte Methode des drug designs einsetzen. Dazu benötigt man allerdings ein Ensemble von Liganden unterschiedlicher chemischer Architektur, die alle an dasselbe Zielprotein binden. Es gibt zahlreiche Methoden des indirekten drug design. Neu an dieser Arbeit ist, dass mit der hier verwendeten Methode des drug design explizit die Koordinaten der Liganden konstruiert werden, die dessen Struktur in der Bindungstasche wiedergeben [BLC (bound ligand conformer)]. Diese Strukturinformation ermöglicht es, neue Wege beim dem design von Medikamenten zu beschreiten. Die Arbeit beruht auf der Annahme, dass ein Ensemble von Liganden, deren chemischer Aufbau sich unterscheidet, die aber an dasselbe Protein binden, versteckte Informationen ihres entsprechenden BLC in sich tragen. Die hier eingesetzte Methode des paarweisen flexiblen Strukturalignment (FSA) nutzt Molekulardynamik (MD) Simulationen von Ligandenpaaren mit für diese Aufgaben maßgeschneiderten Energiefunktionen, um die versteckten Strukturähnlichkeiten der Liganden zu extrahieren. Diese Energiefunktionen sind so konstruiert, dass die Liganden beim Strukturalignment das gleiche Volumen einnehmen können. Die Konformationen verschiedener Ligandenpaare, die über FSA erhalten werden ergeben in Kombination eine Konsensuskonformation für die Liganden (consensus ligand conformers: CLCs), die dann Ähnlichkeit mit dem BLC aufweisen. Der vorgeschlagene Ansatz ist anhand von 44 der PDB entnommenen HIV-1 Proteaseliganden (HIV-P) von ausreichend unterschiedlichem chemischem Aufbau überprüft worden. Die auf Grundlage der Strukturähnlichkeiten gefundenen vier unterschiedlichen Cluster von HIV-P BLCs sind konsistent mit den vier in den Kristallstrukturen vorgefundenen Mustern an Wasserstoffbrücken zwischen den gebundenen Formen der Liganden und der HIV-P. MD Simulation von Biomolekülen mit einer realistischen Repräsentation der zellulären Umgebung ist eine Herausforderung für computergestützte theoretische Methoden. Die Verwendung von explizit modelliertem Wasser (explicit solvent: ES) ist der realistischste Ansatz, erfordern aber enorme Rechenzeiten. MD Simulationsverfahren mit implizit modelliertem Wasser (implicit solvent: IS), bei denen den Effekt des Wassers durch ein angepaßtes Kraftfeld ersetzt wird, haben gegenüber ES- Verfahren einen großen Vorteil, da sie deutlich weniger Rechenzeit benötigen. Der Preis den man für diesen Vorteil zahlt, sind aber Abweichungen in den Simulationsdaten, die durch die Näherungen im IS Kraftfeld bedingt sind. Im Rahmen dieser Arbeit wurde ein Polyether-Kraftfeld für IS MD Simulationen entwickelt, das lokale Konformationen von 1,2-dimethoxy-ethane (DME) und globale Konformationen von Polyethylenglykol (PEG) akkurat beschreibt. Im IS Model fehlen die Wassermoleküle. Dies hat zur Folge, dass die anziehenden van der Waals (vdW) Wechselwirkungen der Atompaare in PEG einerseits und der Atompaare zwischen Wassermolekülen und PEG andererseits nicht ausbalanciert ist. Die im IS Model verbleibenden anziehenden vdW Wechselwirkungen zwischen den Atompaaren in PEG führen, verglichen mit einem ES Model, zu kompakteren Strukturen. Um die Balance der anziehenden vdW Kräfte im IS Model wieder herzustellen, werden entweder die Oberflächenenergieterme abgeschaltet und die vdW Wechselwirkungen auf 90% reduziert oder alternativ die Oberflächenenergieterme leicht negativ gesetzt, während die vdW Wechselwirkungen unverändert bleiben. Jedoch werden, um quantitativ die gleichen PEG Konformationen wie im ES-Model zu erreichen, zusätzliche Kraftfeldänderungen für die Torsionspotentiale und die Coulomb Wechselwirkungen der 1-4 und 1-5 Atompaare vorgeschlagen. Das vorgestellte IS Etherkraftfeld für MD Simulationen ist für Monomere wie für Polyether gleichermaßen gültig.