Enthalpy–Entropy Compensation in Biomolecular Halogen Bonds Measured in DNA Junctions

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Modulation of Rotational Dynamics in Halogen-Bonded Cocrystalline Solids

Dynamic processes are responsible for the functionality of a range of materials, biomolecules, and catalysts. We report a detailed systematic study of the modulation of methyl rotational dynamics via the direct and the indirect influence of noncovalent halogen bonds. For this purpose, a novel series of cocrystalline architectures featuring halogen bonds (XB) to tetramethylpyrazine (TMP) is designed and prepared using gas-phase, solution, and solid-state mechanochemical methods. Single-crystal X-ray diffraction reveals the capacity of molecular bromine as well as weak chloro-XB donors to act as robust directional structure-directing elements. Methyl rotational barriers (Ea) measured using variable-temperature deuterium solid-state NMR range from 3.75 ± 0.04 kJ mol-1 in 1,3,5-trichloro-2,4,6-trifluorobenzene·TMP to 7.08 ± 0.15 kJ mol-1 in 1,4-dichlorotetrafluorobenzene·TMP. Ea data for a larger series of TMP cocrystals featuring chloro-, bromo-, and iodo-XB donors are shown to be governed by a combination of steric and electronic factors. The average number of carbon-carbon close contacts to the methyl group is found to be a key steric metric capable of rationalizing the observed trends within each of the Cl, Br, and I series. Differences between each series are accounted for by considering the strength of the σ-hole on the XB donor. One possible route to modulating dynamics is therefore via designer cocrystals of variable stoichiometry, maintaining the core chemical features of interest between a given donor and acceptor while simultaneously modifying the number of carbon close contacts affecting dynamics. These principles may provide design opportunities to modulate more complex geared or cascade dynamics involving larger functional groups.

On the Importance of Halogen and Chalcogen Bonds in the Solid State of Nucleic Acids: A Combined Crystallographic and Theoretical Perspective

In this work, intra- and intermolecular halogen and chalcogen bonds (HlgBs and ChBs, respectively) present in the solid state of nucleic acids (NAs) have been studied at the RI-MP2/def2-TZVP level of theory. To achieve this, a Protein Data Bank (PDB) survey was carried out, revealing a series of structures in which Br/I or S/Se/Te atoms belonging to nucleobases or pentose rings were involved in noncovalent interactions (NCIs) with electron-rich species. The energetics and directionality of these NCIs were rationalized through a computational study, which included the use of Molecular Electrostatic Potential (MEP) surfaces, the Quantum Theory of Atoms in Molecules (QTAIM), and Non Covalent Interaction plot (NCIplot) and Natural Bonding Orbital (NBO) techniques.

A Comprehensive Ab Initio Study of Halogenated A···U and G···C Base Pair Geometries and Energies

Unraveling the binding preferences involved in the formation of a supramolecular complex is key to properly understand molecular recognition and aggregation phenomena, which are of pivotal importance to biology. The halogenation of nucleic acids has been routinely carried out for decades to assist in their X-ray diffraction analysis. The incorporation of a halogen atom on a DNA/RNA base not only affected its electronic distribution, but also expanded the noncovalent interactions toolbox beyond the classical hydrogen bond (HB) by incorporating the halogen bond (HalB). In this regard, an inspection of the Protein Data Bank (PDB) revealed 187 structures involving halogenated nucleic acids (either unbound or bound to a protein) where at least 1 base pair (BP) exhibited halogenation. Herein, we were interested in disclosing the strength and binding preferences of halogenated A···U and G···C BPs, which are predominant in halogenated nucleic acids. To achieve that, computations at the RI-MP2/def2-TZVP level of theory together with state of the art theoretical modeling tools (including the computation of molecular electrostatic potential (MEP) surfaces, the quantum theory of “Atoms in Molecules” (QTAIM) and the non-covalent interactions plot (NCIplot) analyses) allowed for the characterization of the HB and HalB complexes studied herein.

Aggregation-Induced Electrochemiluminescence and Nitric Oxide Recognition by Halogen Bonding with a Ruthenium(II) Complex

