Halogen Bonding in Nucleic Acid Complexes - American Chemical

Jun 12, 2017 - Michal H. Kolář*,†,§ and Oriana Tabarrini. ‡. †. Institute of Organic Chemistry and Biochemistry of the Czech Academy of Scien...
0 downloads 0 Views 2MB Size
Perspective Cite This: J. Med. Chem. 2017, 60, 8681-8690

pubs.acs.org/jmc

Halogen Bonding in Nucleic Acid Complexes Miniperspective Michal H. Kolár*̌ ,†,§ and Oriana Tabarrini‡ †

Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, 16610 Prague, Czech Republic ‡ Department of Pharmaceutical Sciences, University of Perugia, Via del Liceo 1, I-06123 Perugia, Italy

ABSTRACT: Halogen bonding (X-bonding) has attracted notable attention among noncovalent interactions. This highly directional attraction between a halogen atom and an electron donor has been exploited in knowledge-based drug design. A great deal of information has been gathered about X-bonds in protein−ligand complexes, as opposed to nucleic acid complexes. Here we provide a thorough analysis of nucleic acid complexes containing either halogenated building blocks or halogenated ligands. We analyzed close contacts between halogens and electron-rich moieties. The phosphate backbone oxygen is clearly the most common halogen acceptor. We identified 21 X-bonds within known structures of nucleic acid complexes. A vast majority of the X-bonds is formed by halogenated nucleobases, such as bromouridine, and feature excellent geometries. Noncovalent ligands have been found to form only interactions with suboptimal interaction geometries. Hence, the first X-bonded nucleic acid binder remains to be discovered.



INTRODUCTION Bringing a new drug to the market consumes enormous intellectual and financial resources. It typically takes more than a decade from the discovery of an active compound (lead) to the final approval of a drug, which is based on it. The drug discovery and development declined from the trial-and-error approach when the number of recognized diseases steeply rose. Nowadays, the workflow stems from knowledge gained in a variety of studies believing that this can speed up the whole drug development. The number of known drug targets is slightly higher than 300.1 This amount is not extraordinarily high bearing in mind the number of recognized human genes (less than 20 000 2). Among the targets, there are mostly proteins and only a few nucleic acids (NAs). Given that NAs are ubiquitous biopolymers with a myriad of cellular functions though, they represent a clinically prominent class of targets.1 Many strategies have appeared to optimize the lead compound into a therapeutic substance. The use of halogens is in this sense traditional. Xu et al. estimated that about 25% of approved drugs are halogenated, and the portion is similar in all stages of drug discovery and development.3 Halogen atoms modulate physicochemical properties of the molecular scaffolds; they affect the polarity and hydrophilicity/lipophilicity, which in turn changes membrane and blood−brain barrier permeation of the molecule. Also, carbon−halogen covalent bond is difficult to metabolize, so the halogenation © 2017 American Chemical Society

prolongs the lifetime of the active compound, but at the same time might increase its liver toxicity.4 Apart from nonspecific effects, it was about a decade ago that it was recognized that halogens might partake in a structurally specific and directional noncovalent interaction called a halogen bond (X-bond).5 The X-bond is an interaction between a halogen and a Lewis base or an electron-rich moiety. The electron density donors may be represented by electronegative atoms such as oxygen, nitrogen, sulfur, but also by aromatic rings or conjugated π-systems. The X-bond has been found in many protein−ligand complexes including pharmaceutically relevant ones (for review see refs 6−8). There have been only a few studies on X-bonds, where NAs played a role. The distinguished exceptions are the efforts of Shing Ho and co-workers. They focused on halogen bonding in so-called Holliday junction, a four-stranded (branched) complex of deoxyribonucleic acid (DNA).9,10 Using bromouridine as a building block, they directed a DNA into one of the several nearly isoenergetic conformers.11 The DNA model system also demonstrated that among the halogens it is the bromine that has the most favorable entropy/enthalpy compensation. Consequently, bromine was claimed an optimal element for X-bonding in DNA Holliday junctions.12,13 Received: March 1, 2017 Published: June 12, 2017 8681

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

protein−ligand X-bonds. The X-bonds are mostly established with the protein backbone. Its carbonyl oxygen is the most frequent electron donor.29 No preference for backbones of αhelices, β-sheets, or loop structures has been found.31 Lu et al. reported that about one-third of protein−ligand X-bonds involve an aromatic ring of the amino acid side chains.29 It was also reported that the X-bonds might not disrupt existing networks of H-bonds.32 Instead, an orthogonal pattern of Xand H-bonds is preferred involving the same electron donor atom.33 Another effect of ligand halogenation is a shortening of proximal H-bonds.34

To the best of our knowledge, no ligand has been reported so far to form an X-bond in any NA complex. The lack of information on X-bonding in NAs is somewhat surprising because NAs are naturally rich in electronegative atoms which make them (in theory) prospective X-bond acceptors. It seems that there is a missing relation between the worlds of NAs and X-bonds. This work aims at building such a relation. Hence, we continue with the introductions of the two worlds trying to find some overlap between them and highlight the pharmaceutical significance. Later, we analyze known structures of NAs and reveal main features of their complexes. Finally, we discuss all of the few examples of low-molecular compounds whose halogens are involved in interactions with electron donors.



NUCLEIC ACIDS AS DRUG TARGETS The deoxyribonucleic and ribonucleic acids (DNA and RNA) appear in living organisms in several forms. Whereas the DNA adopts only a few conformational classes, the structural diversity of RNA is much broader. DNA mostly occurs in the B-form helical conformation. A key RNA motif is the A-form double helix. However, loop regions, bulges, and other forms of mismatched nucleobases often disrupt the motif. Such a conformational variability is mirrored in the rapidly expanding variety of RNA functions recognized, especially in the past 2 decades. Apart from the classic roles of ribosomal, transfer, and messenger RNAs, RNA was also shown to store genetic information, regulate gene expression, or act as enzyme.35,36 Some NA-binders are halogenated. The DNA alkylating agents often contain halogen atoms because halogens facilitate the alkylating reactions on NAs. This way, anticancer cisplatin and its analogues37,38 covalently modify DNA. The reaction products block replication or transcription processes. The effect of alkylation is nonspecific in terms of DNA sequence and also cell type. The modifications occur in both normal and cancer cells, but the higher proliferation rate of the cancer cells makes these drugs effective. Alternative to covalent modifications, various classes of compounds bind to double-stranded DNA (dsDNA) in a noncovalent manner. Examples include the anticancer anthracycline type antibiotics daunomycin and adriamycin, and the polypeptide antibiotic dactynomicin, which works mainly by intercalating into DNA. The intercalation, described as an insertion of a planar molecule between consecutive base pairs, interferes with DNA-processing enzymes.39 Other small molecules bind to dsDNA interacting with base pair functional groups on the floor of the minor or major groove. Much of the research has concentrated on the DNA minor groove recognition leading to improved sequence selectivity of the compounds.40,41 Important classes of drugs target ribosomal RNA (rRNA). In fact, most of the RNA-targeting drugs on the market act on ribosomes.42 The ribosome is a biomachine that synthesizes proteins. As such, it represents a critical center of cellular life. The differences between bacterial and eukaryotic ribosomes have allowed developing specific antibacterial agents that are often based on natural products.43−45 There is a variety of binding locations within the ribosome; the most frequent sites are the peptidyl-transferase center on the 50S subunit and the decoding center on the 30S subunit. Aminoglycosides and tetracyclines are the best-known classes of small molecules whose primary target is the 30S subunit. Oxazolidonones instead exert their antibacterial action by binding to the 23S rRNA present within the 50S subunit. Macrolides and related compounds (lincosamides and streptogramins) as well as chloramphenicol and clindamycin also target the 50S subunit.



