Identifying functionally important cis-peptide containing segments in proteins and their utility in molecular function annotation

Cis-peptide embedded segments are rare in proteins but often highlight their important role in molecular function when they do occur. The high evolutionary conservation of these segments illustrates this observation almost universally, although no attempt has been made to systematically use this information for the purpose of function annotation. In the present study, we demonstrate how geometric clustering and level-specific Gene Ontology molecular-function terms (also known as annotations) can be used in a statistically significant manner to identify cis-embedded segments in a protein linked to its molecular function. The present study identifies novel cis-peptide fragments, which are subsequently used for fragment-based function annotation. Annotation recall benchmarks interpreted using the receiver-operator characteristic plot returned an area-under-curve >0.9, corroborating the utility of the annotation method. In addition, we identified cis-peptide fragments occurring in conjunction with functionally important trans-peptide fragments, providing additional insights into molecular function. We further illustrate the applicability of our method in function annotation where homology-based annotation transfer is not possible. The findings of the present study add to the repertoire of function annotation approaches and also facilitate engineering, design and allied studies around the cis-peptide neighborhood of proteins.

Effects of growth temperature on nonpolar a-plane InN grown by molecular beam epitaxy

Nonpolar a-plane InN films were grown on r-plane sapphire substrate by plasma assisted molecular beam epitaxy with GaN underlayer. Effect of growth temperature on structural, morphological, and optical properties has been studied. The growth of nonpolar a-plane (1 1 -2 0) orientation was confirmed by high resolution X-ray diffraction study. The film grown at 500 degrees C shows better crystallinity with the rocking curve FWHM 0.67 degrees and 0.85 degrees along 0 0 0 1] and 1 - 1 0 0] directions, respectively. Scanning electron micrograph shows formation of Indium droplets at higher growth temperature. Room temperature absorption spectra show growth temperature dependent band gap variation from 0.74-0.81 eV, consistent with the expected Burstein-Moss effect. The rectifying behaviour of the I-V curve indicates the existence of Schottky barrier at the InN and GaN interface. (C) 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Salts of hexamethylenetetramine with organic acids: Enhanced anomeric interactions with a lowering of molecular symmetry revealed by crystal structures

The hexamethylenetetramine (HMT) framework displays interesting stereoelectronic interactions of the anomeric type. In the highly symmetrical parent system, the nitrogen centres act as both donors and acceptors. Protonation lowers symmetry and also leads to an enhancement of the anomeric interaction around the protonated centre. X-ray diffraction crystal structures of four derivatives of HMT - with succinic, (DL)-malic, phthalic and 4-hydroxybenzoic acids - reveal significant trends. (The first three form well-defined salts, 4-hydroxybenzoic acid forming a co-crystalline compound.) Each molecular structure is essentially characterised by a major anomeric interaction involving the protonated centre as acceptor. In two cases (succinic and 4-hydroxybenzoic), secondary protonation leads to a weaker anomeric interaction site that apparently competes with the dominant one. Bond length changes indicate that the anomeric interaction decreases as malic > phthalic > succinic > 4-hydroxybenzoic, which correlates with the degree of proton transfer to the nitrogen centre. Along with other bond length and angle changes, the results offer insight into the applicability of the antiperiplanar lone pair hypothesis (ALPH) in a rigid system. (C) 2014 Elsevier B.V. All rights reserved.

Salts of hexamethylenetetramine with organic acids: Enhanced anomeric interactions with a lowering of molecular symmetry revealed by crystal structures

The hexamethylenetetramine (HMT) framework displays interesting stereoelectronic interactions of the anomeric type. In the highly symmetrical parent system, the nitrogen centres act as both donors and acceptors. Protonation lowers symmetry and also leads to an enhancement of the anomeric interaction around the protonated centre. X-ray diffraction crystal structures of four derivatives of HMT - with succinic, (DL)-malic, phthalic and 4-hydroxybenzoic acids - reveal significant trends. (The first three form well-defined salts, 4-hydroxybenzoic acid forming a co-crystalline compound.) Each molecular structure is essentially characterised by a major anomeric interaction involving the protonated centre as acceptor. In two cases (succinic and 4-hydroxybenzoic), secondary protonation leads to a weaker anomeric interaction site that apparently competes with the dominant one. Bond length changes indicate that the anomeric interaction decreases as malic > phthalic > succinic > 4-hydroxybenzoic, which correlates with the degree of proton transfer to the nitrogen centre. Along with other bond length and angle changes, the results offer insight into the applicability of the antiperiplanar lone pair hypothesis (ALPH) in a rigid system. (C) 2014 Elsevier B.V. All rights reserved.

