A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

The Optical Properties, UV-Vis. Absorption and Fluorescence Spectra of 4-Pentylphenyl 4-n-benzoate Derivatives in Different Solvents

In this paper, electronic absorbance and fluorescence spectra of 4-Pentylphenyl 4-n-benzoate derivatives have been measured in 29 solvents, which are non-polar, polar protic and polar aprotic solvents, and electronic transitions that vary depending on the solvent are identified. As the solvent polarity increases, the forbidden energy difference between the frontier orbitals decreases. The statistical models in order to describe the solvent effect were derived using different solvent parameters. Quite complex and multiple absorbance transitions were observed in different solvent environments. Local fluorescence transition and intramolecular charge transfer occurred in the fluorescence spectra. Absorbance transitions are global transitions, and π*←π is the absorbance electronic transition. The frontier molecular orbitals and electrostatic potential surface were founded using quantum chemical calculations. Refractive indices were found with five different methods and forbidden energy gaps were found with the Tauc method. The forbidden energy ranges were found around 4.1 and 4.5, and the forbidden energy gap decreased as the alkyl chain became longer. All compounds can be defined as insulation materials according to the forbidden energy range. Refractive index values close to the E7 liquid crystal mixture used in liquid crystal display panels were found in the investigated liquid crystals.

Piperine solubility enhancement via DES formation: Elucidation of intermolecular interactions and impact of counterpart structure via computational and spectroscopic approaches

The development of new forms of existing APIs with enhanced physicochemical properties is critical for improving their therapeutic potential. In this context, ionic liquids (ILs) and deep eutectic solvents (DESs) have gained significant attention in recent years due to their unique properties and potential for solubility enhancement. In this study, we explore the role of different counterparts in the formation of IL/DESs with piperine (PI), a poorly water-soluble drug. After screening a library of fourteen counterpart molecules, ten liquid PI-counterpart systems were developed and investigated. Thermal analysis confirmed the formation of IL/DES, while computational and spectroscopic studies revealed that hydrogen bonding played a crucial role in the interaction between PI and the counterparts, confirming DES formation. The solubility enhancement of PI in these systems ranged from ∼ 36 % to 294 %, with PI-Oxalic acid (OA) exhibiting the highest saturation solubility (49.71 μg/mL) and PI-Ibuprofen (IB) the lowest (17.23 μg/mL). The presence of hydrogen bonding groups in counterparts was key to successful DES formation. A negative correlation was observed between solubility and logP (r =  - 0.75, p* = 0.0129), while a positive correlation was found between solubility and normalized polar surface area (PSA) (r = 0.68, p* = 0.029). PI-OA and PI-IB were located at the extreme ends of these regression lines, further validating the relationship between these properties and solubility enhancement. These findings highlight essential aspects of rational IL/DES design, optimizing their properties for broader applications.

Investigation of cross-association behavior in water-ethanol solutions: A combined computational-ATR spectroscopy study

The water/ethanol system possesses complexities at the molecular level, which render its description a difficult task. For the elucidation of the system's hydrogen bonding features that are the key factors in its complex behavior, we conduct a Density Functional Theory analysis on relevant water/ethanol clusters inside implicit solvent cavities for the determination of the ethanol donor hydrogen bond strength. We record Attenuated Total Reflectance spectra of water/ethanol-OD solutions and utilize our density and refractive index measurements for post-processing. The application of the Badger-Bauer rule reveals a minimum in the strength of the ethanol donor hydrogen bond for a composition of xwater = 0.74. We attempt to analyze further this result by estimating the effect of the implicit solvent on the ethanol donor hydrogen bond strength, finding it to be incremental. A brief analysis of different cluster conformations is carried out to determine the cooperativity conditions that can potentially explain the observed minimum in the ethanol donor hydrogen bond strength. These observations are related to notions of microheterogeneity in water/alcohol mixtures and provide context toward a more elaborate picture of association in heteroclusters.

Synthesis, docking and biological evaluation of purine-5-N-isosteresas anti-inflammatory agents

An operationally simple one-pot three-component and convenient synthesis method for a series of diverse purine analogues of 5-amino-7-(substituted)-N-(4-sulfamoylphenyl)-4,7-dihydro-[1,2,4]-triazolo[1,5-a][1,3,5]triazine-2-carboxamide derivatives generated in situ via the reaction of 2-hydrazinyl-N-(4-sulfamoylphenyl)-2-thioxoacetamide, cyanoguanidine and a variety of aldehydes was achieved under green conditions. This experiment was conducted to evaluate the anti-inflammatory effect of the newly synthesized compounds using indomethacin as a reference medication; all compounds were tested for in vitro anti-inflammatory activity using the inhibition of albumin denaturation, RBC hemolysis technique and COX inhibition assay. The results showed that all evaluated compounds exhibited significant in vitro anti-inflammatory efficacy leading to excellently effective RBC membrane stabilization, inhibition of protein denaturation, and inhibition of COX enzymes when compared to those of indomethacin. At concentrations of 50, 100, 200, and 300 μg ml-1, these compounds decreased COX-1 and COX-2 activities more than indomethacin and have IC50 values in the range of 40.04-87.29 μg ml-1 for COX-1 and 27.76-42.3 μg ml-1 for COX-2 while indomethacin showed IC50 = 91.57 for COX-1 and 42.66 μg ml-1 for COX-2. The anti-inflammatory findings show the need for more investigation to define the properties underlying the evaluated compounds' anti-inflammatory abilities. The enzyme cyclooxygenase-2 (COX 2) (PDB ID: 5IKT) was docked with ten synthetic substances. With docking scores (S) of -8.82, -7.82, and -7.76 kcal mol-1, 7-furan triazolo-triazine (4), 7-(2-hydroxy phenyl) triazolo-triazine (11), and 7-(4-dimethylamino phenyl) triazolo-triazine (12) had the greatest binding affinities, respectively. Therefore, these substances have COX-2 (PDB ID: 5IKT) inhibitory capabilities and hence may be investigated for COX 2 targeting development. Furthermore, both the top-ranked compounds (4 and 11) and the standard indomethacin were subjected to DFT analysis. The HOMO - LUMO energy difference (ΔE) of the mentioned compounds was found to be less than that of indomethacin.

