Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.

Numerical optimization is a common technique in various areas of computational chemistry, molecular modeling and drug design. It is a key element of 3D techniques, for example, the optimization of protein-ligand poses and small-molecule conformers. Here, often the BFGS algorithm or variants thereof are used. However, the BFGS algorithm tends to make unreasonable large changes to the optimized system under certain circumstances. This behavior has been known for a long time and different solutions have been suggested. Recently, we have analyzed the optimization behavior of our novel JAMDA scoring function in detail and proposed the limited step length (LSL)-BFGS algorithm as a new solution to the problem of excessively large steps during optimization. The LSL-BFGS algorithm allows to control the step sizes during optimization. Its unique feature is the inclusion of arbitrary domain knowledge into the selection of the step sizes. Here, we introduce the open-source LSLOpt C++ library that implements this LSL-BFGS algorithm and demonstrate its usage.

Consistent basis sets of triple-zeta valence with polarization quality for the elements Cs-Po were derived for periodic quantum-chemical solid-state calculations. They are an extension of the pob-TZVP-rev2 [Vilela Oliveira, D.; Laun, J.; Peintinger, M. F. and Bredow, T., J. Comput. Chem., 2019, 40 (27), 2364-2376] basis sets and are based on the fully relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and on the def2-TZVP valence basis of the Ahlrichs group. The basis sets are constructed to minimize the basis set superposition error (BSSE) in crystalline systems. The contraction scheme, the orbital exponents, and contraction coefficients were optimized in order to ensure robust and stable self-consistent-field (SCF) convergence for a set of compounds and metals. For the applied PW1PW hybrid functional, the average deviations of the calculated lattice constants from experimental references are smaller with pob-TZVP-rev2 than with standard basis sets available from the CRYSTAL basis set database.

Computation of the thermodynamic consequences of protein mutations holds great promise in protein biophysics and design. Alchemical free energy methods can give improved estimates of mutational free energies, and are already widely used in calculations of relative and absolute binding free energies in small molecule design problems. In principle, alchemical methods can address any amino acid mutation with an appropriate alchemical pathway, but identifying a strategy that produces such a path for proline and glycine mutations is an ongoing challenge. Most current strategies perturb only side chain atoms, while proline and glycine mutations also alter the backbone parameters and backbone ring topology. Some strategies also perturb backbone parameters and enable glycine mutations. This work presents a strategy that enables both proline and glycine mutations and comprises two key elements: a dual backbone with restraints and scaling of bonded terms, facilitating backbone parameter changes, and a soft bond in the proline ring, enabling ring topology changes in proline mutations. These elements also have utility for core hopping and macrocycle studies in computer-aided drug design. This new strategy shows slight improvements over an alternative side chain perturbation strategy for a set T4 lysozyme mutations lacking proline and glycine, and yields good agreement with experiment for a set of T4 lysozyme proline and glycine mutations not previously studied. To our knowledge this is the first report comparing alchemical predictions of proline mutations with experiment. With this strategy in hand, alchemical methods now have access to the full palette of amino acid mutations.

In the era of exascale supercomputers, large-scale, and long-time molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)-communication in the MPI-parallelized fast multipole method (FMM) combined with MD calculations under three-dimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimum-transferred data (MTD) method, should enable large-scale and long-time MD calculations to be calculated efficiently, under the condition of massive MPI-parallelization on exascale supercomputers.

The combination of a Monte Carlo (MC) sampling of the configurational space with time dependent-density functional theory (TD-DFT) to estimate vertical excitations energies has been applied to compute the absorption spectra of a family of merocyanine dyes in both their monomeric and dimeric forms. These results have been compared to those obtained using a static DFT/TD-DFT approach as well as to the available experimental spectra. Though suffering of the limitations related to the use of DFT and TD-DFT for this type of systems, our data clearly show that the classical MC sampling provides a suitable alternative to classical molecular dynamics to explore the structural flexibility of these donor-acceptor (D-π-A) chromophores enabling a realistic description of the potential energy surface of both their monomers and aggregates (here dimers) and thus of their spectra. Overall, the combination of MC sampling with quantum mechanics (TD-DFT) calculations, carried out in implicit dioxane solvent on random snapshots, provides a workable compromise to solve the combined challenge of accuracy and time-consuming problem not only for merocyanines momers, but also for their dimers, up to now less investigated. Indeed, the simulated absorption spectra fairly agree with the experimental ones, suggesting the general reliability of the method.