HALOGEN BONDING FEATURES The attraction of halogens with other electronegative atoms was observed as early as the 1950s in the crystallographic studies of Hassel et al.,14 although the synthesis of X-bonded complexes dates back to 19th century.15 The puzzling nature of the attraction between two electron-rich chemical groups was attributed merely to charge transfer effects until Politzer et al. came up with a simple model explaining many of the X-bond features.16 On the basis of quantum chemical calculations of the molecular electrostatic potentials (MEPs), they proposed that the surface of halogens contains regions of both positive and negative electrostatic potential (ESP).16 The positive region was labeled a σ-hole (Figure 1), and interestingly enough this

Figure 1. Molecular structure of bromouracil (a), the intersection of a molecular surface (b), and full ESP in atomic units (au) projected onto the molecular surface (c). The blue disk of positive ESP in the forefront is the σ-hole.

label appeared only 15 years after its first evidence.17 The Xbond has been exploited in many areas of chemistry and materials science (e.g., reviewed in refs 5 and 18), and a great amount of work has been done on the theoretical aspects of Xbonds too.6,19 X-bonds are similar in strength to the more common hydrogen bonds (H-bonds). The X-bond stabilization energy typically amounts to 5−25 kJ/mol,20,21 and it increases with the increasing atomic number of the halogen involved.22 The reason for this is the increasing size and magnitude of the halogen σ-hole.23 Another factor is the halogen polarizability,24 which also increases with the halogen atomic number. Contrary to traditional view, modern theoretical studies suggest that the role of charge transfer in the X-bond stabilization is modest.25 The strongest X-bonds are found in complexes of iodinated molecules; brominated and chlorinated molecules form weaker X-bonds. Fluorine is the least polarizable and the most electronegative halogen. It possesses positive σ-holes in rare cases and mostly in inorganic molecules,26 so it is of lower importance for biological applications. The structural trends of X-bonds in biomolecules have been inferred from several Protein Data Bank (PDB)27 surveys and theoretical analyses.7,28−30 All of the studies focused on 8682

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

Figure 2. 2D projections of the electrostatic potential (ESP) of the most common nucleobases A, C, G, T, U onto the planes of their aromatic rings. Note the ESP in atomic units (au) was calculated for the N-methylated analogs.

number of small molecules are able to bind the quadruplex structures. They are characterized by polycyclic heteroaromatic scaffolds or by cyclic/acyclic nonfused aromatic rings. Thus far, only a few molecules have been found to selectively bind the telomeric G4s, although their therapeutic potential appears high.59

Several ribosome-binding molecules have been prepared containing a halogen atom stemming mostly from the pioneering case study of chloramphenicol.46 In the past decade, other noncoding RNAs (ncRNAs) have emerged as prospective drug targets.42,47 Highly conserved ncRNAs provide new opportunities to expand the repertoire of drug targets to treat infections. Viral regulatory elements located in untranslated regions of mRNA often form folded structures that harbor potential binding sites for small molecules. The absence of homologous host cell RNAs makes them attractive for the development of innovative antiviral compounds. There are a plethora of ncRNAs under intense pharmaceutical research due to their ability to interact with low-molecular ligands. For instance, human immunodeficiency virus type 1 (HIV-1) transactivation response (TAR) RNA plays an essential role in HIV-1 replication through its interaction with the viral transactivator of transcription. Such interaction might be disrupted by ligands (reviewed in refs 48 and 49), where some of them contain halogen atoms.50 Another example is an internal ribosome entry site (IRES) from hepatitis C virus (HCV).51 The IRES RNA contains several independently folding domains that are potential targets for the development of selective viral translation inhibitors. A diverse set of ligands including oligonucleotides, peptides, and small molecules have been reported to block IRES function by distinct mechanisms.52 Like in the case of other ncRNAs, the drug candidates occasionally contain halogens. Nevertheless, no drug targeting a ncRNA has been approved for the market. Many lines of evidence are linking mutations and dysregulations of ncRNAs to neurodegenerative disorders.53,54 The presence of expanded CNG repeats in 5′ and 3′ untranslated regions is related to important diseases such as myotonic dystrophy type 1, spinocerebellar ataxia type 3, and fragile X-associated tremor-ataxia syndrome. For example, the pathological expansion of CAG repeats (>35 consecutive CAG codons) in huntingtin exon 1 encodes a mutant protein whose abnormal function determines Huntington’s disease. Finding compounds able to bind pathogenic CNG repeats with high specificity may be a valuable strategy against these devastating diseases. Their design is still limited by the lack of structural information, although some small molecules have emerged through various strategies.55,56 Both DNA and RNA can also adopt a noncanonical higherorder structure called G-quadruplexes (G4s) that are involved in regulating multiple biological pathways such as transcription, replication, translation, and telomere structure.57 The building blocks of G4s are guanosine-rich quartets that self-associate into a square-planar platform through a cyclic Hoogsteen Hbonded arrangement. G4s are found in oncogene promoters, in telomeres, and in introns of mRNAs. These regions have been recognized as potential targets for anticancer drugs.58 A large



NUCLEIC ACIDS AS X-BOND ACCEPTORS In the X-bond, the σ-hole on a halogen represents a Lewis acid that interacts with a Lewis base. NAs seem to offer an abundance of basic chemical groups. The backbone contains phosphate groups with two oxygens carrying a charge of −1. Likewise, ribose and deoxyribose contain oxygens with lone electron pairs that could serve as σ-hole acceptors too. Further, the nucleobases form H-bonds that could be potentially replaced or augmented by X-bonds. Nucleobases themselves show certain electrostatic diversity, where the negative sites may play a role of halogen acceptors. Figure 2 depicts MEPs of five most common nucleobases projected onto the plane of their aromatic systems. Cytosine and guanine contain areas of more negative ESP than the other nucleobases. Thymine and uracil resemble each other having two negative sites on the oxygens separated by a small positive site. Perhaps the most heterogeneous ESP is around adenine which exhibits six areas of zero ESP near the molecular surface, as compared to four (G, T, U) and two (C). Apart from this, all of the nucleobases are aromatic, which allows them to act as electron donors via their π-electrons above and below their rings (not shown). In nucleic acids, the situation is more complicated than in nucleobases: the three-dimensional structure of DNA and RNA is electrostatically diverse due to the presence of the phosphate backbone. For instance, the DNA minor groove was shown to be more electronegative than the major groove.60 Because of its plasticity, RNA may create folds with even more unusual electrostatic characteristics. Indeed, regions of strong negative electrostatic potentials were exploited in designing efficient TAR binders.61 Electrostatic interaction is also the main driving force of aminoglycoside binding to the ribosome.62 Overall, NAs seem to offer favorable electrostatics to attract positive σholes on halogens. How such ability is employed in ligand recognition and NA self-assembly remains elusive.