Effect of functionalization on the adsorption performance of carbon nanotube as a drug delivery system for imatinib: molecular simulation study

In this study, efficiency of functionalized carbon nanotube as a potential delivery system for imatinib anti-cancer drug was investigated. Accordingly, carboxyl and hydroxyl functionalized carbon nanotube were inspected as a notable candidate for the carriage of this drug in aqueous media. For this purpose, possible interactions of imatinib with pure and functionalized carbon nanotube were considered in aqueous media. The compounds were optimized in gas phase using density functional calculations. Solvation free energies and association free energies of the optimized structures were then studied by Monte Carlo simulation and perturbation method in water environment. Outcomes of quantum mechanical calculations presented that pure and functionalized carbon nanotubes can act as imatinib drug adsorbents in gas phase. However, results of association free energy calculations in aqueous solution indicated that only carboxyl and hydroxyl functionalized carbon nanotubes could interact with imatinib. Monte Carlo simulation results revealed that electrostatic interactions play a vital role in the intermolecular interaction energies after binding of drug and nanotube in aqueous solution. Computed solvation free energies in water showed that the interactions with functionalized carbon nanotubes significantly enhance the solubility of imatinib, which could improve its in vivo bioavailability.

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. By performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

Exploring borderline SN1-SN2 mechanisms: the role of explicit solvation protocols in the DFT investigation of isopropyl chloride

Nucleophilic substitution at saturated carbon is a crucial class of organic reactions, playing a pivotal role in various chemical transformations that yield valuable compounds for society. Despite the well-established SN1 and SN2 mechanisms, secondary substrates, particularly in solvolysis reactions, often exhibit a borderline pathway. A molecular-level understanding of these processes is fundamental for developing more efficient chemical transformations. Typically, quantum-chemical simulations of the solvent medium combine explicit and implicit solvation methods. The configuration of explicit molecules can be defined through top-down approaches, such as Monte Carlo (MC) calculations for generating initial configurations, and bottom-up methods that involve user-dependent protocols to add solvent molecules around the substrate. Herein, we investigated the borderline mechanism of the hydrolysis of a secondary substrate, isopropyl chloride (iPrCl), at DFT-M06-2X/aug-cc-pVDZ level, employing explicit and explicit + implicit protocols. Top-down and bottom-up approaches were employed to generate substrate-solvent complexes of varying number (n = 1, 3, 5, 7, 9, and 12) and configurations of H2O molecules. Our findings consistently reveal that regardless of the solvation approach, the hydrolysis of iPrCl follows a loose-SN2-like mechanism with nucleophilic solvent assistance. Increasing the water cluster around the substrate in most cases led to reaction barriers of ΔH‡ ≈ 21 kcal mol-1, with nine water molecules from MC configurations sufficient to describe the reaction. The More O'Ferrall-Jencks plot demonstrates an SN1-like character for all transition state structures, showing a clear merged profile. The fragmentation activation strain analyses indicate that energy barriers are predominantly controlled by solvent-substrate interactions, supported by the leaving group stabilization assessed through CHELPG atomic charges.

The stability of oxygen-centered radicals and its response to hydrogen bonding interactions

The stability of various alkoxy/aryloxy/peroxy radicals, as well as TEMPO and triplet dioxygen (3 O2 ) has been explored at a variety of theoretical levels. Good correlations between RSEtheor and RSEexp are found for hybrid DFT methods, for compound schemes such as G3B3-D3, and also for DLPNO-CCSD(T) calculations. The effects of hydrogen bonding interactions on the stability of oxygen-centered radicals have been probed by addition of a single solvating water molecule. While this water molecule always acts as a H-bond donor to the oxygen-centered radical itself, it can act as a H-bond donor or acceptor to the respective closed-shell parent.

Structure and stability of the sH binary hydrate cavity and host-guest versus guest-guest interactions therein: A DFT approach

