Hydrogen Bonding, Electrostatic Potential, and Molecular Design

Bioisoteres for carboxylic acids: From ionized isosteres to novel unionized replacements

Carboxylic acids are key pharmacophoric elements in many molecules. They can be seen as a problem by some, due to perceived permeability challenges, potential for high plasma protein binding and the risk of forming reactive metabolites due to acyl-glucuronidation. By others they are viewed more favorably as they can decrease lipophilicity by adding an ionizable center which can be beneficial for solubility, and can add enthalpic interactions with the target protein. However, there are many instances where the replacement of a carboxylic acid with a bioisosteric group is required. This has led to the development of a number of ionizable groups which sufficiently mimic the carboxylic acid functionality whilst improving, for example, the metabolic profile of the molecule in question. An alternative strategy involves replacement of the carboxylate by neutral functional groups. This review initially details carefully selected examples whereby tetrazoles, acyl sulfonamides or isoxazolols have been beneficially utilized as carboxylic acid bioisosteres altering physicohemical properties, interactions with the target and metabolism and/or pharmacokinetics, before delving further into the binding mode of carboxylic acid derivatives with their target proteins. This analysis highlights new ways to consider the replacement of carboxylic acids by neutral bioisosteric groups which either rely on hydrogen bonds or cation-π interactions. It should serve as a useful guide for scientists working in drug discovery.

Quantification of Hydrogen-Bond-Donating Ability of Biologically Relevant Compounds

Hydrogen bonding is a key factor in the design of ligands for biological binding, including drug targets. Our group previously developed a method for experimentally assessing the hydrogen-bond-donating ability of an analyte using UV-vis titrations with a colorimetric sensor. Using this method, 79 new titrations were performed on weak hydrogen-bond donors, with a focus on heterocycles and pharmaceutically relevant motifs. The hydrogen-bond donating abilities of drug compounds and the substructures of drug compounds were also measured. These titrations will be used to build a database of hydrogen-bond donors.

Investigations of the reaction mechanism of sodium with hydrogen fluoride to form sodium fluoride and the adsorption of hydrogen fluoride on sodium fluoride monomer and tetramer

CONTEXT

