MolBridge: a program for identifying nonbonded interactions in small molecules and biomolecular structures

Identification and analysis of nonbonded interactions within a molecule and with the surrounding molecules are an essential part of structural studies, given the importance of these interactions in defining the structure and function of any supramolecular entity. MolBridge is an easy to use algorithm based purely on geometric criteria that can identify all possible nonbonded interactions, such as hydrogen bond, halogen bond, cation-pi, pi-pi and van der Waals, in small molecules as well as biomolecules. The user can either upload three-dimensional coordinate files or enter the molecular ID corresponding to the relevant database. The program is available in a standalone form and as an interactive web server with Jmol and JME incorporated into it. The program is freely downloadable and the web server version is also available at http://nucleix.mbu.iisc.ernet.in/molbridge/index.php.

Characterization of Interactions Involving Bromine in 2,2-Dibromo2,3-dihydroinden-1-one via Experimental Charge Density Analysis

Experimental and theoretical charge density analyses on 2,2-dibromo-2,3-dihydroinden-1-one have been carried out to quantify the topological features of a short CBr....O halogen bond with nearly linear geometry (2.922 angstrom, angle CBr....O = 172.7 degrees) and to assess the strength of the interactions using the topological features of the electron density. The electrostatic potential map indicates the presence of the s-hole on bromine, while the interaction energy is comparable to that of a moderate OH....O hydrogen bond. In addition, the energetic contribution of CH.....Br interaction is demonstrated to be on par with that of the CBr....O halogen bond in stabilizing the crystal structure.

Computational Study of the Formation of Short Centrosymmetric N-center dot center dot center dot S Supramolecular Synthon and Related Weak Interactions in Crystalline 1,2,4-Triazoles

A comprehensive analysis of the crystal packing and the energetic features of a series of four biologically active molecules belonging to the family of substituted 4-(benzylideneamino)-3-(4-fluoro-3-phenoxyphenyl)-1H-1,2,4-triazole-5-(4 H)-thione derivatives have been performed based on the molecular conformation and the supramolecular packing. This involves the formation of a short centrosymmetric R-2(2)(8) NH...S supramolecular synthon in the solid state, including the presence of CH...S, CH...O, CH...N, CH...F, CH...Cl, CF...FC, CCl...ClC, and CH...pi intermolecular interactions along with pp stacking to evaluate the role of noncovalent interactions in the crystal. The presence of such synthons has a substantial contribution toward the interaction energy (-18 to -20 kcal/mol) as obtained from the PIXEL calculation, wherein the Coulombic and polarization contribution are more significant than the dispersion contribution. The geometrical characteristics of such synthons favor short distance, and the population of related molecules having these geometries is rare as has been obtained from the Cambridge Structural Database (CSD). Furthermore, their interaction energies have been compared with those present in our molecules in the solid state. The topological characteristics of the NH...S supramolecular synthon, in addition to related weak interactions, CH...N, CH...Cl, CF...FC, and CCl...ClC, have been estimated using the quantum theory of atoms in molecules (QTAIM). In addition, an analysis of the Hirshfeld surface and associated fingerprint plots of these four molecules also have provided a platform for the evaluation of the contribution of different atom...atom contacts, which contribute toward the packing of the molecules in solids.

Quick re-introduction of selective scalar interactions in a pure-shift NMR spectrum