The intrinsic ability of clathrate hydrates to encage gaseous molecules is explored. Encapsulation ability depends on the cavity size and the type of guest gaseous species in addition to the physical parameters, temperature and pressure. Here we have reported the structure, stability and nature of interaction in dissimilar guest occupied sH hydrate cavity. Diatomic gas molecules and small polyatomic hydrocarbons are considered as guests. The irregular icosahedron (5¹² 6⁸ ) cavity of sH hydrate is the host. Different thermodynamic parameters for the guest molecules encapsulation were calculated using three different hybrid DFT functionals, B3LYP, M05-2X, M06, and moreover using dispersion correction (PBE0-D3). With the consideration of large H-bonded systems the 6-31G* and cc-pVTZ basis sets were used for two set of computations. To disclose the nature of interaction between the host-guest systems as well as the interaction between the guest molecules inside the host the non-covalent interaction (NCI) indices and energy decomposition analysis (EDA) were done. Impact of host-guest and guest-guest interactions are discussed.

How Clustered DNA Damage Can Change the Electronic Properties of ds-DNA—Differences between GAG, GAOXOG, and OXOGAOXOG

Every 24 h, roughly 3 × 1017 incidences of DNA damage are generated in the human body as a result of intra- or extra-cellular factors. The structure of the formed lesions is identical to that formed during radio- or chemotherapy. Increases in the clustered DNA damage (CDL) level during anticancer treatment have been observed compared to those found in untreated normal tissues. 7,8-dihydro-8-oxo-2′-deoxyguanosine (OXOG) has been recognized as the most common lesion. In these studies, the influence of OXOG, as an isolated (oligo-OG) or clustered DNA lesion (oligo-OGOG), on charge transfer has been analyzed in comparison to native oligo-G. DNA lesion repair depends on the damage recognition step, probably via charge transfer. Here the electronic properties of short ds-oligonucleotides were calculated and analyzed at the M062x/6-31++G** level of theory in a non-equilibrated and equilibrated solvent state. The rate constant of hole and electron transfer according to Marcus’ theory was also discussed. These studies elucidated that OXOG constitutes the sink for migrated radical cations. However, in the case of oligo-OGOG containing a 5′-OXOGAXOXG-3′ sequence, the 3′-End OXOG becomes predisposed to electron-hole accumulation contrary to the undamaged GAG fragment. Moreover, it was found that the 5′-End OXOG present in an OXOGAOXOG fragment adopts a higher adiabatic ionization potential than the 2′-deoxyguanosine of an undamaged analog if both ds-oligos are present in a cationic form. Because increases in CDL formation have been observed during radio- or chemotherapy, understanding their role in the above processes can be crucial for the efficiency and safety of medical cancer treatment.

The 2Ih and OXOG Proximity Consequences on Charge Transfer through ds-DNA: Theoretical Studies of Clustered DNA Damage

Genetic information is continuously exposed to harmful factors, both intra- and extracellular. Their activity can lead to the formation of different types of DNA damage. Clustered lesions (CDL) are problematic for DNA repair systems. In this study, the short ds-oligos with a CDL containing (R) or (S) 2Ih and OXOG in their structure were chosen as the most frequent in vitro lesions. In the condensed phase, the spatial structure was optimized at the M062x/D95**:M026x/sto-3G level of theory, while the electronic properties were optimized at the M062x/6-31++G** level. The influence of equilibrated and non-equilibrated solvent-solute interactions was then discussed. It was found that the presence of (R)2Ih in the ds-oligo structure causes a greater increase in structure sensitivity towards charge adoption than (S)2Ih, while OXOG shows high stability. Moreover, the analysis of charge and spin distribution reveals the different effects of 2Ih diastereomers. Additionally, the adiabatic ionization potential was found as follows for (R)-2Ih and (S)-2Ih in eV: 7.02 and 6.94. This was in good agreement with the AIP of the investigated ds-oligos. It was found that the presence of (R)-2Ih has a negative influence on excess electron migration through ds-DNA. Finally, according to the Marcus theory, the charge transfer constant was calculated. The results presented in the article show that both diastereomers of 5-carboxamido-5-formamido-2-iminohydantoin should play a significant role in the CDL recognition process via electron transfer. Moreover, it should be pointed out that even though the cellular level of (R and S)-2Ih has been obscured, their mutagenic potential should be at the same level as other similar guanine lesions found in different cancer cells.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

Supramolecular Diversity, Theoretical Investigation and Antibacterial Activity of Cu, Co and Cd Complexes Based on the Tridentate N,N,O-Schiff Base Ligand Formed In Situ

