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.

Optoelectronic Properties of Graphene on Silicon Substrate: Effect of Defects in Graphene

Engineering of electronic energy band structure in graphene based nanostructures has several potential applications. Substrate induced bandgap opening in graphene results several optoelectronic properties due to the inter-band transitions. Various defects like structures, including Stone-Walls and higher-order defects are observed when a graphene sheet is exfoliated from graphite and in many other growth conditions. Existence of defect in graphene based nanostructures may cause changes in optoelectronic properties. Defect engineered graphene on silicon system are considered in this paper to study the tunability of optoelectronic properties. Graphene on silicon atomic system is equilibrated using molecular dynamics simulation scheme. Based on this study, we confirm the existence of a stable super-lattice. Density functional calculations are employed to determine the energy band structure for the super-lattice. Increase in the optical energy bandgap is observed with increasing of order of the complexity in the defect structure. Optical conductivity is computed as a function of incident electromagnetic energy which is also increasing with increase in the defect order. Tunability in optoelectronic properties will be useful in understanding graphene based design of photodetectors, photodiodes and tunnelling transistors.

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.

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.

Multifarious facets of sugar-derived molecular gels: molecular features, mechanisms of self-assembly and emerging applications

The remarkable capability of nature to design and create excellent self-assembled nano-structures, especially in the biological world, has motivated chemists to mimic such systems with synthetic molecular and supramolecular systems. The hierarchically organized self-assembly of low molecular weight gelators (LMWGs) based on non-covalent interactions has been proven to be a useful tool in the development of well-defined nanostructures. Among these, the self-assembly of sugar-derived LMWGs has received immense attention because of their propensity to furnish biocompatible, hierarchical, supramolecular architectures that are macroscopically expressed in gel formation. This review sheds light on various aspects of sugar-derived LMWGs, uncovering their mechanisms of gelation, structural analysis, and tailorable properties, and their diverse applications such as stimuli-responsiveness, sensing, self-healing, environmental problems, and nano and biomaterials synthesis.

Anisotropy in the mechanical properties of organic crystals: temperature dependence

The nanoindentation technique has recently been utilized for quantitative evaluation of the mechanical properties of molecular materials successfully, including their temperature (T) dependence. In this paper, we examine how the mechanical anisotropy varies with T in saccharin and L-alanine single crystals. Our results show that elastic modulus (E) decreases linearly in all the cases examined, with the T-dependence of E being anisotropic. Correspondence between directional dependence of the slopes of the E vs. T plots and the linear thermal expansion coefficients was found. The T-dependence of hardness (H), on the other hand, was found to be nonlinear and significant when (100) of saccharin and (001) of L-alanine are indented. While the anisotropies in E and H of saccharin and E of L-alanine enhance with T, the anisotropy in H of L-alanine was found to reduce with T. Possible mechanistic origins of these variations are discussed.

Graphene scavenges free radicals to synergistically enhance structural properties in a gamma-irradiated polyethylene composite through enhanced interfacial interactions

A unique strategy for scavenging free radicals in situ on exposure to gamma irradiation in polyethylene (PE) nanocomposites is presented. Blends of ultra-high molecular weight PE and linear low-density PE (PEB) and their nanocomposites with graphene (GPEB) were prepared by melt mixing to develop materials for biomedical implants. The effect of gamma irradiation on the microstructure and mechanical properties was systematically investigated. The neat blend and the nanocomposite were subjected to gamma-ray irradiation in order to improve the interfacial adhesion between PE and graphene sheets. Structural and thermal characterization revealed that irradiation induced crosslinking and increased the crystallinity of the polymer blend. The presence of graphene further enhanced the crystallinity via crosslinks between the polymer matrix and the filler on irradiation. Graphene was found to scavenge free radicals as confirmed by electron paramagnetic resonance spectroscopy. Irradiation of graphene-containing polymer composites resulted in the largest increase in modulus and hardness compared to either irradiation or addition of graphene to PEB alone. This study provides new insight into the role of graphene in polymer matrices during irradiation and suggests that irradiated graphene-polymer composites could emerge as promising materials for use as articulating surfaces in biomedical implants.

Synthesis of Eu3+-activated ZnO superstructures: Photoluminescence, Judd-Ofelt analysis and Sunlight photocatalytic properties

Photoluminescence and photocatalytic properties of Eu-doped ZnO nanoparticles (NPs) were synthesized by facile phyto route. XPS results demonstrated the existence of Eu3+ as dopant into ZnO. Morphologies of the NPs were mainly dependent on Eu3+ and Aloe vera gel. Red shift of energy band gap was due to the creation of intermediate energy states of Eu3+ and oxygen vacancies in the band gap. PL emission of ZnO:Eu3+ (1-11 mol%, 8 ml and 7 mol%, 2-12 ml) exhibit characteristic peaks of D-5(0) -> F-7(2) transitions. From the Judd-Ofelt analysis, intensities of transitions between different.' levels dependent on the symmetry of the local environment of Eu3+ ions. CIE chromaticity co-ordinates confirm reddish emission of the phosphor. Further, NPs exhibit excellent photocatalytic activity for the degradation of Rhodamine B (94%) under Sunlight was attributed to crystallite size, band gap, morphology and oxygen vacancies. In addition, photocatalyst reusability studies were conducted and found that Eu-doped catalyst could be reused several times with negligible decrease in catalytic activity. The present work directs new possibilities to provide some new insights into the design of new phyto synthesized nanophosphors for display devices, photocatalysts with high activity for environmental clean-up and solar energy conversion. (C) 2015 Elsevier B.V. All rights reserved.

Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Dispersible crystalline nanobundles of YPO4 and Ln (Eu, Tb)-doped YPO4: rapid synthesis, optical properties and bio-probe applications