STRUCTURAL SURVEY YIELDS TWO DATA SETS To understand interaction preferences of halogens, we analyzed known NA structures. We started with a broad set of X-ray structures from the PDB (September 2016).27 Apart from X-ray structures, the PDB also contains NA structures determined by nuclear magnetic resonance (NMR) techniques. For two reasons, we deliberately omitted those from our analyses. First, we wanted to be consistent with the strategies adopted by previous structural surveys of X-bonds in protein−ligand 8683

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

80 000 protein structures deposited in the PDB at that time. The 21 X-bonds here represent about 0.2% of ∼9000 structures containing NA. Hence to better capture possible geometric properties of halogen interactions, we collected a more extended set of complexes with a longer interaction distance. An arbitrary threshold of 4 Å was chosen such that it is higher than the Xbond length but still short enough to hint for an attractive interaction. Within the article, the interactions are referred to as linear contacts and they comprehend the X-bonds as well. We found 72 linear contacts. Table 2 summarizes various subsets of the PDB query.

complexes. Second, there are indications that it may be difficult to assess the quality of the deposited NMR structural ensembles,63 which could complicate the geometric characterization. Within the X-ray structures, we searched for complexes containing a nucleic acid and a halogen atom (Cl, Br, or I). Fluorine was excluded from the search due to its extremely low ability to form X-bonds in biological systems. The selected structures comprise nucleic acids and their complexes with other nucleic acids, low-molecular ligands, and/or proteins. To get reliable geometric characteristics, only data with the resolution better than 3.0 Å were considered. Note that the previous PDB surveys of X-bonds in protein complexes used the same resolution threshold.28,30 From the PDB, we obtained 672 files that were subsequently filtered. We excluded structures containing halogens in the form of ions. Following the recommendation of the International Union of Pure and Applied Chemistry,64 we selected Xbonds as the contacts between halogens and electronegative atoms (N, O, P) shorter than or equal to the sum of van der Waals (vdW) radii65 (Table 1). The interactions involving at



X-BONDS FAVOR NUCLEOBASE···PHOSPHATE PATTERN Within the X-bond set, the variety of interacting partners is low; 20 X-bonds involve a halogenated nucleobase, one X-bond is formed by cisplatin chlorine. The set contains 19 X-bonds with a phosphate oxygen as the electron donor; further, there is one X-bond with an aromatic ring, and one with cytosine oxygen. Overall, only two X-bonds do not concur with the dominant nucleobase···phosphate interaction pattern. The X-bond geometries are close to ideal. The X-bond lengths are shorter than the sum of the vdW radii by about 8%, which conforms with the contractions reported for the protein X-bonds.7 About 90% of the X-bonds are straighter than 160°. The shortest X-bond in the set belongs to a structure of a Holliday junction; the X-bond is found between a bromodeoxyuridine and a phosphate oxygen of two neighboring residues (Figure 3a) (PDB code 2org, resolution 2.0 Å).11 The X-bond

Table 1. Sums of van der Waals Radii65 of Atomic Pairs in Å Cl Br I

N

O

P

3.30 3.40 3.53

3.27 3.37 3.50

3.55 3.65 3.78

least one of the two interacting atoms with the crystal occupancy lower than 0.5 were omitted. Because of the Xbond directionality,23 we also required the angle of R−X···Y to be higher than 120°. Same or similar geometric criteria were used previously to define biological X-bonds with proteins,28,29 although the wider angular range was used in other studies as well.7,30 We also searched for the X-bonds that are formed with the aromatic systems of the nucleobases. To this aim, we considered only halogen contacts closer than 5 Å to the aromatic plane that make an angle between the plane normal vector and the X−C bond smaller than 60°. In the end, we obtained a set of 21 X-bonds that satisfy the data quality, chemical and geometric criteria. This amount is a rather low. Scholfield et al. reported 760 protein complexes with an X-bond in 2012,7 which stands for about 1% of the

Figure 3. Shortest (a) and the straightest (b) X-bonds found. Hydrogens were omitted for clarity. Residues are labeled as follows: BRU, 5-bromouridine; dBRU, 5-bromo-2-deoxyuridine; U, uridine; dA, deoxyadenosine.

Table 2. Analysis of the PDBa Cl

Br

I

sum

402 43 31 22 2

204 162 135 29 17

66 34 22 2 2

672 239 188 53 21

611 43 22 2

1319 315 48 17

205 41 2 2

2135 399 72 21

files PDB search count contains contact(s) contains contact(s) with NA building block contains linear contacts(s) with NA building block contains X-bond(s) with NA building block interactions number of contacts contacts with NA building block linear contacts with NA building block X-bonds with NA building block a

The contact is defined by the interatomic distance between a halogen and oxygen, nitrogen, or phosphorus shorter than 4 Å. The linear contact has the angle R−X···O/P/N higher than 120°. X-bond is a linear contact shorter than the sum of vdW radii. 8684

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

In the NA complexes, the geometric quality of the halogen interactions increases in the order of Cl < Br < I. The interaction angles are more linear for heavier halogens. Especially the contacts of chlorine are rather bent with the median angle of 141°. The medians of the interaction lengths (Table 3) decrease with increasing atomic number of the halogen. The surveys on protein−ligand X-bonds28,29 revealed the opposite trend, i.e., the increasing length of X-bonds with the increasing atomic number of the halogen involved. In this work (but also in ref 28), the statistical sample might be insufficient to provide reliable statistics, which is true especially for the two iodine X-bonds. Figure 5 shows the histograms of the interaction lengths and angles. In the histogram of lengths, the minor peak near 3.0 Å,

was shown to stabilize particular DNA assembly in competition with H-bond. The straightest X-bond also involves a brominated uridine and phosphate oxygen (Figure 3b) (PDB code 2bu1, resolution 2.2 Å) in a complex of RNA and phage MS2 coat protein,66 and the overall geometry is remarkably similar to the one found in Holliday junction. Whereas the study on Holliday junction fully appreciated the role of Xbonding, in the latter case the specific interaction of the bromine remained unrecognized. What is important, we have identified no noncovalent ligand involved in X-bonding. The only non-nucleobase residue that participates in an X-bond is cisplatin covalently bound to an adenine (PDB code 5j4c, resolution 2.8 Å 67). It is also the only halogen donor that forms an X-bond with a π-system, although the interaction with one of the guanine nitrogens would alone classify the interaction as X-bond. The X-bond features superb geometric characteristics (Figure 4). Unlike proteins, in nucleic

Figure 4. X-bond with an aromatic system between cisplatin (CPT) and a guanine (G). The distance stands for the perpendicular distance from the aromatic plane.

acids, the halogen−π interaction competes with π−π interactions more often. Especially in dsDNA, it is hardly conceivable that there is a space for a halogen to attack the nucleobases from above or below of their aromatic planes. The situation in RNAs might be more favorable, but the single occurrence of such X-bond is hard to generalize.

