Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

Molecular Characterization of Sclerospora graminicola, the Incitant of Pearl Millet Downy Mildew Using ISSR Markers

The DNA polymorphism among 22 isolates of Sclerospora graminicola, the causal agent of downy mildew disease of pearl millet was assessed using 20 inter simple sequence repeats (ISSR) primers. The objective of the study was to examine the effectiveness of using ISSR markers for unravelling the extent and pattern of genetic diversity in 22 S. graminicola isolates collected from different host cultivars in different states of India. The 19 functional ISSR primers generated 410 polymorphic bands and revealed 89% polymorphism and were able to distinguish all the 22 isolates. Polymorphic bands used to construct an unweighted pair group method of averages (UPGMA) dendrogram based on Jaccard's co-efficient of similarity and principal coordinate analysis resulted in the formation of four major clusters of 22 isolates. The standardized Nei genetic distance among the 22 isolates ranged from 0.0050 to 0.0206. The UPGMA clustering using the standardized genetic distance matrix resulted in the identification of four clusters of the 22 isolates with bootstrap values ranging from 15 to 100. The 3D-scale data supported the UPGMA results, which resulted into four clusters amounting to 70% variation among each other. However, comparing the two methods show that sub clustering by dendrogram and multi dimensional scaling plot is slightly different. All the S. graminicola isolates had distinct ISSR genotypes and cluster analysis origin. The results of ISSR fingerprints revealed significant level of genetic diversity among the isolates and that ISSR markers could be a powerful tool for fingerprinting and diversity analysis in fungal pathogens.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Specificity for the Exchange of Phospholipids Through Polymyxin B Mediated Intermembrane Molecular Contacts

Structural specificity for the direct vesicle−vesicle exchange of phospholipids through stable molecular contacts formed by the antibiotic polymyxin B (PxB) is characterized by kinetic and spectroscopic methods. As shown elsewhere [Cajal, Y., Rogers, J., Berg, O. G., & Jain, M. K. (1996) Biochemistry 35, 299−308], intermembrane molecular contacts between anionic vesicles are formed by a small number of PxB molecules, which suggests that a stoichiometric complex may be responsible for the exchange of phospholipids. Larger clusters containing several vesicles are formed where each vesicle can make multiple contacts if sterically allowed. In this paper we show that the overall process can be dissected into three functional steps: binding of PxB to vesicles, formation of stable vesicle−vesicle contacts, and exchange of phospholipids. Polycationic PxB binds to anionic vesicles. Formation of molecular contacts and exchange of monoanionic phospholipids through PxB contacts does not depend on the chain length of the phospholipid. Only monoanionic phospholipids (with methanol, serine, glycol, butanol, or phosphatidylglycerol as the second phosphodiester substituent in the head group) exchange through these contacts, whereas dianionic phosphatidic acid does not. Selectivity for the exchange was also determined with covesicles of phosphatidylmethanol and other phospholipids. PxB does not bind to vesicles of zwitterionic phosphatidylcholine, and its exchange in covesicles is not mediated by PxB. Vesicles of dianionic phospholipids, like phosphatidic acid, bind PxB; however, this phospholipid does not exchange. The structural features of the contacts are characterized by the spectroscopic and chemical properties of PxB at the interface. PxB in intermembrane contacts is readily accessible from the aqueous phase to quenchers and reagents that modify amino groups. Results show that PxB at the interface can exist in two forms depending on the lipid/PxB ratio. Additional studies show that stable PxB-mediated vesicle−vesicle contacts may be structurally and functionally distinct from “stalks”, the putative transient intermediate for membrane fusion. The phenomenon of selective exchange of phospholipids through peptide-mediated contacts could serve as a prototype for intermembrane targeting and sorting of phospholipids during their biosynthesis and trafficking in different compartments of a cell. The protocols and results described here also extend the syllogistic foundations of interfacial equilibria and catalysis.

Electronic structure of square planar CuO46− clusters

The results of spin-polarized MSXagr calculations show that the ground state of the CuO 4 6– cluster is essentially non-magnetic in spite of odd number of electrons in the system for short Cu–O distances (1.90 Å) as found in the highT c superconductors. This arises due to the fact that the unpaired electron resides in a molecular orbital with primarily oxygen 3s character. The stability of this molecular orbital is found to be sensitive to the cluster geometry and thus, increase in Cu–O distance (as well as other changes affecting oxygen-oxygen distance) tend to favour a magnetic state. From these calculations we have also estimated the Coulomb correlation strength within the Cu 3d to be about 5.3 eV.

A study of chemically prepared small metal clusters in the nonmetallic regime

Small gold clusters [mean diameter (d)[less, similar] 1.4 nm], unlike larger clusters, show a higher Au(4f) binding energy relative to the bulk value and the presence of a conductance gap in tunnelling measurements, just as the molecular cluster compound, Au55(PPh3)12Cl6; small platinum clusters show similar nonmetallic features.

Inverse surface melting in confined clusters: Ar-13 in zeolite L

Monte Carlo and molecular dynamics simulations on an Ar-13 cluster in zeolite L have been carried out at a series of temperatures to understand the rigid-nonrigid transition corresponding to the solid-liquid transition exhibited by the free Ar-13 cluster. The icosahedral geometry of the free cluster is no longer preferred when the cluster is confined in the zeolite. The root-mean-squared pair distance fluctuation, delta, exhibits a sharp, well-defined rigid-nonrigid transition at 17 K as compared to 27 K for the free cluster. Multiple peaks in the distribution of short-time averages of the guest-host interaction energy indicate coexistence of two phases.; It is shown that this transition is associated with the inner atoms becoming mobile at 17 K even while the outer layer atoms, which are in close proximity to the zeolitic wall, continue to be comparatively immobile. This may be contrasted with the melting of large free clusters of 40 or more atoms which exhibit surface melting. Guest-host interactions seem to play a predominant role in determining the properties of confined clusters. We demonstrate that the volume of the cluster increases rather sharply at 17 and 27 K respectively for the confined and the free cluster. Power spectra suggest that the motion of the inner atoms is generally parallel to the atoms which form the cage wall.

