Like-Charge Guanidinium Pairing between Ligand and Receptor: An Unusual Interaction for Drug Discovery and Design?

A database survey in this study revealed for the first time that there are 227 counterintuitive like-charge guanidinium pairings (Gdm(+)-Arg pairings) between ligands and receptors in the Protein Data Bank, implying the potential guanidinium-arginine binding between guanidine-containing drugs and their target proteins. Furthermore, there are 145 guanidine-containing molecules in the DrugBank, showing the prevalence of guanidinium groups in drugs. It has also been reported that the introduction of a guanidinium group forming Gdm(+)-Arg pairing improved the potency of the drug by more than 8-fold in a typical case. On the basis of the survey, six ligand-protein complexes with typical Gdm(+)-Arg pairings were chosen for QM/MM calculations. The calculations at the B97-D/6-311++g(d,p) level revealed that the interaction could be as strong as -1.0 to -2.5 kcal/mol in DMSO and water, comparable to common intermolecular interactions. The calculations also unveiled that the Gdm(+)-Arg pairing interactions change from repulsive to attractive with the increase of dielectric constant, suggesting that the dielectric constant has a general stabilization effect on the Gdm(+)-Arg pairing. This study suggested that the like-charge guanidinium pairing interaction could be used not only for tuning the physical and chemical properties of drug leads but also for improving ligand binding affinity.

[1]  H. Scheraga,et al.  Contribution of unusual Arginine-Arginine short-range interactions to stabilization and recognition in proteins , 1994, Journal of protein chemistry.

[2]  A. Moscona,et al.  The use of a quantitative fusion assay to evaluate HN-receptor interaction for human parainfluenza virus type 3. , 1999, Virology.

[3]  K. Houk,et al.  Torsional steering controls the stereoselectivity of epoxidation in the guanacastepene a synthesis. , 2006, Organic letters.

[4]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[5]  B. Maigret,et al.  Monte Carlo simulations on the like-charged guanidinium-guanidinium ion pair in water , 1990 .

[6]  R. Abagyan,et al.  Unusual arginine formations in protein function and assembly: rings, strings, and stacks. , 2012, The journal of physical chemistry. B.

[7]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[8]  Jacopo Tomasi,et al.  The ONIOM-PCM method: Combining the hybrid molecular orbital method and the polarizable continuum model for solvation. Application to the geometry and properties of a merocyanine in solution , 2001 .

[9]  F. Albericio,et al.  Combinatorial approaches towards the discovery of new tryptase inhibitors. , 2005, Bioorganic & medicinal chemistry letters.

[10]  P. Gilbert,et al.  A comparative study of the bactericidal and growth inhibitory activities of the bisbiguanides alexidine and chlorhexidine. , 1989, The Journal of applied bacteriology.

[11]  F Guarnieri,et al.  A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins. , 1999, Biophysical journal.

[12]  R. J. Boyd,et al.  Effect of counterions on the protonation state in a poly(G)-poly(C) radical cation. , 2011, The journal of physical chemistry. B.

[13]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[14]  F. Uhlig,et al.  Like-Charge Ion Pairing in Water: An Ab Initio Molecular Dynamics Study of Aqueous Guanidinium Cations , 2012 .

[15]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[16]  Thom Vreven,et al.  Combining Quantum Mechanics Methods with Molecular Mechanics Methods in ONIOM. , 2006, Journal of chemical theory and computation.

[17]  Keiji Morokuma,et al.  Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants. , 2008, The journal of physical chemistry. B.

[18]  Mario Vazdar,et al.  Like-charge guanidinium pairing from molecular dynamics and ab initio calculations. , 2011, The journal of physical chemistry. A.

[19]  S. Plotkin,et al.  A theory for the anisotropic and inhomogeneous dielectric properties of proteins. , 2011, Physical chemistry chemical physics : PCCP.

[20]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[21]  F. Maseras Hybrid Quantum Mechanics/Molecular Mechanics Methods in Transition Metal Chemistry , 1999 .

[22]  J. Pitera,et al.  Dielectric properties of proteins from simulation: the effects of solvent, ligands, pH, and temperature. , 2001, Biophysical journal.

[23]  V. Gogonea The Qm/mm Method. an Overview , 2002 .

[24]  David S. Wishart,et al.  DrugBank 4.0: shedding new light on drug metabolism , 2013, Nucleic Acids Res..

[25]  Lenwood S. Heath,et al.  H++: a server for estimating pKas and adding missing hydrogens to macromolecules , 2005, Nucleic Acids Res..

[26]  Armando Rossello,et al.  Amber force field implementation, molecular modelling study, synthesis and MMP-1/MMP-2 inhibition profile of (R)- and (S)-N-hydroxy-2-(N-isopropoxybiphenyl-4-ylsulfonamido)-3-methylbutanamides. , 2006, Bioorganic & medicinal chemistry.

[27]  Alexander D. MacKerell,et al.  The structure of aqueous guanidinium chloride solutions. , 2004, Journal of the American Chemical Society.

[28]  C. Chipot,et al.  Effect of Polarizability on the Potential of Mean Force of Two Cations. The Guanidinium−Guanidinium Ion Pair in Water , 1997 .

[29]  M. D. de Jong,et al.  Neuraminidase Inhibitors and their Role in Avian and Pandemic Influenza , 2005, Antiviral therapy.

[30]  Sarah A Harris,et al.  A demonstration of the inhomogeneity of the local dielectric response of proteins by molecular dynamics simulations. , 2010, The Journal of chemical physics.

