Use of Vertex Index in Structure-Activity Analysis and Design of Molecules

The last few decades have witnessed application of graph theory and topological indices derived from molecular graph in structure-activity analysis. Such applications are based on regression and various multivariate analyses. Most of the topological indices are computed for the whole molecule and used as descriptors for explaining properties/activities of chemical compounds. However, some substructural descriptors in the form of topological distance based vertex indices have been found to be useful in identifying activity related substructures and in predicting pharmacological and toxicological activities of bioactive compounds. Another important aspect of drug discovery e. g. designing novel pharmaceutical candidates could also be done from the distance distribution associated with such vertex indices. In this article, we will review the development and applications of this approach both in activity prediction as well as in designing novel compounds.

Wave propagation in single-walled carbon nanotube under longitudinal magnetic field using nonlocal Euler-Bernoulli beam theory

In the present work, the effect of longitudinal magnetic field on wave dispersion characteristics of equivalent continuum structure (ECS) of single-walled carbon nanotubes (SWCNT) embedded in elastic medium is studied. The ECS is modelled as an Euler-Bernoulli beam. The chemical bonds between a SWCNT and the elastic medium are assumed to be formed. The elastic matrix is described by Pasternak foundation model, which accounts for both normal pressure and the transverse shear deformation. The governing equations of motion for the ECS of SWCNT under a longitudinal magnetic field are derived by considering the Lorentz magnetic force obtained from Maxwell's relations within the frame work of nonlocal elasticity theory. The wave propagation analysis is performed using spectral analysis. The results obtained show that the velocity of flexural waves in SWCNTs increases with the increase of longitudinal magnetic field exerted on it in the frequency range: 0-20 THz. The present analysis also shows that the flexural wave dispersion in the ECS of SWCNT obtained by local and nonlocal elasticity theories differ. It is found that the nonlocality reduces the wave velocity irrespective of the presence of the magnetic field and does not influences it in the higher frequency region. Further it is found that the presence of elastic matrix introduces the frequency band gap in flexural wave mode. The band gap in the flexural wave is found to independent of strength of the longitudinal magnetic field. (C) 2011 Elsevier Inc. All rights reserved.

Study of molecular structure and vibrational spectra of poly (beta,-1-malic acid) PMLA] using DFT approach

Poly (beta-L-malic acid) (PMLA) is a biodegradable polymer and it has various important applications in the biomedical field. In the present work the structural and spectral characteristics of PMLA have been studied by methods of infrared. Raman spectroscopy and quantum chemistry. Electrostatic potential surface, optimized geometry, harmonic vibrational frequencies, infrared intensities and activities of Raman scattering were calculated by density functional theory (DFT) using oligomeric approach employing B3LYP with complete relaxation in the potential energy surface using 6-311++G (d, p) basis set. Based on results, we have discussed the correlation between the vibrational modes and the structure of the PMLA. A complete analysis of the experimental infrared and Raman spectra has been reported on the basis of wavenumber of the vibrational bands and potential energy distribution. The calculated HOMO and LUMO energies shows that charge transfer occur within the molecule. The calculated infrared and the Raman spectra of the polymer based on DFT calculations show reasonable agreement with the experimental results. (c) 2012 Elsevier Ltd. All rights reserved.

A Rare Phenoxido/Acetato/Azido Bridged Trinuclear and an Unprecedented Phenoxido/Azido Bridged One-Dimensional Polynuclear Nickel(II) Complexes: Synthesis, Crystal Structure, and Magnetic Properties with Theoretical Investigations on the Exchange Mechanism