Ligand electrochemical parameters, EL , more commonly known as Lever parameters, have played a major research role in understanding redox processes involved in inorganic electrochemistry, enzymatic reactions, catalysis, solar cells, biochemistry, and materials science. Despite their broad usefulness, Lever parameters are not well understood at a first-principles level. Using density functional theory, we demonstrate in this contribution that a ligand's Lever parameter is fundamentally related to the ligand's ability to alter the eigenvalue of the electroactive spin-orbital in an octahedral transition metal complex. Our analysis furthers a first-principles understanding of the nature of Lever parameters.

In the present study, the Divide and Conquer MBAR (DC-MBAR) method is proposed to predict the free energies based on the data sampled by multi-states simulations. For DC-MBAR method, the overlap between any two alchemical states is calculated first and those with sufficient overlap are defined as the adjacent states. Unlike the traditional MBAR method, which calculates the free energy of each state using all the data at once, DC-MBAR focuses on predicting the free energy changes between adjacent states. To estimate the free energy changes accurately, the other states with overlaps with the two adjacent states bigger than the defined threshold are included in the MBAR equation. At a specific threshold, the free energies predicted by DC-MBAR are very close to those calculated by the traditional MABR method. Furthermore, DC-MBAR scheme can reduce both the computation and memory cost. One important characteristic of DC-MBAR method is linear scaling, which means the CPU time with the change of the number of states is a straight-line relation. As the pair-based calculations are mutually independent and parallelizable, all accessible CPU cores on the HPC cluster could be utilized, which makes DC-MBAR strategy more efficient.

There have been significant advances in the biological use of hypervalent selenium and tellurium compounds as cysteine protease inhibitors. However, the full understanding of their reaction mechanisms for and cysteine proteases inhibition is still elusive. Kinetic studies suggest an irreversible inhibition mechanism, which was explained by forming a covalent bond between the enzyme sulfhydryl group and the chalcogen atom at its hypervalent state (+4). In this work, we performed a theoretical investigation using density functional theory to propose the active inhibitor form in an aqueous solution. To this end, we investigated chloride ligand exchange reactions by oxygen and sulfur nucleophiles on hypervalent selenium and tellurium compounds. All tetra- and tri-coordinated chalcogen compounds and distinct protonation states of the nucleophiles were considered, totaling 34 unique species, 7 nucleophiles, and 155 free energies reactions. We discovered that chloride is easily replaced by a nonprotonated nucleophile (SH- or OH- ) in R2 SeCl2 . We also found that tri-coordinate species are more stable than their tetra-coordinate counterparts, with selenoxide (R2 SeO) protonation being strongly exergonic in acid pH. The thermodynamic and kinetic results suggest that the protonated selenoxide (R2 SeOH+ ) is the most probable active chemical species in biological media. The computed energetic profiles paint a possible picture for selenuranes activity, with successive exergonic steps leading to a covalent inhibition of thiol-dependent enzymes, like cysteine proteases. A second pathway has also been uncovered, with a direct reaction to chalcogenonium cation (R2 SeCl+ ) as the inhibition step. Tellurium compounds showed similar trends but formed telluroxide in a pH-independent fashion.

We propose an efficient and accurate non-iterative method, dubbed uESE, for calculating solvation free energies. Apart from a COSMO-like electrostatic term, the model takes into account non-electrostatic contributions, which depend on atomic surfaces, induced surface charge densities, and the molecular volume. uESE is tested on 35 polar and 57 non-polar solvents. The calculated and experimental solvation free energies are compared for 2892 systems. The method exhibits an excellent performance, which is superior to major solvation methods. The mean absolute error of predicted solvation energies is found below 1 kcal/mol for neutral solutes and below 3 kcal/mol for ions. The calculated data are almost independent of the quantum-chemical method or/and basis sets employed.