The four new complexes, [Cu(HL¹)(L²)Cl] (1), [Cu(HL¹)(L¹)]∙Cl∙2H₂O (2), [Co(L¹)₂]∙Cl (3) and [Cd(HL¹)I₂]∙dmso (4), have been prepared by one-pot reactions of the respective chloride or iodide metal salt with a non-aqueous solution of the polydentate Schiff base, HL¹, resulted from in situ condensation of benzhydrazide and 2-pyridinecarboxaldehyde, while a ligand HL², in case of 1, has been formed due to the oxidation of 2-pyridinecarboxaldehyde under reaction conditions. The crystallographic analysis revealed that the molecular building units in 1-4 are linked together into complex structures by hydrogen bonding, resulting in 1D, 2D and 3D supramolecular architectures for 1, 2 and 4, respectively, and the supramolecular trimer for 3. The electronic structures of 1-4 were investigated by the DFT theoretical calculations. The non-covalent interactions in the crystal structures of 1-4 were studied by means of the Hirshfeld surface analysis and the QTAIM theory with a special focus on the C-H⋯Cl bonding. From the DFT/DLPNO-CCSD(T) calculations, using a series of charged model {R₃C-H}⁰⋯Cl^- assemblies, we propose linear regressions for assessment of the interaction enthalpy (ΔH, kcal mol^-1) and the binding energy (BE, kcal mol^-1) between {R₃C-H}⁰ and Cl^- sites starting from the electron density at the bond critical point (ρ(r_BCP), a.u.): ΔH = -678 × ρ(r) + 3 and BE = -726 × ρ(r) + 4. It was also has been found that compounds 1, 3 and 4 during in vitro screening showed an antibacterial activity toward the nine bacteria species, comprising both Gram-positive and Gram-negative, with MIC values ranging from 156.2 to 625 mg/L. The best results have been obtained against Acinetobacter baumannii MβL.

Molecular mechanisms from reaction coordinate graph enabled multidimensional free energies illustrated on water dimer hydrogen bonding

Computing the free energies of molecular mechanisms in multidimensional space relies on combinations of geometrically complex reaction coordinates. We show how a graph theory implementation reduces complexity, and illustrate this on the arrangements of hydrogen bonding of a water dimer. The reaction coordinates and forces are computed using graphs that define the dependencies on the atoms in the Free Energy from Adaptive Reaction Coordinate Forces (FEARCF) library. The library can be interfaced with classical molecular dynamics as well as quantum molecular dynamics packages. Multidimensional interdependent reaction coordinates are constructed to produce complex free energy hypersurfaces. The reaction coordinates are graphed from atomic and molecular components to define points, distances, vectors, angles, planes and combinations thereof. The resultant free energy surfaces that are a function of the distance, angles, planes, and so on, can represent molecular mechanisms in reduced dimensions from the component atomic Cartesian coordinate degrees of freedom. The FEARCF library can be interfaced with any molecular package. Here, we demonstrate the link to NWChem to compute a hyperdimensional DFT (aug‐cc‐pVDZ basis set and X3LYP exchange correlation functionals) free energy space of a water dimer. Analysis of the water dimer free energy hypervolume reveals that while the chain and cyclic hydrogen bonding configurations are located in stable minimum energy wells, the bifurcated hydrogen bond configuration is a gateway to instability and dimer dissociation.

Estimations of OH·N hydrogen bond length from positions and intensities of IR bands

In this computational work applicability of IR spectral parameters for evaluations of OH···N hydrogen bond length is discussed. For a set of 124 complexes with OH···N hydrogen bond formed by combinations of methanol/acetic acid and pyridine (and their fluorine substituted versions) geometries, energies and IR parameters were calculated at MP2/def2-TZVP level of theory. For a number of IR parameters (the shift of proton donor group stretching vibration Δν_s, increase of its intensity I, the low-frequency hydrogen bond stretching vibration ν_σ, bending in-plane δ and out-of-plane γ vibrations) equations linking them with interatomic distances are proposed, the robustness and accuracy of such equations are discussed. The enthalpy of OH···N hydrogen bond formation ΔH was also linked with electron density parameters in (3; -1) critical point.

Weakly Bound Complex Formation between HCN and CH3Cl: A Matrix-Isolation and Computational Study

The matrix-isolated infrared spectrum of a hydrogen cyanide-methyl chloride complex was investigated in a solid argon matrix. HCN and CH₃Cl were co-condensed onto a substrate held at 10 K with an excess of argon gas, and the infrared spectrum was measured using Fourier-transform infrared spectroscopy. Quantum chemical geometry optimization, harmonic frequency, and natural bonding orbital calculations indicate stabilized hydrogen- and halogen-bonded structures. The two resulting weakly bound complexes are both composed of one CH₃Cl molecule bound to a (HCN)₃ subunit, where the three HCN molecules are bound head-to-tail in a ring formation. Our study suggests that─in the presence of CH₃Cl─the formation of (HCN)₃ is promoted through complexation. Since HCN aggregates are an important precursor to prebiotic monomers (amino acids and nucleobases) and other life-bearing polymers, this study has astrophysical implications toward the search for life in space.

Inclusive DFT insight into sensing mechanism of cyclotetrapyrole towards lung irritants