Figure 5. Histograms of interaction lengths (a, c) and interaction angles represented by the R−X···Y angle (b, d); color-coded according to the halogen involved (a, b) or electron donating atom (c, d).



INTERACTIONS LONGER THAN X-BONDS Table 3 summarizes counts of various types of interactions and the statistics of the interaction geometries of the set of 72 halogen linear contacts.

which is comparable with the sum of the vdW radii, stands for almost ideal X-bond length. There is a major peak near 3.5 Å too. Two peaks also appear in the histogram of interaction angles, around 140° and 170°. Each of the halogen subsets contributes to the histograms with a different weight. Chlorine- and bromine-containing complexes span the whole range of interaction lengths. The two iodine complexes feature short interactions. The situation with angles is different. Chlorine complexes appear only in the region of lower interaction angles. The highest angle in the chlorine subset is 154°, so chlorine interactions are likely weak. On the other hand, bromine complexes are scattered across the whole range of angles (Figure 5b) with a cumulation around 170° (Figure 6). The two iodine X-bonds are very straight (Xbond angle higher than 170°), suggesting a strong interaction. We analyzed the electron donors of halogen interactions found in the NA complexes. The backbone oxygen atom is the most common one. We identified 51 linear contacts of halogens with oxygen (71%) and 21 with nitrogen (29%). Although we included phosphorus as an electron donor in the search, no linear contact was found. The actual occurrence of the oxygen

Table 3. Medians (med) and Interquartile Ranges (iqr) of Geometric Characteristics of Interactions Involving Various Types of Partners

halogen chlorine bromine iodine electron donor N nonbackbone O backbone O halogenated residue halogenated nucleotide low-molecular-weight ligand

count

length [Å] (med ± iqr)

angle [deg] (med ± iqr)

22 48 2

3.49 ± 0.31 3.45 ± 0.30 3.07 ± 0.06

141 ± 8 165 ± 22 173 ± 3

21 9 42

3.47 ± 0.28 3.69 ± 0.42 3.41 ± 0.45

140 ± 11 135 ± 25 166 ± 8

51 21

3.43 ± 0.43 3.50 ± 0.29

165 ± 23 140 ± 7 8685

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

Figure 6. Interaction geometries in polar coordinates color-coded according to the halogen (a) or the electron donor atom (b).

and nitrogen in NAs is roughly 2:1 (O/N), which likely contributes to the dominance of the linear contacts with oxygen. Most of the nitrogen interactions were found with chlorine, whereas bromine preferentially interacts with oxygen. 82% of the all interacting oxygens belong to the phosphate backbone. Unlike other oxygens and nitrogens in NAs, the phosphate oxygens carry a negative charge, which explains their higher propensity to halogen σ-holes. According to the electron-donating atoms, the geometric quality of the interactions increases in the order of nonbackbone oxygen < nitrogen < backbone oxygen (Table 3). The interactions that employ a backbone oxygen are typically shorter and straighter than the others. We conclude that different interaction atoms likely occur in different interaction geometries. It is also apparent on a projection of the interaction geometries to the polar coordinates (Figure 6). We observed no correlation between interaction lengths and the corresponding angles.

Figure 7. Structural formulas of several low-molecular-weight ligands or their building blocks that participate in halogen interactions.



HALOGENATED LIGANDS SHOW SUBOPTIMAL INTERACTION GEOMETRIES We did not find any X-bonded noncovalent ligands in NA complexes. Nevertheless, in the set of linear contacts we have found 15 unique ligands: 14 chlorinated and one brominated. The remaining interactions involve halogenated nucleobases as halogen donors. The halogen interactions of ligands are longer and less linear compared to the interactions of halogenated nucleobases (Table 3). The interaction angles deviate notably from linearity, which suggests that such interactions do not play a critical role in the NA−ligand recognition. From the pharmaceutical point of view, the low-molecularweight ligands are of higher interest than the halogenated building blocks. Below we discuss all of the instances among the linear-contact data set which involve a low-molecular ligand. Structural formulas of the ligands discussed are shown in Figure 7. Several halogenated compounds bind to the bacterial ribosome. Chloramphenicol (1)46 is a classic antibiotic compound that targets the ribosomal A-site crevice in the 50S subunit (PDB code 1nji, resolution 3.0 Å;68 PDB code 4v7w, resolution 3.0 Å 69). 1 contains two aliphatic chlorines that are activated by the nearby carbonyl group. Two distinct binding orientations in two different ribosomal system were proposed (T. thermophilus68 and H. marismortui69) (Figure 8). The lengths of both chlorine interactions are beyond the respective sums of the vdW radii (about 3.3 Å, Table 1) and quite bent. One of the chlorines approaches either the ribose

Figure 8. Chloramphenicol 1 (CLM) interactions in the 50S ribosomal subunit with guanosine (G) and uridine (U). Hydrogens were omitted for clarity. The overlay of the two binding poses found in T. thermophilus (PDB code 1nji) and H. marismortui (PDB code 4v7w) into a single reference frame is shown in red and black, respectively. Only residues within 7 Å from each of the halogens are shown.

in-ring oxygen (3.9 Å, 154°) or a guanine nitrogen (3.8 Å, 126°), respectively. The geometries suggest that the interactions of halogens contribute weakly to the complex stability. Clindamycin (2)70 is another halogenated ribosome binder71 (PDB code 1yjn, resolution 3.0 Å) that binds into the 50S subunit (Figure 9a) and forms linear contacts. 2 contains one chlorine that is directed toward the sugar edge of guanosine. There is an interaction with guanine nitrogen (3.5 Å, 148°). At the same time, there is a shorter but more bent contact with sugar O2′ oxygen (3.1 Å, 112°). 7-Chlorotetracycline (3)72 belongs to a class of tetracycline antibiotics that bind into the ribosomal 30S subunit. This class acts by preventing correct processing of aminoacyl-tRNA. Although no complex of 3 with the ribosome satisfies the data set criteria, a complex with an RNA aptamer does. The aptamer was designed to bind 3 with subnanomolar affinity (PDB code 3egz, resolution 2.2 Å 73). It features an interaction between a chlorine and a ribose oxygen of a cytidine (Figure 9b) (3.7 Å, 151°). Due to the suboptimal geometry, the role of halogen interaction in complex stabilization is likely marginal. Moreover, there are many other intermolecular interactions between 8686

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

Figure 9. Selected halogen interactions of ligand complexes. Clindamycin 2 (CLY) bound to the ribosome (a), chlorotetracycline 3 (CTC) bound to an RNA aptamer (b), cisplatin 4 (CPT) bound to ribosome (c), an intasome binder (ZZX) (d), a brominated imidazole analogue 7 (2ZZ) in an RNA aptamer (e). Hydrogens were omitted for clarity. The interaction lengths and R−X···Y angles are provided.

Besides this, it forms a contact with a phosphate oxygen of 3.5 Å under the angle of 137° (Figure 9e).