Synthesis, molecular docking and anti-mycobacterial evaluation of new imidazo1,2-a]pyridine-2-carboxamide derivatives

New anti-tubercular agents, imidazo1,2-a]pyridine-2-carboxamide derivatives (5a-q) have been designed and synthesized. The structural considerations of the designed molecules were further supported by the docking study with a long-chain enoyl-acyl carrier protein reductase (InhA). The chemical structures of the new compounds were characterized by IR, H-1 NMR, C-13 NMR, HRMS and elemental analysis. In addition, single crystal X-ray diffraction has also been recorded for compound 5f. Compounds were evaluated in vitro against Mycobacterium tuberculosis H37Rv, and cytotoxicity against HEK-293T cell line. Amongst the tested compounds 5j, 5l and 5q were emerged as good anti-tubercular agents with low cytotoxicity. The structure-anti TB activity relationship of these derivatives was explained by molecular docking. (C) 2014 Elsevier Masson SAS. All rights reserved.

Modulation of Electronic and Self-Assembly Properties of a Donor-Acceptor-Donor-Based Molecular Materials via Atomistic Approach

The performance of molecular materials in optoelectronic devices critically depends upon their electronic properties and solid-state structure. In this report, we have synthesized sulfur and selenium based (T4BT and T4BSe) donor-acceptor-donor (D-A-D) organic derivatives in order to understand the structure-property correlation in organic semiconductors by selectively tuning the chalcogen atom. The photophysical properties exhibit a significant alteration upon varying a single atom in the molecular structure. A joint theoretical and experimental investigation suggests that replacing sulfur with selenium significantly reduces the band gap and molar absorption coefficient because of lower electronegativity and ionization potential of selenium. Single-crystal X-ray diffraction analysis showed differences in their solid-state packing and intermolecular interactions. Subsequently, difference in the solid-state packing results variation in self-assembly. Micorstructural changes within these materials are correlated to their electrical resistance variation, investigated by conducting probe atomic force microscopy (CP-AFM) measurements. These results provide useful guidelines to understand the fundamental properties of D-A-D materials prepared by atomistic modulation.

Molecular effects of encapsulation of glucose oxidase dimer by graphene

Knowing the nature of the enzyme-graphene interface is critical for a design of graphene-based biosensors. Extensive contacts between graphene and enzyme could be obtained by employing a suitable encapsulation which does not impede its enzymatic reaction. We have performed molecular dynamics simulations to obtain an insight on many forms of contact between glucose oxidase dimer and the single-layer graphene nano-sheets. The unconnected graphene sheets tended to form a flat stack regardless of their initial positions around the enzyme, whereas the same graphene sheets linked together formed a flower-like shape engendering different forms of wrapping of the enzyme. During the encapsulation no core hydrophobic residues of the enzyme were exposed. Since the polar and charged amino acids populated the enzyme's surface we also estimated, using DFT calculations, the interaction energies of individual polar and charged amino acid residues with graphene. It was found that the negatively charged residues can bind to graphene unexpectedly strongly; however, the main effect of encapsulation comes from the overlap of adjacent edges of graphene sheets.

Structural studies on Mycobacterium tuberculosis RecA: Molecular plasticity and interspecies variability

