An efficient protocol for obtaining accurate hydration free energies using quantum chemistry and reweighting from molecular dynamics simulations

The non-Boltzmann Bennett (NBB) free energy estimator method is applied to 21 molecules from the blind subset of the SAMPL4 challenge. When NBB is applied with the SMD implicit solvent model, and the OLYP/DZP level of quantum chemistry, highly accurate hydration free energy calculations are obtained with respect to experiment (RMSD = 0.89 kcal·mol⁻¹). Other quantum chemical methods are also tested, and the effects of solvent model, density functional, basis set are explored in this benchmarking study, providing a framework for improvements in calculating hydration free energies. We provide a practical guide for using the best QM-NBB protocols that are consistently more accurate than either pure QM or pure MM alone. In situations where high accuracy hydration free energy predictions are needed, the QM-NBB method with SMD implicit solvent should be the first choice of quantum chemists.     

Correlation between Li Plating Behavior and Surface Characteristics of Carbon Matrix toward Stable Li Metal Anodes

Carbonaceous materials are widely employed to host Li for stable and safe Li metal batteries while relatively little effort is devoted to tailoring the surface properties of carbon to facilitate uniform Li plating. Herein, the correlation between Li plating behavior and the surface characteristics of electrospun porous carbon nanofibers (PCNFs) is systemically elucidated through experiments and theoretical calculations. It is revealed that the neat carbon surface suffers from severe lattice mismatch with Li metal, hindering uniform Li plating. In contrast, open pores created on the PCNF surface serve as active sites for controlled initial nucleation of Li. The introduction of oxygenated functional groups further facilitates the nucleation of Li on PCNFs through the largely reduced nucleation energy barrier. The Li film uniformly deposited on PCNFs enables efficient use of the whole carbon surface, giving rise to enhanced cyclic stability of the electrode. When used as an anode in lithium–sulfur batteries, the modified electrode delivers an excellent energy density of 385 Wh kg^?1 after 100 cycles. The fundamental correlation established in this study is universal to all types of carbonaceous materials and sheds new light on the rational design of high‐performance Li metal anodes by controlling the initial Li nucleation.     

Comparative analysis of the methods used for finding surface energy to investigate protein interaction behavior on chromatographic supports

