Sulfur as an Acceptor to Bromine in Biomolecular Halogen Bonds

Department of Biochemistry & Molecular Biology, 1870 Campus Delivery, Colorado ...... Boys, S. F.; Bernardi, F., Calculation of Small Molecular Intera...
0 downloads 0 Views 2MB Size
Letter pubs.acs.org/JPCL

Sulfur as an Acceptor to Bromine in Biomolecular Halogen Bonds Melissa Coates Ford, Matthew Saxton, and P. Shing Ho* Department of Biochemistry & Molecular Biology, Colorado State University, 1870 Campus Delivery, Fort Collins, Colorado 80523-1870, United States S Supporting Information *

ABSTRACT: The halogen bond (X-bond) has become an important design element in chemistry, including medicinal chemistry and biomolecular engineering. Although oxygen is the most prevalent and best characterized X-bond acceptor in biomolecules, the interaction is seen with nitrogen, sulfur, and aromatic systems as well. In this study, we characterize the structure and thermodynamics of a Br···S X-bond between a 5bromouracil base and a phosphorothioate in a model DNA junction. The single-crystal structure of the junction shows the geometry of the Br···S to be variable, while calorimetric studies show that the anionic S acceptor is comparable to or slightly more stable than the analogous O acceptor, with a −3.5 kcal/mol difference in ΔΔH25°C and −0.4 kcal/mol ΔΔG25°C (including an entropic penalty ΔΔS25°C of −10 cal/(mol K)). Thus sulfur is shown to be a favorable acceptor for bromine X-bonds, extending the application of this interaction for the design of inhibitors and biological materials.

I

The most common X-bond acceptor in biomolecules is the carbonyl oxygen of the peptide bond, primarily because of its prevalence in proteins,4,27,28 although the amino acid side chains collectively are nearly equivalent in their representation.8,29 The sulfur-containing amino acids, particularly Cys, are underrepresented as X-bond acceptors,8,29 but Met has potential as a target for inhibitor design.14 Experimental and computational studies on small-molecule systems show sulfur as similar to or slightly more favorable than oxygen as an Xbond acceptor30,31 with geometries that deviate significantly from linearity.30,32 A study of iodobenzenes bound in cavities engineered into T4 lysozyme, however, indicated that the I···S or I···Se X-bond was not significantly different in energy compared with standard van der Waals interactions.33 Thus it is not clear how strong S would be as an X-bond acceptor in a biomolecular system. Our lab has developed a DNA Holliday junction model to study X-bonds in a biological context,10,34,35 in which a halogenated uracil (XU in a decanucleotide sequence d(CCGGTA6XU7CCGG), Figure 2) is placed at a nucleotide position to form an X-bond. This X-bond replaces a stabilizing H-bond to the phosphate backbone of nucleotide A6 at the sharp Uturn of the crossing strand of the four-stranded junction in the stacked-X geometry. The single-crystal structures of these constructs reveal the detailed geometries, while melting studies through differential scanning calorimetry (DSC) measures the energetics of the molecular interaction in solution. Finally, quantum-mechanical (QM) analyses on models of the interacting groups within the DNA junction help bridge the

nterest in the halogen bond as a tool in the biological sciences1−3 has grown since its rediscovery as a stabilizing interaction that both conveys specificity in inhibitor recognition4−9 and controls the structure and stability of nucleic acids10 and proteins.11 The term “halogen bond”12 (X-bond) has been adopted to reflect its analogous behavior to the better known hydrogen bond (H-bond). In particular, any H-bond acceptor can potentially serve as an X-bond acceptor.13 Most recent studies have focused on oxygen, nitrogen, and πelectron systems as X-bond acceptors because of their prevalence in biomolecules.8 Sulfur is also commonly seen as an acceptor4,8 and is gaining attention as a potential X-bond target in drug design,14 but there is currently very little understanding of its X-bonding potential in a biological system. Here we adopt a DNA Holliday junction as a model to determine the structure-energy relationship of Br···S X-bonds and find that an anionic sulfur is comparable to or slightly stronger as an X-bond acceptor relative to oxygen, even in geometries that are not as ideal for this type of interaction. The X-bond, similar to the H-bond, is considered to be primarily an electrostatically driven interaction,15−17 although charge transfer, dispersion, and polarization contribute significantly.18−20 In the electrostatic σ-hole model, the Xbond donor is an electropositive crown created at the tip of a halogen substituent as a consequence of depletion of electron density opposite a covalent σ-bond (e.g., a CX bond, Figure 1).15 The resulting “σ-hole” interaction with an electron-rich acceptor is shorter than the sum of the van der Waals radii (∑rvdW) of the respective atoms, with a linear approach of the acceptor to the halogen (θ1 angle). The σ-hole is most pronounced with the heavier halogens (Cl < Br < I; F is not a significant X-bond donor in biology) and is accentuated by electron-withdrawing groups. © XXXX American Chemical Society