Structures of crystals of Mycobacterium tuberculosis RecA, grown and analysed under different conditions, provide insights into hitherto underappreciated details of molecular structure and plasticity. In particular, they yield information on the invariant and variable features of the geometry of the P-loop, whose binding to ATP is central for all the biochemical activities of RecA. The strengths of interaction of the ligands with the P-loop reveal significant differences. This in turn affects the magnitude of the motion of the `switch' residue, Gln195 in M. tuberculosis RecA, which triggers the transmission of ATP-mediated allosteric information to the DNA binding region. M. tuberculosis RecA is substantially rigid compared with its counterparts from M smegmatis and E. coli, which exhibit concerted internal molecular mobility. The interspecies variability in the plasticity of the two mycobacterial proteins is particularly surprising as they have similar sequence and 3D structure. Details of the interactions of ligands with the protein, characterized in the structures reported here, could be useful for design of inhibitors against M. tuberculosis RecA.

Electron transfer statistics and thermal fluctuations in molecular junctions

We derive analytical expressions for probability distribution function (PDF) for electron transport in a simple model of quantum junction in presence of thermal fluctuations. Our approach is based on the large deviation theory combined with the generating function method. For large number of electrons transferred, the PDF is found to decay exponentially in the tails with different rates due to applied bias. This asymmetry in the PDF is related to the fluctuation theorem. Statistics of fluctuations are analyzed in terms of the Fano factor. Thermal fluctuations play a quantitative role in determining the statistics of electron transfer; they tend to suppress the average current while enhancing the fluctuations in particle transfer. This gives rise to both bunching and antibunching phenomena as determined by the Fano factor. The thermal fluctuations and shot noise compete with each other and determine the net (effective) statistics of particle transfer. Exact analytical expression is obtained for delay time distribution. The optimal values of the delay time between successive electron transfers can be lowered below the corresponding shot noise values by tuning the thermal effects. (C) 2015 AIP Publishing LLC.

Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

Response to ``Comment on `Frequency-domain stimulated and spontaneous light emission signals at molecular junctions''' J. Chem. Phys. 142, 137101 (2015)]

In a recent work [U. Harbola, B. K. Agrawalla, and S. Mukamel, J. Chem. Phys. 141, 074107 (2014)], we have presented a superoperator (Liouville space) diagrammatic formulation of spontaneous and stimulated optical signals from current-carrying molecular junctions. We computed the diagrams that contribute to the spontaneous light emission SLE (fluorescence and Raman) signal using a diagrammatic method which clearly distinguishes between the Raman and the fluorescence contributions. We pointed out some discrepancies with the work of Galperin, Ratner and Nitzan (GRN) [M. Galperin, M. A. Ratner and, A. Nitzan, J. Chem. Phys. 130, 144109 (2009)]. In their response [M. Galperin, M. A. Ratner and A. Nitzan, “Comment on‘ Frequency-domain stimulated and spontaneous light emission signals at molecular junctions’” [J. Chem. Phys. 141, 074107 (2014)], J. Chem. Phys. 142, 137101 (2015)] to our work, GRN have argued that there are no differences in the choice of Raman diagrams in both works. Here we reply to their points and show where the differences exist.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

CHEXVIS: a tool for molecular channel extraction and visualization

Background: Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties. Results: We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, CHEXVIS. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that CHEXVIS performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how CHEXVIS can be used to study channels, extract their properties and gain insights into molecular function. Conclusion: CHEXVIS supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The CHEXVIS web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in.

Dynamics of the chemical bond: inter- and intra-molecular hydrogen bond

In this discussion, we show that a static definition of a `bond' is not viable by looking at a few examples for both inter-and intra-molecular hydrogen bonding. This follows from our earlier work (Goswami and Arunan, Phys. Chem. Chem. Phys. 2009, 11, 8974) which showed a practical way to differentiate `hydrogen bonding' from `van der Waals interaction'. We report results from ab initio and atoms in molecules theoretical calculations for a series of Rg center dot center dot center dot HX complexes (Rg = He/Ne/Ar and X = F/Cl/Br) and ethane-1,2-diol. Results for the Rg center dot center dot center dot HX/DX complexes show that Rg center dot center dot center dot DX could have a `deuterium bond' even when Rg center dot center dot center dot HX is not `hydrogen bonded', according to the practical criterion given by Goswami and Arunan. Results for ethane-1,2-diol show that an `intra-molecular hydrogen bond' can appear during a normal mode vibration which is dominated by the O center dot center dot center dot O stretching, though a `bond' is not found in the equilibrium structure. This dynamical `bond' formation may nevertheless be important in ensuring the continuity of electron density across a molecule. In the former case, a vibration `breaks' an existing bond and in the later case, a vibration leads to `bond' formation. In both cases, the molecule/complex stays bound irrespective of what happens to this `hydrogen bond'. Both these cases push the borders on the recent IUPAC recommendation on hydrogen bonding (Arunan et al. Pure. Appl. Chem. 2011, 83 1637) and justify the inclusive nature of the definition.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.