The development of smart sensing devices for toxic analytes detection especially lung irritants is much essential. The cyclic conducting polymers having infinite π-conjugation are proved to be highly sensitive for toxic analytes. Herein, by using the DFT approach, we investigated the sensing mechanism of cyclotetrapyrole (CTPy) for accurate detection of phosgene, diphosgene, chloropicrin and chlorine at the B3LYP-D3/6-31 + G (d, p) level. The calculated interaction energies show the physisorption of analytes over the CTPy surface. Natural bond orbital (NBO) and charge decomposition (CDA) analyses predict charge transfer interactions in the complexes. The reduced density gradient (RDG) approach reveals that hydrogen bonding interactions dominate in the complexes. The sensitivity of CTPy towards lung irritants is further illustrated by the reduction in HOMO-LUMO energy gaps, red shifting of [Formula: see text] in UV-Visible spectra. Density of state (DOS) analysis affirm that enhanced conductivity upon complexation is due to the origination of new energy states in occupied and virtual orbitals nearer to the Fermi level. Moreover, PDOS spectra show that CTPy primarily contributes to the energy of HOMO. The outcome of the current study depicts appreciable sensitivity of CTPy towards lung irritants. Moreover, the competing role of naturally occurring atmospheric water is also investigated. We believe that the upshot of the current findings and their forecasts will provide useful guidelines for an experimentalist to design highly sensitive sensors for toxic analytes using CTPy. HIGHLIGHTS: • The highest QNBO transfer towards the analyte (- 0.121) is seen in the chlorine@CTPy complex. • The highest reduction in Eg (61%) between occupied and virtual orbitals is noticed in chlorine@CTPy. • The orbital overlap results in a 41% red shifting of [Formula: see text] in chlorine@CTPy. • Cyclotetrapyrole is highly sensitive for chlorine.

A Study on the Effect of the Substituent against PAK4 Inhibition Using In Silico Methods

The intrinsic inductive properties of atoms or functional groups depend on the chemical properties of either electron-withdrawing groups (EWGs) or electron-donating groups (EDGs). This study aimed to evaluate in silico methods to determine whether changes in chemical properties of the compound by single atomic substitution affect the biological activity of target proteins and whether the results depend on the properties of the functional groups. We found an imidazo[4,5-b]pyridine-based PAK4 inhibitor, compound 1, as an initial hit compound with the well-defined binding mode for PAK4. In this study, we used both experimental and in silico methods to investigate the effect of atomic substitution on biological activity to optimize the initial hit compound. In biological assays, in the case of EWG, as the size of the halogen atom became smaller and the electronegativity increased, the biological activity IC₅₀ value ranged from 5150 nM to inactive; in the case of EDG, biological activity was inactive. Furthermore, we analyzed the interactions of PAK4 with compounds, focusing on the hinge region residues, L398 and E399, and gatekeeper residues, M395 and K350, of the PAK4 protein using molecular docking studies and fragment molecular orbital (FMO) methods to determine the differences between the effect of EWG and EDG on the activity of target proteins. These results of the docking score and binding energy did not explain the differences in biological activity. However, the pair-interaction energy obtained from the results of the FMO method indicated that there was a difference in the interaction energy between the EWG and EDG in the hinge region residues, L398 and E399, as well as in M395 and K350. The two groups with different properties exhibited opposite electrostatic energy and charge transfer energy between L398 and E399. Additionally, we investigated the electron distribution of the parts interacting with the hinge region by visualizing the molecular electrostatic potential (MEP) surface of the compounds. In conclusion, we described the properties of functional groups that affect biological activity using an in silico method, FMO.

A general justification for hybrid functionals in DFT by means of linear response theory

In the present work, resorting to linear response theory, we examine the plausibility of postulating Kohn-Sham (KS)-type equations which contain, by definition, an effective hybrid potential made up by some arbitrary mixture of local and non-local terms. In this way a general justification for the construction of hybrid functionals is provided without resorting to arguments based on the adiabatic connection, the generalized KS theory or the Levy's constrained search (or its variations). In particular, we examine the cases of single-hybrid functionals, derived from non-local exchange and of double-hybrid functionals, emerging from non-local second-order expressions obtained from the KS perturbation theory. A further generalization for higher-order hybrid functionals is also included.

Synthesis and evaluation of potent yaku'amide A analogs

Two full-length analogs of the anticancer peptide yaku'amide A (1a) and four partial structures have been synthesized. These analogs were identified by computational studies in which the three E- and Z-ΔIle residues of the natural product were replaced by the more accessible dehydroamino acids ΔVal and ΔEnv. Of the eight possible analogs, modeling showed that the targeted structures 2a and 2b most closely resembled the three-dimensional structure of 1a. Synthesis of 2a and 2b followed a convergent route that was streamlined by the absence of ΔIle in the targets. Screening of the compounds against various cancer cell lines revealed that 2a and 2b mimic the potent anticancer activity of 1a, thereby validating the computational studies.

Simplified full-length analogs of yaku'amide A were designed with the aid of computations and then synthesized. The analog that was predicted to most closely resemble the structure of the natural product mimicked its anticancer activity.

All-Atom Molecular Dynamics Investigations on the Interactions between D2 Subunit Dopamine Receptors and Three 11C-Labeled Radiopharmaceutical Ligands

The D2 subunit dopamine receptor represents a key factor in modulating dopamine release. Moreover, the investigated radiopharmaceutical ligands used in positron emission tomography imaging techniques are known to bind D2 receptors, allowing for dopaminergic pathways quantification in the living human brain. Thus, the biophysical characterization of these radioligands is expected to provide additional insights into the interaction mechanisms between the vehicle molecules and their targets. Using molecular dynamics simulations and QM calculations, the present study aimed to investigate the potential positions in which the D2 dopamine receptor would most likely interact with the three distinctive synthetic 11C-labeled compounds (raclopride (3,5-dichloro-N-[[(2S)-1-ethylpyrrolidin-2-yl]methyl]-2-hydroxy-6-methoxybenzamide)—RACL, FLB457 (5-bromo-N-[[(2S)-1-ethylpyrrolidin-2-yl]methyl]-2,3-dimethoxybenzamide)—FLB457 and SCH23390 (R(+)-7-Chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine)—SCH)), as well as to estimate the binding affinities of the ligand-receptor complexes. A docking study was performed prior to multiple 50 ns molecular dynamics productions for the ligands situated at the top and bottom interacting pockets of the receptor. The most prominent motions for the RACL ligand were described by the high fluctuations of the peripheral aliphatic -CH3 groups and by its C-Cl aromatic ring groups. In good agreement with the experimental data, the D2 dopamine receptor-RACL complex showed the highest interacting patterns for ligands docked at the receptor’s top position.