the 3 and the aptamer, such as a stacking of a planar part of 3 with a guanine, and an oxygen−magnesium coordination (not shown). Those reduce the relative importance of the halogen contact even more. Recently, crystal structures of T. thermophilus ribosomes with cisplatin37 (4 in Figure 9c) were resolved (PDB code 5j4b, resolution 2.6 Å; PDB code 5j4c, resolution 2.8 Å 67). Nine molecules of cisplatin are covalently bound to the ribosome in place of one of the chlorines. Two of the cisplatin residues were identified to interfere with the mRNA tunnel and the GTPase site, the places critical for the ribosome functions.67 One cisplatin out of the nine forms an X-bond with a guanine πsystem (Figure 4). Another cisplatin directs with its chlorine in between two stacked adenines with the shortest distance to a nitrogen of 4.0 Å and angle 151°. Another nitrogen is about 4.2 Å far away and forms the angle of 148° (Figure 9c). Our data set contains several ligands that bind into the prototype foamy virus (PFV) intasome, i.e., a complex of a viral integrase and DNA. Such a system serves as a model for HIV-1 integrase inhibition. Several inhibitors were developed with a halogenated phenyl ring (5 and 6) that form X-bonds with NAs.74−76 No difference in the interaction geometries was found between 5 and 6 analogs. 5-analog interactions show the best geometric characteristics among the PFV ligands (but generally still modest); the chlorine on 5 interacts with two adenine nitrogens (Figure 9d). While the contact length of 3.3 Å is comparable with the sum of vdW radii, the angle of 142° is far from the ideal 180°. The other contact is 3.5 Å long and even more bent (angle about 114°). We add that an aromatic fluorine activates the neighboring chlorine by making its σ-hole more positive.23 The effect of fluorination on the ligand binding affinity is not straightforward, however.77 We found only one brominated ligand with a contacts shorter than 4 Å. A brominated imidazole analog (7) binds into a G-quadruplex formed within an RNA aptamer called Spinach (PDB code 4q9q, resolution 2.45 Å).78 The Spinach module allows a green-fluorescence-protein-like functionality by selectively activating fluorescence of 7. 7 binds in a planar conformation with a strong stacking interaction with adenine.



SUMMARY AND OUTLOOK The role of X-bonding in drug development has been emphasized for a decade. Nonetheless, the knowledge-based design of halogenated drugs that feature an X-bond with their biomolecular targets is still a tedious task. Although there have been several published pioneering studies that involved protein targets,77,79,80 the nucleic acids still lack a successful application. In this work, we focused on X-ray structures of biomolecular complexes containing halogen interactions with nucleic acids. The general ability of nucleic acids to provide electron donating groups to X-bonding has been confirmed here. Using criteria recommended by IUPAC for X-bond definition, we found 21 NA complexes with the X-bond. We also analyzed interactions of the halogen atoms longer than X-bonds. All but one of the X-bonds involved a halogenated nucleobase (preferably uridine). Halogenated nucleotides are often utilized in X-ray crystallography, which explains their high occurrence in the NA complexes. The incorporation of heavy atoms such as bromine into the structure helps to overcome the phase problem. Moreover, the halogenated nucleotides alter the physicochemical properties of the NAs and thus modulate the crystallization conditions; in the classical theory this is due to their stacking properties,81 but recently X-bonding has also been shown to play a role.11 What is more, halogenated nucleobases may also be used in radiation anticancer therapy as radiosensitizers.82 Although the exact mechanism of the increased radiosensitivity is not known, the changes in the ionization potential brought by the halogens or global structural changes associated with modified nucleobase pairing have been proposed.83,84 An important feature of X-bonds in NA complexes is their preference toward backbone oxygens. This is true also for the longer halogen interactions. The phosphate oxygens tend to form straighter and shorter X-bonds than other electron donors, possibly due to their negative charge. In NA···ligand noncovalent recognition, a variety of interactions play a role. H-bonding and London dispersion8687

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

Present Address

driven stacking interactions belong to the classic interactions, complemented by salt bridges, water-network disruptions, or metal-ion interactions. Strictly speaking, no noncovalent ligand has been found to employ X-bonding in the NA binding (according to the IUPAC recommendations). No X-bonded NA ligand found is an intriguing fact, which may be explained in two ways: (i) The X-bonding ligands may be inefficient in NA recognition. If so, the halogenated ligands that bind into NAs bind due to other kinds of interactions that are stronger than Xbonds. Our data support this hypothesis. All of the halogen interactions of ligands found in the survey here show suboptimal geometries. Given the high X-bond directionality,85 the stabilization energy of a complex reduces notably when deviating from the ideal geometry. Hence, the contribution to the stabilization of the NA complexes is likely not dominant. Of course, there are exceptions where X-bonding is the driving interaction. For instance, it was proven that a single X-bond/Hbond exchange might transform the DNA conformation completely.11 (ii) Medicinal chemists may underappreciate the role of Xbonding in the drug-design strategies. Our data also support this hypothesis because we found only a few halogenated ligands: 20 chlorinated, 8 brominated, and 1 iodinated. The current study on its own is unable to dissect which reason is more likely. The starting set of PDBs was biased toward the lighter halogens, which are, however, less suitable for strong X-bonds. For a better understanding of the role of Xbonding, structural and functional studies are required, especially of the less-frequent brominated and iodinated compounds. The effects of halogens are not limited to the X-bonding though. The halogenation affects the global electron distribution and consequently many of the molecular properties. For instance, Fanfrlı ́k et al. demonstrated that solvation/desolvation may compensate for favorable σ-hole···lone pair interaction in halogen-to-hydrogen substitution in a protein−ligand complex.86 Also, a PDB survey of two protein families revealed shortening on H-bonds proximal to ligand halogens.34 This corresponds with the notion that halogenation increases the acidity of proximal H-bond donors. Still, these effects in NAs are difficult to inspect with the limited statistical sample of Xbonds here and of halogenated NA ligands in general. Nevertheless, the current analyses may help in designing novel X-bonded NA binders. Such binders should contain a bromine or iodine to form a strong X-bond. The ligand binding should be directed to the vicinity of sugar−phosphate backbone, for example, to groove or bulge regions. Consequently, such a ligand should likely not contain extensive aromatic systems that are prone to intercalation into canonical helices. Systematic halogenation of known pharmacologically active compounds may be an optimal strategy to identify new agents with enhanced activity, presumably supported by the specific X-bonds.



§

M.H.K.: Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, D-37077 Göttingen, Germany.

Notes

The authors declare no competing financial interest. Biographies Michal H. Kolár ̌ received his Ph.D. in 2013 from Charles University in Prague, Czech Republic, and from the Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic. With Pavel Hobza he focused on theoretical and computational description of noncovalent interactions. He is a recipient of the Humboldt Research Fellowship for Postdoctoral Researchers. From 2014 he worked with Paolo Carloni in Forschungszentrum Jülich, Germany, on RNA−ligand recognition. In 2016 he moved to Helmut Grubmüller’s group at the Max Planck Institute for Biophysical Chemistry, Göttingen, Germany, pursuing his interest in noncoding RNAs and computer simulations. Oriana Tabarrini is a Professor of Medicinal Chemistry at the Department of Pharmaceutical Sciences (University of Perugia, Italy). After the degree in Pharmaceutical Chemistry and Technology, she has worked for several years with research grants from pharmaceutical industries. In 2002 she was promoted to Associate Professor. Her research has mainly aimed at developing small molecules as pharmacological tools and potential chemotherapeutics with particular focus on nucleic acid binders as antivirals and more recently for the treatment of neurodegenerative disorders. She has published over 85 research articles in leading peer-reviewed journals, including some invited reviews and patents. She has also received several project grants and is an ad hoc reviewer for several top journals.