The reaction of a tridentate Schiff base ligand HL (2-(3-dimethylaminopropylimino)-methyl]-phenol) with Ni(II) acetate or perchlorate salts in the presence of azide as coligand has led to two new Ni(II) complexes of formulas Ni3L2(OAc)(2)(mu(1,1)-N-3)(2)(H2O)(2)]center dot 2H(2)O (1) and Ni2L2(mu(1,1)-N-3) (mu(1,3)-N-3)](n)(2). Single crystal X-ray structures show that complex 1 is a linear trinuclear Ni(II) compound containing a mu(2)-phenwddo, an end-on (EO) azido and a syn-syn acetato bridge between the terminal and the central Ni(II) ions. Complex 2 can be viewed as a one-dimensional (1D) chain in which the triply bridged (di-mu(2)-phenoxido and EO azido) dimeric Ni-2 units are linked to each other in a zigzag pattern by a single end-to-end (EE) azido bridge. Variable-temperature magnetic susceptibility studies indicate the presence of moderate ferromagnetic exchange coupling in complex 1 with J value of 16.51(6) cm(-1). The magnetic behavior of 2 can be fitted in an alternating ferro- and antiferromagnetic model J(FM) = +34.2(2.8) cm(-1) and J(AF) = -21.6(1.1) cm(-1)] corresponding to the triple bridged dinuclear core and EE azido bridge respectively. Density functional theory (DFT) calculations were performed to corroborate the magnetic results of 1 and 2. The contributions of the different bridges toward magnetic interactions in both compounds have also been calculated.

DFT study of structure and vibrational spectra of ceramide 3: comparison to experimental data

The Fourier transform Raman and infrared (IR) spectra of the Ceramide 3 (CER3) have been recorded in the regions 200-3500 cm(-1) and 680-4000 cm(-1), respectively. We have calculated the equilibrium geometry, harmonic vibrational wavenumbers, electrostatic potential surfaces, absolute Raman scattering activities and IR absorption intensities by the density functional theory with B3LYP functionals having extended basis set 6-311G. This work is undertaken to study the vibrational spectra of CER3 completely and to identify the various normal modes with better wavenumber accuracy. Good consistency is found between the calculated results and experimental data for the IR and Raman spectra.

Activation of Lattice Oxygen of TiO2 by Pd2+ Ion: Correlation of Low-Temperature CO and Hydrocarbon Oxidation with Structure of Ti1-xPdxO2-x (x=0.01-0.03)

Lattice oxygen of TiO2 is activated by the substitution of Pd ion in its lattice. Ti1-xPdxO2-x (x = 0.01-0.03) have been synthesized by solution combustion method crystallizing in anatase TiO2 structure. Pd is in +2 oxidation state and Ti is in +4 oxidation state in the catalyst. Pd is more ionic in TiO2 lattice compared to Pd in PdO. Oxygen storage capacity defined by ``amount of oxygen that is used reversibly to oxidize CO'' is as high as 5100 mu mol/g of Ti0.97Pd0.03O1.97. Oxygen is extracted by CO to CO2 in absence of feed oxygen even at room temperature which is more than 20 times compared to pure TiO2. Rate of CO oxidation is 2.75 mu mol g(-1) s(-1) at 60 degrees C over Ti0.97Pd0.03O1.97 and C2H2 gets oxidized to CO2 and H2O at room temperature. Catalyst is not poisoned on long time operation of the reactor. Such high catalytic activity is due to activated lattice oxygen created by the substitution of Pd ion as seen from first-principles density functional theory (DFT) calculations with 96 atom supercells of Ti32O64, Ti31Pd1O63, Ti30Pd2O62, and Ti29Pd3O61. The compounds crystallize in anatase TiO2 structure with Pd2+ ion in nearly square planar geometry and TiO6 octahedra are distorted by the creation of weakly bound oxygens. Structural analysis of Ti31Pd1O63 which is close to 3% Pd ion substituted TiO2 shows that oxygens associated with both Ti and Pd ions in the lattice show bond valence sum of 1.87, a low value characteristic of weak oxygen in the lattice compared to oxygens with valence 2 and above in the same lattice. Exact positions of activated oxygens have been identified in the lattice from DFT calculations.

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a ``spatial motif'' and several ``fold specific hot spots'' that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.