In this study, a new strategy for NO detection based on the aggregation-induced electrochemical luminescence (AIECL) of a ruthenium-based complex and the halogen bonding effect was developed. First, [Ru(phen)₂ (phen-Br₂ )]²⁺ (phen : 1,10-phenanthroline, phen-Br₂  : 3,8-dibromo-1,10-phenanthroline) was synthesized and exhibited aggregation-induced emission (AIE) and AIECL properties in a poor solvent (H₂ O). [Ru(phen)₂ (phen-Br₂ )]²⁺ exhibited greatly enhanced AIECL properties compared to its AIE intensity. When the volume fraction of water (f_w , v %) in the H₂ O-acetonitrile (MeCN) system was increased from 30 to 90 %, the photoluminescence and electrochemiluminescence (ECL) intensities were three- and 800-fold that of the pure MeCN system, respectively. Dynamic light scattering and scanning electron microscopy results indicated that [Ru(phen)₂ (phen-Br₂ )]²⁺ aggregated into nanoparticles. AIECL is sensitive to NO because of its halogen bonding effect. The C-Br⋅⋅⋅N bond between [Ru(phen)₂ (phen-Br₂ )]²⁺ and NO increased the distance of complex molecules, resulting in ECL quenching. A detection limit of 2 nM was obtained with 5 orders of magnitude linear range. The combination of the AIECL system and the halogen bond effect expands the theoretical research and applications in biomolecular detection, molecular sensors, and stages of medical diagnosis.

Exploring Orthogonality between Halogen and Hydrogen Bonding Involving Benzene

The concept of orthogonality between halogen and hydrogen bonding, brought out by Ho and coworkers some years ago, has become a widely accepted idea within the chemists' community. While the original work was based on a common carbonyl oxygen as acceptor for both interactions, we explore here, by means of M06-2X, M11, ωB97X, and ωB97XD/aug-cc-PVTZ DFT calculations, the interdependence of halogen and hydrogen bonding with a shared π-electron system of benzene. The donor groups (specifically NCBr and H2O) were placed on either or the same side of the ring, according to a double T-shaped or a perpendicular geometry, respectively. The results demonstrate that the two interactions with benzene are not strictly independent on each other, therefore outlining that the orthogonality between halogen and hydrogen bonding, intended as energetical independence between the two interactions, should be carefully evaluated according to the specific acceptor group.

Group 14 Central Atoms and Halogen Bonding in Different Dielectric Environments: How Germanium Outperforms Silicon

The nature of halogen bonding under different dielectric conditions remains underexplored, especially for inorganic systems. The structural and energetic properties of model halogen bonded complexes (R₃ M-I-NH₃ for R=H and F, and M=C, Si, and Ge) are examined computationally for relative permittivities between 1 and 109 using an implicit solvent model. We confirm and assess the exceptionally high maximum potentials at the sigma hole on I (V_s,max ) in F₃ Ge-I relative to cases where M=C or Si. In particular, Ge far outperforms Si in mediating inductive effects. Linear relationships, typically with R² >0.97, are identified between V_s,max , the full point charge on I in R₃ M-I, and the interaction energy, and optimized I-N distance in the complexes. An anomalous trend is identified in which, for each M, F₃ M-I-NH₃ becomes less stable as the optimized I-N distance gets shorter in different dielectric environments; it is explained using the F-I-NH₃ complex as a reference.

A Reduced Generalized Force Field for Biological Halogen Bonds

The halogen bond (or X-bond) is a noncovalent interaction that is increasingly recognized as an important design tool for engineering protein-ligand interactions and controlling the structures of proteins and nucleic acids. In the past decade, there have been significant efforts to characterize the structure-energy relationships of this interaction in macromolecules. Progress in the computational modeling of X-bonds in biological molecules, however, has lagged behind these experimental studies, with most molecular mechanics/dynamics-based simulation methods not properly treating the properties of the X-bond. We had previously derived a force field for biological X-bonds (ffBXB) based on a set of potential energy functions that describe the anisotropic electrostatic and shape properties of halogens participating in X-bonds. Although fairly accurate for reproducing the energies within biomolecular systems, including X-bonds engineered into a DNA junction, the ffBXB with its seven variable parameters was considered to be too unwieldy for general applications. In the current study, we have generalized the ffBXB by reducing the number of variables to just one for each halogen type and show that this remaining electrostatic variable can be estimated for any new halogenated molecule through a standard restricted electrostatic potential calculation of atomic charges. In addition, we have generalized the ffBXB for both nucleic acids and proteins. As a proof of principle, we have parameterized this reduced and more general ffBXB against the AMBER force field. The resulting parameter set was shown to accurately recapitulate the quantum mechanical landscape and experimental interaction energies of X-bonds incorporated into DNA junction and T4 lysozyme model systems. Thus, this reduced and generalized ffBXB is more readily adaptable for incorporation into classical molecular mechanics/dynamics algorithms, including those commonly used to design inhibitors against therapeutic targets in medicinal chemistry and materials in biomolecular engineering.

Enthalpy-Entropy Compensation in Biomolecular Recognition: A Computational Perspective