ACKNOWLEDGMENTS M.H.K. is thankful to N. Leioatts and A. Obarska-Kosinska for valuable discussions on the manuscript and acknowledges the support provided by the Alexander von Humboldt Foundation. Most of the analyses were prospered from NumPy87 and Matplotlib88 Python libraries.



ABBREVIATIONS USED BRU, 5-bromouridine; CLM, chloramphenicol; CLY, clindamycine; CPT, cisplatin; DNA, deoxyribonucleic acid; dsDNA, double-stranded DNA; G4, G-quadruplex; H-bond, hydrogen bond; HCV, hepatitis C virus; HIV-1, human immunodeficiency virus type 1; IRES, internal ribosome entry site; NA, nucleic acid; ncRNA, noncoding RNA; PDB, Protein Data Bank; PFV, prototype foamy virus; RNA, ribonucleic acid; TAR, transactivation response RNA element; X-bond, halogen bond



REFERENCES

(1) Overington, J. P.; Al-Lazikani, B.; Hopkins, A. L. How many drug targets are there? Nat. Rev. Drug Discovery 2006, 5, 993−996. (2) Ezkurdia, I.; Juan, D.; Rodriguez, J. M.; Frankish, A.; Harrow, M.; Diekhans, J.; Vazquez, J.; Valencia, A.; Tress, M. L. Multiple evidence strands suggest that there may be as few as 19 000 human proteincoding genes. Hum. Mol. Genet. 2014, 23, 5866−5878. (3) Xu, Z.; Yang, Z.; Liu, Y.; Lu, Y.; Chen, K.; Zhu, W. Halogen bond: its role beyond drug-target binding affinity for drug discovery and development. J. Chem. Inf. Model. 2014, 54, 69−78. (4) Hernandes, M. Z.; Cavalcanti, S. M. T.; Moreira, D. R. M.; de Azevedo, W. F., Jr.; Leite, A. C. L. Halogen atoms in the modern medicinal chemistry: hints for the drug design. Curr. Drug Targets 2010, 11, 303−314.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Michal H. Kolár:̌ 0000-0003-0841-944X 8688

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

(5) Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The halogen bond. Chem. Rev. 2016, 116, 2478−2601. (6) Ford, M. C.; Ho, P. S. Computational tools to model halogen bonds in medicinal chemistry. J. Med. Chem. 2016, 59, 1655−1670. (7) Scholfield, M. R.; Zanden, C. M. V.; Carter, M.; Ho, P. S. Halogen bonding (X-bonding): a biological perspective. Protein Sci. 2013, 22, 139−152. (8) Wilcken, R.; Zimmermann, M. O.; Lange, A.; Joerger, A. C.; Boeckler, F. M. Principles and applications of halogen bonding in medicinal chemistry and chemical biology. J. Med. Chem. 2013, 56, 1363−1388. (9) Duckett, D. R.; Murchie, A. I.; Diekmann, S.; von Kitzing, E.; Kemper, B.; Lilley, D. M. The structure of the Holliday junction, and its resolution. Cell 1988, 55, 79−89. (10) Lilley, D. M. Structures of helical junctions in nucleic acids. Q. Rev. Biophys. 2000, 33, 109−159. (11) Voth, A. R.; Hays, F. A.; Ho, P. S. Directing macromolecular conformation through halogen bonds. Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 6188−6193. (12) Carter, M.; Ho, P. S. Assaying the energies of biological halogen bonds. Cryst. Growth Des. 2011, 11, 5087−5095. (13) Carter, M.; Voth, A. R.; Scholfield, M. R.; Rummel, B.; Sowers, L. C.; Ho, P. S. Enthalpy-entropy compensation in biomolecular halogen bonds measured in DNA junctions. Biochemistry 2013, 52, 4891−4903. (14) Hassel, O.; Hvoslef, J. The structure of bromine 1,4-dioxanate. Acta Chem. Scand. 1954, 8, 873−873. (15) Guthrie, F. XXVIII.On the iodide of iodammonium. J. Chem. Soc. 1863, 16, 239−244. (16) Brinck, T.; Murray, J. S.; Politzer, P. Surface electrostatic potentials of halogenated methanes as indicators of directional intermolecular interactions. Int. J. Quantum Chem. 1992, 44, 57−64. (17) Clark, T.; Hennemann, M.; Murray, J. S.; Politzer, P. Halogen bonding: the σ-hole. J. Mol. Model. 2007, 13, 291−296. (18) Gilday, L. C.; Robinson, S. W.; Barendt, T. A.; Langton, M. J.; Mullaney, B. R.; Beer, P. D. Halogen bonding in supramolecular chemistry. Chem. Rev. 2015, 115, 7118−7195. (19) Kolár,̌ M. H.; Hobza, P. Computer modeling of halogen bonds and other σ-hole interactions. Chem. Rev. 2016, 116, 5155−5187. (20) Ř ezác,̌ J.; Riley, K. E.; Hobza, P. Benchmark calculations of noncovalent interactions of halogenated molecules. J. Chem. Theory Comput. 2012, 8, 4285−4292. (21) Kolár,̌ M. H.; Deepa, P.; Ajani, H.; Pecina, A.; Hobza, P. Characteristics of a σ-hole and the nature of a halogen bond. In Halogen Bonding II; Metrangolo, P., Resnati, G., Eds.; Springer International Publishing: Cham, Switzerland, 2014; pp 1−25. (22) Politzer, P.; Murray, J. S.; Clark, T. Halogen bonding: an electrostatically-driven highly directional noncovalent interaction. Phys. Chem. Chem. Phys. 2010, 12, 7748−7757. (23) Kolár,̌ M.; Hostaš, J.; Hobza, P. The strength and directionality of a halogen bond are co-determined by the magnitude and size of the σ-hole. Phys. Chem. Chem. Phys. 2014, 16, 9987−9996. (24) Riley, K. E.; Hobza, P. Investigations into the nature of halogen bonding including symmetry adapted perturbation theory analyses. J. Chem. Theory Comput. 2008, 4, 232−242. (25) Ř ezác,̌ J.; de la Lande, A. On the role of charge transfer in halogen bonding. Phys. Chem. Chem. Phys. 2017, 19, 791−803. (26) Metrangolo, P.; Murray, J. S.; Pilati, T.; Politzer, P.; Resnati, G.; Terraneo, G. Fluorine-centered halogen bonding: a factor in recognition phenomena and reactivity. Cryst. Growth Des. 2011, 11, 4238−4246. (27) Berman, H. M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T. N.; Weissig, H.; Shindyalov, I. N.; Bourne, P. E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235−242. (28) Auffinger, P.; Hays, F. A.; Westhof, E.; Ho, P. S. Halogen bonds in biological molecules. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 16789−16794.