K T, Ismail TM, Sajith PK. Theoretical investigation into the effect of water on the N2O decomposition reaction over Cu-ZSM-5 catalyst

Copper exchanged zeolites are an admirable catalyst for the direct decomposition reaction of harmful N2O gas. However, the inhibition of the decomposition reaction in the presence of water vapor greatly...

The influence of oxoG on the electronic properties of ds-DNA. Damage versus mismatch: A theoretical approach

The seed of life is concealed in the base sequence in DNA. This macromolecule is continuously exposed to harmful factors which can cause it damage. The stability of genetic information depends on the protein efficiency of repair systems. Glycosylases are the scouts which recognize and remove damaged bases. Their efficiency depends on how rapidly they recognize DNA lesions. One theory states that charge transfer is involved in protein cross talking through ds-DNA. For these reasons a comparative analysis of ds-oligo containing a mismatched base pair dA:::dG and a damaged dA::dG^OXO is proposed. Additionally, the electronic properties of the short ds-oligo in the context of non-equilibrated and equilibrated solvent modes were taken into theoretical consideration. All energetic calculations were performed at the M062x/6-31++G** level of theory, while for geometry optimized ONIOM methodology was used. The lowest adiabatic ionization potential was assigned for DNA containing a dA:dG^OXO pair. Moreover, the adiabatic electron affinity was assigned at the same level for the mismatched and lesioned ds-oligo. Surprisingly, in the non-equilibrated mode, a significantly higher vertical electro affinity was found for lesioned DNA. The higher VEA in a non-equilibrated solvent state supported faster recognition in the A:G^OXO base pair than A:G by MutY glycosylases under electron transfer mechanism.

A molecular dynamics simulation study on the propensity of Asn-Gly-containing heptapeptides towards β-turn structures: Comparison with ab initio quantum mechanical calculations

Both molecular mechanical and quantum mechanical calculations play an important role in describing the behavior and structure of molecules. In this work, we compare for the same peptide systems the results obtained from folding molecular dynamics simulations with previously reported results from quantum mechanical calculations. More specifically, three molecular dynamics simulations of 5 μs each in explicit water solvent were carried out for three Asn-Gly-containing heptapeptides, in order to study their folding and dynamics. Previous data, based on quantum mechanical calculations within the DFT framework have shown that these peptides adopt β-turn structures in aqueous solution, with type I’ β-turn being the most preferred motif. The results from our analyses indicate that at least for the given systems, force field and simulation protocol, the two methods diverge in their predictions. The possibility of a force field-dependent deficiency is examined as a possible source of the observed discrepancy.

Towards more effective beryllium chelation: an investigation of second-sphere hydrogen bonding

A comparative study between three experimentally known beryllium chelators (EDTA, NTP, and 10-HBQS) and two tetradentate tripodal di-pyridine-based receptors (HL and HL-NH2), specifically designed to bind Be2+ cations, has been undertaken in the aqueous phase at the B3LYP/6-311++G(d,p) computational level. The relative binding energies of these five ligand systems to a variety of first row and pre-transition metal cations have been calculated, specifically to investigate their binding strength to Be2+ and the binding enhancement that a second sphere hydrogen bonding interaction could afford to the pyridyl based systems. The complexes of EDTA were calculated to have the highest average binding energy; followed by those of NTP, HL-NH2, HL, and finally 10-HBQS. The calculated binding energy of the HL-NH2Be complex, which includes second sphere interactions, was found to be almost 9% greater than the HL Be complex, with an average binding energy increase of 13.5% observed across all metals upon inclusion of second sphere hydrogen bonding.

The Electronic Property Differences between dA::dG and dA::dGoxo. A Theoretical Approach

The dA::dG^oxo pair appearing in nucleic ds-DNA can lead to a mutation in the genetic information. Depending on the dG^oxo source, an AT→GC and GC→AC transversion might be observed. As a result, glycosylases are developed during the evolution, i.e., OGG1 and MutY. While the former effectively removes G^oxo from the genome, the second one removes adenine from the dA::dG^oxo and dA:dG pair. However, dA::dG^oxo is recognized by MutY as ~6–10 times faster than dA:dG. In this article, the structural and electronic properties of simple nucleoside pairs dA:dG, dC:::dG^oxo, dC:::dG, dA::dG^oxo in the aqueous phase have been taken into theoretical consideration. The influence of solvent relaxation on the above is also discussed. It can be concluded that the dA::dG^oxo nucleoside pair shows a lower ionization potential and higher electron affinity than the dA:dG pair in both a vertical and adiabatic mode. Therefore, it could be predicted, under electronic properties, that the electron ejected, for instance by a MutY 4[Fe-S]²⁺ cluster, is predisposed to trapping by the ds-DNA part containing the dA::dG^oxo pair rather than by dA::dG.