This mini-review provides an overview of the enthalpy-entropy compensation phenomenon in the simulation of biomacromolecular recognition, with particular emphasis on ligand binding. We approach this complex phenomenon from the point of view of practical computational chemistry. Without providing a detailed description of the plethora of existing methodologies already reviewed in depth elsewhere, we present a series of examples to illustrate different approaches to interpret and predict compensation phenomena at an atomistic level, which is far from trivial to predict using canonical, classic textbook assumptions.

A Biological Take on Halogen Bonding and Other Non-Classical Non-Covalent Interactions

Classical hydrogen bonds have, for many decades, been the dominant non-covalent interaction in the toolbox that chemists and chemical engineers have used to design and control the structures of compounds and molecular assemblies as novel materials. Recently, a set of non-classical non-covalent (NC-NC) interactions have emerged that exploit the properties of the Group IV, V, VI, and VII elements of the periodic table (the tetrel, pnictogen, chalcogen, and halogen bonds, respectively). Our research group has been characterizing the prevalence, geometric constraints, and structure-function relationship specifically of the halogen bond in biological systems. We have been particularly interested in exploiting the biological halogen bonds (or BXBs) to control the structures, stabilities, and activities of biomolecules, including the DNA Holliday junction and enzymes. In this review, we first provide a set of criteria for how to determine whether BXBs or any other NC-NC interactions would have biological relevance. We then navigate the trail of studies that had led us from an initial, very biological question to our current point in the journey to establish BXBs as a tool for biomolecular engineering. Finally, we close with a perspective on future directions for this line of research.

Quantifying Intramolecular Halogen Bonds in Nucleic Acids: A Combined Protein Data Bank and Theoretical Study