Testing concordance in species boundaries using acoustic, morphological, and molecular data in the field cricket genus Itaropsis (Orthoptera: Grylloidea, Gryllidae: Gryllinae)

In most taxa, species boundaries are inferred based on differences in morphology or DNA sequences revealed by taxonomic or phylogenetic analyses. In crickets, acoustic mating signals or calling songs have species-specific structures and provide a third data set to infer species boundaries. We examined the concordance in species boundaries obtained using acoustic, morphological, and molecular data sets in the field cricket genus Itaropsis. This genus is currently described by only one valid species, Itaropsis tenella, with a broad distribution in western peninsular India and Sri Lanka. Calling songs of males sampled from four sites in peninsular India exhibited significant differences in a number of call features, suggesting the existence of multiple species. Cluster analysis of the acoustic data, molecular phylogenetic analyses, and phylogenetic analyses combining all data sets suggested the existence of three clades. Whatever the differences in calling signals, no full congruence was obtained between all the data sets, even though the resultant lineages were largely concordant with the acoustic clusters. The genus Itaropsis could thus be represented by three morphologically cryptic incipient species in peninsular India; their distributions are congruent with usual patterns of endemism in the Western Ghats, India. Song evolution is analysed through the divergence in syllable period, syllable and call duration, and dominant frequency.

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte-Carlo simulation to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB), which, for a specific set of parameters, sustains two solid phases: honeycomb and oblique. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by increasing temperature. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common belief and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as a linear strip, followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions, which enables the dominance of stabilizing energy over destabilizing surface energy. The nucleus of stable oblique phase is wetted by intermediate order particles, which minimizes the surface free energy. In the case of pressure induced transition at low temperature the collapsed state is a disordered solid. The disordered solid phase has diverse local quasi-stable structures along with oblique-solid like domains. (C) 2013 AIP Publishing LLC.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).

A Reactive Intermediate, Ni-5(C6H4N3)(6)(CO)(4)], in the Formation of Nonameric Clusters of Nickel, Ni-9(C6H4N3)(12)(CO)(6)] and Ni-9(C6H4N3)(12)(CO)(6)].2(C3H7NO)

Three new molecular compounds, Ni-5(bta)(6)(CO)(4)], I, Ni-9(bta)(12)(CO)(6)], II, Ni-9(bta)(12)(CO)(6)]. 2(C3H7NO), III, (bta = benzotriazole) were prepared employing solvothermal reactions. Of these, I have pentanuclear nickel, whereas II and III have nonanuclear nickel species. The structures are formed by the connectivity between the nickel and benzotriazole giving rise to the 5- and 9-membered nickel clusters. The structures are stabilised by extensive pi aEuro broken vertical bar pi and C-H... pi interactions. Compound II and III are solvotamorphs as they have the same 9-membered nickel clusters and have different solvent molecules. To the best of our knowledge, the compounds I-III represent the first examples of the same transition element existing in two distinct coordination environment in this class of compounds. The studies reveal that compound I is reactive and could be an intermediate in the preparation of II and III. Thermal studies indicate that the compounds are stable upto 350(a similar to)C and at higher temperatures (similar to 800(a similar to)C) the compounds decompose into NiO. Magnetic studies reveal that II is anti-ferromagnetic.

Electronically nonadiabatic decomposition mechanisms of clusters of zinc and dimethylnitramine

Electronically nonadiabatic decomposition mechanisms of dimethylnitramine (DMNA) in presence of zinc metal clusters are explored. Complete active space self-consistent field (CASSCF) calculation is employed for DMNA-Zn and ONIOM (Our own N-layered integrated molecular orbital and molecular mechanics) methodology is coupled with CASSCF methodology for DMNA-Zn-10 cluster. Present computational results show that DMNA-Zn clusters undergo electronically nonadiabatic reactions, rendering nitro-nitrite isomerization followed by NO elimination. The overall reactions are also found to be highly exothermic in nature. This is the first report on electronically nonadiabatic decomposition pathways of DMNA-Zn-n neutral clusters. (C) 2014 Elsevier B.V. All rights reserved.

High indium non-polar InGaN clusters with infrared sensitivity grown by PAMBE

Studies on the optical properties of InGaN alloy of relatively higher indium content are of potential interest to understand the effect of indium content on the optical band gap of epitaxial InGaN. We report the growth of self assembled non-polar high indium clusters of In0.55Ga0.45N over non-polar (11-20) a-plane In0.17Ga0.83N epilayer grown on a-plane (11-20) GaN/(1-102) r-plane sapphire substrate using plasma assisted molecular beam epitaxy (PAMBE). Such structures are potential candidates for high brightness LEDs emitting in longer wavelengths. The high resolution X-ray diffraction studies revealed the formation of two distinct compositions of InxGa1-xN alloys, which were further confirmed by photoluminescence studies. A possible mechanism for the formation of such structure was postulated which was supported with the results obtained by energy dispersive X-ray analysis. The structure hence grown when investigated for photo-detecting properties, showed sensitivity to both infrared and ultraviolet radiations due to the different composition of InGaN region. (C) 2015 Author(s).

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.