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.

Temperature Dependence of Mechanical Properties in Molecular Crystals

Quantitative evaluation of the mechanical behavior of molecular materials by a nanoindentation technique has gained prominence recently. However, all the reported data have been on room-temperature properties despite many interesting phenomena observed in them with variations in temperature. In this paper, we report the results of nanoindentation experiments conducted as a function of temperature, T, between 283 and 343 K, on the major faces of three organic crystals: saccharin, sulfathiazole (form 2), and L-alanine, which are distinct in terms of the number and strength of intermolecular interactions in them. Results show that elastic modulus, E, and hardness, H, decrease markedly with increasing T. While E decreases linearly with T, the variations in H with T are not so, and were observed to drop by similar to 50% over the range of T investigated. The slope of the linear fits to E vs T for the organic crystals was found to be around 1, which is considerably higher than the values of 0.3-0.5 reported in the literature for metallic, ionic, and covalently bonded crystalline materials. Possible implications of the observed remarkable changes in H for pharmaceutical manufacturing are highlighted.

Potential histopathological and molecular changes in rat vas deferens inhaled by Boswellia papyrifera and Boswellia carterii

Boswellia papyrifera and Boswellia carterii released from smoke contaminate indoor environment and consequently adversely affect humans as evidenced by respiratory disturbances. The aim of this study was to determine the effects of these plants on pathological and biochemical changes in vas deferens of albino rats. Animals were administered 4g/kg body weight B. papyrifera and B. carterii daily for 120days along with controls. Significant changes were observed in epithelial cell types and some cells showed signs of degeneration. The ultrastructural studies revealed marked changes in cytoplasmic organelles. Microvilli were missing and lysosomes were found in the cytoplasm. In addition, all treated groups plasma fructose and other biochemical parameters were decreased indicating reduced energy necessary for motility and contractility of spermatozoa. Many spermatozoa were disorganized and agglomerated. Data suggest that smoke from these plants adversely affects vas deferens.

Molecular and Functional Characterization of RecD, a Novel Member of the SF1 Family of Helicases, from Mycobacterium tuberculosis

Background: The heterotrimeric M. tuberculosis RecBCD complex, or each of its individual subunits, remains uncharacterized. Results: MtRecD exists as a homodimer in solution, catalyzes ssDNA-dependent ATP hydrolysis, unwinding of DNA replication/recombination intermediates, and interacts with RecA. Conclusion: MtRecD possesses strong 5 3- and weak 3 5-helicase activities. Significance: These findings provide insights into the mechanism underlying DSB repair and homologous recombination in mycobacteria. The annotated whole-genome sequence of Mycobacterium tuberculosis revealed the presence of a putative recD gene; however, the biochemical characteristics of its encoded protein product (MtRecD) remain largely unknown. Here, we show that MtRecD exists in solution as a stable homodimer. Protein-DNA binding assays revealed that MtRecD binds efficiently to single-stranded DNA and linear duplexes containing 5 overhangs relative to the 3 overhangs but not to blunt-ended duplex. Furthermore, MtRecD bound more robustly to a variety of Y-shaped DNA structures having 18-nucleotide overhangs but not to a similar substrate containing 5-nucleotide overhangs. MtRecD formed more salt-tolerant complexes with Y-shaped structures compared with linear duplex having 3 overhangs. The intrinsic ATPase activity of MtRecD was stimulated by single-stranded DNA. Site-specific mutagenesis of Lys-179 in motif I abolished the ATPase activity of MtRecD. Interestingly, although MtRecD-catalyzed unwinding showed a markedly higher preference for duplex substrates with 5 overhangs, it could also catalyze significant unwinding of substrates containing 3 overhangs. These results support the notion that MtRecD is a bipolar helicase with strong 5 3 and weak 3 5 unwinding activities. The extent of unwinding of Y-shaped DNA structures was approximate to 3-fold lower compared with duplexes with 5 overhangs. Notably, direct interaction between MtRecD and its cognate RecA led to inhibition of DNA strand exchange promoted by RecA. Altogether, these studies provide the first detailed characterization of MtRecD and present important insights into the type of DNA structure the enzyme is likely to act upon during the processes of DNA repair or homologous recombination.

Molecular dynamics insight to phase transition in n-alkanes with carbon nanofillers

The present work aims to investigate the phase transition, dispersion and diffusion behavior of nanocomposites of carbon nanotube (CNT) and straight chain alkanes. These materials are potential candidates for organic phase change materials(PCMs) and have attracted flurry of research recently. Accurate experimental evaluation of the mass, thermal and transport properties of such composites is both difficult as well as economically taxing. Additionally it is crucial to understand the factors that results in modification or enhancement of their characteristic at atomic or molecular level. Classical molecular dynamics approach has been extended to elucidate the same. Bulk atomistic models have been generated and subjected to rigorous multistage equilibration. To reaffirm the approach, both canonical and constant-temperature, constant-pressure ensembles were employed to simulate the models under consideration. Explicit determination of kinetic, potential, non-bond and total energy assisted in understanding the enhanced thermal and transport property of the nanocomposites from molecular point of view. Crucial parameters including mean square displacement and simulated self diffusion coefficient precisely define the balance of the thermodynamic and hydrodynamic interactions. Radial distribution function also reflected the density variation, strength and mobility of the nanocomposites. It is expected that CNT functionalization could improve the dispersion within n-alkane matrix. This would further ameliorate the mass and thermal properties of the composite. Additionally, the determined density was in good agreement with experimental data. Thus, molecular dynamics can be utilized as a high throughput technique for theoretical investigation of nanocomposites PCMs. (C) 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Inter molecular azide-diisocyanate coupling: new insights for energetic solid propellants