The Influence of Single, Tandem, and Clustered DNA Damage on the Electronic Properties of the Double Helix: A Theoretical Study

Oxidatively generated damage to DNA frequently appears in the human genome as the effect of aerobic metabolism or as the result of exposure to exogenous oxidizing agents, such as ionization radiation. In this paper, the electronic properties of single, tandem, and clustered DNA damage in comparison with native ds-DNA are discussed as a comparative analysis for the first time. A single lesion—8-oxo-7,8-dihydro-2′-deoxyguanosine (G^oxo), a tandem lesion—(5′S) and (5′R) 5′,8-cyclo-2′-deoxyadenosine (cdA), and the presence of both of them in one helix turn as clustered DNA damage were chosen and taken into consideration. The lowest vertical and adiabatic potential (VIP ~ 5.9 and AIP ~ 5.5 eV, respectively) were found for G^oxo, independently of the discussed DNA lesion type and their distribution within the double helix. Moreover, the VIP and AIP were assigned for ds-trimers, ds- dimers and single base pairs isolated from parental ds-hexamers in their neutral and cationic forms. The above results were confirmed by the charge and spin density population, which revealed that G^oxo can be considered as a cation radical point of destination independently of the DNA damage type (single, tandem, or clustered). Additionally, the different influences of cdA on the charge transfer rate were found and discussed in the context of tandem and clustered lesions. Because oligonucleotide lesions are effectively produced as a result of ionization factors, the presented data in this article might be valuable in developing a new scheme of anticancer radiotherapy efficiency.

Structural Characterization, DFT Calculation, NCI, Scan-Rate Analysis and Antifungal Activity against Botrytis cinerea of (E)-2-{[(2-Aminopyridin-2-yl)imino]-methyl}-4,6-di-tert-butylphenol (Pyridine Schiff Base)

Botrytis cinerea is a ubiquitous necrotrophic filamentous fungal phytopathogen that lacks host specificity and can affect more than 1000 different plant species. In this work, we explored L1 [(E)-2-{[(2-aminopyridin-2-yl)imino]-methyl}-4,6-di-tert-butylphenol], a pyridine Schiff base harboring an intramolecular bond (IHB), regarding their antifungal activity against Botrytis cinerea. Moreover, we present a full characterization of the L1 by NMR and powder diffraction, as well as UV–vis, in the presence of previously untested different organic solvents. Complementary time-dependent density functional theory (TD-DFT) calculations were performed, and the noncovalent interaction (NCI) index was determined. Moreover, we obtained a scan-rate study on cyclic voltammetry of L1. Finally, we tested the antifungal activity of L1 against two strains of Botrytis cinerea (B05.10, a standard laboratory strain; and A1, a wild type strains isolated from Chilean blueberries). We found that L1 acts as an efficient antifungal agent against Botrytis cinerea at 26 °C, even better than the commercial antifungal agent fenhexamid. Although the antifungal activity was also observed at 4 °C, the effect was less pronounced. These results show the high versatility of this kind of pyridine Schiff bases in biological applications.

Developing an effective polarizable bond method for small molecules with application to optimized molecular docking

Electrostatic interaction plays an essential role in protein-ligand binding. Due to the polarization effect, electrostatic interactions are largely impacted by their local environments. However, traditional force fields use fixed point charge-charge interactions to describe electrostatic interactions but is unable to include the polarization effect. The lack of the polarization effect in the force field representation can result in substantial error in biomolecular studies, such as molecular dynamics and molecular docking. Docking programs usually employ traditional force fields to estimate the binding energy between a ligand and a protein for pose selection or scoring. The intermolecular interaction energy mainly consists of van der Waals and electrostatic interaction in the force field representation. In the current study, we developed an Effective Polarizable Bond (EPB) method for small organic molecules and applied this EPB method to optimize protein-ligand docking in computational tests for a variety of protein-ligand systems. We tested the method on a set of 38 cocrystallized structures taken from the Protein Data Bank (PDB) and found that the maximum error was reduced from 7.98 Å to 2.03 Å when using EPB Dock, providing strong evidence that the use of EPB charges is important. We found that our optimized docking approach with EPB charges could improve the docking performance, sometimes dramatically, and the maximum error was reduced from 12.88 Å to 1.57 Å in Optimized Docking (in the case of 1fqx). The average RMSD decreased from 2.83 Å to 1.85 Å. Further investigations showed that the use of the EBP method could enhance intermolecular hydrogen bonding, which is a major contributing factor to improved docking performance. Developed tools for the calculation of the polarized ligand charge from a protein-ligand complex structure with the EPB method are freely available on GitHub (https://github.com/Xundrug/EPB).

In silico studies of the inhibition mechanism of dengue with papain

Dengue virus is becoming a major global disease; the envelope protein is the major target for vaccine development against Dengue. Nowadays, the attention has focused on developing inhibitors based on Papain is a promising target for treating Dengue. In the present work, the theoretical studies of E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and His159) complexes are analysed by Density Functional Theory (M06-2X/cc-pVDZ) method. Among the E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and Hys159) complexes, E-protein(Glu76)…Papain(Cys25) complex has the highest interaction value of -352.22 kcal/mol. Moreover, the natural bond orbital analysis also supports the above results. The 100 ns Molecular Dynamics simulation reveals that, E-protein(Ala54-Ile129)…Papain(Cys25) complex had the lowest root mean square deviation value of 1 Å compared to the E-protein(Ala54-Ile129)… Papain(Asn175 & His159) complexes. The salt bridge formation between the Asp103 and Lys110 residues are the important stabilizing factor in E-protein(Ala54-Ile129)…Papain(Cys25) complex. This result can extend our knowledge of the functional behaviour of Papain and provides structural insight to target Envelope protein as forthcoming drug targets in Dengue.

