Hidden Hemibonding in the Aqueous Hydroxyl Radical

The existence of a two-center, three-electron hemibond in the first solvation shell of ^•OH(aq) has long been a matter of debate. The hemibond manifests in ab initio molecular dynamics simulations as a small-r feature in the oxygen radial distribution function (RDF) for H₂O···^•OH, but that feature disappears when semilocal density functionals are replaced with hybrids, suggesting a self-interaction artifact. Using periodic simulations at the PBE0+D3 level, we demonstrate that the hemibond is actually still present (as evidenced by delocalization of the spin density) but is obscured by the hydrogen-bonded feature in the RDF due to a slight elongation of the hemibond. Computed electronic spectra for ^•OH(aq) are in excellent agreement with experiment and confirm that hemibond-like configurations play an outsized role in the spectroscopy due to an intense charge-transfer transition that is strongly attenuated in hydrogen-bonded configurations. Apparently, 25% exact exchange (as in PBE0) is insufficient to eliminate delocalization of unpaired spins.

Infrared Spectroscopy and Anharmonic Vibrational Analysis of (H2O-Krn)+ (n = 1-3): Hemibond Formation of the Water Radical Cation

The hemibond is a nonclassical covalent bond formed between a radical (cation) and a closed shell molecule. The hemibond formation ability of water has attracted great interest, concerning its role in ionization of water. While many computational studies on the water hemibond have been performed, clear experimental evidence has been hardly reported because the hydrogen bond formation overwhelms the hemibond formation. In the present study, infrared photodissociation spectroscopy is applied to (H₂O-Kr_n)⁺ (n = 1-3) radical cation clusters. The observed spectra of (H₂O-Kr_n)⁺ are well reproduced by the anharmonic vibrational simulations based on the hemibonded isomer structures. The firm evidence of the hemibond formation ability of water is revealed.

Shape-based Machine Learning Models for the Potential Novel COVID-19 Protease Inhibitors Assisted by Molecular Dynamics Simulation

BACKGROUND

The vast geographical expansion of novel coronavirus and an increasing number of COVID-19 affected cases have overwhelmed health and public health services. Artificial Intelligence (AI) and Machine Learning (ML) algorithms have extended their major role in tracking disease patterns, and in identifying possible treatments.

OBJECTIVE

This study aims to identify potential COVID-19 protease inhibitors through shape-based Machine Learning assisted by Molecular Docking and Molecular Dynamics simulations.

METHODS

31 Repurposed compounds have been selected targeting the main coronavirus protease (6LU7) and a machine learning approach was employed to generate shape-based molecules starting from the 3D shape to the pharmacophoric features of their seed compound. Ligand-Receptor Docking was performed with Optimized Potential for Liquid Simulations (OPLS) algorithms to identify highaffinity compounds from the list of selected candidates for 6LU7, which were subjected to Molecular Dynamic Simulations followed by ADMET studies and other analyses.

RESULTS

Shape-based Machine learning reported remdesivir, valrubicin, aprepitant, and fulvestrant as the best therapeutic agents with the highest affinity for the target protein. Among the best shape-based compounds, a novel compound identified was not indexed in any chemical databases (PubChem, Zinc, or ChEMBL). Hence, the novel compound was named 'nCorv-EMBS'. Further, toxicity analysis showed nCorv-EMBS to be suitable for further consideration as the main protease inhibitor in COVID-19.

CONCLUSION

Effective ACE-II, GAK, AAK1, and protease 3C blockers can serve as a novel therapeutic approach to block the binding and attachment of the main COVID-19 protease (PDB ID: 6LU7) to the host cell and thus inhibit the infection at AT2 receptors in the lung. The novel compound nCorv- EMBS herein proposed stands as a promising inhibitor to be evaluated further for COVID-19 treatment.

Ab Initio Study of Ionized Water Radical Cation (H2O)8+ in Combination with the Particle Swarm Optimization Method

The structures of cationic water clusters (H₂O)₈⁺ have been globally explored by the particle swarm optimization method in combination with quantum chemical calculations. Geometry optimization and vibrational analysis for the 15 most interesting clusters were computed at the MP2/aug-cc-pVDZ level and infrared spectrum calculation at MPW1K/6-311++G** level. Special attention was paid to the relationships between their configurations and energies. Both MP2 and B3LYP-D3 calculations revealed that the cage-like structure is the most stable, which is different from a five-membered ring lowest energy structure but agrees well with a cage-like structure in the literature. Furthermore, our obtained cage-like structure is more stable by 0.87 and 1.23 kcal/mol than the previously reported structures at MP2 and B3LYP-D3 levels, respectively. Interestingly, on the basis of their relative Gibbs free energies and the temperature dependence of populations, the cage-like structure predominates only at very low temperatures, and the most dominating species transforms into a newfound four-membered ring structure from 100 to 400 K, which can contribute greatly to the experimental infrared spectrum. By topological analysis and reduced density gradient analysis, we also investigated the structural characteristics and bonding strengths of these water cluster radical cations.

