Nonempirical Prediction of the Relative Electrospray Ionization Efficiencies of Nitroanilines by Combined CBS-QB3 and SCC-DFTB Calculations.

Computational study for three isomeric nitroanilines was carried out to find a nonempirical model to predict their relative ionization efficiencies compared to the corresponding experimental logIE values. The CBS-QB3 method was used for the gas-phase calculations of protonated nitroanilines, and the SCC-DFTB method for the calculations of spherical droplets, containing one protonated nitroaniline molecule, one hydronium cation, and 48 acetonitrile molecules, randomly generated by the PACKMOL program. The obtained results show that neither the gas-phase energy of protonated nitroaniline (or gas-phase basicity) nor the average energy of droplets can be used to predict the ionization efficiency ordering of those isomers. However, the difference of gas-phase and droplet energies gives a good correlation with logIE values and, thus, can be used for the prediction of relative ionization efficiencies.

1,3-Thiazolbenzamide Derivatives as Chikungunya Virus nsP2 Protease Inhibitors.

Chikungunya fever results from an infection with Chikungunya virus (CHIKV, genus Alphavirus) that is prevalent in tropical regions and is spreading fast to temperate climates with documented outbreaks in Europe and the Americas. Currently, there are no available vaccines or antiviral drugs for prevention or treatment of Chikungunya fever. The nonstructural proteins (nsPs) of CHIKV responsible for virus replication are promising targets for the development of new antivirals. This study was attempted to find out new potential inhibitors of CHIKV nsP2 protease using the ligand-based drug design. Two compounds 10 and 10c, identified by molecular docking, showed antiviral activity against CHIKV with IC50 of 13.1 and 8.3 µM, respectively. Both compounds demonstrated the ability to inhibit the activity of nsP2 in a cell-free assay, and the impact of compound 10 on virus replication was confirmed by western blot. The molecular dynamics study of the interactions of compounds 10 and 10c with CHIKV nsP2 showed that a possible mechanism of action of these compounds is the blocking of the active site and the catalytic dyad of nsP2.

Novel Analogues of the Chikungunya Virus Protease Inhibitor: Molecular Design, Synthesis, and Biological Evaluation.

The Chikungunya virus (CHIKV) is an arbovirus belonging to the genus Alphavirus of the Togaviridae family. CHIKV is transmitted by the mosquitoes and causes Chikungunya fever. CHIKV outbreaks have occurred in Africa, Asia, Europe, and the countries of Indian and Pacific Oceans. In 2013, CHIKV cases were registered for the first time in the Americas on the Caribbean islands. There is currently no vaccine to prevent or medicines to treat CHIKV infection. The CHIKV nonstructural protease (nsP2) is a promising potential target for the development of drugs against CHIKV infection because this protein is one of the key components of the viral replication complex and is involved in multiple steps of virus infection. In this work, novel analogues of the potential CHIKV nsP2 protease inhibitor, first reported by Das et al. in 2016, were identified using molecular modeling methods, synthesized, and evaluated in vitro. The optimization of the structure of the inhibitor allowed to increase the antiviral activity of the compound 2-10 times. The possible mechanism of action of the identified potential inhibitors of the CHIKV nsP2 protease was studied in detail using molecular dynamics (MD) simulations. According to the MD results, the most probable mechanism of action is the blocking of conformational changes in the nsP2 protease required for substrate recognition and binding.

Computational study on the reactivity of tetrazines toward organometallic reagents.

The possible reaction pathways between methyllithium and disubstituted 1,2,4,5-tetrazines (bearing methyl, methylthio, phenyl, and 3,5-dimethylpyrazolyl groups) were investigated by means of the density functional theory B3LYP/6-31G* method. Solvation was modeled using the supermolecule approach, adding one tetrahydrofuran molecule to the complexes. Comparison of the calculated energies and structures for the alternate azaphilic and nucleophilic addition pathways showed that the azaphilic addition is kinetically favored over nucleophilic addition, while thermodynamically the nucleophilic addition is usually preferred. The coordination of the tetrazine molecule with methyllithium was found to play a crucial role in the process. These findings provide the first rationale for the experimentally observed unique reactivity of tetrazines toward polar organometallic reagents, suggesting the presence of a kinetically controlled process.

Computational study of cesium cation interactions with neutral and anionic compounds related to soil organic matter.

The gas-phase cesium cation affinities (CsCAs) and basicities (CsCBs) for 56 simple neutral compounds (mostly aromatic molecules) and 41 anions (carboxylates and phenolates) were calculated using density functional theory (DFT), in the context of the interaction of Cs(+) with soil organic matter (SOM). The B3LYP/def2-TZVP method gives in general CsCAs and CsCBs in a good agreement with experimental data. The strong deviations in case of NO(3)(-) and CsSO(4)(-) anions need further experimental investigations as the high-level CCSD(T) calculations support B3LYP results. Different cesium cation complexation patterns between Cs(+) and the neutral and anionic systems are discussed. As expected, the strongest CsCAs are observed for anions. The corresponding quantities are approximately by 4-5 times higher than for the neutral counterparts, being in the range 90-118 kcal/mol. The weakest cesium cation bonding is observed in the case of unsubstituted aromatic systems (11-15 kcal/mol).