We present a generalization of a self-consistent electrostatic embedding approach (SC-Ewald) devised to investigate the photophysical properties of 3D periodic materials, to systems in one- or two-dimensional (2D) reduced periodicity. In this approach, calculations are carried out on a small finite molecular cluster extracted from a periodic model, while the crystalline environment is accounted for by an array of point charges which are fitted to reproduce the exact electrostatic potential (at ground or the excited state) of the infinite periodic system. Periodic density functional theory (DFT) calculations are combined with time dependent DFT calculations to simulate absorption and emission properties of the extended system under investigation. We apply this method to compute the UV-Vis. spectra of bulk and quantum-confined 0D quantum dots and 2D extended nanoplatelets of CdSe, due to their relevance as sensitizers in solar cells technologies. The influence of the size and shape of the finite cluster model chosen in the excited state calculations was also investigated and revealed that, although the long-range electrostatics of the environment are important for the calculation of the UV-Vis, a subtle balance between short- and long-range effects exists. These encouraging results demonstrate that this self-consistent electrostatic embedding approach, when applied in different dimensions, can successfully model the photophysical properties of diverse material classes, making it an attractive low-cost alternative to far more computationally demanding electronic structure methods for excited state calculations.

The reactions of γ-dehydronitration of furaxanenitrolic acids have been studied within the density functional theory using molecular electron density theory scheme at the MPWB1K(PCM)/6-311G(d,p) level of theory. The alteration of bonding along the course of the reaction is studied in the topology of the electron density functional within the bonding evolution theory perspective. The characteristics of electron density changes indicate that we can distinguish six different phases in the nitrous acid extrusion from furaxanenitrolic acid 1a. These different phases related to the intrinsic reaction coordinate path of the analyzed reaction denote the non-concerted nature of the molecular mechanism.

The mechanism of 6-π electrocyclization of all-cis, mono-trans, and double-trans [10]annulene to yield 4a,8a-dihydronaphthalene has been explored at various quantum-chemical methods. The mono-trans configuration cyclizes preferentially to trans-4a,8a-dihydronaphthalene, in agreement with the experimental results. The cyclization of the all-cis configuration requires firstly a bond-shifting to the naphthalene-like conformation of double-trans [10]annulene, which is the rate-limiting step, and finally its azulene-like conformation electrocyclizes quickly to cis-4a,8a-dihydronaphthalene. Its experimental rate coefficient is consistent with the computed one for the cyclization of the all-cis configuration, unlike the calculated one for the double-trans configuration. These results confirm the configurations assigned by Masamune et al. to the two isomers which they isolated.