Hydroxyl terminated azide binders can undergo a spurious reaction with diisocyanates to form tetrazoline-5-one via an inter molecular 1,3-dipolar cycloaddition reaction apart from urethane/allophanate groups which has been overlooked. This has serious implications on solid propellants. The computed activation barrier using density functional theory (DFT) for urethane formation reaction is 28.4 kJ mol(-1) and that for tetrazoline-5-one formation reaction is 108.0 kJ mol(-1). DFT studies reveal that the rate limiting step of the reaction is 1,3-dipolar cycloaddition between azide and isocyanate. A dual cure was observed in the temperature ranges 42-77 degrees C and 78-146 degrees C by differential scanning calorimetry (DSC) and rheological studies, confirming multiple reactions. Tetrazoline-5-one formation was confirmed by Fourier transform infrared spectroscopy (FTIR) and solid state nuclear magnetic resonance spectroscopy (NMR).

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

The annotated whole-genome sequence of Mycobacterium tuberculosis indicated that Rv1388 (Mtihf) likely encodes a putative 20 kDa integration host factor (mIHF). However, very little is known about the functional properties of mIHF or organization of 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 Arg170, Arg171, and Arg173, which might be involved in DNA binding, and a conserved proline (P150) in the tight turn. The phenotypic sensitivity of Escherichia coli Delta ihfA and Delta ihfB strains to UV and methylmethanesulfonate 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, bind 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 alpha beta. Furthermore, we reveal novel insights into IHF-mediated DNA compaction depending on the placement of its preferred binding sites; mIHF promotes compaction of DNA into nucleoid-like or higher-order filamentous structures. We hence 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.

Multicomponent Adducts of Pyridoxine: An Evaluation of the Formation of Eutectics and Molecular Salts

Cocrystallization of pyridoxine (vitamin B6) with several biologically important molecules was undertaken with the intent of successfully designing various hydrogen bonded adducts such as salts, cocrystals, and eutectics. Pyridoxine formed eutectics with isoniazid (an antitubercular drug) and nicotinic acid (vitamin B3) and molecular salts with para-aminobenzoic acid (a bioactive) and saccharin (an artificial sweetener), respectively, in accordance to our design strategy. A salt cocrystal, a precisely conjugate acid-base cocrystal, was obtained for the pyridoxine-para-nitrobenzoic acid combination. The role of supramolecular affinity of hydrogen bonding functional groups and Delta pK(a) differences leading to the formation of above diverse adducts was discussed. This study underpins the need for full-fledged supramolecular compatibility studies of multivitamin/drug combinations toward the development of optimal and/or synergistic combination formulations.

Catalytic pathway, substrate binding and stability in SAICAR synthetase: A structure and molecular dynamics study

The de novo purine biosynthesis is one of the highly conserved pathways among all organisms and is essential for the cell viability. A clear understanding of the enzymes in this pathway would pave way for the development of antimicrobial and anticancer drugs. Phosphoribosylaminoimidazole-succinocar boxamide (SAICAR) synthetase is one of the enzymes in this pathway that catalyzes ATP dependent ligation of carboxyaminoimidazole ribotide (CAIR) with L-aspartate (ASP). Here, we describe eight crystal structures of this enzyme, in C222(1) and H3 space groups, bound to various substrates and substrate mimics from a hyperthermophilic archaea Pyrococcus horikoshii along with molecular dynamics simulations of the structures with substrates. Complexes exhibit minimal deviation from its apo structure. The CAIR binding site displays a preference for pyrimidine nucleotides. In the ADP.TMP-ASP complex, the ASP binds at a position equivalent to that found in Saccharomyces cerevisiae structure (PDB: 2CNU) and thus, clears the ambiguity regarding ASP's position. A possible mode for the inhibition of the enzyme by CTP and UTP, observed earlier in the yeast enzyme, is clearly illustrated in the structures bound to CMP and UMP. The ADP.Mg2+.PO4.CD/MP complex having a phosphate ion between the ATP and CAIR sites strengthens one of the two probable pathways (proposed in Escherichia coli study) of catalytic mechanism and suggests the possibility of a phosphorylation taking place before the ASP's attack on CAIR. Molecular dynamic simulations of this enzyme along with its substrates at 90 degrees C reveal the relative strengths of substrate binding, possible antagonism and the role of Mg2+ ions. (C) 2015 Elsevier Inc. All rights reserved.