(29) Lu, Y.; Wang, Y.; Zhu, W. Nonbonding interactions of organic halogens in biological systems: implications for drug discovery and biomolecular design. Phys. Chem. Chem. Phys. 2010, 12, 4543−4551. (30) Sirimulla, S.; Bailey, J. B.; Vegesna, R.; Narayan, M. Halogen interactions in protein-ligand complexes: implications of halogen bonding for rational drug design. J. Chem. Inf. Model. 2013, 53, 2781− 2791. (31) Wilcken, R.; Zimmermann, M. O.; Lange, A.; Zahn, S.; Boeckler, F. M. Using halogen bonds to address the protein backbone: a systematic evaluation. J. Comput.-Aided Mol. Des. 2012, 26, 935−945. (32) Rowe, R. K.; Ho, P. S. Relationships between hydrogen bonds and halogen bonds in biological systems. Acta Crystallogr., Sect. B: Struct. Sci., Cryst. Eng. Mater. 2017, 73, 255−264. (33) Voth, A. R.; Khuu, P.; Oishi, K.; Ho, P. S. Halogen bonds as orthogonal molecular interactions to hydrogen bonds. Nat. Chem. 2009, 1, 74−79. (34) Poznański, J.; Poznańska, A.; Shugar, D. A protein data bank survey reveals shortening of intermolecular hydrogen bonds in ligandprotein complexes when a halogenated ligand is an H-bond donor. PLoS One 2014, 9, e99984. (35) Kiss, T. Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functions. Cell 2002, 109, 145−148. (36) Mercer, T. R.; Dinger, M. E.; Mattick, J. S. Long non-coding RNAs: insights into functions. Nat. Rev. Genet. 2009, 10, 155−159. (37) Rosenberg, B.; Van Camp, L.; Krigas, T. Inhibition of cell division in Escherichia coli by electrolysis products from a platinum electrode. Nature 1965, 205, 698−699. (38) Weiss, R. B.; Christian, M. C. New cisplatin analogues in development. Drugs 1993, 46, 360−377. (39) Baguley, B. DNA intercalating anti-tumour agents. Anti-Cancer Drug Des. 1991, 6, 1−35. (40) Baraldi, P. G.; Bovero, A.; Fruttarolo, F.; Preti, D.; Tabrizi, M. A.; Pavani, M. G.; Romagnoli, R. DNA minor groove binders as potential antitumor and antimicrobial agents. Med. Res. Rev. 2004, 24, 475−528. (41) Barrett, M. P.; Gemmell, C. G.; Suckling, C. J. Minor groove binders as anti-infective agents. Pharmacol. Ther. 2013, 139, 12−23. (42) Thomas, J. R.; Hergenrother, P. J. Targeting RNA with small molecules. Chem. Rev. 2008, 108, 1171−1224. (43) Hermann, T. Drugs targeting the ribosome. Curr. Opin. Struct. Biol. 2005, 15, 355−366. (44) Poehlsgaard, J.; Douthwaite, S. The bacterial ribosome as a target for antibiotics. Nat. Rev. Microbiol. 2005, 3, 870−881. (45) Wilson, D. N. Ribosome-targeting antibiotics and mechanisms of bacterial resistance. Nat. Rev. Microbiol. 2014, 12, 35−48. (46) Ehrlich, J.; Bartz, Q. R.; Smith, R. M.; Joslyn, D. A. Chloromyeetin, a new antibiotic from a soil actinomycete. Science (Washington, DC, U. S.) 1947, 106, 417. (47) Hermann, T. Small molecules targeting viral RNA. WIREs RNA 2016, 7, 726−743. (48) Mousseau, G.; Mediouni, S.; Valente, S. T. Targeting HIV transcription: the quest for a functional cure. In The Future of HIV-1 Therapeutics; Metrangolo, B., Metrangolo, D. S., Metrangolo, D., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp 121− 145. (49) Tabarrini, O.; Desantis, J.; Massari, S. Recent advances in the identification of Tat-mediated transactivation inhibitors: progressing toward a functional cure of HIV. Future Med. Chem. 2016, 8, 421−442. (50) Mayer, M.; James, T. L. NMR-based characterization of phenothiazines as a RNA binding scaffold. J. Am. Chem. Soc. 2004, 126, 4453−4460. (51) Lukavsky, P. J. Structure and function of HCV IRES domains. Virus Res. 2009, 139, 166−171. (52) Dibrov, S. M.; Parsons, J.; Carnevali, M.; Zhou, S.; Rynearson, K. D.; Ding, K.; Garcia Sega, E.; Brunn, N. D.; Boerneke, M. A.; Castaldi, M. P.; Hermann, T. Hepatitis C virus translation inhibitors targeting the internal ribosomal entry site. J. Med. Chem. 2014, 57, 1694−1707. 8689

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690

Journal of Medicinal Chemistry

Perspective