Some Mo-oxo complexes bearing pyridine rings have the capability for dihydrogen production from water. However, energy barrier and overall energy vary depending on the effect exerted by several substituent groups located at different positions around one or more pyridine rings which are ligands of these compounds. Based on the Karunadasa and coworkers investigation where the para-position was experimentally tested in compounds derivatised from the 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine oxo-molybdenum complex synthesized (Karunadasa et al., Nature, 2010, 464, 1329), we tested the combined effect of electron-withdrawing and electron-donating groups simulated as perturbations represented by point-charges. Then, we used the density polarization concept, δρ(r), a local reactivity descriptor corresponding to the partially integrated linear response function, χ(r, r') (a non-local reactivity descriptor), which is able to reveal different displacements of π-electrons on molecular structures. We perturbed the para-positions in the pentadentate ligand 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine in the Mo-based complex by means of point-charges. They were located in three different configurations of the organic ligand (trans, geminal, and cis) which could help to explain energy barriers and overall energy of reactions catalyzed by this type of Mo-complexes. Our results indicate that the trans configuration of point-charges induces the most amount of fraction of electron shifted on the complex. A Mo-based complex bearing the same trans configuration for electron-withdrawing and electron-donating substituent groups (cyano and amino, respectively), leads to a kinetically more favorable H2 release than the cis or geminal configuration of the substituent groups aforementioned.

This study presents the benchmark calculations of proton affinities (PAs) and gas-phase basicities (GBs) of 8-para substituted benzaldehyde compounds using the multilevel model chemistries (G3B3 and G4), density-functional quantum model (B3LYP) and ab initio model (MP2). The results show that the computed properties are strongly correlated with the available experimental data. The PAs and the GBs of other eight para-substituted benzaldehyde compounds, for which the experimental data does not currently exist, have been calculated using G3B3 and B3LYP methods. The correlations between the experimental PAs and GBs with the computed properties such as PA, GB, chemical properties (bond lengths, electron density and δ1 H NMR chemical shift) of the investigated benzaldehydes have been studied and statistically analyzed. The influence of the substituted groups has been discussed in terms of inductive effect and electron donating and electron withdrawing effect. The results obtained show that the chemical properties of the benzaldehyde compounds are controlled by the strong coupling between the CHO group and the nature of the para-substituent groups through the benzene ring as a conducting linkage.

The study of the water-splitting process, which can proceed in 2e- as well as 4e- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.

Quadruple bonding is uncommon for main group elements and the identification of species forming such bonds is remarkably interesting particularly in diatomic anions for which there is a lack of information. Here, it is found that the MX- anions, TcN- , RuC- , RhB- , and PdBe- , present quadruple bonding, as do the corresponding MX neutrals, even though a different type of σ2 bond is involved in ∑+ states of neutral and anions. Specifically, the ground states (X2 Δ or X2 ∑+ ) of the four anions and their first excited states (A2 ∑+ or A2 Δ) of TcN- , RuC- , and RhB- present quadruple bonds consisting of two σ and two π bonds: (4dz2  - 2pz )2 , 5pz 0  ← 2s2 , (4dxz  - 2px )2 , and (4dyz  - 2py )2 . Bond lengths, dissociation energies, spectroscopic data and electron affinities were calculated via high-level multireference and coupled-cluster methodology using the aug-cc-pV5ZX (-PP)M basis set. Strong bonding results in short bond lengths ranging from 1.602 (TcN- ) to 1.944 (PdBe- ) Å. Adiabatic (diabatic) binding energies reach up to 139 (184) kcal/mol. Electron affinities (EA) were calculated at 1.368 (TcN), 1.242 (RuC), 0.873 (RhB), 0.743 (PdBe) eV. Only for RhB has EA been measured experimentally at 0.961 eV, in good agreement with the value reported here.

The stability and bonding of a series of hetero-diatomic molecules with general formula (cAAC)EM(CO)n , where cAAC = cyclic alkyl(amino) carbene; E = group 14 elements (C, Si, and Ge); M = transition metal (Ni, Fe, and Cr) have been studied by quantum chemical calculations using density functional theory (DFT) and energy decomposition analysis-natural orbital chemical valence (EDA-NOCV). The equilibrium geometries were calculated at the BP86/def2-TZVPP level of theory. The tri-coordinated group 14 complex (1a, 4a, and 7a) in which one of the CO groups is migrated to the central group 14 element from adjacent metal is theoretically found to be more stable when the central atom (E) is carbon. On the other hand, the two-coordinate group 14 element containing metal-complexes (2, 5, 8, 3, 6, and 9) are found to be more stable with their corresponding heavier analogues. The electronic structures of all the molecules have been analyzed by molecular orbital, topological analysis of electron density and natural bond orbital (NBO) analysis at the M06/def2-TZVPP//BP86/def2-TZVPP level of theory. The nature of the cAACE and EM bonds has been studied by EDA-NOCV calculations at BP86-D3(BJ)/TZ2P level of theory. The EDA analysis suggests that the bonding of cAACC(CO) can be best represented by electron sharing σ and π interactions, whereas, C(CO)M(CO)n-1 by dative σ and π interactions. On the other hand, EDA-NOCV calculations suggests both dative σ and π interactions for cAACE and EM(CO)n bonds of the corresponding Si and Ge analogues having stronger σ- and relatively weaker π-bonds. The topological analysis of electron density supports the closed-shell interaction for the Si and Ge complexes and open-shell interaction for the carbon complexes. The calculated proton affinity and hydride affinity values corroborated well with the present bonding description. This class of complexes might act as efficient future catalysts for different organic transformations due to the presence of electron rich group 14 element and metal carbonyl.