Received: July 5, 2017 Accepted: August 10, 2017 Published: August 10, 2017 4246

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters

Figure 2. X-bonding interactions with oxygen as the acceptor at the crossover in a DNA junction. The overall junction (PDB code 2ORG10) is shown, with the inset depicting the 5-bromouracil base at Br U7 interacting with phosphate oxygen pA6. Replacement of an oxygen with a sulfur in this phosphate group, as a phosphorothioate in the junction (ps-J), creates a prochiral center at the phosphorus, defining the two prochiral configurations of the sulfur (labeled pro-Rp and pro-Rs), either of which could form a Br···S X-bond.

Figure 1. X-bond. (a) σ-Hole model for X-bonding.15 As a halogen becomes covalently bonded, for example, to a carbon, the pz-orbital of the halogen is depopulated when the valence electron is subsumed by the resulting σ-bond molecular orbital, creating a “σ-hole”. This hole accounts for the electropositive charge and polar flattening in the direction opposite the σ-bond. The electrons that remain in the perpendicular px,y-orbitals create an anisotropic charge distribution across the halogen’s surface. (b) Comparison of H- and X-bond donors and acceptors in biological systems and relationships between the two. Because both H- and X-bonds are primarily electrostatically driven interactions, they share a common set of electron-rich acceptors. When a H- and X-bond share a common acceptor, such as the carbonyl oxygen of the peptide backbone in a protein, the interactions are known to be orthogonally related (in terms of both geometry and their energy independence).21 The anisotropic charge distribution results in halogens being amphoteric, allowing them to serve simultaneously as X-bond donors and H-bond acceptors.13,22−25 Adapted with permission from refs 4 and 26. © 2004 National Academy of Sciences; © 2016 American Chemical Society.

Unlike previous work, all four strands of each construct are identical, forcing an X-bond to form in all structures instead of having one compete against the other. We will, however, use the previously published structure with a Br···O X-bond (PDB code 2ORG10) as a model for the interactions in p-J. The overall structure of the ps-J is that of the four-stranded DNA junction, similar to the stacked-X junctions seen previously (Figure 3a).10,51 The overall junction geometries (Jtwist, which relates the rotational angle between the stacked arms along their helical axes, and Jroll, which relates the stacked arms perpendicular to the helix axes) indicate that the global geometries of the two conformations are very similar (Table 1). The primary difference is that the stacked duplex arms of the sulfur containing junction are related by a much shallower angle (Jtwist ≈ 30° compared with ≥40° in all previous structures,36−40 including those with X-bonds51), making it a more planar junction structure. During refinement of the structure, it was clear that the sulfur of the phosphorothioate at A6 was the pro-Rp enantiomer at the core of the junction and positioned close to the bromine on BrU7. Careful analysis of the residual electron densities in the Fo-Fc difference map (Figure 3b), however, indicated that the ps-linkage assumed two conformations, labeled here as conformer A (ps-JA, with 52% occupancy) and conformer B (ps-JB, with 48% occupancy), and was subsequently confirmed by reductions in the Rwork and Rfree values when both conformations were included in the final model (with Rwork dropping from 22.5 to 21.8% and Rfree dropping from 27.8 to 27.6% with the addition of two conformations). Finally, the presence of these two conformations was further supported by analysis of the structural B-

structural details of the crystals with the associated stabilizing energies in solution. We apply this same DNA junction system here to determine the structure−energy relationship of a Br···S X-bond by introducing a phosphorothioate (ps) in place of the standard phosphate (p) linkage between A6 and XU7 (where X = Br), thereby allowing us for the first time to directly compare sulfur to oxygen as biomolecular X-bond acceptors. The goals of the current studies are to determine the structure−energy relationship of sulfur as an X-bond acceptor and to compare against the more prevalent oxygen acceptor in a biomolecular system. From our previous studies,10,51 bromine was found to be the overall most stabilizing X-bond donor in our model DNA junction to an anionic oxygen acceptor. We have thus designed two DNA junction constructs, the first in which the X-bond acceptor is the anionic oxygen of a standard phosphate (p) diester linkage between T5 and A6 in the sequence (CCGGT5pA6BrUCCGG, called p-J), and the second is the sulfur of a phosphorothioate (ps) in the sequence (CCGGT5psA6BrUCCGG, called ps-J). 4247

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters

The local structures of the T5 and the A6 nucleotides at the point of junction crossover differed among the two Br···S conformers. In particular, the base pair from T5 to its complementary A6 base is slightly stretched in ps-JA compared with ps-JB (Table 1). A more detailed analysis showed the sulfur of ps-JA approaches the Br of BrU7 at a closer distance (rBr···(S/O)) and more linearly (θ1) than that of ps-JB (Figure 3c,d, Table 2). It is clear from the geometries that ps-JA forms an X-bond, which pulls T5 base slightly away from the A6 base on the complementary strand. Table 2. Geometries of the Br···(O/S) X-bond phosphate pJ Junction Construct and the Two Conformers of the Phosphorothioate (ps-JA and ps-JB) Structures, Including the Distance between Interacting Atoms (rBr···(O/S)), the Distance as a Percentage of the Sum of the Respective van der Waals Radii (%∑rvdW), and the Angle of Approach of the Acceptor to the Br−C Bond (θ1)

Figure 3. X-ray crystal structure of phosphorothioate DNA junction (ps-J). (a) Phosphorothioate DNA junction (ps-J) crystallizes as a stacked-X four-stranded junction. The junction was seen to adopt two conformations (ps-JA in blue and ps-JB in orange). (b) Electron density map at the junction crossover. The 2Fo-Fc electron density map calculated at the 1σ level (blue wire) for the psA6 and brominated uracil (BrU7) nucleotides is shown, indicating the presence of the ps-JA (blue carbons) and ps-JB (orange carbons) conformers (c) Geometry of the Br···S interaction from the phosphothioate of A6 (yellow sphere) to Br of BrU7 at the junction core of ps-JA. The short distance and angle of approach are labeled for the Br···S interaction, along with distance to the associated water molecule (WA). (d) Geometry of the Br···S interaction at the junction core for ps-JB. The distance and angle of approach are labeled for the Br···S interaction.

a

p-Ja

ps-JA

Junction Parameters 40.7° 30.1° Jtwist Jroll 138.7° 129.6° Rotational Helix Parametersb helical twist (34.7°)c 36.85° 36.90° propeller twist (−12.0°)c −19.90° −12/20° tilt (−0.62°)c 0.15° 0.50° roll (1.74°)c 3.35° 2.60° Buckle (−0.23°)c −7.80° 8.10° Opening 3.10° 4.70° Translational Helix Parametersb rise (3.30 Å)c 3.29 Å 3.21 Å slide (0.66 Å)c 0.75 Å 0.79 Å shear 0.30 Å 0.59 Å stretch 0.12 Å −0.34 Å stagger −0.06 Å −0.60 Å shift 0.03 Å 0.03 Å

p-Ja

ps-JA

ps-JB

rBr···(O/S) %∑rvdW θ1

2.9 Å 86% 167.2°

3.5 Å 96% 161.3°

4.1 Å 112% 150.1°

Parameters for the p-J construct were from structure 2ORG.10

The geometry of the X-bond in ps-JA, however, is longer (in terms of the %Σrvdw) and less linear compared with the analogous Br···O interaction in the previously studied structure (PDB code 2ORG). The less linear angle of approach of the acceptor to the halogen, however, is consistent with results from previous work on sulfur acceptors in small-molecule assemblies,14 where the approach angles fall between 160 and 170°, as compared with those of oxygen acceptors,43 which fall between 170 and 180°. Furthermore, the X-bond of sulfur acceptors in the small-molecule systems has been shown to be more energetically favorable when compared with oxygen acceptors,30 suggesting that the angle of approach may not be as stringent in this large element compared with O. Finally, a water molecule was seen in ps-JA bridging the sulfur acceptor to the Br, which could confer additional stability to this interaction in the junction (Figure 3d). Dif ferential Scanning Calorimetry Comparison of S and O XBond Energies. The effect of the Br···S X-bond on the thermal stability of the junction in solution was determined through DSC melting studies. We had previously shown that junction formation is concentration-dependent, with the DNA remaining as a duplex at low concentrations and as a four-stranded junction at higher concentrations. At the concentration of DNA (300 μM) for the current study, both junction and duplex are present, allowing the melting parameters of the two forms to be determined simultaneously. The DSC melting and cooling profiles of the ps-J construct showed fully reversible melting and reannealing of the DNA and were best fit using two two-state van’t Hoff models. Subtracting the duplex energies from the junction leaves the energy of stabilization associated with the interactions at the junction crossover.35 Subsequent comparison of the crossover energies between ps-J and p-J allows us to determine the specific thermodynamic parameters (ΔΔH25°C, ΔΔS25°C, and ΔΔG25°C) of a sulfur relative an oxygen as an X-bond acceptor. The DSC melting data show that both the junction and duplex forms of the phosphothioate ps-J construct are more

Table 1. Comparison of Overall Junction Parameters (Jtwist and Jroll) and Helix Parameters for the T5·A6 Base Pair at the Crossover of the Phosphate p-J junction Construct and the Two Conformers of the Phosphorothioate (ps-JA and ps-JB) Structures parameter

parameter

ps-JB 30.2° 124.8° 36.95° −19/10° 1.30° −1.50° −17.10° 5.80° 3.27 Å 0.19 Å −0.93 Å 0.12 Å 0.54 Å 0.03 Å