(53) Castel, A. L.; Cleary, J. D.; Pearson, C. E. Repeat instability as the basis for human diseases and as a potential target for therapy. Nat. Rev. Mol. Cell Biol. 2010, 11, 165−170. (54) Mirkin, S. M. Expandable DNA repeats and human disease. Nature 2007, 447, 932−940. (55) Kumar, A.; Parkesh, R.; Sznajder, L. J.; Childs-Disney, J. L.; Sobczak, K.; Disney, M. D. Chemical correction of pre-mRNA splicing defects associated with sequestration of muscleblind-like 1 protein by expanded r (CAG)-containing transcripts. ACS Chem. Biol. 2012, 7, 496−505. (56) Bochicchio, A.; Rossetti, G.; Tabarrini, O.; Krauss, S.; Carloni, P. Molecular view of ligands specificity for CAG repeats in antiHuntington therapy. J. Chem. Theory Comput. 2015, 11, 4911−4922. (57) Lipps, H. J.; Rhodes, D. G-quadruplex structures: in vivo evidence and function. Trends Cell Biol. 2009, 19, 414−422. (58) Han, H.; Hurley, L. H. G-quadruplex DNA: a potential target for anti-cancer drug design. Trends Pharmacol. Sci. 2000, 21, 136−142. (59) Neidle, S. The structures of quadruplex nucleic acids and their drug complexes. Curr. Opin. Struct. Biol. 2009, 19, 239−250. (60) Chin, K.; Sharp, K. A.; Honig, B.; Pyle, A. M. Calculating the electrostatic properties of RNA provides new insights into molecular interactions and function. Nat. Struct. Biol. 1999, 6, 1055−1061. (61) Davis, B.; Afshar, M.; Varani, G.; Murchie, A. I.; Karn, J.; Lentzen, G.; Drysdale, M.; Bower, J.; Potter, A. J.; Starkey, I. D.; Swarbrick, T.; Aboul-ela, F. Rational design of inhibitors of HIV-1 TAR RNA through the stabilisation of electrostatic hot spots. J. Mol. Biol. 2004, 336, 343−356. (62) Ma, C.; Baker, N. A.; Joseph, S.; McCammon, J. A. Binding of aminoglycoside antibiotics to the small ribosomal subunit: a continuum electrostatics investigation. J. Am. Chem. Soc. 2002, 124, 1438−1442. (63) Nabuurs, S. B.; Spronk, C. A. M.; Vuister, G. W.; Vriend, G. Traditional biomolecular structure determination by NMR spectroscopy allows for major errors. PLoS Comput. Biol. 2006, 2, e9. (64) Desiraju, G. R.; Ho, P. S.; Kloo, L.; Legon, A. C.; Marquardt, R.; Metrangolo, P.; Politzer, P.; Resnati, G.; Rissanen, K. Definition of the halogen bond (IUPAC recommendations 2013). Pure Appl. Chem. 2013, 85, 1711−1713. (65) Bondi, A. Van der Waals volumes and radii. J. Phys. Chem. 1964, 68, 441−451. (66) Grahn, E.; Moss, T.; Helgstrand, C.; Fridborg, K.; Sundaram, M.; Tars, K.; Lago, H.; Stonehouse, N. J.; Davis, D. R.; Stockley, P. G.; Liljas, L. Structural basis of pyrimidine specificity in the MS2 RNA hairpin-coat-protein complex. RNA 2001, 7, 1616−1627. (67) Melnikov, S. V.; Söll, D.; Steitz, T. A.; Polikanov, Y. S. Insights into RNA binding by the anticancer drug cisplatin from the crystal structure of cisplatin-modified ribosome. Nucleic Acids Res. 2016, 44, 4978−4987. (68) Hansen, J. L.; Moore, P. B.; Steitz, T. A. Structures of five antibiotics bound at the peptidyl transferase center of the large ribosomal subunit. J. Mol. Biol. 2003, 330, 1061−1075. (69) Bulkley, D.; Innis, C. A.; Blaha, G.; Steitz, T. A. Revisiting the structures of several antibiotics bound to the bacterial ribosome. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 17158−17163. (70) Birkenmeyer, R. D.; Kagan, F. Lincomycin. XI. Synthesis and structure of clindamycin, a potent antibacterial agent. J. Med. Chem. 1970, 13, 616−619. (71) Tu, D.; Blaha, G.; Moore, P. B.; Steitz, T. A. Structures of MLS B K antibiotics bound to mutated large ribosomal subunits provide a structural explanation for resistance. Cell 2005, 121, 257−270. (72) Duggar, B. M. Aureomycin and Preparation of Same. U.S. Patent 2,482,055, 1949. (73) Xiao, H.; Edwards, T. E.; Ferré-D’Amaré, A. R. Structural basis for specific, high-affinity tetracycline binding by an in vitro evolved aptamer and artificial riboswitch. Chem. Biol. 2008, 15, 1125−1137. (74) Hare, S.; Vos, A. M.; Clayton, R. F.; Thuring, J. W.; Cummings, M. D.; Cherepanov, P. Molecular mechanisms of retroviral integrase inhibition and the evolution of viral resistance. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 20057−20062.

(75) Métifiot, M.; Maddali, K.; Johnson, B. C.; Hare, S.; Smith, S. J.; Zhao, X. Z.; Marchand, C.; Burke, T. R., Jr; Hughes, S. H.; Cherepanov, P.; Pommier, Y. Activities, crystal structures, and molecular dynamics of dihydro-1H-isoindole derivatives, inhibitors of HIV-1 integrase. ACS Chem. Biol. 2013, 8, 209−217. (76) Raheem, I. T.; et al. Discovery of 2-pyridinone aminals: A prodrug strategy to advance a second generation of HIV-1 integrase strand transfer inhibitors. J. Med. Chem. 2015, 58, 8154−8165. (77) Fanfrlík, J.; Kolár,̌ M.; Kamlar, M.; Hurný, D.; Ruiz, F. X.; Cousido-Siah, A.; Mitschler, A.; Ř ezác,̌ J.; Munusamy, E.; Lepšík, M.; Matějíček, P.; Veselý, J.; Podjarny, A.; Hobza, P. Modulation of aldose reductase inhibition by halogen bond tuning. ACS Chem. Biol. 2013, 8, 2484−2492. (78) Huang, H.; Suslov, N. B.; Li, N.-S.; Shelke, S. A.; Evans, M. E.; Koldobskaya, Y.; Rice, P. A.; Piccirilli, J. A. A G-quadruplex-containing RNA activates fluorescence in a GFP-like fluorophore. Nat. Chem. Biol. 2014, 10, 686−691. (79) Llinas-Brunet, M.; et al. Discovery of a potent and selective noncovalent linear inhibitor of the hepatitis C virus NS3 protease (BI 201335). J. Med. Chem. 2010, 53, 6466−6476. (80) Hardegger, L. A.; Kuhn, B.; Spinnler, B.; Anselm, L.; Ecabert, R.; Stihle, M.; Gsell, B.; Thoma, R.; Diez, J.; Benz, J.; Plancher, J.-M.; Hartmann, G.; Banner, D. W.; Haap, W.; Diederich, F. Systematic investigation of halogen bonding in protein-ligand interactions. Angew. Chem., Int. Ed. 2011, 50, 314−318. (81) Bugg, C. E.; Thewalt, U. Effects of halogen substituents on base stacking in nucleic acid components: the crystal structure of 8bromoguanosine. Biochem. Biophys. Res. Commun. 1969, 37, 623−629. (82) Erikson, R. L.; Szybalski, W. Molecular radiobiology of human cell lines: V. Comparative radiosensitizing properties of 5-halodeoxycytidines and 5-halodeoxyuridines. Radiat. Res. 1963, 20, 252−262. (83) Wetmore, S. D.; Boyd, R. J.; Eriksson, L. A. A theoretical study of 5-halouracils: electron affinities, ionization potentials and dissociation of the related anions. Chem. Phys. Lett. 2001, 343, 151− 158. (84) Heshmati, E.; Abdolmaleki, P.; Mozdarani, H.; Sarvestani, A. S. Effects of halogen substitution on Watson-Crick base pairing: A possible mechanism for radiosensitivity. Bioorg. Med. Chem. Lett. 2009, 19, 5256−5260. (85) Kolár,̌ M. H.; Carloni, P.; Hobza, P. Statistical analysis of σholes: a novel complementary view on halogen bonding. Phys. Chem. Chem. Phys. 2014, 16, 19111−19114. (86) Fanfrlík, J.; Ruiz, F. X.; Kadlčíková, A.; Ř ezác,̌ J.; Cousido-Siah, A.; Mitschler, A.; Haldar, S.; Lepšík, M.; Kolár,̌ M. H.; Majer, P.; Podjarny, A. D.; Hobza, P. The effect of halogen-to-hydrogen bond substitution on human aldose reductase inhibition. ACS Chem. Biol. 2015, 10, 1637−1642. (87) Van Der Walt, S.; Colbert, S. C.; Varoquaux, G. The NumPy array: a structure for efficient numerical computation. Comput. Sci. Eng. 2011, 13, 22−30. (88) Hunter, J. D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90−95.

8690

DOI: 10.1021/acs.jmedchem.7b00329 J. Med. Chem. 2017, 60, 8681−8690