Undoped and Ln(3+) (Eu and Tb)-doped crystalline nanobundles of YPO4 were prepared by a facile microwave-assisted route with water as a solvent and without using any surfactant. TEM investigations reveal that the as-prepared powder consists of lenticular-shaped nanobundles (similar to 100 nm in diameter) made of very small nanorods with diameter less than 10 nm and length varying from 20 to 50 nm. Each nanorod in turn is single crystalline, as revealed by HRTEM imaging. The as-prepared nanobundles are easily dispersible in various solvents, especially water, without any surface functionalization, which is critical for various bio-probe applications like cell and tissue imaging. The Eu- and Tb-doped YPO4 nanobundles show good photoluminescence properties and were further evaluated for their use as fluorescent biolabels. Our results show that HeLa cells labelled with Eu- and Tb-doped YPO4 nanobundles show bright red (Eu) and green (Tb) intracellular luminescence under a confocal microscope. Concentration-and time-dependent MTT cell viability assays show that the nanobundles show low toxicity towards cells which makes them promising in bioimaging field.

Bioinspired Reductionistic Peptide Engineering for Exceptional Mechanical Properties

A simple solution-processing and self-assembly approach that exploits the synergistic interactions between multiple hydrogen bonded networks and aromatic interactions was utilized to synthesize molecular crystals of cyclic dipeptides (CDPs), whose molecular weights (similar to 0.2 kDa) are nearly three orders of magnitude smaller than that of natural structural proteins (50-300 kDa). Mechanical properties of these materials, measured using the nanoindentation technique, indicate that the stiffness and strength are comparable and sometimes better than those of natural fibres. The measured mechanical responses were rationalized by recourse to the crystallographic structural analysis and intermolecular interactions in the self-assembled single crystals. With this work we highlight the significance of developing small molecule based bioinspired design strategies to emulate biomechanical properties. A particular advantage of the successfully demonstrated reductionistic strategy of the present work is its amenability for realistic industrial scale manufacturing of designer biomaterials with desired mechanical properties.

Crystal chemistry and photomechanical behavior of 3,4-dimethoxycinnamic acid: correlation between maximum yield in the solid-state topochemical reaction and cooperative molecular motion

A new monoclinic polymorph, form II (P2(1)/c, Z = 4), has been isolated for 3,4-dimethoxycinnamic acid (DMCA). Its solid-state 2 + 2 photoreaction to the corresponding alpha-truxillic acid is different from that of the first polymorph, the triclinic form I (P (1) over bar, Z = 4) that was reported in 1984. The crystal structures of the two forms are rather different. The two polymorphs also exhibit different photomechanical properties. Form I exhibits photosalient behavior but this effect is absent in form II. These properties can be explained on the basis of the crystal packing in the two forms. The nanoindentation technique is used to shed further insights into these structure-property relationships. A faster photoreaction in form I and a higher yield in form II are rationalized on the basis of the mechanical properties of the individual crystal forms. It is suggested that both Schmidt-type and Kaupp-type topochemistry are applicable for the solid-state trans-cinnamic acid photodimerization reaction. Form I of DMCA is more plastic and seems to react under Kaupp-type conditions with maximum molecular movements. Form II is more brittle, and its interlocked structure seems to favor Schmidt-type topochemistry with minimum molecular movement.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).

Targeting human telomeric G-quadruplex DNA with curcumin and its synthesized analogues under molecular crowding conditions

The formation of telomeric G-quadruplexes has been shown to inhibit telomerase activity. Indeed, a number of small molecules capable of p-stacking with G-tetrads have shown the ability to inhibit telomerase activity through the stabilization of G-quadruplexes. Curcumin displays a wide spectrum of medicinal properties ranging from anti-bacterial, anti-viral, anti-protozoal, anti-fungal and anti-inflammatory to anti-cancer activity. We have investigated the interactions of curcumin and its structural analogues with the human telomeric sequence AG(3)(T(2)AG(3))(3) under molecular crowding conditions. Experimental studies indicated the existence of a AG(3)(T(2)AG(3))(3)/curcumin complex with binding affinity of 0.72 x 10(6) M-1 under molecular crowding conditions. The results from UV-visible absorption spectroscopy, a fluorescent TO displacement assay, circular dichroism and molecular docking studies, imply that curcumin and their analogues interact with G-quadruplex DNA via groove binding. While other analogs of curcumin studied here bind to G-quadruplexes in a qualitatively similar manner their affinities are relatively lower in comparison to curcumin. The Knoevenagel condensate, a methoxy-benzylidene derivative of curcumin, also exhibited significant binding to G-quadruplex DNA, although with two times decreased affinity. Our study establishes the potential of curcumin as a promising natural product for G-quadruplex specific ligands.