Kevin E. Riley

Directional Weak Intermolecular Interactions: σ-Hole Bonding

The prototypical directional weak interactions, hydrogen bonding and σ-hole bonding (including the special case of halogen bonding) are reviewed in a united picture that depends on the anisotropic nature of the molecular electrostatic potential around the donor atom. Qualitative descriptions of the effects that lead to these anisotropic distributions are given and examples of the importance of σ-hole bonding in crystal engineering and biological systems are discussed.

Noncovalent interactions in biochemistry

Noncovalent interactions are known to play a key role in biochemistry. The knowledge of stabilization (relative) energies and their components is very important for understanding the nature of these interactions. Accurate and benchmark data on interaction (relative) energies and structures can be obtained from coupled‐cluster with single and double and perturbative triple excitations [CCSD(T)] calculations with extended basis of atomic orbitals or even at the complete basis set limit. These methods cannot be, however, used for systems larger than about 50 atoms. In this contribution, the applicability and performance of various recently introduced wavefunction and density functional methods are examined in detail. It is shown that a very good performance by some of these methods is obtained only by introducing empirical parameters fitted mostly to CCSD(T) benchmark data. Among the methods described, special attention is paid to two techniques. First, the symmetry‐adapted perturbation technique that allows obtaining not only accurate values of total interaction energies but also their components. Results of these calculations reveal a key role of dispersion energy in stabilizing the structures of biomolecular systems. Second, the semiempirical quantum chemical parameterized model 6 method (PM6) augmented by empirical terms describing the dispersion and H‐bonding energies. The method is suitable for much extended systems having several thousands of atoms and can be thus used, e.g., in the drug design.

The relative roles of electrostatics and dispersion in the stabilization of halogen bonds

In this work we highlight recent work aimed at the characterization of halogen bonds. Here we discuss the origins of the σ-hole, the modulation of halogen bond strength by changing of neighboring chemical groups (i.e. halogen bond tuning), the performance of various computational methods in treating halogen bonds, and the strength and character of the halogen bond, the dihalogen bond, and two hydrogen bonds in bromomethanol dimers (which serve as model complexes) are compared. Symmetry adapted perturbation theory analysis of halogen bonding complexes indicates that halogen bonds strongly depend on both dispersion and electrostatics. The electrostatic interaction that occurs between the halogen σ-hole and the electronegative halogen bond donor is responsible for the high degree of directionality exhibited by halogen bonds. Because these noncovalent interactions have a strong dispersion component, it is important that the computational method used to treat a halogen bonding system be chosen very carefully, with correlated methods (such as CCSD(T)) being optimal. It is also noted here that most forcefield-based molecular mechanics methods do not describe the halogen σ-hole, and thus are not suitable for treating systems with halogen bonds. Recent attempts to improve the molecular mechanics description of halogen bonds are also discussed.

MP2.5 and MP2.X: approaching CCSD(T) quality description of noncovalent interaction at the cost of a single CCSD iteration.

The performance of the second-order Møller-Plesset perturbation theory MP2.5 and MP2.X methods, tested on the S22, S66, X40, and other benchmark datasets is briefly reviewed. It is found that both methods produce highly accurate binding energies for the complexes contained in these data sets. Both methods also provide reliable potential energy curves for the complexes in the S66 set. Among the routinely used wavefunction methods, the only other technique that consistently produces lower errors, both for stabilization energies and geometry scans, is the spin-component-scaled coupled-clusters method covering iterative single- and double-electron excitations, which is, however, substantially more computationally intensive. The structures originated from full geometrical gradient optimizations at the MP2.5 and MP2.X level of theory were confirmed to be the closest to the CCSD(T)/CBS (coupled clusters covering iterative single- and double-electron excitations and perturbative triple-electron excitations performed at the complete basis set limit) geometries among all the tested methods (e.g. MP3, SCS(MI)-MP2, MP2, M06-2X, and DFT-D method evaluated with the TPSS functional). The MP2.5 geometries for the tested complexes deviate from the references almost negligibly. Inclusion of the scaled third-order correlation energy results in a substantial improvement of the ability to accurately describe noncovalent interactions. The results shown here serve to support the notion that MP2.5 and MP2.X are reasonable alternative methods for benchmark calculations in cases where system size or (lack of) computational resources preclude the use of CCSD(T)/CBS computations. MP2.X allows for the use of smaller basis sets (i.e. 6-31G*) with results that are nearly identical to those of MP2.5 with larger basis sets, which dramatically decreases computation times and makes calculations on much larger systems possible.

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