Ab initio investigation of structure, stability, thermal behavior, bonding, and infrared spectra of ionized water cluster (H2O)6. J Chem Phys 2016;145:154307. [PMID: 27782468 DOI: 10.1063/1.4964860]

The low-lying isomers of cationic water cluster (H₂O)₆⁺ have been globally explored by using particle swarm optimization algorithm in conjunction with quantum chemical calculations. Compared with previous results, our searching method covers a wide range of structural isomers of (H₂O)₆⁺ and therefore turns out to be more effective. With these local minima, geometry optimization and vibrational analysis are performed for the most interesting clusters at second-order Møller-Plesset (MP2)/aug-cc-pVDZ level, and their energies are further refined at MP2/aug-cc-pVTZ and coupled-cluster theory with single, double, and perturbative triple excitations/aug-cc-pVDZ level. The interaction energies using the complete basis set limits at MP2 level are also reported. The relationships between their structure arrangement and their energies are discussed. Based on the results of thermal simulation, structural change from a four-numbered ring to a tree-like structure occurs at T ≈ 45 K, and the relative population of six lowest-free-energy isomers is found to exceed 4% at some point within the studied temperature range. Studies reveal that, among these six isomers, two new-found isomers constitute 10% of isomer population at 180 K, and the experimental spectra can be better explained with inclusions of the two isomers. The molecular orbitals for six representative cationic water clusters are also studied. Through topological and reduced density gradient analysis, we investigated the structural characteristics and the bonding strengths of these water cluster radical cations.

Potential Paths for the Hydrogen-Bond Relaxing With (H2O)N Cluster Size

Relaxation of the inter- and intra-molecular interactions for the hydrogen bond (O:H-O) between undercoordinated molecules determines the unusual behavior of water nanodroplets and nanobubbles. However, probing such potentials remains unreality. Here we show that the Lagrangian solution [Huang et al., J. Phys. Chem. B, 2013. 117: 13639] transforms the observed H-O bond (x = H) and O:H nonbond (x = L) lengths and phonon frequencies (dx, x) [Sun et al., J. Phys. Chem. Lett., 2013. 4: 2565] into the respective force constants and bond energies (kx, Ex) and hence enables the mapping of the potential paths for the O:H-O bond relaxing with water cluster size. Results show that molecular undercoordination not only reduces the molecular size (dH) with enhanced H-O energy from the bulk value of 3.97 to 5.10 eV for a H2O monomer, but also enlarges the molecular separation (dL) with reduced O:H energy from 95 to 35 meV for a dimer. The H-O energy gain raises the melting point from bulk value 273 to 310 K for the skin and the O:H energy loss lowers the freezing temperature from bulk value 258 to 202 K for 1.4 nm sized droplet, by dispersing the quasisolid phase boundaries.

Proton transfer rates in ionized water clusters (H2O)n (n = 2–4)

A proton transfer process is usually dominant in several biological phenomena such as the energy relaxation of photo-excited DNA base pairs and a charge relay process in Ser-His-Glu.

Structural evolution and solvation of the OH radical in ionized water radical cations (H2O)n+, n = 5–8

Structural evolution of ionized water radical cations (H₂O)_n⁺, n = 5–8, is studied by ab intio methods.

Water desalination through armchair carbon nanotubes: a molecular dynamics study

Separation of ions from water using armchair carbon nanotubes.

Isomer-selective infrared spectroscopy of the cationic trimethylamine dimer to reveal its charge sharing and enhanced acidity of the methyl groups

The isomer-selective infrared spectroscopy revealed the charge-shared (hemibond) and the C⋯HN hydrogen-bond structures of the trimethylamine dimer cation.

Size, separation, structural order, and mass density of molecules packing in water and ice

The structural symmetry and molecular separation in water and ice remain uncertain. We present herewith a solution to unifying the density, the structure order and symmetry, the size (H-O length dH), and the separation (d(OO) = d(L) + d(H) or the O:H length d(L)) of molecules packing in water and ice in terms of statistic mean. This solution reconciles: i) the d(L) and the d(H) symmetrization of the O:H-O bond in compressed ice, ii) the d(OO) relaxation of cooling water and ice and, iii) the d(OO) expansion of a dimer and between molecules at water surface. With any one of the d(OO), the density ρ(g·cm⁻³), the d(L), and the d(H), as a known input, one can resolve the rest quantities using this solution that is probing conditions or methods independent. We clarified that: i) liquid water prefers statistically the mono-phase of tetrahedrally-coordinated structure with fluctuation, ii) the low-density phase (supersolid phase as it is strongly polarized with even lower density) exists only in regions consisting molecules with fewer than four neighbors and, iii) repulsion between electron pairs on adjacent oxygen atoms dictates the cooperative relaxation of the segmented O:H-O bond, which is responsible for the performance of water and ice.