Use of Molecular Dynamics for the Refinement of an Electrostatic Model for the In Silico Design of a Polymer Antidote for the Anticoagulant Fondaparinux
Molecular dynamics (MD) simulations results are herein incorporated into an electrostatic model used to determine the structure of an effective polymer-based antidote to the anticoagulant fondaparinux. In silico data for the polymer or its cationic binding groups has not, up to now, been available, and experimental data on the structure of the polymer-fondaparinux complex is extremely limited. Consequently, the task of optimizing the polymer structure is a daunting challenge. MD simulations provided a means to gain microscopic information on the interactions of the binding groups and fondaparinux that would have otherwise been inaccessible. This was used to refine the electrostatic model and improve the quantitative model predictions of binding affinity. Once refined, the model provided guidelines to improve electrostatic forces between candidate polymers and fondaparinux in order to increase association rate constants. 1. Introduction While anticoagulation therapy is widely used, it has certain undesirable side effects such as the potential to cause life-threatening hemorrhages. Such bleeding complications can be mitigated, in the event of an overdose of anticoagulants, by the administration of antidotes which neutralize the anticoagulants while still avoiding thrombosis [1, 2]. The most commonly used anticoagulants are heparin-derived drugs [3], which include unfractionated heparin (UFH), low molecular weight heparins (LMWHs), and the synthetic pentasaccharide derivatives fondaparinux and idraparinux [4–7]. Because of its predictable dose response, almost complete bioavailability [4, 7], increased half-life [1], and no occurrence of heparin-induced thrombocytopenia [5], fondaparinux is becoming increasingly important in clinical medicine; however, its widespread use is limited by a lack of a specific antidote. Administration of protamine, the antidote for UFH and LMWHs, does not reverse the anticoagulant effect of fondaparinux, and hemodialysis only reduces fondaparinux plasma levels by 20% [1]. Hence, the development of a clinically safe antidote for this anticoagulant has become critical [8]. Currently, only limited experimental work has been reported for the development of an antidote to fondaparinux. It has been shown that heparinase I and the recombinant factor VII (rVIIa) can partially reverse fondaparinux in vitro; however, these studies were limited in scope: there is no clinical data for heparinase I, and there is only one volunteer study and one clinical case for rVIIa [1]. More recently, Borgel et al. have experimentally developed
References
[1]
S. Schulman and N. R. Bijsterveld, “Anticoagulants and their reversal,” Transfusion Medicine Reviews, vol. 21, no. 1, pp. 37–48, 2007.
[2]
J. H. Levy and K. A. Tanaka, “The anticoagulated patient: strategies for effective blood loss management,” Surgery, vol. 142, no. 4, pp. S71–77, 2007.
[3]
I. Capila and R. J. Lindhardt, “Heparin-protein interactions,” Angewandte Chemie International Edition, vol. 41, no. 3, pp. 390–412, 2002.
[4]
J. I. Weitz, “Emerging anticoagulants for the treatment of venous thromboembolism,” Thrombosis and Haemostasis, vol. 96, no. 3, pp. 274–284, 2006.
[5]
J. I. Weitz and L.-A. Linkins, “Beyond heparin and warfarin: the new generation of anticoagulants,” Expert Opinion on Investigational Drugs, vol. 16, no. 3, pp. 271–282, 2007.
[6]
L.-A. Linkins and J. I. Weitz, “New anticoagulants,” Seminars in Thrombosis and Hemostasis, vol. 29, no. 6, pp. 619–631, 2003.
[7]
S. Alban, “From heparins to factor Xa inhibitors and beyond,” European Journal of Clinical Investigation, vol. 35, no. 1, pp. 12–20, 2005.
[8]
B. K. Adler, “Unfractionated heparin and other antithrombin mediated anticoagulants,” Clinical Laboratory Science, vol. 17, no. 2, pp. 113–117, 2004.
[9]
E. P. Bianchini, J. Fazavana, V. Picard, and D. Borgel, “Development of a recombinant antithrombin variant as a potent antidote to fondaparinux and other heparin derivatives,” Blood, vol. 117, no. 6, pp. 2054–2060, 2011.
[10]
J. Fazavana, S. F. Bianchini, E. P. C. Smadja, V. Picard, M. Taverna, and D. Borgel, “A chemically-modified inactive antithrombin as a potent antagonist of fondaparinux and heparins anticoagulant activity,” Journal of Thrombosis and Haemostasis, 2013.
[11]
M. Remko, “Molecular structure, lipophilicity, solubility, absorption, and polar surface area of novel anticoagulant agents,” Journal of Molecular Structure, vol. 916, no. 1–3, pp. 76–85, 2009.
[12]
G. Scapin, “Structural biology and drug discovery,” Current Pharmaceutical Design, vol. 12, no. 17, pp. 2087–2097, 2006.
[13]
A. Cajiao, B. Gopaluni, E. Kwok, and J. N. Kizhakkedathu, “In silico design of polymeric antidote for anticoagulant fondaparinux,” Journal of Medical and Biological Engineering, vol. 31, no. 2, pp. 129–134, 2011.
Accelrys Software, Materials Studio, Online Help, Release 5. 5, 5. 5 edition, San Diego, Calif, USA, 2010.
[17]
H. Sun, “COMPASS: an ab initio force-field optimized for condensed-phase applications: overview with details on alkane and benzene compounds,” Journal of Physical Chemistry B, vol. 102, no. 38, pp. 7338–7364, 1998.
[18]
M. Remko and C.-W. von der Lieth, “Conformational structure of some trimeric and pentameric structural units of heparin,” Journal of Physical Chemistry A, vol. 111, no. 51, pp. 13484–13491, 2007.
[19]
M. Remko, M. Swart, and F. M. Bickelhaupt, “Conformational behavior of basic monomeric building units of glycosaminoglycans: isolated systems and solvent effect,” Journal of Physical Chemistry B, vol. 111, no. 9, pp. 2313–2321, 2007.
[20]
R. Tyler-Cross, M. Sobel, L. E. McAdory, and R. B. Harris, “Structure-function relations of antithrombin III-heparin interactions as assessed by biophysical and biological assays and molecular modeling of peptide-pentasaccharide-docked complexes,” Archives of Biochemistry and Biophysics, vol. 334, no. 2, pp. 206–213, 1996.
[21]
M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids, Oxford University Press, 1991.
[22]
S.-T. Lin, P. K. Maiti, and W. A. Goddard III, “Dynamics and thermodynamics of water in PAMAM dendrimers at subnanosecond time scales,” Journal of Physical Chemistry B, vol. 109, no. 18, pp. 8663–8672, 2005.
[23]
P. K. Maiti, T. ?a?in, S.-T. Lin, and W. A. Goddard III, “Effect of solvent and pH on the structure of PAMAM dendrimers,” Macromolecules, vol. 38, no. 3, pp. 979–991, 2005.
[24]
Y. Cui, “Parallel stacking of caffeine with riboflavin in aqueous solutions: the potential mechanism for hydrotropic solubilization of riboflavin,” International Journal of Pharmaceutics, vol. 397, no. 1-2, pp. 36–43, 2010.
[25]
T. Ba?tu? and S. Kuyucak, “Free energy simulations of single and double ion occupancy in gramicidin A,” Journal of Chemical Physics, vol. 126, no. 10, Article ID 105103, 2007.
[26]
D. Trzesniak, A.-P. E. Kunz, and W. F. van Gunsteren, “A comparison of methods to compute the potential of mean force,” ChemPhysChem, vol. 8, no. 1, pp. 162–169, 2007.
[27]
T. W. Allen, T. Ba?tu?, S. Kuyucak, and S.-H. Chung, “Gramicidin A channel as a test ground for molecular dynamics force fields,” Biophysical Journal, vol. 84, no. 4, pp. 2159–2168, 2003.
[28]
T. Selzer, S. Albeck, and G. Schreiber, “Rational design of faster associating and tighter binding protein complexes,” Nature Structural Biology, vol. 7, no. 7, pp. 537–541, 2000.
[29]
O. G. Berg and P. H. von Hippel, “Diffusion-controlled macromolecular interactions,” Annual Review of Biophysics and Biophysical Chemistry, vol. 14, pp. 131–160, 1985.
[30]
G. Schreiber and A. R. Fersht, “Rapid, electrostatically assisted association of proteins,” Nature Structural Biology, vol. 3, no. 5, pp. 427–431, 1996.
[31]
G. Schreiber, G. Haran, and H.-X. Zhou, “Fundamental aspects of protein: protein association kinetics,” Chemical Reviews, vol. 109, no. 3, pp. 839–860, 2009.
[32]
T. Selzer and G. Schreiber, “New insights into the mechanism of protein-protein association,” Proteins, vol. 45, no. 3, pp. 190–198, 2001.
[33]
H.-X. Zhou, “Disparate ionic-strength dependencies of on and off rates in protein-protein association,” Biopolymers, vol. 59, pp. 427–433, 2001.
[34]
M. Vijayakumar, K.-Y. Wong, G. Schreiber, A. R. Fersht, A. Szabo, and H.-X. Zhou, “Electrostatic enhancement of diffusion-controlled protein-protein association: comparison of theory and experiment on barnase and barstar,” Journal of Molecular Biology, vol. 278, no. 5, pp. 1015–1024, 1998.
[35]
J. T. H.G. Overbeek and B. H. Bijsterbosch, “Electrokinetic separation methods,” in The Electrical Double Layer and the Theory of Electrophoresis, pp. 1–33, Elsevier; North-Holland Biomedical Press, 1979.
[36]
S. Park and K. Schulten, “Calculating potentials of mean force from steered molecular dynamics simulations,” Journal of Chemical Physics, vol. 120, no. 13, pp. 5946–5961, 2004.
[37]
M. A. Kostiainen, J. G. Hardy, and D. K. Smith, “High-affinity multivalent DNA binding by using low-molecular-weight dendrons,” Angewandte Chemie, vol. 44, no. 17, pp. 2556–2559, 2005.
[38]
C. K. Nisha, S. V. Manorama, M. Ganguli, S. Maiti, and J. N. Kizhakkedathu, “Complexes of poly(ethylene glycol)-based cationic random copolymer and calf thymus DNA: a complete biophysical characterization,” Langmuir, vol. 20, no. 6, pp. 2386–2396, 2004.
[39]
J. K?stener, “Umbrella sampling,” Wiley Interdisciplinary Reviews, vol. 1, no. 6, pp. 932–942, 2011.
[40]
J. K?stner, “Umbrella integration in two or more reaction coordinates,” Journal of Chemical Physics, vol. 131, no. 3, Article ID 034109, 2009.