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.

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).