The reaction between Na and HF is a typical harpooning reaction which is of great interest due to its significance in understanding the elementary chemical reaction kinetics. This work aims to investigate the detailed reaction mechanisms of sodium with hydrogen fluoride and the adsorption of HF on the resultant NaF as well as the (NaF)4 tetramer. The results suggest that the reaction between Na and HF leads to the formation of sodium fluoride salt NaF and hydrogen gas. Na interacts with HF to form a complex HF···Na, and then the approaching of F atom of HF to Na results in a transition state H···F···Na. Accompanied by the broken of H-F bond, the bond forms between F and Na atoms as NaF, then the product NaF is yielded due to the removal of H atom. The resultant NaF can further form (NaF)4 tetramer. The interaction of NaF with HF leads to the complex NaF···HF; the form I as well as II of (NaF)4 can interact with HF to produce two complexes (i.e., (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF), but the form III of (NaF)4 can interact with HF to produce only one complex (NaF)4(III)···HF. These complexes were explored in terms of noncovalent interaction (NCI) and quantum theory of atoms in molecules (QTAIM) analyses. NCI analyses confirm the existences of attractive interactions in the complexes HF···Na, NaF···HF, (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF, and (NaF)4(III)···HF. QTAIM analyses suggest that the F···Na interaction forms in the HF···Na complex while the F···H hydrogen bonds form in NaF···HF, (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF, and (NaF)4(III)···HF complexes. Natural bond orbital (NBO) analyses were also applied to analyze the intermolecular donor-acceptor orbital interactions in these complexes. These results would provide valuable insight into the chemical reaction of Na and HF and the adsorption interaction between sodium fluoride salt and HF.

METHODS

The calculations were carried out at the M06-L/6-311++G(2d,2p) level of theory which were performed using the Gaussian16 program. Intrinsic reaction coordinate (IRC) calculations were carried out at the same level of theory to confirm that the obtained transition state was true. The molecular surface electrostatic potential (MSEP) was employed to understand how the complex forms. Quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) analysis was used to know the topology parameters at bond critical points (BCPs) and intermolecular interactions in the complex and intermediate. The topology parameters and the BCP plots were obtained by the Multiwfn software.

Computing the relative binding affinity of ligands based on a pairwise binding comparison network

Structure-based lead optimization is an open challenge in drug discovery, which is still largely driven by hypotheses and depends on the experience of medicinal chemists. Here we propose a pairwise binding comparison network (PBCNet) based on a physics-informed graph attention mechanism, specifically tailored for ranking the relative binding affinity among congeneric ligands. Benchmarking on two held-out sets (provided by Schrödinger and Merck) containing over 460 ligands and 16 targets, PBCNet demonstrated substantial advantages in terms of both prediction accuracy and computational efficiency. Equipped with a fine-tuning operation, the performance of PBCNet reaches that of Schrödinger's FEP+, which is much more computationally intensive and requires substantial expert intervention. A further simulation-based experiment showed that active learning-optimized PBCNet may accelerate lead optimization campaigns by 473%. Finally, for the convenience of users, a web service for PBCNet is established to facilitate complex relative binding affinity prediction through an easy-to-operate graphical interface.

New Thiazole Carboxamide Derivatives as COX Inhibitors: Design, Synthesis, Anticancer Screening, In Silico Molecular Docking, and ADME Profile Studies

The goal of this work was to create and test a new series of thiazole carboxamide derivatives for their cyclooxygenase (COX) suppressor and anticancer effects. The compounds were characterized using 1H, 13C NMR, and HRMS spectrum analysis, and their selectivity toward COX-1 and COX-2 was assessed using an in vitro COX inhibition assay kit. Cytotoxicity was assessed using an MTS assay against a panel of cancer and normal cell lines. The docking studies were aided by the Prime MM-GBSA method for estimating binding affinities. The density functional theory (DFT) analysis was performed to assess compound chemical reactivity, which was calculated by computing the border orbital energy of both HOMO and LUMO orbitals, as well as the HOMO-LUMO energy gap. For ADME-T analysis, the QiKProp module was employed. Furthermore, using human X-ray crystal structures, molecular docking studies were carried out to discover the probable binding patterns of these drugs within both COX-1 and COX-2 isozymes. The results demonstrated that the most effective compound against the COX-1 enzyme was 2b with an IC50 of 0.239 μM. It also showed potent activity against COX-2 with an IC50 value of 0.191 μM and a selectivity ratio of 1.251. The highest selectivity ratio was 2.766 for compound 2a against COX-2 with an IC50 dose of 0.958 μM relating to the celecoxib ratio of 23.8 and its IC50 against COX-2 of 0.002 μM. Compound 2j also showed good selectivity toward COX-2 (1.507) with an IC50 value of 0.957 μM. All compounds showed negligible cytotoxic activity against the evaluated normal cell lines, and the IC50 values were more than 300 μM, except for compound 2b, whose IC50 values were 203.71 ± 1.89 and 116.96 ± 2.05 μM against LX-2 and Hek293t cell lines, respectively. Moreover, compound 2b showed moderate anticancer activity against COLO205 and B16F1 cancer cell lines with IC50 values of 30.79 and 74.15 μM, respectively.

Hydrogen Bonding Parameters by Rapid Colorimetric Assessment: Evaluation of Structural Components Found in Biological Ligands and Organocatalysts

Hydrogen bonding is a key molecular interaction in biological processes, drug delivery, and catalysis. This report describes a high throughput UV-Vis spectroscopic method to measure hydrogen bonding capacity using a pyrazinone sensor. This colormetric sensor reversibly binds to a hydrogen bond donor, resulting in a blue shift as additional equivalents of donor are added. Titration with excess equivalents of donor is used to determine the binding coefficient, ln(Keq ). Over 100 titrations were performed for a variety of biologically relevant compounds. This data enabled development a multiple linear regression model that is capable of predicting 95 % of ln(Keq ) values within 1 unit, allowing for the estimation of hydrogen bonding affinity from a single measurement. To show the effectiveness of the single point measurements, hydrogen bond strengths were obtained for a set of carboxylic acid bioisosteres. The values from the single point measurements were validated with full titrations.

Natural Lipid Extracts as an Artificial Membrane for Drug Permeability Assay: In Vitro and In Silico Characterization

In vitro non-cellular permeability models such as the parallel artificial membrane permeability assay (PAMPA) are widely applied tools for early-phase drug candidate screening. In addition to the commonly used porcine brain polar lipid extract for modeling the blood–brain barrier’s permeability, the total and polar fractions of bovine heart and liver lipid extracts were investigated in the PAMPA model by measuring the permeability of 32 diverse drugs. The zeta potential of the lipid extracts and the net charge of their glycerophospholipid components were also determined. Physicochemical parameters of the 32 compounds were calculated using three independent forms of software (Marvin Sketch, RDKit, and ACD/Percepta). The relationship between the lipid-specific permeabilities and the physicochemical descriptors of the compounds was investigated using linear correlation, Spearman correlation, and PCA analysis. While the results showed only subtle differences between total and polar lipids, permeability through liver lipids highly differed from that of the heart or brain lipid-based models. Correlations between the in silico descriptors (e.g., number of amide bonds, heteroatoms, and aromatic heterocycles, accessible surface area, and H-bond acceptor–donor balance) of drug molecules and permeability values were also found, which provides support for understanding tissue-specific permeability.

Discovery and Optimization of Potent, Selective, and Brain-Penetrant 1-Heteroaryl-1H-Indazole LRRK2 Kinase Inhibitors for the Treatment of Parkinson's Disease

Inhibition of leucine-rich repeat kinase 2 (LRRK2) kinase activity represents a genetically supported, chemically tractable, and potentially disease-modifying mechanism to treat Parkinson's disease. Herein, we describe the optimization of a novel series of potent, selective, central nervous system (CNS)-penetrant 1-heteroaryl-1H-indazole type I (ATP competitive) LRRK2 inhibitors. Type I ATP-competitive kinase physicochemical properties were integrated with CNS drug-like properties through a combination of structure-based drug design and parallel medicinal chemistry enabled by sp³-sp² cross-coupling technologies. This resulted in the discovery of a unique sp³-rich spirocarbonitrile motif that imparted extraordinary potency, pharmacokinetics, and favorable CNS drug-like properties. The lead compound, 25, demonstrated exceptional on-target potency in human peripheral blood mononuclear cells, excellent off-target kinase selectivity, and good brain exposure in rat, culminating in a low projected human dose and a pre-clinical safety profile that warranted advancement toward pre-clinical candidate enabling studies.

Anticancer Activity of Thiophene Carboxamide Derivatives as CA-4 Biomimetics: Synthesis, Biological Potency, 3D Spheroid Model, and Molecular Dynamics Simulation

The present study aimed to synthesize thiophene carboxamide derivatives, which are considered biomimetics of the anticancer medication Combretastatin A-4 (CA-4), and compare the similarity in the polar surface area (PSA) between the novel series and CA-4. Our results showed that the PSA of the most synthesized structures was biomimetic to CA-4, and similar chemical and biological properties were observed against Hep3B cancer cell line. Among the synthesized series 2b and 2e compounds were the most active molecules on Hep3B (IC50 = 5.46 and 12.58 µM, respectively). The 3D results revealed that both 2b and 2e structures confuse the surface of Hep3B cancer cell lines' spheroid formation and force these cells to aggregate into a globular-shaped spheroid. The 2b and 2e showed a comparable interaction pattern to that observed for CA-4 and colchicine within the tubulin-colchicine-binding pocket. The thiophene ring, due to holding a high aromaticity character, participated critically in that observed interaction profile and showed additional advanced interactions over CA-4. The 2b and 2e tubulin complexes showed optimal dynamics trajectories within a time scale of 100 ns at 300 K temperature, which asserts their high stability and compactness. Together, these findings revealed the biomimetic role of 2b and 2e compounds in CA-4 in preventing cancer progression.

The surface site interaction point approach to non-covalent interactions

The functional properties of molecular systems are generally determined by the sum of many weak non-covalent interactions, and therefore methods for predicting the relative magnitudes of these interactions is fundamental to understanding the relationship between function and structure in chemistry, biology and materials science. This review focuses on the Surface Site Interaction Point (SSIP) approach which describes molecules as a set of points that capture the properties of all possible non-covalent interactions that the molecule might make with another molecule. The first half of the review focuses on the empirical non-covalent interaction parameters, α and β, and provides simple rules of thumb to estimate free energy changes for interactions between different types of functional group. These parameters have been used to have been used to establish a quantitative understanding of the role of solvent in solution phase equilibria, and to describe non-covalent interactions at the interface between macroscopic surfaces as well as in the solid state. The second half of the review focuses on a computational approach for obtaining SSIPs and applications in multi-component systems where many different interactions compete. Ab initio calculation of the Molecular Electrostatic Potential (MEP) surface is used to derive an SSIP description of a molecule, where each SSIP is assigned a value equivalent to the corresponding empirical parameter, α or β. By considering the free energies of all possible pairing interactions between all SSIPs in a molecular ensemble, it is possible to calculate the speciation of all intermolecular interactions and hence predict thermodynamic properties using the SSIMPLE algorithm. SSIPs have been used to describe both the solution phase and the solid state and provide accurate predictions of partition coefficients, solvent effects on association constants for formation of intermolecular complexes, and the probability of cocrystal formation. SSIPs represent a simple and intuitive tool for describing the relationship between chemical structure and non-covalent interactions with sufficient accuracy to understand and predict the properties of complex molecular ensembles without the need for computationally expensive simulations.

Quantification of secondary electrostatic interactions in H-bonded complexes

The H-bonding properties of compounds that contain multiple functional groups are difficult to predict, because there are through-bond polarisation effects and long-range secondary electrostatic interactions that have significant effects on the interactions with solvents and other molecules. Here we use experimental measurements of association constants for formation of 1 : 1 H-bonded complexes that contain a single well-defined H-bond and a single well-defined secondary electrostatic interaction to quantify the magnitude of this effect. The results were used to develop a computational method for calculating functional group H-bond parameters that accurately reproduce the magnitudes of both primary H-bonding interaction and secondary electrostatic interactions. The effects of secondary electrostatic interactions are observed in calculations of ab initio Molecular Electrostatic Potential (MEP) values, but at the van der Waals surface, the magnitude of the effect is highly overestimated. MEP values calculated on electron density isosurfaces that lie closer to the nuclei provide a more accurate description of the experimental observations. H-bond parameters calculated using this approach successfully account for the properties of arrays of multiple H-bond donor and acceptor groups in different configurations. The results provide insight into the factors that govern the interaction properties of molecules that contain multiple functional groups and provide an accurate method for prediction of solution phase complexation free energies based on gas phase calculations of individual molecules.

Efficient Screening of Target-Specific Selected Compounds in Mixtures by 19F NMR Binding Assay with Predicted 19F NMR Chemical Shifts

Ligand-based 19 F NMR screening is a highly effective and well-established hit-finding approach. The high sensitivity to protein binding makes it particularly suitable for fragment screening. Different criteria can be considered for generating fluorinated fragment libraries. One common strategy is to assemble a large, diverse, well-designed and characterized fragment library which is screened in mixtures, generated based on experimental 19 F NMR chemical shifts. Here, we introduce a complementary knowledge-based 19 F NMR screening approach, named 19 Focused screening, enabling the efficient screening of putative active molecules selected by computational hit finding methodologies, in mixtures assembled and on-the-fly deconvoluted based on predicted 19 F NMR chemical shifts. In this study, we developed a novel approach, named LEFshift , for 19 F NMR chemical shift prediction using rooted topological fluorine torsion fingerprints in combination with a random forest machine learning method. A demonstration of this approach to a real test case is reported.

Evaluating anti-coronavirus activity of some phosphoramides and their influencing inhibitory factors using molecular docking, DFT, QSAR, and NCI-RDG studies

The recent prevalence of coronavirus disease in 2019 (COVID-19) has triggered widespread global health concerns.Antiviral drugs based on phosphoramides have significant inhibitory activity against the main protease (M^pro) of the virus and prevent transcription and viral replication. Hence, in order to design and introduce a group of inhibitors affecting the coronavirus, 35 phosphoramide compounds based on phospho-guanine and phospho-pyrazine derivatives were selected for molecular docking study. The results showed that most phosphoguanides containing the amino benzimidazole have a high interaction tendency with COVID-19. Among them, compound 19 was identified as the strongest inhibitor with the -9.570 kcal/mol binding energy whereas, the binding energy of Remdesivir is -6.75 kcal/mol. The quantitative structure-activity relationship (QSAR) results demonstrated that the number of aromatic rings, amide's nitrogens and their ability in π-staking, and hydrogen interactions with M^pro active sites are major factors contributing to the inhibitory activity of these compounds.Also, the NCI-RDG and DFT results were in good accordance with those of QSAR and molecular docking. The findings of this investigation can be underlying the synthesis of effective and efficient drugs against COVID-19.

Combining machine learning and quantum mechanics yields more chemically aware molecular descriptors for medicinal chemistry applications

Molecular interaction fields (MIFs), describing molecules in terms of their ability to interact with any chemical entity, are one of the most established and versatile concepts in drug discovery. Improvement of this molecular description is highly desirable for in silico drug discovery and medicinal chemistry applications. In this work, we revised a well‐established molecular mechanics' force field and applied a hybrid quantum mechanics and machine learning approach to parametrize the hydrogen‐bonding (HB) potentials of small molecules, improving this aspect of the molecular description. Approximately 66,000 molecules were chosen from available drug databases and subjected to density functional theory calculations (DFT). For each atom, the molecular electrostatic potential (EP) was extracted and used to derive new HB energy contributions; this was subsequently combined with a fingerprint‐based description of the structural environment via partial least squares modeling, enabling the new potentials to be used for molecules outside of the training set. We demonstrate that parameter prediction for molecules outside of the training set correlates with their DFT‐derived EP, and that there is correlation of the new potentials with hydrogen‐bond acidity and basicity scales. We show the newly derived MIFs vary in strength for various ring substitution in accordance with chemical intuition. Finally, we report that this derived parameter, when extended to non‐HB atoms, can also be used to estimate sites of reaction.

Hydrogen Bond Acceptor Propensity of Different Fluorine Atom Types: An Analysis of Experimentally and Computationally Derived Parameters

The propensity of organic fluorine acting as a weak hydrogen bond acceptor (HBA) in intermolecular and intramolecular interactions has been the subject of many experimental and theoretical studies often reaching different conclusions. Over the last few years, new and stronger evidences have emerged for the direct involvement of fluorine in weak hydrogen bond (HB) formation. However, not all the fluorine atom types can act as weak HBA. In this work, the differential HBA propensity of various types of fluorine atoms was analyzed with a particular emphasis for the different types of alkyl fluorides. This was carried out by evaluating ab initio computed parameters, experimental ¹⁹ F NMR chemical shifts and small molecule crystallographic structures (extracted from the CSD database). According to this analysis, shielded (with reference to the ¹⁹ F NMR chemical shift) alkyl mono-fluorinated motifs display the highest HBA propensity in agreement with solution studies. Although much weaker than other well-characterized HB complexes, the fragile HBs formed by these fluorinated motifs have important implications for the chemical-physical and structural properties of the molecules, chemical reactions, and protein-ligand recognition.

A Chemical Probe for Dark Kinase STK17B Derives Its Potency and High Selectivity through a Unique P-Loop Conformation

STK17B is a member of the death-associated protein kinase family and has been genetically linked to the development of diverse diseases. However, the role of STK17B in normal and disease pathology is poorly defined. Here, we present the discovery of thieno[3,2-d] pyrimidine SGC-STK17B-1 (11s), a high-quality chemical probe for this understudied "dark" kinase. 11s is an ATP-competitive inhibitor that showed remarkable selectivity over other kinases including the closely related STK17A. X-ray crystallography of 11s and related thieno[3,2-d]pyrimidines bound to STK17B revealed a unique P-loop conformation characterized by a salt bridge between R41 and the carboxylic acid of the inhibitor. Molecular dynamic simulations of STK17B revealed the flexibility of the P-loop and a wide range of R41 conformations available to the apo-protein. The isomeric thieno[2,3-d]pyrimidine SGC-STK17B-1N (19g) was identified as a negative control compound. The >100-fold lower activity of 19g on STK17B was attributed to the reduced basicity of its pyrimidine N1.

Discovery of Novel Imidazopyridine GSK-3β Inhibitors Supported by Computational Approaches

The interest of research groups and pharmaceutical companies to discover novel GSK-3β inhibitors has increased over the years considering the involvement of this enzyme in many pathophysiological processes and diseases. Along this line, we recently reported on 1H-indazole-3-carboxamide (INDZ) derivatives 1-6, showing good GSK-3β inhibition activity. However, they suffered from generally poor central nervous system (CNS) permeability. Here, we describe the design, synthesis, and in vitro characterization of novel imidazo[1,5-a]pyridine-1-carboxamide (IMID 1) and imidazo[1,5-a]pyridine-3-carboxamide (IMID 2) compounds (7-18) to overcome such liability. In detail, structure-based approaches and fine-tuning of physicochemical properties guided the design of derivatives 7-18 resulting in ameliorated absorption, distribution, metabolism, and excretion (ADME) properties. A crystal structure of 16 in complex with GSK-3β enzyme (PDB entry 6Y9S) confirmed the in silico models. Despite the nanomolar inhibition activity, the new core compounds showed a reduction in potency with respect to INDZ derivatives 1-6. In this context, Molecular Dynamics (MD) and Quantum Mechanics (QM) based approaches along with NMR investigation helped to rationalize the observed structure activity relationship (SAR). With these findings, the key role of the acidic hydrogen of the central core for a tight interaction within the ATP pocket of the enzyme reflecting in good GSK-3β affinity was demonstrated.

Excited state intramolecular proton transfer (ESIPT) luminescence mechanism for 4-N,N-diethylamino-3-hydroxyflavone in propylene carbonate, acetonitrile and the mixed solvents

In this work, density functional theory (DFT) and time density functional theory (TDDFT) methods were employed to investigate the nature of the double fluorescence emission of DEAHF in these three solvents. We analyzed the geometric structures, vibrational frequencies, frontier molecular orbitals (MOs), molecular electrostatic potential surface (MEPS), calculated absorption and fluorescence spectra and the potential-energy curves for DEAHF. All the results show that the intramolecular hydrogen bond of DEAHF is strengthened from S₀ to S₁ and the electron density redistribution occurs between the proton acceptor and donor, which can facilitate ESIPT. Moreover, the geometric structures, absorption and emission spectra, MEPS and potential-energy curve of DEAHF are identical. It reveals theoretically that ACN and PC can maintain the polarity of the solvent with 1:1 mixing, which is consistent with the experimental results.

Competition and cooperativity of hydrogen-bonding and tetrel-bonding interactions involving triethylene diamine (DABCO), H2O and CO2in air

Triethylene diamine (DABCO) can interact with H₂O and CO₂in air to form dimeric and trimeric complexesviahydrogen bond, tetrel bond as well as van der Waals interactions.

Fragmentation method reveals a wide spectrum of intramolecular hydrogen bond energies in antioxidant natural products

Very strong and weak IHBs in curcumin.

Synthesis, crystal structure, insecticidal activities, molecular docking and QSAR studies of some new phospho guanidines and phospho pyrazines as cholinesterase inhibitors

Novel phospho guanidine and phospho pyrazine derivatives were synthesized and characterized by ³¹P, ¹³C, ¹HNMR and IR spectroscopy to obtain novel and human-safe insecticides. Compound 35 [(C₄H₄N₂NH)₂P(O)(C₆H₆)] was investigated by X-ray crystallography. The inhibitory effects of synthesized compounds were evaluated on human and insect acetylcholinesterase (AChE) using in vitro Ellman method. A few of these compounds, which had low human toxicity, were selected for assessing the killing effects (in vivo) on the elm leaf beetle (X.luteola). The in vitro and in vivo results indicated that compounds bearing both phosphoryl groups and aromatic systems were found to possess a good selectivity for the inhibition of insect AChE over human AChE; up to 550-fold selectivity was achieved for compound 19. Docking studies were performed to explain reasons for the selective behavior of AChE inhibitors. Additionally, the quantitative structure-activity relationship (QSAR) and density functional theory (DFT) results of AChEs demonstrated that the size, shape, dipole moment, and ability to form hydrogen bond played the main role in both models. In addition, the aromatic π - π interactions and charge of the amide nitrogen had a major effect on insecticidal activity of the compounds. The present research can be helpful to gain a better understanding of the interactions between the insect AChE and its inhibitors and introduces compounds which are capable of becoming human-safe insecticides.

Exploration of N5-CAIR Mutase Novel Inhibitors from Pyrococcus horikoshii OT3: A Computational Study

In bacterial and archaeal purine biosynthetic pathways, sixth step involves utilization of enzyme PurE, catalyzing the translation of aminoimidazole ribonucleotide to 4-carboxy-5-aminoimidazole ribonucleotide (CAIR) with carbon dioxide. The formation of CAIR takes place through an unstable intermediate N5-CAIR, played by two enzymes-N5-CAIR synthetase (PurK) and N5-CAIR mutase (PurE) that further catalyzes the reaction of N5-CAIR to CAIR. In this study, N5-CAIR mutase (PH0320) from Pyrococcus horikoshii OT3 (PurE) was considered. The three-dimensional structure of Pyrococcus horikoshii OT3 was modeled based on the structure of PurE from Escherichia coli. The modeled structure was subjected to molecular dynamics simulation up to 100 ns, and least energy structure from the simulation was subjected to virtual screening and induced fit docking to identify the best potent leads. A total of five best antagonists were identified based on their affinity and mode of binding leading with conserved residues Ser18, Ser20, Asp21, Ser45, Ala46, His47, Arg48, Ala72, Gly73, Ala75, and His77 promotes the activity of Ph-N5-CAIR mutase. In addition to molecular dynamics, absorption, digestion, metabolism, and excretion properties, binding free energy and density functional theory calculations of compounds were carried out. Based on analyses, compound from National Cancer Institute (NCI) database, NCI_826 was adjudged as the best potent lead molecule and could be suggested as the suitable inhibitor of N5-CAIR mutase.

3-Fluoro-4-hydroxyprolines: Synthesis, Conformational Analysis, and Stereoselective Recognition by the VHL E3 Ubiquitin Ligase for Targeted Protein Degradation

Hydroxylation and fluorination of proline alters the pyrrolidine ring pucker and the trans:cis amide bond ratio in a stereochemistry-dependent fashion, affecting molecular recognition of proline-containing molecules by biological systems. While hydroxyprolines and fluoroprolines are common motifs in medicinal and biological chemistry, the synthesis and molecular properties of prolines containing both modifications, i.e., fluoro-hydroxyprolines, have not been described. Here we present a practical and facile synthesis of all four diastereoisomers of 3-fluoro-4-hydroxyprolines (F-Hyps), starting from readily available 4-oxo-l-proline derivatives. Small-molecule X-ray crystallography, NMR spectroscopy, and quantum mechanical calculations are consistent with fluorination at C3 having negligible effects on the hydrogen bond donor capacity of the C4 hydroxyl, but inverting the natural preference of Hyp from C4-exo to C4-endo pucker. In spite of this, F-Hyps still bind to the von Hippel-Lindau (VHL) E3 ligase, which naturally recognizes C4-exo Hyp in a stereoselective fashion. Co-crystal structures and electrostatic potential calculations support and rationalize the observed preferential recognition for (3 R,4 S)-F-Hyp over the corresponding (3 S,4 S) epimer by VHL. We show that (3 R,4 S)-F-Hyp provides bioisosteric Hyp substitution in both hypoxia-inducible factor 1 alpha (HIF-1α) substrate peptides and peptidomimetic ligands that form part of PROTAC (proteolysis targeting chimera) conjugates for targeted protein degradation. Despite a weakened affinity, Hyp substitution with (3 S,4 S)-F-Hyp within the PROTAC MZ1 led to Brd4-selective cellular degradation at concentrations >100-fold lower than the binary Kd for VHL. We anticipate that the disclosed chemistry of 3-fluoro-4-hydroxyprolines and their application as VHL ligands for targeted protein degradation will be of wide interest to medicinal organic chemists, chemical biologists, and drug discoverers alike.

Intermolecular interactions between σ- and π-holes of bromopentafluorobenzene and pyridine: computational and experimental investigations

The characters of σ- and π-holes of bromopentafluorobenzene (C6F5Br) enable it to interact with an electron-rich atom or group like pyridine which possesses an electron lone-pair N atom and a π ring. Theoretical studies of intermolecular interactions between C6F5Br and C5H5N have been carried out at the M06-2X/aug-cc-pVDZ level without and with the counterpoise method, together with single point calculations at M06-2X/TZVP, wB97-XD/aug-cc-pVDZ and CCSD(T)/aug-cc-pVDZ levels. The σ- and π-holes of C6F5Br exhibiting positive electrostatic potentials make these sites favorably interact with the N atom and the π ring of C5H5N with negative electrostatic potentials, leading to five different dimers connected by a σ-holen bond, a σ-holeπ bond or a π-holeπ bond. Their geometrical structures, characteristics, nature and spectroscopy behaviors were systematically investigated. EDA analyses reveal that the driving forces in these dimers are different. NCI, QTAIM and NBO analyses confirm the existence of intermolecular interactions formed via σ- and π-holes of C6F5Br and the N atom and the π ring of C5H5N. The experimental IR and Raman spectra gave us important information about the formation of molecular complexes between C6F5Br and C5H5N. We expect that the results could provide valuable insights into the investigation of intermolecular interactions involving σ- and π-holes.

Effects of the Hofmeister series of sodium salts on the solvent properties of water

The solvent features of water (solvent dipolarity/polarizability, π*, hydrogen bond donor acidity, α, and hydrogen bond acceptor basicity, β) were examined in aqueous solutions of Na₂SO₄, NaF, CH₃COONa, NaCl, NaBr, NaI, and NaClO₄ at concentrations of each salt from 0 to 1.0 M (up to 2.0 M for NaClO₄). The solvent features of water in solutions of different concentrations for each salt were found to be linearly related as π* = z + aα + bβ. The coefficients of this relationship were suggested to represent the signature of the salt effect on the solvent features of water. The normalized distances for each salt were calculated using glucose as a reference compound. These distances may be used as the relative measures of the salt-water interactions. It is demonstrated that the distances for all salts examined are interrelated with structural water entropies and static polarizabilities of anions. It is shown that the distance may be used as a measure of the relative effects of salts on precipitation of ferric oxide, excessive chemical potential of propanol in salt solutions, surface tension, and viscosity. The distance represents the relative measure of the salt effect on the solvent features of water in a salt solution. The examples presented confirm that the approach used does enable us to characterize the differences between the effects of salts in the Hofmeister series on the properties of water.

Photoexcitation effect on the adsorption of hazardous gases on silica surface

There is very little scientific understanding of photoexcitation effect on the adsorption properties of adsorbent. The adsorption of four hazardous gases (SARIN (propan-2-ylmethylphospho-nofluoridate), methyl dichlorophosphate (MDCP), trimethyl phosphate (TMP) and hydrogen sulfide (H₂S)) on silica surface is taken as target sample in this work. The adsorption energy order (MDCP<SARIN<TMP) in the ground state is consistent with the strength order of intermolecular hydrogen bond (inter-HB) between hydroxyl group of silica surface and hazardous gas, and the desorption order of the three gases in previous reports. However, with the adsorption energy increase of MDCP and the decrease of SARIN and TMP, this order changes remarkably to SARIN<TMP<MDCP after photoexcitation to excited state by absorbing shortwave ultraviolet irradiation. This change is opposite to the inter-HB weakening of MDCP in the first excited (S₁) state and the strengthening of TMP and SARIN in the second excited (S₂) state. This opposite change is caused by formation of intermolecular charge transfer state of MDCP and local excitation of SARIN and TMP. The H₂S is dissociated after photoexcitation to the S₁ state. This work presents photoexcitation as a new standard for the design and detection of adsorption properties of adsorbent for its striking effect on adsorption behaviors.

Electronic Structures of LNA Phosphorothioate Oligonucleotides

Important oligonucleotides in anti-sense research have been investigated in silico and experimentally. This involves quantum mechanical (QM) calculations and chromatography experiments on locked nucleic acid (LNA) phosphorothioate (PS) oligonucleotides. iso-potential electrostatic surfaces are essential in this study and have been calculated from the wave functions derived from the QM calculations that provide binding information and other properties of these molecules. The QM calculations give details of the electronic structures in terms of e.g., energy and bonding, which make them distinguish or differentiate between the individual PS diastereoisomers determined by the position of sulfur atoms. Rules are derived from the electronic calculations of these molecules and include the effects of the phosphorothioate chirality and formation of electrostatic potential surfaces. Physical and electrochemical descriptors of the PS oligonucleotides are compared to the experiments in which chiral states on these molecules can be distinguished. The calculations demonstrate that electronic structure, electrostatic potential, and topology are highly sensitive to single PS configuration changes and can give a lead to understanding the activity of the molecules.

The theoretical study of excited-state intramolecular proton transfer of 2,5-bis(benzoxazol-2-yl)thiophene-3,4-diol

The symmetrical structures 2,5-bis(benzoxazol-2-yl)thiophene-3,4-diol (BBTD) can take shape two intramolecular hydrogen bonds in chloroform. In order to research the molecular dynamic behavior of BBTD upon photo-induced process, we utilize density functional theory (DFT) and time-dependent density functional theory (TDDFT) to complete theoretical calculation. Through the comparison of bond length, bond angle, IR spectra, and frontier molecular orbitals between ground state (S₀) and first excited state (S₁), it clearly indicates that photoexcitation have slightly influence for intensity of hydrogen bond. For the sake of understanding the mechanism of excited state intramolecular proton transfer (ESIPT) of BBTD in chloroform, potential energy surfaces have been scanned along with the orientation of O₁-H₂ and O₄-H₅ in S₀ and S₁ state, respectively. A intrigued hydrogen bond dynamic phenomenon has been found that ESIPT of BBTD is not a synergetic double proton transfer process, but a stepwise single proton transfer process BBTD→BBTD-S→BBTD-D. Moreover, the proton transfer process of BBTD-S→BBTD-D is easier to occur than that of BBTD→BBTD-S in S₁ state.