Small-Molecule Ligands as Potential GDNF Family Receptor Agonists.

To find out potential GDNF family receptor α1 (GFRα1) agonists, small molecules were built up by molecular fragments according to the structure-based drug design approach. Molecular docking was used to identify their binding modes to the biological target GFRα1 in GDNF-binding pocket. Thereafter, commercially available compounds based on the best predicted structures were searched from ZINC and MolPort databases (similarity ≥ 80%). Five compounds from the ZINC library were tested in phosphorylation and luciferase assays to study their ability to activate GFRα1-RET. A bidental compound with two carboxyl groups showed the highest activity in molecular modeling and biological studies. However, the relative position of these groups was important. The meta-substituted structure otherwise identical to the most active compound 2-[4-(5-carboxy-1H-1,3-benzodiazol-2-yl)phenyl]-1H-1,3-benzodiazole-5-carboxylic acid was inactive. A weaker activity was detected for a compound with a single carboxyl group, that is, 4-(1,3-benzoxazol-2-yl)benzoic acid. The substitution of the carboxyl group by the amino or acetamido group also led to the loss of the activity.

Molecular Dynamics Simulations of the Interactions between Glial Cell Line-Derived Neurotrophic Factor Family Receptor GFRα1 and Small-Molecule Ligands.

The glial cell line-derived neurotrophic factor (GDNF) family ligands (GFLs) support the survival and functioning of various neuronal populations. Thus, they could be attractive therapeutic agents against a multitude of neurodegenerative diseases caused by progressive death of GFLs responsive neurons. Small-molecule ligands BT13 and BT18 show an effect on GDNF family receptor GFRα1 and RET receptor tyrosine kinase RetA function. Thus, their potential binding sites and interactions were explored in the GDNF-GFRα1-RetA complex using molecular docking calculations as well as molecular dynamics (MD) simulations. Three possible regions were examined: the interface between GDNF and GFRα1 (region A), the RetA interface with GFRα1 (region B), and a possible allosteric site in GFRα1 (region C). The results obtained by the docking calculations and the MD simulations indicate that the preferable binding occurs at the allosteric site. A less preferable binding site was detected on the RetA surface interfacing GFRα1. In the membrane-bound state of RetA this can enable compounds BT13 and BT18 to act as direct RetA agonists. The analysis of the MD simulations shows hydrogen bonds for BT13 and significant hydrophobic interactions with GFRα1 for BT13 and BT18 at the allosteric site.

Indole-like Trk receptor antagonists.

The virtual screening for new scaffolds for TrkA receptor antagonists resulted in potential low molecular weight drug candidates for the treatment of neuropathic pain and cancer. In particular, the compound (Z)-3-((5-methoxy-1H-indol-3-yl)methylene)-2-oxindole and its derivatives were assessed for their inhibitory activity against Trk receptors. The IC50 values were computationally predicted in combination of molecular and fragment-based QSAR. Thereafter, based on the structure-activity relationships (SAR), a series of new compounds were designed and synthesized. Among the final selection of 13 compounds, (Z)-3-((5-methoxy-1-methyl-1H-indol-3-yl)methylene)-N-methyl-2-oxindole-5-sulfonamide showed the best TrkA inhibitory activity using both biochemical and cellular assays and (Z)-3-((5-methoxy-1-methyl-1H-indol-3-yl)methylene)-2-oxindole-5-sulfonamide was the most potent inhibitor of TrkB and TrkC.

Interaction between the Cesium Cation and Cesium Carboxylates: An Extended Cs+ Basicity Scale.

The interaction between the cesium cation and fulvic or humic acids is supposed to play a role in cesium mobility in the environment, which is of importance in the context of geographical dispersion or concentration of the corresponding radionuclides. Among the singly charged positive clusters generated by electrospray ionization of mixtures of carboxylic acids (AH) and cesium salts (nitrate, iodide, or trifluoroacetate), the cluster [A- Cs+ ]Cs+ [A'- Cs+ ] was subjected to collision-induced dissociation. The fragmentation into Cs+ [A- Cs+ ] and Cs+ [A'- Cs+ ] was treated using the kinetic method. A gas-phase basicity ladder was built by a step-by-step addition of the relative basicity data. The relative basicity scale deduced from the kinetic method was calibrated using the affinity and basicity (the enthalpy and Gibbs energy scales, respectively) obtained from DFT calculations. The enhanced basicity of the polyacid salts, relative to their monoacid salts, as well as the substituent effects on aliphatic and aromatic structures, are discussed.

Computational study of the Sonogashira cross-coupling reaction in the gas phase and in dichloromethane solution.

