Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Herringbone to cofacial solid state packing via H-bonding in diketopyrrolopyrrole (DPP) based molecular crystals: influence on charge transport

Themono-alkylation of DPP derivatives leads to cofacial pi-pi stacking via H-bonding unlike their di-alkylated counterparts, which exhibit a classical herringbone packing pattern. Single crystal organic field-effect transistor (OFET) measurements reveal a significant enhancement of charge carrier mobility for mono-hexyl DPP derivatives.

Molecular Dissection of Mycobacterium tuberculosis Integration Host Factor Reveals Novel Insights into the Mode of DNA Binding and Nucleoid Compaction

Background: mIHF belongs to a subfamily of proteins, distinct from E. coli IHF. Results: Functionally important amino acids of mIHF and the mechanism(s) underlying DNA binding, DNA bending, and site-specific recombination are distinct from that of E. coli IHF. Conclusion: mIHF functions could contribute beyond nucleoid compaction. Significance: Because mIHF is essential for growth, the molecular mechanisms identified here can be exploited in drug screening efforts. The annotated whole-genome sequence of Mycobacterium tuberculosis revealed that Rv1388 (Mtihf) is likely to encode for a putative 20-kDa integration host factor (mIHF). However, very little is known about the functional properties of mIHF or the organization of the mycobacterial nucleoid. Molecular modeling of the mIHF three-dimensional structure, based on the cocrystal structure of Streptomyces coelicolor IHF duplex DNA, a bona fide relative of mIHF, revealed the presence of Arg-170, Arg-171, and Arg-173, which might be involved in DNA binding, and a conserved proline (Pro-150) in the tight turn. The phenotypic sensitivity of Escherichia coli ihfA and ihfB strains to UV and methyl methanesulfonate could be complemented with the wild-type Mtihf but not its alleles bearing mutations in the DNA-binding residues. Protein-DNA interaction assays revealed that wild-type mIHF, but not its DNA-binding variants, binds with high affinity to fragments containing attB and attP sites and curved DNA. Strikingly, the functionally important amino acid residues of mIHF and the mechanism(s) underlying its binding to DNA, DNA bending, and site-specific recombination are fundamentally different from that of E. coli IHF. Furthermore, we reveal novel insights into IHF-mediated DNA compaction depending on the placement of its preferred binding sites; mIHF promotes DNA compaction into nucleoid-like or higher order filamentous structures. We therefore propose that mIHF is a distinct member of a subfamily of proteins that serve as essential cofactors in site-specific recombination and nucleoid organization and that these findings represent a significant advance in our understanding of the role(s) of nucleoid-associated proteins.

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.