In this study, we report experimental (Protein Data Bank (PDB) search) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the nature, stability, and directionality properties of intramolecular halogen bonding interactions (HaBs) between 5-bromo/5-iodoracil bases and backbone phosphate groups in nucleic acids (NAs). A PDB survey revealed relevant examples where intramolecular HaBs are undertaken and serve as a structural source of stability in RNA and DNA molecules. In order to develop suitable energy predictors, we started this investigation by calculating the interaction energy values and both the potential V(r) and kinetic G(r) energy densities (using Bader's "atoms in molecules" theory) of several halogen bond complexes involving 5-bromo/5-iodoracil molecules and biologically relevant electron donors. Once the energy predictors based on V(r)/G(r) energy densities were developed, we analyzed the HaBs observed in the biological examples retrieved from the PDB search in order to estimate the strength of the interaction. As far as our knowledge extends, intramolecular halogen bonds in NAs have not been previously quantified in the literature using this methodology and may be of great importance in understanding their structural properties as well as in the construction of molecular materials with DNA and other biological macromolecules.

Halogen Bonds in Protein Nucleic Acid Recognition

Identifying new binding forces between electron donor and acceptor entities is key to properly understanding molecular recognition and aggregation phenomena, which are of inmense importance to biology. For decades, the halogenation of DNA/RNA bases has been routinely carried out to solve solid state structures of nucleic acids (NA). However, the effects of this modification might be deeper than just a simple atom substitution since halogens are also known to undergo noncovalent binding (halogen bonding). Herein we show that halogenated NAs with either Br or I atoms are able to establish halogen bonds with properly disposed protein residues. An inspection of the Protein Data Bank (PDB) reveals several examples involving 5-iodo/5-bromouracil, 8-bromoadenine, and 5-iodocytosine bases that are consistent with the halogen bond geometry features. Computations reveal the favorable and moderately strong nature of this interaction, thus confirming the ability of halogenated bases to actively participate in protein-NA binding.

Cyclic Automated Model Building (CAB) Applied to Nucleic Acids

Obtaining high-quality models for nucleic acid structures by automated model building programs (AMB) is still a challenge. The main reasons are the rather low resolution of the diffraction data and the large number of rotatable bonds in the main chains. The application of the most popular and documented AMB programs (e.g., PHENIX.AUTOBUILD, NAUTILUS and ARP/wARP) may provide a good assessment of the state of the art. Quite recently, a cyclic automated model building (CAB) package was described; it is a new AMB approach that makes the use of BUCCANEER for protein model building cyclic without modifying its basic algorithms. The applications showed that CAB improves the efficiency of BUCCANEER. The success suggested an extension of CAB to nucleic acids—in particular, to check if cyclically including NAUTILUS in CAB may improve its effectiveness. To accomplish this task, CAB algorithms designed for protein model building were modified to adapt them to the nucleic acid crystallochemistry. CAB was tested using 29 nucleic acids (DNA and RNA fragments). The phase estimates obtained via molecular replacement (MR) techniques were automatically submitted to phase refinement and then used as input for CAB. The experimental results from CAB were compared with those obtained by NAUTILUS, ARP/wARP and PHENIX.AUTOBUILD.

Hydrogen Bond Enhanced Halogen Bonds: A Synergistic Interaction in Chemistry and Biochemistry

The halogen bond (XB) has become an important tool for molecular design in all areas of chemistry, including crystal and materials engineering and medicinal chemistry. Its similarity to the hydrogen bond (HB) makes the relationship between these interactions complex, at times competing against and other times orthogonal to each other. Recently, our two laboratories have independently reported and characterized a synergistic relationship, in which the XB is enhanced through direct intramolecular HBing to the electron-rich belt of the halogen. In one study, intramolecular HBing from an amine polarizes the iodopyridinium XB donors of a bidentate anion receptor. The resulting HB enhanced XB (or HBeXB) preorganizes and further augments the XB donors. Consequently, the affinity of the receptor for halogen anions was significantly increased. In a parallel study, a meta-chlorotyrosine was engineered into T4 lysozyme, resulting in a HBeXB that increased the thermal stability and activity of the enzyme at elevated temperatures. The crystal structure showed that the chlorine of the noncanonical amino acid formed a XB to the protein backbone, which augmented the HB of the wild-type enzyme. Calorimetric analysis resulted in an enthalpic contribution of this Cl-XB to the stability of the protein that was an order of magnitude greater than previously determined in biomolecules. Quantum mechanical (QM) calculations showed that rotating the hydroxyl group of the tyrosine to point toward rather than away from the halogen greatly increased its potential to serve as a XB donor, equivalent to what was observed experimentally. In sum, the two systems described here show that the HBeXB concept extends the range of interaction energies and geometries to be significantly greater than that of the XB alone. Additionally, surveys of structural databases indicate that the components for this interaction are already present in many existing molecular systems. The confluence of the independent studies from our two laboratories demonstrates the reach of the HBeXB across both chemistry and biochemistry and that intentional engineering of this enhanced interaction will extend the applications of XBs beyond these two initial examples.

Instantaneous Iodine-Assisted DNAzyme Cleavage of Phosphorothioate RNA

Metal ions play a critical role in the RNA-cleavage reaction by interacting with the scissile phosphate and stabilizing the highly negatively charged transition state. Many metal-dependent DNAzymes have been selected for RNA cleavage. Herein, we report that the Ce13d DNAzyme can use nonmetallic iodine (I₂) to cleave a phosphorothioate (PS)-modified substrate. The cleavage yield exceeded 60% for both the R_p and S_p stereoisomers in 10 s, while the yield without the enzyme strand was only ∼10%. The Ce13d cleavage with I₂ also required Na⁺, consistent with the property of Ce13d and confirming the similar role of I₂ as a metal ion. Ce13d had the highest yield among eight tested DNAzymes, with the second highest DNAzyme showing only 20% cleavage. The incomplete cleavage was due to competition from desulfurization and isomerization reactions. This DNAzyme was engineered for fluorescence-based I₂ detection. With EDTA for masking metal ions, I₂ was selectively detected down to 4.7 nM. Oxidation of I^- with Fe³⁺ produced I₂ in situ, allowing detection of Fe³⁺ down to 78 nM. By harnessing nonelectrostatic interactions, such as the I₂/sulfur interaction observed here, more nonmetal species might be discovered to assist DNAzyme-based RNA cleavage.

Biomolecular Simulations of Halogen Bonds with a GROMOS Force Field

Halogen bonds (XBs) are non-covalent interactions in which halogens (X), acting as electrophiles, interact with Lewis bases. XBs are able to mediate protein-ligand recognition and therefore play an important role in rational drug design. In this context, the development of molecular modeling tools that can tackle XBs is paramount. XBs are predominantly explained by the existence of a positive region on the electrostatic potential of X named the σ-hole. Typically, with molecular mechanics force fields, this region is modeled using a charged extra point (EP) linked to X along the R-X covalent bond axis. In this work, we developed the first EP-based strategy for GROMOS force fields (specifically GROMOS 54A7) using bacteriophage T4 lysozyme in complex with both iodobenzene and iodopentafluorobenzene as a prototype system. Several EP parametrization schemes were tested by adding a virtual interaction site to ligand topologies retrieved from the Automated Topology Builder (ATB) and Repository. Contrary to previous approaches using other force fields, our analysis is based on the capability of each parametrization scheme to sample XBs during MD simulations. Our results indicate that the implementation of an EP at a distance from iodine corresponding to R_min provides a good qualitative description of XBs in MD simulations, supporting the compatibility of our approach with the GROMOS 54A7 force field.

Increasing Enzyme Stability and Activity through Hydrogen Bond-Enhanced Halogen Bonds

The construction of more stable proteins is important in biomolecular engineering, particularly in the design of biologics-based therapeutics. We show here that replacing the tyrosine at position 18 (Y18) of T4 lysozyme with the unnatural amino acid m-chlorotyrosine (^mClY) increases both the thermal stability (increasing the melting temperature by ∼1 °C and the melting enthalpy by 3 kcal/mol) and the enzymatic activity at elevated temperatures (15% higher than that of the parent enzyme at 40 °C) of this classic enzyme. The chlorine of ^mClY forms a halogen bond (XB) to the carbonyl oxygen of the peptide bond at glycine 28 (G28) in a tight loop near the active site. In this case, the XB potential of the typically weak XB donor Cl is shown from quantum chemical calculations to be significantly enhanced by polarization via an intramolecular hydrogen bond (HB) from the adjacent hydroxyl substituent of the tyrosyl side chain, resulting in a distinctive synergistic HB-enhanced XB (or HeX-B for short) interaction. The larger halogens (bromine and iodine) are not well accommodated within this same loop and, consequently, do not exhibit the effects on protein stability or function associated with the HeX-B interaction. Thus, we have for the first time demonstrated that an XB can be engineered to stabilize and increase the activity of an enzyme, with the increased stabilizing potential of the HeX-B further extending the application of halogenated amino acids in the design of more stable protein therapeutics.

The Molecular Origin of Enthalpy/Entropy Compensation in Biomolecular Recognition

Biomolecular recognition can be stubborn; changes in the structures of associating molecules, or the environments in which they associate, often yield compensating changes in enthalpies and entropies of binding and no net change in affinities. This phenomenon—termed enthalpy/entropy (H/S) compensation—hinders efforts in biomolecular design, and its incidence—often a surprise to experimentalists—makes interactions between biomolecules difficult to predict. Although characterizing H/S compensation requires experimental care, it is unquestionably a real phenomenon that has, from an engineering perspective, useful physical origins. Studying H/S compensation can help illuminate the still-murky roles of water and dynamics in biomolecular recognition and self-assembly. This review summarizes known sources of H/ S compensation (real and perceived) and lays out a conceptual framework for understanding and dissecting—and, perhaps, avoiding or exploiting—this phenomenon in biophysical systems.

Aggregation-Induced Emission-Active Ruthenium(II) Complex of 4,7-Dichloro Phenanthroline for Selective Luminescent Detection and Ribosomal RNA Imaging

The development of red emissive aggregation-induced emission (AIE) active probes for organelle-specific imaging is of great importance. Construction of metal complex-based AIE-active materials with metal-to-ligand charge transfer (MLCT), ligand-to-metal charge transfer (LMCT) emission together with the ligand-centered and intraligand (LC/ILCT) emission is a challenging task. We developed a red emissive ruthenium(II) complex, 1[PF₆]₂, and its perchlorate analogues of the 4,7-dichloro phenanthroline ligand. 1[PF₆]₂ has been characterized by spectroscopic and single-crystal X-ray diffraction. Complex 1 showed AIE enhancement in water, highly dense polyethylene glycol media, and also in the solid state. The possible reason behind the AIE property may be the weak supramolecular π···π, C-H···π, and C-Cl···H interactions between neighboring phen ligands as well as C-Cl···O halogen bonding (XB). The crystal structures of the two perchlorate analogues revealed C-Cl···O distances shorter than the sum of the van der Waals radii, which confirmed the XB interaction. The AIE property was supported by scanning electron microscopy, transmission electron microscopy, dynamic light scattering, and atomic force microscopy studies. Most importantly, the probe was found to be low cytotoxicity and to efficiently permeate the cell membrane. The cell-imaging experiments revealed rapid staining of the nucleolus in HeLa cells via the interaction with nucleolar ribosomal ribonucleic acid (rRNA). It is expected that the supramolecular interactions as well as C-Cl···O XB interaction with rRNA is the origin of aggregation and possible photoluminescence enhancement. To the best of our knowledge, this is the first report of red emissive ruthenium(II) complex-based probes with AIE characteristics for selective rRNA detection and nucleolar imaging.

Structural Examination of Halogen-Bonded Co-Crystals of Tritopic Acceptors

A series of tritopic N-heterocyclic compounds containing electrostatically and geometrically equivalent binding sites were synthesized and subjected to systematic co-crystallizations with selected perfluoroiodoarenes in order to map out their structural landscapes. More than 70% of the attempted reactions produced a co-crystal as indicated by IR spectroscopy. Four new crystal structures are reported and in all of them, at least one potential binding site on the acceptor is left vacant. The absence of halogen bonds to all sites can be ascribed primarily due to deactivation of the σ-hole on the iodo-arene donors and partially due to steric hindrance. The tritopic acceptors containing 5,6-dimethylbenzimidazole derivatives yield discrete tetrameric aggregates in the solid state, whereas the pyrazole and imidazole analogues assemble into halogen-bonded 1-D chains.

Structure-based domain assignment in Leishmania infantum EndoG: characterization of a pH-dependent regulatory switch and a C-terminal extension that largely dictates DNA substrate preferences

Mitochondrial endonuclease G from Leishmania infantum (LiEndoG) participates in the degradation of double-stranded DNA (dsDNA) during parasite cell death and is catalytically inactive at a pH of 8.0 or above. The presence, in the primary sequence, of an acidic amino acid-rich insertion exclusive to trypanosomatids and its spatial position in a homology-built model of LiEndoG led us to postulate that this peptide stretch might act as a pH sensor for self-inhibition. We found that a LiEndoG variant lacking residues 145–180 is indeed far more active than its wild-type counterpart at pH values >7.0. In addition, we discovered that (i) LiEndoG exists as a homodimer, (ii) replacement of Ser211 in the active-site SRGH motif with the canonical aspartate from the DRGH motif of other nucleases leads to a catalytically deficient enzyme, (iii) the activity of the S211D variant can be restored upon the concomitant replacement of Ala247 with Arg and (iv) a C-terminal extension is responsible for the observed preferential cleavage of single-stranded DNA (ssDNA) and ssDNA–dsDNA junctions. Taken together, our results support the view that LiEndoG is a multidomain molecular machine whose nuclease activity can be subtly modulated or even abrogated through architectural changes brought about by environmental conditions and interaction with other binding partners.

Sulfur as an Acceptor to Bromine in Biomolecular Halogen Bonds

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 5-bromouracil 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 ΔΔH^25°C and -0.4 kcal/mol ΔΔG^25°C (including an entropic penalty ΔΔS^25°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.

Structure of the Holliday junction: applications beyond recombination

The Holliday junction (HJ) is an essential element in recombination and related mechanisms. The structure of this four-stranded DNA assembly, which is now well-defined alone and in complex with proteins, has led to its applications in areas well outside of molecular recombination, including nanotechnology and biophysics. This minireview explores some interesting recent research on the HJ, as it has been adapted to design regular two- or three-dimensional lattices for crystal engineering, and more complex systems through DNA origami. In addition, the sequence dependence of the structure is discussed in terms how it can be applied to characterize the geometries and energies of various noncovalent interactions, including halogen bonds in oxidatively damaged (halogenated) bases and hydrogen bonds associated with the epigenetic 5-hydroxylmethylcytosine base.

Halogen Bonding in Nucleic Acid Complexes

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.

Halogen-π Interactions in the Cytochrome P450 Active Site: Structural Insights into Human CYP2B6 Substrate Selectivity

Numerous cytochrome P450 (CYP) 2B6 substrates including drugs and environmental chemicals are halogenated. To assess the role of halogen-π bonds in substrate selectivity and orientation in the active site, structures of four CYP2B6 monoterpenoid complexes were solved by X-ray crystallography. Bornyl bromide exhibited dual orientations in the active site with the predominant orientation revealing a bromine-π bond with the Phe108 side chain. Bornane demonstrated two orientations with equal occupancy; in both, the C2 atom that bears the bromine in bornyl bromide was displaced by more than 2.5 Å compared with the latter complex. The bromine in myrtenyl bromide π-bonded with Phe297 in CYP2B6, whereas the two major orientations in the active site mutant I114V exhibited bromine-π interactions with two additional residues, Phe108 and Phe115. Analysis of existing structures suggests that halogen-π interactions may be unique to the CYP2B enzymes within CYP family 2 but are also important for CYP3A enzymes.

Relationships between hydrogen bonds and halogen bonds in biological systems

The recent recognition that halogen bonding (XB) plays important roles in the recognition and assembly of biological molecules has led to new approaches in medicinal chemistry and biomolecular engineering. When designing XBs into strategies for rational drug design or into a biomolecule to affect its structure and function, we must consider the relationship between this interaction and the more ubiquitous hydrogen bond (HB). In this review, we explore these relationships by asking whether and how XBs can replace, compete against or behave independently of HBs in various biological systems. The complex relationships between the two interactions inform us of the challenges we face in fully utilizing XBs to control the affinity and recognition of inhibitors against their therapeutic targets, and to control the structure and function of proteins, nucleic acids and other biomolecular scaffolds.

Halogen-bonding contacts determining the crystal structure and fluorescence properties of organic salts

Different inorganic anions tuning halogen-bonding contacts to form different 3D networks with various absorption and emission properties.

Effect of Hydroxymethylcytosine on the Structure and Stability of Holliday Junctions

5-Hydroxymethylcytosine (5hmC) is an epigenetic marker that has recently been shown to promote homologous recombination (HR). In this study, we determine the effects of 5hmC on the structure, thermodynamics, and conformational dynamics of the Holliday junction (the four-stranded DNA intermediate associated with HR) in its native stacked-X form. The hydroxymethyl and the control methyl substituents are placed in the context of an amphimorphic GxCC trinucleotide core sequence (where xC is C, 5hmC, or the methylated 5mC), which is part of a sequence also recognized by endonuclease G to promote HR. The hydroxymethyl group of the 5hmC junction adopts two distinct rotational conformations, with an in-base-plane form being dominant over the competing out-of-plane rotamer that has typically been seen in duplex structures. The in-plane rotamer is seen to be stabilized by a more stable intramolecular hydrogen bond to the junction backbone. Stabilizing hydrogen bonds (H-bonds) formed by the hydroxyl substituent in 5hmC or from a bridging water in the 5mC structure provide approximately 1.5-2 kcal/mol per interaction of stability to the junction, which is mostly offset by entropy compensation, thereby leaving the overall stability of the G5hmCC and G5mCC constructs similar to that of the GCC core. Thus, both methyl and hydroxymethyl modifications are accommodated without disrupting the structure or stability of the Holliday junction. Both 5hmC and 5mC are shown to open the structure to make the junction core more accessible. The overall consequences of incorporating 5hmC into a DNA junction are thus discussed in the context of the specificity in protein recognition of the hydroxymethyl substituent through direct and indirect readout mechanisms.

Insights into halogen bond-driven enantioseparations

Although the halogen bond (XB) has been so far mainly studied in silico and in the solid state, its potential impact in solution is yet to be fully understood. In this study, we describe the first systematic investigation on the halogen bond in solvated environment by high-performance liquid chromatography (HPLC). Thirty three atropisomeric polyhalogenated-4,4'-bipyridines (HBipys), containing Cl, Br and I as substituents, were selected and used as potential XB donors (XBDs) on two cellulose-based chiral stationary phases (CSPs) containing potential XB acceptors (XBAs). The impact of the halogens on the enantiodiscrimination mechanism was investigated and iodine showed a pivotal role on the enantioseparation in non-polar medium. Electrostatic potentials (EPs) were computed to understand the electrostatic component of CSP-analyte interaction. Moreover, van't Hoff studies for ten HBipys were performed and the thermodynamic parameters governing the halogen-dependent enantioseparations are discussed. Finally, a molecular dynamic (MD) simulation is proposed to model halogen bond in polysaccharide-analyte complexes by inclusion of a charged extra point to represent the positive 'σ-hole' on the halogen atom. On the basis of both experimental results and theoretical data, we have profiled the halogen bond as a chemo-, regio-, site- and stereoselective interaction which can work in HPLC environment besides other known interactions based on the complementarity between selector and selectand.

The Halogen Bond

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

A PGSE NMR approach to the characterization of single and multi-site halogen-bonded adducts in solution

The diffusional PGSE NMR technique allows a detailed characterization of single- and multi-site halogen-bonded systems.

Evaluating the potential of halogen bonding in molecular design: automated scaffold decoration using the new scoring function XBScore

We present a QM-derived empirical scoring function for the interaction between aromatic halogenated ligands and the carbonyl oxygen atom of the protein backbone. Applying this scoring function, we developed an algorithm that evaluates the potential of protein-bound ligands to form favorable halogen-bonding contacts upon scaffold decoration with chlorine, bromine, or iodine. Full recovery of all existing halogen bonds in the PDB involving the protein backbone was achieved with our protocol. Interestingly, the potential for introducing halogen bonds through scaffold decoration of unsubstituted aromatic carbon atoms appears to easily match the number of previously known halogen bonds. Our approach can thus be used as a blueprint for integration of halogen bonding into general empirical scoring functions, which at present ignore this interaction. Most importantly, we were able to identify a substantial number of protein-ligand complexes where the benefits and challenges of introducing a halogen bond by molecular design can be studied experimentally.

Force Field Model of Periodic Trends in Biomolecular Halogen Bonds

The study of the noncovalent interaction now defined as a halogen bond (X-bond) has become one of the fastest growing areas in experimental and theoretical chemistry--its applications as a design tool are highly extensive. The significance of the interaction in biology has only recently been recognized, but has now become important in medicinal chemistry. We had previously derived a set of empirical potential energy functions to model the structure-energy relationships for bromines in biomolecular X-bonds (BXBs). Here, we have extended this force field for BXBs (ffBXB) to the halogens (Cl, Br, and I) that are commonly seen to form stable X-bonds. The ffBXB calculated energies show a remarkable one-to-one linear relationship to explicit BXB energies determined from an experimental DNA junction system, thereby validating the approach and the model. The resulting parameters allow us to interpret the stabilizing effects of BXBs in terms of well-defined physical properties of the halogen atoms, including their size, shape, and charge, showing periodic trends that are predictable along the Group VII column of elements. Consequently, we have established the ffBXB as an accurate computational tool that can be applied, for example, for the design of new therapeutic compounds against clinically important targets and new biomolecular-based materials.

Halogen- and hydrogen-bonded salts and co-crystals formed from 4-halo-2,3,5,6-tetrafluorophenol and cyclic secondary and tertiary amines: orthogonal and non-orthogonal halogen and hydrogen bonding, and synthetic analogues of halogen-bonded biological systems

Co-crystallisation of, in particular, 4-iodotetrafluorophenol with a series of secondary and tertiary cyclic amines results in deprotonation of the phenol and formation of the corresponding ammonium phenate. Careful examination of the X-ray single-crystal structures shows that the phenate anion develops a C=O double bond and that the C-C bond lengths in the ring suggest a Meissenheimer-like delocalisation. This delocalisation is supported by the geometry of the phenate anion optimised at the MP2(Full) level of theory within the aug-cc-pVDZ basis (aug-cc-pVDZ-PP on I) and by natural bond orbital (NBO) analyses. With sp(2) hybridisation at the phenate oxygen atom, there is strong preference for the formation of two non-covalent interactions with the oxygen sp(2) lone pairs and, in the case of secondary amines, this occurs through hydrogen bonding to the ammonium hydrogen atoms. However, where tertiary amines are concerned, there are insufficient hydrogen atoms available and so an electrophilic iodine atom from a neighbouring 4-iodotetrafluorophenate group forms an I⋅⋅⋅O halogen bond to give the second interaction. However, in some co-crystals with secondary amines, it is also found that in addition to the two hydrogen bonds forming with the phenate oxygen sp(2) lone pairs, there is an additional intermolecular I⋅⋅⋅O halogen bond in which the electrophilic iodine atom interacts with the C=O π-system. All attempts to reproduce this behaviour with 4-bromotetrafluorophenol were unsuccessful. These structural motifs are significant as they reproduce extremely well, in low-molar-mass synthetic systems, motifs found by Ho and co-workers when examining halogen-bonding interactions in biological systems. The analogy is cemented through the structures of co-crystals of 1,4-diiodotetrafluorobenzene with acetamide and with N-methylbenzamide, which, as designed models, demonstrate the orthogonality of hydrogen and halogen bonding proposed in Ho's biological study.

Thermodynamic and structural characterization of halogen bonding in protein-ligand interactions: a case study of PDE5 and its inhibitors

The significance of halogen bonding in protein-ligand interactions has been recognized recently. We present here the first comprehensive thermodynamic and structural characterization of halogen bonding in PDE5-inhibitor interactions. ITC studies reveal that binding strength of the halogen bonding between chlorine, bromine, and iodine of inhibitor and the protein is -1.57, -3.09, and -5.59 kJ/mol, respectively. The halogens interact with the designed residue Y612 and an unexpected buried water molecule.

Halogen bonding controls selectivity of FRET substrate probes for MMP-9

Matrix metalloproteinases (MMPs) are a large family of zinc-dependent endoproteases that catalyze cleavage of extracellular matrix and nonmatrix proteins. MMPs play a role in tissue remodeling, and their uncontrolled activity is associated with number of diseases, including tumor metastasis. Thus, there is a need to develop methods to monitor MMP activity, and number of probes has been previously described. The key problem many probes encounter is the issue of selectivity, since 23 human MMPs, despite playing different physiological roles, have structurally similar active sites. Here, we introduce the halogen bonding concept into the probe design and show that the probe containing iodine exhibits an unprecedented selectivity for MMP-9. We provide structure-based explanation for the selectivity, confirming that it is due to formation of the halogen bond that supports catalysis, and we highlight the value of exploring halogen bonding in the context of selective probe design.

Probing halogen bonds with solid-state NMR spectroscopy: observation and interpretation of J(77Se,31P) coupling in halogen-bonded PSe⋯I motifs

Intra-halogen bond J couplings measured via NMR spectroscopy and interpreted using natural localized molecular orbitals offer novel insights into this class of non-covalent interaction.