Understanding the quantum mechanical properties of hydrogen bonds in solvated biomolecules from cluster calculations

Many-body analysis and spectroscopic characterization of diazene oligomers: A theoretical study

The present study reports the many-body analysis and spectroscopic characterization of linear and cyclic diazene oligomers in gas and water solvent states. The oligomers of diazene from monomer to pentamer have been considered for the study. The spectroscopic studies such as geometrical parameters, infrared spectra, electronic absorption spectra, and natural transition orbitals (NTOs) were reported. Many-body analysis techniques have been implemented to study the interactions among the diazene oligomers. These calculations have been performed using exchange and correlation functional (B3LYP) and 6-311++G (d,p) basis set. The geometrical parameters and infrared modes of monomer diazene in the gas state are well-matched with the available experimental determinations at this level of theory. A significant change in vibrational modes of linear and cyclic diazene oligomers has been observed in the gas phase-to-water solvent state. The time-dependent density functional theory (TD-DFT) has been used to calculate the electronic absorption spectra of diazene oligomers. The Wavelength of electronic transitions, oscillator strength, and HOMO to LUMO gap has been reported. Many-body analysis shows that two-, three-, four-, and five-body energies have a remarkable contribution to the binding energy in addition to relaxation energies. All these calculations have been performed using Gaussian 16 program package.

Direct and Reliable Method for Estimating the Hydrogen Bond Energies and Cooperativity in Water Clusters, Wn, n = 3 to 8

No direct method for estimating the individual O-H···O hydrogen bond (H-bond) energies in water clusters (W_n) exists in the literature. In this work, we propose such a direct method based on the molecular tailoring approach, which also enables the estimation of the cooperativity contributions. The calculated H-bond energies at MP2(full)/aug-cc-pVTZ and CCSD(T)/aug-cc-pVDZ levels for W_n, n = 3 to 8, agree well with one another and fall between 0.3 and 11.6 kcal mol^-1 with the cooperativity contributions in the range of -1.2 and 7.0 kcal mol^-1. For gauging the accuracy of our H-bond energies for a cluster, the H-bond energy sum is added to the sum of monomer energies, and the results are compared with the respective total energy. These two values agree with each other to within 8.3 mH (∼5 kcal mol^-1), testifying the accuracy of our estimated H-bond energies. Further, these H-bond strengths show a good correlation with the respective O-H stretching frequencies and the molecular electron density values at the (3, -1) O-H···O H-bond critical point.

A DFT study of structure and stability of pleated and rippled cross-β sheets with hydrophobic sidechains

The rippled cross-β sheet, a topography, in which mirror-image peptides are arranged with alternating chirality into a periodic two-dimensional network, is burgeoning as a new design principle for materials and biomedical applications. Experiments by the Schneider, Nilsson, and Raskatov labs have independently shown diverse racemic mixtures of aggregation-prone peptide of different sizes to favor the rippled over the pleated topography. Yet, systematic ab initio studies are lacking, and the field is yet to develop rules that would enable the design of new rippled cross-β frameworks from first principles. Here, DFT calculations were performed on a set of model systems, designed to begin understanding the impact that bulky, hydrophobic sidechains have upon the formation of pleated and rippled cross-β frameworks. It is hoped that this study will help stimulate the development of a predictive, general framework to enable rational design of rippled cross-β sheets in the future.

Molecular Dynamics Simulation Studies on Structure, Dynamics, and Thermodynamics of Uranyl Nitrate Solution at Various Acid Concentrations

The structural and dynamical characteristics of uranyl ions in an aqueous acidic environment are of immense importance in the field of nuclear fuel reprocessing. In view of that, the structural and dynamical behavior of the uranyl ion in water has been investigated by performing molecular dynamics (MD) simulations using different force fields. All the force fields have depicted similar structural and dynamical properties except the free energy of hydration where the Guilbaud-Wipff (GW) model performs well over the others. The calculated density using MD simulations is found to be in excellent agreement with the measured experimental density, which ensures the accuracy of the adopted GW force field. The calculated surface tension and shear viscosity are seen to be increased with uranyl nitrate concentrations. At a higher concentration of about 4.0 mol/L, the supersaturation effect has been captured by an inflection in the plot of surface tension and shear viscosity against concentration because of the solution heterogeneity, which was correlated by an inflection in the scattering intensity observed by performing the dynamic light scattering experiment. The binding mode of nitrate ions with the uranyl ion is found to be concentration-dependent, and at higher concentration, it is predominantly monodentate.