[31]  Weiliang Zhu,et al.  Halogen bonding--a novel interaction for rational drug design? , 2009, Journal of medicinal chemistry.

[32]  A. W. Johnston,et al.  Bethanidine, Guanethidine, and Methyldopa in Treatment of Hypertension: a Within-patient Comparison , 1968, British medical journal.

[33]  G. Vriend,et al.  Optimizing the hydrogen‐bond network in Poisson–Boltzmann equation‐based pKa calculations , 2001, Proteins.

[34]  D. E. Resetarits,et al.  Errors in chlorothiazide bioavailability estimates based on a Bratton--Marshall colorimetric method for chlorothiazide in urine. , 1979, Journal of pharmaceutical sciences.

[35]  R. Adamson,et al.  Biological Sciences: Hydroxyguanidine—a New Antitumour Drug , 1972, Nature.

[36]  Weiliang Zhu,et al.  The stabilization effect of dielectric constant and acidic amino acids on arginine-arginine (Arg-Arg) pairings: database survey and computational studies. , 2013, The journal of physical chemistry. B.

[37]  J Andrew McCammon,et al.  Optimized Radii for Poisson-Boltzmann Calculations with the AMBER Force Field. , 2005, Journal of chemical theory and computation.

[38]  Nobuaki Koga,et al.  Role of structural water molecule in HIV protease‐inhibitor complexes: A QM/MM study , 2008, J. Comput. Chem..

[39]  Vasiliy P. Mishin,et al.  Susceptibilities of Antiviral-Resistant Influenza Viruses to Novel Neuraminidase Inhibitors , 2005, Antimicrobial Agents and Chemotherapy.

[40]  Erik L. G. Wernersson,et al.  Guanidinium Cations Pair with Positively Charged Arginine Side Chains in Water , 2011 .

[41]  Thom Vreven,et al.  Geometry optimization with QM/MM, ONIOM, and other combined methods. I. Microiterations and constraints , 2003, J. Comput. Chem..

[42]  R. Saf,et al.  Structure-activity relationships of oligoguanidines-influence of counterion, diamine, and average molecular weight on biocidal activities. , 2003, Biomacromolecules.

[43]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[44]  Z. Maksić,et al.  Basicity of some organic superbases in acetonitrile. , 2001, Organic letters.

[45]  J. Esko,et al.  On Guanidinium and Cellular Uptake , 2014, The Journal of organic chemistry.

[46]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[47]  Xiaowu Chen,et al.  Neuraminidase Inhibitors as Anti-Influenza Virus Agents , 1999, Antiviral chemistry & chemotherapy.

[48]  T. Yumura,et al.  The Role of Weak Bonding in Determining the Structure of Thiophene Oligomers inside Carbon Nanotubes , 2012 .

[49]  F. Sa̧czewski,et al.  Biological activities of guanidine compounds, 2008 – 2012 update , 2013, Expert opinion on therapeutic patents.

[50]  M. Lawrence,et al.  Structure of the haemagglutinin-neuraminidase from human parainfluenza virus type III. , 2004, Journal of molecular biology.

[51]  Abhijit R. Tendulkar,et al.  Electrostatics‐defying interaction between arginine termini as a thermodynamic driving force in protein–protein interaction , 2009, Proteins.

[52]  Leming,et al.  Efficacy and safety of the neuraminidase inhibitor zanamivir in the treatment of influenzavirus infections , 1997 .

[53]  Rebecca C. Wade,et al.  Improving the Continuum Dielectric Approach to Calculating pKas of Ionizable Groups in Proteins , 1996 .

[54]  C David Sherrill,et al.  An Assessment of Density Functional Methods for Potential Energy Curves of Nonbonded Interactions: The XYG3 and B97-D Approximations. , 2010, Journal of chemical theory and computation.

[55]  Jirí Vondrásek,et al.  The molecular origin of like-charge arginine-arginine pairing in water. , 2009, The journal of physical chemistry. B.

[56]  B C Finzel,et al.  Three-dimensional structure of an antibody-antigen complex. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[57]  K. Baldridge,et al.  Implementation and Performance of DFT-D with Respect to Basis Set and Functional for Study of Dispersion Interactions in Nanoscale Aromatic Hydrocarbons. , 2008, Journal of chemical theory and computation.

[58]  Claude Lecomte,et al.  Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities , 1998 .

[59]  David Quiñonero,et al.  Structure and binding energy of anion-pi and cation-pi complexes: a comparison of MP2, RI-MP2, DFT, and DF-DFT methods. , 2005, The journal of physical chemistry. A.

[60]  H. Scheraga,et al.  Stability of Like and Oppositely Charged Organic Ion Pairs in Aqueous Solution , 1997 .

[61]  P. Colman,et al.  The structure of the complex between influenza virus neuraminidase and sialic acid, the viral receptor , 1992, Proteins.

[62]  F. Sa̧czewski,et al.  Biological activities of guanidine compounds , 2009, Expert opinion on therapeutic patents.

[63]  J. Zimmerberg,et al.  The Anti-Influenza Virus Agent 4-GU-DANA (Zanamivir) Inhibits Cell Fusion Mediated by Human Parainfluenza Virus and Influenza Virus HA , 2000, Journal of Virology.

[64]  L. Chavas,et al.  The catalytic mechanism of human parainfluenza virus type 3 haemagglutinin-neuraminidase revealed. , 2015, Angewandte Chemie.