A new 1D NMR experiment cited as `Quick G-SERF', which re-introduces selective proton-proton scalar interactions in a pure shift spectrum during real time data acquisition, is reported. The method provides information on multiple proton-proton couplings from a single experiment, analogous to the 2D G-SERF technique, while significantly shortening the experimental time by 1-2 orders of magnitude due to reduced dimension and enhanced sensitivity.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

Heterogeneity in Dye-TiO2 Interactions Dictate Charge Transfer Efficiencies for Diketopyrrolopyrrole-Based Polymer Sensitized Solar Cells

Power conversion efficiency of a solar cell is a complex parameter which usually hides the molecular details of the charge generation process. For rationally tailoring the overall device efficiency of the dye-sensitized solar cell, detailed molecular understanding of photoinduced reactions at the dye-TiO2 interface has to be achieved. Recently, near-IR absorbing diketopyrrolopyrrole-based (DPP) low bandgap polymeric dyes with enhanced photostabilities have been used for TiO2 sensitization with moderate efficiencies. To improve the reported device performances, a critical analysis of the polymerTiO(2) interaction and electron transfer dynamics is imperative. Employing a combination of time-resolved optical measurements complemented by low temperature EPR and steady-state Raman spectroscopy on polymerTiO(2) conjugates, we provide direct evidence for photoinduced electron injection from the TDPP-BBT polymer singlet state into TiO2 through the C-O group of the DPP-core. A detailed excited state description of the electron transfer process in films reveals instrument response function (IRF) limited (<110 fs) charge injection from a minor polymer fraction followed by a picosecond recombination. The major fraction of photoexcited polymers, however, does not show injection indicating pronounced ground state heterogeneity induced due to nonspecific polymerTiO(2) interactions. Our work therefore underscores the importance of gathering molecular-level insight into the competitive pathways of ultrafast charge generation along with probing the chemical heterogeneity at the nanoscale within the polymerTiO2 films for optimizing photovoltaic device efficiencies.

Cocrystals of Hydrochlorothiazide: Solubility and Diffusion/Permeability Enhancements through Drug-Coformer Interactions

Hydrochlorothiazide (HCT) is a diuretic and a BCS class IV drug with low solubility and low permeability, exhibiting poor oral absorption. The present study attempts to improve the physicochemical properties of the drug using a crystal engineering approach with cocrystals. Such multicomponent crystals of HCT with nicotinic acid (NIC), nicotinamide (NCT), 4-aminobenzoic acid (PABA), succinamide (SAM), and resorcinol (RES) were prepared using liquid-assisted grinding, and their solubilities in pH 7.4 buffer were evaluated. Diffusion and membrane permeability were studied using a Franz diffusion cell. Except for the SAM and NIC cocrystals, all other binary systems exhibited improved solubility. All of the cocrystals showed improved diffusion/membrane permeability compared to that of HCT with the exception of the SAM cocrystal. When the solubility was high, as in the case of PABA, NCT, and RES cocrystals, the flux/permeability dropped slightly. This is in agreement with the expected interplay between solubility and permeability. Improved solubility/permeability is attributed to new drug-coformer interactions. Cocrystals of SAM, however, showed poor solubility and flux This cocrystal contains a primary sulfonamide dimer synthon similar to that of HCT polymorphs, which may be a reason for its unusual behavior. Hirshfeld surface analysis was carried out in all cases to determine whether a correlation exists between cocrystal permeability and drug-coformer interactions.

Polyfunctional Lewis Acids: Intriguing Solid-State Structure and Selective Detection and Discrimination of Nitroaromatic Explosives

Synthesis and crystal structures of three porphyrin-based polyfunctional Lewis acids 1-3 are reported. Intermolecular HgClHgCl (linear and -type) interactions in the solid state of the peripherally ArHgCl-decorated compound 3 lead to a fascinating 3D supramolecular architecture. Compound3 shows a selective fluorescence quenching response to picric acid and discriminates other nitroaromatic-based explosives. For the first time, an electron-deficient polyfunctional Lewis acid is shown to be useful for the selective detection and discrimination of nitroaromatic explosives. The Stern-Volmer quenching constant and detection limits of compound3 for picric acid are the best among the reported small-molecular receptors for nitroaromatic explosives. The electronic structure, Lewis acidity, and selective sensing characteristics of 3 are well corroborated by DFT calculations.

An ultrastable conjugate of silver nanoparticles and protein formed through weak interactions

In recent years, silver nanoparticles (AgNPs) have attracted significant attention owing to their unique physicochemical, optical, conductive and antimicrobial properties. One of the properties of AgNPs which is crucial for all applications is their stability. In the present study we unravel a mechanism through which silver nanoparticles are rendered ultrastable in an aqueous solution in complex with the protein ubiquitin (Ubq). This involves a dynamic and reversible association and dissociation of ubiquitin from the surface of AgNP. The exchange occurs at a rate much greater than 25 s(-1) implying a residence time of <40 ms for the protein. The AgNP-Ubq complex remains stable for months due to steric stabilization over a wide pH range compared to unconjugated AgNPs. NMR studies reveal that the protein molecules bind reversibly to AgNP with an approximate dissociation constant of 55 mu M and undergo fast exchange. At pH > 4 the positively charged surface of the protein comes in contact with the citrate capped AgNP surface. Further, NMR relaxation-based experiments suggest that in addition to the dynamic exchange, a conformational rearrangement of the protein takes place upon binding to AgNP. The ultrastability of the AgNP-Ubq complex was found to be useful for its anti-microbial activity, which allowed the recycling of this complex multiple times without the loss of stability. Altogether, the study provides new insights into the mechanism of protein-silver nanoparticle interactions and opens up new avenues for its application in a wide range of systems.

Energy-savvy solid-state and sonochemical synthesis of lithium sodium titanate as an anode active material for Li-ion batteries

Lithium sodium titanate insertion-type anode has been synthesized by classical solid-state (dry) and an alternate solution-assisted (wet) sonochemical synthesis routes. Successful synthesis of the target compound has been realized using simple Na- and Li-hydroxide salts along with titania. In contrast to the previous reports, these energy-savvy synthesis routes can yield the final product by calcination at 650 -750 degrees C for limited duration of 1-10 h. Owing to the restricted calcination duration (dry route for 1-2 h and wet route for 1-5 h), they yield homogeneous nanoscale lithium sodium titanate particles. Sono-chemical synthesis reduces the lithium sodium titanate particle size down to 80-100 nm vis-a-vis solid-state method delivering larger (200-500 nm) particles. Independent of the synthetic methods, the end products deliver reversible electrochemical performance with reversible capacity exceeding 80 mAh.g(-1) acting as a 1.3 V anode for Li-ion batteries. (C) 2015 Elsevier B.V. All rights reserved.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

Origin of the Spin-Orbital Liquid State in a Nearly J=0 Iridate Ba3ZnIr2O9

We show using detailed magnetic and thermodynamic studies and theoretical calculations that the ground state of Ba3ZnIr2O9 is a realization of a novel spin-orbital liquid state. Our results reveal that Ba3ZnIr2O9 with Ir5+ (5d(4)) ions and strong spin-orbit coupling (SOC) arrives very close to the elusive J = 0 state but each Ir ion still possesses a weak moment. Ab initio density functional calculations indicate that this moment is developed due to superexchange, mediated by a strong intradimer hopping mechanism. While the Ir spins within the structural Ir2O9 dimer are expected to form a spin-orbit singlet state (SOS) with no resultant moment, substantial frustration arising from interdimer exchange interactions induce quantum fluctuations in these possible SOS states favoring a spin-orbital liquid phase down to at least 100 mK.

Triplet excited electronic state switching induced by hydrogen bonding: A transient absorption spectroscopy and time-dependent DFT study

The solvent plays a decisive role in the photochemistry and photophysics of aromatic ketones. Xanthone (XT) is one such aromatic ketone and its triplet-triplet (T-T) absorption spectra show intriguing solvatochromic behavior. Also, the reactivity of XT towards H-atom abstraction shows an unprecedented decrease in protic solvents relative to aprotic solvents. Therefore, a comprehensive solvatochromic analysis of the triplet-triplet absorption spectra of XT was carried out in conjunction with time dependent density functional theory using the ad hoc explicit solvent model approach. A detailed solvatochromic analysis of the T-T absorption bands of XT suggests that the hydrogen bonding interactions are different in the corresponding triplet excited states. Furthermore, the contributions of non-specific and hydrogen bonding interactions towards differential solvation of the triplet states in protic solvents were found to be of equal magnitude. The frontier molecular orbital and electron density difference analysis of the T-1 and T-2 states of XT indicates that the charge redistribution in these states leads to intermolecular hydrogen bond strengthening and weakening, respectively, relative to the S-0 state. This is further supported by the vertical excitation energy calculations of the XT-methanol supra-molecular complex. The intermolecular hydrogen bonding potential energy curves obtained for this complex in the S-0, T-1, and T-2 states support the model. In summary, we propose that the different hydrogen bonding mechanisms exhibited by the two lowest triplet excited states of XT result in a decreasing role of the n pi* triplet state, and are thus responsible for its reduced reactivity towards H-atom abstraction in protic solvents. (C) 2016 AIP Publishing LLC.