Hydrophobic interaction chromatography, an important and effective purification strategy, is generally used for the purification of variety of biomolecules. A basic understanding of the protein interaction behavior is required to effectively separate these biomolecules. A colloidal type extended Derjaguin, Landau, Verwey, and Overbeek calculations were utilized to study the interactions behavior of model proteins to commercially available hydrophobic chromatographic materials that is, Toyopearl Phenyl 650C and Toyopearl Butyl 650C. Physicochemical properties of selected model proteins were achieved by contact angle and zeta potential measurements. The contact angle of chromatographic materials used was achieved through sessile drop method on disrupted beads and capillary penetration method (CPM) on intact beads. The surface properties were further used to calculate the interactions of the proteins to chromatographic supports. The calculated secondary energy minimum of the proteins with the chromatographic materials (from the contact angle values determined through both methods can be correlated with the retention volumes from the real chromatography. The secondary energy minimum values are higher for each protein to the chromatographic materials calculated from the inputs derived through sessile drop method compared to CPM. For instance, immunoglobulin G has secondary energy minimum value of 0.17 kT compared to 0.11 kT, obtained through sessile drop method and CPM, respectively. Average relative values of the energy minimum calculated for all proteins are as 1.51 kT and 1.29 kT for Toyopearl Butyl 650C and Toyopearl Phenyl 650C, respectively, as a conversion factor for estimation of secondary energy minimum for both methods.     

B3LYP/6-311++G* * study of structure and spin-spin coupling constant in heparin disaccharide

Structures of heparin disaccharide have been analyzed by DFT using the B3LYP/6-311++G( * *) method. The optimized geometries of two forms of this disaccharide, differing in the conformation ((1)C(4) and (2)S(0)) of the IdoA2S residue, confirmed considerable influences of the sulfate and the carboxylate groups upon the pyranose ring geometries. The computed energies showed that disaccharide having the (1)C(4) form of the IdoA2S residue is more stable than that with the (2)S(0) form. Interatomic distances, bond and torsion angles showed that interconversion of the IdoA2S residue results in geometry changes in the GlcN,6S residue as well. Three-bond proton-proton and proton-carbon spin-spin coupling constants computed for both forms agree with the experimental data and indicate that only two chair forms contribute to the conformational equilibrium in disaccharide. Influences of the charged groups upon the magnitudes of spin-spin coupling constants are also discussed.     

Two new antibacterial sesquiterpenoids from Centipeda minima

Two new guaiane-type sesquiterpene lactones, compounds 1 and 2, along with three known guaianolide- or pseudoguaianolides, were isolated from Centipeda minima (whole plant). Their structures were identified by spectroscopic and mass-spectrometric analyses. The configuration at C5 of the guaiane framework of 1 was rationalized by quantum-mechanical calculations (Table 2). All compounds were found to be active against eight different microbial pathogens (Table 3), with MIC values in the range of 6.25-100 microg/ml.     

Identification and characterization of carbapenem binding sites within the RND-transporter AcrB

Understanding the molecular determinants for recognition, binding and transport of antibiotics by multidrug efflux systems is important for basic research and useful for the design of more effective antimicrobial compounds. Imipenem and meropenem are two carbapenems whose antibacterial activity is known to be poorly and strongly affected by MexAB-OprM, the major efflux pump transporter in Pseudomonas aeruginosa. However, not much is known regarding recognition and transport of these compounds by AcrAB-TolC, which is the MexAB-OprM homologue in Escherichia coli and by definition the paradigm model for structural studies on efflux pumps. Prompted by this motivation, we unveiled the molecular details of the interaction of imipenem and meropenem with the transporter AcrB by combining computer simulations with biophysical experiments. Regarding the interaction with the two main substrate binding regions of AcrB, the so-called access and deep binding pockets, molecular dynamics simulations revealed imipenem to be more mobile than meropenem in the former, while comparable mobilities were observed in the latter. This result is in line with isothermal titration calorimetry, differential scanning experiments, and binding free energy calculations, indicating a higher affinity for meropenem than imipenem at the deep binding pocket, while both sharing similar affinities at the access pocket. Our findings rationalize how different physico-chemical properties of compounds reflect on their interactions with AcrB. As such, they constitute precious information to be exploited for the rational design of antibiotics able to evade efflux pumps.     

Exploring the antimalarial potential of the methoxy-thiazinoquinone scaffold: Identification of a new lead candidate

A small library of antiplasmodial methoxy-thiazinoquinones, rationally designed on the model of the previously identified hit 1, has been prepared by a simple and inexpensive procedure. The synthetic derivatives have been subjected to in vitro pharmacological screening, including antiplasmodial and toxicity assays. These studies afforded a new lead candidate, compound 9, endowed with higher antiplasmodial potency compared to 1, a good selectivity index when tested against a panel of mammalian cells, no toxicity against RBCs, a synergistic antiplasmodial action in combination with dihydroartemisinin, and a promising inhibitory activity on stage V gametocyte growth. Computational studies provided useful insights into the structural requirements needed for the antiplasmodial activity of thiazinoquinone compounds and on their putative mechanism of action.     

First-principles investigations on structural stability,elastic and electronic properties of Co7M6 (M= W,Mo, Nb) µ phases

The structural, elastic and electronic properties of Co₇M₆ (M?=?W, Mo, Nb) μ phases were investigated by first-principles calculations based on the density functional theory (DFT). The calculated cohesive energy indicates that Co₇M₆ (M?=?W, Mo, Nb) μ phases are thermodynamically stable. Besides, Co₇W₆ owns a higher structural stability than that of Co₇Mo₆ and Co₇Nb₆. The obtained elastic constant demonstrates that Co₇M₆ (M?=?W, Mo, Nb) are mechanically stable. With Voigt-Reuss-Hill (VRH) approximation, the elastic bulk modulus (B), shear modulus (G), Young's modulus (E) and Poisson's ratio (ν) were derived. The ductility and plasticity as well as the elastic anisotropy of the three phases were discussed in details. Finally, the density of states and charge density difference were also analysed to reveal the underlying mechanism of structural stability and the elastic properties.     

Molecular modelling and vibrational investigations of ammonium-based ionic liquid (CLTOAB)

CLTOAB is an ammonium-based ionic liquid composed of ε-Caprolactam (CL) C₆H₁₁NO and tetraoctylammonium bromide (TOAB) (C₃₂H₆₈BrN). In this study, experimental IR and Raman spectra of CLTOAB ionic liquid together with the computational results of the compound have been reported. The optimized geometry, vibrational frequencies, IR intensities and Raman activities of the CLTOAB were calculated using the wb97xd and B3LYP density functional methods combined with the 6-31G(d,p) basis set using Gaussian 03 program. The complete assignment of the bands was performed based on the potential energy distributions (PED%). The HOMO and LUMO analysis is used to determine the charge transfer within the molecule. The Gauge-including atomic orbital ¹H-NMR and ¹³C-NMR chemical shifts calculations were carried out and compared with the experimental data. Furthermore to evaluate interaction between CLTOAB and DNA, molecular docking study was carried out.

Communicated by Ramaswamy H. Sarma     

Ab initio vibrational calculations on ara-T molecule: application to analysis of IR and Raman spectra

The FTIR and FT-Raman spectra are reported for the arabinonucleoside ara-T (1-beta-D-arabinofuranosylthymine), which shows antiviral activity. The accurate knowledge of the vibrational modes is a prerequisite for the elucidation of drug-nucleotide and drug-enzyme interactions. The FTIR and FT-Raman spectra of ara-T were recorded from 4000 to 30 cm(-1). A tetradeuterated derivative (deuteration at N3, and hydroxyl groups O'2, O'3, and O'5) was synthesized and the observed isotopic shifts in its spectra were used for the vibrational analysis of ara-T. The theoretical frequencies and the potential energy distribution (PED) of the vibrational modes of ara-T were calculated using the ab initio Hartree-Fock/3-21G method. An assignment of the vibrational spectra of ara-T is proposed considering the scaled PED and the observed band shifts under deuteration. The scaled ab initio frequencies were in reasonable agreement with the experimental data.