AbstractO-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7xa0kcal mol-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3xa0kcal mol-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.n FigureCompetition of hydrogen, halogen, and dihalogen bonding in the bromomethanol dimer are investigated

Ligands of Therapeutic Utility for the Liver X Receptors

Liver X receptors (LXRs) have been increasingly recognized as a potential therapeutic target to treat pathological conditions ranging from vascular and metabolic diseases, neurological degeneration, to cancers that are driven by lipid metabolism. Amidst intensifying efforts to discover ligands that act through LXRs to achieve the sought-after pharmacological outcomes, several lead compounds are already being tested in clinical trials for a variety of disease interventions. While more potent and selective LXR ligands continue to emerge from screening of small molecule libraries, rational design, and empirical medicinal chemistry approaches, challenges remain in minimizing undesirable effects of LXR activation on lipid metabolism. This review provides a summary of known endogenous, naturally occurring, and synthetic ligands. The review also offers considerations from a molecular modeling perspective with which to design more specific LXRβ ligands based on the interaction energies of ligands and the important amino acid residues in the LXRβ ligand binding domain.

Exploring the (Very Flat) Potential Energy Landscape of R−Br⋅⋅⋅π Interactions with Accurate CCSD(T) and SAPT Techniques

Halogen bonds involving an aromatic moiety as an acceptor, otherwise known as R-X⋅⋅⋅π interactions, have increasingly been recognized as being important in materials and in protein-ligand complexes. These types of interactions have been the subject of many recent investigations, but little is known about the ways in which the strengths of R-X⋅⋅⋅π interactions vary as a function of the relative geometries of the interacting pairs. Here we use the accurate CCSD(T) and SAPT2+3δMP2 methods to investigate the potential energy landscapes for systems of HBr, HCCBr, and NCBr complexed with benzene. It is found that only the separation between the complexed molecules have a strong effect on interaction strength while other geometric parameters, such as tilting and shifting R-Br⋅⋅⋅π donor relative to the benzene plane, affect these interactions only mildly. Importantly, it is found that the C6v (T-shaped) configuration is not the global minimum for any of the dimers investigated.

Osteoinductive effects of glyceollins on adult mesenchymal stromal/stem cells from adipose tissue and bone marrow

BACKGROUNDnWhile current therapies for osteoporosis focus on reducing bone resorption, the development of therapies to regenerate bone may also be beneficial. Promising anabolic therapy candidates include phytoestrogens, such as daidzein, which effectively induce osteogenesis of adipose-derived stromal cells (ASCs) and bone marrow stromal cells (BMSCs).nnnPURPOSEnTo investigate the effects of glyceollins, structural derivatives of daidzein, on osteogenesis of ASCs and BMSCs.nnnSTUDY DESIGNnHerein, the osteoinductive effects of glyceollin I and glyceollin II were assessed and compared to estradiol in ASCs and BMSCs. The mechanism by which glyceollin II induces osteogenesis was further examined.nnnMETHODSnThe ability of glyceollins to promote osteogenesis of ASCs and BMSCs was evaluated in adherent and scaffold cultures. Relative deposition of calcium was analyzed using Alizarin Red staining, Bichinchoninic acid Protein Assay, and Alamar Blue Assay. To further explore the mechanism by which glyceollin II exerts its osteoinductive effects, docking studies of glyceollin II, RNA isolation, cDNA synthesis, and quantitative RT-PCR (qPCR) were performed.nnnRESULTSnIn adherent cultures, ASCs and BMSCs treated with estradiol, glyceollin I, or glyceollin II demonstrated increased calcium deposition relative to vehicle-treated cells. During evaluation on PLGA scaffolds seeded with ASCs and BMSCs, glyceollin II was the most efficacious in inducing ASC and BMSC osteogenesis compared to estradiol and glyceollin I. Dose-response analysis in ASCs and BMSCs revealed that glyceollin II has the highest potency at 10nM in adherent cultures and 1µM in tissue scaffold cultures. At all doses, osteoinductive effects were attenuated by fulvestrant, suggesting that glyceollin II acts at least in part through estrogen receptor-mediated pathways to induce osteogenesis. Analysis of gene expression demonstrated that, similar to estradiol, glyceollin II induces upregulation of genes involved in osteogenic differentiation.nnnCONCLUSIONnThe ability of glyceollin II to induce osteogenic differentiation in ASCs and BMSCs indicates that glyceollins hold the potential for the development of pharmacological interventions to improve clinical outcomes of patients with osteoporosis.