Parameters for the p-J construct were from structure 2ORG.10 bHelix parameters were calculated using CURVES+.41 cValues in parentheses are the averages for all B-type structures from Hays et al.42 a

factors, which were significantly higher for the T5 and A6 residues with only a single model, but became consistent with the remainder of the structure with two conformers in the phasing model. 4248

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters

Table 3. Melting Temperatures (Tm), Enthalpies (ΔHm), and Entropies (ΔSm) Measured by Differential Scanning Calorimetry (DSC) for the p-J and the Phosphorothioate (ps-J) DNA constructs duplex

junction

DNA

Tm (°C)

ΔHm (kcal/mol)

ΔSm (cal/(mol K))

Tm (°C)

ΔHm (kcal/mol)

ΔSm (cal/(mol K))

p-J ps-J

51.9 ± 0.12 52.6 ± 0.10

53.2 ± 0.5 58.4 ± 0.3

163.5 ± 1.6 179 ± 0.9

57.2 ± 0.08 58.2 ± 0.13

70.7 ± 1.0 83.8 ± 0.6

214 ± 3 253.1 ± 1.8

thermally stable than the native p-J, with melting temperatures (Tm) that are higher by ∼0.7 °C for the duplex and by 1 °C for the junction (Table 3). Furthermore, ΔHm for ps-J duplex and junction is ∼6 and ∼13 kcal/mol, respectively, higher than the comparable oxygen-containing DNA. To directly compare the stabilizing potential of the sulfur Xbond acceptor to that of oxygen, the DSC melting parameters were translated to stabilizing enthalpies and entropies at a reference temperature of 25 °C (ΔH25°C and ΔS25°C, respectively), applying the change in heat capacity (ΔCp) in the standard relationships.44 Because the core structures are nearly identical for the oxygen and sulfur constructs, the difference in ΔH25°C and ΔS25°C reflects the difference in specific X-bonding potential for the two types of acceptors once the effects of the modifications on the duplexes have been accounted for. The resulting thermodynamic values (Table 4) show that the sulfur−bromine X-bond is enthalpi-

Table 5. Quantum-Mechanical Energies Calculated at the MP2 Level (EMP2), Applying the aug-cc-PVTZ Basis Set with BSSE Correction, for a 5-Bromouracil Complex with Dimethylphosphate (for 2ORG) or with Dimethylphosphorothioate (for ps-JA and ps-JB Conformers of the ps-J Construct) EMP2 (kcal/mol)

ΔS25°C(J−D), and ΔG25°C(J−D)), and Associated Differences between the Phosphorothioate and Phosphate Junction Constructs (ps-J−p-J) ΔH25°C(J−D) (kcal/mol)

ΔS25°C(J−D) (cal/(mol K))

ΔG25°C(J−D) (kcal/mol)

DNA construct

−5.6 ± 0.5 −9.1 ± 0.3 ΔΔH25°C (kcal/mol)

−15.6 ± 1.5 −25.7 ± 0.9 ΔΔS25°C (cal/(mol K))

−0.97 ± 0.06 −1.40 ± 0.04 ΔΔG25°C (kcal/mol)

ps-J−p-J

−3.5 ± 0.6

−10.2 ± 1.8

−0.43 ± 0.07

p-J ps-J

Br···S (ps-JA)

Br···S (ps-JB)

−2.22

−1.75

−0.94

the QM calculated and the DSC measured energies. The MP2 energies show that the Br···S X-bond of the ps-JA conformer is about twice as stabilizing as that of ps-JB. This is not surprising, given the shorter and more linear geometry of the ps-JA conformer. In contrast with the DSC results, however, the QM energies for the Br···S X-bond of ps-JA were slightly less stabilizing by ∼0.5 kcal/mol than that of Br···O. Thus the QM energies support the expectation that the more ideal geometry of ps-JA conformer is energetically more stable than the ps-JB conformer and the observation from the DSC studies that the sulfur is comparable to oxygen as an X-bond acceptor. In the current work, we have shown that an anionic sulfur can serve as a favorable acceptor for bromine X-bonds. Although previous studies have characterized the structure− energy relationship of sulfur X-bonds in small-molecule assemblies both experimentally and computationally, there has been a dearth of knowledge about sulfur acceptors in a biological context. In our DNA model system, a bromine Xbond with a sulfur acceptor was shown to be enthalpically more stabilizing than a comparable oxygen acceptor, despite what appears to be less favorable geometries. Both the energetics and broader range of geometries are consistent with prior computational studies on sulfur−halogen interactions.14,30 The stronger interaction, however, comes with an entropic cost, which renders sulfur overall only slightly more stabilizing as an X-bond acceptor. Sulfur-containing amino acids have been observed to comprise only 5% of X-bond acceptors in protein−ligand complexes,8 which may reflect the overall lower occurrence of the element at accessible regions in protein structures. Cysteine and methionine, however, are both prevalent in active sites of enzymes.45 Cysteines, for example, are key reactive nucleophiles for many hydrolases, deubiquitinases, caspases, and enzymes involved in redox reactions.46 Methionine is a particularly interesting in drug design due to its flexibility and lipophilicity. Similar to cysteine, it is found at higher rates in active sites of many enzymes, including nuclear hormone receptors, catechol-O-methyltransferase, and herpes simplex virus type 1 thymidine kinase.47 With the understanding now that sulfur is comparable to or more potent than oxygen as an X-bond acceptor, even in geometries that are not as short or linear, this underrepresented element becomes an attractive target when designing halogenated compounds as therapeutic agents.