Effect of chlorine substitution on triplet state structure of thioxanthone: A time-resolved resonance Raman study

Substitution plays an important role in determining the triplet state reactivity. In this paper, we have studied the effect of chlorine substitution on the triplet state structure and the reactivity of thioxanthone (TX). We have employed time-resolved resonance Raman technique to understand the structure of the lowest triplet excited state of 2-chlorothioxanthone (CTX). The experimental findings have been corroborated with the computational results using density functional theory. Akin to the parent compound (TX), coexistence of two lowest triplet states has been observed in case of CTX, which has been substantiated using resonant probe wavelength dependence study. The relative contribution of 3n-pi* to 3 pi-pi* to the equilibrated triplet state has been found to be more for CTX compared to TX suggesting increase in the triplet state reactivity after the substitution. The above observation has been further supported by the flash photolysis experiments. Copyright (C) 2013 John Wiley & Sons, Ltd.

Exploring the color of transition metal ions in irregular coordination geometries: new colored inorganic oxides based on the spiroffite structure, Zn2-xMxTe3O8 (M = Co, Ni, Cu)

We describe the synthesis, crystal structures, and optical absorption spectra of transition metal substituted spiroffite derivatives, Zn2-xMxTe3O8 (M-II = Co, Ni, Cu; 0 < x <= 1.0). The oxides are readily synthesized by solid state reaction of stoichiometric mixtures of the constituent binaries at 620 degrees C. Reitveld refinement of the crystal structures from powder X-ray diffraction (XRD) data shows that the Zn/MO6 octahedra are strongly distorted, as in the parent Zn2Te3O8 structure, consisting of five relatively short Zn/M-II-O bonds (1.898-2.236 angstrom) and one longer Zn/M-II-O bond (2.356-2.519 angstrom). We have interpreted the unique colors and the optical absorption/diffuse reflectance spectra of Zn2-xMxTe3O8 in the visible, in terms of the observed/irregular coordination geometry of the Zn/M-II-O chromophores. We could not however prepare the fully substituted M2Te3O8 (M-II = Co, Ni, Cu) by the direct solid state reaction method. Density Functional Theory (DFT) modeling of the electronic structure of both the parent and the transition metal substituted derivatives provides new insights into the bonding and the role of transition metals toward the origin of color in these materials. We believe that transition metal substituted spiroffites Zn2-xMxTe3O8 reported here suggest new directions for the development of colored inorganic materials/pigments featuring irregular/distorted oxygen coordination polyhedra around transition metal ions.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Structure constants of beta deformed super Yang-Mills

We study the structure constants of the N = 1 beta deformed theory perturbatively and at strong coupling. We show that the planar one loop corrections to the structure constants of single trace gauge invariant operators in the scalar sector is determined by the anomalous dimension Hamiltonian. This result implies that 3 point functions of the chiral primaries of the theory do not receive corrections at one loop. We then study the structure constants at strong coupling using the Lunin-Maldacena geometry. We explicitly construct the supergravity mode dual to the chiral primary with three equal U(1) R-charges in the Lunin-Maldacena geometry. We show that the 3 point function of this supergravity mode with semi-classical states representing two other similar chiral primary states but with large U(1) charges to be independent of the beta deformation and identical to that found in the AdS(5) x S-5 geometry. This together with the one-loop result indicate that these structure constants are protected by a non-renormalization theorem. We also show that three point function of U(1) R-currents with classical massive strings is proportional to the R-charge carried by the string solution. This is in accordance with the prediction of the R-symmetry Ward identity.

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Spatial variations of effective elastic thickness over the Ninetyeast Ridge and implications for its structure and tectonic evolution