The Sonogashira cross-couplig reaction, consisting of oxidative addition, cis-trans isomerization, transmetalation, and reductive elimination, was computationally modeled using the DFT B3LYP/cc-pVDZ method for reaction between bromobenzene and phenylacetylene. Palladium diphosphane was used as a catalyst, copper(I) bromide as a co-catalyst and trimethylamine as a base. The reaction mechanism was studied both in the gas phase and in dichloromethane solution using PCM method. The complete catalytic cycle is thermodynamically strongly shifted toward products (diphenylacetylene and regenerated palladium catalyst) and is exothermic being in accordance with experimental data. The rate-determining step is the oxidative addition, since the highest point on the Gibbs energy graph of the complete reaction is the transition state of this step. This conclusion is also supported by recent experimental data. The computed energy profile suggests that the transmetalation step is initiated by the dissociation of neutral ligand, while the activation Gibbs energy of this step is 0.1 kcal mol(-1) in the gas phase.

A study of the cesium cation bonding to carboxylate anions by the kinetic method and quantum chemical calculations.

Collision-induced dissociation (CID) of the Cs(+) heterodimer adducts of the nitrate anion (NO(3)(-)) and a variety of substituted benzoates (XBenz(-)) [(XBenz(-))(Cs(+))(NO(3)(-))](-) produces essentially nitrate and benzoate ions. A plot of the natural logarithm of their intensity ratio, ln[I (NO(3)(-))/I(XBenz(-))], versus the calculated cesium cation affinity (DFT B3LYP) of the substituted benzoate ions (equivalent to the enthalpy of heterolytic dissociation of the salt) is reasonably linear. This suggests that the kinetic method can be used as a source of data on the intrinsic interaction between the anionic and the cationic moieties in a salt.

Interaction of the cesium cation with mono-, di-, and tricarboxylic acids in the gas phase. A Cs+ affinity scale for cesium carboxylates ion pairs.

Humic substances (HS), including humic and fulvic acids, play a significant role in the fate of metals in soils. The interaction of metal cations with HS occurs predominantly through the ionized (anionic) acidic functions. In the context of the effect of HS on transport of radioactive cesium isotopes in soils, a study of the interaction between the cesium cation and model carboxylic acids was undertaken. Structure and energetics of the adducts formed between Cs+ and cesium carboxylate salts [Cs+RCOO-] were studied by the kinetic method and density functional theory (DFT). Clusters generated by electrospray ionization mass spectrometry from mixtures of a cesium salt (nitrate, iodide, trifluoroacetate) and carboxylic acids were quantitatively studied by CID. By combining the results of the kinetic method and the energetic data from DFT calculations, a scale of cesium cation affinity, CsCA, was built for 33 cesium carboxylates representing the first scale of cation affinity of molecular salts. The structural effects on the CsCA values are discussed.

Bonding energetics in clusters formed by cesium salts: a study by collision-induced dissociation and density functional theory.

In relation to the interaction between (137)Cs and soil organic matter, electrospray mass spectrometry experiments and density functional theory (DFT) calculations were carried out on the dissociation of positively charged adducts formed by cesium nitrate and cesium organic salts attached to a cesium cation [Cs(CsNO(3))(CsA)](+) (A = benzoate, salicylate, hydrogen phthalate, hydrogen maleate, hydrogen fumarate, hydrogen oxalate, and hydrogen malonate ion). These mixed clusters were generated by electrospray from methanol solutions containing cesium nitrate and an organic acid. Collision-induced dissociation of [Cs(CsNO(3))(CsA)](+) in a quadrupole ion trap gave [Cs(CsNO(3))](+) and [Cs(CsA)](+) as major product ions. Loss of HNO(3) was observed, and also CO(2) loss in the case of A = hydrogen malonate. Branching ratios for the dissociation into [Cs(CsNO(3))](+) and [Cs(CsA)](+) were treated by the Cooks' kinetic method to obtain a quantitative order of bonding energetics (enthalpies and Gibbs free energies) between Cs(+) and the molecular salt (ion pair) CsA, and were correlated with the corresponding values calculated using DFT. The kinetic method leads to relative scales of Cs(+) affinities and basicities that are consistent with the DFT-calculated values. This study brings new data on the strong interaction between the cesium cation and molecular salts CsA.

Investigations of cluster ions formed between cesium cations and benzoic, salicylic and phthalic acids by electrospray mass spectrometry and density-functional theory calculations. Toward a modeling of the interaction of Cs+ with humic substances.

A concerted theoretical (density-functional theory) and experimental electrospray mass spectrometry study was conducted on the formation of cesium cation adducts with small molecules taken as models of specific interactions sites in humic substances. Electrospray experiments with phenol, benzoic acid, salicylic acid, and phthalic acid, in methanolic solution containing cesium nitrate, were performed using a quadrupole ion trap. The formation of positively charged mixed clusters, [Cs(CsNO3)(n)(CsA1)(m)(Cs2A2)(p)]+ (A1 = benzoate, salicylate, and hydrogenophthalate, A2 = phthalate), was observed. Calculations of structures and bonding energetics of Cs+ in simple adducts formed with NO3-, CsNO3, A-, AH, and CsA are reported. The observation of variable cluster stoichiometry (n, m and p values) was interpreted in terms of more or less favorable interaction energies between Cs+ and the neutral species constituting the clusters. Phenol did not form clusters in significant abundances, despite a strong calculated interaction between Cs+ and cesium phenolate. This was attributed to its weak acid dissociation in the electrospray solution.