Table 4. Stabilization Enthalpies, Entropies, and Free Energies Derived from DSC Melting Studies (Table 3), Extrapolated to 25 °C, for the Junction Minus Duplex Forms of the p-J and ps-J DNA Constructs (ΔH25°C(J−D),

DNA construct

Br···O (2ORG)

cally more stabilizing by −3.5 kcal/mol compared with oxygen. As had been previously seen,51 this more thermally stable interaction comes at an entropic cost of 10 cal/(mol K), which is reflected in the significantly lower crystallographic B-factors of the nucleotides at the junction crossing-overthose that participate in the X-bond (Figure S2). The resulting ΔΔG25°C, however, shows that sulfur remains overall a slightly stronger X-bond acceptor by 0.4 kcal/mol in this DNA system. QM Calculations on X-Bonding Energies. To associate the geometries of the X-bonds seen in the crystal structures with the stabilizing potentials determined by DSC, we performed QM calculations at the MP2 level on a set of model compounds (Figure S3) that mimic the Br···O interaction of the interaction in the previous 2ORG and the Br···S interaction of the current ps-J constructs (Table 5). We note that these MP2 calculations do not take into account the explicit solvent interactions, including the water seen in ps-JA that bridges the sulfur to the Br atom, or base stacking effects, all of which may account for the significant deviations between 4249

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters



METHODS DNA Synthesis and Purif ication. The DNA constructs for this study were designed as self-complementary sequences of the form 5′-CCGATXBrUCGG-3′ (where X has either a standard 2′-deoxyadenosine 5′-monophosphate, pA, or a 2-deoxyadenosine-5′-thiomonophosphate, psA). Chemically synthesized DNA oligonucleotides are purchased from Midland Certified Reagent Company attached to the solid controlled-pore glass (CPG) support, with the terminal dimethoxytrityl protecting group intact to allow for purification. The full-length product sequences were HPLC-purified and subsequently detritylated, as previously described.10,34 Phosphorothioate constructs required an additional HPLC step to separate the pro-Rp from pro-Rs prochiral enantiomers (Figure 2) using a protocol derived from that as described in Liu et al.48 Crystallization and Structure Determination. DNA constructs were crystallized by hanging drop vapor diffusion under conditions similar to those used previously to crystallize brominated DNA junctions.10 The initial 6 μL crystallization drops contained 0.75 mM DNA, 25 mM sodium cacodylate buffer (pH 7.0), 10−20 mM calcium chloride, and 0.9−1.2 mM spermine. Thin diamond-shaped crystals of the p-J and the pro-Rp ps-J enantiomer grew after 7 days of equilibration against a reservoir containing 30% 2-methyl-2,4-pentanediol for both constructs and 1% ß-mercaptoethanol for the psA construct. X-ray diffraction data on DNA crystals were collected using a Rigaku Compact Home Lab equipped with a copper sealedtube microfocus source and a PILATUS detector. Data were indexed, integrated, and scaled using the HKL2000 program (Table S1). Crystal structures were solved by molecular replacement using PHENIX,49 with a brominated junction (Protein Data Bank, PDB50 code 2ORG10) in which the bromine had been removed, serving as the initial phasing model. The final model (PDB code 5VBJ) was refined to Rcryst = 21.8% and Rfree = 27.6%. Dif ferential Scanning Calorimetry Studies. DSC was performed to determine energies of the sulfur-containing and non-sulfur-containing constructs, as previously described.34,35 DNA constructs (300 μM concentration in a solution containing 50 mM sodium cacodylate buffer, pH 7.0, and 15 mM calcium chloride) were annealed in the DSC instrument (TA Instruments Nano DSC) by heating to 90 °C for 20 min then allowing to cool to 10 °C at a rate of 0.9 °C/min. DNA melting energies were collected in the DSC, starting with equilibration for 400 s, and scanning from 10 to 95 °C at a rate of 0.9 °C/min at a pressure of 3.0 atm. Melting temperatures (Tm) and enthalpies of melting (ΔHm) were determined by fitting the data with TA Nano Analyze software. Previous work34,51 showed the presence of both junction and duplex for 300 μM DNA concentrations; therefore, a two-component (junction and duplex), two-state scaled model was used to analyze the data (Figure S1). Data for each construct were measured for at least seven replicates. Quantum-Mechanical Calculations. Models for QM energy calculations of the X-bonded junctions consisted of a 5brominouracil base interacting with either (dimethylphosphate) for the p-J construct or dimethylphosphorothioate (for the ps-J construct), with atomic position placed according to their respective crystal structures (Figure S2). The Møller− Plesset second-order (MP2) QM energies were calculated using Gaussian 09,52 with cyclohexane as the solvent (D = 2,

relative to a vacuum) to mimic the hydrophobic environment of the junction core, applying an appropriate polarizable basis set that includes dispersion (aug-cc-PVTZ) from the EMSL Basis Set Exchange.53 Basis set superposition errors (BSSEs) were determined from a separate counterpoise gas-phase calculation and directly summed into the calculated solvent phase energy.54



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.7b01725. Table S1. Parameters from Crystallization and Structure Solution of ps-J (PDB code 5VBJ). Figure S1. Sample differential scanning calorimetry (DSC) melting profile. Figure S2. Molecular models for QM calculations of Xbond interaction energies. Table S2. Atomic coordinates of the 5-bromouracil-dimethylphosphate complex for quantum-chemical calculations of Br···O X-bonding interaction. Table S3. Atomic coordinates of 5bromouracil-dimethylphosphorothioate complex for quantum-chemical calculations of Br···S X-bonding interaction in ps-JA conformation. Table S4. Atomic coordinates of 5-bromouracil-dimethylphosphorothioate complex for quantum -chemical calculations of Br···S Xbonding interaction in ps-JB conformation. (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Tel: 970-491-0569. Fax: 970491-0494. ORCID

Melissa Coates Ford: 0000-0002-0253-5389 P. Shing Ho: 0000-0002-8082-4311 Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS This work was supported by a grant from the National Science Foundation (CHE-1608146). REFERENCES

(1) 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 (4), 1363−1388. (2) 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. (3) Ho, P. S. Halogen Bonding in Medicinal Chemistry: From Observation to Prediction. Future Med. Chem. 2017, 9 (7), 637−640. (4) Auffinger, P.; Hays, F. A.; Westhof, E.; Ho, P. S. Halogen bonds in biological molecules. Proc. Natl. Acad. Sci. U. S. A. 2004, 101 (48), 16789−94. (5) Voth, A. R.; Ho, P. S. The Role of Halogen Bonding in Inhibitor Recognition and Binding by Protein Kinases. Curr. Top. Med. Chem. 2007, 7, 1336−1348. (6) Lam, P. Y. S.; Clark, C. G.; Smallwood, A. M.; Alexander, R. S. Structure-Based Drug Design Utilizing Halogen Bonding: Factor Xa Inhibitors. In The 238th ACS National Meeting; Metrangolo, P., Resnati, G., Eds.; ACS: Washington, DC, 2009; Vol. ORGN, p 58.

4250

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters (7) Lu, Y.; Liu, Y.; Xu, Z.; Li, H.; Liu, H.; Zhu, W. Halogen Bonding for Rational Drug Design and New Drug Discovery. Expert Opin. Drug Discovery 2012, 7 (5), 375−83. (8) 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 (11), 2781−2791. (9) Li, Y. S.; Guo, B. B.; Xu, Z. J.; Li, B.; Cai, T. T.; Zhang, X. B.; Yu, Y. Q.; Wang, H. Y.; Shi, J. Y.; Zhu, W. L. Repositioning Organohalogen Drugs: A Case Study for Identification of Potent BRaf V600E Inhibitors Via Docking And Bioassay. Sci. Rep. 2016, 6, 6. (10) Voth, A. R.; Hays, F. A.; Ho, P. S. Directing Macromolecular Conformation By Halogen Bonds. Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 6188−6193. (11) Scholfield, M. R.; Ford, M. C.; Carlsson, A. C.; Butta, H.; Mehl, R. A.; Ho, P. S. Structure-Energy Relationships of Halogen Bonds in Proteins. Biochemistry 2017, 56 (22), 2794−2802. (12) 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 (8), 1711−1713. (13) 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 (2), 255−264. (14) Wilcken, R.; Zimmermann, M. O.; Lange, A.; Zahn, S.; Kirchner, B.; Boeckler, F. M. Addressing Methionine in Molecular Design through Directed Sulfur-Halogen Bonds. J. Chem. Theory Comput. 2011, 7 (7), 2307−2315. (15) Clark, T.; Hennemann, M.; Murray, J. S.; Politzer, P. Halogen Bonding: The Sigma-Hole. J. Mol. Model. 2007, 13 (2), 291−296. (16) Murray, J. S.; Lane, P.; Politzer, P. Expansion of the SigmaHole Concept. J. Mol. Model. 2009, 15 (6), 723−729. (17) Murray, J. S.; Lane, P.; Clark, T.; Riley, K. E.; Politzer, P. Sigma-Holes, Pi-Holes and Electrostatically-Driven Interactions. J. Mol. Model. 2012, 18 (2), 541−548. (18) Eskandari, K.; Zariny, H. Halogen Bonding: A Lump-Hole Interaction. Chem. Phys. Lett. 2010, 492 (1−3), 9−13. (19) Rosokha, S. V.; Neretin, I. S.; Rosokha, T. Y.; Hecht, J.; Kochi, J. K. Charge-Transfer Character of Halogen Bonding: Molecular Structures and Electronic Spectroscopy of Carbon Tetrabromide and Bromoform Complexes with Organic Sigma- And Pi-Donors. Heteroat. Chem. 2006, 17 (5), 449−459. (20) Wang, C. W.; Danovich, D.; Mo, Y. R.; Shaik, S. On The Nature of the Halogen Bond. J. Chem. Theory Comput. 2014, 10 (9), 3726−3737. (21) 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. (22) Brammer, L.; Bruton, E. A.; Sherwood, P. Understanding the Behavior of Halogens as Hydrogen Bond Acceptors. Cryst. Growth Des. 2001, 1 (4), 277−290. (23) Aakeröy, C.; Desper, J.; Helfrich, B. A.; Metrangolo, P.; Pilati, T.; Resnati, G.; Stevenazzi, A. Combining Halogen Bonds and Hydrogen Bonds in the Modular Assembly of Heteromeric Infinite 1D Chains. Chem. Commun. 2007, 4236−4238. (24) Nelyubina, Y. V.; Antipin, M. Y.; Dunin, D. S.; Kotov, V. Y.; Lyssenko, K. A. Unexpected ″Amphoteric″ Character of the Halogen Bond: The Charge Density Study of the Co-Crystal of NMethylpyrazine Iodide with I-2. Chem. Commun. 2010, 46 (29), 5325−5327. (25) Zhou, P. P.; Qiu, W. Y.; Liu, S.; Jin, N. Z. Halogen as HalogenBonding Donor and Hydrogen-Bonding Acceptor Simultaneously in Ring-Shaped H3N.X(Y).HF (X = Cl, Br and Y = F, Cl, Br) Complexes. Phys. Chem. Chem. Phys. 2011, 13 (16), 7408−18. (26) Ford, M. C.; Ho, P. S. Computational Tools To Model Halogen Bonds in Medicinal Chemistry. J. Med. Chem. 2016, 59 (5), 1655−1670. (27) 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 (1), 314−318. (28) 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 (8), 935−945. (29) Zhang, Q.; Xu, Z.; Zhu, W. The Underestimated Halogen Bonds Forming with Protein Side Chains in Drug Discovery and Design. J. Chem. Inf. Model. 2017, 57 (1), 22−26. (30) Hauchecorne, D.; Moiana, A.; van der Veken, B. J.; Herrebout, W. A. Halogen Bonding to a Divalent Sulfur Atom: An Experimental Study of the Interactions of CF3X (X = Cl, Br, I) with Dimethyl Sulfide. Phys. Chem. Chem. Phys. 2011, 13 (21), 10204−13. (31) Li, Q. Z.; Jing, B.; Li, R.; Liu, Z. B.; Li, W. Z.; Luan, F.; Cheng, J. B.; Gong, B. A.; Sun, J. Z. Some Measures for Making Halogen Bonds Stronger than Hydrogen Bonds in H2CS-HOX (X = F, Cl, And Br) Complexes. Phys. Chem. Chem. Phys. 2011, 13 (6), 2266− 2271. (32) Lu, Y. X.; Zou, J. W.; Wang, Y. H.; Jiang, Y. J.; Yu, Q. S. Ab Initio Investigation of the Complexes Between Bromobenzene and Several Electron Donors: Some Insights into the Magnitude and Nature of Halogen Bonding Interactions. J. Phys. Chem. A 2007, 111 (42), 10781−8. (33) Liu, L.; Baase, W. A.; Matthews, B. W. Halogenated Benzenes Bound within a Non-Polar Cavity in T4 Lysozyme Provide Examples of I···S And I···Se Halogen-Bonding. J. Mol. Biol. 2009, 385 (2), 595− 605. (34) Carter, M.; Ho, P. S. Assaying the Energies of Biological Halogen Bonds. Cryst. Growth Des. 2011, 11 (11), 5087−5095. (35) Vander Zanden, C. M.; Carter, M.; Ho, P. S. Determining Thermodynamic Properties of Molecular Interactions from Single Crystal Studies. Methods 2013, 64, 12−18. (36) Eichman, B. F.; Vargason, J. M.; Mooers, B. H. M.; Ho, P. S. The Holliday Junction in an Inverted Repeat Sequence: Sequence Effects on the Structure of Four-Way Junctions. Proc. Natl. Acad. Sci. U. S. A. 2000, 97, 3971−3976. (37) Eichman, B. F.; Ortiz-Lombardia, M.; Aymami, J.; Coll, M.; Ho, P. S. The Inherent Properties of DNA Four-Way Junctions: Comparing the Crystal Structures of Holliday Junctions. J. Mol. Biol. 2002, 320 (5), 1037−51. (38) Hays, F. A.; Watson, J.; Ho, P. S. Caution! DNA Crossing: Crystal Structures of Holliday Junctions. J. Biol. Chem. 2003, 278 (50), 49663−6. (39) Khuu, P.; Ho, P. S. A Rare Nucleotide Tautomer in the Structure of an Asymmetric DNA Junction. Biochemistry 2009, 48, 7824−7832. (40) Vander Zanden, C. M.; Rowe, R. K.; Broad, A. J.; Robertson, A. B.; Ho, P. S. Effect of Hydroxymethylcytosine on the Structure and Stability of Holliday Junctions. Biochemistry 2016, 55 (41), 5781− 5789. (41) Lavery, R.; Moakher, M.; Maddocks, J. H.; Petkeviciute, D.; Zakrzewska, K. Conformational Analysis of Nucleic Acids Revisited: Curves+. Nucleic Acids Res. 2009, 37 (17), 5917−29. (42) Hays, F. A.; Teegarden, A.; Jones, Z. J.; Harms, M.; Raup, D.; Watson, J.; Cavaliere, E.; Ho, P. S. How Sequence Defines Structure: A Crystallographic Map of DNA Structure and Conformation. Proc. Natl. Acad. Sci. U. S. A. 2005, 102 (20), 7157−62. (43) Scholfield, M. R.; Ford, M. C.; Vander Zanden, C. M.; Billman, M. M.; Ho, P. S.; Rappe, A. K. Force Field Model of Periodic Trends in Biomolecular Halogen Bonds. J. Phys. Chem. B 2015, 119 (29), 9140−9149. (44) van Holde, K.; Johnson, W. C.; Ho, P. S. Principles of Physical Biochemistry, 2nd ed.; Pearson Prentice Hall: Upper Saddle River, NJ, 2006. (45) Khazanov, N. A.; Carlson, H. A. Exploring the Composition of Protein-Ligand Binding Sites on a Large Scale. PLoS Comput. Biol. 2013, 9 (11), e1003321. 4251

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252

Letter

The Journal of Physical Chemistry Letters (46) Hallenbeck, K. K.; Turner, D. M.; Renslo, A. R.; Arkin, M. R. Targeting Non-Catalytic Cysteine Residues Through StructureGuided Drug Discovery. Curr. Top. Med. Chem. 2016, 17 (1), 4−15. (47) Bissantz, C.; Kuhn, B.; Stahl, M. A Medicinal Chemist’s Guide to Molecular Interactions. J. Med. Chem. 2010, 53 (14), 5061−5084. (48) Liu, J.; Declais, A. C.; Lilley, D. M. Electrostatic Interactions and the Folding of the Four-Way DNA Junction: Analysis by Selective Methyl Phosphonate Substitution. J. Mol. Biol. 2004, 343 (4), 851−64. (49) Adams, P. D.; Afonine, P. V.; Bunkoczi, G.; Chen, V. B.; Davis, I. W.; Echols, N.; Headd, J. J.; Hung, L. W.; Kapral, G. J.; GrosseKunstleve, R. W.; McCoy, A. J.; Moriarty, N. W.; Oeffner, R.; Read, R. J.; Richardson, D. C.; Richardson, J. S.; Terwilliger, T. C.; Zwart, P. H. PHENIX: A Comprehensive Python-Based System for Macromolecular Structure Solution. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2010, 66, 213−221. (50) 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 (1), 235−242. (51) 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 (29), 4891−4903. (52) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F., et al. Gaussian 09, revision E.01; Gaussian, Inc.: 2009. (53) Schuchardt, K. L.; Didier, B. T.; Elsethagen, T.; Sun, L. S.; Gurumoorthi, V.; Chase, J.; Li, J.; Windus, T. L. Basis Set Exchange: A Community Database for Computational Sciences. J. Chem. Inf. Model. 2007, 47 (3), 1045−1052. (54) Boys, S. F.; Bernardi, F. Calculation of Small Molecular Interactions by Differences of Separate Total Energies - Some Procedures with Reduced Errors. Mol. Phys. 1970, 19 (4), 553−566.

4252

DOI: 10.1021/acs.jpclett.7b01725 J. Phys. Chem. Lett. 2017, 8, 4246−4252