We present new data on the strength of oceanic lithosphere along the Ninetyeast Ridge (NER) from two independent methods: spectral analysis (Bouguer coherence) using the fan wavelet transform technique, and spatial analysis (flexure inversion) with the convolution method. The two methods provide effective elastic thickness (T-e) patterns that broadly complement each other, and correlate well with known surface structures and regional-scale features. Furthermore, our study presents a new high resolution database on the Moho configuration, which obeys flexural isostasy, and exhibit regional correlations with the T-e variations. A continuous ridge structure with a much lower T-e value than that of normal oceanic lithosphere provides strong support for the hotspot theory. The derived T-e values vary over the northern (higher T-e similar to 10-20 km), central (anomalously low T-e similar to 0-5 km), and southern (low T-e similar to 5 km) segments of the NER. The lack of correlation of the T-e value with the progressive aging of the lithosphere implies differences in thermo-mechanical setting of the crust and underlying mantle in different parts of the NER, again indicating diversity in their evolution. The anomalously low T-e and deeper Moho (similar to 22 km) estimates of the central NER (between 0.5 degrees N and 17 degrees S) are attributed to the interaction of a hotspot with the Wharton spreading ridge that caused significant thermal rejuvenation and hence weakening of the lithosphere. The higher mechanical strength values in the northern NER (north of 0.5 degrees N) may support the idea of off-ridge emplacement and a relatively large plate motion at the time of volcanism. The low T-e and deeper Moho (similar to 22 km) estimates in the southern part (south of 17 degrees S) suggest that the lithosphere was weak and therefore younger at the time of volcanism, and this supports the idea that the southern NER was emplaced on the edge of the Indian plate. (C) 2013 Elsevier B.V. All rights reserved.

Electronic structure of Nd1-xYxMnO3 from Mn K edge absorption spectroscopy and DFT methods

The electronic structure of Nd1-xYxMnO3 (x-0-0.5) is studied using x-ray absorption near-edge structure (XANES) spectroscopy at the Mn K-edge along with the DFT-based LSDA+U and real space cluster calculations. The main edge of the spectra does not show any variation with doping. The pre-edge shows two distinct features which appear well-separated with doping. The intensity of the pre-edge decreases with doping. The theoretical XANES were calculated using real space multiple scattering methods which reproduces the entire experimental spectra at the main edge as well as the pre-edge. Density functional theory calculations are used to obtain the Mn 4p, Mn 3d and O 2p density of states. For x=0, the site-projected density of states at 1.7 eV above Fermi energy shows a singular peak of unoccupied e(g) (spin-up) states which is hybridized Mn 4p and O 2p states. For x=0.5, this feature develops at a higher energy and is highly delocalized and overlaps with the 3d spin-down states which changes the pre-edge intensity. The Mn 4p DOS for both compositions, show considerable difference between the individual p(x), p(y) and p(z)), states. For x=0.5, there is a considerable change in the 4p orbital polarization suggesting changes in the Jahn-Teller effect with doping. (C) 2013 Elsevier Ltd. All rights reserved.

Enhanced Gas Adsorption on Graphitic Substrates via Defects and Local Curvature: A Density Functional Theory Study

Using van-der-Waals-corrected density functional theory calculations, we explore the possibility of engineering the local structure and morphology of high-surface-area graphene-derived materials to improve the uptake of methane and carbon dioxide for gas storage and sensing. We test the sensitivity of the gas adsorption energy to the introduction of native point defects, curvature, and the application of strain. The binding energy at topological point defect sites is inversely correlated with the number of missing carbon atoms, causing Stone-Wales defects to show the largest enhancement with respect to pristine graphene (similar to 20%). Improvements of similar magnitude are observed at concavely curved surfaces in buckled graphene sheets under compressive strain, whereas tensile strain tends to weaken gas binding. Trends for CO2 and CH4 are, similar, although CO2 binding is generally stronger by similar to 4 to 5 kJ mol(-1). However, the differential between the adsorption of CO2 and CH4 is much higher on folded graphene sheets and at concave curvatures; this could possibly be leveraged for CH4/CO2 flow separation and gasselective sensors.