Point-of-Care Determination of Acetaminophen Levels with Multi-Hydrogen Bond Manipulated Single-Molecule Recognition (eMuHSiR)

This work aims to face the challenge of monitoring small molecule drugs accurately and rapidly for point-of-care (POC) diagnosis in current clinical settings. Overdose of acetaminophen (AP), a commonly used over the counter (OTC) analgesic drug, has been determined to be a major cause of acute liver failure in the US and the UK. However, there is no rapid and accurate detection method available for this drug in the emergency room. The present study examined an AP sensing strategy that relies on a previously unexplored strong interaction between AP and the arginine (Arg) molecule. It was found that as many as 4 hydrogen bonds can be formed between one Arg molecule and one AP molecule. By taking advantages of this structural selectivity and high tenability of hydrogen bonds, Arg, immobilized on a graphene surface via electrostatic interactions, was utilized to structurally capture AP. Interestingly, bonded AP still remained the perfect electrochemical activities. The extent of Arg-AP bonds was quantified using a newly designed electrochemical (EC) sensor. To verify the feasibility of this novel assay, based on multihydrogen bond manipulated single-molecule recognition (eMuHSiR), both pharmaceutical and serum sample were examined. In commercial tablet measurement, no significant difference was seen between the results of eMuHSiR and other standard methods. For measuring AP concentration in the mice blood, the substances in serum, such as sugars and fats, would not bring any interference to the eMuHSiR in a wide concentration range. This eMuHSiR method opens the way for future development of small molecule detection for the POC testing.

Ortho-Methylarylamines as Time-Dependent Inhibitors of Cytochrome P450 1A1 Enzyme.

BACKGROUNDnMembers of the cytochrome P450 1A family metabolize many procarcinogens such as polycyclicaromatic hydrocarbons and heterocyclic amines. Inactivation of these enzymes is a prerequisite for cancer prevention and treatment in certain cases. Mechanism-based inhibition (time and co-factor dependent) is an effective method for the inactivation of these enzymes. Our recent study on emodin analogs revealed an anthraquinone with ortho-methylarylamine moiety that exhibited timedependent inhibition of P450 enzymes 1A1 and 1A2.nnnMETHODSnTo determine whether the amino group or the methyl group or both were responsible for the time-dependent inhibition of these enzymes, a set of eleven compounds containing the orthomethylarylamine moiety were identified through a database search, and studied for the inhibition of the P450 enzymes 1A1, 1A2, 2A6 and 2B1. Our earlier studies on carbazole derivatives provided us with highly selective P450 1A2 inhibitors. Glycine scanning studies were performed on the docked proteinligand complexes of compounds 1-20 in order to understand the contribution of different protein residues towards the ligand binding.nnnRESULTSnFour compounds were found to cause selective time-dependent inhibition of P450 1A1 with KI values ranging from 0.24 to 8.25 mM. These compounds exhibited only direct inhibition of P450 1A2. Molecular modeling studies of these molecules indicated that the shapes of the molecules, their binding modes, and the methyl substituent in close proximity (4.5-5.7 Å) to the heme-Fe all contributed to their selective time-dependent inhibition activity on P450 1A1. Glycine scanning studies for P450 1A1 indicated that ligand interaction with Phe123 was the strongest binding contributor and similar studies for P450 1A2 indicated that ligand interactions with the phenylalanine residues 226 and 260 were the largest binding contributors.nnnCONCLUSIONnFour compounds have been identified that exhibit selective time-dependent inhibition of P450 1A1. Modeling studies have indicated that the proximity of the aromatic methyl group to the heme-Fe could be the main contributor for time-dependent inhibition. Future studies will focus on the confirmation of the involvement of the aromatic methyl group in enzyme inactivation.