Methodologies in Spectral Tuning of DSSC Chromophores through Rational Design and Chemical-Structure Engineering

The investigation of new photosensitizers for Grätzel-type organic dye-sensitized solar cells (DSSCs) remains a topic of interest for researchers of alternative solar cell materials. Over the past 20 years, considerable and increasing research efforts have been devoted to the design and synthesis of new materials, based on "donor, π-conjugated bridge, acceptor" (D-π-A) organic dye photosensitizers. In this paper, the computational chemistry methods are outlined and the design of organic sensitizers (compounds, dyes) is discussed. With reference to recent literature reports, rational molecular design is demonstrated as an effective process to study structure-property relationships. Examples from established organic dye sensitizer structures, such as TA-St-CA, Carbz-PAHTDDT (S9), and metalloporphyrin (PZn-EDOT), are used as reference structures for an examination of this concept applied to generate systematically modified structural derivatives and hence new photosensitizers (i.e., dyes). Using computer-aided rational design (CARD), the in silico design of new chromophores targeted an improvement in spectral properties via the tuning of electronic structures by substitution of molecular fragments, as evaluated by the calculation of absorption profiles. This mini review provides important rational design strategies for engineering new organic light-absorbing compounds towards improved spectral absorption and related optoelectronic properties of chromophores for photovoltaic applications, including the dye-sensitized solar cell (DSSC).

Functional and Basis Set Dependence for Time-Dependent Density Functional Theory Trajectory Surface Hopping Molecular Dynamics: Cis-Azobenzene Photoisomerization

Within three functionals (TD-B3LYP, TD-BHandHLYP, and TD-CAM-B3LYP) in combination with four basis sets (3-21g, 6-31g, 6-31g(d), and cc-pvdz), global switching (GS) trajectory surface hopping molecular dynamics has been performed for cis-to-trans azobenzene photoisomerization up to the S₁ (nπ*) excitation. Although all the combinations show artificial double-cone structure of conical intersection between ground and first excited states, simulated quantum yields and lifetimes are in good agreement with one another; 0.6 (±5%) and 40.5 fs (±10%) by TD-B3LYP, 0.5 (±10%) and 35.5 fs (±4%) by TD-BHandHLYP, and 0.44 (±9%) and 35.2 fs (±10%) by TD-CAM-B3LYP. By analyzing distributions of excited-state population decays, hopping spots, and typical trajectories with performance of 12 functional/basis set combinations, it has been concluded that functional dependence for given basis set is slightly more sensitive than basis set dependence for given functional. The present GS on-the-fly time-dependent density functional theory (TDDFT) trajectory surface hopping simulation can provide practical benchmark guidelines for conical intersection driven excited-state molecular dynamics simulation involving in large complex system within ordinary TDDFT framework. © 2019 Wiley Periodicals, Inc.

The Driving Force for the Acylation of β-Lactam Antibiotics by L,D-Transpeptidase 2: Quantum Mechanics/Molecular Mechanics (QM/MM) Study

β-lactam antibiotics, which are used to treat infectious diseases, are currently the most widely used class of antibiotics. This study focused on the chemical reactivity of five- and six-membered ring systems attached to the β-lactam ring. The ring strain energy (RSE), force constant (FC) of amide (C-N), acylation transition states and second-order perturbation stabilization energies of 13 basic structural units of β-lactam derivatives were computed using the M06-2X and G3/B3LYP multistep method. In the ring strain calculations, an isodesmic reaction scheme was used to obtain the total energies. RSE is relatively greater in the five-(1a-2c) compared to the six-membered ring systems except for 4b, which gives a RSE that is comparable to five-membered ring lactams. These variations were also observed in the calculated inter-atomic amide bond distances (C-N), which is why the six-membered ring lactams C-N bond are more rigid than those with five-membered ring lactams. The calculated ΔG^# values from the acylation reaction of the lactams (involving the S-H group of the cysteine active residue from L,D transpeptidase 2) revealed a faster rate of C-N cleavage in the five-membered ring lactams especially in the 1-2 derivatives (17.58 kcal mol^-1 ). This observation is also reflected in the calculated amide bond force constant (1.26 mDyn/A) indicating a weaker bond strength, suggesting that electronic factors (electron delocalization) play more of a role on reactivity of the β-lactam ring, than ring strain.